Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 19:44:48 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 Hg.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 6S 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 51 21 6 1659 435 76 Max 52 22 7 1664 449 79 Sum 1861 761 241 59813 15895 2777 bravais-lattice index = 14 lattice parameter (alat) = 7.4739 a.u. unit-cell volume = 417.4803 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 22.00 number of Kohn-Sham states= 30 kinetic-energy cutoff = 43.0000 Ry charge density cutoff = 416.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.473866 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 Hg read from file: /users/gautes/Pseudo/Hg.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: a7fd98bfb2640362c6e2ed76594f1960 Pseudo is Ultrasoft + core correction, Zval = 12.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1279 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 Sr read from file: /users/gautes/Pseudo/Sr.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 971034ad0a3914f6282d12395b0824f8 Pseudo is Ultrasoft + core correction, Zval = 10.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1221 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 Hg 12.00 200.59000 Hg( 1.00) Sr 10.00 87.62000 Sr( 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.0046296 k( 2) = ( 0.0000000 0.0000000 0.1666667), wk = 0.0277778 k( 3) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0277778 k( 4) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0138889 k( 5) = ( 0.0000000 0.1666667 0.1666667), wk = 0.0555556 k( 6) = ( 0.0000000 0.1666667 0.3333333), wk = 0.1111111 k( 7) = ( 0.0000000 0.1666667 -0.5000000), wk = 0.0555556 k( 8) = ( 0.0000000 0.3333333 0.3333333), wk = 0.0555556 k( 9) = ( 0.0000000 0.3333333 -0.5000000), wk = 0.0555556 k( 10) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0138889 k( 11) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0370370 k( 12) = ( 0.1666667 0.1666667 0.3333333), wk = 0.1111111 k( 13) = ( 0.1666667 0.1666667 -0.5000000), wk = 0.0555556 k( 14) = ( 0.1666667 0.3333333 0.3333333), wk = 0.1111111 k( 15) = ( 0.1666667 0.3333333 -0.5000000), wk = 0.1111111 k( 16) = ( 0.1666667 -0.5000000 -0.5000000), wk = 0.0277778 k( 17) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0370370 k( 18) = ( 0.3333333 0.3333333 -0.5000000), wk = 0.0555556 k( 19) = ( 0.3333333 -0.5000000 -0.5000000), wk = 0.0277778 k( 20) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0046296 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0046296 k( 2) = ( 0.0000000 0.0000000 0.1666667), wk = 0.0277778 k( 3) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0277778 k( 4) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0138889 k( 5) = ( 0.0000000 0.1666667 0.1666667), wk = 0.0555556 k( 6) = ( 0.0000000 0.1666667 0.3333333), wk = 0.1111111 k( 7) = ( 0.0000000 0.1666667 -0.5000000), wk = 0.0555556 k( 8) = ( 0.0000000 0.3333333 0.3333333), wk = 0.0555556 k( 9) = ( 0.0000000 0.3333333 -0.5000000), wk = 0.0555556 k( 10) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0138889 k( 11) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0370370 k( 12) = ( 0.1666667 0.1666667 0.3333333), wk = 0.1111111 k( 13) = ( 0.1666667 0.1666667 -0.5000000), wk = 0.0555556 k( 14) = ( 0.1666667 0.3333333 0.3333333), wk = 0.1111111 k( 15) = ( 0.1666667 0.3333333 -0.5000000), wk = 0.1111111 k( 16) = ( 0.1666667 -0.5000000 -0.5000000), wk = 0.0277778 k( 17) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0370370 k( 18) = ( 0.3333333 0.3333333 -0.5000000), wk = 0.0555556 k( 19) = ( 0.3333333 -0.5000000 -0.5000000), wk = 0.0277778 k( 20) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0046296 Dense grid: 59813 G-vectors FFT dimensions: ( 50, 50, 50) Smooth grid: 15895 G-vectors FFT dimensions: ( 32, 32, 32) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.05 Mb ( 120, 30) NL pseudopotentials 0.06 Mb ( 60, 68) Each V/rho on FFT grid 0.08 Mb ( 5000) Each G-vector array 0.01 Mb ( 1661) G-vector shells 0.00 Mb ( 383) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.22 Mb ( 120, 120) Each subspace H/S matrix 0.01 Mb ( 30, 30) Each matrix 0.06 Mb ( 68, 2, 30) Arrays for rho mixing 0.61 Mb ( 5000, 8) Initial potential from superposition of free atoms starting charge 21.99232, renormalised to 22.00000 Starting wfc are 34 randomized atomic wfcs total cpu time spent up to now is 2.5 secs per-process dynamical memory: 37.3 Mb Self-consistent Calculation iteration # 1 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 6.0 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.71E-04, avg # of iterations = 4.5 total cpu time spent up to now is 4.3 secs total energy = -169.94309793 Ry Harris-Foulkes estimate = -169.97730112 Ry estimated scf accuracy < 0.04692174 Ry iteration # 2 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.13E-04, avg # of iterations = 3.6 total cpu time spent up to now is 5.2 secs total energy = -169.95800535 Ry Harris-Foulkes estimate = -169.97240536 Ry estimated scf accuracy < 0.02419009 Ry iteration # 3 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.10E-04, avg # of iterations = 2.0 total cpu time spent up to now is 5.9 secs total energy = -169.96331842 Ry Harris-Foulkes estimate = -169.96332243 Ry estimated scf accuracy < 0.00024211 Ry iteration # 4 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 2 eigenvalues not converged c_bands: 1 eigenvalues not converged ethr = 1.10E-06, avg # of iterations = 6.9 total cpu time spent up to now is 7.1 secs total energy = -169.96372297 Ry Harris-Foulkes estimate = -169.96382085 Ry estimated scf accuracy < 0.00023747 Ry iteration # 5 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.08E-06, avg # of iterations = 1.4 total cpu time spent up to now is 7.7 secs total energy = -169.96370187 Ry Harris-Foulkes estimate = -169.96373812 Ry estimated scf accuracy < 0.00005022 Ry iteration # 6 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.28E-07, avg # of iterations = 4.3 total cpu time spent up to now is 8.6 secs total energy = -169.96372812 Ry Harris-Foulkes estimate = -169.96374304 Ry estimated scf accuracy < 0.00002698 Ry iteration # 7 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.23E-07, avg # of iterations = 2.4 total cpu time spent up to now is 9.3 secs total energy = -169.96373261 Ry Harris-Foulkes estimate = -169.96373274 Ry estimated scf accuracy < 0.00000021 Ry iteration # 8 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.69E-10, avg # of iterations = 5.0 total cpu time spent up to now is 10.4 secs total energy = -169.96373328 Ry Harris-Foulkes estimate = -169.96373338 Ry estimated scf accuracy < 0.00000019 Ry iteration # 9 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.85E-10, avg # of iterations = 2.5 total cpu time spent up to now is 11.0 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 1935 PWs) bands (ev): -26.4321 -26.4321 -9.1575 -9.1575 -8.0020 -8.0020 -8.0020 -8.0020 1.1948 1.1948 1.1948 1.1948 2.6087 2.6087 2.9316 2.9316 2.9316 2.9316 3.1588 3.1588 11.0051 11.0051 11.0051 11.0051 12.5323 12.5323 13.9242 13.9242 13.9242 13.9242 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.1667 ( 1983 PWs) bands (ev): -26.4286 -26.4286 -9.1752 -9.1752 -8.0482 -8.0482 -7.9973 -7.9973 1.2006 1.2006 1.2218 1.2218 2.8924 2.8924 2.9426 2.9426 3.0154 3.0154 3.1495 3.1495 10.2159 10.2159 11.0208 11.0208 12.2726 12.2726 12.9811 12.9811 13.0796 13.0796 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.3333 ( 1984 PWs) bands (ev): -26.4215 -26.4215 -9.2155 -9.2155 -8.1319 -8.1319 -7.9881 -7.9881 1.2119 1.2119 1.2668 1.2668 2.9627 2.9627 3.0045 3.0045 3.1343 3.1343 3.9009 3.9009 8.1894 8.1894 11.0546 11.0546 11.6616 11.6616 12.1091 12.1091 14.0902 14.0902 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000-0.5000 ( 1984 PWs) bands (ev): -26.4179 -26.4179 -9.2381 -9.2381 -8.1693 -8.1693 -7.9834 -7.9834 1.2174 1.2174 1.2854 1.2854 2.9711 2.9711 3.0388 3.0388 3.1291 3.1291 4.7770 4.7770 6.8595 6.8595 11.0727 11.0727 11.3840 11.3840 11.7957 11.7957 14.6941 14.6941 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.1667 ( 1991 PWs) bands (ev): -26.4251 -26.4251 -9.1906 -9.1906 -8.0683 -8.0683 -8.0160 -8.0160 1.2111 1.2111 1.2327 1.2327 2.9417 2.9417 2.9813 2.9813 3.1061 3.1061 3.2958 3.2958 10.1749 10.1749 10.5363 10.5363 11.7329 11.7329 12.5377 12.5377 14.0698 14.0699 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.3333 ( 1989 PWs) bands (ev): -26.4181 -26.4181 -9.2260 -9.2260 -8.1408 -8.1408 -8.0125 -8.0125 1.2161 1.2161 1.2617 1.2617 2.9588 2.9588 3.0274 3.0274 3.0992 3.0992 4.1127 4.1127 8.4357 8.4357 10.7197 10.7197 11.6034 11.6034 12.1716 12.1716 14.3023 14.3023 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667-0.5000 ( 1984 PWs) bands (ev): -26.4145 -26.4145 -9.2460 -9.2460 -8.1763 -8.1763 -8.0074 -8.0074 1.2172 1.2172 1.2736 1.2736 2.9624 2.9624 3.0441 3.0441 3.0994 3.0994 4.8321 4.8321 7.2978 7.2978 10.9800 10.9800 11.3229 11.3229 11.9702 11.9702 14.3433 14.3433 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333 0.3333 ( 2003 PWs) bands (ev): -26.4113 -26.4113 -9.2512 -9.2512 -8.1809 -8.1809 -8.0317 -8.0317 1.2096 1.2096 1.2547 1.2547 2.9395 2.9395 2.9835 2.9835 3.1444 3.1444 4.6478 4.6478 8.8119 8.8119 9.7113 9.7113 11.7405 11.7405 12.6492 12.6492 13.8243 13.8243 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333-0.5000 ( 1996 PWs) bands (ev): -26.4079 -26.4079 -9.2661 -9.2661 -8.2071 -8.2071 -8.0314 -8.0314 1.2023 1.2023 1.2505 1.2505 2.9245 2.9245 2.9535 2.9535 3.1800 3.1800 5.0401 5.0401 8.3810 8.3810 9.9220 9.9220 11.4779 11.4779 12.8405 12.8405 13.8153 13.8153 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000-0.5000-0.5000 ( 2000 PWs) bands (ev): -26.4045 -26.4045 -9.2783 -9.2783 -8.2285 -8.2285 -8.0337 -8.0337 1.1902 1.1902 1.2376 1.2376 2.9018 2.9018 2.9065 2.9065 3.2277 3.2277 5.2205 5.2205 9.0794 9.0794 9.5700 9.5700 10.9509 10.9509 13.5816 13.5816 13.6231 13.6231 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.0019 0.0019 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.1667 ( 1980 PWs) bands (ev): -26.4216 -26.4216 -9.2044 -9.2044 -8.0667 -8.0667 -8.0545 -8.0545 1.2309 1.2309 1.2311 1.2311 2.9835 2.9835 2.9984 2.9984 3.0862 3.0862 3.5440 3.5440 10.2519 10.2519 10.6678 10.6678 10.8389 10.8389 12.7267 12.7267 12.9237 12.9237 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.3333 ( 1992 PWs) bands (ev): -26.4147 -26.4147 -9.2362 -9.2362 -8.1415 -8.1415 -8.0459 -8.0459 1.2252 1.2252 1.2617 1.2617 2.9864 2.9864 3.0416 3.0416 3.0747 3.0747 4.2866 4.2866 8.7017 8.7017 10.5815 10.5815 11.5164 11.5164 11.9403 11.9403 13.5286 13.5286 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667-0.5000 ( 1974 PWs) bands (ev): -26.4112 -26.4112 -9.2545 -9.2545 -8.1776 -8.1776 -8.0389 -8.0389 1.2212 1.2212 1.2732 1.2732 2.9768 2.9768 3.0402 3.0402 3.0970 3.0970 4.9150 4.9150 7.6477 7.6477 10.6014 10.6014 11.6334 11.6334 12.0881 12.0881 13.7246 13.7246 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333 0.3333 ( 1990 PWs) bands (ev): -26.4079 -26.4079 -9.2608 -9.2608 -8.1729 -8.1729 -8.0755 -8.0755 1.2268 1.2268 1.2647 1.2647 2.9734 2.9734 3.0133 3.0133 3.1417 3.1417 4.7895 4.7895 9.0370 9.0370 9.6254 9.6254 11.6715 11.6715 12.1565 12.1565 13.3919 13.3919 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333-0.5000 ( 1996 PWs) bands (ev): -26.4045 -26.4045 -9.2752 -9.2752 -8.2005 -8.2005 -8.0747 -8.0747 1.2220 1.2220 1.2659 1.2659 2.9437 2.9437 3.0089 3.0089 3.1822 3.1822 5.1909 5.1909 8.4635 8.4635 9.6866 9.6866 11.8681 11.8681 12.4574 12.4574 13.4062 13.4062 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667-0.5000-0.5000 ( 2008 PWs) bands (ev): -26.4012 -26.4012 -9.2877 -9.2877 -8.2221 -8.2221 -8.0790 -8.0790 1.2230 1.2230 1.2558 1.2558 2.9168 2.9168 2.9916 2.9916 3.2339 3.2339 5.4231 5.4231 8.8355 8.8355 9.3255 9.3255 11.4903 11.4903 12.6391 12.6391 13.8392 13.8406 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 k = 0.3333 0.3333 0.3333 ( 1980 PWs) bands (ev): -26.4012 -26.4012 -9.2837 -9.2837 -8.1587 -8.1587 -8.1546 -8.1546 1.2492 1.2492 1.2822 1.2822 2.9533 2.9533 3.1362 3.1362 3.1499 3.1499 5.2533 5.2533 8.9873 8.9873 9.5113 9.5113 10.6513 10.6513 13.1318 13.1318 13.1427 13.1428 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.1242 0.1242 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333-0.5000 ( 1992 PWs) bands (ev): -26.3978 -26.3978 -9.2970 -9.2970 -8.1938 -8.1938 -8.1485 -8.1485 1.2567 1.2567 1.2882 1.2882 2.9382 2.9382 3.1444 3.1444 3.2031 3.2031 5.6701 5.6701 8.5177 8.5177 9.2071 9.2071 10.8238 10.8238 13.0423 13.0423 13.1157 13.1157 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 k = 0.3333-0.5000-0.5000 ( 1964 PWs) bands (ev): -26.3945 -26.3945 -9.3097 -9.3097 -8.2092 -8.2092 -8.1635 -8.1635 1.2791 1.2791 1.2884 1.2884 2.9323 2.9323 3.1730 3.1730 3.2507 3.2507 6.1085 6.1085 8.3866 8.3866 8.8949 8.8949 10.0779 10.0779 13.0178 13.0178 14.3824 14.4086 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 k =-0.5000-0.5000-0.5000 ( 2008 PWs) bands (ev): -26.3911 -26.3911 -9.3224 -9.3224 -8.2027 -8.2027 -8.2027 -8.2027 1.3030 1.3030 1.3030 1.3030 2.9376 2.9376 3.2610 3.2610 3.2610 3.2610 6.9681 6.9681 8.0257 8.0257 8.7013 8.7013 8.7013 8.7013 14.6506 14.6506 14.7504 14.7505 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 the Fermi energy is 9.4847 ev ! total energy = -169.96373329 Ry Harris-Foulkes estimate = -169.96373330 Ry estimated scf accuracy < 8.1E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -29.36200817 Ry hartree contribution = 28.11176802 Ry xc contribution = -50.33227048 Ry ewald contribution = -118.38119440 Ry smearing contrib. (-TS) = -0.00002826 Ry convergence has been achieved in 9 iterations Writing output data file SrHg.save init_run : 0.37s CPU 0.43s WALL ( 1 calls) electrons : 8.30s CPU 8.59s WALL ( 1 calls) Called by init_run: wfcinit : 0.17s CPU 0.19s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 6.36s CPU 6.58s WALL ( 10 calls) sum_band : 1.38s CPU 1.42s WALL ( 10 calls) v_of_rho : 0.04s CPU 0.04s WALL ( 10 calls) v_h : 0.00s CPU 0.00s WALL ( 10 calls) v_xc : 0.03s CPU 0.03s WALL ( 10 calls) newd : 0.52s CPU 0.54s WALL ( 10 calls) mix_rho : 0.02s CPU 0.02s WALL ( 10 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.02s WALL ( 420 calls) cegterg : 6.11s CPU 6.28s WALL ( 200 calls) Called by sum_band: sum_band:bec : 0.38s CPU 0.40s WALL ( 200 calls) addusdens : 0.42s CPU 0.43s WALL ( 10 calls) Called by *egterg: h_psi : 2.79s CPU 2.98s WALL ( 992 calls) s_psi : 0.15s CPU 0.17s WALL ( 992 calls) g_psi : 0.02s CPU 0.01s WALL ( 772 calls) cdiaghg : 2.79s CPU 2.76s WALL ( 952 calls) cegterg:over : 0.19s CPU 0.18s WALL ( 772 calls) cegterg:upda : 0.15s CPU 0.14s WALL ( 772 calls) cegterg:last : 0.06s CPU 0.06s WALL ( 213 calls) cdiaghg:chol : 0.12s CPU 0.14s WALL ( 952 calls) cdiaghg:inve : 0.04s CPU 0.05s WALL ( 952 calls) cdiaghg:para : 0.19s CPU 0.18s WALL ( 1904 calls) Called by h_psi: h_psi:vloc : 2.43s CPU 2.58s WALL ( 992 calls) h_psi:vnl : 0.34s CPU 0.39s WALL ( 992 calls) add_vuspsi : 0.17s CPU 0.22s WALL ( 992 calls) General routines calbec : 0.23s CPU 0.21s WALL ( 1192 calls) fft : 0.10s CPU 0.08s WALL ( 304 calls) ffts : 0.00s CPU 0.01s WALL ( 80 calls) fftw : 2.58s CPU 2.72s WALL ( 83520 calls) interpolate : 0.03s CPU 0.03s WALL ( 80 calls) Parallel routines fft_scatter : 1.74s CPU 1.75s WALL ( 83904 calls) PWSCF : 10.90s CPU 11.83s WALL This run was terminated on: 19:45: 0 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=