
=== q1: q = (0.01500, 0.02598, 0.00000) 2pi/a, |q| = 0.0592 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.108   Z          0.000     0.000     0.000         0.000
     2    0.383   T         34.520     6.427     0.380         4.886
     3    0.617   L         11.119     2.626     0.155         6.304
     4    8.295   T          0.000     0.000     0.000         0.000
     5    8.311   L          0.000     0.000     0.000         0.000
     6   11.196   T         31.074    31.267     1.850         8.418
     7   11.238   L         99.805   100.617     5.955        12.836
     8   11.895   Z         90.039    93.387     5.527        24.901
     9   13.727   Z          0.000     0.000     0.000         0.000

=== q2: q = (0.03000, 0.05196, 0.00000) 2pi/a, |q| = 0.1184 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.172   Z          0.000     0.000     0.000         0.000
     2    0.769   T         25.551     3.369     0.399         3.119
     3    1.232   L         16.201     2.704     0.320         7.741
     4    8.311   T          0.000     0.000     0.000         0.000
     5    8.372   L          0.000     0.000     0.000         0.000
     6   11.186   T         13.970     7.026     0.832        17.839
     7   11.202   L         41.022    20.644     2.443        11.671
     8   11.902   Z         87.349    45.311     5.363        13.236
     9   13.674   Z          0.000     0.000     0.000         0.000

=== q3: q = (0.06000, 0.10392, 0.00000) 2pi/a, |q| = 0.2367 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.580   Z          0.000     0.000     0.000         0.000
     2    1.540   T         19.105     1.783     0.422         9.997
     3    2.411   L         25.296     2.953     0.699         1.434
     4    8.374   T          0.000     0.000     0.000         0.000
     5    8.592   L          0.000     0.000     0.000         0.000
     6   11.046   L         24.392     6.095     1.443        13.125
     7   11.166   T         20.512     5.153     1.220         6.786
     8   11.925   Z         69.158    17.954     4.250        30.403
     9   13.468   Z          0.000     0.000     0.000         0.000

=== q4: q = (0.10000, 0.17321, 0.00000) 2pi/a, |q| = 0.3945 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    1.418   Z          0.000     0.000     0.000         0.000
     2    2.547   T          7.325     0.527     0.208        24.928
     3    3.813   L         30.038     2.646     1.044         9.223
     4    8.503   T          0.000     0.000     0.000         0.000
     5    8.987   L          0.000     0.000     0.000         0.000
     6   10.753   L         36.986     5.471     2.159        15.460
     7   11.111   T         14.672     2.206     0.870         9.393
     8   11.951   Z         23.236     3.623     1.430        56.678
     9   13.028   Z          0.000     0.000     0.000         0.000

=== q5: q = (0.15000, 0.25981, 0.00000) 2pi/a, |q| = 0.5918 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    2.675   Z          0.000     0.000     0.000         5.399
     2    3.738   T          4.734     0.275     0.163         0.000
     3    5.156   L         29.977     2.047     1.211         0.000
     4    8.723   T          0.000     0.000     0.000         7.283
     5    9.457   L          0.000     0.000     0.000        98.619
     6   10.381   L         28.363     2.748     1.626         0.000
     7   11.022   T          9.976     0.996     0.589         0.000
     8   11.932   Z         11.203     1.164     0.689         0.000
     9   12.299   Z          0.000     0.000     0.000        66.740

=== q6: q = (-0.05196, 0.03000, 0.00000) 2pi/a, |q| = 0.1184 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.147   Z         45.088     2.600     0.308         0.551
     2    0.766   T          2.201     0.290     0.034        14.374
     3    1.240   L         19.201     3.215     0.381        12.958
     4    8.311   T          0.048     0.021     0.002         0.001
     5    8.372   L          1.102     0.479     0.057         0.014
     6   11.188   T          0.960     0.483     0.057         7.529
     7   11.199   L         55.567    27.960     3.309         6.337
     8   11.902   Z         89.675    46.517     5.506        30.776
     9   13.673   Z          4.546     2.528     0.299         0.058

=== q7: q = (-0.10392, 0.06000, 0.00000) 2pi/a, |q| = 0.2367 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    0.568   Z         14.899     0.844     0.200         1.306
     2    1.503   T          6.392     0.589     0.139        19.940
     3    2.438   L         10.575     1.241     0.294        20.955
     4    8.371   T          0.039     0.008     0.002         0.006
     5    8.593   L          1.322     0.291     0.069         0.121
     6   11.053   L         33.911     8.476     2.006        10.907
     7   11.158   T          3.981     1.000     0.237         9.654
     8   11.925   Z         82.563    21.435     5.074        36.490
     9   13.469   Z          3.362     0.928     0.220         0.269

=== q8: q = (-0.17321, 0.10000, 0.00000) 2pi/a, |q| = 0.3945 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    1.431   Z          1.071     0.058     0.023         1.161
     2    2.424   T         22.497     1.580     0.623        14.927
     3    3.892   L         28.142     2.504     0.988        20.034
     4    8.491   T          0.005     0.001     0.000         0.006
     5    8.993   L          0.211     0.028     0.011         0.236
     6   10.780   L         19.962     2.956     1.166        26.007
     7   11.082   T         13.616     2.045     0.807         8.176
     8   11.958   Z         27.893     4.351     1.717        79.000
     9   13.028   Z          0.459     0.075     0.030         0.479

=== q9: q = (-0.34333, 0.00000, 0.00000) 2pi/a, |q| = 0.6773 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    3.210   Z          0.000     0.000     0.000         0.000
     2    4.215   T         21.378     1.153     0.781        16.165
     3    5.562   L         92.228     5.715     3.871        23.042
     4    8.837   T          0.000     0.000     0.000         0.000
     5    9.610   L          0.000     0.000     0.000         0.000
     6   10.261   L         13.018     1.096     0.742        42.534
     7   10.982   T         23.975     2.088     1.414        27.875
     8   11.896   Z         33.713     3.055     2.069        51.643
     9   11.965   Z          0.000     0.000     0.000         0.000

=== q10: q = (-0.51500, -0.29734, 0.00000) 2pi/a, |q| = 1.1731 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    4.578   T          0.002     0.000     0.000        70.170
     2    5.042   Z          0.174     0.006     0.007         0.137
     3    6.885   L         41.393     1.648     1.933        91.219
     4    8.906   T          0.000     0.000     0.000         0.002
     5    9.672   Z          0.013     0.001     0.001         0.003
     6   10.347   L         55.014     2.685     3.149        31.053
     7   10.612   T          0.002     0.000     0.000        28.515
     8   11.245   L          0.133     0.007     0.008         0.107
     9   11.871   Z         48.375     2.529     2.966        20.584

=== q11: q = (-0.34333, -0.59467, 0.00000) 2pi/a, |q| = 1.3546 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    5.328   Z          0.000     0.000     0.000         0.000
     2    5.441   T         70.266     2.153     2.917        12.358
     3    6.889   T        109.043     3.760     5.094        51.372
     4    9.307   Z          0.000     0.000     0.000         0.000
     5    9.817   T         25.790     1.062     1.438        34.269
     6    9.820   T          0.000     0.000     0.000         0.000
     7   10.832   T          0.000     0.000     0.000         0.000
     8   11.039   T          7.050     0.308     0.417         9.303
     9   11.473   Z         40.328     1.795     2.431        37.580

=== q12: q = (-0.50500, -0.87469, 0.00000) 2pi/a, |q| = 1.9924 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    4.585   L         40.891     0.782     1.558         0.000
     2    5.055   Z          0.000     0.000     0.000        36.790
     3    6.888   T         99.409     2.330     4.643         0.000
     4    8.909   L          0.000     0.000     0.000        12.662
     5    9.667   Z          0.000     0.000     0.000        13.527
     6   10.350   T          4.526     0.130     0.259         0.000
     7   10.612   L         49.335     1.436     2.860         0.000
     8   11.242   T          0.000     0.000     0.000        68.986
     9   11.868   Z         80.730     2.484     4.950         0.000

=== q13: q = (0.16167, 0.28001, 0.00000) 2pi/a, |q| = 0.6378 1/A
  mode  f[THz]  pol   |g_CB->CB|[meV]  Xi[eV]   D[eV/A]  |g_CB->CB+1|[meV]
     1    2.966   Z          0.000     0.000     0.000         2.625
     2    3.999   T          6.154     0.343     0.219         0.000
     3    5.389   L         27.852     1.804     1.151         0.000
     4    8.783   T          0.000     0.000     0.000         6.343
     5    9.544   L          0.000     0.000     0.000       104.583
     6   10.311   L         25.162     2.254     1.438         0.000
     7   11.000   T          9.762     0.903     0.576         0.000
     8   11.915   Z         16.795     1.618     1.032         0.000
     9   12.118   Z          0.000     0.000     0.000        58.970

Xi_LA(q) [eV] from the small-q points: |q|=0.118/A: 2.70, |q|=0.118/A: 3.22, |q|=0.237/A: 2.95, |q|=0.237/A: 1.24
  -> Xi_LA: q -> 0 extrapolation 3.82 eV, mean over the window 2.53 eV   [Kaasbjerg: LA 2.8, TA 1.6]

Xi_TA(q) [eV] from the small-q points: |q|=0.118/A: 3.37, |q|=0.118/A: 0.29, |q|=0.237/A: 1.78, |q|=0.237/A: 0.59
  -> Xi_TA: q -> 0 extrapolation 2.47 eV, mean over the window 1.51 eV   [Kaasbjerg: LA 2.8, TA 1.6]
