Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 21:17:15 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 Tl.rel-pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 6S 5D 5D renormalized file Y.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4S 4P 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 28 20 6 690 390 68 Max 29 21 7 693 413 71 Sum 1041 725 221 24885 14531 2469 bravais-lattice index = 14 lattice parameter (alat) = 7.0884 a.u. unit-cell volume = 356.1539 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 24.00 number of Kohn-Sham states= 32 kinetic-energy cutoff = 45.0000 Ry charge density cutoff = 258.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)= 7.088362 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 Tl read from file: /users/gautes/Pseudo/Tl.rel-pbe-dn-rrkjus_psl.1.0.0.UPF MD5 check sum: cd969aa200d2b9b89aab382fbbb9d684 Pseudo is Ultrasoft + core correction, Zval = 13.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 Y read from file: /users/gautes/Pseudo/Y.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 25a826a599bbcba7ebf619f9598d90b1 Pseudo is Ultrasoft + core correction, Zval = 11.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1223 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 Tl 13.00 204.38330 Tl( 1.00) Y 11.00 88.90590 Y( 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= 20 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 Dense grid: 24885 G-vectors FFT dimensions: ( 40, 40, 40) Smooth grid: 14531 G-vectors FFT dimensions: ( 32, 32, 32) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.06 Mb ( 114, 32) NL pseudopotentials 0.06 Mb ( 57, 68) Each V/rho on FFT grid 0.05 Mb ( 3200) Each G-vector array 0.01 Mb ( 693) G-vector shells 0.00 Mb ( 175) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.22 Mb ( 114, 128) Each subspace H/S matrix 0.02 Mb ( 32, 32) Each matrix 0.07 Mb ( 68, 2, 32) Arrays for rho mixing 0.39 Mb ( 3200, 8) Check: negative/imaginary core charge= -0.000002 0.000000 Initial potential from superposition of free atoms starting charge 23.99395, renormalised to 24.00000 Starting wfc are 44 randomized atomic wfcs total cpu time spent up to now is 2.7 secs per-process dynamical memory: 24.7 Mb Self-consistent Calculation iteration # 1 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 2.0 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.22E-04, avg # of iterations = 3.0 total cpu time spent up to now is 3.9 secs total energy = -200.03201041 Ry Harris-Foulkes estimate = -200.05498138 Ry estimated scf accuracy < 0.03379670 Ry iteration # 2 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.41E-04, avg # of iterations = 2.0 total cpu time spent up to now is 4.6 secs total energy = -200.03855932 Ry Harris-Foulkes estimate = -200.04890043 Ry estimated scf accuracy < 0.01761227 Ry iteration # 3 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.34E-05, avg # of iterations = 2.0 total cpu time spent up to now is 5.2 secs total energy = -200.04355681 Ry Harris-Foulkes estimate = -200.04405582 Ry estimated scf accuracy < 0.00103151 Ry iteration # 4 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.30E-06, avg # of iterations = 3.4 total cpu time spent up to now is 5.9 secs total energy = -200.04374959 Ry Harris-Foulkes estimate = -200.04375191 Ry estimated scf accuracy < 0.00002410 Ry iteration # 5 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-07, avg # of iterations = 3.0 total cpu time spent up to now is 6.7 secs total energy = -200.04375670 Ry Harris-Foulkes estimate = -200.04375769 Ry estimated scf accuracy < 0.00000246 Ry iteration # 6 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.02E-08, avg # of iterations = 2.4 total cpu time spent up to now is 7.3 secs total energy = -200.04375708 Ry Harris-Foulkes estimate = -200.04375710 Ry estimated scf accuracy < 0.00000007 Ry iteration # 7 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.07E-10, avg # of iterations = 3.6 total cpu time spent up to now is 8.2 secs total energy = -200.04375711 Ry Harris-Foulkes estimate = -200.04375713 Ry estimated scf accuracy < 0.00000002 Ry iteration # 8 ecut= 45.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.04E-10, avg # of iterations = 2.6 total cpu time spent up to now is 8.9 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 1839 PWs) bands (ev): -30.6744 -30.6744 -11.5937 -11.5937 -10.1441 -10.1441 -10.1441 -10.1441 -0.4200 -0.4200 -0.4200 -0.4200 1.6733 1.6733 1.6733 1.6733 1.8465 1.8465 3.3894 3.3894 12.2510 12.2510 12.2510 12.2510 14.3011 14.3011 15.8578 15.8578 15.8578 15.8578 15.9694 15.9694 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.1429 ( 1823 PWs) bands (ev): -30.6719 -30.6719 -11.6046 -11.6046 -10.1741 -10.1741 -10.1412 -10.1412 -0.4175 -0.4175 -0.4039 -0.4039 1.6788 1.6788 1.7010 1.7010 1.8402 1.8402 3.6457 3.6457 11.8038 11.8038 12.2744 12.2744 13.8600 13.8600 14.6691 14.6691 14.8451 14.8451 15.8756 15.8756 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 0.9998 0.9998 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.2857 ( 1811 PWs) bands (ev): -30.6664 -30.6664 -11.6310 -11.6310 -10.2380 -10.2380 -10.1348 -10.1348 -0.4122 -0.4122 -0.3711 -0.3711 1.6898 1.6898 1.7633 1.7633 1.8280 1.8280 4.3742 4.3742 10.1580 10.1580 12.3296 12.3296 13.2495 13.2495 13.7320 13.7320 15.1474 15.1474 15.8323 15.8323 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 0.9894 0.9894 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4286 ( 1823 PWs) bands (ev): -30.6619 -30.6619 -11.6541 -11.6541 -10.2858 -10.2858 -10.1296 -10.1296 -0.4079 -0.4079 -0.3479 -0.3479 1.6965 1.6965 1.8138 1.8138 1.8206 1.8206 5.3274 5.3274 8.5051 8.5051 12.3769 12.3769 12.8043 12.8043 13.2364 13.2364 15.7979 15.7979 16.5012 16.5012 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 0.7439 0.7439 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.1429 ( 1833 PWs) bands (ev): -30.6695 -30.6695 -11.6147 -11.6147 -10.1873 -10.1873 -10.1538 -10.1538 -0.4107 -0.4107 -0.3966 -0.3966 1.6919 1.6919 1.7163 1.7163 1.8292 1.8292 3.8895 3.8895 11.7809 11.7809 12.0302 12.0302 13.2416 13.2416 14.1695 14.1696 15.6142 15.6142 15.6929 15.6929 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.2857 ( 1823 PWs) bands (ev): -30.6640 -30.6640 -11.6391 -11.6391 -10.2437 -10.2437 -10.1517 -10.1517 -0.4072 -0.4072 -0.3719 -0.3719 1.7003 1.7003 1.7696 1.7696 1.8100 1.8100 4.5715 4.5715 10.3097 10.3097 12.1457 12.1457 13.0499 13.0499 13.6733 13.6733 15.5212 15.5212 15.6148 15.6148 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.4286 ( 1820 PWs) bands (ev): -30.6596 -30.6596 -11.6604 -11.6604 -10.2895 -10.2895 -10.1462 -10.1462 -0.4059 -0.4059 -0.3538 -0.3538 1.7019 1.7019 1.7885 1.7885 1.8262 1.8262 5.4275 5.4275 8.7974 8.7974 12.3972 12.3972 12.7206 12.7206 13.3206 13.3206 15.5464 15.5464 16.5572 16.5572 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 0.3958 0.3958 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.2857 ( 1827 PWs) bands (ev): -30.6587 -30.6587 -11.6588 -11.6588 -10.2738 -10.2738 -10.1692 -10.1692 -0.4045 -0.4045 -0.3692 -0.3692 1.7103 1.7103 1.7477 1.7477 1.8262 1.8262 5.1007 5.1007 10.3989 10.3989 11.3820 11.3820 13.1144 13.1144 13.8538 13.8538 15.1217 15.1217 16.1764 16.1764 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.4286 ( 1816 PWs) bands (ev): -30.6544 -30.6544 -11.6763 -11.6763 -10.3087 -10.3087 -10.1693 -10.1693 -0.4081 -0.4081 -0.3651 -0.3651 1.7004 1.7004 1.7253 1.7253 1.8669 1.8669 5.6955 5.6955 9.5840 9.5840 11.6151 11.6151 12.9451 12.9451 13.8911 13.8911 15.0856 15.0856 16.6227 16.6227 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4286 0.4286 ( 1808 PWs) bands (ev): -30.6502 -30.6502 -11.6907 -11.6907 -10.3340 -10.3340 -10.1751 -10.1751 -0.4155 -0.4155 -0.3728 -0.3728 1.6656 1.6656 1.7041 1.7041 1.9097 1.9097 5.9821 5.9821 10.3732 10.3732 11.0408 11.0408 12.6598 12.6598 14.5353 14.5353 14.8073 14.8073 16.6910 16.6910 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.1429 ( 1827 PWs) bands (ev): -30.6671 -30.6671 -11.6241 -11.6241 -10.1869 -10.1869 -10.1789 -10.1789 -0.3991 -0.3991 -0.3965 -0.3965 1.7182 1.7182 1.7184 1.7184 1.8146 1.8146 4.1178 4.1178 11.7817 11.7817 12.1018 12.1018 12.6191 12.6191 14.2293 14.2293 14.4921 14.4921 16.6087 16.6087 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.2857 ( 1820 PWs) bands (ev): -30.6616 -30.6616 -11.6469 -11.6469 -10.2447 -10.2447 -10.1734 -10.1734 -0.4000 -0.4000 -0.3728 -0.3728 1.7187 1.7187 1.7757 1.7757 1.7912 1.7912 4.7536 4.7536 10.4955 10.4955 12.0880 12.0880 12.9155 12.9155 13.4331 13.4331 14.8537 14.8537 16.2247 16.2247 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.4286 ( 1814 PWs) bands (ev): -30.6573 -30.6573 -11.6669 -11.6669 -10.2906 -10.2906 -10.1661 -10.1661 -0.4022 -0.4022 -0.3560 -0.3560 1.7125 1.7125 1.7726 1.7726 1.8308 1.8308 5.5342 5.5342 9.0638 9.0638 12.1747 12.1747 12.8687 12.8687 13.4356 13.4356 15.1527 15.1527 16.6519 16.6519 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.2857 0.2857 ( 1824 PWs) bands (ev): -30.6563 -30.6563 -11.6661 -11.6661 -10.2688 -10.2688 -10.1971 -10.1971 -0.3949 -0.3949 -0.3682 -0.3682 1.7332 1.7332 1.7479 1.7479 1.8214 1.8214 5.2535 5.2535 10.6452 10.6452 11.3464 11.3464 13.0627 13.0627 13.4425 13.4425 14.5732 14.5732 16.3908 16.3908 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.2857 0.4286 ( 1812 PWs) bands (ev): -30.6520 -30.6520 -11.6832 -11.6832 -10.3051 -10.3051 -10.1961 -10.1961 -0.3990 -0.3990 -0.3605 -0.3605 1.7007 1.7007 1.7510 1.7510 1.8673 1.8673 5.8300 5.8300 9.7406 9.7406 11.4174 11.4174 13.1928 13.1928 13.8084 13.8084 14.7612 14.7612 16.0945 16.0945 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.4286 0.4286 ( 1809 PWs) bands (ev): -30.6478 -30.6478 -11.6978 -11.6978 -10.3300 -10.3300 -10.2041 -10.2041 -0.4001 -0.4001 -0.3647 -0.3647 1.6663 1.6663 1.7472 1.7472 1.9121 1.9121 6.1590 6.1590 10.2690 10.2690 10.8753 10.8753 12.9370 12.9370 14.3878 14.3878 15.0098 15.0098 15.3351 15.3351 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.2857 0.2857 ( 1809 PWs) bands (ev): -30.6510 -30.6510 -11.6840 -11.6840 -10.2598 -10.2598 -10.2548 -10.2548 -0.3825 -0.3825 -0.3615 -0.3615 1.7095 1.7095 1.8130 1.8130 1.8229 1.8229 5.6913 5.6913 10.7454 10.7454 11.3471 11.3471 12.1850 12.1850 14.1216 14.1216 14.1433 14.1433 15.6750 15.6750 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.2857 0.4286 ( 1809 PWs) bands (ev): -30.6468 -30.6468 -11.7000 -11.7000 -10.3025 -10.3025 -10.2481 -10.2481 -0.3816 -0.3816 -0.3515 -0.3515 1.6835 1.6835 1.8182 1.8182 1.8791 1.8791 6.2200 6.2200 10.0067 10.0067 10.9562 10.9562 13.0464 13.0464 14.0103 14.0103 14.5016 14.5016 15.2887 15.2887 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.4286 0.4286 ( 1805 PWs) bands (ev): -30.6426 -30.6426 -11.7151 -11.7151 -10.3210 -10.3210 -10.2665 -10.2665 -0.3685 -0.3685 -0.3476 -0.3476 1.6663 1.6663 1.8458 1.8458 1.9196 1.9196 6.6932 6.6932 9.9582 9.9582 10.5386 10.5386 12.4658 12.4658 14.0287 14.0287 15.3787 15.3787 15.4314 15.4315 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 0.0042 0.0042 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4286 0.4286 0.4286 ( 1814 PWs) bands (ev): -30.6384 -30.6384 -11.7304 -11.7304 -10.3139 -10.3139 -10.3138 -10.3138 -0.3448 -0.3448 -0.3362 -0.3362 1.6641 1.6641 1.9251 1.9251 1.9287 1.9287 7.4409 7.4409 9.6311 9.6311 10.2530 10.2530 10.8945 10.8945 15.6870 15.6870 15.7644 15.7644 15.7670 15.7670 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 12.3914 ev ! total energy = -200.04375712 Ry Harris-Foulkes estimate = -200.04375712 Ry estimated scf accuracy < 5.6E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -33.08409152 Ry hartree contribution = 31.65312672 Ry xc contribution = -50.17641477 Ry ewald contribution = -148.43625778 Ry smearing contrib. (-TS) = -0.00011976 Ry convergence has been achieved in 8 iterations Writing output data file YTl.save init_run : 0.36s CPU 0.41s WALL ( 1 calls) electrons : 6.06s CPU 6.24s WALL ( 1 calls) Called by init_run: wfcinit : 0.20s CPU 0.23s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 4.73s CPU 4.85s WALL ( 9 calls) sum_band : 1.03s CPU 1.05s WALL ( 9 calls) v_of_rho : 0.02s CPU 0.02s WALL ( 9 calls) v_h : 0.00s CPU 0.00s WALL ( 9 calls) v_xc : 0.02s CPU 0.02s WALL ( 9 calls) newd : 0.30s CPU 0.30s WALL ( 9 calls) mix_rho : 0.01s CPU 0.01s WALL ( 9 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.01s WALL ( 380 calls) cegterg : 4.52s CPU 4.59s WALL ( 180 calls) Called by sum_band: sum_band:bec : 0.34s CPU 0.36s WALL ( 180 calls) addusdens : 0.15s CPU 0.15s WALL ( 9 calls) Called by *egterg: h_psi : 2.40s CPU 2.50s WALL ( 682 calls) s_psi : 0.14s CPU 0.14s WALL ( 682 calls) g_psi : 0.00s CPU 0.01s WALL ( 482 calls) cdiaghg : 1.81s CPU 1.75s WALL ( 642 calls) cegterg:over : 0.11s CPU 0.12s WALL ( 482 calls) cegterg:upda : 0.11s CPU 0.09s WALL ( 482 calls) cegterg:last : 0.04s CPU 0.04s WALL ( 180 calls) cdiaghg:chol : 0.09s CPU 0.09s WALL ( 642 calls) cdiaghg:inve : 0.02s CPU 0.02s WALL ( 642 calls) cdiaghg:para : 0.14s CPU 0.12s WALL ( 1284 calls) Called by h_psi: h_psi:vloc : 2.07s CPU 2.14s WALL ( 682 calls) h_psi:vnl : 0.33s CPU 0.36s WALL ( 682 calls) add_vuspsi : 0.14s CPU 0.19s WALL ( 682 calls) General routines calbec : 0.26s CPU 0.22s WALL ( 862 calls) fft : 0.03s CPU 0.04s WALL ( 273 calls) ffts : 0.01s CPU 0.00s WALL ( 72 calls) fftw : 2.24s CPU 2.36s WALL ( 72880 calls) interpolate : 0.02s CPU 0.02s WALL ( 72 calls) Parallel routines fft_scatter : 1.46s CPU 1.50s WALL ( 73225 calls) PWSCF : 8.66s CPU 9.77s WALL This run was terminated on: 21:17:25 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=