Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 16:32:54 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 Ca.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 3S 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 79 35 9 3201 931 141 Max 80 36 10 3206 952 146 Sum 2869 1265 357 115313 33833 5185 bravais-lattice index = 14 lattice parameter (alat) = 9.2975 a.u. unit-cell volume = 803.6960 (a.u.)^3 number of atoms/cell = 4 number of atomic types = 2 number of electrons = 42.00 number of Kohn-Sham states= 50 kinetic-energy cutoff = 46.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)= 9.297452 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 Ca read from file: /users/gautes/Pseudo/Ca.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 60bd286d514994a6a955285ae841e6ef 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 1169 points, 10 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 l(7) = 2 l(8) = 2 l(9) = 2 l(10) = 2 Q(r) pseudized with 0 coefficients atomic species valence mass pseudopotential Hg 12.00 200.59000 Hg( 1.00) Ca 10.00 40.07800 Ca( 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= 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: 115313 G-vectors FFT dimensions: ( 64, 64, 64) Smooth grid: 33833 G-vectors FFT dimensions: ( 45, 45, 45) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.19 Mb ( 244, 50) NL pseudopotentials 0.25 Mb ( 122, 136) Each V/rho on FFT grid 0.12 Mb ( 8192) Each G-vector array 0.02 Mb ( 3206) G-vector shells 0.00 Mb ( 617) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.74 Mb ( 244, 200) Each subspace H/S matrix 0.04 Mb ( 50, 50) Each matrix 0.21 Mb ( 136, 2, 50) Arrays for rho mixing 1.00 Mb ( 8192, 8) Initial potential from superposition of free atoms starting charge 41.98536, renormalised to 42.00000 Starting wfc are 66 randomized atomic wfcs total cpu time spent up to now is 3.8 secs per-process dynamical memory: 56.3 Mb Self-consistent Calculation iteration # 1 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 9.3 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 2.26E-04, avg # of iterations = 10.1 total cpu time spent up to now is 8.3 secs total energy = -324.87345600 Ry Harris-Foulkes estimate = -324.98956130 Ry estimated scf accuracy < 0.15163142 Ry iteration # 2 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.61E-04, avg # of iterations = 5.2 total cpu time spent up to now is 10.0 secs total energy = -324.90348511 Ry Harris-Foulkes estimate = -325.02224238 Ry estimated scf accuracy < 0.25853730 Ry iteration # 3 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.61E-04, avg # of iterations = 3.5 total cpu time spent up to now is 11.5 secs total energy = -324.95123834 Ry Harris-Foulkes estimate = -324.95229735 Ry estimated scf accuracy < 0.00257650 Ry iteration # 4 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 1 eigenvalues not converged c_bands: 4 eigenvalues not converged c_bands: 2 eigenvalues not converged ethr = 6.13E-06, avg # of iterations = 11.5 total cpu time spent up to now is 13.9 secs total energy = -324.95452881 Ry Harris-Foulkes estimate = -324.95473625 Ry estimated scf accuracy < 0.00059454 Ry iteration # 5 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 1 eigenvalues not converged ethr = 1.42E-06, avg # of iterations = 8.0 total cpu time spent up to now is 15.7 secs total energy = -324.95461634 Ry Harris-Foulkes estimate = -324.95462215 Ry estimated scf accuracy < 0.00002403 Ry iteration # 6 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.72E-08, avg # of iterations = 5.4 total cpu time spent up to now is 17.8 secs total energy = -324.95464611 Ry Harris-Foulkes estimate = -324.95464701 Ry estimated scf accuracy < 0.00000194 Ry iteration # 7 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.63E-09, avg # of iterations = 4.0 total cpu time spent up to now is 19.3 secs total energy = -324.95464645 Ry Harris-Foulkes estimate = -324.95464657 Ry estimated scf accuracy < 0.00000022 Ry iteration # 8 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.29E-10, avg # of iterations = 4.9 total cpu time spent up to now is 21.4 secs total energy = -324.95464666 Ry Harris-Foulkes estimate = -324.95464670 Ry estimated scf accuracy < 0.00000009 Ry iteration # 9 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.03E-10, avg # of iterations = 2.0 total cpu time spent up to now is 22.6 secs total energy = -324.95464666 Ry Harris-Foulkes estimate = -324.95464667 Ry estimated scf accuracy < 0.00000002 Ry iteration # 10 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.88E-11, avg # of iterations = 5.3 total cpu time spent up to now is 24.5 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 4169 PWs) bands (ev): -33.4549 -33.4549 -33.3925 -33.3925 -33.3925 -33.3925 -14.7726 -14.7726 -14.7726 -14.7726 -14.5969 -14.5969 -14.4435 -14.4435 -14.3593 -14.3593 -14.3593 -14.3593 -14.2824 -14.2824 -14.1821 -14.1821 -14.1821 -14.1821 0.1372 0.1372 0.1372 0.1372 1.9035 1.9035 1.9035 1.9035 1.9461 1.9461 2.0004 2.0004 6.3780 6.3780 8.7718 8.7718 8.7756 8.7756 8.7756 8.7757 9.1975 9.1988 9.2202 9.4025 9.4040 9.4041 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 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.0000 0.0000 0.2000 ( 4212 PWs) bands (ev): -33.4496 -33.4496 -33.3978 -33.3978 -33.3925 -33.3925 -14.7684 -14.7684 -14.7539 -14.7539 -14.6042 -14.6042 -14.4367 -14.4367 -14.3596 -14.3596 -14.3526 -14.3526 -14.2873 -14.2873 -14.2078 -14.2078 -14.1916 -14.1916 0.1428 0.1428 0.1509 0.1509 1.9157 1.9157 1.9269 1.9269 1.9533 1.9533 2.2519 2.2519 6.2395 6.2395 7.7136 7.7136 8.6179 8.6179 8.8810 8.8810 9.0252 9.0252 9.2686 9.2686 9.5447 9.5447 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 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.4000 ( 4220 PWs) bands (ev): -33.4375 -33.4375 -33.4099 -33.4099 -33.3925 -33.3925 -14.7611 -14.7611 -14.7083 -14.7083 -14.6343 -14.6343 -14.4238 -14.4238 -14.3390 -14.3390 -14.3329 -14.3329 -14.3030 -14.3030 -14.2687 -14.2687 -14.2076 -14.2076 0.1519 0.1519 0.1723 0.1723 1.9351 1.9351 1.9613 1.9613 1.9652 1.9652 2.8499 2.8499 5.3054 5.3054 7.6944 7.6944 8.3111 8.3111 8.5918 8.5918 9.3803 9.3803 9.6164 9.6164 9.9001 9.9003 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.2000 ( 4222 PWs) bands (ev): -33.4448 -33.4448 -33.3997 -33.3997 -33.3954 -33.3954 -14.7601 -14.7601 -14.7465 -14.7465 -14.6150 -14.6150 -14.4275 -14.4275 -14.3588 -14.3588 -14.3520 -14.3520 -14.2833 -14.2833 -14.2315 -14.2315 -14.1959 -14.1959 0.1496 0.1496 0.1588 0.1588 1.9291 1.9291 1.9480 1.9480 1.9561 1.9561 2.4798 2.4798 6.2449 6.2449 7.3878 7.3878 8.1138 8.1138 8.6352 8.6352 9.1513 9.1513 9.4844 9.4845 9.5214 9.5214 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 0.2306 0.2306 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.4000 ( 4228 PWs) bands (ev): -33.4334 -33.4334 -33.4099 -33.4099 -33.3963 -33.3963 -14.7527 -14.7527 -14.7154 -14.7154 -14.6484 -14.6484 -14.4131 -14.4131 -14.3549 -14.3549 -14.3341 -14.3341 -14.2886 -14.2886 -14.2629 -14.2629 -14.2146 -14.2146 0.1583 0.1583 0.1731 0.1731 1.9475 1.9475 1.9569 1.9569 1.9840 1.9840 3.0096 3.0096 5.5557 5.5557 7.3413 7.3413 7.9380 7.9380 8.4640 8.4640 9.3574 9.3574 9.4283 9.4283 9.8969 9.8990 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4000 0.4000 ( 4252 PWs) bands (ev): -33.4233 -33.4233 -33.4112 -33.4112 -33.4050 -33.4050 -14.7348 -14.7348 -14.7228 -14.7228 -14.6792 -14.6792 -14.3903 -14.3903 -14.3716 -14.3716 -14.3256 -14.3256 -14.3040 -14.3040 -14.2414 -14.2414 -14.2248 -14.2248 0.1676 0.1676 0.1745 0.1745 1.9407 1.9407 1.9666 1.9666 2.0191 2.0191 3.3686 3.3686 6.0608 6.0608 6.4360 6.4360 7.5697 7.5697 8.4291 8.4291 9.2820 9.2820 9.4447 9.4447 9.7114 9.7114 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.2000 ( 4244 PWs) bands (ev): -33.4404 -33.4404 -33.3997 -33.3997 -33.3997 -33.3997 -14.7604 -14.7604 -14.7385 -14.7385 -14.6292 -14.6292 -14.4151 -14.4151 -14.3696 -14.3696 -14.3441 -14.3441 -14.2759 -14.2759 -14.2504 -14.2504 -14.1949 -14.1949 0.1530 0.1530 0.1650 0.1650 1.9362 1.9362 1.9559 1.9559 1.9651 1.9651 2.6889 2.6889 6.2616 6.2616 7.4530 7.4530 7.7226 7.7226 8.2171 8.2171 9.3893 9.3893 9.4023 9.4023 9.5642 9.5642 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.0002 0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.4000 ( 4234 PWs) bands (ev): -33.4300 -33.4300 -33.4100 -33.4100 -33.3996 -33.3996 -14.7552 -14.7552 -14.7195 -14.7195 -14.6650 -14.6650 -14.3999 -14.3999 -14.3774 -14.3774 -14.3225 -14.3225 -14.2859 -14.2859 -14.2538 -14.2538 -14.2095 -14.2095 0.1599 0.1599 0.1732 0.1732 1.9410 1.9410 1.9628 1.9628 1.9923 1.9923 3.1480 3.1480 5.7881 5.7881 7.2396 7.2396 7.9241 7.9241 8.1052 8.1052 9.1550 9.1550 9.3240 9.3241 9.8251 9.8252 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.3597 0.3597 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.4000 0.4000 ( 4248 PWs) bands (ev): -33.4215 -33.4215 -33.4111 -33.4111 -33.4067 -33.4067 -14.7560 -14.7560 -14.7236 -14.7236 -14.6960 -14.6960 -14.4040 -14.4040 -14.3745 -14.3745 -14.3102 -14.3102 -14.2919 -14.2919 -14.2298 -14.2298 -14.2094 -14.2094 0.1638 0.1638 0.1714 0.1714 1.9228 1.9228 1.9762 1.9762 2.0097 2.0097 3.4297 3.4297 6.2626 6.2626 6.7056 6.7056 7.7425 7.7425 8.0914 8.0914 8.9313 8.9313 9.2366 9.2366 9.3352 9.3352 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.6081 0.6081 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4000 0.4000 0.4000 ( 4236 PWs) bands (ev): -33.4170 -33.4170 -33.4110 -33.4110 -33.4110 -33.4110 -14.7835 -14.7835 -14.7286 -14.7286 -14.7210 -14.7210 -14.4340 -14.4340 -14.3586 -14.3586 -14.2907 -14.2907 -14.2843 -14.2843 -14.2040 -14.2040 -14.1924 -14.1924 0.1594 0.1594 0.1640 0.1640 1.8920 1.8920 1.9917 1.9917 1.9999 1.9999 3.5386 3.5386 6.7661 6.7661 7.3269 7.3269 7.5220 7.5220 8.1956 8.1956 8.3922 8.3922 8.4487 8.4487 8.9343 8.9343 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.0007 0.0007 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 8.0974 ev ! total energy = -324.95464667 Ry Harris-Foulkes estimate = -324.95464667 Ry estimated scf accuracy < 4.0E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -107.15667722 Ry hartree contribution = 68.72050491 Ry xc contribution = -68.31952823 Ry ewald contribution = -218.19833346 Ry smearing contrib. (-TS) = -0.00061268 Ry convergence has been achieved in 10 iterations Writing output data file Ca3Hg.save init_run : 1.08s CPU 1.58s WALL ( 1 calls) electrons : 19.44s CPU 21.01s WALL ( 1 calls) Called by init_run: wfcinit : 0.39s CPU 0.64s WALL ( 1 calls) potinit : 0.02s CPU 0.03s WALL ( 1 calls) Called by electrons: c_bands : 15.52s CPU 16.48s WALL ( 11 calls) sum_band : 2.81s CPU 2.85s WALL ( 11 calls) v_of_rho : 0.06s CPU 0.06s WALL ( 11 calls) v_h : 0.00s CPU 0.00s WALL ( 11 calls) v_xc : 0.06s CPU 0.06s WALL ( 11 calls) newd : 1.06s CPU 1.10s WALL ( 11 calls) mix_rho : 0.04s CPU 0.04s WALL ( 11 calls) Called by c_bands: init_us_2 : 0.03s CPU 0.03s WALL ( 230 calls) cegterg : 15.02s CPU 15.27s WALL ( 110 calls) Called by sum_band: sum_band:bec : 0.45s CPU 0.46s WALL ( 110 calls) addusdens : 0.91s CPU 0.92s WALL ( 11 calls) Called by *egterg: h_psi : 8.64s CPU 8.84s WALL ( 812 calls) s_psi : 0.50s CPU 0.48s WALL ( 812 calls) g_psi : 0.02s CPU 0.02s WALL ( 692 calls) cdiaghg : 4.84s CPU 4.94s WALL ( 792 calls) cegterg:over : 0.42s CPU 0.43s WALL ( 692 calls) cegterg:upda : 0.44s CPU 0.44s WALL ( 692 calls) cegterg:last : 0.14s CPU 0.14s WALL ( 129 calls) cdiaghg:chol : 0.26s CPU 0.28s WALL ( 792 calls) cdiaghg:inve : 0.13s CPU 0.18s WALL ( 792 calls) cdiaghg:para : 0.32s CPU 0.30s WALL ( 1584 calls) Called by h_psi: h_psi:vloc : 7.48s CPU 7.70s WALL ( 812 calls) h_psi:vnl : 1.12s CPU 1.11s WALL ( 812 calls) add_vuspsi : 0.56s CPU 0.56s WALL ( 812 calls) General routines calbec : 0.69s CPU 0.68s WALL ( 922 calls) fft : 0.12s CPU 0.13s WALL ( 335 calls) ffts : 0.03s CPU 0.02s WALL ( 88 calls) fftw : 8.08s CPU 8.25s WALL ( 83804 calls) interpolate : 0.06s CPU 0.06s WALL ( 88 calls) Parallel routines fft_scatter : 2.66s CPU 2.69s WALL ( 84227 calls) PWSCF : 22.83s CPU 27.39s WALL This run was terminated on: 16:33:22 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=