Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 19:45:28 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 101 48 13 4613 1533 225 Max 102 49 14 4616 1553 230 Sum 3657 1757 489 166101 55467 8217 bravais-lattice index = 14 lattice parameter (alat) = 10.5037 a.u. unit-cell volume = 1158.8373 (a.u.)^3 number of atoms/cell = 8 number of atomic types = 2 number of electrons = 48.00 number of Kohn-Sham states= 58 kinetic-energy cutoff = 50.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)= 10.503664 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 Zr read from file: /users/gautes/Pseudo/Zr.rel-pbe-n-nc.UPF MD5 check sum: 453101d5cee66cf3691e2371f73c5d7d Pseudo is Norm-conserving + core correction, Zval = 4.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1225 points, 4 beta functions with: l(1) = 1 l(2) = 1 l(3) = 2 l(4) = 2 atomic species valence mass pseudopotential Hg 12.00 200.59000 Hg( 1.00) Zr 4.00 91.22400 Zr( 1.00) 48 Sym. Ops., with inversion, found (24 have fractional translation) 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 ) f =( -0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) cart. s( 5) = ( 0.0000000 1.0000000 -0.0000000 ) f =( -0.5000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.5000000 ) isym = 6 180 deg rotation - cart. axis [1,-1,0] cryst. s( 6) = ( 0 -1 0 ) f =( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) cart. s( 6) = ( -0.0000000 -1.0000000 0.0000000 ) f =( 0.5000000 ) ( -1.0000000 0.0000000 -0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.5000000 ) isym = 7 90 deg rotation - cart. axis [0,0,-1] cryst. s( 7) = ( 0 -1 0 ) f =( -0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) cart. s( 7) = ( -0.0000000 1.0000000 0.0000000 ) f =( -0.5000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.5000000 ) isym = 8 90 deg rotation - cart. axis [0,0,1] cryst. s( 8) = ( 0 1 0 ) f =( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) cart. s( 8) = ( 0.0000000 -1.0000000 0.0000000 ) f =( 0.5000000 ) ( 1.0000000 -0.0000000 0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.5000000 ) isym = 9 180 deg rotation - cart. axis [1,0,1] cryst. s( 9) = ( 0 0 1 ) f =( -0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) cart. s( 9) = ( 0.0000000 -0.0000000 1.0000000 ) f =( -0.5000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.5000000 ) isym = 10 180 deg rotation - cart. axis [-1,0,1] cryst. s(10) = ( 0 0 -1 ) f =( 0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) cart. s(10) = ( -0.0000000 -0.0000000 -1.0000000 ) f =( 0.5000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.5000000 ) isym = 11 90 deg rotation - cart. axis [0,1,0] cryst. s(11) = ( 0 0 -1 ) f =( -0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) cart. s(11) = ( -0.0000000 0.0000000 1.0000000 ) f =( -0.5000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.5000000 ) isym = 12 90 deg rotation - cart. axis [0,-1,0] cryst. s(12) = ( 0 0 1 ) f =( 0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) cart. s(12) = ( 0.0000000 0.0000000 -1.0000000 ) f =( 0.5000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.5000000 ) isym = 13 180 deg rotation - cart. axis [0,1,1] cryst. s(13) = ( -1 0 0 ) f =( 0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) cart. s(13) = ( -1.0000000 0.0000000 0.0000000 ) f =( 0.5000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.5000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) isym = 14 180 deg rotation - cart. axis [0,1,-1] cryst. s(14) = ( -1 0 0 ) f =( 0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) cart. s(14) = ( -1.0000000 0.0000000 0.0000000 ) f =( 0.5000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.5000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) isym = 15 90 deg rotation - cart. axis [-1,0,0] cryst. s(15) = ( 1 0 0 ) f =( -0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) cart. s(15) = ( 1.0000000 -0.0000000 0.0000000 ) f =( -0.5000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.5000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) isym = 16 90 deg rotation - cart. axis [1,0,0] cryst. s(16) = ( 1 0 0 ) f =( -0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) cart. s(16) = ( 1.0000000 0.0000000 -0.0000000 ) f =( -0.5000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.5000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) 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 ) f =( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) cart. s(29) = ( -0.0000000 -1.0000000 0.0000000 ) f =( 0.5000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.5000000 ) isym = 30 inv. 180 deg rotation - cart. axis [1,-1,0] cryst. s(30) = ( 0 1 0 ) f =( -0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) cart. s(30) = ( 0.0000000 1.0000000 -0.0000000 ) f =( -0.5000000 ) ( 1.0000000 -0.0000000 0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.5000000 ) isym = 31 inv. 90 deg rotation - cart. axis [0,0,-1] cryst. s(31) = ( 0 1 0 ) f =( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) cart. s(31) = ( 0.0000000 -1.0000000 -0.0000000 ) f =( 0.5000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.5000000 ) isym = 32 inv. 90 deg rotation - cart. axis [0,0,1] cryst. s(32) = ( 0 -1 0 ) f =( -0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) cart. s(32) = ( -0.0000000 1.0000000 -0.0000000 ) f =( -0.5000000 ) ( -1.0000000 0.0000000 -0.0000000 ) ( 0.5000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.5000000 ) isym = 33 inv. 180 deg rotation - cart. axis [1,0,1] cryst. s(33) = ( 0 0 -1 ) f =( 0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) cart. s(33) = ( -0.0000000 0.0000000 -1.0000000 ) f =( 0.5000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.5000000 ) isym = 34 inv. 180 deg rotation - cart. axis [-1,0,1] cryst. s(34) = ( 0 0 1 ) f =( -0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) cart. s(34) = ( 0.0000000 0.0000000 1.0000000 ) f =( -0.5000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.5000000 ) isym = 35 inv. 90 deg rotation - cart. axis [0,1,0] cryst. s(35) = ( 0 0 1 ) f =( 0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) ( -1 0 0 ) ( 0.5000000 ) cart. s(35) = ( 0.0000000 -0.0000000 -1.0000000 ) f =( 0.5000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.5000000 ) isym = 36 inv. 90 deg rotation - cart. axis [0,-1,0] cryst. s(36) = ( 0 0 -1 ) f =( -0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) ( 1 0 0 ) ( 0.5000000 ) cart. s(36) = ( -0.0000000 -0.0000000 1.0000000 ) f =( -0.5000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.5000000 ) isym = 37 inv. 180 deg rotation - cart. axis [0,1,1] cryst. s(37) = ( 1 0 0 ) f =( -0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) cart. s(37) = ( 1.0000000 -0.0000000 -0.0000000 ) f =( -0.5000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.5000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) isym = 38 inv. 180 deg rotation - cart. axis [0,1,-1] cryst. s(38) = ( 1 0 0 ) f =( -0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) cart. s(38) = ( 1.0000000 0.0000000 0.0000000 ) f =( -0.5000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.5000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) isym = 39 inv. 90 deg rotation - cart. axis [-1,0,0] cryst. s(39) = ( -1 0 0 ) f =( 0.5000000 ) ( 0 0 1 ) ( 0.5000000 ) ( 0 -1 0 ) ( 0.5000000 ) cart. s(39) = ( -1.0000000 0.0000000 -0.0000000 ) f =( 0.5000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.5000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.5000000 ) isym = 40 inv. 90 deg rotation - cart. axis [1,0,0] cryst. s(40) = ( -1 0 0 ) f =( 0.5000000 ) ( 0 0 -1 ) ( 0.5000000 ) ( 0 1 0 ) ( 0.5000000 ) cart. s(40) = ( -1.0000000 0.0000000 0.0000000 ) f =( 0.5000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.5000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.5000000 ) isym = 41 inv. 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(41) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(41) = ( -0.0000000 -1.0000000 0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 42 inv. 120 deg rotation - cart. axis [-1,1,1] cryst. s(42) = ( 0 0 1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(42) = ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 43 inv. 120 deg rotation - cart. axis [1,1,-1] cryst. s(43) = ( 0 0 1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(43) = ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 44 inv. 120 deg rotation - cart. axis [1,-1,1] cryst. s(44) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(44) = ( -0.0000000 1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 45 inv. 120 deg rotation - cart. axis [1,1,1] cryst. s(45) = ( 0 -1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(45) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 46 inv. 120 deg rotation - cart. axis [-1,1,-1] cryst. s(46) = ( 0 1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(46) = ( 0.0000000 0.0000000 -1.0000000 ) ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 47 inv. 120 deg rotation - cart. axis [1,-1,-1] cryst. s(47) = ( 0 1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(47) = ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 48 inv. 120 deg rotation - cart. axis [-1,-1,1] cryst. s(48) = ( 0 -1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(48) = ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) double point group O_h (m-3m) there are 16 classes and 6 irreducible representations the character table: E -E 8C3 -8C3 3C2 6C4 -6C4 6C2' i -i 8S6 -8S6 -3C2 -6C2' G_6+ 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 2.00 -2.00 1.00 -1.00 G_7+ 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 2.00 -2.00 1.00 -1.00 G_8+ 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 4.00 -4.00 -1.00 1.00 G_6- 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 -2.00 2.00 -1.00 1.00 G_7- 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 -2.00 2.00 -1.00 1.00 G_8- 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 -4.00 4.00 1.00 -1.00 3s_h 6S4 -6S4 6s_d -3s_h -6s_d G_6+ 0.00 1.41 -1.41 0.00 G_7+ 0.00 -1.41 1.41 0.00 G_8+ 0.00 0.00 0.00 0.00 G_6- 0.00 -1.41 1.41 0.00 G_7- 0.00 1.41 -1.41 0.00 G_8- 0.00 0.00 0.00 0.00 the symmetry operations in each class and the name of the first element: E 1 identity -E -1 identity E 8C3 17 19 20 18 24 21 22 23 120 deg rotation - cart. axis [-1,-1,-1] -8C3 -17 -19 -20 -18 -24 -21 -22 -23 120 deg rotation - cart. axis [-1,-1,-1] E 3C2 -3C2 2 -2 4 -4 3 -3 180 deg rotation - cart. axis [0,0,1] 6C4 7 8 15 16 12 11 90 deg rotation - cart. axis [0,0,-1] -6C4 -7 -8 -15 -16 -12 -11 90 deg rotation - cart. axis [0,0,-1] E 6C2'-6C2' 5 -5 6 -6 14 -13 -14 13 -10 -9 10 9 180 deg rotation - cart. axis [1,1,0] i 25 inversion -i -25 inversion E 8S6 41 43 44 42 48 45 46 47 inv. 120 deg rotation - cart. axis [-1,-1,-1] -8S6 -41 -43 -44 -42 -48 -45 -46 -47 inv. 120 deg rotation - cart. axis [-1,-1,-1] E 3s_h-3s_h 26 -26 28 -28 27 -27 inv. 180 deg rotation - cart. axis [0,0,1] 6S4 31 32 39 40 36 35 inv. 90 deg rotation - cart. axis [0,0,-1] -6S4 -31 -32 -39 -40 -36 -35 inv. 90 deg rotation - cart. axis [0,0,-1] E 6s_d-6s_d 29 -29 30 -30 38 -37 -38 37 -34 -33 34 33 inv. 180 deg rotation - cart. axis [1,1,0] Cartesian axes number of k points= 10 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 Dense grid: 166101 G-vectors FFT dimensions: ( 72, 72, 72) Smooth grid: 55467 G-vectors FFT dimensions: ( 48, 48, 48) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.34 Mb ( 388, 58) NL pseudopotentials 0.49 Mb ( 194, 164) Each V/rho on FFT grid 0.16 Mb ( 10368) Each G-vector array 0.04 Mb ( 4616) G-vector shells 0.01 Mb ( 818) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 1.37 Mb ( 388, 232) Each subspace H/S matrix 0.05 Mb ( 58, 58) Each matrix 0.29 Mb ( 164, 2, 58) Arrays for rho mixing 1.27 Mb ( 10368, 8) Initial potential from superposition of free atoms starting charge 47.98686, renormalised to 48.00000 Starting wfc are 144 randomized atomic wfcs total cpu time spent up to now is 2.5 secs per-process dynamical memory: 33.6 Mb Self-consistent Calculation iteration # 1 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 2.0 total cpu time spent up to now is 4.1 secs total energy = -257.84847873 Ry Harris-Foulkes estimate = -258.47075876 Ry estimated scf accuracy < 0.77997807 Ry iteration # 2 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.62E-03, avg # of iterations = 3.8 total cpu time spent up to now is 6.3 secs total energy = -258.11015207 Ry Harris-Foulkes estimate = -258.58368684 Ry estimated scf accuracy < 0.89835128 Ry iteration # 3 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.62E-03, avg # of iterations = 2.0 total cpu time spent up to now is 8.0 secs total energy = -258.26742145 Ry Harris-Foulkes estimate = -258.27688645 Ry estimated scf accuracy < 0.01650003 Ry iteration # 4 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.44E-05, avg # of iterations = 4.4 total cpu time spent up to now is 10.8 secs total energy = -258.28942507 Ry Harris-Foulkes estimate = -258.29448580 Ry estimated scf accuracy < 0.01134381 Ry iteration # 5 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.36E-05, avg # of iterations = 2.0 total cpu time spent up to now is 12.3 secs total energy = -258.29011006 Ry Harris-Foulkes estimate = -258.29071542 Ry estimated scf accuracy < 0.00104205 Ry iteration # 6 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 2 eigenvalues not converged c_bands: 1 eigenvalues not converged ethr = 2.17E-06, avg # of iterations = 8.6 total cpu time spent up to now is 15.3 secs total energy = -258.29089256 Ry Harris-Foulkes estimate = -258.29109099 Ry estimated scf accuracy < 0.00038084 Ry iteration # 7 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.93E-07, avg # of iterations = 2.6 total cpu time spent up to now is 17.1 secs total energy = -258.29097886 Ry Harris-Foulkes estimate = -258.29097699 Ry estimated scf accuracy < 0.00000150 Ry iteration # 8 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.12E-09, avg # of iterations = 9.9 total cpu time spent up to now is 21.2 secs total energy = -258.29098736 Ry Harris-Foulkes estimate = -258.29099042 Ry estimated scf accuracy < 0.00000593 Ry iteration # 9 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.12E-09, avg # of iterations = 4.3 total cpu time spent up to now is 23.7 secs total energy = -258.29098877 Ry Harris-Foulkes estimate = -258.29098877 Ry estimated scf accuracy < 0.00000007 Ry iteration # 10 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.39E-10, avg # of iterations = 5.6 total cpu time spent up to now is 26.3 secs total energy = -258.29098879 Ry Harris-Foulkes estimate = -258.29098880 Ry estimated scf accuracy < 0.00000007 Ry iteration # 11 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.39E-10, avg # of iterations = 1.0 total cpu time spent up to now is 27.6 secs total energy = -258.29098875 Ry Harris-Foulkes estimate = -258.29098879 Ry estimated scf accuracy < 0.00000006 Ry iteration # 12 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.17E-10, avg # of iterations = 4.0 total cpu time spent up to now is 29.5 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 6931 PWs) bands (ev): 1.7554 1.7554 1.7554 1.7554 1.8102 1.8102 1.8102 1.8102 2.6720 2.6720 3.4705 3.4705 3.4705 3.4705 3.5034 3.5034 3.5578 3.5578 3.5578 3.5578 3.6980 3.6980 4.6278 4.6278 7.7300 7.7300 8.0185 8.0185 8.0185 8.0185 8.2299 8.2299 8.6206 8.6206 8.6206 8.6206 8.6963 8.6963 8.6963 8.6963 8.9054 8.9054 9.0072 9.0072 9.1797 9.1797 9.1797 9.1797 9.5549 9.5549 10.7235 10.7235 10.7236 10.7237 10.8930 10.8941 10.8943 10.8943 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.2000 ( 6920 PWs) bands (ev): 1.7512 1.7512 1.7575 1.7575 1.8004 1.8004 1.8032 1.8032 2.9005 2.9005 3.4524 3.4524 3.4865 3.4865 3.5071 3.5071 3.5341 3.5341 3.5517 3.5517 3.6648 3.6648 4.5226 4.5226 7.7568 7.7568 7.7713 7.7713 8.0400 8.0400 8.1362 8.1362 8.4872 8.4872 8.5439 8.5439 8.7546 8.7546 9.1166 9.1166 9.1487 9.1487 9.1953 9.1953 9.2437 9.2437 9.2538 9.2538 9.7103 9.7103 10.6933 10.6933 10.7959 10.7959 10.8192 10.8192 10.9357 10.9357 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4000 ( 6920 PWs) bands (ev): 1.7401 1.7401 1.7452 1.7452 1.7967 1.7967 1.7987 1.7987 3.4020 3.4020 3.4199 3.4199 3.4956 3.4956 3.5345 3.5345 3.5486 3.5486 3.5845 3.5845 3.6113 3.6113 4.1693 4.1693 7.1692 7.1692 7.2039 7.2039 8.1522 8.1522 8.2641 8.2641 8.5312 8.5312 8.6592 8.6592 9.2619 9.2619 9.3020 9.3020 9.5236 9.5236 9.6402 9.6402 9.7086 9.7086 10.0770 10.0770 10.1131 10.1131 10.4853 10.4853 10.5703 10.5703 10.5804 10.5804 10.8803 10.8803 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.6669 0.6669 0.1235 0.1235 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.2000 ( 6915 PWs) bands (ev): 1.7497 1.7497 1.7567 1.7567 1.7908 1.7908 1.7950 1.7950 3.1147 3.1147 3.4487 3.4487 3.4834 3.4834 3.5062 3.5062 3.5236 3.5236 3.5455 3.5455 3.6344 3.6344 4.4868 4.4868 7.4704 7.4704 7.6482 7.6482 8.0916 8.0916 8.2082 8.2082 8.5480 8.5480 8.7617 8.7617 8.7798 8.7798 9.0342 9.0342 9.1620 9.1620 9.3213 9.3213 9.3832 9.3832 9.6267 9.6267 9.8319 9.8319 10.7454 10.7454 10.7982 10.7983 10.8314 10.8315 10.9150 10.9150 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.4000 ( 6921 PWs) bands (ev): 1.7376 1.7376 1.7430 1.7430 1.7875 1.7875 1.7909 1.7909 3.3962 3.3962 3.4529 3.4529 3.4890 3.4890 3.5186 3.5186 3.5516 3.5516 3.5838 3.5838 3.7327 3.7327 4.2666 4.2666 7.0169 7.0169 7.0834 7.0834 8.1349 8.1349 8.4110 8.4110 8.7086 8.7086 8.7863 8.7863 9.1671 9.1671 9.4412 9.4412 9.4952 9.4952 9.5570 9.5570 9.7619 9.7619 10.0436 10.0436 10.3423 10.3423 10.5541 10.5541 10.6356 10.6356 10.7439 10.7439 11.0319 11.0319 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9587 0.9587 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4000 0.4000 ( 6922 PWs) bands (ev): 1.7251 1.7251 1.7297 1.7297 1.7818 1.7818 1.7868 1.7868 3.3632 3.3632 3.4477 3.4477 3.4765 3.4765 3.5119 3.5119 3.5376 3.5376 3.5944 3.5944 4.1402 4.1402 4.3858 4.3858 6.5924 6.5924 6.7219 6.7219 8.4746 8.4746 8.8167 8.8167 8.8966 8.8966 8.9950 8.9950 9.1760 9.1760 9.3251 9.3251 9.4245 9.4245 9.5699 9.5699 9.7955 9.7955 9.8685 9.8685 10.9379 10.9379 11.0097 11.0097 11.0374 11.0374 11.0580 11.0580 11.2976 11.2976 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 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.2000 0.2000 0.2000 ( 6929 PWs) bands (ev): 1.7526 1.7526 1.7546 1.7546 1.7808 1.7808 1.7834 1.7834 3.3045 3.3045 3.4586 3.4586 3.4625 3.4625 3.4928 3.4928 3.5322 3.5322 3.5524 3.5524 3.6031 3.6031 4.4853 4.4853 7.5062 7.5062 7.6114 7.6114 7.7840 7.7840 8.3465 8.3465 8.5345 8.5345 8.8234 8.8234 8.8838 8.8838 8.9528 8.9528 9.2875 9.2875 9.5204 9.5204 9.6627 9.6627 9.7234 9.7234 9.9579 9.9579 10.8334 10.8334 10.8468 10.8468 10.8779 10.8779 10.9677 10.9677 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9999 0.9999 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.4000 ( 6910 PWs) bands (ev): 1.7396 1.7396 1.7427 1.7427 1.7759 1.7759 1.7792 1.7792 3.4096 3.4096 3.4521 3.4521 3.4808 3.4808 3.5127 3.5127 3.5343 3.5343 3.5502 3.5502 3.8970 3.8970 4.3544 4.3544 7.0926 7.0926 7.1076 7.1076 8.0469 8.0469 8.3291 8.3291 8.7663 8.7663 8.9389 8.9389 9.1566 9.1566 9.2760 9.2760 9.6447 9.6447 9.6961 9.6961 9.7584 9.7584 10.3287 10.3287 10.3829 10.3829 10.4452 10.4452 10.7288 10.7288 10.8337 10.8337 11.2659 11.2670 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.4000 0.4000 ( 6919 PWs) bands (ev): 1.7298 1.7298 1.7340 1.7340 1.7661 1.7661 1.7722 1.7722 3.3906 3.3906 3.4434 3.4434 3.4674 3.4674 3.5066 3.5066 3.5123 3.5123 3.5423 3.5423 4.2706 4.2706 4.4685 4.4685 6.8731 6.8731 6.9688 6.9688 8.2052 8.2052 8.5141 8.5141 8.8022 8.8022 9.1312 9.1312 9.3042 9.3042 9.3470 9.3470 9.4627 9.4627 9.7309 9.7309 9.9289 9.9289 9.9743 9.9743 10.8720 10.8720 10.9067 10.9067 10.9908 10.9908 11.1936 11.1936 11.2820 11.2820 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9997 0.9997 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4000 0.4000 0.4000 ( 6923 PWs) bands (ev): 1.7353 1.7353 1.7361 1.7361 1.7482 1.7482 1.7483 1.7483 3.4334 3.4334 3.4391 3.4391 3.4414 3.4414 3.4691 3.4691 3.4873 3.4873 3.4961 3.4961 4.5303 4.5303 4.6075 4.6075 7.3006 7.3006 7.5569 7.5569 7.7281 7.7281 8.3063 8.3063 8.4982 8.4982 8.8799 8.8799 9.4705 9.4705 9.5579 9.5579 9.5787 9.5787 9.9643 9.9643 10.0314 10.0314 10.0524 10.0524 10.9097 10.9097 10.9123 10.9123 11.2561 11.2562 11.3358 11.3360 11.3727 11.3846 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9999 0.9999 0.9828 0.9828 0.9244 0.9244 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 10.0864 ev ! total energy = -258.29098876 Ry Harris-Foulkes estimate = -258.29098877 Ry estimated scf accuracy < 9.4E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -19.18857225 Ry hartree contribution = 48.77807866 Ry xc contribution = -103.07225978 Ry ewald contribution = -184.80802563 Ry smearing contrib. (-TS) = -0.00020975 Ry convergence has been achieved in 12 iterations Writing output data file Zr3Hg.save init_run : 1.26s CPU 1.36s WALL ( 1 calls) electrons : 25.69s CPU 27.05s WALL ( 1 calls) Called by init_run: wfcinit : 1.00s CPU 1.04s WALL ( 1 calls) potinit : 0.03s CPU 0.03s WALL ( 1 calls) Called by electrons: c_bands : 21.85s CPU 22.15s WALL ( 12 calls) sum_band : 2.88s CPU 3.36s WALL ( 12 calls) v_of_rho : 0.10s CPU 0.09s WALL ( 13 calls) v_h : 0.00s CPU 0.01s WALL ( 13 calls) v_xc : 0.09s CPU 0.08s WALL ( 13 calls) newd : 0.78s CPU 1.37s WALL ( 13 calls) mix_rho : 0.07s CPU 0.07s WALL ( 12 calls) Called by c_bands: init_us_2 : 0.08s CPU 0.07s WALL ( 250 calls) cegterg : 21.06s CPU 21.30s WALL ( 120 calls) Called by sum_band: sum_band:bec : 0.24s CPU 0.28s WALL ( 120 calls) addusdens : 0.45s CPU 0.91s WALL ( 12 calls) Called by *egterg: h_psi : 13.72s CPU 13.93s WALL ( 632 calls) s_psi : 0.67s CPU 0.65s WALL ( 632 calls) g_psi : 0.03s CPU 0.04s WALL ( 502 calls) cdiaghg : 4.75s CPU 4.83s WALL ( 622 calls) cegterg:over : 0.93s CPU 0.96s WALL ( 502 calls) cegterg:upda : 0.74s CPU 0.79s WALL ( 502 calls) cegterg:last : 0.37s CPU 0.35s WALL ( 151 calls) cdiaghg:chol : 0.26s CPU 0.29s WALL ( 622 calls) cdiaghg:inve : 0.16s CPU 0.19s WALL ( 622 calls) cdiaghg:para : 0.37s CPU 0.32s WALL ( 1244 calls) Called by h_psi: h_psi:vloc : 11.76s CPU 11.95s WALL ( 632 calls) h_psi:vnl : 1.89s CPU 1.92s WALL ( 632 calls) add_vuspsi : 0.94s CPU 0.91s WALL ( 632 calls) General routines calbec : 1.22s CPU 1.27s WALL ( 752 calls) fft : 0.20s CPU 0.20s WALL ( 387 calls) ffts : 0.02s CPU 0.02s WALL ( 100 calls) fftw : 12.41s CPU 12.51s WALL ( 116552 calls) interpolate : 0.07s CPU 0.07s WALL ( 100 calls) Parallel routines fft_scatter : 4.61s CPU 4.64s WALL ( 117039 calls) PWSCF : 28.73s CPU 30.71s WALL This run was terminated on: 19:45:58 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=