Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 15:56:45 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 Br.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4S 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 109 83 23 5156 3460 509 Max 111 84 24 5158 3492 512 Sum 3945 3017 845 185639 125039 18349 bravais-lattice index = 14 lattice parameter (alat) = 10.9037 a.u. unit-cell volume = 1296.3550 (a.u.)^3 number of atoms/cell = 5 number of atomic types = 3 number of electrons = 42.00 number of Kohn-Sham states= 50 kinetic-energy cutoff = 80.0000 Ry charge density cutoff = 416.0000 Ry convergence threshold = 1.0E-08 mixing beta = 0.7000 number of iterations used = 8 plain mixing Exchange-correlation = PBE ( 1 4 3 4 0 0) Non magnetic calculation with spin-orbit celldm(1)= 10.903719 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 Br read from file: /users/gautes/Pseudo/Br.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 5d639767687de88634f0d235e6526498 Pseudo is Ultrasoft + core correction, Zval = 7.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1213 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 PseudoPot. # 3 for Cs read from file: /users/gautes/Pseudo/Cs.rel-pbe-oncvpsp.UPF MD5 check sum: b26d75608c6aa7e713949ad9f7a49707 Pseudo is Norm-conserving, Zval = 9.0 Generated using ONCVPSP code by D. R. Hamann Using radial grid of 2414 points, 14 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 l(11) = 3 l(12) = 3 l(13) = 3 l(14) = 3 atomic species valence mass pseudopotential Hg 12.00 200.59000 Hg( 1.00) Br 7.00 79.90400 Br( 1.00) Cs 9.00 132.90550 Cs( 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.0156250 k( 2) = ( 0.0000000 0.0000000 0.2500000), wk = 0.0937500 k( 3) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0468750 k( 4) = ( 0.0000000 0.2500000 0.2500000), wk = 0.1875000 k( 5) = ( 0.0000000 0.2500000 -0.5000000), wk = 0.1875000 k( 6) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0468750 k( 7) = ( 0.2500000 0.2500000 0.2500000), wk = 0.1250000 k( 8) = ( 0.2500000 0.2500000 -0.5000000), wk = 0.1875000 k( 9) = ( 0.2500000 -0.5000000 -0.5000000), wk = 0.0937500 k( 10) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0156250 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0156250 k( 2) = ( 0.0000000 0.0000000 0.2500000), wk = 0.0937500 k( 3) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0468750 k( 4) = ( 0.0000000 0.2500000 0.2500000), wk = 0.1875000 k( 5) = ( 0.0000000 0.2500000 -0.5000000), wk = 0.1875000 k( 6) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0468750 k( 7) = ( 0.2500000 0.2500000 0.2500000), wk = 0.1250000 k( 8) = ( 0.2500000 0.2500000 -0.5000000), wk = 0.1875000 k( 9) = ( 0.2500000 -0.5000000 -0.5000000), wk = 0.0937500 k( 10) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0156250 Dense grid: 185639 G-vectors FFT dimensions: ( 72, 72, 72) Smooth grid: 125039 G-vectors FFT dimensions: ( 64, 64, 64) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.68 Mb ( 888, 50) NL pseudopotentials 0.93 Mb ( 444, 138) Each V/rho on FFT grid 0.16 Mb ( 10368) Each G-vector array 0.04 Mb ( 5158) G-vector shells 0.01 Mb ( 895) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 2.71 Mb ( 888, 200) Each subspace H/S matrix 0.04 Mb ( 50, 50) Each matrix 0.21 Mb ( 138, 2, 50) Arrays for rho mixing 1.27 Mb ( 10368, 8) Initial potential from superposition of free atoms starting charge 41.96166, renormalised to 42.00000 Starting wfc are 52 randomized atomic wfcs total cpu time spent up to now is 3.3 secs per-process dynamical memory: 87.9 Mb Self-consistent Calculation iteration # 1 ecut= 80.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 = 5.95E-04, avg # of iterations = 3.8 total cpu time spent up to now is 8.5 secs total energy = -228.84650158 Ry Harris-Foulkes estimate = -229.04170296 Ry estimated scf accuracy < 0.25622516 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.10E-04, avg # of iterations = 3.8 total cpu time spent up to now is 11.6 secs total energy = -228.91763391 Ry Harris-Foulkes estimate = -228.98805162 Ry estimated scf accuracy < 0.11664334 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.78E-04, avg # of iterations = 2.0 total cpu time spent up to now is 14.2 secs total energy = -228.94901381 Ry Harris-Foulkes estimate = -228.94950217 Ry estimated scf accuracy < 0.00164337 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.91E-06, avg # of iterations = 6.3 total cpu time spent up to now is 18.2 secs total energy = -228.94990590 Ry Harris-Foulkes estimate = -228.94992239 Ry estimated scf accuracy < 0.00010560 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.51E-07, avg # of iterations = 3.1 total cpu time spent up to now is 21.2 secs total energy = -228.94993590 Ry Harris-Foulkes estimate = -228.94993119 Ry estimated scf accuracy < 0.00000813 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.93E-08, avg # of iterations = 3.6 total cpu time spent up to now is 24.3 secs total energy = -228.94993893 Ry Harris-Foulkes estimate = -228.94993849 Ry estimated scf accuracy < 0.00000039 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.34E-10, avg # of iterations = 3.7 total cpu time spent up to now is 27.4 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 15515 PWs) bands (ev): -16.1290 -16.1290 -10.9098 -10.9098 -10.5532 -10.5532 -10.5532 -10.5532 -4.8724 -4.8724 -3.8238 -3.8238 -3.8238 -3.8238 -3.4261 -3.4261 -3.4261 -3.4261 -2.0490 -2.0490 -1.9617 -1.9617 -1.9617 -1.9617 -0.0516 -0.0516 0.0131 0.0131 0.0131 0.0131 1.4265 1.4265 1.5614 1.5614 1.6609 1.6609 1.6609 1.6609 1.8769 1.8769 2.1023 2.1023 2.1023 2.1023 7.2593 7.2593 7.9814 7.9814 7.9814 7.9814 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.2500 ( 15675 PWs) bands (ev): -16.1278 -16.1278 -10.8621 -10.8621 -10.5732 -10.5732 -10.5520 -10.5520 -4.8674 -4.8674 -3.9266 -3.9266 -3.8275 -3.8275 -3.4164 -3.4164 -3.4135 -3.4135 -2.2655 -2.2655 -2.0654 -2.0654 -1.9641 -1.9641 -0.2231 -0.2231 0.0393 0.0393 0.0448 0.0448 1.2025 1.2025 1.5840 1.5840 1.6229 1.6229 1.7519 1.7519 1.8423 1.8423 2.0350 2.0350 3.2261 3.2261 7.8610 7.8610 8.4041 8.4041 8.4487 8.4488 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.0000-0.5000 ( 15624 PWs) bands (ev): -16.1265 -16.1265 -10.8053 -10.8053 -10.6017 -10.6017 -10.5509 -10.5509 -4.8625 -4.8625 -4.0450 -4.0450 -3.8311 -3.8311 -3.4066 -3.4066 -3.4009 -3.4009 -2.4963 -2.4963 -2.0820 -2.0820 -1.9657 -1.9657 -0.3035 -0.3035 0.0731 0.0731 0.1074 0.1074 1.0767 1.0767 1.5243 1.5243 1.5796 1.5796 1.7006 1.7006 1.8716 1.8716 1.9172 1.9172 4.6664 4.6664 8.5582 8.5582 9.0294 9.0295 9.0487 9.0488 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.2500 0.2500 ( 15680 PWs) bands (ev): -16.1266 -16.1266 -10.8131 -10.8131 -10.5790 -10.5790 -10.5631 -10.5631 -4.8673 -4.8673 -3.9678 -3.9678 -3.8850 -3.8850 -3.4239 -3.4239 -3.4044 -3.4044 -2.3153 -2.3153 -2.1978 -2.1978 -2.0851 -2.0851 -0.9117 -0.9117 0.1112 0.1112 0.9592 0.9592 1.1003 1.1003 1.5140 1.5140 1.6347 1.6347 1.7550 1.7550 1.8926 1.8926 2.1890 2.1890 3.8126 3.8126 8.3058 8.3058 8.6805 8.6805 8.8256 8.8256 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.2500-0.5000 ( 15656 PWs) bands (ev): -16.1255 -16.1255 -10.7524 -10.7524 -10.6031 -10.6031 -10.5667 -10.5667 -4.8671 -4.8671 -4.0635 -4.0635 -3.9050 -3.9050 -3.4243 -3.4243 -3.4018 -3.4018 -2.5090 -2.5090 -2.2543 -2.2543 -2.1043 -2.1043 -1.0666 -1.0666 0.1836 0.1836 1.0207 1.0207 1.2077 1.2077 1.5004 1.5004 1.6226 1.6226 1.6946 1.6946 2.0689 2.0689 2.1365 2.1365 5.0066 5.0066 8.4948 8.4948 9.1451 9.1452 9.2832 9.2833 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.5000-0.5000 ( 15676 PWs) bands (ev): -16.1246 -16.1246 -10.6773 -10.6773 -10.6090 -10.6090 -10.5994 -10.5994 -4.8715 -4.8715 -4.1093 -4.1093 -3.9703 -3.9703 -3.4423 -3.4423 -3.4006 -3.4006 -2.5236 -2.5236 -2.5176 -2.5176 -2.1290 -2.1290 -1.2325 -1.2325 0.2869 0.2869 0.9505 0.9505 1.3762 1.3762 1.5475 1.5475 1.6330 1.6330 1.6448 1.6448 2.2516 2.2516 2.6854 2.6854 5.8881 5.8881 7.8539 7.8539 8.6805 8.6805 9.1583 9.1583 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.9688 0.9688 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2500 0.2500 0.2500 ( 15675 PWs) bands (ev): -16.1257 -16.1257 -10.7655 -10.7655 -10.5794 -10.5794 -10.5744 -10.5744 -4.8735 -4.8735 -3.9720 -3.9720 -3.9719 -3.9719 -3.4404 -3.4404 -3.4070 -3.4070 -2.3208 -2.3208 -2.3160 -2.3160 -2.1063 -2.1063 -1.6398 -1.6398 1.0171 1.0171 1.0458 1.0458 1.1034 1.1034 1.4889 1.4889 1.5287 1.5287 1.9252 1.9252 2.0309 2.0309 2.2542 2.2542 4.2383 4.2383 8.4031 8.4031 8.7849 8.7849 8.8182 8.8182 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.2500 0.2500-0.5000 ( 15674 PWs) bands (ev): -16.1247 -16.1247 -10.7054 -10.7054 -10.6026 -10.6026 -10.5745 -10.5745 -4.8795 -4.8795 -4.0754 -4.0754 -3.9763 -3.9763 -3.4463 -3.4463 -3.4220 -3.4220 -2.5131 -2.5131 -2.3217 -2.3217 -2.1258 -2.1258 -1.8615 -1.8615 1.0509 1.0509 1.1209 1.1209 1.3070 1.3070 1.4050 1.4050 1.8251 1.8251 1.9841 1.9841 2.1316 2.1316 2.2052 2.2052 5.2845 5.2845 8.0339 8.0339 8.8018 8.8018 8.9351 8.9351 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.2500-0.5000-0.5000 ( 15644 PWs) bands (ev): -16.1238 -16.1238 -10.6394 -10.6394 -10.6057 -10.6057 -10.5967 -10.5967 -4.8917 -4.8917 -4.1145 -4.1145 -4.0374 -4.0374 -3.4649 -3.4649 -3.4474 -3.4474 -2.5215 -2.5215 -2.5150 -2.5150 -2.1451 -2.1451 -2.0975 -2.0975 1.0586 1.0586 1.1907 1.1907 1.3979 1.3979 1.6481 1.6481 1.8646 1.8646 2.0590 2.0590 2.2867 2.2867 2.6967 2.6967 6.1120 6.1120 7.4283 7.4283 8.0596 8.0596 8.5190 8.5190 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.9312 0.9312 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k =-0.5000-0.5000-0.5000 ( 15624 PWs) bands (ev): -16.1230 -16.1230 -10.6088 -10.6088 -10.5939 -10.5939 -10.5939 -10.5939 -4.9118 -4.9118 -4.1202 -4.1202 -4.1202 -4.1202 -3.4909 -3.4909 -3.4909 -3.4909 -2.5646 -2.5646 -2.5085 -2.5085 -2.5085 -2.5085 -2.1602 -2.1602 1.1496 1.1496 1.1496 1.1496 1.3593 1.3593 2.0478 2.0478 2.3002 2.3002 2.3002 2.3002 2.7170 2.7170 2.7170 2.7170 6.7869 6.7869 7.0081 7.0081 7.5755 7.5755 7.5755 7.5755 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.7531 0.7531 0.7531 0.7531 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 2.7321 ev ! total energy = -228.94993914 Ry Harris-Foulkes estimate = -228.94993912 Ry estimated scf accuracy < 4.6E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -62.53281546 Ry hartree contribution = 50.88693023 Ry xc contribution = -64.28303197 Ry ewald contribution = -153.02092703 Ry smearing contrib. (-TS) = -0.00009491 Ry convergence has been achieved in 7 iterations Writing output data file CsHgBr3.save init_run : 1.32s CPU 1.68s WALL ( 1 calls) electrons : 22.54s CPU 24.37s WALL ( 1 calls) Called by init_run: wfcinit : 0.68s CPU 0.81s WALL ( 1 calls) potinit : 0.04s CPU 0.04s WALL ( 1 calls) Called by electrons: c_bands : 18.39s CPU 18.98s WALL ( 8 calls) sum_band : 2.86s CPU 3.25s WALL ( 8 calls) v_of_rho : 0.06s CPU 0.06s WALL ( 8 calls) v_h : 0.01s CPU 0.00s WALL ( 8 calls) v_xc : 0.06s CPU 0.06s WALL ( 8 calls) newd : 1.19s CPU 1.64s WALL ( 8 calls) mix_rho : 0.05s CPU 0.05s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.10s CPU 0.12s WALL ( 170 calls) cegterg : 16.70s CPU 16.87s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.12s CPU 0.13s WALL ( 80 calls) addusdens : 0.47s CPU 0.82s WALL ( 8 calls) Called by *egterg: h_psi : 11.97s CPU 12.11s WALL ( 405 calls) s_psi : 0.50s CPU 0.48s WALL ( 405 calls) g_psi : 0.03s CPU 0.04s WALL ( 315 calls) cdiaghg : 2.27s CPU 2.32s WALL ( 385 calls) cegterg:over : 0.77s CPU 0.78s WALL ( 315 calls) cegterg:upda : 0.86s CPU 0.81s WALL ( 315 calls) cegterg:last : 0.20s CPU 0.21s WALL ( 80 calls) cdiaghg:chol : 0.16s CPU 0.14s WALL ( 385 calls) cdiaghg:inve : 0.08s CPU 0.09s WALL ( 385 calls) cdiaghg:para : 0.15s CPU 0.14s WALL ( 770 calls) Called by h_psi: h_psi:vloc : 10.63s CPU 10.72s WALL ( 405 calls) h_psi:vnl : 1.28s CPU 1.33s WALL ( 405 calls) add_vuspsi : 0.68s CPU 0.69s WALL ( 405 calls) General routines calbec : 0.81s CPU 0.84s WALL ( 485 calls) fft : 0.21s CPU 0.20s WALL ( 242 calls) ffts : 0.02s CPU 0.03s WALL ( 64 calls) fftw : 11.90s CPU 11.88s WALL ( 55016 calls) interpolate : 0.06s CPU 0.06s WALL ( 64 calls) Parallel routines fft_scatter : 4.02s CPU 4.08s WALL ( 55322 calls) PWSCF : 26.54s CPU 31.28s WALL This run was terminated on: 15:57:16 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=