Program EPW v.2.3.2 starts ... Today is 29Sep2010 at 14:12:54 Ultrasoft (Vanderbilt) Pseudopotentials bravais-lattice index = 2 lattice parameter (a_0) = 9.2000 a.u. unit-cell volume = 194.6720 (a.u.)^3 number of atoms/cell = 1 number of atomic types = 1 kinetic-energy cut-off = 40.0000 Ry charge density cut-off = 160.0000 Ry convergence threshold = 0.0E+00 beta = 0.0000 number of iterations used = 0 Exchange-correlation = SLA PZ NOGX NOGC (1100) celldm(1)= 9.20000 celldm(2)= 0.00000 celldm(3)= 0.00000 celldm(4)= 0.00000 celldm(5)= 0.00000 celldm(6)= 0.00000 crystal axes: (cart. coord. in units of a_0) a(1) = ( -0.5000 0.0000 0.5000 ) a(2) = ( 0.0000 0.5000 0.5000 ) a(3) = ( -0.5000 0.5000 0.0000 ) reciprocal axes: (cart. coord. in units 2 pi/a_0) b(1) = ( -1.0000 -1.0000 1.0000 ) b(2) = ( 1.0000 1.0000 1.0000 ) b(3) = ( -1.0000 1.0000 -1.0000 ) Atoms inside the unit cell: Cartesian axes site n. atom mass positions (a_0 units) 1 Pb 207.2175 tau( 1) = ( 0.00000 0.00000 0.00000 ) 49 Sym.Ops. (with q -> -q+G ) G cutoff = 343.0330 ( 6615 G-vectors) FFT grid: ( 27, 27, 27) number of k points= 216 gaussian broad. (Ry)= 0.0500 ngauss = 0 cart. coord. in units 2pi/a_0 k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0092593 k( 2) = ( -0.1666667 0.1666667 -0.1666667), wk = 0.0092593 k( 3) = ( -0.3333333 0.3333333 -0.3333333), wk = 0.0092593 k( 4) = ( -0.5000000 0.5000000 -0.5000000), wk = 0.0092593 k( 5) = ( -0.6666667 0.6666667 -0.6666667), wk = 0.0092593 k( 6) = ( -0.8333333 0.8333333 -0.8333333), wk = 0.0092593 k( 7) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0092593 k( 8) = ( 0.0000000 0.3333333 0.0000000), wk = 0.0092593 k( 9) = ( -0.1666667 0.5000000 -0.1666667), wk = 0.0092593 k( 10) = ( -0.3333333 0.6666667 -0.3333333), wk = 0.0092593 k( 11) = ( -0.5000000 0.8333333 -0.5000000), wk = 0.0092593 k( 12) = ( -0.6666667 1.0000000 -0.6666667), wk = 0.0092593 k( 13) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0092593 k( 14) = ( 0.1666667 0.5000000 0.1666667), wk = 0.0092593 k( 15) = ( 0.0000000 0.6666667 0.0000000), wk = 0.0092593 k( 16) = ( -0.1666667 0.8333333 -0.1666667), wk = 0.0092593 k( 17) = ( -0.3333333 1.0000000 -0.3333333), wk = 0.0092593 k( 18) = ( -0.5000000 1.1666667 -0.5000000), wk = 0.0092593 k( 19) = ( 0.5000000 0.5000000 0.5000000), wk = 0.0092593 k( 20) = ( 0.3333333 0.6666667 0.3333333), wk = 0.0092593 k( 21) = ( 0.1666667 0.8333333 0.1666667), wk = 0.0092593 k( 22) = ( 0.0000000 1.0000000 0.0000000), wk = 0.0092593 k( 23) = ( -0.1666667 1.1666667 -0.1666667), wk = 0.0092593 k( 24) = ( -0.3333333 1.3333333 -0.3333333), wk = 0.0092593 k( 25) = ( 0.6666667 0.6666667 0.6666667), wk = 0.0092593 k( 26) = ( 0.5000000 0.8333333 0.5000000), wk = 0.0092593 k( 27) = ( 0.3333333 1.0000000 0.3333333), wk = 0.0092593 k( 28) = ( 0.1666667 1.1666667 0.1666667), wk = 0.0092593 k( 29) = ( 0.0000000 1.3333333 0.0000000), wk = 0.0092593 k( 30) = ( -0.1666667 1.5000000 -0.1666667), wk = 0.0092593 k( 31) = ( 0.8333333 0.8333333 0.8333333), wk = 0.0092593 k( 32) = ( 0.6666667 1.0000000 0.6666667), wk = 0.0092593 k( 33) = ( 0.5000000 1.1666667 0.5000000), wk = 0.0092593 k( 34) = ( 0.3333333 1.3333333 0.3333333), wk = 0.0092593 k( 35) = ( 0.1666667 1.5000000 0.1666667), wk = 0.0092593 k( 36) = ( 0.0000000 1.6666667 0.0000000), wk = 0.0092593 k( 37) = ( -0.1666667 -0.1666667 0.1666667), wk = 0.0092593 k( 38) = ( -0.3333333 0.0000000 0.0000000), wk = 0.0092593 k( 39) = ( -0.5000000 0.1666667 -0.1666667), wk = 0.0092593 k( 40) = ( -0.6666667 0.3333333 -0.3333333), wk = 0.0092593 k( 41) = ( -0.8333333 0.5000000 -0.5000000), wk = 0.0092593 k( 42) = ( -1.0000000 0.6666667 -0.6666667), wk = 0.0092593 k( 43) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0092593 k( 44) = ( -0.1666667 0.1666667 0.1666667), wk = 0.0092593 k( 45) = ( -0.3333333 0.3333333 0.0000000), wk = 0.0092593 k( 46) = ( -0.5000000 0.5000000 -0.1666667), wk = 0.0092593 k( 47) = ( -0.6666667 0.6666667 -0.3333333), wk = 0.0092593 k( 48) = ( -0.8333333 0.8333333 -0.5000000), wk = 0.0092593 k( 49) = ( 0.1666667 0.1666667 0.5000000), wk = 0.0092593 k( 50) = ( -0.0000000 0.3333333 0.3333333), wk = 0.0092593 k( 51) = ( -0.1666667 0.5000000 0.1666667), wk = 0.0092593 k( 52) = ( -0.3333333 0.6666667 -0.0000000), wk = 0.0092593 k( 53) = ( -0.5000000 0.8333333 -0.1666667), wk = 0.0092593 k( 54) = ( -0.6666667 1.0000000 -0.3333333), wk = 0.0092593 k( 55) = ( 0.3333333 0.3333333 0.6666667), wk = 0.0092593 k( 56) = ( 0.1666667 0.5000000 0.5000000), wk = 0.0092593 k( 57) = ( 0.0000000 0.6666667 0.3333333), wk = 0.0092593 k( 58) = ( -0.1666667 0.8333333 0.1666667), wk = 0.0092593 k( 59) = ( -0.3333333 1.0000000 -0.0000000), wk = 0.0092593 k( 60) = ( -0.5000000 1.1666667 -0.1666667), wk = 0.0092593 k( 61) = ( 0.5000000 0.5000000 0.8333333), wk = 0.0092593 k( 62) = ( 0.3333333 0.6666667 0.6666667), wk = 0.0092593 k( 63) = ( 0.1666667 0.8333333 0.5000000), wk = 0.0092593 k( 64) = ( 0.0000000 1.0000000 0.3333333), wk = 0.0092593 k( 65) = ( -0.1666667 1.1666667 0.1666667), wk = 0.0092593 k( 66) = ( -0.3333333 1.3333333 0.0000000), wk = 0.0092593 k( 67) = ( 0.6666667 0.6666667 1.0000000), wk = 0.0092593 k( 68) = ( 0.5000000 0.8333333 0.8333333), wk = 0.0092593 k( 69) = ( 0.3333333 1.0000000 0.6666667), wk = 0.0092593 k( 70) = ( 0.1666667 1.1666667 0.5000000), wk = 0.0092593 k( 71) = ( -0.0000000 1.3333333 0.3333333), wk = 0.0092593 k( 72) = ( -0.1666667 1.5000000 0.1666667), wk = 0.0092593 k( 73) = ( -0.3333333 -0.3333333 0.3333333), wk = 0.0092593 k( 74) = ( -0.5000000 -0.1666667 0.1666667), wk = 0.0092593 k( 75) = ( -0.6666667 0.0000000 0.0000000), wk = 0.0092593 k( 76) = ( -0.8333333 0.1666667 -0.1666667), wk = 0.0092593 k( 77) = ( -1.0000000 0.3333333 -0.3333333), wk = 0.0092593 k( 78) = ( -1.1666667 0.5000000 -0.5000000), wk = 0.0092593 k( 79) = ( -0.1666667 -0.1666667 0.5000000), wk = 0.0092593 k( 80) = ( -0.3333333 0.0000000 0.3333333), wk = 0.0092593 k( 81) = ( -0.5000000 0.1666667 0.1666667), wk = 0.0092593 k( 82) = ( -0.6666667 0.3333333 -0.0000000), wk = 0.0092593 k( 83) = ( -0.8333333 0.5000000 -0.1666667), wk = 0.0092593 k( 84) = ( -1.0000000 0.6666667 -0.3333333), wk = 0.0092593 k( 85) = ( 0.0000000 0.0000000 0.6666667), wk = 0.0092593 k( 86) = ( -0.1666667 0.1666667 0.5000000), wk = 0.0092593 k( 87) = ( -0.3333333 0.3333333 0.3333333), wk = 0.0092593 k( 88) = ( -0.5000000 0.5000000 0.1666667), wk = 0.0092593 k( 89) = ( -0.6666667 0.6666667 -0.0000000), wk = 0.0092593 k( 90) = ( -0.8333333 0.8333333 -0.1666667), wk = 0.0092593 k( 91) = ( 0.1666667 0.1666667 0.8333333), wk = 0.0092593 k( 92) = ( 0.0000000 0.3333333 0.6666667), wk = 0.0092593 k( 93) = ( -0.1666667 0.5000000 0.5000000), wk = 0.0092593 k( 94) = ( -0.3333333 0.6666667 0.3333333), wk = 0.0092593 k( 95) = ( -0.5000000 0.8333333 0.1666667), wk = 0.0092593 k( 96) = ( -0.6666667 1.0000000 0.0000000), wk = 0.0092593 k( 97) = ( 0.3333333 0.3333333 1.0000000), wk = 0.0092593 k( 98) = ( 0.1666667 0.5000000 0.8333333), wk = 0.0092593 k( 99) = ( 0.0000000 0.6666667 0.6666667), wk = 0.0092593 k( 100) = ( -0.1666667 0.8333333 0.5000000), wk = 0.0092593 k( 101) = ( -0.3333333 1.0000000 0.3333333), wk = 0.0092593 k( 102) = ( -0.5000000 1.1666667 0.1666667), wk = 0.0092593 k( 103) = ( 0.5000000 0.5000000 1.1666667), wk = 0.0092593 k( 104) = ( 0.3333333 0.6666667 1.0000000), wk = 0.0092593 k( 105) = ( 0.1666667 0.8333333 0.8333333), wk = 0.0092593 k( 106) = ( 0.0000000 1.0000000 0.6666667), wk = 0.0092593 k( 107) = ( -0.1666667 1.1666667 0.5000000), wk = 0.0092593 k( 108) = ( -0.3333333 1.3333333 0.3333333), wk = 0.0092593 k( 109) = ( -0.5000000 -0.5000000 0.5000000), wk = 0.0092593 k( 110) = ( -0.6666667 -0.3333333 0.3333333), wk = 0.0092593 k( 111) = ( -0.8333333 -0.1666667 0.1666667), wk = 0.0092593 k( 112) = ( -1.0000000 0.0000000 0.0000000), wk = 0.0092593 k( 113) = ( -1.1666667 0.1666667 -0.1666667), wk = 0.0092593 k( 114) = ( -1.3333333 0.3333333 -0.3333333), wk = 0.0092593 k( 115) = ( -0.3333333 -0.3333333 0.6666667), wk = 0.0092593 k( 116) = ( -0.5000000 -0.1666667 0.5000000), wk = 0.0092593 k( 117) = ( -0.6666667 -0.0000000 0.3333333), wk = 0.0092593 k( 118) = ( -0.8333333 0.1666667 0.1666667), wk = 0.0092593 k( 119) = ( -1.0000000 0.3333333 -0.0000000), wk = 0.0092593 k( 120) = ( -1.1666667 0.5000000 -0.1666667), wk = 0.0092593 k( 121) = ( -0.1666667 -0.1666667 0.8333333), wk = 0.0092593 k( 122) = ( -0.3333333 -0.0000000 0.6666667), wk = 0.0092593 k( 123) = ( -0.5000000 0.1666667 0.5000000), wk = 0.0092593 k( 124) = ( -0.6666667 0.3333333 0.3333333), wk = 0.0092593 k( 125) = ( -0.8333333 0.5000000 0.1666667), wk = 0.0092593 k( 126) = ( -1.0000000 0.6666667 0.0000000), wk = 0.0092593 k( 127) = ( 0.0000000 0.0000000 1.0000000), wk = 0.0092593 k( 128) = ( -0.1666667 0.1666667 0.8333333), wk = 0.0092593 k( 129) = ( -0.3333333 0.3333333 0.6666667), wk = 0.0092593 k( 130) = ( -0.5000000 0.5000000 0.5000000), wk = 0.0092593 k( 131) = ( -0.6666667 0.6666667 0.3333333), wk = 0.0092593 k( 132) = ( -0.8333333 0.8333333 0.1666667), wk = 0.0092593 k( 133) = ( 0.1666667 0.1666667 1.1666667), wk = 0.0092593 k( 134) = ( 0.0000000 0.3333333 1.0000000), wk = 0.0092593 k( 135) = ( -0.1666667 0.5000000 0.8333333), wk = 0.0092593 k( 136) = ( -0.3333333 0.6666667 0.6666667), wk = 0.0092593 k( 137) = ( -0.5000000 0.8333333 0.5000000), wk = 0.0092593 k( 138) = ( -0.6666667 1.0000000 0.3333333), wk = 0.0092593 k( 139) = ( 0.3333333 0.3333333 1.3333333), wk = 0.0092593 k( 140) = ( 0.1666667 0.5000000 1.1666667), wk = 0.0092593 k( 141) = ( 0.0000000 0.6666667 1.0000000), wk = 0.0092593 k( 142) = ( -0.1666667 0.8333333 0.8333333), wk = 0.0092593 k( 143) = ( -0.3333333 1.0000000 0.6666667), wk = 0.0092593 k( 144) = ( -0.5000000 1.1666667 0.5000000), wk = 0.0092593 k( 145) = ( -0.6666667 -0.6666667 0.6666667), wk = 0.0092593 k( 146) = ( -0.8333333 -0.5000000 0.5000000), wk = 0.0092593 k( 147) = ( -1.0000000 -0.3333333 0.3333333), wk = 0.0092593 k( 148) = ( -1.1666667 -0.1666667 0.1666667), wk = 0.0092593 k( 149) = ( -1.3333333 0.0000000 0.0000000), wk = 0.0092593 k( 150) = ( -1.5000000 0.1666667 -0.1666667), wk = 0.0092593 k( 151) = ( -0.5000000 -0.5000000 0.8333333), wk = 0.0092593 k( 152) = ( -0.6666667 -0.3333333 0.6666667), wk = 0.0092593 k( 153) = ( -0.8333333 -0.1666667 0.5000000), wk = 0.0092593 k( 154) = ( -1.0000000 -0.0000000 0.3333333), wk = 0.0092593 k( 155) = ( -1.1666667 0.1666667 0.1666667), wk = 0.0092593 k( 156) = ( -1.3333333 0.3333333 0.0000000), wk = 0.0092593 k( 157) = ( -0.3333333 -0.3333333 1.0000000), wk = 0.0092593 k( 158) = ( -0.5000000 -0.1666667 0.8333333), wk = 0.0092593 k( 159) = ( -0.6666667 -0.0000000 0.6666667), wk = 0.0092593 k( 160) = ( -0.8333333 0.1666667 0.5000000), wk = 0.0092593 k( 161) = ( -1.0000000 0.3333333 0.3333333), wk = 0.0092593 k( 162) = ( -1.1666667 0.5000000 0.1666667), wk = 0.0092593 k( 163) = ( -0.1666667 -0.1666667 1.1666667), wk = 0.0092593 k( 164) = ( -0.3333333 -0.0000000 1.0000000), wk = 0.0092593 k( 165) = ( -0.5000000 0.1666667 0.8333333), wk = 0.0092593 k( 166) = ( -0.6666667 0.3333333 0.6666667), wk = 0.0092593 k( 167) = ( -0.8333333 0.5000000 0.5000000), wk = 0.0092593 k( 168) = ( -1.0000000 0.6666667 0.3333333), wk = 0.0092593 k( 169) = ( 0.0000000 0.0000000 1.3333333), wk = 0.0092593 k( 170) = ( -0.1666667 0.1666667 1.1666667), wk = 0.0092593 k( 171) = ( -0.3333333 0.3333333 1.0000000), wk = 0.0092593 k( 172) = ( -0.5000000 0.5000000 0.8333333), wk = 0.0092593 k( 173) = ( -0.6666667 0.6666667 0.6666667), wk = 0.0092593 k( 174) = ( -0.8333333 0.8333333 0.5000000), wk = 0.0092593 k( 175) = ( 0.1666667 0.1666667 1.5000000), wk = 0.0092593 k( 176) = ( -0.0000000 0.3333333 1.3333333), wk = 0.0092593 k( 177) = ( -0.1666667 0.5000000 1.1666667), wk = 0.0092593 k( 178) = ( -0.3333333 0.6666667 1.0000000), wk = 0.0092593 k( 179) = ( -0.5000000 0.8333333 0.8333333), wk = 0.0092593 k( 180) = ( -0.6666667 1.0000000 0.6666667), wk = 0.0092593 k( 181) = ( -0.8333333 -0.8333333 0.8333333), wk = 0.0092593 k( 182) = ( -1.0000000 -0.6666667 0.6666667), wk = 0.0092593 k( 183) = ( -1.1666667 -0.5000000 0.5000000), wk = 0.0092593 k( 184) = ( -1.3333333 -0.3333333 0.3333333), wk = 0.0092593 k( 185) = ( -1.5000000 -0.1666667 0.1666667), wk = 0.0092593 k( 186) = ( -1.6666667 0.0000000 0.0000000), wk = 0.0092593 k( 187) = ( -0.6666667 -0.6666667 1.0000000), wk = 0.0092593 k( 188) = ( -0.8333333 -0.5000000 0.8333333), wk = 0.0092593 k( 189) = ( -1.0000000 -0.3333333 0.6666667), wk = 0.0092593 k( 190) = ( -1.1666667 -0.1666667 0.5000000), wk = 0.0092593 k( 191) = ( -1.3333333 0.0000000 0.3333333), wk = 0.0092593 k( 192) = ( -1.5000000 0.1666667 0.1666667), wk = 0.0092593 k( 193) = ( -0.5000000 -0.5000000 1.1666667), wk = 0.0092593 k( 194) = ( -0.6666667 -0.3333333 1.0000000), wk = 0.0092593 k( 195) = ( -0.8333333 -0.1666667 0.8333333), wk = 0.0092593 k( 196) = ( -1.0000000 0.0000000 0.6666667), wk = 0.0092593 k( 197) = ( -1.1666667 0.1666667 0.5000000), wk = 0.0092593 k( 198) = ( -1.3333333 0.3333333 0.3333333), wk = 0.0092593 k( 199) = ( -0.3333333 -0.3333333 1.3333333), wk = 0.0092593 k( 200) = ( -0.5000000 -0.1666667 1.1666667), wk = 0.0092593 k( 201) = ( -0.6666667 0.0000000 1.0000000), wk = 0.0092593 k( 202) = ( -0.8333333 0.1666667 0.8333333), wk = 0.0092593 k( 203) = ( -1.0000000 0.3333333 0.6666667), wk = 0.0092593 k( 204) = ( -1.1666667 0.5000000 0.5000000), wk = 0.0092593 k( 205) = ( -0.1666667 -0.1666667 1.5000000), wk = 0.0092593 k( 206) = ( -0.3333333 0.0000000 1.3333333), wk = 0.0092593 k( 207) = ( -0.5000000 0.1666667 1.1666667), wk = 0.0092593 k( 208) = ( -0.6666667 0.3333333 1.0000000), wk = 0.0092593 k( 209) = ( -0.8333333 0.5000000 0.8333333), wk = 0.0092593 k( 210) = ( -1.0000000 0.6666667 0.6666667), wk = 0.0092593 k( 211) = ( 0.0000000 0.0000000 1.6666667), wk = 0.0092593 k( 212) = ( -0.1666667 0.1666667 1.5000000), wk = 0.0092593 k( 213) = ( -0.3333333 0.3333333 1.3333333), wk = 0.0092593 k( 214) = ( -0.5000000 0.5000000 1.1666667), wk = 0.0092593 k( 215) = ( -0.6666667 0.6666667 1.0000000), wk = 0.0092593 k( 216) = ( -0.8333333 0.8333333 0.8333333), wk = 0.0092593 PseudoPot. # 1 for Pb read from file Pb.pz-bhs.UPF Pseudo is Norm-conserving, Zval = 4.0 Generated by new atomic code, or converted to UPF format Using radial grid of 331 points, 3 beta functions with: l(1) = 0 l(2) = 1 l(3) = 2 EPW : 2.68s CPU time, 4.22s wall time EPW : 3.30s CPU time, 5.16s wall time EPW : 3.30s CPU time, 5.16s wall time No wavefunction gauge setting applied ------------------------------------------------------------------- Wannierization on 6 x 6 x 6 electronic grid ------------------------------------------------------------------- Spin CASE ( default = unpolarized ) Initializing Wannier90 Initial Wannier projections ( 0.00000 0.00000 0.00000) : l = -3 mr = 1 ( 0.00000 0.00000 0.00000) : l = -3 mr = 2 ( 0.00000 0.00000 0.00000) : l = -3 mr = 3 ( 0.00000 0.00000 0.00000) : l = -3 mr = 4 - Number of bands is ( 6) - Number of wannier functions is ( 4) - All guiding functions are given - All neighbours are found AMN 1 of 216 on ionode 2 of 216 on ionode 3 of 216 on ionode 4 of 216 on ionode 5 of 216 on ionode 6 of 216 on ionode 7 of 216 on ionode 8 of 216 on ionode 9 of 216 on ionode 10 of 216 on ionode 11 of 216 on ionode 12 of 216 on ionode 13 of 216 on ionode 14 of 216 on ionode 15 of 216 on ionode 16 of 216 on ionode 17 of 216 on ionode 18 of 216 on ionode 19 of 216 on ionode 20 of 216 on ionode 21 of 216 on ionode 22 of 216 on ionode 23 of 216 on ionode 24 of 216 on ionode 25 of 216 on ionode 26 of 216 on ionode 27 of 216 on ionode 28 of 216 on ionode 29 of 216 on ionode 30 of 216 on ionode 31 of 216 on ionode 32 of 216 on ionode 33 of 216 on ionode 34 of 216 on ionode 35 of 216 on ionode 36 of 216 on ionode 37 of 216 on ionode 38 of 216 on ionode 39 of 216 on ionode 40 of 216 on ionode 41 of 216 on ionode 42 of 216 on ionode 43 of 216 on ionode 44 of 216 on ionode 45 of 216 on ionode 46 of 216 on ionode 47 of 216 on ionode 48 of 216 on ionode 49 of 216 on ionode 50 of 216 on ionode 51 of 216 on ionode 52 of 216 on ionode 53 of 216 on ionode 54 of 216 on ionode 55 of 216 on ionode 56 of 216 on ionode 57 of 216 on ionode 58 of 216 on ionode 59 of 216 on ionode 60 of 216 on ionode 61 of 216 on ionode 62 of 216 on ionode 63 of 216 on ionode 64 of 216 on ionode 65 of 216 on ionode 66 of 216 on ionode 67 of 216 on ionode 68 of 216 on ionode 69 of 216 on ionode 70 of 216 on ionode 71 of 216 on ionode 72 of 216 on ionode 73 of 216 on ionode 74 of 216 on ionode 75 of 216 on ionode 76 of 216 on ionode 77 of 216 on ionode 78 of 216 on ionode 79 of 216 on ionode 80 of 216 on ionode 81 of 216 on ionode 82 of 216 on ionode 83 of 216 on ionode 84 of 216 on ionode 85 of 216 on ionode 86 of 216 on ionode 87 of 216 on ionode 88 of 216 on ionode 89 of 216 on ionode 90 of 216 on ionode 91 of 216 on ionode 92 of 216 on ionode 93 of 216 on ionode 94 of 216 on ionode 95 of 216 on ionode 96 of 216 on ionode 97 of 216 on ionode 98 of 216 on ionode 99 of 216 on ionode 100 of 216 on ionode 101 of 216 on ionode 102 of 216 on ionode 103 of 216 on ionode 104 of 216 on ionode 105 of 216 on ionode 106 of 216 on ionode 107 of 216 on ionode 108 of 216 on ionode 109 of 216 on ionode 110 of 216 on ionode 111 of 216 on ionode 112 of 216 on ionode 113 of 216 on ionode 114 of 216 on ionode 115 of 216 on ionode 116 of 216 on ionode 117 of 216 on ionode 118 of 216 on ionode 119 of 216 on ionode 120 of 216 on ionode 121 of 216 on ionode 122 of 216 on ionode 123 of 216 on ionode 124 of 216 on ionode 125 of 216 on ionode 126 of 216 on ionode 127 of 216 on ionode 128 of 216 on ionode 129 of 216 on ionode 130 of 216 on ionode 131 of 216 on ionode 132 of 216 on ionode 133 of 216 on ionode 134 of 216 on ionode 135 of 216 on ionode 136 of 216 on ionode 137 of 216 on ionode 138 of 216 on ionode 139 of 216 on ionode 140 of 216 on ionode 141 of 216 on ionode 142 of 216 on ionode 143 of 216 on ionode 144 of 216 on ionode 145 of 216 on ionode 146 of 216 on ionode 147 of 216 on ionode 148 of 216 on ionode 149 of 216 on ionode 150 of 216 on ionode 151 of 216 on ionode 152 of 216 on ionode 153 of 216 on ionode 154 of 216 on ionode 155 of 216 on ionode 156 of 216 on ionode 157 of 216 on ionode 158 of 216 on ionode 159 of 216 on ionode 160 of 216 on ionode 161 of 216 on ionode 162 of 216 on ionode 163 of 216 on ionode 164 of 216 on ionode 165 of 216 on ionode 166 of 216 on ionode 167 of 216 on ionode 168 of 216 on ionode 169 of 216 on ionode 170 of 216 on ionode 171 of 216 on ionode 172 of 216 on ionode 173 of 216 on ionode 174 of 216 on ionode 175 of 216 on ionode 176 of 216 on ionode 177 of 216 on ionode 178 of 216 on ionode 179 of 216 on ionode 180 of 216 on ionode 181 of 216 on ionode 182 of 216 on ionode 183 of 216 on ionode 184 of 216 on ionode 185 of 216 on ionode 186 of 216 on ionode 187 of 216 on ionode 188 of 216 on ionode 189 of 216 on ionode 190 of 216 on ionode 191 of 216 on ionode 192 of 216 on ionode 193 of 216 on ionode 194 of 216 on ionode 195 of 216 on ionode 196 of 216 on ionode 197 of 216 on ionode 198 of 216 on ionode 199 of 216 on ionode 200 of 216 on ionode 201 of 216 on ionode 202 of 216 on ionode 203 of 216 on ionode 204 of 216 on ionode 205 of 216 on ionode 206 of 216 on ionode 207 of 216 on ionode 208 of 216 on ionode 209 of 216 on ionode 210 of 216 on ionode 211 of 216 on ionode 212 of 216 on ionode 213 of 216 on ionode 214 of 216 on ionode 215 of 216 on ionode 216 of 216 on ionode AMN calculated MMN 1 of 216 on ionode 2 of 216 on ionode 3 of 216 on ionode 4 of 216 on ionode 5 of 216 on ionode 6 of 216 on ionode 7 of 216 on ionode 8 of 216 on ionode 9 of 216 on ionode 10 of 216 on ionode 11 of 216 on ionode 12 of 216 on ionode 13 of 216 on ionode 14 of 216 on ionode 15 of 216 on ionode 16 of 216 on ionode 17 of 216 on ionode 18 of 216 on ionode 19 of 216 on ionode 20 of 216 on ionode 21 of 216 on ionode 22 of 216 on ionode 23 of 216 on ionode 24 of 216 on ionode 25 of 216 on ionode 26 of 216 on ionode 27 of 216 on ionode 28 of 216 on ionode 29 of 216 on ionode 30 of 216 on ionode 31 of 216 on ionode 32 of 216 on ionode 33 of 216 on ionode 34 of 216 on ionode 35 of 216 on ionode 36 of 216 on ionode 37 of 216 on ionode 38 of 216 on ionode 39 of 216 on ionode 40 of 216 on ionode 41 of 216 on ionode 42 of 216 on ionode 43 of 216 on ionode 44 of 216 on ionode 45 of 216 on ionode 46 of 216 on ionode 47 of 216 on ionode 48 of 216 on ionode 49 of 216 on ionode 50 of 216 on ionode 51 of 216 on ionode 52 of 216 on ionode 53 of 216 on ionode 54 of 216 on ionode 55 of 216 on ionode 56 of 216 on ionode 57 of 216 on ionode 58 of 216 on ionode 59 of 216 on ionode 60 of 216 on ionode 61 of 216 on ionode 62 of 216 on ionode 63 of 216 on ionode 64 of 216 on ionode 65 of 216 on ionode 66 of 216 on ionode 67 of 216 on ionode 68 of 216 on ionode 69 of 216 on ionode 70 of 216 on ionode 71 of 216 on ionode 72 of 216 on ionode 73 of 216 on ionode 74 of 216 on ionode 75 of 216 on ionode 76 of 216 on ionode 77 of 216 on ionode 78 of 216 on ionode 79 of 216 on ionode 80 of 216 on ionode 81 of 216 on ionode 82 of 216 on ionode 83 of 216 on ionode 84 of 216 on ionode 85 of 216 on ionode 86 of 216 on ionode 87 of 216 on ionode 88 of 216 on ionode 89 of 216 on ionode 90 of 216 on ionode 91 of 216 on ionode 92 of 216 on ionode 93 of 216 on ionode 94 of 216 on ionode 95 of 216 on ionode 96 of 216 on ionode 97 of 216 on ionode 98 of 216 on ionode 99 of 216 on ionode 100 of 216 on ionode 101 of 216 on ionode 102 of 216 on ionode 103 of 216 on ionode 104 of 216 on ionode 105 of 216 on ionode 106 of 216 on ionode 107 of 216 on ionode 108 of 216 on ionode 109 of 216 on ionode 110 of 216 on ionode 111 of 216 on ionode 112 of 216 on ionode 113 of 216 on ionode 114 of 216 on ionode 115 of 216 on ionode 116 of 216 on ionode 117 of 216 on ionode 118 of 216 on ionode 119 of 216 on ionode 120 of 216 on ionode 121 of 216 on ionode 122 of 216 on ionode 123 of 216 on ionode 124 of 216 on ionode 125 of 216 on ionode 126 of 216 on ionode 127 of 216 on ionode 128 of 216 on ionode 129 of 216 on ionode 130 of 216 on ionode 131 of 216 on ionode 132 of 216 on ionode 133 of 216 on ionode 134 of 216 on ionode 135 of 216 on ionode 136 of 216 on ionode 137 of 216 on ionode 138 of 216 on ionode 139 of 216 on ionode 140 of 216 on ionode 141 of 216 on ionode 142 of 216 on ionode 143 of 216 on ionode 144 of 216 on ionode 145 of 216 on ionode 146 of 216 on ionode 147 of 216 on ionode 148 of 216 on ionode 149 of 216 on ionode 150 of 216 on ionode 151 of 216 on ionode 152 of 216 on ionode 153 of 216 on ionode 154 of 216 on ionode 155 of 216 on ionode 156 of 216 on ionode 157 of 216 on ionode 158 of 216 on ionode 159 of 216 on ionode 160 of 216 on ionode 161 of 216 on ionode 162 of 216 on ionode 163 of 216 on ionode 164 of 216 on ionode 165 of 216 on ionode 166 of 216 on ionode 167 of 216 on ionode 168 of 216 on ionode 169 of 216 on ionode 170 of 216 on ionode 171 of 216 on ionode 172 of 216 on ionode 173 of 216 on ionode 174 of 216 on ionode 175 of 216 on ionode 176 of 216 on ionode 177 of 216 on ionode 178 of 216 on ionode 179 of 216 on ionode 180 of 216 on ionode 181 of 216 on ionode 182 of 216 on ionode 183 of 216 on ionode 184 of 216 on ionode 185 of 216 on ionode 186 of 216 on ionode 187 of 216 on ionode 188 of 216 on ionode 189 of 216 on ionode 190 of 216 on ionode 191 of 216 on ionode 192 of 216 on ionode 193 of 216 on ionode 194 of 216 on ionode 195 of 216 on ionode 196 of 216 on ionode 197 of 216 on ionode 198 of 216 on ionode 199 of 216 on ionode 200 of 216 on ionode 201 of 216 on ionode 202 of 216 on ionode 203 of 216 on ionode 204 of 216 on ionode 205 of 216 on ionode 206 of 216 on ionode 207 of 216 on ionode 208 of 216 on ionode 209 of 216 on ionode 210 of 216 on ionode 211 of 216 on ionode 212 of 216 on ionode 213 of 216 on ionode 214 of 216 on ionode 215 of 216 on ionode 216 of 216 on ionode MMN calculated Running Wannier90 Wannier Function centers (cartesian, alat) and spreads (ang): ( 0.07840 0.07840 0.07840) : 2.27123 ( 0.07840 -0.07840 -0.07840) : 2.27123 ( -0.07840 0.07840 -0.07840) : 2.27123 ( -0.07840 -0.07840 0.07840) : 2.27123 ------------------------------------------------------------------- WANNIER : 112.01s CPU ------------------------------------------------------------------- Dipole matrix elements calculated Calculating kmap and kgmap Progress kmap: ######################################## Progress kgmap: ######################################## kmaps : 9.39s CPU =================================================================== irreducible q point # 1 =================================================================== Imposing acoustic sum rule on the dynamical matrix q( 1 ) = ( 0.0000000 0.0000000 0.0000000 ) =================================================================== irreducible q point # 2 =================================================================== q( 2 ) = ( -0.1666667 0.1666667 -0.1666667 ) q( 3 ) = ( 0.1666667 -0.1666667 -0.1666667 ) q( 4 ) = ( 0.1666667 0.1666667 0.1666667 ) q( 5 ) = ( -0.1666667 -0.1666667 0.1666667 ) q( 6 ) = ( 0.1666667 -0.1666667 0.1666667 ) q( 7 ) = ( -0.1666667 0.1666667 0.1666667 ) q( 8 ) = ( 0.1666667 0.1666667 -0.1666667 ) q( 9 ) = ( -0.1666667 -0.1666667 -0.1666667 ) =================================================================== irreducible q point # 3 =================================================================== q( 10 ) = ( -0.3333333 0.3333333 -0.3333333 ) q( 11 ) = ( 0.3333333 -0.3333333 -0.3333333 ) q( 12 ) = ( 0.3333333 0.3333333 0.3333333 ) q( 13 ) = ( -0.3333333 -0.3333333 0.3333333 ) q( 14 ) = ( 0.3333333 -0.3333333 0.3333333 ) q( 15 ) = ( -0.3333333 0.3333333 0.3333333 ) q( 16 ) = ( 0.3333333 0.3333333 -0.3333333 ) q( 17 ) = ( -0.3333333 -0.3333333 -0.3333333 ) =================================================================== irreducible q point # 4 =================================================================== q( 18 ) = ( 0.5000000 -0.5000000 0.5000000 ) q( 19 ) = ( -0.5000000 0.5000000 0.5000000 ) q( 20 ) = ( -0.5000000 -0.5000000 -0.5000000 ) q( 21 ) = ( 0.5000000 0.5000000 -0.5000000 ) =================================================================== irreducible q point # 5 =================================================================== q( 22 ) = ( 0.0000000 0.3333333 0.0000000 ) q( 23 ) = ( 0.0000000 -0.3333333 0.0000000 ) q( 24 ) = ( 0.3333333 0.0000000 0.0000000 ) q( 25 ) = ( -0.3333333 0.0000000 0.0000000 ) q( 26 ) = ( 0.0000000 0.0000000 0.3333333 ) q( 27 ) = ( 0.0000000 0.0000000 -0.3333333 ) =================================================================== irreducible q point # 6 =================================================================== q( 28 ) = ( -0.1666667 0.5000000 -0.1666667 ) q( 29 ) = ( 0.1666667 -0.5000000 -0.1666667 ) q( 30 ) = ( 0.1666667 0.5000000 0.1666667 ) q( 31 ) = ( -0.1666667 -0.5000000 0.1666667 ) q( 32 ) = ( 0.5000000 -0.1666667 0.1666667 ) q( 33 ) = ( -0.5000000 0.1666667 0.1666667 ) q( 34 ) = ( 0.5000000 0.1666667 -0.1666667 ) q( 35 ) = ( -0.5000000 -0.1666667 -0.1666667 ) q( 36 ) = ( -0.1666667 -0.5000000 -0.1666667 ) q( 37 ) = ( 0.1666667 -0.5000000 0.1666667 ) q( 38 ) = ( -0.1666667 0.5000000 0.1666667 ) q( 39 ) = ( 0.1666667 0.5000000 -0.1666667 ) q( 40 ) = ( 0.1666667 -0.1666667 0.5000000 ) q( 41 ) = ( 0.1666667 0.1666667 -0.5000000 ) q( 42 ) = ( -0.1666667 -0.1666667 -0.5000000 ) q( 43 ) = ( -0.1666667 0.1666667 0.5000000 ) q( 44 ) = ( 0.5000000 -0.1666667 -0.1666667 ) q( 45 ) = ( -0.5000000 -0.1666667 0.1666667 ) q( 46 ) = ( 0.5000000 0.1666667 0.1666667 ) q( 47 ) = ( -0.5000000 0.1666667 -0.1666667 ) q( 48 ) = ( -0.1666667 -0.1666667 0.5000000 ) q( 49 ) = ( -0.1666667 0.1666667 -0.5000000 ) q( 50 ) = ( 0.1666667 0.1666667 0.5000000 ) q( 51 ) = ( 0.1666667 -0.1666667 -0.5000000 ) =================================================================== irreducible q point # 7 =================================================================== q( 52 ) = ( 0.6666667 -0.3333333 0.6666667 ) q( 53 ) = ( -0.6666667 0.3333333 0.6666667 ) q( 54 ) = ( -0.6666667 -0.3333333 -0.6666667 ) q( 55 ) = ( 0.6666667 0.3333333 -0.6666667 ) q( 56 ) = ( -0.3333333 0.6666667 -0.6666667 ) q( 57 ) = ( 0.3333333 -0.6666667 -0.6666667 ) q( 58 ) = ( -0.3333333 -0.6666667 0.6666667 ) q( 59 ) = ( 0.3333333 0.6666667 0.6666667 ) q( 60 ) = ( 0.6666667 0.3333333 0.6666667 ) q( 61 ) = ( -0.6666667 0.3333333 -0.6666667 ) q( 62 ) = ( 0.6666667 -0.3333333 -0.6666667 ) q( 63 ) = ( -0.6666667 -0.3333333 0.6666667 ) q( 64 ) = ( -0.6666667 0.6666667 -0.3333333 ) q( 65 ) = ( -0.6666667 -0.6666667 0.3333333 ) q( 66 ) = ( 0.6666667 0.6666667 0.3333333 ) q( 67 ) = ( 0.6666667 -0.6666667 -0.3333333 ) q( 68 ) = ( -0.3333333 0.6666667 0.6666667 ) q( 69 ) = ( 0.3333333 0.6666667 -0.6666667 ) q( 70 ) = ( -0.3333333 -0.6666667 -0.6666667 ) q( 71 ) = ( 0.3333333 -0.6666667 0.6666667 ) q( 72 ) = ( 0.6666667 0.6666667 -0.3333333 ) q( 73 ) = ( 0.6666667 -0.6666667 0.3333333 ) q( 74 ) = ( -0.6666667 -0.6666667 -0.3333333 ) q( 75 ) = ( -0.6666667 0.6666667 0.3333333 ) =================================================================== irreducible q point # 8 =================================================================== q( 76 ) = ( 0.5000000 -0.1666667 0.5000000 ) q( 77 ) = ( -0.5000000 0.1666667 0.5000000 ) q( 78 ) = ( -0.5000000 -0.1666667 -0.5000000 ) q( 79 ) = ( 0.5000000 0.1666667 -0.5000000 ) q( 80 ) = ( -0.1666667 0.5000000 -0.5000000 ) q( 81 ) = ( 0.1666667 -0.5000000 -0.5000000 ) q( 82 ) = ( -0.1666667 -0.5000000 0.5000000 ) q( 83 ) = ( 0.1666667 0.5000000 0.5000000 ) q( 84 ) = ( 0.5000000 0.1666667 0.5000000 ) q( 85 ) = ( -0.5000000 0.1666667 -0.5000000 ) q( 86 ) = ( 0.5000000 -0.1666667 -0.5000000 ) q( 87 ) = ( -0.5000000 -0.1666667 0.5000000 ) q( 88 ) = ( -0.5000000 0.5000000 -0.1666667 ) q( 89 ) = ( -0.5000000 -0.5000000 0.1666667 ) q( 90 ) = ( 0.5000000 0.5000000 0.1666667 ) q( 91 ) = ( 0.5000000 -0.5000000 -0.1666667 ) q( 92 ) = ( -0.1666667 0.5000000 0.5000000 ) q( 93 ) = ( 0.1666667 0.5000000 -0.5000000 ) q( 94 ) = ( -0.1666667 -0.5000000 -0.5000000 ) q( 95 ) = ( 0.1666667 -0.5000000 0.5000000 ) q( 96 ) = ( 0.5000000 0.5000000 -0.1666667 ) q( 97 ) = ( 0.5000000 -0.5000000 0.1666667 ) q( 98 ) = ( -0.5000000 -0.5000000 -0.1666667 ) q( 99 ) = ( -0.5000000 0.5000000 0.1666667 ) =================================================================== irreducible q point # 9 =================================================================== q( 100 ) = ( 0.3333333 0.0000000 0.3333333 ) q( 101 ) = ( -0.3333333 0.0000000 0.3333333 ) q( 102 ) = ( -0.3333333 0.0000000 -0.3333333 ) q( 103 ) = ( 0.3333333 0.0000000 -0.3333333 ) q( 104 ) = ( 0.0000000 0.3333333 -0.3333333 ) q( 105 ) = ( 0.0000000 -0.3333333 -0.3333333 ) q( 106 ) = ( 0.0000000 -0.3333333 0.3333333 ) q( 107 ) = ( 0.0000000 0.3333333 0.3333333 ) q( 108 ) = ( -0.3333333 0.3333333 0.0000000 ) q( 109 ) = ( -0.3333333 -0.3333333 0.0000000 ) q( 110 ) = ( 0.3333333 0.3333333 0.0000000 ) q( 111 ) = ( 0.3333333 -0.3333333 0.0000000 ) =================================================================== irreducible q point # 10 =================================================================== q( 112 ) = ( 0.0000000 0.6666667 0.0000000 ) q( 113 ) = ( 0.0000000 -0.6666667 0.0000000 ) q( 114 ) = ( 0.6666667 0.0000000 0.0000000 ) q( 115 ) = ( -0.6666667 0.0000000 0.0000000 ) q( 116 ) = ( 0.0000000 0.0000000 0.6666667 ) q( 117 ) = ( 0.0000000 0.0000000 -0.6666667 ) =================================================================== irreducible q point # 11 =================================================================== q( 118 ) = ( 0.8333333 -0.1666667 0.8333333 ) q( 119 ) = ( -0.8333333 0.1666667 0.8333333 ) q( 120 ) = ( -0.8333333 -0.1666667 -0.8333333 ) q( 121 ) = ( 0.8333333 0.1666667 -0.8333333 ) q( 122 ) = ( -0.1666667 0.8333333 -0.8333333 ) q( 123 ) = ( 0.1666667 -0.8333333 -0.8333333 ) q( 124 ) = ( -0.1666667 -0.8333333 0.8333333 ) q( 125 ) = ( 0.1666667 0.8333333 0.8333333 ) q( 126 ) = ( 0.8333333 0.1666667 0.8333333 ) q( 127 ) = ( -0.8333333 0.1666667 -0.8333333 ) q( 128 ) = ( 0.8333333 -0.1666667 -0.8333333 ) q( 129 ) = ( -0.8333333 -0.1666667 0.8333333 ) q( 130 ) = ( -0.8333333 0.8333333 -0.1666667 ) q( 131 ) = ( -0.8333333 -0.8333333 0.1666667 ) q( 132 ) = ( 0.8333333 0.8333333 0.1666667 ) q( 133 ) = ( 0.8333333 -0.8333333 -0.1666667 ) q( 134 ) = ( -0.1666667 0.8333333 0.8333333 ) q( 135 ) = ( 0.1666667 0.8333333 -0.8333333 ) q( 136 ) = ( -0.1666667 -0.8333333 -0.8333333 ) q( 137 ) = ( 0.1666667 -0.8333333 0.8333333 ) q( 138 ) = ( 0.8333333 0.8333333 -0.1666667 ) q( 139 ) = ( 0.8333333 -0.8333333 0.1666667 ) q( 140 ) = ( -0.8333333 -0.8333333 -0.1666667 ) q( 141 ) = ( -0.8333333 0.8333333 0.1666667 ) =================================================================== irreducible q point # 12 =================================================================== q( 142 ) = ( 0.6666667 0.0000000 0.6666667 ) q( 143 ) = ( -0.6666667 0.0000000 0.6666667 ) q( 144 ) = ( -0.6666667 0.0000000 -0.6666667 ) q( 145 ) = ( 0.6666667 0.0000000 -0.6666667 ) q( 146 ) = ( 0.0000000 0.6666667 -0.6666667 ) q( 147 ) = ( 0.0000000 -0.6666667 -0.6666667 ) q( 148 ) = ( 0.0000000 -0.6666667 0.6666667 ) q( 149 ) = ( 0.0000000 0.6666667 0.6666667 ) q( 150 ) = ( -0.6666667 0.6666667 0.0000000 ) q( 151 ) = ( -0.6666667 -0.6666667 0.0000000 ) q( 152 ) = ( 0.6666667 0.6666667 0.0000000 ) q( 153 ) = ( 0.6666667 -0.6666667 0.0000000 ) =================================================================== irreducible q point # 13 =================================================================== q( 154 ) = ( 0.0000000 -1.0000000 0.0000000 ) q( 155 ) = ( -1.0000000 0.0000000 0.0000000 ) q( 156 ) = ( 0.0000000 0.0000000 -1.0000000 ) =================================================================== irreducible q point # 14 =================================================================== q( 157 ) = ( 0.6666667 -0.3333333 1.0000000 ) q( 158 ) = ( -0.6666667 0.3333333 1.0000000 ) q( 159 ) = ( -0.6666667 -0.3333333 -1.0000000 ) q( 160 ) = ( 0.6666667 0.3333333 -1.0000000 ) q( 161 ) = ( -0.3333333 0.6666667 -1.0000000 ) q( 162 ) = ( 0.3333333 -0.6666667 -1.0000000 ) q( 163 ) = ( -0.3333333 -0.6666667 1.0000000 ) q( 164 ) = ( 0.3333333 0.6666667 1.0000000 ) q( 165 ) = ( 1.0000000 0.3333333 0.6666667 ) q( 166 ) = ( -1.0000000 0.3333333 -0.6666667 ) q( 167 ) = ( 1.0000000 -0.3333333 -0.6666667 ) q( 168 ) = ( -1.0000000 -0.3333333 0.6666667 ) q( 169 ) = ( -0.6666667 1.0000000 -0.3333333 ) q( 170 ) = ( -0.6666667 -1.0000000 0.3333333 ) q( 171 ) = ( 0.6666667 1.0000000 0.3333333 ) q( 172 ) = ( 0.6666667 -1.0000000 -0.3333333 ) q( 173 ) = ( -0.3333333 1.0000000 0.6666667 ) q( 174 ) = ( 0.3333333 1.0000000 -0.6666667 ) q( 175 ) = ( -0.3333333 -1.0000000 -0.6666667 ) q( 176 ) = ( 0.3333333 -1.0000000 0.6666667 ) q( 177 ) = ( 1.0000000 0.6666667 -0.3333333 ) q( 178 ) = ( 1.0000000 -0.6666667 0.3333333 ) q( 179 ) = ( -1.0000000 -0.6666667 -0.3333333 ) q( 180 ) = ( -1.0000000 0.6666667 0.3333333 ) =================================================================== irreducible q point # 15 =================================================================== q( 181 ) = ( 0.5000000 -0.1666667 0.8333333 ) q( 182 ) = ( -0.5000000 0.1666667 0.8333333 ) q( 183 ) = ( -0.5000000 -0.1666667 -0.8333333 ) q( 184 ) = ( 0.5000000 0.1666667 -0.8333333 ) q( 185 ) = ( -0.1666667 0.5000000 -0.8333333 ) q( 186 ) = ( 0.1666667 -0.5000000 -0.8333333 ) q( 187 ) = ( -0.1666667 -0.5000000 0.8333333 ) q( 188 ) = ( 0.1666667 0.5000000 0.8333333 ) q( 189 ) = ( 0.8333333 0.1666667 0.5000000 ) q( 190 ) = ( -0.8333333 0.1666667 -0.5000000 ) q( 191 ) = ( 0.8333333 -0.1666667 -0.5000000 ) q( 192 ) = ( -0.8333333 -0.1666667 0.5000000 ) q( 193 ) = ( -0.5000000 0.1666667 -0.8333333 ) q( 194 ) = ( 0.5000000 -0.1666667 -0.8333333 ) q( 195 ) = ( 0.5000000 0.1666667 0.8333333 ) q( 196 ) = ( -0.5000000 -0.1666667 0.8333333 ) q( 197 ) = ( 0.1666667 -0.5000000 0.8333333 ) q( 198 ) = ( -0.1666667 0.5000000 0.8333333 ) q( 199 ) = ( 0.1666667 0.5000000 -0.8333333 ) q( 200 ) = ( -0.1666667 -0.5000000 -0.8333333 ) q( 201 ) = ( -0.8333333 -0.1666667 -0.5000000 ) q( 202 ) = ( 0.8333333 -0.1666667 0.5000000 ) q( 203 ) = ( -0.8333333 0.1666667 0.5000000 ) q( 204 ) = ( 0.8333333 0.1666667 -0.5000000 ) =================================================================== irreducible q point # 16 =================================================================== q( 205 ) = ( -0.3333333 -1.0000000 -0.0000000 ) q( 206 ) = ( 0.3333333 1.0000000 0.0000000 ) q( 207 ) = ( -1.0000000 -0.3333333 0.0000000 ) q( 208 ) = ( 1.0000000 0.3333333 0.0000000 ) q( 209 ) = ( -0.0000000 1.0000000 -0.3333333 ) q( 210 ) = ( 0.0000000 1.0000000 0.3333333 ) q( 211 ) = ( 0.3333333 0.0000000 -1.0000000 ) q( 212 ) = ( -0.3333333 -0.0000000 1.0000000 ) q( 213 ) = ( -1.0000000 0.0000000 -0.3333333 ) q( 214 ) = ( 1.0000000 -0.0000000 0.3333333 ) q( 215 ) = ( -0.0000000 -0.3333333 -1.0000000 ) q( 216 ) = ( 0.0000000 0.3333333 1.0000000 ) Writing epmatq on .epb files band disentanglement is used: nbndsub = 4 Using uniform MP q-mesh: 4 4 4 Size of q point mesh for interpolation: 8 Using uniform MP k-mesh: 30 30 30 Size of k point mesh for interpolation: 1504 Max number of k points per pool: 1504 Number of ep-matrix elements per pool : 36096 ~= 282.00 Kb (@ 8 bytes/ DP) Loading Wannier rep into memory: 57.01 MB =================================================================== Phonon (Imaginary) Self-Energy in the Migdal Approximation (on the fly) =================================================================== Fermi Surface thickness = 0.500000 Ry Golden Rule strictly enforced with T = 0.00000029 Ry Gaussian Broadening: 0.045 Ry, ngauss= 1 DOS = 3.398804 states/spin/Ry/Unit Cell at Ef= 5.702811 eV DOS = 3.398804 states/spin/Ry/Unit Cell at Ef= 5.812463 eV iq = 1 coord.: 0.00000 0.00000 0.00000 wt: 0.03125 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.000 gamma= 0.000 meV omega= -0.003 meV lambda( 2 )= 0.000 gamma= 0.000 meV omega= -0.003 meV lambda( 3 )= 0.000 gamma= 0.000 meV omega= -0.003 meV lambda( sum )= 0.000 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 745 out of 752 iq = 2 coord.: 0.00000 0.00000 0.25000 wt: 0.25000 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.466 gamma= 0.004 meV omega= 3.377 meV lambda( 2 )= 0.123 gamma= 0.001 meV omega= 3.377 meV lambda( 3 )= 0.687 gamma= 0.021 meV omega= 6.287 meV lambda( sum )= 1.276 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 738 out of 752 iq = 3 coord.: 0.00000 0.00000 -0.50000 wt: 0.12500 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.199 gamma= 0.003 meV omega= 4.050 meV lambda( 2 )= 0.110 gamma= 0.001 meV omega= 4.050 meV lambda( 3 )= 0.391 gamma= 0.021 meV omega= 8.342 meV lambda( sum )= 0.701 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 733 out of 752 iq = 4 coord.: 0.00000 0.25000 0.25000 wt: 0.18750 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.011 gamma= 0.000 meV omega= 5.094 meV lambda( 2 )= 0.169 gamma= 0.003 meV omega= 5.094 meV lambda( 3 )= 0.875 gamma= 0.035 meV omega= 7.173 meV lambda( sum )= 1.056 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 745 out of 752 iq = 5 coord.: 0.00000 0.25000 -0.50000 wt: 0.75000 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.004 gamma= 0.000 meV omega= 3.780 meV lambda( 2 )= 0.096 gamma= 0.003 meV omega= 5.891 meV lambda( 3 )= 0.046 gamma= 0.002 meV omega= 7.636 meV lambda( sum )= 0.145 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 745 out of 752 iq = 6 coord.: 0.00000 0.25000 -0.25000 wt: 0.37500 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.114 gamma= 0.001 meV omega= 4.011 meV lambda( 2 )= 0.604 gamma= 0.017 meV omega= 6.000 meV lambda( 3 )= 3.347 gamma= 0.145 meV omega= 7.442 meV lambda( sum )= 4.064 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 745 out of 752 iq = 7 coord.: 0.00000 -0.50000 -0.50000 wt: 0.09375 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.008 gamma= 0.000 meV omega= 6.040 meV lambda( 2 )= 0.236 gamma= 0.007 meV omega= 6.040 meV lambda( 3 )= 0.864 gamma= 0.037 meV omega= 7.436 meV lambda( sum )= 1.108 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 734 out of 752 iq = 8 coord.: 0.25000 -0.50000 -0.25000 wt: 0.18750 ismear: 1 ------------------------------------------------------------------- lambda( 1 )= 0.375 gamma= 0.010 meV omega= 5.922 meV lambda( 2 )= 1.473 gamma= 0.050 meV omega= 6.602 meV lambda( 3 )= 0.975 gamma= 0.033 meV omega= 6.602 meV lambda( sum )= 2.823 ------------------------------------------------------------------- Number of (k,k+q) pairs on the Fermi surface: 745 out of 752 =================================================================== Eliashberg Spectral Function in the Migdal Approximation =================================================================== a2F : 0.02s CPU ep-interp : 9.72s CPU ( 8 calls, 1.215 s avg) ephwann : 44.26s CPU INITIALIZATION: init_vloc : 0.04s CPU ( 18 calls, 0.002 s avg) init_us_1 : 0.77s CPU ( 18 calls, 0.043 s avg) Electron-Phonon interpolation ephwann : 44.26s CPU ep-interp : 9.72s CPU ( 8 calls, 1.215 s avg) PH SELF-ENER : 0.00s CPU ( 8 calls, 0.000 s avg) Ham: step 1 : -0.00s CPU Ham: step 2 : 0.02s CPU ep: step 1 : 0.40s CPU ( 648 calls, 0.001 s avg) ep: step 2 : 9.50s CPU ( 648 calls, 0.015 s avg) DynW2B : -0.00s CPU ( 16 calls, -0.000 s avg) HamW2B : 5.98s CPU ( 61664 calls, 0.000 s avg) ephW2Bp : 0.99s CPU ( 8 calls, 0.124 s avg) Total program execution EPW : 1h 7m CPU time, 1h 8m wall time