Program PWSCF v.5.3.0 (svn rev. 11974) starts on 3Jan2017 at 13:57:34 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 Ba.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 5S 6S renormalized 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 56 23 7 1878 489 83 Max 57 24 8 1884 509 87 Sum 2017 829 253 67715 17941 3071 bravais-lattice index = 14 lattice parameter (alat) = 7.7951 a.u. unit-cell volume = 473.6618 (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.795120 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 Ba read from file: /users/gautes/Pseudo/Ba.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 3006f3397849ee49953ff32be8e81b23 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 1251 points, 11 beta functions with: l(1) = 0 l(2) = 0 l(3) = 0 l(4) = 1 l(5) = 1 l(6) = 1 l(7) = 1 l(8) = 2 l(9) = 2 l(10) = 2 l(11) = 2 Q(r) pseudized with 0 coefficients PseudoPot. # 2 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 atomic species valence mass pseudopotential Ba 10.00 137.32700 Ba( 1.00) Hg 12.00 200.59000 Hg( 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: 67715 G-vectors FFT dimensions: ( 54, 54, 54) Smooth grid: 17941 G-vectors FFT dimensions: ( 36, 36, 36) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.06 Mb ( 128, 30) NL pseudopotentials 0.07 Mb ( 64, 69) Each V/rho on FFT grid 0.09 Mb ( 5832) Each G-vector array 0.01 Mb ( 1884) G-vector shells 0.00 Mb ( 416) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.23 Mb ( 128, 120) Each subspace H/S matrix 0.01 Mb ( 30, 30) Each matrix 0.06 Mb ( 69, 2, 30) Arrays for rho mixing 0.71 Mb ( 5832, 8) Initial potential from superposition of free atoms starting charge 21.98288, renormalised to 22.00000 Starting wfc are 28 randomized atomic wfcs + 2 random wfc total cpu time spent up to now is 2.6 secs per-process dynamical memory: 42.9 Mb Self-consistent Calculation iteration # 1 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 5.2 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 2.30E-04, avg # of iterations = 3.4 total cpu time spent up to now is 4.5 secs total energy = -162.52327457 Ry Harris-Foulkes estimate = -162.56388910 Ry estimated scf accuracy < 0.05310416 Ry iteration # 2 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.41E-04, avg # of iterations = 3.3 total cpu time spent up to now is 5.4 secs total energy = -162.54084243 Ry Harris-Foulkes estimate = -162.55695995 Ry estimated scf accuracy < 0.02676878 Ry iteration # 3 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.22E-04, avg # of iterations = 2.3 total cpu time spent up to now is 6.2 secs total energy = -162.54691033 Ry Harris-Foulkes estimate = -162.54692915 Ry estimated scf accuracy < 0.00023755 Ry iteration # 4 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 1 eigenvalues not converged ethr = 1.08E-06, avg # of iterations = 6.0 total cpu time spent up to now is 7.6 secs total energy = -162.54727088 Ry Harris-Foulkes estimate = -162.54733834 Ry estimated scf accuracy < 0.00015576 Ry iteration # 5 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.08E-07, avg # of iterations = 2.0 total cpu time spent up to now is 8.3 secs total energy = -162.54727604 Ry Harris-Foulkes estimate = -162.54728657 Ry estimated scf accuracy < 0.00001659 Ry iteration # 6 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.54E-08, avg # of iterations = 4.0 total cpu time spent up to now is 9.4 secs total energy = -162.54728545 Ry Harris-Foulkes estimate = -162.54728785 Ry estimated scf accuracy < 0.00000440 Ry iteration # 7 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.00E-08, avg # of iterations = 2.3 total cpu time spent up to now is 10.2 secs total energy = -162.54728621 Ry Harris-Foulkes estimate = -162.54728622 Ry estimated scf accuracy < 0.00000003 Ry iteration # 8 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.31E-10, avg # of iterations = 5.3 total cpu time spent up to now is 11.5 secs total energy = -162.54728633 Ry Harris-Foulkes estimate = -162.54728634 Ry estimated scf accuracy < 0.00000002 Ry iteration # 9 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.46E-11, avg # of iterations = 1.0 total cpu time spent up to now is 12.2 secs total energy = -162.54728633 Ry Harris-Foulkes estimate = -162.54728633 Ry estimated scf accuracy < 0.00000001 Ry iteration # 10 ecut= 43.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.51E-11, avg # of iterations = 3.0 total cpu time spent up to now is 13.0 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 2301 PWs) bands (ev): -18.7952 -18.7952 -5.2613 -5.2613 -3.1924 -3.1924 -3.1924 -3.1924 1.4261 1.4261 1.4261 1.4261 3.1675 3.1675 3.1675 3.1675 3.3411 3.3411 3.4043 3.4043 10.4390 10.4390 10.4390 10.4390 12.5283 12.5283 14.0274 14.0274 14.0274 14.0274 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 ( 2244 PWs) bands (ev): -18.7852 -18.7852 -5.2997 -5.2997 -3.3147 -3.3147 -3.1830 -3.1830 1.4332 1.4332 1.4580 1.4580 3.1802 3.1802 3.2102 3.2102 3.3444 3.3444 3.7566 3.7566 10.1736 10.1736 10.4665 10.4665 12.0182 12.0182 12.7797 12.7797 12.9131 12.9131 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 ( 2260 PWs) bands (ev): -18.7652 -18.7652 -5.3881 -5.3881 -3.5195 -3.5195 -3.1641 -3.1641 1.4475 1.4475 1.5167 1.5167 3.2062 3.2062 3.2953 3.2953 3.3523 3.3523 4.6492 4.6492 8.6623 8.6623 10.5261 10.5261 11.2996 11.2996 11.7417 11.7417 13.8107 13.8107 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 ( 2240 PWs) bands (ev): -18.7551 -18.7551 -5.4379 -5.4379 -3.6048 -3.6048 -3.1546 -3.1546 1.4549 1.4549 1.5438 1.5438 3.2195 3.2195 3.3378 3.3378 3.3569 3.3569 5.3989 5.3989 7.5664 7.5664 10.5584 10.5584 10.9651 10.9651 11.3695 11.3695 14.4605 14.4605 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 ( 2261 PWs) bands (ev): -18.7755 -18.7755 -5.3316 -5.3316 -3.3662 -3.3662 -3.2307 -3.2307 1.4459 1.4459 1.4677 1.4677 3.1959 3.1959 3.2314 3.2314 3.3279 3.3279 4.0728 4.0728 10.0965 10.0965 10.4721 10.4721 11.3483 11.3483 12.2854 12.2854 13.7827 13.7827 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 ( 2260 PWs) bands (ev): -18.7559 -18.7559 -5.4064 -5.4064 -3.5379 -3.5379 -3.2224 -3.2224 1.4509 1.4509 1.5051 1.5051 3.2095 3.2095 3.2930 3.2930 3.3212 3.3212 4.7981 4.7981 8.8541 8.8541 10.6046 10.6046 11.2393 11.2393 11.8011 11.8011 13.9449 13.9449 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 ( 2250 PWs) bands (ev): -18.7460 -18.7460 -5.4494 -5.4494 -3.6167 -3.6167 -3.2097 -3.2097 1.4520 1.4520 1.5241 1.5241 3.2158 3.2158 3.3023 3.3023 3.3431 3.3431 5.3180 5.3180 8.0393 8.0393 10.8380 10.8380 10.8995 10.8995 11.5408 11.5408 14.0139 14.0139 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 ( 2254 PWs) bands (ev): -18.7371 -18.7371 -5.4536 -5.4536 -3.6223 -3.6223 -3.2625 -3.2625 1.4446 1.4446 1.4918 1.4918 3.1931 3.1931 3.2351 3.2351 3.3546 3.3546 5.0614 5.0614 9.1533 9.1533 10.3437 10.3437 11.4996 11.4996 12.2188 12.2188 13.3752 13.3752 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 ( 2246 PWs) bands (ev): -18.7277 -18.7277 -5.4826 -5.4826 -3.6758 -3.6758 -3.2576 -3.2576 1.4350 1.4350 1.4900 1.4900 3.1771 3.1771 3.2064 3.2064 3.3817 3.3817 5.1929 5.1929 9.1559 9.1559 10.4248 10.4248 11.3585 11.3585 12.3575 12.3575 13.3783 13.3783 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 ( 2216 PWs) bands (ev): -18.7184 -18.7184 -5.5046 -5.5046 -3.7166 -3.7166 -3.2571 -3.2571 1.4200 1.4200 1.4764 1.4764 3.1470 3.1470 3.1621 3.1621 3.4086 3.4086 5.1524 5.1524 10.0242 10.0242 10.3842 10.3842 10.8829 10.8829 13.1274 13.1274 13.1314 13.1314 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.1667 ( 2280 PWs) bands (ev): -18.7659 -18.7659 -5.3591 -5.3591 -3.3662 -3.3662 -3.3227 -3.3227 1.4592 1.4592 1.4732 1.4732 3.2315 3.2315 3.2387 3.2387 3.3118 3.3118 4.3186 4.3186 10.2451 10.2451 10.5910 10.5910 10.8327 10.8327 12.2337 12.2337 12.4539 12.4539 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 ( 2250 PWs) bands (ev): -18.7466 -18.7466 -5.4247 -5.4247 -3.5366 -3.5366 -3.3038 -3.3038 1.4547 1.4547 1.5097 1.5097 3.2369 3.2369 3.2701 3.2701 3.3286 3.3286 4.9060 4.9060 9.1322 9.1322 10.3689 10.3689 11.4150 11.4150 11.6674 11.6674 13.0866 13.0866 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 ( 2236 PWs) bands (ev): -18.7369 -18.7369 -5.4628 -5.4628 -3.6163 -3.6163 -3.2851 -3.2851 1.4501 1.4501 1.5279 1.5279 3.2372 3.2372 3.2601 3.2601 3.3636 3.3636 5.3156 5.3156 8.3876 8.3876 10.2253 10.2253 11.2776 11.2776 12.1839 12.1839 13.3578 13.3578 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 ( 2248 PWs) bands (ev): -18.7279 -18.7279 -5.4707 -5.4707 -3.6022 -3.6022 -3.3703 -3.3703 1.4567 1.4567 1.5177 1.5177 3.2237 3.2237 3.2633 3.2633 3.3660 3.3660 5.1616 5.1616 9.2729 9.2729 9.9497 9.9497 11.4676 11.4676 12.0813 12.0813 13.1230 13.1230 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 ( 2260 PWs) bands (ev): -18.7185 -18.7185 -5.4988 -5.4988 -3.6604 -3.6604 -3.3650 -3.3650 1.4543 1.4543 1.5220 1.5220 3.2045 3.2045 3.2550 3.2550 3.3982 3.3982 5.3249 5.3249 9.0739 9.0739 9.7376 9.7376 11.7305 11.7305 12.4375 12.4375 13.2277 13.2277 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.0062 0.0062 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667-0.5000-0.5000 ( 2268 PWs) bands (ev): -18.7093 -18.7093 -5.5216 -5.5216 -3.7036 -3.7036 -3.3713 -3.3713 1.4629 1.4629 1.5101 1.5101 3.1914 3.1914 3.2339 3.2339 3.4270 3.4270 5.3416 5.3416 9.3502 9.3502 9.7156 9.7156 11.5383 11.5383 12.0109 12.0109 13.8967 13.8970 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.0302 0.0302 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333 0.3333 ( 2266 PWs) bands (ev): -18.7094 -18.7094 -5.5131 -5.5131 -3.5653 -3.5653 -3.5564 -3.5564 1.4894 1.4894 1.5523 1.5523 3.2283 3.2283 3.3716 3.3716 3.3785 3.3785 5.4596 5.4596 8.9347 8.9347 9.3931 9.3931 10.9171 10.9171 13.0493 13.0493 13.1884 13.1884 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.5000 ( 2274 PWs) bands (ev): -18.7001 -18.7001 -5.5387 -5.5387 -3.6439 -3.6439 -3.5402 -3.5402 1.5085 1.5085 1.5635 1.5635 3.2309 3.2309 3.3721 3.3721 3.4293 3.4293 5.6919 5.6919 8.6250 8.6250 9.0488 9.0488 10.6276 10.6276 13.1588 13.1588 13.3665 13.3665 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 ( 2272 PWs) bands (ev): -18.6909 -18.6909 -5.5627 -5.5627 -3.6774 -3.6774 -3.5743 -3.5743 1.5496 1.5496 1.5682 1.5682 3.2460 3.2460 3.4027 3.4027 3.4643 3.4643 5.9638 5.9638 8.4163 8.4163 8.8302 8.8302 9.6814 9.6814 13.3539 13.3539 15.0312 15.0312 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.2784 0.2784 0.0000 0.0000 0.0000 0.0000 k =-0.5000-0.5000-0.5000 ( 2272 PWs) bands (ev): -18.6816 -18.6816 -5.5869 -5.5869 -3.6642 -3.6642 -3.6642 -3.6642 1.5937 1.5937 1.5937 1.5937 3.2726 3.2726 3.4828 3.4828 3.4828 3.4828 6.6804 6.6804 7.9232 7.9232 8.4810 8.4810 8.4810 8.4810 15.3653 15.3653 15.9221 15.9221 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.6685 ev ! total energy = -162.54728633 Ry Harris-Foulkes estimate = -162.54728633 Ry estimated scf accuracy < 1.2E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -22.38120280 Ry hartree contribution = 25.60364023 Ry xc contribution = -52.26722074 Ry ewald contribution = -113.50245427 Ry smearing contrib. (-TS) = -0.00004875 Ry convergence has been achieved in 10 iterations Writing output data file BaHg.save init_run : 0.40s CPU 0.47s WALL ( 1 calls) electrons : 10.10s CPU 10.41s WALL ( 1 calls) Called by init_run: wfcinit : 0.16s CPU 0.19s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 7.56s CPU 7.82s WALL ( 11 calls) sum_band : 1.82s CPU 1.84s WALL ( 11 calls) v_of_rho : 0.04s CPU 0.05s WALL ( 11 calls) v_h : 0.01s CPU 0.00s WALL ( 11 calls) v_xc : 0.04s CPU 0.04s WALL ( 11 calls) newd : 0.68s CPU 0.70s WALL ( 11 calls) mix_rho : 0.03s CPU 0.03s WALL ( 11 calls) Called by c_bands: init_us_2 : 0.02s CPU 0.02s WALL ( 460 calls) cegterg : 7.25s CPU 7.46s WALL ( 220 calls) Called by sum_band: sum_band:bec : 0.42s CPU 0.45s WALL ( 220 calls) addusdens : 0.54s CPU 0.56s WALL ( 11 calls) Called by *egterg: h_psi : 3.98s CPU 4.11s WALL ( 995 calls) s_psi : 0.20s CPU 0.19s WALL ( 995 calls) g_psi : 0.00s CPU 0.01s WALL ( 755 calls) cdiaghg : 2.58s CPU 2.70s WALL ( 955 calls) cegterg:over : 0.20s CPU 0.20s WALL ( 755 calls) cegterg:upda : 0.15s CPU 0.15s WALL ( 755 calls) cegterg:last : 0.06s CPU 0.07s WALL ( 243 calls) cdiaghg:chol : 0.17s CPU 0.15s WALL ( 955 calls) cdiaghg:inve : 0.04s CPU 0.04s WALL ( 955 calls) cdiaghg:para : 0.13s CPU 0.18s WALL ( 1910 calls) Called by h_psi: h_psi:vloc : 3.50s CPU 3.65s WALL ( 995 calls) h_psi:vnl : 0.48s CPU 0.46s WALL ( 995 calls) add_vuspsi : 0.28s CPU 0.25s WALL ( 995 calls) General routines calbec : 0.27s CPU 0.27s WALL ( 1215 calls) fft : 0.13s CPU 0.13s WALL ( 335 calls) ffts : 0.00s CPU 0.01s WALL ( 88 calls) fftw : 3.78s CPU 3.90s WALL ( 90520 calls) interpolate : 0.03s CPU 0.03s WALL ( 88 calls) Parallel routines fft_scatter : 2.05s CPU 2.04s WALL ( 90943 calls) PWSCF : 12.86s CPU 14.05s WALL This run was terminated on: 13:57:48 3Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=