
-----------------------------
------- calculate fc2 -------
-----------------------------

        _
  _ __ | |__   ___  _ __   ___   _ __  _   _
 | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | |
 | |_) | | | | (_) | | | | (_) || |_) | |_| |
 | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, |
 |_|                            |_|    |___/
                                      2.47.1

-------------------------[time 2026-01-08 10:21:19]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phonopy.load mode.
Python version 3.14.2
Spglib version 2.6.1

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [2 2 1]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
Number of symmetry operations in supercell: 96
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   10.921106497213408    0.000000000000000    0.000000000000000
  b   -5.460553248606704    9.457955664022098    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
   *1 Cs  0.20929583065401  0.79070416934600  0.45863866870204 132.905
    2 Cs  0.79070416934600  0.58140833869199  0.95863866870204 132.905
    3 Cs  0.41859166130801  0.20929583065401  0.95863866870204 132.905
    4 Cs  0.58140833869199  0.79070416934600  0.45863866870204 132.905
    5 Cs  0.20929583065401  0.41859166130801  0.45863866870204 132.905
    6 Cs  0.79070416934600  0.20929583065401  0.54136133129796 132.905
    7 Cs  0.79070416934600  0.58140833869199  0.54136133129796 132.905
    8 Cs  0.79070416934600  0.20929583065401  0.95863866870204 132.905
    9 Cs  0.41859166130801  0.20929583065401  0.54136133129796 132.905
   10 Cs  0.58140833869199  0.79070416934600  0.04136133129796 132.905
   11 Cs  0.20929583065401  0.79070416934600  0.04136133129796 132.905
   12 Cs  0.20929583065401  0.41859166130801  0.04136133129796 132.905
  *13 Ga  0.33333333333333  0.66666666666667  0.75000000000000  69.723
   14 Ga  0.66666666666667  0.33333333333333  0.25000000000000  69.723
  *15 Sb  0.47341102757325  0.94682205514651  0.75000000000000 121.760
   16 Sb  0.94682205514651  0.47341102757325  0.25000000000000 121.760
   17 Sb  0.52658897242674  0.47341102757325  0.25000000000000 121.760
   18 Sb  0.47341102757325  0.52658897242674  0.75000000000000 121.760
   19 Sb  0.05317794485349  0.52658897242674  0.75000000000000 121.760
   20 Sb  0.52658897242674  0.05317794485349  0.25000000000000 121.760
  *21 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760
   22 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760
  *23 K   0.11904674310150  0.88095325689850  0.75000000000000  39.098
   24 K   0.88095325689850  0.76190651379700  0.25000000000000  39.098
   25 K   0.23809348620300  0.11904674310150  0.25000000000000  39.098
   26 K   0.76190651379700  0.88095325689850  0.75000000000000  39.098
   27 K   0.11904674310150  0.23809348620300  0.75000000000000  39.098
   28 K   0.88095325689850  0.11904674310150  0.25000000000000  39.098
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a   10.921106497213408    0.000000000000000    0.000000000000000
  b   -5.460553248606704    9.457955664022098    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
   *1 Cs  0.20929583065401  0.79070416934600  0.45863866870204 132.905 > 1
    2 Cs  0.79070416934600  0.58140833869199  0.95863866870204 132.905 > 2
    3 Cs  0.41859166130801  0.20929583065401  0.95863866870204 132.905 > 3
    4 Cs  0.58140833869199  0.79070416934600  0.45863866870204 132.905 > 4
    5 Cs  0.20929583065401  0.41859166130801  0.45863866870204 132.905 > 5
    6 Cs  0.79070416934600  0.20929583065401  0.54136133129796 132.905 > 6
    7 Cs  0.79070416934600  0.58140833869199  0.54136133129796 132.905 > 7
    8 Cs  0.79070416934600  0.20929583065401  0.95863866870204 132.905 > 8
    9 Cs  0.41859166130801  0.20929583065401  0.54136133129796 132.905 > 9
   10 Cs  0.58140833869199  0.79070416934600  0.04136133129796 132.905 > 10
   11 Cs  0.20929583065401  0.79070416934600  0.04136133129796 132.905 > 11
   12 Cs  0.20929583065401  0.41859166130801  0.04136133129796 132.905 > 12
  *13 Ga  0.33333333333333  0.66666666666667  0.75000000000000  69.723 > 13
   14 Ga  0.66666666666667  0.33333333333333  0.25000000000000  69.723 > 14
  *15 Sb  0.47341102757325  0.94682205514651  0.75000000000000 121.760 > 15
   16 Sb  0.94682205514651  0.47341102757325  0.25000000000000 121.760 > 16
   17 Sb  0.52658897242674  0.47341102757325  0.25000000000000 121.760 > 17
   18 Sb  0.47341102757325  0.52658897242674  0.75000000000000 121.760 > 18
   19 Sb  0.05317794485349  0.52658897242674  0.75000000000000 121.760 > 19
   20 Sb  0.52658897242674  0.05317794485349  0.25000000000000 121.760 > 20
  *21 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760 > 21
   22 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760 > 22
  *23 K   0.11904674310150  0.88095325689850  0.75000000000000  39.098 > 23
   24 K   0.88095325689850  0.76190651379700  0.25000000000000  39.098 > 24
   25 K   0.23809348620300  0.11904674310150  0.25000000000000  39.098 > 25
   26 K   0.76190651379700  0.88095325689850  0.75000000000000  39.098 > 26
   27 K   0.11904674310150  0.23809348620300  0.75000000000000  39.098 > 27
   28 K   0.88095325689850  0.11904674310150  0.25000000000000  39.098 > 28
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   21.842212994426816    0.000000000000000    0.000000000000000
  b  -10.921106497213408   18.915911328044196    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
   *1 Cs  0.10464791532700  0.39535208467300  0.45863866870204 132.905 > 1
    2 Cs  0.60464791532700  0.39535208467300  0.45863866870204 132.905 > 1
    3 Cs  0.10464791532700  0.89535208467300  0.45863866870204 132.905 > 1
    4 Cs  0.60464791532700  0.89535208467300  0.45863866870204 132.905 > 1
    5 Cs  0.39535208467300  0.29070416934600  0.95863866870204 132.905 > 2
    6 Cs  0.89535208467300  0.29070416934600  0.95863866870204 132.905 > 2
    7 Cs  0.39535208467300  0.79070416934600  0.95863866870204 132.905 > 2
    8 Cs  0.89535208467300  0.79070416934600  0.95863866870204 132.905 > 2
    9 Cs  0.20929583065400  0.10464791532700  0.95863866870204 132.905 > 3
   10 Cs  0.70929583065400  0.10464791532700  0.95863866870204 132.905 > 3
   11 Cs  0.20929583065400  0.60464791532700  0.95863866870204 132.905 > 3
   12 Cs  0.70929583065400  0.60464791532700  0.95863866870204 132.905 > 3
   13 Cs  0.29070416934600  0.39535208467300  0.45863866870204 132.905 > 4
   14 Cs  0.79070416934600  0.39535208467300  0.45863866870204 132.905 > 4
   15 Cs  0.29070416934600  0.89535208467300  0.45863866870204 132.905 > 4
   16 Cs  0.79070416934600  0.89535208467300  0.45863866870204 132.905 > 4
   17 Cs  0.10464791532700  0.20929583065400  0.45863866870204 132.905 > 5
   18 Cs  0.60464791532700  0.20929583065400  0.45863866870204 132.905 > 5
   19 Cs  0.10464791532700  0.70929583065400  0.45863866870204 132.905 > 5
   20 Cs  0.60464791532700  0.70929583065400  0.45863866870204 132.905 > 5
   21 Cs  0.39535208467300  0.10464791532700  0.54136133129796 132.905 > 6
   22 Cs  0.89535208467300  0.10464791532700  0.54136133129796 132.905 > 6
   23 Cs  0.39535208467300  0.60464791532700  0.54136133129796 132.905 > 6
   24 Cs  0.89535208467300  0.60464791532700  0.54136133129796 132.905 > 6
   25 Cs  0.39535208467300  0.29070416934600  0.54136133129796 132.905 > 7
   26 Cs  0.89535208467300  0.29070416934600  0.54136133129796 132.905 > 7
   27 Cs  0.39535208467300  0.79070416934600  0.54136133129796 132.905 > 7
   28 Cs  0.89535208467300  0.79070416934600  0.54136133129796 132.905 > 7
   29 Cs  0.39535208467300  0.10464791532700  0.95863866870204 132.905 > 8
   30 Cs  0.89535208467300  0.10464791532700  0.95863866870204 132.905 > 8
   31 Cs  0.39535208467300  0.60464791532700  0.95863866870204 132.905 > 8
   32 Cs  0.89535208467300  0.60464791532700  0.95863866870204 132.905 > 8
   33 Cs  0.20929583065400  0.10464791532700  0.54136133129796 132.905 > 9
   34 Cs  0.70929583065400  0.10464791532700  0.54136133129796 132.905 > 9
   35 Cs  0.20929583065400  0.60464791532700  0.54136133129796 132.905 > 9
   36 Cs  0.70929583065400  0.60464791532700  0.54136133129796 132.905 > 9
   37 Cs  0.29070416934600  0.39535208467300  0.04136133129796 132.905 > 10
   38 Cs  0.79070416934600  0.39535208467300  0.04136133129796 132.905 > 10
   39 Cs  0.29070416934600  0.89535208467300  0.04136133129796 132.905 > 10
   40 Cs  0.79070416934600  0.89535208467300  0.04136133129796 132.905 > 10
   41 Cs  0.10464791532700  0.39535208467300  0.04136133129796 132.905 > 11
   42 Cs  0.60464791532700  0.39535208467300  0.04136133129796 132.905 > 11
   43 Cs  0.10464791532700  0.89535208467300  0.04136133129796 132.905 > 11
   44 Cs  0.60464791532700  0.89535208467300  0.04136133129796 132.905 > 11
   45 Cs  0.10464791532700  0.20929583065400  0.04136133129796 132.905 > 12
   46 Cs  0.60464791532700  0.20929583065400  0.04136133129796 132.905 > 12
   47 Cs  0.10464791532700  0.70929583065400  0.04136133129796 132.905 > 12
   48 Cs  0.60464791532700  0.70929583065400  0.04136133129796 132.905 > 12
  *49 Ga  0.16666666666667  0.33333333333333  0.75000000000000  69.723 > 13
   50 Ga  0.66666666666667  0.33333333333333  0.75000000000000  69.723 > 13
   51 Ga  0.16666666666667  0.83333333333333  0.75000000000000  69.723 > 13
   52 Ga  0.66666666666667  0.83333333333333  0.75000000000000  69.723 > 13
   53 Ga  0.33333333333333  0.16666666666667  0.25000000000000  69.723 > 14
   54 Ga  0.83333333333333  0.16666666666667  0.25000000000000  69.723 > 14
   55 Ga  0.33333333333333  0.66666666666667  0.25000000000000  69.723 > 14
   56 Ga  0.83333333333333  0.66666666666667  0.25000000000000  69.723 > 14
  *57 Sb  0.23670551378663  0.47341102757325  0.75000000000000 121.760 > 15
   58 Sb  0.73670551378663  0.47341102757325  0.75000000000000 121.760 > 15
   59 Sb  0.23670551378663  0.97341102757325  0.75000000000000 121.760 > 15
   60 Sb  0.73670551378663  0.97341102757325  0.75000000000000 121.760 > 15
   61 Sb  0.47341102757325  0.23670551378663  0.25000000000000 121.760 > 16
   62 Sb  0.97341102757325  0.23670551378663  0.25000000000000 121.760 > 16
   63 Sb  0.47341102757325  0.73670551378663  0.25000000000000 121.760 > 16
   64 Sb  0.97341102757325  0.73670551378663  0.25000000000000 121.760 > 16
   65 Sb  0.26329448621337  0.23670551378663  0.25000000000000 121.760 > 17
   66 Sb  0.76329448621337  0.23670551378663  0.25000000000000 121.760 > 17
   67 Sb  0.26329448621337  0.73670551378663  0.25000000000000 121.760 > 17
   68 Sb  0.76329448621337  0.73670551378663  0.25000000000000 121.760 > 17
   69 Sb  0.23670551378663  0.26329448621337  0.75000000000000 121.760 > 18
   70 Sb  0.73670551378663  0.26329448621337  0.75000000000000 121.760 > 18
   71 Sb  0.23670551378663  0.76329448621337  0.75000000000000 121.760 > 18
   72 Sb  0.73670551378663  0.76329448621337  0.75000000000000 121.760 > 18
   73 Sb  0.02658897242675  0.26329448621337  0.75000000000000 121.760 > 19
   74 Sb  0.52658897242675  0.26329448621337  0.75000000000000 121.760 > 19
   75 Sb  0.02658897242675  0.76329448621337  0.75000000000000 121.760 > 19
   76 Sb  0.52658897242675  0.76329448621337  0.75000000000000 121.760 > 19
   77 Sb  0.26329448621337  0.02658897242675  0.25000000000000 121.760 > 20
   78 Sb  0.76329448621337  0.02658897242675  0.25000000000000 121.760 > 20
   79 Sb  0.26329448621337  0.52658897242675  0.25000000000000 121.760 > 20
   80 Sb  0.76329448621337  0.52658897242675  0.25000000000000 121.760 > 20
  *81 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760 > 21
   82 Sb  0.50000000000000  0.00000000000000  0.00000000000000 121.760 > 21
   83 Sb  0.00000000000000  0.50000000000000  0.00000000000000 121.760 > 21
   84 Sb  0.50000000000000  0.50000000000000  0.00000000000000 121.760 > 21
   85 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760 > 22
   86 Sb  0.50000000000000  0.00000000000000  0.50000000000000 121.760 > 22
   87 Sb  0.00000000000000  0.50000000000000  0.50000000000000 121.760 > 22
   88 Sb  0.50000000000000  0.50000000000000  0.50000000000000 121.760 > 22
  *89 K   0.05952337155075  0.44047662844925  0.75000000000000  39.098 > 23
   90 K   0.55952337155075  0.44047662844925  0.75000000000000  39.098 > 23
   91 K   0.05952337155075  0.94047662844925  0.75000000000000  39.098 > 23
   92 K   0.55952337155075  0.94047662844925  0.75000000000000  39.098 > 23
   93 K   0.44047662844925  0.38095325689850  0.25000000000000  39.098 > 24
   94 K   0.94047662844925  0.38095325689850  0.25000000000000  39.098 > 24
   95 K   0.44047662844925  0.88095325689850  0.25000000000000  39.098 > 24
   96 K   0.94047662844925  0.88095325689850  0.25000000000000  39.098 > 24
   97 K   0.11904674310150  0.05952337155075  0.25000000000000  39.098 > 25
   98 K   0.61904674310150  0.05952337155075  0.25000000000000  39.098 > 25
   99 K   0.11904674310150  0.55952337155075  0.25000000000000  39.098 > 25
  100 K   0.61904674310150  0.55952337155075  0.25000000000000  39.098 > 25
  101 K   0.38095325689850  0.44047662844925  0.75000000000000  39.098 > 26
  102 K   0.88095325689850  0.44047662844925  0.75000000000000  39.098 > 26
  103 K   0.38095325689850  0.94047662844925  0.75000000000000  39.098 > 26
  104 K   0.88095325689850  0.94047662844925  0.75000000000000  39.098 > 26
  105 K   0.05952337155075  0.11904674310150  0.75000000000000  39.098 > 27
  106 K   0.55952337155075  0.11904674310150  0.75000000000000  39.098 > 27
  107 K   0.05952337155075  0.61904674310150  0.75000000000000  39.098 > 27
  108 K   0.55952337155075  0.61904674310150  0.75000000000000  39.098 > 27
  109 K   0.44047662844925  0.05952337155075  0.25000000000000  39.098 > 28
  110 K   0.94047662844925  0.05952337155075  0.25000000000000  39.098 > 28
  111 K   0.44047662844925  0.55952337155075  0.25000000000000  39.098 > 28
  112 K   0.94047662844925  0.55952337155075  0.25000000000000  39.098 > 28
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            9.9027445    0.0000000    0.0000000
           -0.0000000    9.9027445    0.0000000
            0.0000000    0.0000000    6.0926774
-------------------------- Born effective charges --------------------------
    1 Cs    1.1405182    0.0439356    0.2763771
            0.0439356    1.1912507   -0.1595664
            0.0811556   -0.0468552    1.1503096
    2 Cs    1.2166170    0.0000000    0.0000000
           -0.0000000    1.1151519   -0.3191328
           -0.0000000   -0.0937104    1.1503096
    3 Cs    1.1405182   -0.0439356    0.2763771
           -0.0439356    1.1912507    0.1595664
            0.0811556    0.0468552    1.1503096
    4 Cs    1.1405182   -0.0439356   -0.2763771
           -0.0439356    1.1912507   -0.1595664
           -0.0811556   -0.0468552    1.1503096
    5 Cs    1.2166170    0.0000000   -0.0000000
           -0.0000000    1.1151519    0.3191328
           -0.0000000    0.0937104    1.1503096
    6 Cs    1.1405182    0.0439356    0.2763771
            0.0439356    1.1912507   -0.1595664
            0.0811556   -0.0468552    1.1503096
    7 Cs    1.2166170    0.0000000   -0.0000000
           -0.0000000    1.1151519    0.3191328
           -0.0000000    0.0937104    1.1503096
    8 Cs    1.1405182    0.0439356   -0.2763771
            0.0439356    1.1912507    0.1595664
           -0.0811556    0.0468552    1.1503096
    9 Cs    1.1405182   -0.0439356   -0.2763771
           -0.0439356    1.1912507   -0.1595664
           -0.0811556   -0.0468552    1.1503096
   10 Cs    1.1405182   -0.0439356    0.2763771
           -0.0439356    1.1912507    0.1595664
            0.0811556    0.0468552    1.1503096
   11 Cs    1.1405182    0.0439356   -0.2763771
            0.0439356    1.1912507    0.1595664
           -0.0811556    0.0468552    1.1503096
   12 Cs    1.2166170    0.0000000    0.0000000
           -0.0000000    1.1151519   -0.3191328
           -0.0000000   -0.0937104    1.1503096
   13 Ga    2.7985137    0.0000000    0.0000000
           -0.0000000    2.7985137    0.0000000
           -0.0000000   -0.0000000   -1.0476192
   14 Ga    2.7985137    0.0000000    0.0000000
           -0.0000000    2.7985137    0.0000000
           -0.0000000   -0.0000000   -1.0476192
   15 Sb   -2.2548971    0.0000000    0.0000000
            0.0000000   -4.1494014    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   16 Sb   -3.6757753   -0.8203444    0.0000000
           -0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   17 Sb   -3.6757753    0.8203444    0.0000000
            0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   18 Sb   -3.6757753    0.8203444    0.0000000
            0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   19 Sb   -3.6757753   -0.8203444    0.0000000
           -0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   20 Sb   -2.2548971    0.0000000    0.0000000
            0.0000000   -4.1494014    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   21 Sb   -3.0985017   -0.0000000    0.0000000
            0.0000000   -3.0985017    0.0000000
           -0.0000000   -0.0000000   -2.5818845
   22 Sb   -3.0985017   -0.0000000    0.0000000
            0.0000000   -3.0985017    0.0000000
           -0.0000000   -0.0000000   -2.5818845
   23 K     1.1474204   -0.3066493    0.0000000
           -0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   24 K     0.6162882   -0.0000000    0.0000000
           -0.0000000    1.3244644    0.0000000
           -0.0000000   -0.0000000    0.9490339
   25 K     1.1474204    0.3066493    0.0000000
            0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   26 K     1.1474204    0.3066493    0.0000000
            0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   27 K     0.6162882   -0.0000000    0.0000000
           -0.0000000    1.3244644    0.0000000
           -0.0000000   -0.0000000    0.9490339
   28 K     1.1474204   -0.3066493    0.0000000
           -0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
----------------------------------------------------------------------------
Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
-------------------------------- Symfc start -------------------------------
Symfc version 1.5.4 (https://github.com/symfc/symfc)
Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024)
Computing [2] order force constants.
Permutation basis: 336/336
Permutation basis: 14448/14448
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 701
Number of blocks in projector: 388
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 375
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 196
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 130
Use standard eigh solver.
Tree of FC basis block matrices:
- (701, 686), data: False
|-- (130, 125), data: True
|-- (196, 194), data: True
|-- (375, 367), data: True
-----
Solver_atoms: 1 -- 112 / 112
Time (Solver_compr_matrix_reshape): 0.008
Solver_block: 100 / 160
 - Time: 0.646
Solver_block: 160 / 160
 - Time: 0.425
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.129
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 336/336
Permutation basis: 14448/14448
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 701
Number of blocks in projector: 388
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 375
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 196
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 130
Use standard eigh solver.
Tree of FC basis block matrices:
- (701, 686), data: False
|-- (130, 125), data: True
|-- (196, 194), data: True
|-- (375, 367), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) 
Force constants are written into "force_constants.hdf5".

----------------------------------------------------------------------------
 One of the following run modes may be specified for phonon calculations.
 - Mesh sampling (MESH, --mesh)
 - Q-points (QPOINTS, --qpoints)
 - Band structure (BAND, --band)
 - Animation (ANIME, --anime)
 - Modulation (MODULATION, --modulation)
 - Characters of Irreps (IRREPS, --irreps)
 - Create displacements (CREATE_DISPLACEMENTS, -d)
----------------------------------------------------------------------------

Summary of calculation was written in "phonopy.yaml".
-------------------------[time 2026-01-08 10:21:31]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


----------------------------
------ calculate fc3 -------
----------------------------

        _                      _____
  _ __ | |__   ___  _ __   ___|___ / _ __  _   _
 | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | |
 | |_) | | | | (_) | | | | (_) |__) | |_) | |_| |
 | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, |
 |_|                                |_|    |___/ 
                                      3.23.0

-------------------------[time 2026-01-08 10:21:32]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [2 2 1]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   10.921106497213408    0.000000000000000    0.000000000000000
  b   -5.460553248606704    9.457955664022098    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
    1 Cs  0.20929583065401  0.79070416934600  0.45863866870204 132.905
    2 Cs  0.79070416934600  0.58140833869199  0.95863866870204 132.905
    3 Cs  0.41859166130801  0.20929583065401  0.95863866870204 132.905
    4 Cs  0.58140833869199  0.79070416934600  0.45863866870204 132.905
    5 Cs  0.20929583065401  0.41859166130801  0.45863866870204 132.905
    6 Cs  0.79070416934600  0.20929583065401  0.54136133129796 132.905
    7 Cs  0.79070416934600  0.58140833869199  0.54136133129796 132.905
    8 Cs  0.79070416934600  0.20929583065401  0.95863866870204 132.905
    9 Cs  0.41859166130801  0.20929583065401  0.54136133129796 132.905
   10 Cs  0.58140833869199  0.79070416934600  0.04136133129796 132.905
   11 Cs  0.20929583065401  0.79070416934600  0.04136133129796 132.905
   12 Cs  0.20929583065401  0.41859166130801  0.04136133129796 132.905
   13 Ga  0.33333333333333  0.66666666666667  0.75000000000000  69.723
   14 Ga  0.66666666666667  0.33333333333333  0.25000000000000  69.723
   15 Sb  0.47341102757325  0.94682205514651  0.75000000000000 121.760
   16 Sb  0.94682205514651  0.47341102757325  0.25000000000000 121.760
   17 Sb  0.52658897242674  0.47341102757325  0.25000000000000 121.760
   18 Sb  0.47341102757325  0.52658897242674  0.75000000000000 121.760
   19 Sb  0.05317794485349  0.52658897242674  0.75000000000000 121.760
   20 Sb  0.52658897242674  0.05317794485349  0.25000000000000 121.760
   21 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760
   22 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760
   23 K   0.11904674310150  0.88095325689850  0.75000000000000  39.098
   24 K   0.88095325689850  0.76190651379700  0.25000000000000  39.098
   25 K   0.23809348620300  0.11904674310150  0.25000000000000  39.098
   26 K   0.76190651379700  0.88095325689850  0.75000000000000  39.098
   27 K   0.11904674310150  0.23809348620300  0.75000000000000  39.098
   28 K   0.88095325689850  0.11904674310150  0.25000000000000  39.098
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   21.842212994426816    0.000000000000000    0.000000000000000
  b  -10.921106497213408   18.915911328044196    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
    1 Cs  0.10464791532700  0.39535208467300  0.45863866870204 132.905 > 1
    2 Cs  0.60464791532700  0.39535208467300  0.45863866870204 132.905 > 1
    3 Cs  0.10464791532700  0.89535208467300  0.45863866870204 132.905 > 1
    4 Cs  0.60464791532700  0.89535208467300  0.45863866870204 132.905 > 1
    5 Cs  0.39535208467300  0.29070416934600  0.95863866870204 132.905 > 5
    6 Cs  0.89535208467300  0.29070416934600  0.95863866870204 132.905 > 5
    7 Cs  0.39535208467300  0.79070416934600  0.95863866870204 132.905 > 5
    8 Cs  0.89535208467300  0.79070416934600  0.95863866870204 132.905 > 5
    9 Cs  0.20929583065400  0.10464791532700  0.95863866870204 132.905 > 9
   10 Cs  0.70929583065400  0.10464791532700  0.95863866870204 132.905 > 9
   11 Cs  0.20929583065400  0.60464791532700  0.95863866870204 132.905 > 9
   12 Cs  0.70929583065400  0.60464791532700  0.95863866870204 132.905 > 9
   13 Cs  0.29070416934600  0.39535208467300  0.45863866870204 132.905 > 13
   14 Cs  0.79070416934600  0.39535208467300  0.45863866870204 132.905 > 13
   15 Cs  0.29070416934600  0.89535208467300  0.45863866870204 132.905 > 13
   16 Cs  0.79070416934600  0.89535208467300  0.45863866870204 132.905 > 13
   17 Cs  0.10464791532700  0.20929583065400  0.45863866870204 132.905 > 17
   18 Cs  0.60464791532700  0.20929583065400  0.45863866870204 132.905 > 17
   19 Cs  0.10464791532700  0.70929583065400  0.45863866870204 132.905 > 17
   20 Cs  0.60464791532700  0.70929583065400  0.45863866870204 132.905 > 17
   21 Cs  0.39535208467300  0.10464791532700  0.54136133129796 132.905 > 21
   22 Cs  0.89535208467300  0.10464791532700  0.54136133129796 132.905 > 21
   23 Cs  0.39535208467300  0.60464791532700  0.54136133129796 132.905 > 21
   24 Cs  0.89535208467300  0.60464791532700  0.54136133129796 132.905 > 21
   25 Cs  0.39535208467300  0.29070416934600  0.54136133129796 132.905 > 25
   26 Cs  0.89535208467300  0.29070416934600  0.54136133129796 132.905 > 25
   27 Cs  0.39535208467300  0.79070416934600  0.54136133129796 132.905 > 25
   28 Cs  0.89535208467300  0.79070416934600  0.54136133129796 132.905 > 25
   29 Cs  0.39535208467300  0.10464791532700  0.95863866870204 132.905 > 29
   30 Cs  0.89535208467300  0.10464791532700  0.95863866870204 132.905 > 29
   31 Cs  0.39535208467300  0.60464791532700  0.95863866870204 132.905 > 29
   32 Cs  0.89535208467300  0.60464791532700  0.95863866870204 132.905 > 29
   33 Cs  0.20929583065400  0.10464791532700  0.54136133129796 132.905 > 33
   34 Cs  0.70929583065400  0.10464791532700  0.54136133129796 132.905 > 33
   35 Cs  0.20929583065400  0.60464791532700  0.54136133129796 132.905 > 33
   36 Cs  0.70929583065400  0.60464791532700  0.54136133129796 132.905 > 33
   37 Cs  0.29070416934600  0.39535208467300  0.04136133129796 132.905 > 37
   38 Cs  0.79070416934600  0.39535208467300  0.04136133129796 132.905 > 37
   39 Cs  0.29070416934600  0.89535208467300  0.04136133129796 132.905 > 37
   40 Cs  0.79070416934600  0.89535208467300  0.04136133129796 132.905 > 37
   41 Cs  0.10464791532700  0.39535208467300  0.04136133129796 132.905 > 41
   42 Cs  0.60464791532700  0.39535208467300  0.04136133129796 132.905 > 41
   43 Cs  0.10464791532700  0.89535208467300  0.04136133129796 132.905 > 41
   44 Cs  0.60464791532700  0.89535208467300  0.04136133129796 132.905 > 41
   45 Cs  0.10464791532700  0.20929583065400  0.04136133129796 132.905 > 45
   46 Cs  0.60464791532700  0.20929583065400  0.04136133129796 132.905 > 45
   47 Cs  0.10464791532700  0.70929583065400  0.04136133129796 132.905 > 45
   48 Cs  0.60464791532700  0.70929583065400  0.04136133129796 132.905 > 45
   49 Ga  0.16666666666667  0.33333333333333  0.75000000000000  69.723 > 49
   50 Ga  0.66666666666667  0.33333333333333  0.75000000000000  69.723 > 49
   51 Ga  0.16666666666667  0.83333333333333  0.75000000000000  69.723 > 49
   52 Ga  0.66666666666667  0.83333333333333  0.75000000000000  69.723 > 49
   53 Ga  0.33333333333333  0.16666666666667  0.25000000000000  69.723 > 53
   54 Ga  0.83333333333333  0.16666666666667  0.25000000000000  69.723 > 53
   55 Ga  0.33333333333333  0.66666666666667  0.25000000000000  69.723 > 53
   56 Ga  0.83333333333333  0.66666666666667  0.25000000000000  69.723 > 53
   57 Sb  0.23670551378663  0.47341102757325  0.75000000000000 121.760 > 57
   58 Sb  0.73670551378663  0.47341102757325  0.75000000000000 121.760 > 57
   59 Sb  0.23670551378663  0.97341102757325  0.75000000000000 121.760 > 57
   60 Sb  0.73670551378663  0.97341102757325  0.75000000000000 121.760 > 57
   61 Sb  0.47341102757325  0.23670551378663  0.25000000000000 121.760 > 61
   62 Sb  0.97341102757325  0.23670551378663  0.25000000000000 121.760 > 61
   63 Sb  0.47341102757325  0.73670551378663  0.25000000000000 121.760 > 61
   64 Sb  0.97341102757325  0.73670551378663  0.25000000000000 121.760 > 61
   65 Sb  0.26329448621337  0.23670551378663  0.25000000000000 121.760 > 65
   66 Sb  0.76329448621337  0.23670551378663  0.25000000000000 121.760 > 65
   67 Sb  0.26329448621337  0.73670551378663  0.25000000000000 121.760 > 65
   68 Sb  0.76329448621337  0.73670551378663  0.25000000000000 121.760 > 65
   69 Sb  0.23670551378663  0.26329448621337  0.75000000000000 121.760 > 69
   70 Sb  0.73670551378663  0.26329448621337  0.75000000000000 121.760 > 69
   71 Sb  0.23670551378663  0.76329448621337  0.75000000000000 121.760 > 69
   72 Sb  0.73670551378663  0.76329448621337  0.75000000000000 121.760 > 69
   73 Sb  0.02658897242675  0.26329448621337  0.75000000000000 121.760 > 73
   74 Sb  0.52658897242675  0.26329448621337  0.75000000000000 121.760 > 73
   75 Sb  0.02658897242675  0.76329448621337  0.75000000000000 121.760 > 73
   76 Sb  0.52658897242675  0.76329448621337  0.75000000000000 121.760 > 73
   77 Sb  0.26329448621337  0.02658897242675  0.25000000000000 121.760 > 77
   78 Sb  0.76329448621337  0.02658897242675  0.25000000000000 121.760 > 77
   79 Sb  0.26329448621337  0.52658897242675  0.25000000000000 121.760 > 77
   80 Sb  0.76329448621337  0.52658897242675  0.25000000000000 121.760 > 77
   81 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760 > 81
   82 Sb  0.50000000000000  0.00000000000000  0.00000000000000 121.760 > 81
   83 Sb  0.00000000000000  0.50000000000000  0.00000000000000 121.760 > 81
   84 Sb  0.50000000000000  0.50000000000000  0.00000000000000 121.760 > 81
   85 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760 > 85
   86 Sb  0.50000000000000  0.00000000000000  0.50000000000000 121.760 > 85
   87 Sb  0.00000000000000  0.50000000000000  0.50000000000000 121.760 > 85
   88 Sb  0.50000000000000  0.50000000000000  0.50000000000000 121.760 > 85
   89 K   0.05952337155075  0.44047662844925  0.75000000000000  39.098 > 89
   90 K   0.55952337155075  0.44047662844925  0.75000000000000  39.098 > 89
   91 K   0.05952337155075  0.94047662844925  0.75000000000000  39.098 > 89
   92 K   0.55952337155075  0.94047662844925  0.75000000000000  39.098 > 89
   93 K   0.44047662844925  0.38095325689850  0.25000000000000  39.098 > 93
   94 K   0.94047662844925  0.38095325689850  0.25000000000000  39.098 > 93
   95 K   0.44047662844925  0.88095325689850  0.25000000000000  39.098 > 93
   96 K   0.94047662844925  0.88095325689850  0.25000000000000  39.098 > 93
   97 K   0.11904674310150  0.05952337155075  0.25000000000000  39.098 > 97
   98 K   0.61904674310150  0.05952337155075  0.25000000000000  39.098 > 97
   99 K   0.11904674310150  0.55952337155075  0.25000000000000  39.098 > 97
  100 K   0.61904674310150  0.55952337155075  0.25000000000000  39.098 > 97
  101 K   0.38095325689850  0.44047662844925  0.75000000000000  39.098 > 101
  102 K   0.88095325689850  0.44047662844925  0.75000000000000  39.098 > 101
  103 K   0.38095325689850  0.94047662844925  0.75000000000000  39.098 > 101
  104 K   0.88095325689850  0.94047662844925  0.75000000000000  39.098 > 101
  105 K   0.05952337155075  0.11904674310150  0.75000000000000  39.098 > 105
  106 K   0.55952337155075  0.11904674310150  0.75000000000000  39.098 > 105
  107 K   0.05952337155075  0.61904674310150  0.75000000000000  39.098 > 105
  108 K   0.55952337155075  0.61904674310150  0.75000000000000  39.098 > 105
  109 K   0.44047662844925  0.05952337155075  0.25000000000000  39.098 > 109
  110 K   0.94047662844925  0.05952337155075  0.25000000000000  39.098 > 109
  111 K   0.44047662844925  0.55952337155075  0.25000000000000  39.098 > 109
  112 K   0.94047662844925  0.55952337155075  0.25000000000000  39.098 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            9.9027445    0.0000000    0.0000000
           -0.0000000    9.9027445    0.0000000
            0.0000000    0.0000000    6.0926774
-------------------------- Born effective charges --------------------------
    1 Cs    1.1405182    0.0439356    0.2763771
            0.0439356    1.1912507   -0.1595664
            0.0811556   -0.0468552    1.1503096
    2 Cs    1.2166170    0.0000000    0.0000000
           -0.0000000    1.1151519   -0.3191328
           -0.0000000   -0.0937104    1.1503096
    3 Cs    1.1405182   -0.0439356    0.2763771
           -0.0439356    1.1912507    0.1595664
            0.0811556    0.0468552    1.1503096
    4 Cs    1.1405182   -0.0439356   -0.2763771
           -0.0439356    1.1912507   -0.1595664
           -0.0811556   -0.0468552    1.1503096
    5 Cs    1.2166170    0.0000000   -0.0000000
           -0.0000000    1.1151519    0.3191328
           -0.0000000    0.0937104    1.1503096
    6 Cs    1.1405182    0.0439356    0.2763771
            0.0439356    1.1912507   -0.1595664
            0.0811556   -0.0468552    1.1503096
    7 Cs    1.2166170    0.0000000   -0.0000000
           -0.0000000    1.1151519    0.3191328
           -0.0000000    0.0937104    1.1503096
    8 Cs    1.1405182    0.0439356   -0.2763771
            0.0439356    1.1912507    0.1595664
           -0.0811556    0.0468552    1.1503096
    9 Cs    1.1405182   -0.0439356   -0.2763771
           -0.0439356    1.1912507   -0.1595664
           -0.0811556   -0.0468552    1.1503096
   10 Cs    1.1405182   -0.0439356    0.2763771
           -0.0439356    1.1912507    0.1595664
            0.0811556    0.0468552    1.1503096
   11 Cs    1.1405182    0.0439356   -0.2763771
            0.0439356    1.1912507    0.1595664
           -0.0811556    0.0468552    1.1503096
   12 Cs    1.2166170    0.0000000    0.0000000
           -0.0000000    1.1151519   -0.3191328
           -0.0000000   -0.0937104    1.1503096
   13 Ga    2.7985137    0.0000000    0.0000000
           -0.0000000    2.7985137    0.0000000
           -0.0000000   -0.0000000   -1.0476192
   14 Ga    2.7985137    0.0000000    0.0000000
           -0.0000000    2.7985137    0.0000000
           -0.0000000   -0.0000000   -1.0476192
   15 Sb   -2.2548971    0.0000000    0.0000000
            0.0000000   -4.1494014    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   16 Sb   -3.6757753   -0.8203444    0.0000000
           -0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   17 Sb   -3.6757753    0.8203444    0.0000000
            0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   18 Sb   -3.6757753    0.8203444    0.0000000
            0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   19 Sb   -3.6757753   -0.8203444    0.0000000
           -0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   20 Sb   -2.2548971    0.0000000    0.0000000
            0.0000000   -4.1494014    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   21 Sb   -3.0985017   -0.0000000    0.0000000
            0.0000000   -3.0985017    0.0000000
           -0.0000000   -0.0000000   -2.5818845
   22 Sb   -3.0985017   -0.0000000    0.0000000
            0.0000000   -3.0985017    0.0000000
           -0.0000000   -0.0000000   -2.5818845
   23 K     1.1474204   -0.3066493    0.0000000
           -0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   24 K     0.6162882   -0.0000000    0.0000000
           -0.0000000    1.3244644    0.0000000
           -0.0000000   -0.0000000    0.9490339
   25 K     1.1474204    0.3066493    0.0000000
            0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   26 K     1.1474204    0.3066493    0.0000000
            0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   27 K     0.6162882   -0.0000000    0.0000000
           -0.0000000    1.3244644    0.0000000
           -0.0000000   -0.0000000    0.9490339
   28 K     1.1474204   -0.3066493    0.0000000
           -0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 49, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 57, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [-0.0050  0.0087  0.0000]
    [ 0.0050 -0.0087  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 81, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 89, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000012 (xxy) 0.00000012 (xxy) 0.00000012 (xyx)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000025 (zy) -0.00000025 (zy) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 10:21:46]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


----------------------------
------ calculate LTC -------
----------------------------

        _                      _____
  _ __ | |__   ___  _ __   ___|___ / _ __  _   _
 | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | |
 | |_) | | | | (_) | | | | (_) |__) | |_) | |_| |
 | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, |
 |_|                                |_|    |___/ 
                                      3.23.0

-------------------------[time 2026-01-08 10:21:47]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [2 2 1]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   10.921106497213408    0.000000000000000    0.000000000000000
  b   -5.460553248606704    9.457955664022098    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
    1 Cs  0.20929583065401  0.79070416934600  0.45863866870204 132.905
    2 Cs  0.79070416934600  0.58140833869199  0.95863866870204 132.905
    3 Cs  0.41859166130801  0.20929583065401  0.95863866870204 132.905
    4 Cs  0.58140833869199  0.79070416934600  0.45863866870204 132.905
    5 Cs  0.20929583065401  0.41859166130801  0.45863866870204 132.905
    6 Cs  0.79070416934600  0.20929583065401  0.54136133129796 132.905
    7 Cs  0.79070416934600  0.58140833869199  0.54136133129796 132.905
    8 Cs  0.79070416934600  0.20929583065401  0.95863866870204 132.905
    9 Cs  0.41859166130801  0.20929583065401  0.54136133129796 132.905
   10 Cs  0.58140833869199  0.79070416934600  0.04136133129796 132.905
   11 Cs  0.20929583065401  0.79070416934600  0.04136133129796 132.905
   12 Cs  0.20929583065401  0.41859166130801  0.04136133129796 132.905
   13 Ga  0.33333333333333  0.66666666666667  0.75000000000000  69.723
   14 Ga  0.66666666666667  0.33333333333333  0.25000000000000  69.723
   15 Sb  0.47341102757325  0.94682205514651  0.75000000000000 121.760
   16 Sb  0.94682205514651  0.47341102757325  0.25000000000000 121.760
   17 Sb  0.52658897242674  0.47341102757325  0.25000000000000 121.760
   18 Sb  0.47341102757325  0.52658897242674  0.75000000000000 121.760
   19 Sb  0.05317794485349  0.52658897242674  0.75000000000000 121.760
   20 Sb  0.52658897242674  0.05317794485349  0.25000000000000 121.760
   21 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760
   22 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760
   23 K   0.11904674310150  0.88095325689850  0.75000000000000  39.098
   24 K   0.88095325689850  0.76190651379700  0.25000000000000  39.098
   25 K   0.23809348620300  0.11904674310150  0.25000000000000  39.098
   26 K   0.76190651379700  0.88095325689850  0.75000000000000  39.098
   27 K   0.11904674310150  0.23809348620300  0.75000000000000  39.098
   28 K   0.88095325689850  0.11904674310150  0.25000000000000  39.098
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   21.842212994426816    0.000000000000000    0.000000000000000
  b  -10.921106497213408   18.915911328044196    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.404273330000001
Atomic positions (fractional):
    1 Cs  0.10464791532700  0.39535208467300  0.45863866870204 132.905 > 1
    2 Cs  0.60464791532700  0.39535208467300  0.45863866870204 132.905 > 1
    3 Cs  0.10464791532700  0.89535208467300  0.45863866870204 132.905 > 1
    4 Cs  0.60464791532700  0.89535208467300  0.45863866870204 132.905 > 1
    5 Cs  0.39535208467300  0.29070416934600  0.95863866870204 132.905 > 5
    6 Cs  0.89535208467300  0.29070416934600  0.95863866870204 132.905 > 5
    7 Cs  0.39535208467300  0.79070416934600  0.95863866870204 132.905 > 5
    8 Cs  0.89535208467300  0.79070416934600  0.95863866870204 132.905 > 5
    9 Cs  0.20929583065400  0.10464791532700  0.95863866870204 132.905 > 9
   10 Cs  0.70929583065400  0.10464791532700  0.95863866870204 132.905 > 9
   11 Cs  0.20929583065400  0.60464791532700  0.95863866870204 132.905 > 9
   12 Cs  0.70929583065400  0.60464791532700  0.95863866870204 132.905 > 9
   13 Cs  0.29070416934600  0.39535208467300  0.45863866870204 132.905 > 13
   14 Cs  0.79070416934600  0.39535208467300  0.45863866870204 132.905 > 13
   15 Cs  0.29070416934600  0.89535208467300  0.45863866870204 132.905 > 13
   16 Cs  0.79070416934600  0.89535208467300  0.45863866870204 132.905 > 13
   17 Cs  0.10464791532700  0.20929583065400  0.45863866870204 132.905 > 17
   18 Cs  0.60464791532700  0.20929583065400  0.45863866870204 132.905 > 17
   19 Cs  0.10464791532700  0.70929583065400  0.45863866870204 132.905 > 17
   20 Cs  0.60464791532700  0.70929583065400  0.45863866870204 132.905 > 17
   21 Cs  0.39535208467300  0.10464791532700  0.54136133129796 132.905 > 21
   22 Cs  0.89535208467300  0.10464791532700  0.54136133129796 132.905 > 21
   23 Cs  0.39535208467300  0.60464791532700  0.54136133129796 132.905 > 21
   24 Cs  0.89535208467300  0.60464791532700  0.54136133129796 132.905 > 21
   25 Cs  0.39535208467300  0.29070416934600  0.54136133129796 132.905 > 25
   26 Cs  0.89535208467300  0.29070416934600  0.54136133129796 132.905 > 25
   27 Cs  0.39535208467300  0.79070416934600  0.54136133129796 132.905 > 25
   28 Cs  0.89535208467300  0.79070416934600  0.54136133129796 132.905 > 25
   29 Cs  0.39535208467300  0.10464791532700  0.95863866870204 132.905 > 29
   30 Cs  0.89535208467300  0.10464791532700  0.95863866870204 132.905 > 29
   31 Cs  0.39535208467300  0.60464791532700  0.95863866870204 132.905 > 29
   32 Cs  0.89535208467300  0.60464791532700  0.95863866870204 132.905 > 29
   33 Cs  0.20929583065400  0.10464791532700  0.54136133129796 132.905 > 33
   34 Cs  0.70929583065400  0.10464791532700  0.54136133129796 132.905 > 33
   35 Cs  0.20929583065400  0.60464791532700  0.54136133129796 132.905 > 33
   36 Cs  0.70929583065400  0.60464791532700  0.54136133129796 132.905 > 33
   37 Cs  0.29070416934600  0.39535208467300  0.04136133129796 132.905 > 37
   38 Cs  0.79070416934600  0.39535208467300  0.04136133129796 132.905 > 37
   39 Cs  0.29070416934600  0.89535208467300  0.04136133129796 132.905 > 37
   40 Cs  0.79070416934600  0.89535208467300  0.04136133129796 132.905 > 37
   41 Cs  0.10464791532700  0.39535208467300  0.04136133129796 132.905 > 41
   42 Cs  0.60464791532700  0.39535208467300  0.04136133129796 132.905 > 41
   43 Cs  0.10464791532700  0.89535208467300  0.04136133129796 132.905 > 41
   44 Cs  0.60464791532700  0.89535208467300  0.04136133129796 132.905 > 41
   45 Cs  0.10464791532700  0.20929583065400  0.04136133129796 132.905 > 45
   46 Cs  0.60464791532700  0.20929583065400  0.04136133129796 132.905 > 45
   47 Cs  0.10464791532700  0.70929583065400  0.04136133129796 132.905 > 45
   48 Cs  0.60464791532700  0.70929583065400  0.04136133129796 132.905 > 45
   49 Ga  0.16666666666667  0.33333333333333  0.75000000000000  69.723 > 49
   50 Ga  0.66666666666667  0.33333333333333  0.75000000000000  69.723 > 49
   51 Ga  0.16666666666667  0.83333333333333  0.75000000000000  69.723 > 49
   52 Ga  0.66666666666667  0.83333333333333  0.75000000000000  69.723 > 49
   53 Ga  0.33333333333333  0.16666666666667  0.25000000000000  69.723 > 53
   54 Ga  0.83333333333333  0.16666666666667  0.25000000000000  69.723 > 53
   55 Ga  0.33333333333333  0.66666666666667  0.25000000000000  69.723 > 53
   56 Ga  0.83333333333333  0.66666666666667  0.25000000000000  69.723 > 53
   57 Sb  0.23670551378663  0.47341102757325  0.75000000000000 121.760 > 57
   58 Sb  0.73670551378663  0.47341102757325  0.75000000000000 121.760 > 57
   59 Sb  0.23670551378663  0.97341102757325  0.75000000000000 121.760 > 57
   60 Sb  0.73670551378663  0.97341102757325  0.75000000000000 121.760 > 57
   61 Sb  0.47341102757325  0.23670551378663  0.25000000000000 121.760 > 61
   62 Sb  0.97341102757325  0.23670551378663  0.25000000000000 121.760 > 61
   63 Sb  0.47341102757325  0.73670551378663  0.25000000000000 121.760 > 61
   64 Sb  0.97341102757325  0.73670551378663  0.25000000000000 121.760 > 61
   65 Sb  0.26329448621337  0.23670551378663  0.25000000000000 121.760 > 65
   66 Sb  0.76329448621337  0.23670551378663  0.25000000000000 121.760 > 65
   67 Sb  0.26329448621337  0.73670551378663  0.25000000000000 121.760 > 65
   68 Sb  0.76329448621337  0.73670551378663  0.25000000000000 121.760 > 65
   69 Sb  0.23670551378663  0.26329448621337  0.75000000000000 121.760 > 69
   70 Sb  0.73670551378663  0.26329448621337  0.75000000000000 121.760 > 69
   71 Sb  0.23670551378663  0.76329448621337  0.75000000000000 121.760 > 69
   72 Sb  0.73670551378663  0.76329448621337  0.75000000000000 121.760 > 69
   73 Sb  0.02658897242675  0.26329448621337  0.75000000000000 121.760 > 73
   74 Sb  0.52658897242675  0.26329448621337  0.75000000000000 121.760 > 73
   75 Sb  0.02658897242675  0.76329448621337  0.75000000000000 121.760 > 73
   76 Sb  0.52658897242675  0.76329448621337  0.75000000000000 121.760 > 73
   77 Sb  0.26329448621337  0.02658897242675  0.25000000000000 121.760 > 77
   78 Sb  0.76329448621337  0.02658897242675  0.25000000000000 121.760 > 77
   79 Sb  0.26329448621337  0.52658897242675  0.25000000000000 121.760 > 77
   80 Sb  0.76329448621337  0.52658897242675  0.25000000000000 121.760 > 77
   81 Sb  0.00000000000000  0.00000000000000  0.00000000000000 121.760 > 81
   82 Sb  0.50000000000000  0.00000000000000  0.00000000000000 121.760 > 81
   83 Sb  0.00000000000000  0.50000000000000  0.00000000000000 121.760 > 81
   84 Sb  0.50000000000000  0.50000000000000  0.00000000000000 121.760 > 81
   85 Sb  0.00000000000000  0.00000000000000  0.50000000000000 121.760 > 85
   86 Sb  0.50000000000000  0.00000000000000  0.50000000000000 121.760 > 85
   87 Sb  0.00000000000000  0.50000000000000  0.50000000000000 121.760 > 85
   88 Sb  0.50000000000000  0.50000000000000  0.50000000000000 121.760 > 85
   89 K   0.05952337155075  0.44047662844925  0.75000000000000  39.098 > 89
   90 K   0.55952337155075  0.44047662844925  0.75000000000000  39.098 > 89
   91 K   0.05952337155075  0.94047662844925  0.75000000000000  39.098 > 89
   92 K   0.55952337155075  0.94047662844925  0.75000000000000  39.098 > 89
   93 K   0.44047662844925  0.38095325689850  0.25000000000000  39.098 > 93
   94 K   0.94047662844925  0.38095325689850  0.25000000000000  39.098 > 93
   95 K   0.44047662844925  0.88095325689850  0.25000000000000  39.098 > 93
   96 K   0.94047662844925  0.88095325689850  0.25000000000000  39.098 > 93
   97 K   0.11904674310150  0.05952337155075  0.25000000000000  39.098 > 97
   98 K   0.61904674310150  0.05952337155075  0.25000000000000  39.098 > 97
   99 K   0.11904674310150  0.55952337155075  0.25000000000000  39.098 > 97
  100 K   0.61904674310150  0.55952337155075  0.25000000000000  39.098 > 97
  101 K   0.38095325689850  0.44047662844925  0.75000000000000  39.098 > 101
  102 K   0.88095325689850  0.44047662844925  0.75000000000000  39.098 > 101
  103 K   0.38095325689850  0.94047662844925  0.75000000000000  39.098 > 101
  104 K   0.88095325689850  0.94047662844925  0.75000000000000  39.098 > 101
  105 K   0.05952337155075  0.11904674310150  0.75000000000000  39.098 > 105
  106 K   0.55952337155075  0.11904674310150  0.75000000000000  39.098 > 105
  107 K   0.05952337155075  0.61904674310150  0.75000000000000  39.098 > 105
  108 K   0.55952337155075  0.61904674310150  0.75000000000000  39.098 > 105
  109 K   0.44047662844925  0.05952337155075  0.25000000000000  39.098 > 109
  110 K   0.94047662844925  0.05952337155075  0.25000000000000  39.098 > 109
  111 K   0.44047662844925  0.55952337155075  0.25000000000000  39.098 > 109
  112 K   0.94047662844925  0.55952337155075  0.25000000000000  39.098 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            9.9027445    0.0000000    0.0000000
           -0.0000000    9.9027445    0.0000000
            0.0000000    0.0000000    6.0926774
-------------------------- Born effective charges --------------------------
    1 Cs    1.1405182    0.0439356    0.2763771
            0.0439356    1.1912507   -0.1595664
            0.0811556   -0.0468552    1.1503096
    2 Cs    1.2166170    0.0000000    0.0000000
           -0.0000000    1.1151519   -0.3191328
           -0.0000000   -0.0937104    1.1503096
    3 Cs    1.1405182   -0.0439356    0.2763771
           -0.0439356    1.1912507    0.1595664
            0.0811556    0.0468552    1.1503096
    4 Cs    1.1405182   -0.0439356   -0.2763771
           -0.0439356    1.1912507   -0.1595664
           -0.0811556   -0.0468552    1.1503096
    5 Cs    1.2166170    0.0000000   -0.0000000
           -0.0000000    1.1151519    0.3191328
           -0.0000000    0.0937104    1.1503096
    6 Cs    1.1405182    0.0439356    0.2763771
            0.0439356    1.1912507   -0.1595664
            0.0811556   -0.0468552    1.1503096
    7 Cs    1.2166170    0.0000000   -0.0000000
           -0.0000000    1.1151519    0.3191328
           -0.0000000    0.0937104    1.1503096
    8 Cs    1.1405182    0.0439356   -0.2763771
            0.0439356    1.1912507    0.1595664
           -0.0811556    0.0468552    1.1503096
    9 Cs    1.1405182   -0.0439356   -0.2763771
           -0.0439356    1.1912507   -0.1595664
           -0.0811556   -0.0468552    1.1503096
   10 Cs    1.1405182   -0.0439356    0.2763771
           -0.0439356    1.1912507    0.1595664
            0.0811556    0.0468552    1.1503096
   11 Cs    1.1405182    0.0439356   -0.2763771
            0.0439356    1.1912507    0.1595664
           -0.0811556    0.0468552    1.1503096
   12 Cs    1.2166170    0.0000000    0.0000000
           -0.0000000    1.1151519   -0.3191328
           -0.0000000   -0.0937104    1.1503096
   13 Ga    2.7985137    0.0000000    0.0000000
           -0.0000000    2.7985137    0.0000000
           -0.0000000   -0.0000000   -1.0476192
   14 Ga    2.7985137    0.0000000    0.0000000
           -0.0000000    2.7985137    0.0000000
           -0.0000000   -0.0000000   -1.0476192
   15 Sb   -2.2548971    0.0000000    0.0000000
            0.0000000   -4.1494014    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   16 Sb   -3.6757753   -0.8203444    0.0000000
           -0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   17 Sb   -3.6757753    0.8203444    0.0000000
            0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   18 Sb   -3.6757753    0.8203444    0.0000000
            0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   19 Sb   -3.6757753   -0.8203444    0.0000000
           -0.8203444   -2.7285232    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   20 Sb   -2.2548971    0.0000000    0.0000000
            0.0000000   -4.1494014    0.0000000
           -0.0000000   -0.0000000   -2.0398186
   21 Sb   -3.0985017   -0.0000000    0.0000000
            0.0000000   -3.0985017    0.0000000
           -0.0000000   -0.0000000   -2.5818845
   22 Sb   -3.0985017   -0.0000000    0.0000000
            0.0000000   -3.0985017    0.0000000
           -0.0000000   -0.0000000   -2.5818845
   23 K     1.1474204   -0.3066493    0.0000000
           -0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   24 K     0.6162882   -0.0000000    0.0000000
           -0.0000000    1.3244644    0.0000000
           -0.0000000   -0.0000000    0.9490339
   25 K     1.1474204    0.3066493    0.0000000
            0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   26 K     1.1474204    0.3066493    0.0000000
            0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
   27 K     0.6162882   -0.0000000    0.0000000
           -0.0000000    1.3244644    0.0000000
           -0.0000000   -0.0000000    0.9490339
   28 K     1.1474204   -0.3066493    0.0000000
           -0.3066493    0.7933323    0.0000000
           -0.0000000   -0.0000000    0.9490339
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000012 (xxy) 0.00000012 (xxy) 0.00000012 (xyx)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 5 5 4 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.39, Number of G-points: 315, Lambda: 0.12
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/15) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 15
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   0.753   (   0.000    0.000    0.000)    0.000
   0.786   (   0.000    0.000    0.000)    0.000
   0.863   (   0.000    0.000    0.000)    0.000
   0.863   (   0.000    0.000    0.000)    0.000
   1.028   (   0.000    0.000    0.000)    0.000
   1.028   (   0.000    0.000    0.000)    0.000
   1.059   (   0.000    0.000    0.000)    0.000
   1.059   (   0.000    0.000    0.000)    0.000
   1.072   (   0.000    0.000    0.000)    0.000
   1.072   (   0.000    0.000    0.000)    0.000
   1.164   (   0.000    0.000    0.000)    0.000
   1.206   (   0.000    0.000    0.000)    0.000
   1.206   (   0.000    0.000    0.000)    0.000
   1.249   (   0.000    0.000    0.000)    0.000
   1.249   (   0.000    0.000    0.000)    0.000
   1.283   (   0.000    0.000    0.000)    0.000
   1.283   (   0.000    0.000    0.000)    0.000
   1.299   (   0.000    0.000    0.000)    0.000
   1.315   (   0.000    0.000    0.000)    0.000
   1.379   (   0.000    0.000    0.000)    0.000
   1.379   (   0.000    0.000    0.000)    0.000
   1.414   (   0.000    0.000    0.000)    0.000
   1.414   (   0.000    0.000    0.000)    0.000
   1.428   (   0.000    0.000    0.000)    0.000
   1.428   (   0.000    0.000    0.000)    0.000
   1.443   (   0.000    0.000    0.000)    0.000
   1.443   (   0.000    0.000    0.000)    0.000
   1.479   (   0.000    0.000    0.000)    0.000
   1.509   (   0.000    0.000    0.000)    0.000
   1.509   (   0.000    0.000    0.000)    0.000
   1.519   (   0.000    0.000    0.000)    0.000
   1.530   (   0.000    0.000    0.000)    0.000
   1.533   (   0.000    0.000    0.000)    0.000
   1.556   (   0.000    0.000    0.000)    0.000
   1.556   (   0.000    0.000    0.000)    0.000
   1.724   (   0.000    0.000    0.000)    0.000
   1.747   (   0.000    0.000    0.000)    0.000
   1.747   (   0.000    0.000    0.000)    0.000
   1.751   (   0.000    0.000    0.000)    0.000
   1.793   (   0.000    0.000    0.000)    0.000
   1.796   (   0.000    0.000    0.000)    0.000
   1.796   (   0.000    0.000    0.000)    0.000
   1.802   (   0.000    0.000    0.000)    0.000
   1.860   (   0.000    0.000    0.000)    0.000
   1.913   (   0.000    0.000    0.000)    0.000
   1.913   (   0.000    0.000    0.000)    0.000
   1.928   (   0.000    0.000    0.000)    0.000
   1.928   (   0.000    0.000    0.000)    0.000
   1.956   (   0.000    0.000    0.000)    0.000
   2.057   (   0.000    0.000    0.000)    0.000
   2.089   (   0.000    0.000    0.000)    0.000
   2.089   (   0.000    0.000    0.000)    0.000
   2.440   (   0.000    0.000    0.000)    0.000
   2.440   (   0.000    0.000    0.000)    0.000
   2.473   (   0.000    0.000    0.000)    0.000
   3.005   (   0.000    0.000    0.000)    0.000
   3.031   (   0.000    0.000    0.000)    0.000
   3.051   (   0.000    0.000    0.000)    0.000
   3.095   (   0.000    0.000    0.000)    0.000
   3.185   (   0.000    0.000    0.000)    0.000
   3.185   (   0.000    0.000    0.000)    0.000
   3.234   (   0.000    0.000    0.000)    0.000
   3.234   (   0.000    0.000    0.000)    0.000
   3.319   (   0.000    0.000    0.000)    0.000
   3.319   (   0.000    0.000    0.000)    0.000
   3.331   (   0.000    0.000    0.000)    0.000
   3.331   (   0.000    0.000    0.000)    0.000
   3.343   (   0.000    0.000    0.000)    0.000
   3.343   (   0.000    0.000    0.000)    0.000
   3.397   (   0.000    0.000    0.000)    0.000
   3.397   (   0.000    0.000    0.000)    0.000
   3.630   (   0.000    0.000    0.000)    0.000
   3.651   (   0.000    0.000    0.000)    0.000
   3.683   (   0.000    0.000    0.000)    0.000
   3.739   (   0.000    0.000    0.000)    0.000
   3.842   (   0.000    0.000    0.000)    0.000
   3.911   (   0.000    0.000    0.000)    0.000
   5.662   (   0.000    0.000    0.000)    0.000
   5.662   (   0.000    0.000    0.000)    0.000
   5.760   (   0.000    0.000    0.000)    0.000
   5.760   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/15) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 27
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.287   (  11.398    6.581    0.000)   13.162
   0.307   (  12.279    7.089    0.000)   14.178
   0.528   (  20.923   12.080    0.000)   24.160
   0.800   (   3.436    1.984    0.000)    3.968
   0.829   (   2.696    1.557    0.000)    3.113
   0.877   (   1.158    0.668    0.000)    1.337
   0.912   (   2.882    1.664    0.000)    3.328
   0.972   (  -3.773   -2.178    0.000)    4.357
   1.002   (  -1.424   -0.822    0.000)    1.644
   1.065   (   0.706    0.408    0.000)    0.815
   1.071   (   0.352    0.203    0.000)    0.406
   1.091   (   1.035    0.597    0.000)    1.195
   1.117   (   4.172    2.409    0.000)    4.817
   1.118   (  -3.958   -2.285    0.000)    4.570
   1.205   (  -1.650   -0.953    0.000)    1.906
   1.206   (  -2.907   -1.678    0.000)    3.357
   1.248   (  -2.193   -1.266    0.000)    2.533
   1.250   (   1.719    0.993    0.000)    1.985
   1.258   (   2.807    1.620    0.000)    3.241
   1.264   (  -1.431   -0.826    0.000)    1.652
   1.269   (  -1.944   -1.122    0.000)    2.245
   1.319   (   2.128    1.229    0.000)    2.457
   1.383   (  -1.676   -0.967    0.000)    1.935
   1.405   (  -1.238   -0.715    0.000)    1.429
   1.418   (  -3.137   -1.811    0.000)    3.623
   1.422   (  -0.294   -0.170    0.000)    0.340
   1.424   (  -0.673   -0.388    0.000)    0.777
   1.455   (   0.633    0.366    0.000)    0.731
   1.463   (   3.564    2.058    0.000)    4.115
   1.467   (  -1.743   -1.007    0.000)    2.013
   1.472   (   4.443    2.565    0.000)    5.130
   1.496   (  -1.936   -1.118    0.000)    2.235
   1.517   (   2.503    1.445    0.000)    2.890
   1.553   (   0.304    0.176    0.000)    0.351
   1.555   (   0.461    0.266    0.000)    0.533
   1.569   (   3.480    2.009    0.000)    4.019
   1.585   (   4.525    2.612    0.000)    5.225
   1.684   (  -5.080   -2.933    0.000)    5.866
   1.711   (  -2.497   -1.442    0.000)    2.884
   1.735   (   0.297    0.171    0.000)    0.343
   1.745   (  -0.182   -0.105    0.000)    0.210
   1.752   (   0.680    0.393    0.000)    0.786
   1.781   (   1.850    1.068    0.000)    2.136
   1.801   (   2.923    1.688    0.000)    3.376
   1.821   (   1.958    1.130    0.000)    2.260
   1.826   (   1.050    0.606    0.000)    1.213
   1.845   (  -1.103   -0.637    0.000)    1.273
   1.868   (  -2.067   -1.194    0.000)    2.387
   1.891   (  -2.278   -1.315    0.000)    2.631
   1.919   (   0.502    0.290    0.000)    0.580
   1.937   (  -1.516   -0.875    0.000)    1.751
   2.015   (   4.357    2.516    0.000)    5.032
   2.054   (  -0.196   -0.113    0.000)    0.226
   2.108   (   1.220    0.705    0.000)    1.409
   2.131   (   1.612    0.931    0.000)    1.861
   2.415   (  -2.306   -1.331    0.000)    2.663
   2.458   (  -3.651   -2.108    0.000)    4.216
   2.479   (   1.380    0.797    0.000)    1.594
   2.917   (  -8.891   -5.133    0.000)   10.266
   2.970   (  -7.619   -4.399    0.000)    8.798
   2.999   (  -0.426   -0.246    0.000)    0.491
   3.020   (  -0.788   -0.455    0.000)    0.909
   3.093   (  -6.699   -3.868    0.000)    7.736
   3.160   (  -0.202   -0.117    0.000)    0.233
   3.200   (  -1.203   -0.695    0.000)    1.390
   3.223   (  -5.186   -2.994    0.000)    5.988
   3.320   (   0.052    0.030    0.000)    0.061
   3.321   (   0.160    0.092    0.000)    0.184
   3.367   (   2.105    1.215    0.000)    2.431
   3.396   (  -0.134   -0.077    0.000)    0.154
   3.397   (  -0.056   -0.032    0.000)    0.065
   3.432   (   5.516    3.185    0.000)    6.370
   3.456   (  -1.153   -0.666    0.000)    1.331
   3.463   (   5.272    3.044    0.000)    6.088
   3.631   (   0.118    0.068    0.000)    0.136
   3.672   (   1.016    0.587    0.000)    1.173
   3.681   (  -0.153   -0.088    0.000)    0.176
   3.733   (  -0.472   -0.273    0.000)    0.546
   3.865   (   1.831    1.057    0.000)    2.114
   3.970   (   3.888    2.245    0.000)    4.489
   5.653   (  -0.579   -0.334    0.000)    0.669
   5.744   (  -1.162   -0.671    0.000)    1.342
   5.798   (   2.588    1.494    0.000)    2.988
   5.870   (   0.673    0.388    0.000)    0.777
======================= Grid point 2 (3/15) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 27
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.541   (   8.655    4.997    0.000)    9.994
   0.591   (  10.658    6.154    0.000)   12.307
   0.813   (  -6.434   -3.715    0.000)    7.429
   0.814   (  -5.783   -3.339    0.000)    6.677
   0.887   (   2.754    1.590    0.000)    3.180
   0.910   (   1.059    0.611    0.000)    1.223
   0.975   (  -0.660   -0.381    0.000)    0.762
   0.994   (  15.580    8.995    0.000)   17.990
   1.004   (   1.331    0.769    0.000)    1.537
   1.086   (   2.296    1.325    0.000)    2.651
   1.090   (   1.478    0.853    0.000)    1.707
   1.115   (   3.742    2.161    0.000)    4.321
   1.121   (   1.752    1.011    0.000)    2.023
   1.142   (  -1.765   -1.019    0.000)    2.038
   1.161   (  -0.229   -0.132    0.000)    0.264
   1.201   (  -1.867   -1.078    0.000)    2.156
   1.216   (  -2.714   -1.567    0.000)    3.134
   1.223   (   3.074    1.775    0.000)    3.550
   1.229   (  -3.774   -2.179    0.000)    4.358
   1.237   (   0.365    0.211    0.000)    0.421
   1.295   (  -3.332   -1.924    0.000)    3.848
   1.316   (   2.261    1.306    0.000)    2.611
   1.324   (  -5.657   -3.266    0.000)    6.532
   1.350   (  -0.391   -0.226    0.000)    0.452
   1.391   (  -1.061   -0.612    0.000)    1.225
   1.401   (  -5.540   -3.199    0.000)    6.397
   1.405   (   0.025    0.014    0.000)    0.029
   1.407   (   2.356    1.360    0.000)    2.720
   1.449   (  -0.406   -0.234    0.000)    0.469
   1.483   (  -6.515   -3.762    0.000)    7.523
   1.493   (   1.457    0.841    0.000)    1.682
   1.557   (   0.997    0.576    0.000)    1.151
   1.561   (   0.239    0.138    0.000)    0.276
   1.562   (   2.372    1.369    0.000)    2.739
   1.595   (   1.629    0.941    0.000)    1.881
   1.602   (  -1.981   -1.144    0.000)    2.288
   1.631   (   1.871    1.080    0.000)    2.160
   1.659   (   2.340    1.351    0.000)    2.702
   1.699   (   3.681    2.125    0.000)    4.250
   1.716   (  -0.870   -0.502    0.000)    1.005
   1.779   (   1.268    0.732    0.000)    1.464
   1.785   (   1.550    0.895    0.000)    1.789
   1.811   (  -0.808   -0.466    0.000)    0.933
   1.824   (   0.657    0.379    0.000)    0.758
   1.827   (  -0.916   -0.529    0.000)    1.057
   1.842   (  -0.828   -0.478    0.000)    0.957
   1.853   (   1.118    0.645    0.000)    1.291
   1.877   (   0.754    0.436    0.000)    0.871
   1.889   (   2.391    1.380    0.000)    2.760
   1.894   (  -1.469   -0.848    0.000)    1.696
   1.932   (   0.398    0.230    0.000)    0.460
   2.079   (   1.081    0.624    0.000)    1.248
   2.085   (   1.376    0.795    0.000)    1.589
   2.152   (   2.796    1.614    0.000)    3.228
   2.154   (   2.957    1.707    0.000)    3.415
   2.314   (  -5.367   -3.099    0.000)    6.197
   2.325   (  -6.606   -3.814    0.000)    7.628
   2.565   (   3.957    2.285    0.000)    4.569
   2.721   (  -5.401   -3.118    0.000)    6.237
   2.822   (  -3.340   -1.929    0.000)    3.857
   2.947   (  -4.358   -2.516    0.000)    5.032
   2.991   (  -0.146   -0.084    0.000)    0.168
   3.000   (  -0.604   -0.349    0.000)    0.697
   3.101   (  -3.353   -1.936    0.000)    3.871
   3.162   (   0.282    0.163    0.000)    0.326
   3.180   (  -0.567   -0.327    0.000)    0.655
   3.321   (   0.031    0.018    0.000)    0.036
   3.325   (   0.099    0.057    0.000)    0.115
   3.392   (  -0.084   -0.048    0.000)    0.097
   3.395   (   0.447    0.258    0.000)    0.516
   3.395   (  -0.033   -0.019    0.000)    0.038
   3.435   (  -0.482   -0.278    0.000)    0.556
   3.529   (   2.128    1.229    0.000)    2.457
   3.554   (   1.933    1.116    0.000)    2.232
   3.635   (   0.146    0.084    0.000)    0.168
   3.678   (  -0.092   -0.053    0.000)    0.107
   3.690   (   0.474    0.273    0.000)    0.547
   3.721   (  -0.392   -0.226    0.000)    0.453
   3.909   (   1.204    0.695    0.000)    1.390
   4.050   (   1.962    1.133    0.000)    2.265
   5.642   (  -0.230   -0.133    0.000)    0.266
   5.715   (  -0.865   -0.500    0.000)    0.999
   5.853   (   1.344    0.776    0.000)    1.552
   5.887   (   0.514    0.297    0.000)    0.593
======================= Grid point 6 (4/15) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 27
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.481   (   5.804   10.053    0.000)   11.608
   0.538   (   7.120   12.332    0.000)   14.240
   0.867   (   9.676   16.759    0.000)   19.352
   0.874   (   2.409    4.172    0.000)    4.818
   0.881   (   1.187    2.056    0.000)    2.374
   0.890   (  -0.956   -1.657    0.000)    1.913
   0.898   (  -1.267   -2.195    0.000)    2.535
   0.915   (  -1.704   -2.952    0.000)    3.409
   1.026   (   2.867    4.966    0.000)    5.734
   1.062   (  -2.043   -3.539    0.000)    4.086
   1.097   (   0.957    1.658    0.000)    1.915
   1.099   (   1.338    2.317    0.000)    2.675
   1.100   (   0.026    0.045    0.000)    0.052
   1.125   (  -1.508   -2.613    0.000)    3.017
   1.166   (   1.301    2.254    0.000)    2.602
   1.219   (  -0.412   -0.714    0.000)    0.825
   1.223   (  -1.524   -2.640    0.000)    3.048
   1.223   (  -0.470   -0.813    0.000)    0.939
   1.234   (   0.410    0.710    0.000)    0.820
   1.238   (  -1.001   -1.734    0.000)    2.002
   1.291   (   0.418    0.723    0.000)    0.835
   1.332   (   0.129    0.224    0.000)    0.258
   1.336   (  -2.019   -3.497    0.000)    4.038
   1.379   (  -2.031   -3.518    0.000)    4.062
   1.409   (  -0.921   -1.596    0.000)    1.843
   1.414   (  -0.072   -0.124    0.000)    0.144
   1.417   (   0.377    0.652    0.000)    0.753
   1.432   (  -0.138   -0.239    0.000)    0.276
   1.455   (  -0.646   -1.119    0.000)    1.292
   1.494   (  -1.020   -1.766    0.000)    2.040
   1.504   (   0.056    0.097    0.000)    0.112
   1.513   (   1.406    2.434    0.000)    2.811
   1.558   (   2.851    4.937    0.000)    5.701
   1.563   (   1.230    2.130    0.000)    2.460
   1.564   (  -0.844   -1.463    0.000)    1.689
   1.595   (  -2.125   -3.681    0.000)    4.250
   1.648   (   1.937    3.356    0.000)    3.875
   1.656   (   0.916    1.586    0.000)    1.831
   1.701   (   0.502    0.869    0.000)    1.004
   1.717   (  -2.527   -4.376    0.000)    5.053
   1.745   (   0.206    0.357    0.000)    0.412
   1.757   (   0.277    0.480    0.000)    0.554
   1.795   (   0.162    0.281    0.000)    0.324
   1.795   (   0.796    1.378    0.000)    1.591
   1.822   (  -0.722   -1.251    0.000)    1.445
   1.826   (   0.248    0.429    0.000)    0.496
   1.856   (   1.370    2.372    0.000)    2.739
   1.895   (   1.549    2.684    0.000)    3.099
   1.904   (  -1.100   -1.905    0.000)    2.200
   1.919   (   1.765    3.057    0.000)    3.529
   1.955   (   1.442    2.497    0.000)    2.884
   2.049   (   0.432    0.749    0.000)    0.865
   2.062   (   0.744    1.288    0.000)    1.487
   2.098   (  -1.197   -2.073    0.000)    2.394
   2.182   (   2.420    4.192    0.000)    4.841
   2.335   (  -4.816   -8.341    0.000)    9.631
   2.388   (  -1.213   -2.100    0.000)    2.425
   2.526   (   1.555    2.694    0.000)    3.110
   2.781   (  -2.932   -5.078    0.000)    5.863
   2.879   (  -2.230   -3.862    0.000)    4.459
   2.992   (  -0.188   -0.325    0.000)    0.375
   3.005   (  -0.471   -0.816    0.000)    0.942
   3.023   (  -3.033   -5.253    0.000)    6.066
   3.051   (  -2.415   -4.183    0.000)    4.830
   3.084   (  -2.198   -3.808    0.000)    4.397
   3.197   (  -2.684   -4.649    0.000)    5.369
   3.321   (   0.055    0.095    0.000)    0.110
   3.324   (   0.061    0.106    0.000)    0.123
   3.393   (  -0.054   -0.093    0.000)    0.107
   3.395   (  -0.048   -0.084    0.000)    0.097
   3.435   (   1.473    2.552    0.000)    2.946
   3.474   (   0.899    1.557    0.000)    1.798
   3.500   (   1.735    3.005    0.000)    3.470
   3.524   (   1.738    3.010    0.000)    3.475
   3.634   (   0.122    0.210    0.000)    0.243
   3.678   (  -0.081   -0.140    0.000)    0.162
   3.687   (   0.448    0.776    0.000)    0.896
   3.723   (  -0.327   -0.567    0.000)    0.654
   3.898   (   0.914    1.584    0.000)    1.829
   4.034   (   1.679    2.908    0.000)    3.357
   5.659   (   0.467    0.810    0.000)    0.935
   5.733   (  -0.199   -0.344    0.000)    0.397
   5.833   (   0.754    1.305    0.000)    1.507
   5.869   (  -0.261   -0.453    0.000)    0.523
======================= Grid point 7 (5/15) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 27
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.638   (  -1.420    2.459    0.000)    2.840
   0.751   (  -2.458    4.258    0.000)    4.917
   0.783   (  -1.212    2.099    0.000)    2.423
   0.821   (  -3.751    6.498    0.000)    7.503
   0.887   (   3.018   -5.227    0.000)    6.036
   0.901   (   0.505   -0.874    0.000)    1.010
   0.951   (  -2.303    3.988    0.000)    4.605
   1.019   (   3.963   -6.865    0.000)    7.927
   1.073   (  -0.314    0.545    0.000)    0.629
   1.090   (   0.718   -1.244    0.000)    1.437
   1.092   (   0.303   -0.525    0.000)    0.606
   1.104   (   0.887   -1.536    0.000)    1.774
   1.136   (   0.214   -0.371    0.000)    0.429
   1.183   (  -0.138    0.240    0.000)    0.277
   1.189   (   0.032   -0.056    0.000)    0.064
   1.191   (  -2.669    4.623    0.000)    5.338
   1.195   (  -0.769    1.332    0.000)    1.538
   1.235   (   0.798   -1.382    0.000)    1.596
   1.238   (  -1.657    2.870    0.000)    3.313
   1.255   (  -0.627    1.086    0.000)    1.254
   1.274   (  -0.191    0.331    0.000)    0.382
   1.300   (  -0.691    1.197    0.000)    1.382
   1.305   (  -0.612    1.060    0.000)    1.224
   1.327   (   0.784   -1.357    0.000)    1.567
   1.370   (   0.451   -0.781    0.000)    0.902
   1.408   (  -0.960    1.662    0.000)    1.919
   1.418   (   0.209   -0.362    0.000)    0.418
   1.418   (  -0.684    1.185    0.000)    1.369
   1.453   (  -0.099    0.171    0.000)    0.198
   1.472   (   2.164   -3.748    0.000)    4.328
   1.479   (  -1.377    2.385    0.000)    2.754
   1.502   (   2.645   -4.581    0.000)    5.289
   1.547   (  -0.496    0.859    0.000)    0.992
   1.550   (  -1.466    2.540    0.000)    2.933
   1.595   (   0.560   -0.970    0.000)    1.120
   1.614   (   3.570   -6.183    0.000)    7.139
   1.663   (   0.385   -0.667    0.000)    0.771
   1.680   (  -2.529    4.381    0.000)    5.059
   1.698   (  -1.140    1.975    0.000)    2.281
   1.704   (   2.537   -4.395    0.000)    5.074
   1.760   (   1.109   -1.921    0.000)    2.218
   1.769   (   1.130   -1.958    0.000)    2.261
   1.800   (   0.087   -0.150    0.000)    0.174
   1.802   (   0.643   -1.114    0.000)    1.287
   1.809   (  -2.550    4.416    0.000)    5.099
   1.860   (   0.327   -0.566    0.000)    0.654
   1.877   (   0.324   -0.561    0.000)    0.648
   1.883   (  -0.631    1.093    0.000)    1.262
   1.920   (  -1.334    2.310    0.000)    2.667
   1.974   (  -2.669    4.623    0.000)    5.338
   1.989   (  -2.158    3.737    0.000)    4.315
   2.063   (   0.882   -1.528    0.000)    1.764
   2.072   (   1.144   -1.982    0.000)    2.289
   2.101   (   2.249   -3.895    0.000)    4.498
   2.183   (   2.042   -3.536    0.000)    4.083
   2.272   (  -2.165    3.750    0.000)    4.330
   2.331   (  -0.743    1.287    0.000)    1.487
   2.534   (   1.770   -3.065    0.000)    3.539
   2.773   (  -3.443    5.964    0.000)    6.886
   2.830   (  -1.351    2.340    0.000)    2.702
   2.919   (   0.465   -0.806    0.000)    0.931
   2.989   (   0.021   -0.036    0.000)    0.041
   2.991   (   4.039   -6.995    0.000)    8.077
   2.993   (   0.141   -0.244    0.000)    0.282
   3.080   (   2.103   -3.643    0.000)    4.207
   3.122   (   2.141   -3.708    0.000)    4.282
   3.322   (  -0.063    0.109    0.000)    0.126
   3.325   (   0.032   -0.055    0.000)    0.064
   3.392   (  -0.023    0.040    0.000)    0.046
   3.394   (   0.047   -0.081    0.000)    0.093
   3.459   (  -1.832    3.174    0.000)    3.665
   3.482   (  -1.258    2.178    0.000)    2.515
   3.542   (  -0.266    0.461    0.000)    0.533
   3.563   (   0.181   -0.314    0.000)    0.363
   3.637   (  -0.043    0.074    0.000)    0.086
   3.677   (   0.019   -0.033    0.000)    0.038
   3.699   (  -0.294    0.509    0.000)    0.588
   3.714   (   0.212   -0.367    0.000)    0.424
   3.917   (  -0.023    0.040    0.000)    0.046
   4.068   (  -0.241    0.418    0.000)    0.482
   5.673   (  -1.245    2.157    0.000)    2.490
   5.729   (  -0.730    1.264    0.000)    1.460
   5.845   (   0.552   -0.956    0.000)    1.103
   5.862   (   1.169   -2.025    0.000)    2.338
======================= Grid point 31 (6/15) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 10
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.277   (  -0.000   -0.000   12.440)   12.440
   0.277   (   0.000   -0.000   12.440)   12.440
   0.552   (   0.000    0.000   24.385)   24.385
   0.756   (  -0.000   -0.000   -8.227)    8.227
   0.756   (  -0.000    0.000   -8.227)    8.227
   0.758   (   0.000    0.000    0.413)    0.413
   0.781   (   0.000    0.000   -0.418)    0.418
   1.064   (   0.000   -0.000    0.190)    0.190
   1.064   (  -0.000    0.000    0.190)    0.190
   1.069   (   0.000   -0.000    2.995)    2.995
   1.069   (  -0.000   -0.000    2.995)    2.995
   1.101   (   0.000    0.000    2.554)    2.554
   1.101   (   0.000   -0.000    2.554)    2.554
   1.185   (   0.000    0.000    1.792)    1.792
   1.224   (  -0.000    0.000   -0.072)    0.072
   1.224   (   0.000   -0.000   -0.072)    0.072
   1.243   (   0.000    0.000   -0.326)    0.326
   1.243   (   0.000    0.000   -0.326)    0.326
   1.269   (   0.000    0.000   -0.883)    0.883
   1.269   (   0.000    0.000   -0.883)    0.883
   1.291   (   0.000   -0.000   -2.073)    2.073
   1.323   (   0.000    0.000    1.968)    1.968
   1.379   (   0.000   -0.000   -2.282)    2.282
   1.379   (   0.000   -0.000   -2.282)    2.282
   1.380   (  -0.000   -0.000   -0.887)    0.887
   1.380   (   0.000   -0.000   -0.887)    0.887
   1.425   (   0.000    0.000    2.110)    2.110
   1.425   (   0.000   -0.000    2.110)    2.110
   1.472   (   0.000    0.000    1.006)    1.006
   1.472   (  -0.000    0.000    1.006)    1.006
   1.473   (   0.000    0.000  -16.437)   16.437
   1.478   (   0.000    0.000   -3.699)    3.699
   1.510   (  -0.000    0.000    0.264)    0.264
   1.510   (   0.000   -0.000    0.264)    0.264
   1.524   (  -0.000    0.000    0.322)    0.322
   1.531   (   0.000    0.000    0.011)    0.011
   1.559   (   0.000   -0.000    0.006)    0.006
   1.559   (   0.000    0.000    0.006)    0.006
   1.641   (   0.000    0.000    4.373)    4.373
   1.742   (  -0.000    0.000   -0.622)    0.622
   1.742   (   0.000   -0.000   -0.622)    0.622
   1.751   (   0.000    0.000   -2.662)    2.662
   1.751   (   0.000   -0.000   -2.662)    2.662
   1.754   (   0.000    0.000   -0.820)    0.820
   1.809   (   0.000   -0.000    2.581)    2.581
   1.835   (   0.000    0.000    1.408)    1.408
   1.918   (   0.000    0.000    0.051)    0.051
   1.918   (   0.000   -0.000    0.051)    0.051
   1.920   (   0.000    0.000   -0.520)    0.520
   1.920   (   0.000   -0.000   -0.520)    0.520
   2.011   (   0.000   -0.000    5.672)    5.672
   2.045   (   0.000    0.000   -0.745)    0.745
   2.162   (   0.000   -0.000    5.678)    5.678
   2.162   (   0.000    0.000    5.678)    5.678
   2.390   (   0.000   -0.000   -1.216)    1.216
   2.401   (   0.000    0.000   -3.460)    3.460
   2.401   (   0.000   -0.000   -3.460)    3.460
   2.416   (   0.000   -0.000   -4.111)    4.111
   3.014   (   0.000    0.000    0.713)    0.713
   3.049   (   0.000   -0.000   -0.571)    0.571
   3.058   (  -0.000   -0.000    0.567)    0.567
   3.088   (   0.000    0.000   -0.542)    0.542
   3.189   (   0.000    0.000    0.442)    0.442
   3.189   (   0.000    0.000    0.442)    0.442
   3.226   (   0.000    0.000   -0.764)    0.764
   3.226   (   0.000   -0.000   -0.764)    0.764
   3.329   (   0.000    0.000   -0.094)    0.094
   3.329   (   0.000    0.000   -0.094)    0.094
   3.332   (  -0.000   -0.000    0.485)    0.485
   3.332   (  -0.000    0.000    0.485)    0.485
   3.344   (   0.000    0.000    0.633)    0.633
   3.344   (   0.000    0.000    0.633)    0.633
   3.392   (   0.000   -0.000   -0.548)    0.548
   3.392   (   0.000    0.000   -0.548)    0.548
   3.662   (   0.000   -0.000    2.622)    2.622
   3.663   (   0.000    0.000    1.091)    1.091
   3.726   (  -0.000   -0.000   -1.041)    1.041
   3.810   (   0.000    0.000   -2.921)    2.921
   3.860   (   0.000   -0.000    1.426)    1.426
   3.912   (   0.000    0.000   -0.253)    0.253
   5.676   (   0.000    0.000    1.227)    1.227
   5.676   (   0.000    0.000    1.227)    1.227
   5.746   (   0.000   -0.000   -1.250)    1.250
   5.746   (   0.000   -0.000   -1.250)    1.250
======================= Grid point 32 (7/15) =======================
q-point: ( 0.20  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 30
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.403   (   8.471    4.890    7.568)   12.367
   0.442   (   9.305    5.372    8.972)   13.998
   0.714   (  11.994    6.925   15.635)   20.887
   0.774   (   1.853    1.070   -5.878)    6.256
   0.812   (   3.437    1.984   -0.790)    4.046
   0.822   (   4.858    2.805   -7.165)    9.100
   0.825   (   3.020    1.744   -0.626)    3.543
   1.004   (  -3.440   -1.986    2.960)    4.954
   1.018   (  -1.867   -1.078    1.473)    2.611
   1.088   (  -0.602   -0.348    0.575)    0.902
   1.104   (   0.935    0.540    2.338)    2.576
   1.111   (   1.268    0.732    2.220)    2.659
   1.146   (   3.186    1.839    1.338)    3.914
   1.150   (  -3.391   -1.958    2.118)    4.451
   1.192   (  -2.400   -1.386   -0.881)    2.908
   1.204   (  -2.335   -1.348   -0.255)    2.708
   1.220   (  -0.166   -0.096   -1.954)    1.963
   1.239   (  -1.695   -0.978   -0.904)    2.156
   1.244   (  -1.384   -0.799   -1.557)    2.231
   1.271   (   2.603    1.503    0.711)    3.089
   1.288   (  -2.165   -1.250    1.537)    2.935
   1.317   (   0.975    0.563   -0.497)    1.230
   1.369   (  -1.118   -0.646   -3.180)    3.432
   1.372   (  -0.998   -0.576   -0.453)    1.238
   1.386   (   0.352    0.203   -1.670)    1.719
   1.422   (   2.153    1.243   -1.543)    2.926
   1.423   (  -1.657   -0.957   -3.795)    4.250
   1.437   (  -2.961   -1.709    1.358)    3.679
   1.453   (  -0.312   -0.180    0.316)    0.479
   1.469   (  -0.309   -0.178   -4.403)    4.417
   1.482   (   2.892    1.670   -0.046)    3.339
   1.493   (   1.558    0.900    0.735)    1.944
   1.523   (   3.106    1.793   -3.628)    5.101
   1.547   (   0.767    0.443   -0.471)    1.003
   1.569   (   3.812    2.201    0.040)    4.402
   1.570   (   0.638    0.368    1.253)    1.453
   1.588   (   3.618    2.089    0.164)    4.181
   1.594   (   2.051    1.184    4.656)    5.224
   1.667   (   2.821    1.629   -2.198)    3.930
   1.715   (  -2.526   -1.458   -0.809)    3.027
   1.750   (  -3.438   -1.985    1.020)    4.099
   1.769   (   1.824    1.053    0.173)    2.114
   1.783   (  -1.024   -0.591    0.009)    1.182
   1.783   (   2.431    1.404   -1.554)    3.209
   1.793   (   3.034    1.752   -1.526)    3.822
   1.808   (   2.051    1.184   -0.994)    2.568
   1.876   (  -1.401   -0.809    1.277)    2.061
   1.891   (  -0.953   -0.550    0.724)    1.317
   1.934   (   0.198    0.114   -0.094)    0.247
   1.945   (  -3.577   -2.065    3.299)    5.287
   1.978   (   2.307    1.332   -0.295)    2.680
   2.044   (  -0.804   -0.464   -0.577)    1.093
   2.122   (  -6.055   -3.496    8.209)   10.783
   2.163   (  -0.001   -0.001    4.656)    4.656
   2.204   (  -0.474   -0.274    5.323)    5.351
   2.378   (  -2.061   -1.190   -3.250)    4.029
   2.410   (   1.152    0.665   -5.245)    5.411
   2.426   (  -2.223   -1.284   -3.231)    4.127
   2.925   (  -8.775   -5.066    0.679)   10.155
   2.962   (  -7.894   -4.557   -0.662)    9.139
   3.004   (  -0.645   -0.373    0.428)    0.859
   3.019   (  -1.188   -0.686   -0.115)    1.377
   3.097   (  -6.049   -3.493    0.480)    7.002
   3.161   (  -0.316   -0.182    0.226)    0.429
   3.175   (  -4.191   -2.420   -3.126)    5.761
   3.190   (  -1.034   -0.597   -0.796)    1.435
   3.337   (   0.018    0.011    1.341)    1.341
   3.337   (   0.193    0.111    1.184)    1.205
   3.371   (   1.929    1.114    0.370)    2.258
   3.389   (  -0.054   -0.031   -0.634)    0.637
   3.390   (   0.078    0.045   -0.582)    0.588
   3.417   (   1.502    0.867   -1.786)    2.489
   3.436   (   5.474    3.160    0.370)    6.332
   3.458   (   5.269    3.042   -0.414)    6.098
   3.665   (  -0.016   -0.009    2.711)    2.711
   3.679   (   0.643    0.371    0.703)    1.022
   3.723   (  -0.541   -0.312   -0.715)    0.950
   3.734   (  -3.145   -1.816    1.836)    4.069
   3.876   (   1.716    0.991    0.917)    2.184
   3.941   (   3.782    2.184   -2.062)    4.830
   5.667   (  -0.667   -0.385    1.138)    1.374
   5.731   (  -1.080   -0.623   -1.165)    1.706
   5.789   (   2.944    1.700   -0.721)    3.475
   5.809   (   4.435    2.560   -2.961)    5.915
======================= Grid point 33 (8/15) =======================
q-point: ( 0.40  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 30
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.626   (   9.082    5.243    3.143)   10.948
   0.657   (   7.860    4.538    8.332)   12.321
   0.802   (  -3.207   -1.852   -1.375)    3.951
   0.857   (  -4.936   -2.850    3.550)    6.715
   0.860   (   1.538    0.888   -2.113)    2.760
   0.890   (   1.200    0.693   -1.122)    1.783
   0.981   (   5.329    3.077   -2.629)    6.691
   1.001   (   0.251    0.145    1.815)    1.838
   1.038   (   8.896    5.136    3.333)   10.799
   1.073   (  -0.219   -0.126   -1.660)    1.679
   1.099   (   0.992    0.573    1.221)    1.674
   1.112   (   0.199    0.115    0.993)    1.019
   1.143   (   1.795    1.037    1.807)    2.750
   1.153   (  -0.489   -0.282   -1.881)    1.964
   1.159   (  -0.152   -0.087    1.424)    1.435
   1.189   (   0.315    0.182    1.172)    1.227
   1.209   (  -0.828   -0.478   -4.233)    4.339
   1.210   (  -1.690   -0.976   -0.870)    2.137
   1.237   (  -1.319   -0.762    0.275)    1.548
   1.237   (  -4.176   -2.411    0.620)    4.861
   1.269   (  -3.627   -2.094   -1.940)    4.615
   1.272   (   3.591    2.073    0.533)    4.181
   1.338   (  -1.135   -0.655   -1.678)    2.130
   1.358   (  -3.180   -1.836   -1.627)    4.016
   1.373   (  -0.387   -0.224    1.547)    1.610
   1.387   (   0.784    0.452   -0.218)    0.931
   1.405   (   0.301    0.174   -0.805)    0.877
   1.416   (  -1.743   -1.006   -0.114)    2.016
   1.445   (   0.018    0.010   -0.563)    0.563
   1.453   (  -5.158   -2.978   -1.692)    6.192
   1.515   (   2.179    1.258    0.440)    2.554
   1.552   (   1.155    0.667   -0.478)    1.416
   1.565   (   0.983    0.568   -0.183)    1.150
   1.567   (  -0.334   -0.193    0.130)    0.407
   1.607   (   1.457    0.841    1.703)    2.394
   1.633   (   0.322    0.186    3.105)    3.127
   1.643   (   0.780    0.450    1.362)    1.633
   1.671   (   2.936    1.695    0.872)    3.501
   1.695   (   3.726    2.151   -0.287)    4.312
   1.701   (   0.846    0.489   -0.504)    1.100
   1.735   (  -1.553   -0.896   -1.862)    2.585
   1.768   (   2.742    1.583   -1.676)    3.582
   1.827   (   1.272    0.734    0.402)    1.523
   1.836   (  -1.121   -0.647    1.141)    1.725
   1.836   (   1.172    0.677    0.146)    1.362
   1.857   (   1.002    0.578   -0.051)    1.158
   1.860   (  -1.138   -0.657   -0.064)    1.316
   1.874   (  -0.963   -0.556    1.060)    1.536
   1.896   (   1.957    1.130    0.694)    2.364
   1.939   (   0.509    0.294    0.042)    0.590
   1.960   (  -1.337   -0.772    4.214)    4.488
   2.056   (   1.971    1.138   -0.996)    2.484
   2.116   (  -0.031   -0.018    1.609)    1.610
   2.149   (   0.521    0.301   -0.045)    0.604
   2.168   (   0.956    0.552    1.710)    2.035
   2.290   (  -4.663   -2.692   -2.162)    5.802
   2.323   (  -4.916   -2.838    0.256)    5.683
   2.498   (   3.479    2.009   -5.523)    6.830
   2.733   (  -5.174   -2.987    1.126)    6.080
   2.805   (  -3.712   -2.143   -1.467)    4.531
   2.955   (  -3.540   -2.044    0.835)    4.172
   2.995   (  -0.239   -0.138    0.143)    0.310
   3.002   (  -0.935   -0.540    0.261)    1.110
   3.068   (  -3.094   -1.786   -2.499)    4.360
   3.161   (   0.158    0.091   -0.012)    0.183
   3.174   (  -0.399   -0.230   -0.451)    0.645
   3.337   (   0.005    0.003    1.279)    1.279
   3.341   (   0.105    0.061    1.249)    1.255
   3.386   (  -0.071   -0.041   -0.574)    0.580
   3.389   (   0.028    0.016   -0.583)    0.584
   3.399   (   0.365    0.211    0.316)    0.527
   3.425   (  -0.162   -0.093   -0.768)    0.790
   3.532   (   2.104    1.215    0.289)    2.447
   3.550   (   1.962    1.133   -0.340)    2.290
   3.661   (  -0.143   -0.083    1.860)    1.867
   3.685   (  -0.081   -0.047    0.419)    0.429
   3.696   (  -0.247   -0.142    0.446)    0.529
   3.716   (  -0.240   -0.139   -0.380)    0.470
   3.923   (   1.362    0.786    1.281)    2.028
   4.022   (   2.022    1.167   -2.294)    3.273
   5.653   (  -0.325   -0.188    0.907)    0.981
   5.704   (  -0.774   -0.447   -0.944)    1.300
   5.849   (   1.449    0.837   -0.211)    1.686
   5.871   (   1.149    0.663   -1.111)    1.730
======================= Grid point 37 (9/15) =======================
q-point: ( 0.20  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 30
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.582   (   5.970   10.340    3.938)   12.573
   0.599   (   4.904    8.494    8.168)   12.764
   0.828   (   2.067    3.580   -4.759)    6.303
   0.868   (   1.047    1.813   -0.929)    2.290
   0.879   (   2.023    3.503   -0.079)    4.046
   0.912   (   1.649    2.857   -2.936)    4.416
   0.914   (   1.754    3.038    5.031)    6.133
   0.970   (   2.342    4.057    5.368)    7.125
   0.995   (   1.030    1.784    0.703)    2.176
   1.081   (  -1.382   -2.394    2.068)    3.452
   1.089   (   0.457    0.791   -0.598)    1.092
   1.100   (  -0.991   -1.717   -0.630)    2.080
   1.124   (  -1.267   -2.194    0.151)    2.538
   1.146   (   0.554    0.960    3.124)    3.315
   1.195   (  -0.483   -0.837   -0.885)    1.310
   1.197   (   1.036    1.794    1.151)    2.369
   1.206   (  -0.370   -0.641   -2.004)    2.137
   1.220   (  -0.810   -1.403   -4.273)    4.569
   1.229   (   0.561    0.972    0.105)    1.127
   1.240   (  -1.429   -2.475    0.677)    2.937
   1.259   (  -1.279   -2.216    1.246)    2.846
   1.300   (  -0.998   -1.729   -1.664)    2.599
   1.360   (  -1.387   -2.403    0.936)    2.928
   1.365   (  -0.688   -1.192   -2.737)    3.064
   1.377   (  -1.418   -2.456   -0.261)    2.848
   1.401   (  -0.146   -0.252   -2.599)    2.616
   1.420   (   1.090    1.887    0.183)    2.187
   1.429   (   0.183    0.317   -0.603)    0.705
   1.462   (   0.479    0.830    0.273)    0.996
   1.486   (  -0.054   -0.093   -2.137)    2.139
   1.495   (  -0.717   -1.242    0.253)    1.457
   1.511   (   1.311    2.271   -0.089)    2.623
   1.566   (   2.060    3.569    0.532)    4.155
   1.570   (  -0.557   -0.965    0.345)    1.167
   1.593   (   1.267    2.194   -0.625)    2.609
   1.602   (   0.227    0.394    2.179)    2.226
   1.632   (  -0.130   -0.225    1.054)    1.085
   1.655   (   2.090    3.621    0.895)    4.275
   1.703   (   0.894    1.549    0.417)    1.836
   1.709   (  -1.918   -3.321   -0.570)    3.877
   1.720   (   0.829    1.435   -1.020)    1.946
   1.748   (   1.687    2.922   -1.113)    3.553
   1.794   (   0.770    1.333    0.462)    1.607
   1.809   (  -0.967   -1.674    1.352)    2.359
   1.811   (  -0.785   -1.359   -0.404)    1.621
   1.850   (   1.722    2.983   -0.303)    3.457
   1.877   (   1.974    3.420   -1.210)    4.130
   1.886   (  -1.449   -2.509    1.884)    3.456
   1.930   (   1.651    2.860    0.894)    3.422
   1.949   (   0.052    0.090    0.948)    0.954
   1.983   (  -0.889   -1.540    2.811)    3.326
   2.024   (   1.048    1.816   -1.755)    2.734
   2.090   (   0.142    0.246    3.586)    3.597
   2.143   (  -2.391   -4.142    2.440)    5.369
   2.195   (   1.495    2.590    1.841)    3.512
   2.325   (  -3.988   -6.907   -0.660)    8.002
   2.355   (  -1.045   -1.810   -2.932)    3.600
   2.459   (   1.490    2.581   -5.423)    6.187
   2.795   (  -2.811   -4.868    1.240)    5.757
   2.865   (  -2.279   -3.948   -1.218)    4.719
   2.992   (  -0.477   -0.827   -0.032)    0.955
   2.996   (  -1.183   -2.049   -0.455)    2.410
   3.032   (  -2.288   -3.963    0.495)    4.603
   3.045   (  -2.414   -4.180   -0.388)    4.843
   3.086   (  -2.079   -3.601    0.373)    4.175
   3.161   (  -2.279   -3.947   -2.592)    5.244
   3.337   (   0.046    0.079    1.217)    1.221
   3.341   (   0.005    0.008    1.283)    1.283
   3.387   (  -0.037   -0.064   -0.639)    0.643
   3.389   (  -0.038   -0.066   -0.584)    0.589
   3.440   (   1.436    2.488    0.439)    2.906
   3.466   (   1.092    1.891   -0.585)    2.261
   3.503   (   1.852    3.208    0.293)    3.716
   3.520   (   1.797    3.113   -0.343)    3.611
   3.662   (  -0.121   -0.209    2.064)    2.078
   3.685   (   0.038    0.066    0.380)    0.388
   3.700   (  -0.375   -0.650    0.917)    1.185
   3.716   (  -0.260   -0.450   -0.621)    0.810
   3.911   (   1.050    1.818    1.191)    2.414
   4.006   (   1.726    2.989   -2.302)    4.149
   5.669   (   0.356    0.616    0.906)    1.152
   5.722   (  -0.115   -0.200   -0.966)    0.993
   5.828   (   0.828    1.434   -0.312)    1.685
   5.848   (   0.437    0.758   -1.391)    1.643
======================= Grid point 38 (10/15) =======================
q-point: ( 0.40  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 30
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.742   (  -1.295    2.244    7.331)    7.776
   0.769   (  -1.963    3.399    1.611)    4.243
   0.815   (  -2.929    5.073    1.093)    5.959
   0.828   (  -1.016    1.760    3.806)    4.315
   0.886   (   0.120   -0.208   -1.169)    1.194
   0.927   (  -2.347    4.064   -2.192)    5.180
   0.944   (   1.896   -3.285    3.517)    5.173
   1.018   (   3.041   -5.266   -0.084)    6.082
   1.067   (   2.401   -4.159   -1.559)    5.049
   1.080   (   0.053   -0.091   -0.915)    0.921
   1.083   (  -1.970    3.412   -0.775)    4.015
   1.102   (   1.041   -1.802    0.135)    2.086
   1.117   (   0.092   -0.160   -4.079)    4.083
   1.169   (   0.722   -1.250    1.588)    2.146
   1.187   (  -0.664    1.151   -0.209)    1.345
   1.192   (  -1.177    2.039    0.215)    2.364
   1.214   (   0.118   -0.968   -2.436)    2.624
   1.214   (  -0.720    2.012    1.792)    2.789
   1.224   (  -0.286    0.496   -1.185)    1.316
   1.236   (  -0.984    1.704   -0.508)    2.032
   1.247   (  -0.451    0.781   -2.915)    3.051
   1.274   (   0.546   -0.945    0.090)    1.095
   1.321   (   1.023   -1.772    1.002)    2.279
   1.329   (   0.788   -1.365   -0.209)    1.590
   1.342   (   0.184   -0.319   -1.234)    1.288
   1.392   (  -0.462    0.800   -1.868)    2.084
   1.438   (  -1.149    1.990    0.851)    2.450
   1.442   (  -1.858    3.218    1.978)    4.210
   1.454   (  -0.782    1.355   -0.196)    1.577
   1.491   (   1.332   -2.306    1.120)    2.889
   1.495   (   1.623   -2.811    0.506)    3.285
   1.496   (  -1.612    2.793    0.839)    3.332
   1.567   (  -0.754    1.305    1.151)    1.897
   1.573   (  -0.272    0.470    1.925)    2.000
   1.615   (   0.933   -1.616    1.230)    2.235
   1.625   (   3.172   -5.493    0.700)    6.382
   1.656   (   0.605   -1.047   -0.471)    1.298
   1.661   (   0.143   -0.248   -1.507)    1.534
   1.701   (  -1.033    1.789    0.589)    2.148
   1.714   (  -0.423    0.733    0.800)    1.164
   1.744   (   0.376   -0.652   -0.619)    0.974
   1.774   (  -0.166    0.288   -0.914)    0.972
   1.794   (   1.411   -2.444    1.150)    3.048
   1.796   (   0.990   -1.714   -0.159)    1.986
   1.824   (  -1.335    2.313    0.234)    2.681
   1.863   (   0.113   -0.196   -0.362)    0.427
   1.883   (  -0.511    0.884    0.079)    1.024
   1.904   (  -0.705    1.221   -0.288)    1.439
   1.935   (  -0.238    0.412    2.862)    2.901
   1.985   (  -1.989    3.445    1.155)    4.142
   1.998   (  -1.881    3.258    0.361)    3.780
   2.050   (   0.992   -1.717   -1.158)    2.296
   2.085   (   0.807   -1.397    1.165)    1.990
   2.105   (   1.650   -2.857    0.493)    3.336
   2.204   (   1.510   -2.615    2.304)    3.798
   2.257   (  -2.142    3.709   -0.956)    4.389
   2.314   (  -0.435    0.754   -1.316)    1.578
   2.474   (   1.674   -2.899   -5.272)    6.245
   2.772   (  -3.093    5.357    0.201)    6.189
   2.817   (  -1.647    2.853   -1.096)    3.472
   2.923   (   0.783   -1.356    0.502)    1.644
   2.972   (   2.699   -4.674   -1.159)    5.520
   2.995   (   0.005   -0.009    0.426)    0.426
   3.003   (   0.856   -1.483    0.269)    1.734
   3.073   (   2.461   -4.263   -0.425)    4.941
   3.100   (   2.333   -4.041   -1.598)    4.932
   3.338   (  -0.059    0.102    1.207)    1.213
   3.341   (   0.061   -0.106    1.241)    1.247
   3.386   (  -0.026    0.045   -0.595)    0.597
   3.388   (   0.050   -0.086   -0.566)    0.575
   3.463   (  -1.793    3.106    0.306)    3.599
   3.479   (  -1.429    2.475   -0.290)    2.873
   3.545   (  -0.224    0.387    0.311)    0.545
   3.560   (   0.090   -0.157   -0.217)    0.283
   3.660   (   0.019   -0.034    1.630)    1.630
   3.685   (   0.004   -0.007    0.491)    0.491
   3.698   (  -0.182    0.315   -0.053)    0.368
   3.709   (   0.226   -0.392   -0.363)    0.580
   3.933   (  -0.076    0.132    1.491)    1.498
   4.040   (  -0.237    0.411   -2.310)    2.359
   5.680   (  -1.149    1.991    0.647)    2.388
   5.719   (  -0.780    1.351   -0.771)    1.740
   5.842   (   0.585   -1.012   -0.208)    1.187
   5.854   (   1.023   -1.772   -0.640)    2.144
======================= Grid point 62 (11/15) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 10
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.539   (   0.000   -0.000  -11.219)   11.219
   0.539   (   0.000    0.000  -11.219)   11.219
   0.539   (   0.000    0.000   11.219)   11.219
   0.539   (   0.000    0.000   11.219)   11.219
   0.770   (  -0.000    0.000   -0.589)    0.589
   0.770   (   0.000   -0.000    0.589)    0.589
   1.054   (   0.000   -0.000  -21.096)   21.096
   1.054   (   0.000   -0.000   21.096)   21.096
   1.087   (  -0.000   -0.000   -1.257)    1.257
   1.087   (  -0.000    0.000   -1.257)    1.257
   1.087   (  -0.000   -0.000    1.257)    1.257
   1.087   (   0.000    0.000    1.257)    1.257
   1.175   (   0.000   -0.000   -3.727)    3.727
   1.175   (   0.000    0.000   -3.727)    3.727
   1.175   (  -0.000   -0.000    3.727)    3.727
   1.175   (  -0.000    0.000    3.727)    3.727
   1.236   (  -0.000   -0.000   -2.663)    2.663
   1.236   (  -0.000    0.000    2.663)    2.663
   1.250   (   0.000    0.000   -0.932)    0.932
   1.250   (   0.000   -0.000   -0.932)    0.932
   1.250   (   0.000   -0.000    0.932)    0.932
   1.250   (   0.000   -0.000    0.932)    0.932
   1.358   (   0.000    0.000   -0.408)    0.408
   1.358   (   0.000    0.000   -0.408)    0.408
   1.358   (  -0.000   -0.000    0.408)    0.408
   1.358   (  -0.000   -0.000    0.408)    0.408
   1.392   (  -0.000    0.000   -3.727)    3.727
   1.392   (  -0.000    0.000    3.727)    3.727
   1.474   (   0.000    0.000   -1.094)    1.094
   1.474   (   0.000    0.000   -1.094)    1.094
   1.474   (  -0.000    0.000    1.094)    1.094
   1.474   (   0.000    0.000    1.094)    1.094
   1.529   (  -0.000    0.000   -0.154)    0.154
   1.529   (   0.000   -0.000    0.154)    0.154
   1.535   (   0.000    0.000   -2.046)    2.046
   1.535   (   0.000    0.000   -2.046)    2.046
   1.535   (  -0.000    0.000    2.046)    2.046
   1.535   (  -0.000    0.000    2.046)    2.046
   1.721   (   0.000   -0.000   -2.465)    2.465
   1.721   (  -0.000    0.000    2.465)    2.465
   1.724   (   0.000    0.000   -0.398)    0.398
   1.724   (   0.000    0.000   -0.398)    0.398
   1.724   (  -0.000    0.000    0.398)    0.398
   1.724   (   0.000   -0.000    0.398)    0.398
   1.855   (   0.000    0.000   -0.465)    0.465
   1.855   (   0.000    0.000    0.465)    0.465
   1.915   (   0.000    0.000   -0.103)    0.103
   1.915   (   0.000    0.000   -0.103)    0.103
   1.915   (   0.000   -0.000    0.103)    0.103
   1.915   (   0.000   -0.000    0.103)    0.103
   2.050   (   0.000   -0.000   -1.085)    1.085
   2.050   (   0.000   -0.000    1.085)    1.085
   2.295   (   0.000   -0.000   -5.861)    5.861
   2.295   (   0.000    0.000   -5.861)    5.861
   2.295   (   0.000   -0.000    5.861)    5.861
   2.295   (   0.000   -0.000    5.861)    5.861
   2.365   (   0.000    0.000   -0.222)    0.222
   2.365   (  -0.000    0.000    0.222)    0.222
   3.033   (   0.000    0.000   -0.904)    0.904
   3.033   (   0.000    0.000    0.904)    0.904
   3.073   (   0.000    0.000   -0.784)    0.784
   3.073   (   0.000   -0.000    0.784)    0.784
   3.205   (   0.000    0.000   -0.979)    0.979
   3.205   (   0.000    0.000   -0.979)    0.979
   3.205   (  -0.000    0.000    0.979)    0.979
   3.205   (   0.000    0.000    0.979)    0.979
   3.330   (   0.000   -0.000   -0.195)    0.195
   3.330   (   0.000    0.000   -0.195)    0.195
   3.330   (   0.000    0.000    0.195)    0.195
   3.330   (  -0.000   -0.000    0.195)    0.195
   3.372   (   0.000    0.000   -1.297)    1.297
   3.372   (   0.000    0.000   -1.297)    1.297
   3.372   (  -0.000   -0.000    1.297)    1.297
   3.372   (   0.000   -0.000    1.297)    1.297
   3.693   (   0.000    0.000   -1.397)    1.397
   3.693   (   0.000   -0.000    1.397)    1.397
   3.736   (   0.000    0.000   -3.467)    3.467
   3.736   (   0.000    0.000    3.467)    3.467
   3.893   (   0.000   -0.000   -1.391)    1.391
   3.893   (   0.000    0.000    1.391)    1.391
   5.711   (   0.000    0.000   -1.752)    1.752
   5.711   (   0.000    0.000   -1.752)    1.752
   5.711   (   0.000    0.000    1.752)    1.752
   5.711   (   0.000    0.000    1.752)    1.752
======================= Grid point 64 (12/15) =======================
q-point: ( 0.20  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 24
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.594   (   4.226    2.440   -9.153)   10.372
   0.594   (   4.226    2.440    9.153)   10.372
   0.639   (   7.139    4.121   -9.007)   12.209
   0.639   (   7.139    4.121    9.007)   12.209
   0.816   (   3.229    1.864   -0.388)    3.749
   0.816   (   3.229    1.864    0.388)    3.749
   1.064   (  -1.871   -1.080   -2.331)    3.178
   1.064   (  -1.871   -1.080    2.331)    3.178
   1.079   (  -0.565   -0.326   -6.420)    6.453
   1.079   (  -0.565   -0.326    6.420)    6.453
   1.093   (   2.109    1.217  -10.256)   10.541
   1.093   (   2.109    1.217   10.256)   10.541
   1.162   (  -0.586   -0.339   -1.939)    2.054
   1.162   (  -0.586   -0.339    1.939)    2.054
   1.185   (   0.486    0.281   -2.140)    2.212
   1.185   (   0.486    0.281    2.140)    2.212
   1.202   (  -2.030   -1.172   -2.140)    3.174
   1.202   (  -2.030   -1.172    2.140)    3.174
   1.226   (  -0.477   -0.276   -1.294)    1.407
   1.226   (  -0.477   -0.276    1.294)    1.407
   1.295   (   1.845    1.065   -1.447)    2.575
   1.295   (   1.845    1.065    1.447)    2.575
   1.342   (  -2.687   -1.551   -2.867)    4.225
   1.342   (  -2.687   -1.551    2.867)    4.225
   1.363   (   0.272    0.157   -0.042)    0.316
   1.363   (   0.272    0.157    0.042)    0.316
   1.401   (   2.659    1.535   -0.504)    3.111
   1.401   (   2.659    1.535    0.504)    3.111
   1.458   (  -2.368   -1.367   -0.504)    2.781
   1.458   (  -2.368   -1.367    0.504)    2.781
   1.510   (   2.110    1.218   -0.851)    2.581
   1.510   (   2.110    1.218    0.851)    2.581
   1.523   (   1.382    0.798   -1.753)    2.370
   1.523   (   1.382    0.798    1.753)    2.370
   1.580   (   3.394    1.960   -0.804)    4.001
   1.580   (   3.394    1.960    0.804)    4.001
   1.608   (   3.333    1.924   -1.769)    4.236
   1.608   (   3.333    1.924    1.769)    4.236
   1.686   (  -1.339   -0.773   -2.022)    2.545
   1.686   (  -1.339   -0.773    2.022)    2.545
   1.768   (   3.220    1.859   -0.532)    3.756
   1.768   (   3.220    1.859    0.532)    3.756
   1.769   (   0.920    0.531   -0.343)    1.117
   1.769   (   0.920    0.531    0.343)    1.117
   1.789   (  -1.536   -0.887   -0.658)    1.892
   1.789   (  -1.536   -0.887    0.658)    1.892
   1.915   (  -0.093   -0.054   -1.260)    1.265
   1.915   (  -0.093   -0.054    1.260)    1.265
   1.923   (  -0.676   -0.390   -2.391)    2.515
   1.923   (  -0.676   -0.390    2.391)    2.515
   2.020   (  -1.251   -0.722   -1.363)    1.986
   2.020   (  -1.251   -0.722    1.363)    1.986
   2.264   (  -5.384   -3.109   -4.647)    7.762
   2.264   (  -5.384   -3.109    4.647)    7.762
   2.280   (  -1.351   -0.780   -5.379)    5.601
   2.280   (  -1.351   -0.780    5.379)    5.601
   2.338   (   1.438    0.830   -3.854)    4.196
   2.338   (   1.438    0.830    3.854)    4.196
   2.943   (  -8.413   -4.857   -0.932)    9.759
   2.943   (  -8.413   -4.857    0.932)    9.759
   3.014   (  -1.165   -0.672   -0.348)    1.389
   3.014   (  -1.165   -0.672    0.348)    1.389
   3.120   (  -4.868   -2.811   -1.732)    5.882
   3.120   (  -4.868   -2.811    1.732)    5.882
   3.172   (  -0.651   -0.376   -0.722)    1.042
   3.172   (  -0.651   -0.376    0.722)    1.042
   3.364   (   0.486    0.281   -1.116)    1.249
   3.364   (   0.486    0.281    1.116)    1.249
   3.368   (   0.004    0.002   -1.275)    1.275
   3.368   (   0.004    0.002    1.275)    1.275
   3.391   (   1.473    0.850   -0.816)    1.886
   3.391   (   1.473    0.850    0.816)    1.886
   3.447   (   5.355    3.092   -0.582)    6.211
   3.447   (   5.355    3.092    0.582)    6.211
   3.700   (   0.042    0.024   -1.083)    1.084
   3.700   (   0.042    0.024    1.083)    1.084
   3.727   (  -1.317   -0.760   -2.190)    2.667
   3.727   (  -1.317   -0.760    2.190)    2.667
   3.902   (   2.190    1.265   -1.394)    2.888
   3.902   (   2.190    1.265    1.394)    2.888
   5.699   (  -0.876   -0.506   -1.628)    1.917
   5.699   (  -0.876   -0.506    1.628)    1.917
   5.779   (   3.975    2.295   -0.219)    4.595
   5.779   (   3.975    2.295    0.219)    4.595
======================= Grid point 66 (13/15) =======================
q-point: ( 0.40  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 24
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.720   (   5.202    3.003   -4.963)    7.792
   0.720   (   5.202    3.003    4.963)    7.792
   0.826   (   6.834    3.946   -6.369)   10.141
   0.826   (   6.834    3.946    6.369)   10.141
   0.860   (  -1.243   -0.717   -1.429)    2.025
   0.860   (  -1.243   -0.717    1.429)    2.025
   1.027   (  -3.176   -1.833   -6.559)    7.515
   1.027   (  -3.176   -1.833    6.559)    7.515
   1.037   (   0.560    0.323   -1.484)    1.619
   1.037   (   0.560    0.323    1.484)    1.619
   1.109   (  -1.084   -0.626   -2.526)    2.819
   1.109   (  -1.084   -0.626    2.526)    2.819
   1.145   (   0.102    0.059   -2.441)    2.444
   1.145   (   0.102    0.059    2.441)    2.444
   1.183   (   0.339    0.196   -0.748)    0.844
   1.183   (   0.339    0.196    0.748)    0.844
   1.213   (   2.117    1.222   -1.184)    2.716
   1.213   (   2.117    1.222    1.184)    2.716
   1.249   (  -4.379   -2.528   -0.320)    5.066
   1.249   (  -4.379   -2.528    0.320)    5.066
   1.265   (   0.445    0.257   -1.748)    1.822
   1.265   (   0.445    0.257    1.748)    1.822
   1.296   (   1.446    0.835   -1.598)    2.311
   1.296   (   1.446    0.835    1.598)    2.311
   1.387   (   0.879    0.507   -0.682)    1.223
   1.387   (   0.879    0.507    0.682)    1.223
   1.417   (  -0.579   -0.335   -1.725)    1.851
   1.417   (  -0.579   -0.335    1.726)    1.851
   1.419   (  -2.661   -1.536   -1.306)    3.339
   1.419   (  -2.661   -1.536    1.306)    3.339
   1.547   (   2.403    1.387   -1.966)    3.400
   1.547   (   2.403    1.387    1.966)    3.400
   1.553   (   0.605    0.349   -0.618)    0.933
   1.553   (   0.605    0.349    0.618)    0.933
   1.653   (   0.417    0.240   -1.228)    1.319
   1.653   (   0.417    0.240    1.228)    1.319
   1.687   (   3.763    2.173   -0.494)    4.374
   1.687   (   3.763    2.173    0.494)    4.374
   1.697   (  -0.573   -0.331   -1.724)    1.846
   1.697   (  -0.573   -0.331    1.724)    1.846
   1.719   (   0.409    0.236   -2.271)    2.319
   1.719   (   0.409    0.236    2.271)    2.319
   1.840   (   1.402    0.809   -0.664)    1.750
   1.840   (   1.402    0.809    0.664)    1.750
   1.849   (   1.824    1.053   -0.663)    2.208
   1.849   (   1.824    1.053    0.663)    2.208
   1.873   (  -2.381   -1.375   -1.203)    3.001
   1.873   (  -2.381   -1.375    1.203)    3.001
   1.921   (   1.341    0.774   -1.408)    2.093
   1.921   (   1.341    0.774    1.408)    2.093
   2.027   (   1.196    0.690   -1.900)    2.349
   2.027   (   1.196    0.690    1.900)    2.349
   2.148   (  -2.946   -1.701   -0.758)    3.485
   2.148   (  -2.946   -1.701    0.758)    3.485
   2.227   (  -2.518   -1.454   -3.247)    4.358
   2.227   (  -2.518   -1.454    3.247)    4.358
   2.371   (   0.839    0.484   -4.679)    4.778
   2.371   (   0.839    0.484    4.679)    4.778
   2.767   (  -4.518   -2.609   -1.796)    5.518
   2.767   (  -4.518   -2.609    1.796)    5.518
   2.983   (  -1.456   -0.841   -1.328)    2.143
   2.983   (  -1.456   -0.841    1.328)    2.143
   3.019   (  -2.160   -1.247   -1.440)    2.880
   3.019   (  -2.160   -1.247    1.440)    2.880
   3.165   (  -0.101   -0.058   -0.338)    0.358
   3.165   (  -0.101   -0.058    0.338)    0.358
   3.367   (  -0.037   -0.022   -1.220)    1.221
   3.367   (  -0.037   -0.022    1.220)    1.221
   3.369   (   0.094    0.054   -1.161)    1.166
   3.369   (   0.094    0.054    1.161)    1.166
   3.409   (   0.192    0.111   -0.594)    0.634
   3.409   (   0.192    0.111    0.594)    0.634
   3.541   (   2.035    1.175   -0.446)    2.392
   3.541   (   2.035    1.175    0.446)    2.392
   3.689   (  -0.537   -0.310   -0.416)    0.747
   3.689   (  -0.537   -0.310    0.416)    0.747
   3.706   (  -0.186   -0.107   -0.623)    0.659
   3.706   (  -0.186   -0.107    0.623)    0.659
   3.966   (   1.806    1.043   -2.504)    3.258
   3.966   (   1.806    1.043    2.504)    3.258
   5.678   (  -0.551   -0.318   -1.308)    1.455
   5.678   (  -0.551   -0.318    1.308)    1.455
   5.851   (   1.593    0.920   -0.529)    1.914
   5.851   (   1.593    0.920    0.529)    1.914
======================= Grid point 74 (14/15) =======================
q-point: ( 0.20  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 24
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.696   (   3.726    6.453   -5.692)    9.376
   0.696   (   3.726    6.453    5.692)    9.376
   0.772   (   4.464    7.732   -6.534)   11.063
   0.772   (   4.464    7.732    6.534)   11.063
   0.877   (   1.587    2.750   -0.027)    3.175
   0.877   (   1.587    2.750    0.027)    3.175
   1.040   (  -1.217   -2.107   -2.944)    3.819
   1.040   (  -1.217   -2.107    2.944)    3.819
   1.052   (  -0.577   -1.000   -0.238)    1.179
   1.052   (  -0.577   -1.000    0.238)    1.179
   1.093   (  -1.177   -2.038   -3.703)    4.388
   1.093   (  -1.177   -2.038    3.703)    4.388
   1.151   (  -0.817   -1.415   -0.724)    1.787
   1.151   (  -0.817   -1.415    0.724)    1.787
   1.171   (   0.167    0.290   -2.573)    2.595
   1.171   (   0.167    0.290    2.573)    2.595
   1.205   (   0.839    1.453   -1.318)    2.133
   1.205   (   0.839    1.453    1.318)    2.133
   1.241   (   0.574    0.994   -1.821)    2.152
   1.241   (   0.574    0.994    1.821)    2.152
   1.267   (  -1.708   -2.958   -1.406)    3.694
   1.267   (  -1.708   -2.958    1.406)    3.694
   1.330   (   0.139    0.242   -0.025)    0.280
   1.330   (   0.139    0.242    0.025)    0.280
   1.378   (  -1.298   -2.248   -0.054)    2.596
   1.378   (  -1.298   -2.248    0.054)    2.596
   1.416   (   1.535    2.659   -0.058)    3.071
   1.416   (   1.535    2.659    0.058)    3.071
   1.461   (   0.890    1.542   -0.070)    1.782
   1.461   (   0.890    1.542    0.070)    1.782
   1.506   (   0.052    0.089   -0.567)    0.576
   1.506   (   0.052    0.089    0.567)    0.576
   1.573   (   1.057    1.830   -0.091)    2.116
   1.573   (   1.057    1.830    0.091)    2.116
   1.602   (   0.278    0.482   -1.229)    1.349
   1.602   (   0.278    0.482    1.229)    1.349
   1.682   (  -0.016   -0.027   -1.957)    1.957
   1.682   (  -0.016   -0.027    1.957)    1.957
   1.693   (   0.676    1.171   -1.458)    1.988
   1.693   (   0.676    1.171    1.458)    1.988
   1.721   (   1.490    2.581   -1.103)    3.177
   1.721   (   1.490    2.581    1.103)    3.178
   1.803   (  -0.090   -0.156   -0.347)    0.390
   1.803   (  -0.090   -0.156    0.347)    0.390
   1.845   (   1.593    2.760   -0.383)    3.209
   1.845   (   1.593    2.760    0.383)    3.209
   1.874   (  -1.225   -2.123   -1.914)    3.110
   1.874   (  -1.225   -2.123    1.914)    3.110
   1.951   (   1.315    2.278   -0.674)    2.716
   1.951   (   1.315    2.278    0.674)    2.716
   2.007   (   0.070    0.122   -0.085)    0.164
   2.007   (   0.070    0.122    0.085)    0.164
   2.167   (  -2.681   -4.644   -1.347)    5.529
   2.167   (  -2.681   -4.644    1.347)    5.529
   2.269   (  -0.425   -0.736   -4.363)    4.445
   2.269   (  -0.425   -0.736    4.363)    4.445
   2.342   (  -0.222   -0.384   -3.513)    3.540
   2.342   (  -0.222   -0.384    3.513)    3.540
   2.830   (  -2.496   -4.324   -1.793)    5.305
   2.830   (  -2.496   -4.324    1.793)    5.305
   2.991   (  -1.113   -1.928   -0.012)    2.227
   2.991   (  -1.113   -1.928    0.012)    2.227
   3.039   (  -2.160   -3.741   -0.185)    4.323
   3.039   (  -2.160   -3.741    0.185)    4.323
   3.109   (  -2.032   -3.519   -1.796)    4.443
   3.109   (  -2.032   -3.519    1.796)    4.443
   3.365   (   0.008    0.015   -0.982)    0.982
   3.365   (   0.008    0.015    0.982)    0.982
   3.370   (   0.002    0.003   -0.890)    0.890
   3.370   (   0.002    0.003    0.890)    0.890
   3.452   (   1.308    2.266   -0.642)    2.694
   3.452   (   1.308    2.266    0.642)    2.694
   3.511   (   1.882    3.259   -0.421)    3.787
   3.511   (   1.882    3.259    0.421)    3.787
   3.694   (  -0.373   -0.646   -0.040)    0.747
   3.694   (  -0.373   -0.646    0.040)    0.747
   3.706   (  -0.324   -0.562   -0.273)    0.704
   3.706   (  -0.324   -0.562    0.273)    0.704
   3.951   (   1.474    2.553   -2.352)    3.771
   3.951   (   1.474    2.553    2.352)    3.771
   5.695   (   0.105    0.181   -1.317)    1.333
   5.695   (   0.105    0.181    1.317)    1.333
   5.827   (   0.911    1.577   -0.376)    1.860
   5.827   (   0.911    1.577    0.376)    1.860
======================= Grid point 76 (15/15) =======================
q-point: ( 0.40  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.97e-04 1.97e-04 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 
Number of triplets: 24
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.807   (  -1.851    3.206   -0.843)    3.797
   0.807   (  -1.851    3.206    0.843)    3.797
   0.869   (  -1.501    2.600   -1.702)    3.451
   0.869   (  -1.501    2.600    1.702)    3.451
   0.915   (  -1.742    3.017   -0.091)    3.485
   0.915   (  -1.742    3.017    0.091)    3.485
   0.979   (   0.994   -1.722   -2.439)    3.146
   0.979   (   0.994   -1.722    2.439)    3.146
   1.035   (   1.736   -3.006   -0.587)    3.521
   1.035   (   1.736   -3.006    0.587)    3.521
   1.079   (  -0.096    0.166   -0.259)    0.323
   1.079   (  -0.096    0.166    0.259)    0.323
   1.124   (   1.494   -2.588   -2.615)    3.971
   1.124   (   1.494   -2.588    2.615)    3.971
   1.194   (  -0.530    0.917   -0.500)    1.171
   1.194   (  -0.530    0.917    0.500)    1.171
   1.204   (  -0.314    0.544   -0.436)    0.765
   1.204   (  -0.314    0.544    0.436)    0.765
   1.228   (   0.046   -0.080   -0.553)    0.561
   1.228   (   0.046   -0.080    0.553)    0.561
   1.263   (  -0.157    0.273   -1.488)    1.521
   1.263   (  -0.157    0.273    1.488)    1.521
   1.328   (   1.339   -2.319   -0.088)    2.679
   1.328   (   1.339   -2.319    0.088)    2.679
   1.344   (  -0.234    0.404   -1.651)    1.715
   1.344   (  -0.234    0.404    1.651)    1.716
   1.446   (  -1.270    2.200   -0.045)    2.541
   1.446   (  -1.270    2.200    0.045)    2.541
   1.484   (  -2.406    4.167   -1.422)    5.018
   1.484   (  -2.406    4.167    1.422)    5.018
   1.511   (   0.287   -0.497   -0.163)    0.597
   1.511   (   0.287   -0.497    0.163)    0.597
   1.594   (   0.312   -0.541   -1.246)    1.394
   1.594   (   0.312   -0.541    1.246)    1.394
   1.630   (   1.221   -2.116   -0.107)    2.445
   1.630   (   1.221   -2.116    0.107)    2.445
   1.643   (   1.028   -1.780   -0.359)    2.087
   1.643   (   1.028   -1.780    0.359)    2.087
   1.720   (  -1.036    1.794   -0.629)    2.164
   1.720   (  -1.036    1.794    0.629)    2.164
   1.750   (   0.798   -1.383   -1.005)    1.887
   1.750   (   0.798   -1.383    1.005)    1.887
   1.804   (   0.662   -1.146   -0.709)    1.501
   1.804   (   0.662   -1.146    0.709)    1.501
   1.838   (  -0.142    0.245   -1.504)    1.530
   1.838   (  -0.142    0.245    1.504)    1.530
   1.890   (  -0.669    1.158   -0.593)    1.463
   1.890   (  -0.669    1.158    0.593)    1.463
   1.986   (  -1.136    1.968   -1.475)    2.710
   1.986   (  -1.136    1.968    1.475)    2.710
   2.019   (   0.335   -0.580   -1.584)    1.720
   2.019   (   0.335   -0.580    1.584)    1.720
   2.109   (   0.679   -1.176   -0.579)    1.477
   2.109   (   0.679   -1.176    0.579)    1.477
   2.250   (  -1.669    2.891   -0.043)    3.338
   2.250   (  -1.669    2.891    0.043)    3.338
   2.337   (   1.756   -3.041   -5.544)    6.562
   2.337   (   1.756   -3.041    5.544)    6.562
   2.789   (  -2.331    4.038   -1.250)    4.827
   2.789   (  -2.331    4.038    1.250)    4.827
   2.944   (   1.766   -3.058   -1.269)    3.753
   2.944   (   1.766   -3.058    1.269)    3.753
   3.003   (   0.167   -0.289   -0.225)    0.402
   3.003   (   0.167   -0.289    0.225)    0.402
   3.073   (   2.757   -4.776   -0.600)    5.547
   3.073   (   2.757   -4.776    0.600)    5.547
   3.365   (  -0.015    0.027   -0.776)    0.777
   3.365   (  -0.015    0.027    0.776)    0.777
   3.370   (   0.036   -0.062   -0.727)    0.731
   3.370   (   0.036   -0.062    0.727)    0.731
   3.471   (  -1.661    2.877   -0.392)    3.345
   3.471   (  -1.661    2.877    0.392)    3.345
   3.554   (  -0.096    0.166   -0.376)    0.422
   3.554   (  -0.096    0.166    0.376)    0.422
   3.686   (  -0.028    0.049   -0.496)    0.499
   3.686   (  -0.028    0.049    0.496)    0.499
   3.701   (   0.144   -0.250   -0.367)    0.466
   3.701   (   0.144   -0.250    0.367)    0.466
   3.981   (  -0.184    0.318   -2.729)    2.754
   3.981   (  -0.184    0.318    2.729)    2.754
   5.699   (  -0.942    1.631   -0.996)    2.131
   5.699   (  -0.942    1.631    0.996)    2.131
   5.842   (   0.750   -1.300   -0.292)    1.529
   5.842   (   0.750   -1.300    0.292)    1.529
=================== End of collection of collisions ===================
----------- Thermal conductivity (W/m-k) with tetrahedron method -----------
#  T(K)        xx         yy         zz         yz         xz         xy        #ipm
    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/8400
   10.0     13.393     13.393     14.087      0.000     -0.000     -0.000 3/8400
   20.0      2.950      2.950      3.638      0.000     -0.000     -0.000 3/8400
   30.0      1.693      1.693      2.092      0.000     -0.000     -0.000 3/8400
   40.0      1.228      1.228      1.486      0.000     -0.000     -0.000 3/8400
   50.0      0.978      0.978      1.158      0.000     -0.000     -0.000 3/8400
   60.0      0.818      0.818      0.952      0.000     -0.000     -0.000 3/8400
   70.0      0.705      0.705      0.809      0.000     -0.000     -0.000 3/8400
   80.0      0.620      0.620      0.705      0.000     -0.000     -0.000 3/8400
   90.0      0.554      0.554      0.625      0.000     -0.000     -0.000 3/8400
  100.0      0.501      0.501      0.561      0.000     -0.000     -0.000 3/8400
  110.0      0.457      0.457      0.509      0.000     -0.000     -0.000 3/8400
  120.0      0.420      0.420      0.467      0.000     -0.000     -0.000 3/8400
  130.0      0.389      0.389      0.430      0.000     -0.000     -0.000 3/8400
  140.0      0.362      0.362      0.400      0.000     -0.000     -0.000 3/8400
  150.0      0.339      0.339      0.373      0.000     -0.000     -0.000 3/8400
  160.0      0.318      0.318      0.349      0.000     -0.000     -0.000 3/8400
  170.0      0.300      0.300      0.329      0.000     -0.000     -0.000 3/8400
  180.0      0.284      0.284      0.311      0.000     -0.000     -0.000 3/8400
  190.0      0.269      0.269      0.294      0.000     -0.000     -0.000 3/8400
  200.0      0.256      0.256      0.280      0.000     -0.000     -0.000 3/8400
  210.0      0.244      0.244      0.266      0.000     -0.000     -0.000 3/8400
  220.0      0.233      0.233      0.254      0.000     -0.000     -0.000 3/8400
  230.0      0.223      0.223      0.243      0.000     -0.000     -0.000 3/8400
  240.0      0.214      0.214      0.233      0.000     -0.000     -0.000 3/8400
  250.0      0.206      0.206      0.224      0.000     -0.000     -0.000 3/8400
  260.0      0.198      0.198      0.215      0.000     -0.000     -0.000 3/8400
  270.0      0.191      0.191      0.207      0.000     -0.000     -0.000 3/8400
  280.0      0.184      0.184      0.200      0.000     -0.000     -0.000 3/8400
  290.0      0.178      0.178      0.193      0.000     -0.000     -0.000 3/8400
  300.0      0.172      0.172      0.187      0.000     -0.000     -0.000 3/8400
  310.0      0.167      0.167      0.181      0.000     -0.000     -0.000 3/8400
  320.0      0.162      0.162      0.175      0.000     -0.000     -0.000 3/8400
  330.0      0.157      0.157      0.170      0.000     -0.000     -0.000 3/8400
  340.0      0.152      0.152      0.165      0.000     -0.000     -0.000 3/8400
  350.0      0.148      0.148      0.160      0.000     -0.000     -0.000 3/8400
  360.0      0.144      0.144      0.156      0.000     -0.000     -0.000 3/8400
  370.0      0.140      0.140      0.151      0.000     -0.000     -0.000 3/8400
  380.0      0.136      0.136      0.147      0.000     -0.000     -0.000 3/8400
  390.0      0.133      0.133      0.144      0.000     -0.000     -0.000 3/8400
  400.0      0.130      0.130      0.140      0.000     -0.000     -0.000 3/8400
  410.0      0.127      0.127      0.137      0.000     -0.000     -0.000 3/8400
  420.0      0.124      0.124      0.133      0.000     -0.000     -0.000 3/8400
  430.0      0.121      0.121      0.130      0.000     -0.000     -0.000 3/8400
  440.0      0.118      0.118      0.127      0.000     -0.000     -0.000 3/8400
  450.0      0.115      0.115      0.125      0.000     -0.000     -0.000 3/8400
  460.0      0.113      0.113      0.122      0.000     -0.000     -0.000 3/8400
  470.0      0.111      0.111      0.119      0.000     -0.000     -0.000 3/8400
  480.0      0.108      0.108      0.117      0.000     -0.000     -0.000 3/8400
  490.0      0.106      0.106      0.114      0.000     -0.000     -0.000 3/8400
  500.0      0.104      0.104      0.112      0.000     -0.000     -0.000 3/8400
  510.0      0.102      0.102      0.110      0.000     -0.000     -0.000 3/8400
  520.0      0.100      0.100      0.108      0.000     -0.000     -0.000 3/8400
  530.0      0.098      0.098      0.106      0.000     -0.000     -0.000 3/8400
  540.0      0.096      0.096      0.104      0.000     -0.000     -0.000 3/8400
  550.0      0.095      0.095      0.102      0.000     -0.000     -0.000 3/8400
  560.0      0.093      0.093      0.100      0.000     -0.000     -0.000 3/8400
  570.0      0.091      0.091      0.098      0.000     -0.000     -0.000 3/8400
  580.0      0.090      0.090      0.097      0.000     -0.000     -0.000 3/8400
  590.0      0.088      0.088      0.095      0.000     -0.000     -0.000 3/8400
  600.0      0.087      0.087      0.093      0.000     -0.000     -0.000 3/8400
  610.0      0.085      0.085      0.092      0.000     -0.000     -0.000 3/8400
  620.0      0.084      0.084      0.090      0.000     -0.000     -0.000 3/8400
  630.0      0.083      0.083      0.089      0.000     -0.000     -0.000 3/8400
  640.0      0.081      0.081      0.088      0.000     -0.000     -0.000 3/8400
  650.0      0.080      0.080      0.086      0.000     -0.000     -0.000 3/8400
  660.0      0.079      0.079      0.085      0.000     -0.000     -0.000 3/8400
  670.0      0.078      0.078      0.084      0.000     -0.000     -0.000 3/8400
  680.0      0.077      0.077      0.083      0.000     -0.000     -0.000 3/8400
  690.0      0.076      0.076      0.081      0.000     -0.000     -0.000 3/8400
  700.0      0.074      0.074      0.080      0.000     -0.000     -0.000 3/8400
  710.0      0.073      0.073      0.079      0.000     -0.000     -0.000 3/8400
  720.0      0.072      0.072      0.078      0.000     -0.000     -0.000 3/8400
  730.0      0.071      0.071      0.077      0.000     -0.000     -0.000 3/8400
  740.0      0.070      0.070      0.076      0.000     -0.000     -0.000 3/8400
  750.0      0.070      0.070      0.075      0.000     -0.000     -0.000 3/8400
  760.0      0.069      0.069      0.074      0.000     -0.000     -0.000 3/8400
  770.0      0.068      0.068      0.073      0.000     -0.000     -0.000 3/8400
  780.0      0.067      0.067      0.072      0.000     -0.000     -0.000 3/8400
  790.0      0.066      0.066      0.071      0.000     -0.000     -0.000 3/8400
  800.0      0.065      0.065      0.070      0.000     -0.000     -0.000 3/8400
  810.0      0.064      0.064      0.069      0.000     -0.000     -0.000 3/8400
  820.0      0.064      0.064      0.068      0.000     -0.000     -0.000 3/8400
  830.0      0.063      0.063      0.068      0.000     -0.000     -0.000 3/8400
  840.0      0.062      0.062      0.067      0.000     -0.000     -0.000 3/8400
  850.0      0.061      0.061      0.066      0.000     -0.000     -0.000 3/8400
  860.0      0.061      0.061      0.065      0.000     -0.000     -0.000 3/8400
  870.0      0.060      0.060      0.065      0.000     -0.000     -0.000 3/8400
  880.0      0.059      0.059      0.064      0.000     -0.000     -0.000 3/8400
  890.0      0.059      0.059      0.063      0.000     -0.000     -0.000 3/8400
  900.0      0.058      0.058      0.062      0.000     -0.000     -0.000 3/8400
  910.0      0.057      0.057      0.062      0.000     -0.000     -0.000 3/8400
  920.0      0.057      0.057      0.061      0.000     -0.000     -0.000 3/8400
  930.0      0.056      0.056      0.060      0.000     -0.000     -0.000 3/8400
  940.0      0.056      0.056      0.060      0.000     -0.000     -0.000 3/8400
  950.0      0.055      0.055      0.059      0.000     -0.000     -0.000 3/8400
  960.0      0.054      0.054      0.059      0.000     -0.000     -0.000 3/8400
  970.0      0.054      0.054      0.058      0.000     -0.000     -0.000 3/8400
  980.0      0.053      0.053      0.057      0.000     -0.000     -0.000 3/8400
  990.0      0.053      0.053      0.057      0.000     -0.000     -0.000 3/8400
 1000.0      0.052      0.052      0.056      0.000     -0.000     -0.000 3/8400

Thermal conductivity related properties were written into 
"kappa-m554.hdf5".
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 10:22:46]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

