Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 19:50:33 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 In.rel-pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 5S 4D 4D 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 18 5 698 352 57 Max 29 19 6 703 361 62 Sum 1041 665 193 25269 12893 2103 bravais-lattice index = 14 lattice parameter (alat) = 7.1148 a.u. unit-cell volume = 360.1567 (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 = 41.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.114818 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 In read from file: /users/gautes/Pseudo/In.rel-pbe-dn-rrkjus_psl.1.0.0.UPF MD5 check sum: 52ff8a85f9af3f41268f91f861ab85bf 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 1241 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 In 13.00 114.81800 In( 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: 25269 G-vectors FFT dimensions: ( 40, 40, 40) Smooth grid: 12893 G-vectors FFT dimensions: ( 30, 30, 30) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.05 Mb ( 94, 32) NL pseudopotentials 0.05 Mb ( 47, 68) Each V/rho on FFT grid 0.05 Mb ( 3200) Each G-vector array 0.01 Mb ( 703) G-vector shells 0.00 Mb ( 188) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.18 Mb ( 94, 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.000003 0.000000 Initial potential from superposition of free atoms starting charge 23.99300, renormalised to 24.00000 Starting wfc are 44 randomized atomic wfcs total cpu time spent up to now is 2.5 secs per-process dynamical memory: 24.7 Mb Self-consistent Calculation iteration # 1 ecut= 41.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 = 8.42E-05, avg # of iterations = 3.2 total cpu time spent up to now is 3.8 secs total energy = -223.74233249 Ry Harris-Foulkes estimate = -223.75595886 Ry estimated scf accuracy < 0.02264751 Ry iteration # 2 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.44E-05, avg # of iterations = 2.0 total cpu time spent up to now is 4.5 secs total energy = -223.74694960 Ry Harris-Foulkes estimate = -223.75000925 Ry estimated scf accuracy < 0.00484359 Ry iteration # 3 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.02E-05, avg # of iterations = 2.0 total cpu time spent up to now is 5.1 secs total energy = -223.74827507 Ry Harris-Foulkes estimate = -223.74849632 Ry estimated scf accuracy < 0.00042032 Ry iteration # 4 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.75E-06, avg # of iterations = 3.9 total cpu time spent up to now is 6.0 secs total energy = -223.74838470 Ry Harris-Foulkes estimate = -223.74839321 Ry estimated scf accuracy < 0.00002841 Ry iteration # 5 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.18E-07, avg # of iterations = 2.6 total cpu time spent up to now is 6.6 secs total energy = -223.74838973 Ry Harris-Foulkes estimate = -223.74838952 Ry estimated scf accuracy < 0.00000017 Ry iteration # 6 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.18E-10, avg # of iterations = 4.0 total cpu time spent up to now is 7.5 secs total energy = -223.74838985 Ry Harris-Foulkes estimate = -223.74838984 Ry estimated scf accuracy < 0.00000002 Ry iteration # 7 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.43E-11, avg # of iterations = 2.9 total cpu time spent up to now is 8.2 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 1575 PWs) bands (ev): -31.6415 -31.6415 -12.5625 -12.5625 -11.1141 -11.1141 -11.1141 -11.1141 -3.5607 -3.5607 -3.5607 -3.5607 -2.7265 -2.7265 -2.7265 -2.7265 -2.6719 -2.6719 3.3142 3.3142 11.0079 11.0079 11.0079 11.0079 13.9056 13.9056 14.4672 14.4672 14.4672 14.4672 14.7408 14.7408 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 ( 1599 PWs) bands (ev): -31.6393 -31.6393 -12.5728 -12.5728 -11.1421 -11.1421 -11.1111 -11.1111 -3.5598 -3.5598 -3.5524 -3.5524 -2.7250 -2.7250 -2.7141 -2.7141 -2.6740 -2.6740 3.5712 3.5712 10.6455 10.6455 11.0818 11.0818 13.1918 13.1918 13.4378 13.4379 13.8056 13.8056 14.7783 14.7783 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.2857 ( 1598 PWs) bands (ev): -31.6341 -31.6341 -12.5976 -12.5976 -11.2021 -11.2021 -11.1045 -11.1045 -3.5577 -3.5577 -3.5354 -3.5354 -2.7224 -2.7224 -2.6845 -2.6845 -2.6781 -2.6781 4.3041 4.3041 9.2251 9.2251 11.2633 11.2633 12.4004 12.4004 12.5665 12.5665 14.2219 14.2219 14.8027 14.8027 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.9997 0.9997 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4286 ( 1602 PWs) bands (ev): -31.6299 -31.6299 -12.6192 -12.6192 -11.2470 -11.2470 -11.0992 -11.0992 -3.5561 -3.5561 -3.5234 -3.5234 -2.7209 -2.7209 -2.6807 -2.6807 -2.6590 -2.6590 5.3343 5.3343 7.6080 7.6080 11.4284 11.4284 11.9556 11.9556 12.1013 12.1013 14.8226 14.8226 15.5145 15.5145 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.0177 0.0177 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.1429 ( 1604 PWs) bands (ev): -31.6370 -31.6370 -12.5824 -12.5824 -11.1544 -11.1544 -11.1229 -11.1229 -3.5562 -3.5562 -3.5489 -3.5489 -2.7192 -2.7192 -2.7078 -2.7078 -2.6775 -2.6775 3.8174 3.8174 10.7172 10.7172 10.8184 10.8184 12.3949 12.3949 13.1228 13.1228 14.4600 14.4600 14.5826 14.5826 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 ( 1596 PWs) bands (ev): -31.6319 -31.6319 -12.6056 -12.6056 -11.2076 -11.2076 -11.1207 -11.1207 -3.5548 -3.5548 -3.5360 -3.5360 -2.7165 -2.7165 -2.6872 -2.6872 -2.6793 -2.6793 4.5126 4.5126 9.3753 9.3753 11.0366 11.0366 12.1051 12.1051 12.6018 12.6018 14.3940 14.3940 14.6544 14.6544 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 ( 1595 PWs) bands (ev): -31.6277 -31.6277 -12.6258 -12.6258 -11.2510 -11.2510 -11.1152 -11.1152 -3.5545 -3.5545 -3.5265 -3.5265 -2.7168 -2.7168 -2.6902 -2.6902 -2.6565 -2.6565 5.4529 5.4529 7.8656 7.8656 11.3737 11.3737 11.7623 11.7623 12.2311 12.2311 14.4731 14.4731 15.5714 15.5714 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.5005 0.5005 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.2857 ( 1602 PWs) bands (ev): -31.6268 -31.6268 -12.6249 -12.6249 -11.2368 -11.2368 -11.1374 -11.1374 -3.5525 -3.5525 -3.5350 -3.5350 -2.7098 -2.7098 -2.7022 -2.7022 -2.6635 -2.6635 5.0781 5.0781 9.4819 9.4819 10.2346 10.2346 12.0907 12.0907 12.7785 12.7785 13.8952 13.8952 15.1503 15.1504 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 ( 1597 PWs) bands (ev): -31.6227 -31.6227 -12.6420 -12.6420 -11.2706 -11.2706 -11.1378 -11.1378 -3.5539 -3.5539 -3.5327 -3.5327 -2.7185 -2.7185 -2.7044 -2.7044 -2.6441 -2.6441 5.7479 5.7479 8.5865 8.5865 10.5216 10.5216 11.8543 11.8543 12.7461 12.7461 13.8723 13.8723 15.6129 15.6129 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 ( 1608 PWs) bands (ev): -31.6187 -31.6187 -12.6566 -12.6566 -11.2957 -11.2957 -11.1438 -11.1438 -3.5564 -3.5564 -3.5369 -3.5369 -2.7345 -2.7345 -2.7041 -2.7041 -2.6279 -2.6279 6.0387 6.0387 9.4389 9.4389 9.8912 9.8912 11.3966 11.3966 13.2928 13.2928 13.5222 13.5222 15.8031 15.8032 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.1570 0.1570 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.1429 ( 1602 PWs) bands (ev): -31.6348 -31.6348 -12.5914 -12.5914 -11.1533 -11.1533 -11.1474 -11.1474 -3.5499 -3.5499 -3.5492 -3.5492 -2.7082 -2.7082 -2.7068 -2.7068 -2.6824 -2.6824 4.0521 4.0521 10.7526 10.7526 10.8398 10.8398 11.7737 11.7737 13.2393 13.2393 13.3333 13.3333 15.8201 15.8201 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 ( 1592 PWs) bands (ev): -31.6296 -31.6296 -12.6133 -12.6133 -11.2087 -11.2087 -11.1414 -11.1414 -3.5510 -3.5510 -3.5365 -3.5365 -2.7084 -2.7084 -2.6918 -2.6918 -2.6789 -2.6789 4.7103 4.7103 9.5501 9.5501 10.9615 10.9615 11.9204 11.9204 12.3752 12.3752 13.7286 13.7286 15.2716 15.2716 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 ( 1594 PWs) bands (ev): -31.6255 -31.6255 -12.6324 -12.6324 -11.2524 -11.2524 -11.1344 -11.1344 -3.5523 -3.5523 -3.5277 -3.5277 -2.7124 -2.7124 -2.6953 -2.6953 -2.6544 -2.6544 5.5757 5.5757 8.1125 8.1125 11.1070 11.1070 11.8765 11.8765 12.3268 12.3268 13.9802 13.9802 15.6942 15.6942 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 ( 1593 PWs) bands (ev): -31.6245 -31.6245 -12.6321 -12.6321 -11.2322 -11.2322 -11.1641 -11.1641 -3.5479 -3.5479 -3.5343 -3.5343 -2.7089 -2.7089 -2.6947 -2.6947 -2.6641 -2.6641 5.2501 5.2501 9.7424 9.7424 10.1731 10.1731 11.9986 11.9986 12.3565 12.3565 13.4135 13.4135 15.6893 15.6893 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 ( 1601 PWs) bands (ev): -31.6205 -31.6205 -12.6487 -12.6487 -11.2672 -11.2672 -11.1634 -11.1634 -3.5496 -3.5496 -3.5302 -3.5302 -2.7209 -2.7209 -2.6921 -2.6921 -2.6431 -2.6431 5.8955 5.8955 8.7630 8.7630 10.3151 10.3151 12.0681 12.0681 12.6784 12.6784 13.5092 13.5092 15.3650 15.3650 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 ( 1600 PWs) bands (ev): -31.6165 -31.6165 -12.6632 -12.6632 -11.2917 -11.2917 -11.1713 -11.1713 -3.5491 -3.5491 -3.5326 -3.5326 -2.7342 -2.7342 -2.6883 -2.6883 -2.6265 -2.6265 6.2301 6.2301 9.4152 9.4152 9.6922 9.6922 11.6414 11.6414 13.2806 13.2806 13.7007 13.7007 14.4629 14.4629 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 ( 1598 PWs) bands (ev): -31.6195 -31.6195 -12.6495 -12.6495 -11.2235 -11.2235 -11.2197 -11.2197 -3.5419 -3.5419 -3.5304 -3.5304 -2.7190 -2.7190 -2.6657 -2.6657 -2.6634 -2.6634 5.7273 5.7273 9.9026 9.9026 10.0626 10.0626 11.2041 11.2041 12.9607 12.9607 13.0128 13.0128 14.8623 14.8624 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 ( 1588 PWs) bands (ev): -31.6155 -31.6155 -12.6650 -12.6650 -11.2649 -11.2649 -11.2131 -11.2131 -3.5409 -3.5409 -3.5253 -3.5253 -2.7278 -2.7278 -2.6645 -2.6645 -2.6380 -2.6380 6.3067 6.3067 9.0963 9.0963 9.8192 9.8192 12.1105 12.1105 12.8584 12.8584 13.2591 13.2591 14.1993 14.1993 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 ( 1599 PWs) bands (ev): -31.6115 -31.6115 -12.6795 -12.6795 -11.2828 -11.2828 -11.2307 -11.2307 -3.5338 -3.5338 -3.5234 -3.5234 -2.7338 -2.7338 -2.6512 -2.6512 -2.6230 -2.6230 6.7942 6.7942 9.1946 9.1946 9.3313 9.3313 11.8111 11.8111 12.7087 12.7087 13.9989 13.9989 14.1038 14.1038 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.4286 0.4286 0.4286 ( 1598 PWs) bands (ev): -31.6075 -31.6075 -12.6939 -12.6939 -11.2758 -11.2758 -11.2757 -11.2757 -3.5218 -3.5218 -3.5173 -3.5173 -2.7339 -2.7339 -2.6206 -2.6206 -2.6194 -2.6194 7.5941 7.5941 8.9234 8.9234 9.0667 9.0667 10.2901 10.2901 14.1536 14.1536 14.3592 14.3592 14.4330 14.4330 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 11.3737 ev ! total energy = -223.74838985 Ry Harris-Foulkes estimate = -223.74838985 Ry estimated scf accuracy < 1.4E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -53.52886589 Ry hartree contribution = 41.77280206 Ry xc contribution = -64.10789182 Ry ewald contribution = -147.88430357 Ry smearing contrib. (-TS) = -0.00013063 Ry convergence has been achieved in 7 iterations Writing output data file YIn.save init_run : 0.36s CPU 0.41s WALL ( 1 calls) electrons : 5.61s CPU 5.85s WALL ( 1 calls) Called by init_run: wfcinit : 0.22s CPU 0.24s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 4.39s CPU 4.53s WALL ( 8 calls) sum_band : 0.96s CPU 0.97s WALL ( 8 calls) v_of_rho : 0.02s CPU 0.02s WALL ( 8 calls) v_h : 0.00s CPU 0.00s WALL ( 8 calls) v_xc : 0.02s CPU 0.02s WALL ( 8 calls) newd : 0.24s CPU 0.27s WALL ( 8 calls) mix_rho : 0.00s CPU 0.01s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.01s WALL ( 340 calls) cegterg : 4.21s CPU 4.31s WALL ( 160 calls) Called by sum_band: sum_band:bec : 0.34s CPU 0.32s WALL ( 160 calls) addusdens : 0.14s CPU 0.14s WALL ( 8 calls) Called by *egterg: h_psi : 2.35s CPU 2.43s WALL ( 631 calls) s_psi : 0.12s CPU 0.12s WALL ( 631 calls) g_psi : 0.01s CPU 0.00s WALL ( 451 calls) cdiaghg : 1.60s CPU 1.64s WALL ( 591 calls) cegterg:over : 0.13s CPU 0.11s WALL ( 451 calls) cegterg:upda : 0.07s CPU 0.07s WALL ( 451 calls) cegterg:last : 0.03s CPU 0.04s WALL ( 160 calls) cdiaghg:chol : 0.07s CPU 0.09s WALL ( 591 calls) cdiaghg:inve : 0.01s CPU 0.02s WALL ( 591 calls) cdiaghg:para : 0.12s CPU 0.11s WALL ( 1182 calls) Called by h_psi: h_psi:vloc : 2.04s CPU 2.12s WALL ( 631 calls) h_psi:vnl : 0.30s CPU 0.31s WALL ( 631 calls) add_vuspsi : 0.18s CPU 0.17s WALL ( 631 calls) General routines calbec : 0.18s CPU 0.18s WALL ( 791 calls) fft : 0.05s CPU 0.04s WALL ( 242 calls) ffts : 0.01s CPU 0.00s WALL ( 64 calls) fftw : 2.26s CPU 2.36s WALL ( 65228 calls) interpolate : 0.01s CPU 0.01s WALL ( 64 calls) Parallel routines fft_scatter : 1.26s CPU 1.31s WALL ( 65534 calls) PWSCF : 8.16s CPU 9.12s WALL This run was terminated on: 19:50:42 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=