Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 19:40: 7 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 file Mg.rel-pbe-spnl-rrkjus_psl.1.0.0.UPF: wavefunction(s) 2P 2P 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 39 22 6 1089 458 76 Max 40 23 7 1093 472 79 Sum 1405 793 241 39223 16831 2777 bravais-lattice index = 14 lattice parameter (alat) = 6.5007 a.u. unit-cell volume = 274.7083 (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 = 59.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)= 6.500657 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 Mg read from file: /users/gautes/Pseudo/Mg.rel-pbe-spnl-rrkjus_psl.1.0.0.UPF MD5 check sum: 99618f2d1dea131ec68d28ce241ef835 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 1129 points, 6 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 Q(r) pseudized with 0 coefficients atomic species valence mass pseudopotential Hg 12.00 200.59000 Hg( 1.00) Mg 10.00 24.30500 Mg( 1.00) 48 Sym. Ops., with inversion, found s frac. trans. isym = 1 identity cryst. s( 1) = ( 1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) cart. s( 1) = ( 1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 2 180 deg rotation - cart. axis [0,0,1] cryst. s( 2) = ( -1 0 0 ) ( 0 -1 0 ) ( 0 0 1 ) cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 3 180 deg rotation - cart. axis [0,1,0] cryst. s( 3) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 -1 ) cart. s( 3) = ( -1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 4 180 deg rotation - cart. axis [1,0,0] cryst. s( 4) = ( 1 0 0 ) ( 0 -1 0 ) ( 0 0 -1 ) cart. s( 4) = ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 5 180 deg rotation - cart. axis [1,1,0] cryst. s( 5) = ( 0 1 0 ) ( 1 0 0 ) ( 0 0 -1 ) cart. s( 5) = ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 6 180 deg rotation - cart. axis [1,-1,0] cryst. s( 6) = ( 0 -1 0 ) ( -1 0 0 ) ( 0 0 -1 ) cart. s( 6) = ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 7 90 deg rotation - cart. axis [0,0,-1] cryst. s( 7) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 1 ) cart. s( 7) = ( -0.0000000 1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 8 90 deg rotation - cart. axis [0,0,1] cryst. s( 8) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 1 ) cart. s( 8) = ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 9 180 deg rotation - cart. axis [1,0,1] cryst. s( 9) = ( 0 0 1 ) ( 0 -1 0 ) ( 1 0 0 ) cart. s( 9) = ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 10 180 deg rotation - cart. axis [-1,0,1] cryst. s(10) = ( 0 0 -1 ) ( 0 -1 0 ) ( -1 0 0 ) cart. s(10) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 11 90 deg rotation - cart. axis [0,1,0] cryst. s(11) = ( 0 0 -1 ) ( 0 1 0 ) ( 1 0 0 ) cart. s(11) = ( -0.0000000 0.0000000 1.0000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 12 90 deg rotation - cart. axis [0,-1,0] cryst. s(12) = ( 0 0 1 ) ( 0 1 0 ) ( -1 0 0 ) cart. s(12) = ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 13 180 deg rotation - cart. axis [0,1,1] cryst. s(13) = ( -1 0 0 ) ( 0 0 1 ) ( 0 1 0 ) cart. s(13) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 14 180 deg rotation - cart. axis [0,1,-1] cryst. s(14) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) cart. s(14) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 15 90 deg rotation - cart. axis [-1,0,0] cryst. s(15) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) cart. s(15) = ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 16 90 deg rotation - cart. axis [1,0,0] cryst. s(16) = ( 1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) cart. s(16) = ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 17 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(17) = ( 0 0 1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(17) = ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 18 120 deg rotation - cart. axis [-1,1,1] cryst. s(18) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(18) = ( -0.0000000 -1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 19 120 deg rotation - cart. axis [1,1,-1] cryst. s(19) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(19) = ( -0.0000000 1.0000000 -0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 20 120 deg rotation - cart. axis [1,-1,1] cryst. s(20) = ( 0 0 1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(20) = ( 0.0000000 -1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 21 120 deg rotation - cart. axis [1,1,1] cryst. s(21) = ( 0 1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(21) = ( 0.0000000 0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 22 120 deg rotation - cart. axis [-1,1,-1] cryst. s(22) = ( 0 -1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(22) = ( -0.0000000 -0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 23 120 deg rotation - cart. axis [1,-1,-1] cryst. s(23) = ( 0 -1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(23) = ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 24 120 deg rotation - cart. axis [-1,-1,1] cryst. s(24) = ( 0 1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(24) = ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 25 inversion cryst. s(25) = ( -1 0 0 ) ( 0 -1 0 ) ( 0 0 -1 ) cart. s(25) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 26 inv. 180 deg rotation - cart. axis [0,0,1] cryst. s(26) = ( 1 0 0 ) ( 0 1 0 ) ( 0 0 -1 ) cart. s(26) = ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 27 inv. 180 deg rotation - cart. axis [0,1,0] cryst. s(27) = ( 1 0 0 ) ( 0 -1 0 ) ( 0 0 1 ) cart. s(27) = ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 28 inv. 180 deg rotation - cart. axis [1,0,0] cryst. s(28) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) cart. s(28) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 29 inv. 180 deg rotation - cart. axis [1,1,0] cryst. s(29) = ( 0 -1 0 ) ( -1 0 0 ) ( 0 0 1 ) cart. s(29) = ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 30 inv. 180 deg rotation - cart. axis [1,-1,0] cryst. s(30) = ( 0 1 0 ) ( 1 0 0 ) ( 0 0 1 ) cart. s(30) = ( 0.0000000 1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 31 inv. 90 deg rotation - cart. axis [0,0,-1] cryst. s(31) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 -1 ) cart. s(31) = ( 0.0000000 -1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 32 inv. 90 deg rotation - cart. axis [0,0,1] cryst. s(32) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 -1 ) cart. s(32) = ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 33 inv. 180 deg rotation - cart. axis [1,0,1] cryst. s(33) = ( 0 0 -1 ) ( 0 1 0 ) ( -1 0 0 ) cart. s(33) = ( -0.0000000 0.0000000 -1.0000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 34 inv. 180 deg rotation - cart. axis [-1,0,1] cryst. s(34) = ( 0 0 1 ) ( 0 1 0 ) ( 1 0 0 ) cart. s(34) = ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 35 inv. 90 deg rotation - cart. axis [0,1,0] cryst. s(35) = ( 0 0 1 ) ( 0 -1 0 ) ( -1 0 0 ) cart. s(35) = ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 36 inv. 90 deg rotation - cart. axis [0,-1,0] cryst. s(36) = ( 0 0 -1 ) ( 0 -1 0 ) ( 1 0 0 ) cart. s(36) = ( -0.0000000 -0.0000000 1.0000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 37 inv. 180 deg rotation - cart. axis [0,1,1] cryst. s(37) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) cart. s(37) = ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 38 inv. 180 deg rotation - cart. axis [0,1,-1] cryst. s(38) = ( 1 0 0 ) ( 0 0 1 ) ( 0 1 0 ) cart. s(38) = ( 1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 39 inv. 90 deg rotation - cart. axis [-1,0,0] cryst. s(39) = ( -1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) cart. s(39) = ( -1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 40 inv. 90 deg rotation - cart. axis [1,0,0] cryst. s(40) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) cart. s(40) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 41 inv. 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(41) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(41) = ( -0.0000000 -1.0000000 0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 42 inv. 120 deg rotation - cart. axis [-1,1,1] cryst. s(42) = ( 0 0 1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(42) = ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 43 inv. 120 deg rotation - cart. axis [1,1,-1] cryst. s(43) = ( 0 0 1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(43) = ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 44 inv. 120 deg rotation - cart. axis [1,-1,1] cryst. s(44) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(44) = ( -0.0000000 1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 45 inv. 120 deg rotation - cart. axis [1,1,1] cryst. s(45) = ( 0 -1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(45) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 46 inv. 120 deg rotation - cart. axis [-1,1,-1] cryst. s(46) = ( 0 1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(46) = ( 0.0000000 0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 47 inv. 120 deg rotation - cart. axis [1,-1,-1] cryst. s(47) = ( 0 1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(47) = ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 48 inv. 120 deg rotation - cart. axis [-1,-1,1] cryst. s(48) = ( 0 -1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(48) = ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) double point group O_h (m-3m) there are 16 classes and 6 irreducible representations the character table: E -E 8C3 -8C3 3C2 6C4 -6C4 6C2' i -i 8S6 -8S6 -3C2 -6C2' G_6+ 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 2.00 -2.00 1.00 -1.00 G_7+ 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 2.00 -2.00 1.00 -1.00 G_8+ 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 4.00 -4.00 -1.00 1.00 G_6- 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 -2.00 2.00 -1.00 1.00 G_7- 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 -2.00 2.00 -1.00 1.00 G_8- 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 -4.00 4.00 1.00 -1.00 3s_h 6S4 -6S4 6s_d -3s_h -6s_d G_6+ 0.00 1.41 -1.41 0.00 G_7+ 0.00 -1.41 1.41 0.00 G_8+ 0.00 0.00 0.00 0.00 G_6- 0.00 -1.41 1.41 0.00 G_7- 0.00 1.41 -1.41 0.00 G_8- 0.00 0.00 0.00 0.00 the symmetry operations in each class and the name of the first element: E 1 identity -E -1 identity E 8C3 17 19 20 18 24 21 22 23 120 deg rotation - cart. axis [-1,-1,-1] -8C3 -17 -19 -20 -18 -24 -21 -22 -23 120 deg rotation - cart. axis [-1,-1,-1] E 3C2 -3C2 2 -2 4 -4 3 -3 180 deg rotation - cart. axis [0,0,1] 6C4 7 8 15 16 12 11 90 deg rotation - cart. axis [0,0,-1] -6C4 -7 -8 -15 -16 -12 -11 90 deg rotation - cart. axis [0,0,-1] E 6C2'-6C2' 5 -5 6 -6 14 -13 -14 13 -10 -9 10 9 180 deg rotation - cart. axis [1,1,0] i 25 inversion -i -25 inversion E 8S6 41 43 44 42 48 45 46 47 inv. 120 deg rotation - cart. axis [-1,-1,-1] -8S6 -41 -43 -44 -42 -48 -45 -46 -47 inv. 120 deg rotation - cart. axis [-1,-1,-1] E 3s_h-3s_h 26 -26 28 -28 27 -27 inv. 180 deg rotation - cart. axis [0,0,1] 6S4 31 32 39 40 36 35 inv. 90 deg rotation - cart. axis [0,0,-1] -6S4 -31 -32 -39 -40 -36 -35 inv. 90 deg rotation - cart. axis [0,0,-1] E 6s_d-6s_d 29 -29 30 -30 38 -37 -38 37 -34 -33 34 33 inv. 180 deg rotation - cart. axis [1,1,0] Cartesian axes number of k points= 20 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 Dense grid: 39223 G-vectors FFT dimensions: ( 45, 45, 45) Smooth grid: 16831 G-vectors FFT dimensions: ( 32, 32, 32) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.06 Mb ( 126, 30) NL pseudopotentials 0.05 Mb ( 63, 48) Each V/rho on FFT grid 0.06 Mb ( 4050) Each G-vector array 0.01 Mb ( 1091) G-vector shells 0.00 Mb ( 262) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.23 Mb ( 126, 120) Each subspace H/S matrix 0.01 Mb ( 30, 30) Each matrix 0.04 Mb ( 48, 2, 30) Arrays for rho mixing 0.49 Mb ( 4050, 8) Check: negative/imaginary core charge= -0.000009 0.000000 Initial potential from superposition of free atoms starting charge 21.99900, renormalised to 22.00000 Starting wfc are 28 randomized atomic wfcs + 2 random wfc total cpu time spent up to now is 1.8 secs per-process dynamical memory: 19.9 Mb Self-consistent Calculation iteration # 1 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 4.1 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.67E-04, avg # of iterations = 4.0 total cpu time spent up to now is 3.2 secs total energy = -226.40922435 Ry Harris-Foulkes estimate = -226.43604605 Ry estimated scf accuracy < 0.04276809 Ry iteration # 2 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.94E-04, avg # of iterations = 3.1 total cpu time spent up to now is 3.9 secs total energy = -226.41923076 Ry Harris-Foulkes estimate = -226.43003945 Ry estimated scf accuracy < 0.01751495 Ry iteration # 3 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.96E-05, avg # of iterations = 2.0 total cpu time spent up to now is 4.5 secs total energy = -226.42343760 Ry Harris-Foulkes estimate = -226.42347361 Ry estimated scf accuracy < 0.00030865 Ry iteration # 4 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.40E-06, avg # of iterations = 4.3 total cpu time spent up to now is 5.4 secs total energy = -226.42372213 Ry Harris-Foulkes estimate = -226.42377894 Ry estimated scf accuracy < 0.00011831 Ry iteration # 5 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.38E-07, avg # of iterations = 1.9 total cpu time spent up to now is 5.9 secs total energy = -226.42372486 Ry Harris-Foulkes estimate = -226.42373374 Ry estimated scf accuracy < 0.00001246 Ry iteration # 6 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.66E-08, avg # of iterations = 4.0 total cpu time spent up to now is 6.7 secs total energy = -226.42373406 Ry Harris-Foulkes estimate = -226.42373785 Ry estimated scf accuracy < 0.00000917 Ry iteration # 7 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.17E-08, avg # of iterations = 1.1 total cpu time spent up to now is 7.1 secs total energy = -226.42373358 Ry Harris-Foulkes estimate = -226.42373450 Ry estimated scf accuracy < 0.00000206 Ry iteration # 8 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.36E-09, avg # of iterations = 3.1 total cpu time spent up to now is 7.8 secs total energy = -226.42373422 Ry Harris-Foulkes estimate = -226.42373423 Ry estimated scf accuracy < 0.00000002 Ry iteration # 9 ecut= 59.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.00E-11, avg # of iterations = 3.0 total cpu time spent up to now is 8.4 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 2103 PWs) bands (ev): -65.8586 -65.8586 -32.7953 -32.7953 -32.5071 -32.5071 -32.5071 -32.5071 1.3370 1.3370 1.8857 1.8857 1.8857 1.8857 3.6915 3.6915 3.6915 3.6915 4.1919 4.1919 13.6579 13.6579 14.7882 14.7882 14.7882 14.7882 17.2446 17.2446 17.4912 17.4912 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.1429 ( 2088 PWs) bands (ev): -65.8585 -65.8585 -32.7957 -32.7957 -32.5080 -32.5080 -32.5070 -32.5070 1.5957 1.5957 1.8928 1.8928 1.9508 1.9508 3.6885 3.6885 3.6917 3.6917 4.1452 4.1452 12.0875 12.0875 14.4190 14.4190 15.1388 15.1388 17.4743 17.4743 17.5254 17.5254 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.2857 ( 2084 PWs) bands (ev): -65.8584 -65.8584 -32.7965 -32.7965 -32.5100 -32.5100 -32.5068 -32.5068 1.8501 1.8501 1.9067 1.9067 2.5392 2.5392 3.6669 3.6669 3.8128 3.8128 4.0650 4.0650 9.2419 9.2419 15.2588 15.2588 16.1008 16.1008 17.4311 17.4311 18.0863 18.0863 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.4286 ( 2088 PWs) bands (ev): -65.8583 -65.8583 -32.7972 -32.7972 -32.5116 -32.5116 -32.5066 -32.5066 1.8839 1.8839 1.9153 1.9153 3.1498 3.1498 3.6324 3.6324 4.0272 4.0272 4.7008 4.7008 6.7320 6.7320 16.3206 16.3206 17.3636 17.3636 17.3922 17.3922 19.0379 19.0379 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.1429 0.1429 ( 2085 PWs) bands (ev): -65.8584 -65.8584 -32.7960 -32.7960 -32.5084 -32.5084 -32.5074 -32.5074 1.8083 1.8083 1.9273 1.9273 2.0245 2.0245 3.6656 3.6656 3.7182 3.7182 4.0954 4.0954 12.1811 12.1811 12.6769 12.6769 15.1375 15.1375 17.5578 17.5578 19.9103 19.9104 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.1429 0.2857 ( 2092 PWs) bands (ev): -65.8583 -65.8583 -32.7969 -32.7969 -32.5102 -32.5102 -32.5074 -32.5074 1.8926 1.8926 1.9441 1.9441 2.7330 2.7330 3.6462 3.6462 3.8588 3.8588 4.0248 4.0248 9.5447 9.5447 13.2718 13.2718 16.0525 16.0525 17.5181 17.5181 19.6641 19.6641 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.1429 0.4286 ( 2107 PWs) bands (ev): -65.8582 -65.8582 -32.7976 -32.7976 -32.5118 -32.5118 -32.5072 -32.5072 1.9053 1.9053 1.9424 1.9424 3.2564 3.2564 3.6033 3.6033 3.9604 3.9604 4.8701 4.8701 7.0554 7.0554 14.4233 14.4233 16.5190 16.5190 17.3005 17.3005 19.1394 19.1394 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.2857 0.2857 ( 2101 PWs) bands (ev): -65.8582 -65.8582 -32.7977 -32.7977 -32.5113 -32.5113 -32.5080 -32.5080 1.9116 1.9116 1.9583 1.9583 3.2866 3.2866 3.6403 3.6403 3.9557 3.9557 4.1746 4.1746 10.1731 10.1731 10.6165 10.6165 16.1771 16.1771 17.0384 17.0384 19.1312 19.1312 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.0259 0.0259 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.4286 ( 2110 PWs) bands (ev): -65.8581 -65.8581 -32.7984 -32.7984 -32.5125 -32.5125 -32.5081 -32.5081 1.9104 1.9104 1.9264 1.9264 3.5095 3.5095 3.5940 3.5940 3.9849 3.9849 5.3744 5.3744 7.9112 7.9112 11.5502 11.5502 14.3872 14.3872 18.3126 18.3126 19.6339 19.6339 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.4286 0.4286 ( 2104 PWs) bands (ev): -65.8580 -65.8580 -32.7990 -32.7990 -32.5134 -32.5134 -32.5084 -32.5084 1.8406 1.8406 1.9223 1.9223 3.5594 3.5594 3.5791 3.5791 4.2135 4.2135 6.5058 6.5058 8.6730 8.6730 9.3376 9.3376 12.6314 12.6314 19.6388 19.6388 20.3911 20.3911 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.1429 0.1429 0.1429 ( 2103 PWs) bands (ev): -65.8584 -65.8584 -32.7964 -32.7964 -32.5084 -32.5084 -32.5083 -32.5083 1.9222 1.9222 1.9376 1.9376 2.2106 2.2106 3.6671 3.6671 3.7271 3.7271 4.0436 4.0436 12.3100 12.3100 12.8795 12.8795 13.0949 13.0949 19.5838 19.5838 19.8605 19.8605 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.1429 0.1429 0.2857 ( 2108 PWs) bands (ev): -65.8583 -65.8583 -32.7973 -32.7973 -32.5103 -32.5103 -32.5081 -32.5081 1.9415 1.9415 1.9516 1.9516 2.9322 2.9322 3.6551 3.6551 3.8517 3.8517 4.0240 4.0240 9.8484 9.8484 13.3591 13.3591 13.8672 13.8672 18.4700 18.4700 19.4968 19.4968 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.1429 0.1429 0.4286 ( 2095 PWs) bands (ev): -65.8582 -65.8582 -32.7979 -32.7979 -32.5118 -32.5118 -32.5079 -32.5079 1.9363 1.9363 1.9494 1.9494 3.3572 3.3572 3.5957 3.5957 3.8920 3.8920 5.0336 5.0336 7.3828 7.3828 14.4587 14.4587 15.1334 15.1334 16.7139 16.7139 17.5291 17.5291 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.1429 0.2857 0.2857 ( 2102 PWs) bands (ev): -65.8581 -65.8581 -32.7981 -32.7981 -32.5111 -32.5111 -32.5089 -32.5089 1.9497 1.9497 1.9590 1.9590 3.3939 3.3939 3.6558 3.6558 3.8819 3.8819 4.3373 4.3373 10.4984 10.4984 10.9073 10.9073 14.4751 14.4751 16.6433 16.6433 19.6581 19.6581 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.9937 0.9937 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.2857 0.4286 ( 2104 PWs) bands (ev): -65.8580 -65.8580 -32.7987 -32.7987 -32.5124 -32.5124 -32.5089 -32.5089 1.9178 1.9178 1.9507 1.9507 3.5259 3.5259 3.6106 3.6106 3.9459 3.9459 5.5766 5.5766 8.2454 8.2454 11.8715 11.8715 14.6000 14.6000 15.8119 15.8119 17.8150 17.8150 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.1429 0.4286 0.4286 ( 2103 PWs) bands (ev): -65.8579 -65.8579 -32.7993 -32.7993 -32.5133 -32.5133 -32.5093 -32.5093 1.8494 1.8494 1.9461 1.9461 3.4909 3.4909 3.6334 3.6334 4.1863 4.1863 6.7363 6.7363 9.0295 9.0295 9.6661 9.6661 12.8146 12.8146 17.0837 17.0837 18.7043 18.7044 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.2857 0.2857 0.2857 ( 2100 PWs) bands (ev): -65.8580 -65.8580 -32.7988 -32.7988 -32.5109 -32.5109 -32.5108 -32.5108 1.9504 1.9504 1.9574 1.9574 3.5084 3.5084 3.6991 3.6991 3.9061 3.9061 4.9660 4.9660 11.2118 11.2118 11.7444 11.7444 11.7683 11.7683 17.0056 17.0056 17.5290 17.5292 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.2857 0.2857 0.4286 ( 2116 PWs) bands (ev): -65.8579 -65.8579 -32.7995 -32.7995 -32.5124 -32.5124 -32.5107 -32.5107 1.9077 1.9077 1.9456 1.9456 3.4291 3.4291 3.7980 3.7980 3.9859 3.9859 6.2335 6.2335 9.1244 9.1244 12.5867 12.5867 13.0015 13.0015 15.1571 15.1571 15.4481 15.4481 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.2857 0.4286 0.4286 ( 2110 PWs) bands (ev): -65.8578 -65.8578 -32.8002 -32.8002 -32.5131 -32.5131 -32.5113 -32.5113 1.8598 1.8598 1.9242 1.9242 3.3311 3.3311 3.8859 3.8859 4.1678 4.1678 7.4341 7.4341 9.9467 9.9467 10.5796 10.5796 13.3517 13.3517 14.1247 14.1247 16.2440 16.2440 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.2851 0.2851 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4286 0.4286 0.4286 ( 2128 PWs) bands (ev): -65.8577 -65.8577 -32.8009 -32.8009 -32.5130 -32.5130 -32.5129 -32.5129 1.8553 1.8553 1.8604 1.8604 3.2200 3.2200 4.1515 4.1515 4.1848 4.1848 8.6341 8.6341 10.7971 10.7971 11.6963 11.6963 11.9450 11.9450 13.9141 13.9141 14.7064 14.7065 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 the Fermi energy is 10.5671 ev ! total energy = -226.42373422 Ry Harris-Foulkes estimate = -226.42373422 Ry estimated scf accuracy < 2.0E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -102.69081218 Ry hartree contribution = 61.95618871 Ry xc contribution = -49.58505291 Ry ewald contribution = -136.10396042 Ry smearing contrib. (-TS) = -0.00009742 Ry convergence has been achieved in 9 iterations Writing output data file MgHg.save init_run : 0.29s CPU 0.35s WALL ( 1 calls) electrons : 6.45s CPU 6.71s WALL ( 1 calls) Called by init_run: wfcinit : 0.15s CPU 0.16s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 5.28s CPU 5.48s WALL ( 10 calls) sum_band : 0.89s CPU 0.92s WALL ( 10 calls) v_of_rho : 0.04s CPU 0.03s WALL ( 10 calls) v_h : 0.01s CPU 0.00s WALL ( 10 calls) v_xc : 0.03s CPU 0.03s WALL ( 10 calls) newd : 0.23s CPU 0.24s WALL ( 10 calls) mix_rho : 0.01s CPU 0.02s WALL ( 10 calls) Called by c_bands: init_us_2 : 0.02s CPU 0.01s WALL ( 420 calls) cegterg : 5.12s CPU 5.27s WALL ( 200 calls) Called by sum_band: sum_band:bec : 0.22s CPU 0.21s WALL ( 200 calls) addusdens : 0.16s CPU 0.15s WALL ( 10 calls) Called by *egterg: h_psi : 2.50s CPU 2.66s WALL ( 835 calls) s_psi : 0.08s CPU 0.10s WALL ( 835 calls) g_psi : 0.01s CPU 0.01s WALL ( 615 calls) cdiaghg : 2.28s CPU 2.19s WALL ( 795 calls) cegterg:over : 0.13s CPU 0.16s WALL ( 615 calls) cegterg:upda : 0.08s CPU 0.12s WALL ( 615 calls) cegterg:last : 0.05s CPU 0.05s WALL ( 200 calls) cdiaghg:chol : 0.11s CPU 0.12s WALL ( 795 calls) cdiaghg:inve : 0.02s CPU 0.03s WALL ( 795 calls) cdiaghg:para : 0.12s CPU 0.15s WALL ( 1590 calls) Called by h_psi: h_psi:vloc : 2.20s CPU 2.41s WALL ( 835 calls) h_psi:vnl : 0.29s CPU 0.25s WALL ( 835 calls) add_vuspsi : 0.12s CPU 0.13s WALL ( 835 calls) General routines calbec : 0.18s CPU 0.15s WALL ( 1035 calls) fft : 0.07s CPU 0.06s WALL ( 304 calls) ffts : 0.00s CPU 0.01s WALL ( 80 calls) fftw : 2.44s CPU 2.59s WALL ( 78960 calls) interpolate : 0.02s CPU 0.02s WALL ( 80 calls) Parallel routines fft_scatter : 1.59s CPU 1.67s WALL ( 79344 calls) PWSCF : 8.41s CPU 9.38s WALL This run was terminated on: 19:40:16 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=