
=== q1: q = (0.01500, 0.00866, 0.00000) 2pi/a, |q| = 0.0342 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.078   Z       1903.537   276.884     9.461         7.394
     2    0.237   T          3.481     0.882     0.030         0.162
     3    0.336   L         36.774    11.098     0.379         0.137
     4    8.291   T          2.027     3.041     0.104        39.223
     5    8.297   L         10.079    15.122     0.517        39.220
     6   11.197   T          1.515     2.641     0.090         0.013
     7   11.247   L        198.487   346.728    11.847         0.013
     8   11.893   Z         55.038    98.866     3.378         0.001
     9   13.738   Z        130.406   251.763     8.602         1.323

=== q2: q = (0.03000, 0.01732, 0.00000) 2pi/a, |q| = 0.0683 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.053   Z        107.842     6.463     0.442        18.114
     2    0.450   T          2.212     0.386     0.026         0.021
     3    0.714   L         37.881     8.336     0.570         0.017
     4    8.297   T          0.103     0.077     0.005        39.321
     5    8.318   L          0.930     0.698     0.048        39.313
     6   11.195   T          4.196     3.656     0.250         0.002
     7   11.234   L        110.215    96.206     6.574         0.002
     8   11.895   Z         52.245    46.928     3.207         0.000
     9   13.720   Z          6.266     6.045     0.413         2.693

=== q3: q = (0.06000, 0.03464, 0.00000) 2pi/a, |q| = 0.1367 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.176   Z         35.250     1.928     0.263        20.428
     2    0.881   T          5.779     0.706     0.097         0.030
     3    1.430   L         42.080     6.554     0.896         0.027
     4    8.318   T          0.078     0.029     0.004        39.683
     5    8.399   L          1.063     0.401     0.055        39.666
     6   11.182   L         73.562    32.031     4.378         0.004
     7   11.185   T          7.964     3.468     0.474         0.005
     8   11.905   Z         42.788    19.224     2.627         0.001
     9   13.650   Z          4.054     1.950     0.267         5.758

=== q4: q = (0.10000, 0.05774, 0.00000) 2pi/a, |q| = 0.2278 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.507   Z         15.417     0.858     0.195         2.553
     2    1.449   T         10.006     0.941     0.214        25.069
     3    2.350   L         48.384     5.795     1.320        22.324
     4    8.366   T          0.078     0.018     0.004         0.012
     5    8.572   L          1.220     0.279     0.064         0.198
     6   11.067   L         61.461    15.974     3.639        29.887
     7   11.161   T         12.387     3.233     0.737        21.636
     8   11.923   Z         27.198     7.338     1.671         5.512
     9   13.489   Z          3.350     0.961     0.219         0.561

=== q5: q = (0.15000, 0.08660, 0.00000) 2pi/a, |q| = 0.3417 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    1.116   Z          7.820     0.430     0.147         2.935
     2    2.126   T         12.587     0.956     0.327        35.720
     3    3.428   L         51.898     5.005     1.710        42.216
     4    8.448   T          0.077     0.012     0.004         0.024
     5    8.854   L          1.224     0.190     0.065         0.448
     6   10.875   L         49.077     8.430     2.880        41.875
     7   11.112   T         17.838     3.097     1.058        28.549
     8   11.948   Z         12.740     2.294     0.784        10.665
     9   13.194   Z          2.831     0.536     0.183         1.071

=== q6: q = (0.03333, 0.05774, 0.00000) 2pi/a, |q| = 0.1315 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.207   Z          0.000     0.000     0.000        24.000
     2    0.855   T         26.486     3.315     0.436         0.000
     3    1.368   L         31.792     5.032     0.662         0.000
     4    8.316   T          0.000     0.000     0.000        37.199
     5    8.391   L          0.000     0.000     0.000        41.812
     6   11.179   T         43.905    19.865     2.613         0.000
     7   11.194   T         60.882    27.566     3.625         0.000
     8   11.904   Z         44.307    20.687     2.721         0.000
     9   13.657   Z          0.000     0.000     0.000         3.601

=== q7: q = (0.06667, 0.11547, 0.00000) 2pi/a, |q| = 0.2630 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.698   Z          0.000     0.000     0.000         0.000
     2    1.710   T         27.086     2.397     0.630        27.131
     3    2.659   L         38.210     4.216     1.109         0.281
     4    8.392   T          0.000     0.000     0.000         0.000
     5    8.653   L          0.000     0.000     0.000         0.000
     6   11.002   L         52.331    11.745     3.089        20.925
     7   11.159   T         32.484     7.342     1.931        28.109
     8   11.931   Z         23.327     5.452     1.434         8.236
     9   13.408   Z          0.000     0.000     0.000         0.000

=== q8: q = (0.33333, 0.57735, 0.00000) 2pi/a, |q| = 1.3151 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    5.334   Z          2.271     0.071     0.093         0.004
     2    5.438   T          0.041     0.001     0.002         0.121
     3    6.904   T          0.042     0.001     0.002         0.003
     4    9.302   Z          0.019     0.001     0.001        57.551
     5    9.810   T         29.695     1.259     1.655         0.001
     6    9.827   T          0.004     0.000     0.000         0.004
     7   10.828   T          1.831     0.082     0.107         0.002
     8   11.045   T          0.034     0.002     0.002         0.003
     9   11.469   Z          0.023     0.001     0.001         0.018

=== q9: q = (-0.16100, -0.27886, 0.00000) 2pi/a, |q| = 0.6352 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    2.950   Z          0.000     0.000     0.000         0.000
     2    3.985   T          6.205     0.347     0.220        14.929
     3    5.376   L         27.816     1.807     1.148       119.718
     4    8.779   T          0.000     0.000     0.000         0.000
     5    9.540   L          0.000     0.000     0.000         0.000
     6   10.315   L         25.119     2.260     1.436        60.787
     7   11.001   T          9.597     0.892     0.567         2.048
     8   11.916   Z         16.877     1.632     1.037        61.214
     9   12.128   Z          0.000     0.000     0.000         0.000

=== q10: q = (0.06699, 0.01732, 0.00000) 2pi/a, |q| = 0.1365 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.188   Z          0.000     0.000     0.000        13.134
     2    0.882   T         18.959     2.322     0.317         0.000
     3    1.424   L         42.156     6.558     0.895         0.000
     4    8.319   T          0.000     0.000     0.000        41.389
     5    8.399   L          0.000     0.000     0.000        37.993
     6   11.178   L         62.906    27.423     3.743         0.000
     7   11.190   T         39.205    17.100     2.334         0.000
     8   11.904   Z         42.226    18.997     2.593         0.000
     9   13.651   Z          0.000     0.000     0.000         6.838

=== q11: q = (0.03464, 0.06000, 0.00000) 2pi/a, |q| = 0.1367 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.217   Z          0.000     0.000     0.000        24.573
     2    0.888   T         26.243     3.220     0.440         0.000
     3    1.422   L         32.184     4.997     0.683         0.000
     4    8.318   T          0.000     0.000     0.000        37.140
     5    8.399   L          0.000     0.000     0.000        41.914
     6   11.175   T         49.920    21.731     2.970         0.000
     7   11.192   T         54.503    23.743     3.245         0.000
     8   11.905   Z         43.514    19.551     2.672         0.000
     9   13.650   Z          0.000     0.000     0.000         3.686

Xi_LA(q) [eV] from the small-q points: |q|=0.068/A: 8.34, |q|=0.132/A: 5.03, |q|=0.136/A: 6.56, |q|=0.137/A: 5.00, |q|=0.137/A: 6.55, |q|=0.228/A: 5.79, |q|=0.263/A: 4.22, |q|=0.342/A: 5.00
  -> Xi_LA: q -> 0 extrapolation 7.55 eV, mean over the window 5.81 eV   [Kaasbjerg: LA 2.8, TA 1.6]

Xi_TA(q) [eV] from the small-q points: |q|=0.068/A: 0.39, |q|=0.132/A: 3.31, |q|=0.136/A: 2.32, |q|=0.137/A: 3.22, |q|=0.137/A: 0.71, |q|=0.228/A: 0.94, |q|=0.263/A: 2.40, |q|=0.342/A: 0.96
  -> Xi_TA: q -> 0 extrapolation 2.04 eV, mean over the window 1.78 eV   [Kaasbjerg: LA 2.8, TA 1.6]
