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

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

-------------------------[time 2026-01-08 17:47:48]-------------------------
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:
    [-0.5  0.5  0.5]
    [ 0.5 -0.5  0.5]
    [ 0.5  0.5 -0.5]
Spacegroup: I4/mcm (140)
Number of symmetry operations in supercell: 128
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   -4.103592924999999    4.103592924999999    7.419846740000000
  b    4.103592924999999   -4.103592924999999    7.419846740000000
  c    4.103592924999999    4.103592924999999   -7.419846740000000
Atomic positions (fractional):
   *1 Rb  0.25000000000000  0.25000000000000  0.00000000000000  85.468
    2 Rb  0.75000000000000  0.75000000000000  0.00000000000000  85.468
   *3 Sn  0.16800355217870  0.33199644782130  0.50000000000000 118.710
    4 Sn  0.33199644782130  0.83199644782130  0.16399289564260 118.710
    5 Sn  0.83199644782130  0.66800355217870  0.50000000000000 118.710
    6 Sn  0.66800355217870  0.16800355217870  0.83600710435740 118.710
   *7 Br  0.47354753945124  0.29008362008363  0.50000000000000  79.904
    8 Br  0.29008362008363  0.79008362008363  0.81653608063239  79.904
    9 Br  0.52645246054876  0.70991637991637  0.50000000000000  79.904
   10 Br  0.79008362008363  0.97354753945124  0.50000000000000  79.904
   11 Br  0.97354753945124  0.47354753945124  0.18346391936761  79.904
   12 Br  0.70991637991637  0.20991637991637  0.18346391936761  79.904
   13 Br  0.20991637991637  0.02645246054876  0.50000000000000  79.904
  *14 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904
   15 Br  0.02645246054876  0.52645246054876  0.81653608063239  79.904
   16 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    8.207185849999998    0.000000000000000    0.000000000000000
  b    0.000000000000000    8.207185849999998    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.839693479999999
Atomic positions (fractional):
   *1 Rb  0.00000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    2 Rb  0.50000000000000  0.50000000000000  0.25000000000000  85.468 > 2
    3 Rb  0.50000000000000  0.50000000000000  0.75000000000000  85.468 > 1
    4 Rb  0.00000000000000  0.00000000000000  0.75000000000000  85.468 > 2
   *5 Sn  0.33199644782130  0.16800355217870  0.00000000000000 118.710 > 3
    6 Sn  0.83199644782130  0.33199644782130  0.00000000000000 118.710 > 4
    7 Sn  0.66800355217870  0.83199644782130  0.00000000000000 118.710 > 5
    8 Sn  0.16800355217870  0.66800355217870  0.00000000000000 118.710 > 6
    9 Sn  0.83199644782130  0.66800355217870  0.50000000000000 118.710 > 3
   10 Sn  0.33199644782130  0.83199644782130  0.50000000000000 118.710 > 4
   11 Sn  0.16800355217870  0.33199644782130  0.50000000000000 118.710 > 5
   12 Sn  0.66800355217870  0.16800355217870  0.50000000000000 118.710 > 6
  *13 Br  0.15826804031619  0.34173195968381  0.13181557976744  79.904 > 7
   14 Br  0.65826804031619  0.15826804031619  0.13181557976744  79.904 > 8
   15 Br  0.34173195968381  0.15826804031619  0.36818442023256  79.904 > 9
   16 Br  0.84173195968381  0.65826804031619  0.13181557976744  79.904 > 10
   17 Br  0.34173195968381  0.84173195968381  0.13181557976744  79.904 > 11
   18 Br  0.84173195968381  0.34173195968381  0.36818442023256  79.904 > 12
   19 Br  0.65826804031619  0.84173195968381  0.36818442023256  79.904 > 13
  *20 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 14
   21 Br  0.15826804031619  0.65826804031619  0.36818442023256  79.904 > 15
   22 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 16
   23 Br  0.65826804031619  0.84173195968381  0.63181557976744  79.904 > 7
   24 Br  0.15826804031619  0.65826804031619  0.63181557976744  79.904 > 8
   25 Br  0.84173195968381  0.65826804031619  0.86818442023256  79.904 > 9
   26 Br  0.34173195968381  0.15826804031619  0.63181557976744  79.904 > 10
   27 Br  0.84173195968381  0.34173195968381  0.63181557976744  79.904 > 11
   28 Br  0.34173195968381  0.84173195968381  0.86818442023256  79.904 > 12
   29 Br  0.15826804031619  0.34173195968381  0.86818442023256  79.904 > 13
   30 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 14
   31 Br  0.65826804031619  0.15826804031619  0.86818442023256  79.904 > 15
   32 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 16
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   16.414371699999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   16.414371699999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.839693479999999
Atomic positions (fractional):
   *1 Rb  0.00000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    2 Rb  0.50000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    3 Rb  0.00000000000000  0.50000000000000  0.25000000000000  85.468 > 1
    4 Rb  0.50000000000000  0.50000000000000  0.25000000000000  85.468 > 1
    5 Rb  0.25000000000000  0.25000000000000  0.25000000000000  85.468 > 2
    6 Rb  0.75000000000000  0.25000000000000  0.25000000000000  85.468 > 2
    7 Rb  0.25000000000000  0.75000000000000  0.25000000000000  85.468 > 2
    8 Rb  0.75000000000000  0.75000000000000  0.25000000000000  85.468 > 2
    9 Rb  0.25000000000000  0.25000000000000  0.75000000000000  85.468 > 1
   10 Rb  0.75000000000000  0.25000000000000  0.75000000000000  85.468 > 1
   11 Rb  0.25000000000000  0.75000000000000  0.75000000000000  85.468 > 1
   12 Rb  0.75000000000000  0.75000000000000  0.75000000000000  85.468 > 1
   13 Rb  0.00000000000000  0.00000000000000  0.75000000000000  85.468 > 2
   14 Rb  0.50000000000000  0.00000000000000  0.75000000000000  85.468 > 2
   15 Rb  0.00000000000000  0.50000000000000  0.75000000000000  85.468 > 2
   16 Rb  0.50000000000000  0.50000000000000  0.75000000000000  85.468 > 2
  *17 Sn  0.16599822391065  0.08400177608935  0.00000000000000 118.710 > 3
   18 Sn  0.66599822391065  0.08400177608935  0.00000000000000 118.710 > 3
   19 Sn  0.16599822391065  0.58400177608935  0.00000000000000 118.710 > 3
   20 Sn  0.66599822391065  0.58400177608935  0.00000000000000 118.710 > 3
   21 Sn  0.41599822391065  0.16599822391065  0.00000000000000 118.710 > 4
   22 Sn  0.91599822391065  0.16599822391065  0.00000000000000 118.710 > 4
   23 Sn  0.41599822391065  0.66599822391065  0.00000000000000 118.710 > 4
   24 Sn  0.91599822391065  0.66599822391065  0.00000000000000 118.710 > 4
   25 Sn  0.33400177608935  0.41599822391065  0.00000000000000 118.710 > 5
   26 Sn  0.83400177608935  0.41599822391065  0.00000000000000 118.710 > 5
   27 Sn  0.33400177608935  0.91599822391065  0.00000000000000 118.710 > 5
   28 Sn  0.83400177608935  0.91599822391065  0.00000000000000 118.710 > 5
   29 Sn  0.08400177608935  0.33400177608935  0.00000000000000 118.710 > 6
   30 Sn  0.58400177608935  0.33400177608935  0.00000000000000 118.710 > 6
   31 Sn  0.08400177608935  0.83400177608935  0.00000000000000 118.710 > 6
   32 Sn  0.58400177608935  0.83400177608935  0.00000000000000 118.710 > 6
   33 Sn  0.41599822391065  0.33400177608935  0.50000000000000 118.710 > 3
   34 Sn  0.91599822391065  0.33400177608935  0.50000000000000 118.710 > 3
   35 Sn  0.41599822391065  0.83400177608935  0.50000000000000 118.710 > 3
   36 Sn  0.91599822391065  0.83400177608935  0.50000000000000 118.710 > 3
   37 Sn  0.16599822391065  0.41599822391065  0.50000000000000 118.710 > 4
   38 Sn  0.66599822391065  0.41599822391065  0.50000000000000 118.710 > 4
   39 Sn  0.16599822391065  0.91599822391065  0.50000000000000 118.710 > 4
   40 Sn  0.66599822391065  0.91599822391065  0.50000000000000 118.710 > 4
   41 Sn  0.08400177608935  0.16599822391065  0.50000000000000 118.710 > 5
   42 Sn  0.58400177608935  0.16599822391065  0.50000000000000 118.710 > 5
   43 Sn  0.08400177608935  0.66599822391065  0.50000000000000 118.710 > 5
   44 Sn  0.58400177608935  0.66599822391065  0.50000000000000 118.710 > 5
   45 Sn  0.33400177608935  0.08400177608935  0.50000000000000 118.710 > 6
   46 Sn  0.83400177608935  0.08400177608935  0.50000000000000 118.710 > 6
   47 Sn  0.33400177608935  0.58400177608935  0.50000000000000 118.710 > 6
   48 Sn  0.83400177608935  0.58400177608935  0.50000000000000 118.710 > 6
  *49 Br  0.07913402015810  0.17086597984190  0.13181557976744  79.904 > 7
   50 Br  0.57913402015810  0.17086597984190  0.13181557976744  79.904 > 7
   51 Br  0.07913402015810  0.67086597984190  0.13181557976744  79.904 > 7
   52 Br  0.57913402015810  0.67086597984190  0.13181557976744  79.904 > 7
   53 Br  0.32913402015810  0.07913402015810  0.13181557976744  79.904 > 8
   54 Br  0.82913402015810  0.07913402015810  0.13181557976744  79.904 > 8
   55 Br  0.32913402015810  0.57913402015810  0.13181557976744  79.904 > 8
   56 Br  0.82913402015810  0.57913402015810  0.13181557976744  79.904 > 8
   57 Br  0.17086597984190  0.07913402015810  0.36818442023256  79.904 > 9
   58 Br  0.67086597984190  0.07913402015810  0.36818442023256  79.904 > 9
   59 Br  0.17086597984190  0.57913402015810  0.36818442023256  79.904 > 9
   60 Br  0.67086597984190  0.57913402015810  0.36818442023256  79.904 > 9
   61 Br  0.42086597984190  0.32913402015810  0.13181557976744  79.904 > 10
   62 Br  0.92086597984190  0.32913402015810  0.13181557976744  79.904 > 10
   63 Br  0.42086597984190  0.82913402015810  0.13181557976744  79.904 > 10
   64 Br  0.92086597984190  0.82913402015810  0.13181557976744  79.904 > 10
   65 Br  0.17086597984190  0.42086597984190  0.13181557976744  79.904 > 11
   66 Br  0.67086597984190  0.42086597984190  0.13181557976744  79.904 > 11
   67 Br  0.17086597984190  0.92086597984190  0.13181557976744  79.904 > 11
   68 Br  0.67086597984190  0.92086597984190  0.13181557976744  79.904 > 11
   69 Br  0.42086597984190  0.17086597984190  0.36818442023256  79.904 > 12
   70 Br  0.92086597984190  0.17086597984190  0.36818442023256  79.904 > 12
   71 Br  0.42086597984190  0.67086597984190  0.36818442023256  79.904 > 12
   72 Br  0.92086597984190  0.67086597984190  0.36818442023256  79.904 > 12
   73 Br  0.32913402015810  0.42086597984190  0.36818442023256  79.904 > 13
   74 Br  0.82913402015810  0.42086597984190  0.36818442023256  79.904 > 13
   75 Br  0.32913402015810  0.92086597984190  0.36818442023256  79.904 > 13
   76 Br  0.82913402015810  0.92086597984190  0.36818442023256  79.904 > 13
  *77 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 14
   78 Br  0.50000000000000  0.00000000000000  0.00000000000000  79.904 > 14
   79 Br  0.00000000000000  0.50000000000000  0.00000000000000  79.904 > 14
   80 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 14
   81 Br  0.07913402015810  0.32913402015810  0.36818442023256  79.904 > 15
   82 Br  0.57913402015810  0.32913402015810  0.36818442023256  79.904 > 15
   83 Br  0.07913402015810  0.82913402015810  0.36818442023256  79.904 > 15
   84 Br  0.57913402015810  0.82913402015810  0.36818442023256  79.904 > 15
   85 Br  0.25000000000000  0.25000000000000  0.00000000000000  79.904 > 16
   86 Br  0.75000000000000  0.25000000000000  0.00000000000000  79.904 > 16
   87 Br  0.25000000000000  0.75000000000000  0.00000000000000  79.904 > 16
   88 Br  0.75000000000000  0.75000000000000  0.00000000000000  79.904 > 16
   89 Br  0.32913402015810  0.42086597984190  0.63181557976744  79.904 > 7
   90 Br  0.82913402015810  0.42086597984190  0.63181557976744  79.904 > 7
   91 Br  0.32913402015810  0.92086597984190  0.63181557976744  79.904 > 7
   92 Br  0.82913402015810  0.92086597984190  0.63181557976744  79.904 > 7
   93 Br  0.07913402015810  0.32913402015810  0.63181557976744  79.904 > 8
   94 Br  0.57913402015810  0.32913402015810  0.63181557976744  79.904 > 8
   95 Br  0.07913402015810  0.82913402015810  0.63181557976744  79.904 > 8
   96 Br  0.57913402015810  0.82913402015810  0.63181557976744  79.904 > 8
   97 Br  0.42086597984190  0.32913402015810  0.86818442023256  79.904 > 9
   98 Br  0.92086597984190  0.32913402015810  0.86818442023256  79.904 > 9
   99 Br  0.42086597984190  0.82913402015810  0.86818442023256  79.904 > 9
  100 Br  0.92086597984190  0.82913402015810  0.86818442023256  79.904 > 9
  101 Br  0.17086597984190  0.07913402015810  0.63181557976744  79.904 > 10
  102 Br  0.67086597984190  0.07913402015810  0.63181557976744  79.904 > 10
  103 Br  0.17086597984190  0.57913402015810  0.63181557976744  79.904 > 10
  104 Br  0.67086597984190  0.57913402015810  0.63181557976744  79.904 > 10
  105 Br  0.42086597984190  0.17086597984190  0.63181557976744  79.904 > 11
  106 Br  0.92086597984190  0.17086597984190  0.63181557976744  79.904 > 11
  107 Br  0.42086597984190  0.67086597984190  0.63181557976744  79.904 > 11
  108 Br  0.92086597984190  0.67086597984190  0.63181557976744  79.904 > 11
  109 Br  0.17086597984190  0.42086597984190  0.86818442023256  79.904 > 12
  110 Br  0.67086597984190  0.42086597984190  0.86818442023256  79.904 > 12
  111 Br  0.17086597984190  0.92086597984190  0.86818442023256  79.904 > 12
  112 Br  0.67086597984190  0.92086597984190  0.86818442023256  79.904 > 12
  113 Br  0.07913402015810  0.17086597984190  0.86818442023256  79.904 > 13
  114 Br  0.57913402015810  0.17086597984190  0.86818442023256  79.904 > 13
  115 Br  0.07913402015810  0.67086597984190  0.86818442023256  79.904 > 13
  116 Br  0.57913402015810  0.67086597984190  0.86818442023256  79.904 > 13
  117 Br  0.25000000000000  0.25000000000000  0.50000000000000  79.904 > 14
  118 Br  0.75000000000000  0.25000000000000  0.50000000000000  79.904 > 14
  119 Br  0.25000000000000  0.75000000000000  0.50000000000000  79.904 > 14
  120 Br  0.75000000000000  0.75000000000000  0.50000000000000  79.904 > 14
  121 Br  0.32913402015810  0.07913402015810  0.86818442023256  79.904 > 15
  122 Br  0.82913402015810  0.07913402015810  0.86818442023256  79.904 > 15
  123 Br  0.32913402015810  0.57913402015810  0.86818442023256  79.904 > 15
  124 Br  0.82913402015810  0.57913402015810  0.86818442023256  79.904 > 15
  125 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 16
  126 Br  0.50000000000000  0.00000000000000  0.50000000000000  79.904 > 16
  127 Br  0.00000000000000  0.50000000000000  0.50000000000000  79.904 > 16
  128 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 16
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            5.9284900    0.0000000    0.0000000
            0.0000000    5.9284900    0.0000000
            0.0000000    0.0000000    4.0652419
-------------------------- Born effective charges --------------------------
    1 Rb    1.2367491    0.0000000   -0.0000000
            0.0000000    1.2367491    0.0000000
           -0.0000000    0.0000000    1.7464545
    2 Rb    1.2367491    0.0000000   -0.0000000
            0.0000000    1.2367491    0.0000000
           -0.0000000    0.0000000    1.7464545
    3 Sn    4.3725586    0.1008536   -0.0000000
            0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    4 Sn    4.3725586   -0.1008536   -0.0000000
           -0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    5 Sn    4.3725586    0.1008536   -0.0000000
            0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    6 Sn    4.3725586   -0.1008536   -0.0000000
           -0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    7 Br   -1.7462628    0.2622448    0.2123810
            0.2622448   -1.7462628   -0.2123810
            0.0800262   -0.0800262   -1.3168415
    8 Br   -1.7462628   -0.2622448    0.2123810
           -0.2622448   -1.7462628    0.2123810
            0.0800262    0.0800262   -1.3168415
    9 Br   -1.7462628    0.2622448    0.2123810
            0.2622448   -1.7462628   -0.2123810
            0.0800262   -0.0800262   -1.3168415
   10 Br   -1.7462628    0.2622448   -0.2123810
            0.2622448   -1.7462628    0.2123810
           -0.0800262    0.0800262   -1.3168415
   11 Br   -1.7462628   -0.2622448   -0.2123810
           -0.2622448   -1.7462628   -0.2123810
           -0.0800262   -0.0800262   -1.3168415
   12 Br   -1.7462628   -0.2622448    0.2123810
           -0.2622448   -1.7462628    0.2123810
            0.0800262    0.0800262   -1.3168415
   13 Br   -1.7462628    0.2622448   -0.2123810
            0.2622448   -1.7462628    0.2123810
           -0.0800262    0.0800262   -1.3168415
   14 Br   -2.9968150    0.2846536   -0.0000000
           -0.2846536   -2.9968150    0.0000000
           -0.0000000    0.0000000   -0.4476406
   15 Br   -1.7462628   -0.2622448   -0.2123810
           -0.2622448   -1.7462628   -0.2123810
           -0.0800262   -0.0800262   -1.3168415
   16 Br   -2.9968150   -0.2846536   -0.0000000
            0.2846536   -2.9968150    0.0000000
           -0.0000000    0.0000000   -0.4476406
----------------------------------------------------------------------------
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: 384/384
Permutation basis: 11568/11568
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 689
Number of blocks in projector: 689
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 4
--- Eigsh_solver_block: 1 / 4 ---
Block_size: 214
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 4 ---
Block_size: 193
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 4 ---
Block_size: 163
Use standard eigh solver.
--- Eigsh_solver_block: 4 / 4 ---
Block_size: 119
Use standard eigh solver.
Tree of FC basis block matrices:
- (689, 676), data: False
|-- (119, 115), data: True
|-- (163, 160), data: True
|-- (193, 189), data: True
|-- (214, 212), data: True
-----
Solver_atoms: 1 -- 128 / 128
Time (Solver_compr_matrix_reshape): 0.004
Solver_block: 100 / 160
 - Time: 0.652
Solver_block: 160 / 160
 - Time: 0.439
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.105
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 384/384
Permutation basis: 11568/11568
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 689
Number of blocks in projector: 689
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 4
--- Eigsh_solver_block: 1 / 4 ---
Block_size: 214
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 4 ---
Block_size: 193
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 4 ---
Block_size: 163
Use standard eigh solver.
--- Eigsh_solver_block: 4 / 4 ---
Block_size: 119
Use standard eigh solver.
Tree of FC basis block matrices:
- (689, 676), data: False
|-- (119, 115), data: True
|-- (163, 160), data: True
|-- (193, 189), data: True
|-- (214, 212), 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 17:47:58]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 17:47:58]-------------------------
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:
  [-0.5  0.5  0.5]
  [ 0.5 -0.5  0.5]
  [ 0.5  0.5 -0.5]
Spacegroup: I4/mcm (140)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   -4.103592924999999    4.103592924999999    7.419846740000000
  b    4.103592924999999   -4.103592924999999    7.419846740000000
  c    4.103592924999999    4.103592924999999   -7.419846740000000
Atomic positions (fractional):
    1 Rb  0.25000000000000  0.25000000000000  0.00000000000000  85.468
    2 Rb  0.75000000000000  0.75000000000000  0.00000000000000  85.468
    3 Sn  0.16800355217870  0.33199644782130  0.50000000000000 118.710
    4 Sn  0.33199644782130  0.83199644782130  0.16399289564260 118.710
    5 Sn  0.83199644782130  0.66800355217870  0.50000000000000 118.710
    6 Sn  0.66800355217870  0.16800355217870  0.83600710435740 118.710
    7 Br  0.47354753945124  0.29008362008363  0.50000000000000  79.904
    8 Br  0.29008362008363  0.79008362008363  0.81653608063239  79.904
    9 Br  0.52645246054876  0.70991637991637  0.50000000000000  79.904
   10 Br  0.79008362008363  0.97354753945124  0.50000000000000  79.904
   11 Br  0.97354753945124  0.47354753945124  0.18346391936761  79.904
   12 Br  0.70991637991637  0.20991637991637  0.18346391936761  79.904
   13 Br  0.20991637991637  0.02645246054876  0.50000000000000  79.904
   14 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904
   15 Br  0.02645246054876  0.52645246054876  0.81653608063239  79.904
   16 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.414371699999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   16.414371699999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.839693479999999
Atomic positions (fractional):
    1 Rb  0.00000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    2 Rb  0.50000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    3 Rb  0.00000000000000  0.50000000000000  0.25000000000000  85.468 > 1
    4 Rb  0.50000000000000  0.50000000000000  0.25000000000000  85.468 > 1
    5 Rb  0.25000000000000  0.25000000000000  0.25000000000000  85.468 > 5
    6 Rb  0.75000000000000  0.25000000000000  0.25000000000000  85.468 > 5
    7 Rb  0.25000000000000  0.75000000000000  0.25000000000000  85.468 > 5
    8 Rb  0.75000000000000  0.75000000000000  0.25000000000000  85.468 > 5
    9 Rb  0.25000000000000  0.25000000000000  0.75000000000000  85.468 > 1
   10 Rb  0.75000000000000  0.25000000000000  0.75000000000000  85.468 > 1
   11 Rb  0.25000000000000  0.75000000000000  0.75000000000000  85.468 > 1
   12 Rb  0.75000000000000  0.75000000000000  0.75000000000000  85.468 > 1
   13 Rb  0.00000000000000  0.00000000000000  0.75000000000000  85.468 > 5
   14 Rb  0.50000000000000  0.00000000000000  0.75000000000000  85.468 > 5
   15 Rb  0.00000000000000  0.50000000000000  0.75000000000000  85.468 > 5
   16 Rb  0.50000000000000  0.50000000000000  0.75000000000000  85.468 > 5
   17 Sn  0.16599822391065  0.08400177608935  0.00000000000000 118.710 > 17
   18 Sn  0.66599822391065  0.08400177608935  0.00000000000000 118.710 > 17
   19 Sn  0.16599822391065  0.58400177608935  0.00000000000000 118.710 > 17
   20 Sn  0.66599822391065  0.58400177608935  0.00000000000000 118.710 > 17
   21 Sn  0.41599822391065  0.16599822391065  0.00000000000000 118.710 > 21
   22 Sn  0.91599822391065  0.16599822391065  0.00000000000000 118.710 > 21
   23 Sn  0.41599822391065  0.66599822391065  0.00000000000000 118.710 > 21
   24 Sn  0.91599822391065  0.66599822391065  0.00000000000000 118.710 > 21
   25 Sn  0.33400177608935  0.41599822391065  0.00000000000000 118.710 > 25
   26 Sn  0.83400177608935  0.41599822391065  0.00000000000000 118.710 > 25
   27 Sn  0.33400177608935  0.91599822391065  0.00000000000000 118.710 > 25
   28 Sn  0.83400177608935  0.91599822391065  0.00000000000000 118.710 > 25
   29 Sn  0.08400177608935  0.33400177608935  0.00000000000000 118.710 > 29
   30 Sn  0.58400177608935  0.33400177608935  0.00000000000000 118.710 > 29
   31 Sn  0.08400177608935  0.83400177608935  0.00000000000000 118.710 > 29
   32 Sn  0.58400177608935  0.83400177608935  0.00000000000000 118.710 > 29
   33 Sn  0.41599822391065  0.33400177608935  0.50000000000000 118.710 > 17
   34 Sn  0.91599822391065  0.33400177608935  0.50000000000000 118.710 > 17
   35 Sn  0.41599822391065  0.83400177608935  0.50000000000000 118.710 > 17
   36 Sn  0.91599822391065  0.83400177608935  0.50000000000000 118.710 > 17
   37 Sn  0.16599822391065  0.41599822391065  0.50000000000000 118.710 > 21
   38 Sn  0.66599822391065  0.41599822391065  0.50000000000000 118.710 > 21
   39 Sn  0.16599822391065  0.91599822391065  0.50000000000000 118.710 > 21
   40 Sn  0.66599822391065  0.91599822391065  0.50000000000000 118.710 > 21
   41 Sn  0.08400177608935  0.16599822391065  0.50000000000000 118.710 > 25
   42 Sn  0.58400177608935  0.16599822391065  0.50000000000000 118.710 > 25
   43 Sn  0.08400177608935  0.66599822391065  0.50000000000000 118.710 > 25
   44 Sn  0.58400177608935  0.66599822391065  0.50000000000000 118.710 > 25
   45 Sn  0.33400177608935  0.08400177608935  0.50000000000000 118.710 > 29
   46 Sn  0.83400177608935  0.08400177608935  0.50000000000000 118.710 > 29
   47 Sn  0.33400177608935  0.58400177608935  0.50000000000000 118.710 > 29
   48 Sn  0.83400177608935  0.58400177608935  0.50000000000000 118.710 > 29
   49 Br  0.07913402015810  0.17086597984190  0.13181557976744  79.904 > 49
   50 Br  0.57913402015810  0.17086597984190  0.13181557976744  79.904 > 49
   51 Br  0.07913402015810  0.67086597984190  0.13181557976744  79.904 > 49
   52 Br  0.57913402015810  0.67086597984190  0.13181557976744  79.904 > 49
   53 Br  0.32913402015810  0.07913402015810  0.13181557976744  79.904 > 53
   54 Br  0.82913402015810  0.07913402015810  0.13181557976744  79.904 > 53
   55 Br  0.32913402015810  0.57913402015810  0.13181557976744  79.904 > 53
   56 Br  0.82913402015810  0.57913402015810  0.13181557976744  79.904 > 53
   57 Br  0.17086597984190  0.07913402015810  0.36818442023256  79.904 > 57
   58 Br  0.67086597984190  0.07913402015810  0.36818442023256  79.904 > 57
   59 Br  0.17086597984190  0.57913402015810  0.36818442023256  79.904 > 57
   60 Br  0.67086597984190  0.57913402015810  0.36818442023256  79.904 > 57
   61 Br  0.42086597984190  0.32913402015810  0.13181557976744  79.904 > 61
   62 Br  0.92086597984190  0.32913402015810  0.13181557976744  79.904 > 61
   63 Br  0.42086597984190  0.82913402015810  0.13181557976744  79.904 > 61
   64 Br  0.92086597984190  0.82913402015810  0.13181557976744  79.904 > 61
   65 Br  0.17086597984190  0.42086597984190  0.13181557976744  79.904 > 65
   66 Br  0.67086597984190  0.42086597984190  0.13181557976744  79.904 > 65
   67 Br  0.17086597984190  0.92086597984190  0.13181557976744  79.904 > 65
   68 Br  0.67086597984190  0.92086597984190  0.13181557976744  79.904 > 65
   69 Br  0.42086597984190  0.17086597984190  0.36818442023256  79.904 > 69
   70 Br  0.92086597984190  0.17086597984190  0.36818442023256  79.904 > 69
   71 Br  0.42086597984190  0.67086597984190  0.36818442023256  79.904 > 69
   72 Br  0.92086597984190  0.67086597984190  0.36818442023256  79.904 > 69
   73 Br  0.32913402015810  0.42086597984190  0.36818442023256  79.904 > 73
   74 Br  0.82913402015810  0.42086597984190  0.36818442023256  79.904 > 73
   75 Br  0.32913402015810  0.92086597984190  0.36818442023256  79.904 > 73
   76 Br  0.82913402015810  0.92086597984190  0.36818442023256  79.904 > 73
   77 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 77
   78 Br  0.50000000000000  0.00000000000000  0.00000000000000  79.904 > 77
   79 Br  0.00000000000000  0.50000000000000  0.00000000000000  79.904 > 77
   80 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 77
   81 Br  0.07913402015810  0.32913402015810  0.36818442023256  79.904 > 81
   82 Br  0.57913402015810  0.32913402015810  0.36818442023256  79.904 > 81
   83 Br  0.07913402015810  0.82913402015810  0.36818442023256  79.904 > 81
   84 Br  0.57913402015810  0.82913402015810  0.36818442023256  79.904 > 81
   85 Br  0.25000000000000  0.25000000000000  0.00000000000000  79.904 > 85
   86 Br  0.75000000000000  0.25000000000000  0.00000000000000  79.904 > 85
   87 Br  0.25000000000000  0.75000000000000  0.00000000000000  79.904 > 85
   88 Br  0.75000000000000  0.75000000000000  0.00000000000000  79.904 > 85
   89 Br  0.32913402015810  0.42086597984190  0.63181557976744  79.904 > 49
   90 Br  0.82913402015810  0.42086597984190  0.63181557976744  79.904 > 49
   91 Br  0.32913402015810  0.92086597984190  0.63181557976744  79.904 > 49
   92 Br  0.82913402015810  0.92086597984190  0.63181557976744  79.904 > 49
   93 Br  0.07913402015810  0.32913402015810  0.63181557976744  79.904 > 53
   94 Br  0.57913402015810  0.32913402015810  0.63181557976744  79.904 > 53
   95 Br  0.07913402015810  0.82913402015810  0.63181557976744  79.904 > 53
   96 Br  0.57913402015810  0.82913402015810  0.63181557976744  79.904 > 53
   97 Br  0.42086597984190  0.32913402015810  0.86818442023256  79.904 > 57
   98 Br  0.92086597984190  0.32913402015810  0.86818442023256  79.904 > 57
   99 Br  0.42086597984190  0.82913402015810  0.86818442023256  79.904 > 57
  100 Br  0.92086597984190  0.82913402015810  0.86818442023256  79.904 > 57
  101 Br  0.17086597984190  0.07913402015810  0.63181557976744  79.904 > 61
  102 Br  0.67086597984190  0.07913402015810  0.63181557976744  79.904 > 61
  103 Br  0.17086597984190  0.57913402015810  0.63181557976744  79.904 > 61
  104 Br  0.67086597984190  0.57913402015810  0.63181557976744  79.904 > 61
  105 Br  0.42086597984190  0.17086597984190  0.63181557976744  79.904 > 65
  106 Br  0.92086597984190  0.17086597984190  0.63181557976744  79.904 > 65
  107 Br  0.42086597984190  0.67086597984190  0.63181557976744  79.904 > 65
  108 Br  0.92086597984190  0.67086597984190  0.63181557976744  79.904 > 65
  109 Br  0.17086597984190  0.42086597984190  0.86818442023256  79.904 > 69
  110 Br  0.67086597984190  0.42086597984190  0.86818442023256  79.904 > 69
  111 Br  0.17086597984190  0.92086597984190  0.86818442023256  79.904 > 69
  112 Br  0.67086597984190  0.92086597984190  0.86818442023256  79.904 > 69
  113 Br  0.07913402015810  0.17086597984190  0.86818442023256  79.904 > 73
  114 Br  0.57913402015810  0.17086597984190  0.86818442023256  79.904 > 73
  115 Br  0.07913402015810  0.67086597984190  0.86818442023256  79.904 > 73
  116 Br  0.57913402015810  0.67086597984190  0.86818442023256  79.904 > 73
  117 Br  0.25000000000000  0.25000000000000  0.50000000000000  79.904 > 77
  118 Br  0.75000000000000  0.25000000000000  0.50000000000000  79.904 > 77
  119 Br  0.25000000000000  0.75000000000000  0.50000000000000  79.904 > 77
  120 Br  0.75000000000000  0.75000000000000  0.50000000000000  79.904 > 77
  121 Br  0.32913402015810  0.07913402015810  0.86818442023256  79.904 > 81
  122 Br  0.82913402015810  0.07913402015810  0.86818442023256  79.904 > 81
  123 Br  0.32913402015810  0.57913402015810  0.86818442023256  79.904 > 81
  124 Br  0.82913402015810  0.57913402015810  0.86818442023256  79.904 > 81
  125 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 85
  126 Br  0.50000000000000  0.00000000000000  0.50000000000000  79.904 > 85
  127 Br  0.00000000000000  0.50000000000000  0.50000000000000  79.904 > 85
  128 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 85
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            5.9284900    0.0000000    0.0000000
            0.0000000    5.9284900    0.0000000
            0.0000000    0.0000000    4.0652419
-------------------------- Born effective charges --------------------------
    1 Rb    1.2367491    0.0000000   -0.0000000
            0.0000000    1.2367491    0.0000000
           -0.0000000    0.0000000    1.7464545
    2 Rb    1.2367491    0.0000000   -0.0000000
            0.0000000    1.2367491    0.0000000
           -0.0000000    0.0000000    1.7464545
    3 Sn    4.3725586    0.1008536   -0.0000000
            0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    4 Sn    4.3725586   -0.1008536   -0.0000000
           -0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    5 Sn    4.3725586    0.1008536   -0.0000000
            0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    6 Sn    4.3725586   -0.1008536   -0.0000000
           -0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    7 Br   -1.7462628    0.2622448    0.2123810
            0.2622448   -1.7462628   -0.2123810
            0.0800262   -0.0800262   -1.3168415
    8 Br   -1.7462628   -0.2622448    0.2123810
           -0.2622448   -1.7462628    0.2123810
            0.0800262    0.0800262   -1.3168415
    9 Br   -1.7462628    0.2622448    0.2123810
            0.2622448   -1.7462628   -0.2123810
            0.0800262   -0.0800262   -1.3168415
   10 Br   -1.7462628    0.2622448   -0.2123810
            0.2622448   -1.7462628    0.2123810
           -0.0800262    0.0800262   -1.3168415
   11 Br   -1.7462628   -0.2622448   -0.2123810
           -0.2622448   -1.7462628   -0.2123810
           -0.0800262   -0.0800262   -1.3168415
   12 Br   -1.7462628   -0.2622448    0.2123810
           -0.2622448   -1.7462628    0.2123810
            0.0800262    0.0800262   -1.3168415
   13 Br   -1.7462628    0.2622448   -0.2123810
            0.2622448   -1.7462628    0.2123810
           -0.0800262    0.0800262   -1.3168415
   14 Br   -2.9968150    0.2846536   -0.0000000
           -0.2846536   -2.9968150    0.0000000
           -0.0000000    0.0000000   -0.4476406
   15 Br   -1.7462628   -0.2622448   -0.2123810
           -0.2622448   -1.7462628   -0.2123810
           -0.0800262   -0.0800262   -1.3168415
   16 Br   -2.9968150   -0.2846536   -0.0000000
            0.2846536   -2.9968150    0.0000000
           -0.0000000    0.0000000   -0.4476406
----------------------------------------------------------------------------
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[ 17, 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[ 77, 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.00000005 (yyy) -0.00000005 (yyy) -0.00000005 (yyy)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (xy) 0.00000000 (xy) 
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 17:48:09]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 17:48:10]-------------------------
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:
  [-0.5  0.5  0.5]
  [ 0.5 -0.5  0.5]
  [ 0.5  0.5 -0.5]
Spacegroup: I4/mcm (140)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   -4.103592924999999    4.103592924999999    7.419846740000000
  b    4.103592924999999   -4.103592924999999    7.419846740000000
  c    4.103592924999999    4.103592924999999   -7.419846740000000
Atomic positions (fractional):
    1 Rb  0.25000000000000  0.25000000000000  0.00000000000000  85.468
    2 Rb  0.75000000000000  0.75000000000000  0.00000000000000  85.468
    3 Sn  0.16800355217870  0.33199644782130  0.50000000000000 118.710
    4 Sn  0.33199644782130  0.83199644782130  0.16399289564260 118.710
    5 Sn  0.83199644782130  0.66800355217870  0.50000000000000 118.710
    6 Sn  0.66800355217870  0.16800355217870  0.83600710435740 118.710
    7 Br  0.47354753945124  0.29008362008363  0.50000000000000  79.904
    8 Br  0.29008362008363  0.79008362008363  0.81653608063239  79.904
    9 Br  0.52645246054876  0.70991637991637  0.50000000000000  79.904
   10 Br  0.79008362008363  0.97354753945124  0.50000000000000  79.904
   11 Br  0.97354753945124  0.47354753945124  0.18346391936761  79.904
   12 Br  0.70991637991637  0.20991637991637  0.18346391936761  79.904
   13 Br  0.20991637991637  0.02645246054876  0.50000000000000  79.904
   14 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904
   15 Br  0.02645246054876  0.52645246054876  0.81653608063239  79.904
   16 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.414371699999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   16.414371699999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.839693479999999
Atomic positions (fractional):
    1 Rb  0.00000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    2 Rb  0.50000000000000  0.00000000000000  0.25000000000000  85.468 > 1
    3 Rb  0.00000000000000  0.50000000000000  0.25000000000000  85.468 > 1
    4 Rb  0.50000000000000  0.50000000000000  0.25000000000000  85.468 > 1
    5 Rb  0.25000000000000  0.25000000000000  0.25000000000000  85.468 > 5
    6 Rb  0.75000000000000  0.25000000000000  0.25000000000000  85.468 > 5
    7 Rb  0.25000000000000  0.75000000000000  0.25000000000000  85.468 > 5
    8 Rb  0.75000000000000  0.75000000000000  0.25000000000000  85.468 > 5
    9 Rb  0.25000000000000  0.25000000000000  0.75000000000000  85.468 > 1
   10 Rb  0.75000000000000  0.25000000000000  0.75000000000000  85.468 > 1
   11 Rb  0.25000000000000  0.75000000000000  0.75000000000000  85.468 > 1
   12 Rb  0.75000000000000  0.75000000000000  0.75000000000000  85.468 > 1
   13 Rb  0.00000000000000  0.00000000000000  0.75000000000000  85.468 > 5
   14 Rb  0.50000000000000  0.00000000000000  0.75000000000000  85.468 > 5
   15 Rb  0.00000000000000  0.50000000000000  0.75000000000000  85.468 > 5
   16 Rb  0.50000000000000  0.50000000000000  0.75000000000000  85.468 > 5
   17 Sn  0.16599822391065  0.08400177608935  0.00000000000000 118.710 > 17
   18 Sn  0.66599822391065  0.08400177608935  0.00000000000000 118.710 > 17
   19 Sn  0.16599822391065  0.58400177608935  0.00000000000000 118.710 > 17
   20 Sn  0.66599822391065  0.58400177608935  0.00000000000000 118.710 > 17
   21 Sn  0.41599822391065  0.16599822391065  0.00000000000000 118.710 > 21
   22 Sn  0.91599822391065  0.16599822391065  0.00000000000000 118.710 > 21
   23 Sn  0.41599822391065  0.66599822391065  0.00000000000000 118.710 > 21
   24 Sn  0.91599822391065  0.66599822391065  0.00000000000000 118.710 > 21
   25 Sn  0.33400177608935  0.41599822391065  0.00000000000000 118.710 > 25
   26 Sn  0.83400177608935  0.41599822391065  0.00000000000000 118.710 > 25
   27 Sn  0.33400177608935  0.91599822391065  0.00000000000000 118.710 > 25
   28 Sn  0.83400177608935  0.91599822391065  0.00000000000000 118.710 > 25
   29 Sn  0.08400177608935  0.33400177608935  0.00000000000000 118.710 > 29
   30 Sn  0.58400177608935  0.33400177608935  0.00000000000000 118.710 > 29
   31 Sn  0.08400177608935  0.83400177608935  0.00000000000000 118.710 > 29
   32 Sn  0.58400177608935  0.83400177608935  0.00000000000000 118.710 > 29
   33 Sn  0.41599822391065  0.33400177608935  0.50000000000000 118.710 > 17
   34 Sn  0.91599822391065  0.33400177608935  0.50000000000000 118.710 > 17
   35 Sn  0.41599822391065  0.83400177608935  0.50000000000000 118.710 > 17
   36 Sn  0.91599822391065  0.83400177608935  0.50000000000000 118.710 > 17
   37 Sn  0.16599822391065  0.41599822391065  0.50000000000000 118.710 > 21
   38 Sn  0.66599822391065  0.41599822391065  0.50000000000000 118.710 > 21
   39 Sn  0.16599822391065  0.91599822391065  0.50000000000000 118.710 > 21
   40 Sn  0.66599822391065  0.91599822391065  0.50000000000000 118.710 > 21
   41 Sn  0.08400177608935  0.16599822391065  0.50000000000000 118.710 > 25
   42 Sn  0.58400177608935  0.16599822391065  0.50000000000000 118.710 > 25
   43 Sn  0.08400177608935  0.66599822391065  0.50000000000000 118.710 > 25
   44 Sn  0.58400177608935  0.66599822391065  0.50000000000000 118.710 > 25
   45 Sn  0.33400177608935  0.08400177608935  0.50000000000000 118.710 > 29
   46 Sn  0.83400177608935  0.08400177608935  0.50000000000000 118.710 > 29
   47 Sn  0.33400177608935  0.58400177608935  0.50000000000000 118.710 > 29
   48 Sn  0.83400177608935  0.58400177608935  0.50000000000000 118.710 > 29
   49 Br  0.07913402015810  0.17086597984190  0.13181557976744  79.904 > 49
   50 Br  0.57913402015810  0.17086597984190  0.13181557976744  79.904 > 49
   51 Br  0.07913402015810  0.67086597984190  0.13181557976744  79.904 > 49
   52 Br  0.57913402015810  0.67086597984190  0.13181557976744  79.904 > 49
   53 Br  0.32913402015810  0.07913402015810  0.13181557976744  79.904 > 53
   54 Br  0.82913402015810  0.07913402015810  0.13181557976744  79.904 > 53
   55 Br  0.32913402015810  0.57913402015810  0.13181557976744  79.904 > 53
   56 Br  0.82913402015810  0.57913402015810  0.13181557976744  79.904 > 53
   57 Br  0.17086597984190  0.07913402015810  0.36818442023256  79.904 > 57
   58 Br  0.67086597984190  0.07913402015810  0.36818442023256  79.904 > 57
   59 Br  0.17086597984190  0.57913402015810  0.36818442023256  79.904 > 57
   60 Br  0.67086597984190  0.57913402015810  0.36818442023256  79.904 > 57
   61 Br  0.42086597984190  0.32913402015810  0.13181557976744  79.904 > 61
   62 Br  0.92086597984190  0.32913402015810  0.13181557976744  79.904 > 61
   63 Br  0.42086597984190  0.82913402015810  0.13181557976744  79.904 > 61
   64 Br  0.92086597984190  0.82913402015810  0.13181557976744  79.904 > 61
   65 Br  0.17086597984190  0.42086597984190  0.13181557976744  79.904 > 65
   66 Br  0.67086597984190  0.42086597984190  0.13181557976744  79.904 > 65
   67 Br  0.17086597984190  0.92086597984190  0.13181557976744  79.904 > 65
   68 Br  0.67086597984190  0.92086597984190  0.13181557976744  79.904 > 65
   69 Br  0.42086597984190  0.17086597984190  0.36818442023256  79.904 > 69
   70 Br  0.92086597984190  0.17086597984190  0.36818442023256  79.904 > 69
   71 Br  0.42086597984190  0.67086597984190  0.36818442023256  79.904 > 69
   72 Br  0.92086597984190  0.67086597984190  0.36818442023256  79.904 > 69
   73 Br  0.32913402015810  0.42086597984190  0.36818442023256  79.904 > 73
   74 Br  0.82913402015810  0.42086597984190  0.36818442023256  79.904 > 73
   75 Br  0.32913402015810  0.92086597984190  0.36818442023256  79.904 > 73
   76 Br  0.82913402015810  0.92086597984190  0.36818442023256  79.904 > 73
   77 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 77
   78 Br  0.50000000000000  0.00000000000000  0.00000000000000  79.904 > 77
   79 Br  0.00000000000000  0.50000000000000  0.00000000000000  79.904 > 77
   80 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 77
   81 Br  0.07913402015810  0.32913402015810  0.36818442023256  79.904 > 81
   82 Br  0.57913402015810  0.32913402015810  0.36818442023256  79.904 > 81
   83 Br  0.07913402015810  0.82913402015810  0.36818442023256  79.904 > 81
   84 Br  0.57913402015810  0.82913402015810  0.36818442023256  79.904 > 81
   85 Br  0.25000000000000  0.25000000000000  0.00000000000000  79.904 > 85
   86 Br  0.75000000000000  0.25000000000000  0.00000000000000  79.904 > 85
   87 Br  0.25000000000000  0.75000000000000  0.00000000000000  79.904 > 85
   88 Br  0.75000000000000  0.75000000000000  0.00000000000000  79.904 > 85
   89 Br  0.32913402015810  0.42086597984190  0.63181557976744  79.904 > 49
   90 Br  0.82913402015810  0.42086597984190  0.63181557976744  79.904 > 49
   91 Br  0.32913402015810  0.92086597984190  0.63181557976744  79.904 > 49
   92 Br  0.82913402015810  0.92086597984190  0.63181557976744  79.904 > 49
   93 Br  0.07913402015810  0.32913402015810  0.63181557976744  79.904 > 53
   94 Br  0.57913402015810  0.32913402015810  0.63181557976744  79.904 > 53
   95 Br  0.07913402015810  0.82913402015810  0.63181557976744  79.904 > 53
   96 Br  0.57913402015810  0.82913402015810  0.63181557976744  79.904 > 53
   97 Br  0.42086597984190  0.32913402015810  0.86818442023256  79.904 > 57
   98 Br  0.92086597984190  0.32913402015810  0.86818442023256  79.904 > 57
   99 Br  0.42086597984190  0.82913402015810  0.86818442023256  79.904 > 57
  100 Br  0.92086597984190  0.82913402015810  0.86818442023256  79.904 > 57
  101 Br  0.17086597984190  0.07913402015810  0.63181557976744  79.904 > 61
  102 Br  0.67086597984190  0.07913402015810  0.63181557976744  79.904 > 61
  103 Br  0.17086597984190  0.57913402015810  0.63181557976744  79.904 > 61
  104 Br  0.67086597984190  0.57913402015810  0.63181557976744  79.904 > 61
  105 Br  0.42086597984190  0.17086597984190  0.63181557976744  79.904 > 65
  106 Br  0.92086597984190  0.17086597984190  0.63181557976744  79.904 > 65
  107 Br  0.42086597984190  0.67086597984190  0.63181557976744  79.904 > 65
  108 Br  0.92086597984190  0.67086597984190  0.63181557976744  79.904 > 65
  109 Br  0.17086597984190  0.42086597984190  0.86818442023256  79.904 > 69
  110 Br  0.67086597984190  0.42086597984190  0.86818442023256  79.904 > 69
  111 Br  0.17086597984190  0.92086597984190  0.86818442023256  79.904 > 69
  112 Br  0.67086597984190  0.92086597984190  0.86818442023256  79.904 > 69
  113 Br  0.07913402015810  0.17086597984190  0.86818442023256  79.904 > 73
  114 Br  0.57913402015810  0.17086597984190  0.86818442023256  79.904 > 73
  115 Br  0.07913402015810  0.67086597984190  0.86818442023256  79.904 > 73
  116 Br  0.57913402015810  0.67086597984190  0.86818442023256  79.904 > 73
  117 Br  0.25000000000000  0.25000000000000  0.50000000000000  79.904 > 77
  118 Br  0.75000000000000  0.25000000000000  0.50000000000000  79.904 > 77
  119 Br  0.25000000000000  0.75000000000000  0.50000000000000  79.904 > 77
  120 Br  0.75000000000000  0.75000000000000  0.50000000000000  79.904 > 77
  121 Br  0.32913402015810  0.07913402015810  0.86818442023256  79.904 > 81
  122 Br  0.82913402015810  0.07913402015810  0.86818442023256  79.904 > 81
  123 Br  0.32913402015810  0.57913402015810  0.86818442023256  79.904 > 81
  124 Br  0.82913402015810  0.57913402015810  0.86818442023256  79.904 > 81
  125 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 85
  126 Br  0.50000000000000  0.00000000000000  0.50000000000000  79.904 > 85
  127 Br  0.00000000000000  0.50000000000000  0.50000000000000  79.904 > 85
  128 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 85
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            5.9284900    0.0000000    0.0000000
            0.0000000    5.9284900    0.0000000
            0.0000000    0.0000000    4.0652419
-------------------------- Born effective charges --------------------------
    1 Rb    1.2367491    0.0000000   -0.0000000
            0.0000000    1.2367491    0.0000000
           -0.0000000    0.0000000    1.7464545
    2 Rb    1.2367491    0.0000000   -0.0000000
            0.0000000    1.2367491    0.0000000
           -0.0000000    0.0000000    1.7464545
    3 Sn    4.3725586    0.1008536   -0.0000000
            0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    4 Sn    4.3725586   -0.1008536   -0.0000000
           -0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    5 Sn    4.3725586    0.1008536   -0.0000000
            0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    6 Sn    4.3725586   -0.1008536   -0.0000000
           -0.1008536    4.3725586    0.0000000
           -0.0000000    0.0000000    1.9842761
    7 Br   -1.7462628    0.2622448    0.2123810
            0.2622448   -1.7462628   -0.2123810
            0.0800262   -0.0800262   -1.3168415
    8 Br   -1.7462628   -0.2622448    0.2123810
           -0.2622448   -1.7462628    0.2123810
            0.0800262    0.0800262   -1.3168415
    9 Br   -1.7462628    0.2622448    0.2123810
            0.2622448   -1.7462628   -0.2123810
            0.0800262   -0.0800262   -1.3168415
   10 Br   -1.7462628    0.2622448   -0.2123810
            0.2622448   -1.7462628    0.2123810
           -0.0800262    0.0800262   -1.3168415
   11 Br   -1.7462628   -0.2622448   -0.2123810
           -0.2622448   -1.7462628   -0.2123810
           -0.0800262   -0.0800262   -1.3168415
   12 Br   -1.7462628   -0.2622448    0.2123810
           -0.2622448   -1.7462628    0.2123810
            0.0800262    0.0800262   -1.3168415
   13 Br   -1.7462628    0.2622448   -0.2123810
            0.2622448   -1.7462628    0.2123810
           -0.0800262    0.0800262   -1.3168415
   14 Br   -2.9968150    0.2846536   -0.0000000
           -0.2846536   -2.9968150    0.0000000
           -0.0000000    0.0000000   -0.4476406
   15 Br   -1.7462628   -0.2622448   -0.2123810
           -0.2622448   -1.7462628   -0.2123810
           -0.0800262   -0.0800262   -1.3168415
   16 Br   -2.9968150   -0.2846536   -0.0000000
            0.2846536   -2.9968150    0.0000000
           -0.0000000    0.0000000   -0.4476406
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000005 (yyy) -0.00000005 (yyy) -0.00000005 (yyy)
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 ... 
Generalized regular grid: [ 3 6 12 ]
Grid generation matrix:
  [ 0 6 6 ]
  [ 6 0 6 ]
  [ 3 3 0 ]
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.52, Number of G-points: 291, Lambda: 0.13
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/32) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 32
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.720   (   0.000    0.000    0.000)    0.000
   0.720   (   0.000    0.000    0.000)    0.000
   1.078   (   0.000    0.000    0.000)    0.000
   1.098   (   0.000    0.000    0.000)    0.000
   1.098   (   0.000    0.000    0.000)    0.000
   1.141   (   0.000    0.000    0.000)    0.000
   1.463   (   0.000    0.000    0.000)    0.000
   1.468   (   0.000    0.000    0.000)    0.000
   1.468   (   0.000    0.000    0.000)    0.000
   1.482   (   0.000    0.000    0.000)    0.000
   1.661   (   0.000    0.000    0.000)    0.000
   1.752   (   0.000    0.000    0.000)    0.000
   1.768   (   0.000    0.000    0.000)    0.000
   1.768   (   0.000    0.000    0.000)    0.000
   1.945   (   0.000    0.000    0.000)    0.000
   2.064   (   0.000    0.000    0.000)    0.000
   2.064   (   0.000    0.000    0.000)    0.000
   2.084   (   0.000    0.000    0.000)    0.000
   2.136   (   0.000    0.000    0.000)    0.000
   2.136   (   0.000    0.000    0.000)    0.000
   2.188   (   0.000    0.000    0.000)    0.000
   2.188   (   0.000    0.000    0.000)    0.000
   2.215   (   0.000    0.000    0.000)    0.000
   2.262   (   0.000    0.000    0.000)    0.000
   2.313   (   0.000    0.000    0.000)    0.000
   2.329   (   0.000    0.000    0.000)    0.000
   2.343   (   0.000    0.000    0.000)    0.000
   2.343   (   0.000    0.000    0.000)    0.000
   2.360   (   0.000    0.000    0.000)    0.000
   2.623   (   0.000    0.000    0.000)    0.000
   2.661   (   0.000    0.000    0.000)    0.000
   2.956   (   0.000    0.000    0.000)    0.000
   2.956   (   0.000    0.000    0.000)    0.000
   2.993   (   0.000    0.000    0.000)    0.000
   3.287   (   0.000    0.000    0.000)    0.000
   3.553   (   0.000    0.000    0.000)    0.000
   3.553   (   0.000    0.000    0.000)    0.000
   3.561   (   0.000    0.000    0.000)    0.000
   3.808   (   0.000    0.000    0.000)    0.000
   3.808   (   0.000    0.000    0.000)    0.000
   4.081   (   0.000    0.000    0.000)    0.000
   4.180   (   0.000    0.000    0.000)    0.000
   4.180   (   0.000    0.000    0.000)    0.000
   4.684   (   0.000    0.000    0.000)    0.000
   4.993   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/32) =======================
q-point: ( 0.00  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 67
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.525   (  10.432    0.000   -0.262)   10.436
   0.725   (  13.539    0.000    2.791)   13.824
   0.820   (  -6.511   -0.000    3.280)    7.291
   1.000   (  11.554    0.000    3.447)   12.057
   1.031   (   3.403    0.000   -1.903)    3.899
   1.099   (  -2.067   -0.000    3.112)    3.736
   1.180   (  11.793    0.000    3.658)   12.347
   1.309   (  -2.746   -0.000    0.973)    2.913
   1.328   (  -0.576   -0.000   -0.720)    0.922
   1.387   (  -0.242    0.000    1.458)    1.478
   1.471   (   1.088    0.000    0.515)    1.204
   1.472   (   0.413    0.000    0.264)    0.490
   1.568   (   1.741    0.000    2.321)    2.901
   1.618   (   2.486    0.000    4.163)    4.849
   1.640   (  -6.507   -0.000    1.040)    6.590
   1.758   (  -4.142   -0.000   -1.570)    4.430
   1.835   (  -4.564   -0.000   -0.727)    4.621
   1.855   (  -2.033   -0.000   -4.798)    5.212
   1.855   (  -7.823   -0.000   -2.265)    8.145
   1.951   (  -2.529   -0.000   -2.259)    3.391
   2.024   (  -1.975   -0.000   -1.668)    2.585
   2.077   (  -2.627   -0.000   -1.634)    3.094
   2.127   (  -1.608   -0.000   -0.915)    1.850
   2.129   (  -1.526   -0.000   -0.747)    1.699
   2.247   (   1.315    0.000    1.219)    1.793
   2.248   (   3.136    0.000   -1.613)    3.527
   2.302   (   1.281    0.000   -1.349)    1.860
   2.371   (  -0.366   -0.000   -0.584)    0.689
   2.376   (  -0.650   -0.000    0.370)    0.748
   2.437   (   0.423    0.000   -0.747)    0.859
   2.503   (  -1.065   -0.000    1.318)    1.694
   2.544   (  -0.294    0.000    1.798)    1.822
   2.668   (   5.598    0.000    1.484)    5.791
   2.911   (  -3.883   -0.000   -0.094)    3.884
   3.047   (   3.694    0.000   -0.847)    3.789
   3.076   (   0.645    0.000    0.570)    0.861
   3.268   (   4.852    0.000   -0.230)    4.857
   3.451   (   1.961    0.000   -1.544)    2.496
   3.535   (   1.347    0.000    0.017)    1.347
   3.607   (   0.925   -0.000   -7.020)    7.081
   3.690   (   0.747    0.000    1.378)    1.568
   3.728   (  -2.058   -0.000   -2.257)    3.054
   4.067   (   0.877    0.000   -1.401)    1.653
   4.104   (   2.159    0.000   -2.182)    3.070
   4.208   (  10.688    0.000   -2.001)   10.874
   4.507   (  -1.606   -0.000    2.695)    3.137
   4.879   (  -4.201   -0.000   -2.175)    4.731
   4.924   (  -4.500   -0.000    1.512)    4.748
======================= Grid point 3 (3/32) =======================
q-point: ( 0.00  0.00  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 48
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.379   (   8.024    8.024    0.000)   11.348
   0.394   (   9.574    9.574    0.000)   13.539
   0.761   (  16.446   16.446    0.000)   23.258
   0.794   (   3.082    3.082    0.000)    4.359
   0.912   (  -0.065   -0.065    0.000)    0.093
   0.976   (  -2.012   -2.012    0.000)    2.845
   1.270   (   2.297    2.297    0.000)    3.248
   1.279   (   1.579    1.579    0.000)    2.233
   1.329   (  -4.611   -4.611    0.000)    6.520
   1.372   (   0.853    0.853    0.000)    1.207
   1.405   (  -4.618   -4.618    0.000)    6.531
   1.483   (  -6.176   -6.176    0.000)    8.734
   1.510   (   0.340    0.340    0.000)    0.480
   1.527   (  -1.502   -1.502    0.000)    2.124
   1.748   (  -1.206   -1.206    0.000)    1.705
   1.772   (  -5.938   -5.938    0.000)    8.397
   1.925   (   0.100    0.100    0.000)    0.141
   1.999   (   0.110    0.110    0.000)    0.156
   2.023   (  -2.558   -2.558    0.000)    3.618
   2.037   (   2.406    2.406    0.000)    3.403
   2.090   (  -2.515   -2.515    0.000)    3.557
   2.101   (  -1.177   -1.177    0.000)    1.665
   2.129   (  -2.104   -2.104    0.000)    2.975
   2.150   (  -6.721   -6.721    0.000)    9.505
   2.186   (   1.342    1.342    0.000)    1.898
   2.274   (   2.391    2.391    0.000)    3.382
   2.300   (   1.714    1.714    0.000)    2.424
   2.350   (   0.680    0.680    0.000)    0.962
   2.358   (   1.480    1.480    0.000)    2.093
   2.414   (   0.395    0.395    0.000)    0.559
   2.486   (   2.799    2.799    0.000)    3.958
   2.545   (  -1.795   -1.795    0.000)    2.539
   2.616   (  -0.940   -0.940    0.000)    1.329
   2.949   (  -2.078   -2.078    0.000)    2.939
   3.030   (   2.273    2.273    0.000)    3.214
   3.035   (   5.473    5.473    0.000)    7.740
   3.214   (  -1.664   -1.664    0.000)    2.354
   3.512   (   4.049    4.049    0.000)    5.727
   3.530   (  -0.183   -0.183    0.000)    0.258
   3.619   (   1.089    1.089    0.000)    1.540
   3.804   (  -0.036   -0.036    0.000)    0.051
   3.817   (   0.396    0.396    0.000)    0.560
   3.981   (   6.534    6.534    0.000)    9.240
   4.083   (   0.095    0.095    0.000)    0.135
   4.117   (  -2.853   -2.853    0.000)    4.035
   4.567   (  -5.482   -5.482    0.000)    7.752
   4.900   (  -1.226   -1.226    0.000)    1.734
   4.980   (  -0.667   -0.667    0.000)    0.943
======================= Grid point 4 (4/32) =======================
q-point: ( 0.00  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.726   (  -7.484    4.847    1.829)    9.102
   0.785   (   9.603    3.955   -3.101)   10.838
   0.894   (   6.021   -7.253    2.076)    9.652
   0.950   (  -5.064   -6.898   -2.475)    8.908
   1.082   (  -1.672    3.208    0.613)    3.669
   1.114   (   1.516    0.917   -1.988)    2.663
   1.217   (  -0.655   -4.968    1.149)    5.141
   1.228   (   3.719   -5.435    0.330)    6.594
   1.349   (  12.330   -7.590    3.329)   14.856
   1.428   (  -7.090   -4.077   -2.579)    8.576
   1.484   (  -4.646   -2.687    1.139)    5.487
   1.515   (   0.528    1.546   -0.615)    1.746
   1.574   (  -0.021   -1.014   -0.402)    1.091
   1.598   (   1.214    1.072   -0.973)    1.889
   1.621   (   3.394    0.731    2.901)    4.524
   1.711   (   1.539    2.999   -3.682)    4.992
   1.742   (   4.385    2.731    2.483)    5.732
   1.808   (  -3.955    0.431    1.819)    4.374
   1.860   (  -2.808    2.476   -2.160)    4.322
   1.888   (   0.387    0.273    3.072)    3.108
   1.959   (  -1.232    0.219   -3.907)    4.103
   2.003   (  -1.887   -1.216   -4.083)    4.660
   2.084   (  -2.156   -0.994    3.264)    4.036
   2.115   (   1.422   -2.425    0.432)    2.844
   2.292   (   0.655    1.162    1.039)    1.691
   2.308   (  -6.013    0.464    0.501)    6.052
   2.310   (  -0.591    1.817    0.204)    1.921
   2.322   (  -2.440    1.197    0.719)    2.811
   2.405   (   3.478   -1.125   -1.073)    3.809
   2.418   (   1.834    2.290   -1.466)    3.280
   2.510   (   1.985   -0.578    1.722)    2.691
   2.531   (  -0.602   -0.290   -0.170)    0.690
   2.790   (   5.801   -1.218    0.118)    5.929
   2.802   (  -4.715   -1.702   -0.852)    5.084
   3.094   (  -0.547    2.255    1.153)    2.591
   3.118   (   0.408    2.015   -0.865)    2.230
   3.374   (  -7.728   -0.184    0.566)    7.751
   3.391   (   7.608   -0.100   -0.775)    7.648
   3.589   (  -0.146    2.158    1.261)    2.503
   3.632   (  -3.379    3.265   -1.902)    5.069
   3.662   (  -0.830    1.231    3.495)    3.797
   3.742   (  -0.461   -0.100   -3.289)    3.322
   4.047   (  -1.071   -2.151    0.582)    2.472
   4.069   (   0.371   -1.405   -1.524)    2.106
   4.300   ( -11.339   -3.322    7.006)   13.737
   4.451   (  13.362   -4.299   -5.715)   15.155
   4.808   (  -5.091   -0.014    3.886)    6.404
   4.902   (   2.376    0.552   -3.873)    4.577
======================= Grid point 5 (5/32) =======================
q-point: ( 0.00 -0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.432   (  -6.470    6.470   -4.362)   10.136
   0.467   (  -5.900    5.900   -5.900)   10.220
   0.817   (  -3.521    3.521   -2.518)    5.580
   0.845   ( -13.194   13.194   -5.505)   19.454
   0.951   (   0.233   -0.233   -3.418)    3.434
   1.009   (   1.689   -1.689   -2.657)    3.573
   1.211   (  -0.397    0.397    2.111)    2.185
   1.277   (  -1.162    1.162   -0.005)    1.643
   1.340   (   5.177   -5.177   -0.627)    7.349
   1.402   (  -1.878    1.878   -2.587)    3.707
   1.472   (   2.736   -2.736   -1.594)    4.185
   1.481   (   3.942   -3.942   -0.256)    5.581
   1.550   (  -1.384    1.384   -4.709)    5.100
   1.571   (   1.512   -1.512   -3.590)    4.179
   1.752   (   2.623   -2.623   -0.346)    3.726
   1.774   (   5.079   -5.079   -0.163)    7.184
   1.887   (  -1.186    1.186    2.132)    2.712
   1.923   (   1.269   -1.269    7.625)    7.833
   1.971   (   0.689   -0.689    3.764)    3.888
   1.984   (   1.329   -1.329    5.540)    5.850
   1.999   (  -0.294    0.294    2.543)    2.576
   2.059   (  -3.586    3.586    2.822)    5.804
   2.139   (   6.556   -6.556    0.107)    9.273
   2.193   (   1.863   -1.863   -0.941)    2.798
   2.203   (   1.522   -1.522   -1.538)    2.646
   2.284   (  -2.027    2.027    1.014)    3.040
   2.295   (  -2.012    2.012    2.327)    3.676
   2.360   (  -1.351    1.351   -1.656)    2.528
   2.372   (  -0.204    0.204    2.762)    2.777
   2.400   (  -0.380    0.380   -3.322)    3.365
   2.493   (  -3.876    3.876   -0.329)    5.491
   2.531   (   1.828   -1.828    0.381)    2.613
   2.637   (  -0.130    0.130   -2.091)    2.099
   2.951   (   2.137   -2.137   -0.222)    3.030
   3.021   (  -5.534    5.534    1.362)    7.944
   3.033   (  -2.303    2.303   -0.284)    3.269
   3.218   (   1.728   -1.728   -0.352)    2.468
   3.419   ( -10.928   10.928   10.731)   18.815
   3.526   (  -0.208    0.208    0.331)    0.443
   3.620   (  -0.896    0.896   -0.208)    1.284
   3.650   (  -1.751    1.751    5.592)    6.116
   3.793   (  -0.406    0.406    1.949)    2.032
   3.954   (  -6.140    6.140    3.348)    9.306
   4.057   (   0.269   -0.269    1.066)    1.132
   4.126   (   2.637   -2.637   -0.734)    3.801
   4.466   (  -6.109    6.109    4.043)    9.539
   4.963   (   4.481   -4.481   -2.964)    6.997
   4.985   (   4.955   -4.955   -1.078)    7.089
======================= Grid point 6 (6/32) =======================
q-point: ( 0.00  0.00  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 50
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.543   (  -1.658   -1.658    0.000)    2.345
   0.768   (   8.392    8.392    0.000)   11.868
   0.885   (  -5.125   -5.125    0.000)    7.248
   0.960   (   0.893    0.893    0.000)    1.263
   1.026   (  -5.284   -5.284    0.000)    7.473
   1.032   (   3.557    3.557    0.000)    5.030
   1.116   (   2.645    2.645    0.000)    3.741
   1.188   (  -2.927   -2.927    0.000)    4.140
   1.237   (  -2.265   -2.265    0.000)    3.203
   1.323   (   1.057    1.057    0.000)    1.495
   1.481   (  -5.636   -5.636    0.000)    7.970
   1.494   (   4.804    4.804    0.000)    6.794
   1.527   (   0.208    0.208    0.000)    0.295
   1.599   (   9.892    9.892    0.000)   13.990
   1.628   (  -4.702   -4.702    0.000)    6.650
   1.656   (   5.362    5.362    0.000)    7.583
   1.770   (  -7.905   -7.905    0.000)   11.179
   1.896   (  -0.619   -0.619    0.000)    0.875
   1.903   (  -3.903   -3.903    0.000)    5.520
   1.978   (  -1.964   -1.964    0.000)    2.777
   2.002   (  -1.882   -1.882    0.000)    2.662
   2.062   (  -1.322   -1.322    0.000)    1.870
   2.075   (  -0.039   -0.039    0.000)    0.055
   2.117   (  -0.585   -0.585    0.000)    0.828
   2.250   (   1.843    1.843    0.000)    2.606
   2.338   (  -1.642   -1.642    0.000)    2.322
   2.355   (   1.106    1.106    0.000)    1.564
   2.397   (   2.399    2.399    0.000)    3.393
   2.406   (  -3.639   -3.639    0.000)    5.146
   2.417   (   1.258    1.258    0.000)    1.779
   2.446   (   4.144    4.144    0.000)    5.861
   2.526   (  -0.871   -0.871    0.000)    1.232
   2.691   (   2.203    2.203    0.000)    3.115
   2.835   (  -3.437   -3.437    0.000)    4.860
   3.096   (   1.195    1.195    0.000)    1.690
   3.175   (  -0.708   -0.708    0.000)    1.002
   3.282   (   6.347    6.347    0.000)    8.976
   3.548   (  -4.214   -4.214    0.000)    5.959
   3.578   (  -0.719   -0.719    0.000)    1.016
   3.584   (   2.795    2.795    0.000)    3.953
   3.778   (  -2.414   -2.414    0.000)    3.415
   3.834   (   0.361    0.361    0.000)    0.510
   4.022   (  -1.037   -1.037    0.000)    1.466
   4.094   (   0.706    0.706    0.000)    0.998
   4.223   (   4.453    4.453    0.000)    6.297
   4.282   (  -7.629   -7.629    0.000)   10.789
   4.813   (  -2.880   -2.880    0.000)    4.073
   4.962   (   0.000    0.000    0.000)    0.001
======================= Grid point 7 (7/32) =======================
q-point: ( 0.00 -0.67  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 67
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.615   (   3.201   -3.201   -0.788)    4.595
   0.705   (  -6.682    6.682    1.247)    9.532
   0.890   (  10.418  -10.418   -1.234)   14.785
   1.002   (   2.209   -2.209   -0.554)    3.173
   1.074   (  -1.535    1.535   -0.258)    2.186
   1.145   (  -2.868    2.868   -1.694)    4.395
   1.149   (   0.286   -0.286   -3.129)    3.155
   1.215   (  -0.388    0.388   -3.856)    3.895
   1.307   (   1.656   -1.656   -4.255)    4.857
   1.326   (   4.703   -4.703   -1.142)    6.748
   1.418   (   8.811   -8.811    4.339)   13.194
   1.514   (  -0.979    0.979   -0.342)    1.427
   1.562   (   1.332   -1.332   -1.700)    2.538
   1.663   (  -0.970    0.970    0.106)    1.376
   1.693   (  -0.383    0.383    8.837)    8.854
   1.725   (  -0.118    0.118    1.039)    1.052
   1.753   (  -1.401    1.401   -4.713)    5.112
   1.804   (  -2.534    2.534   -2.495)    4.366
   1.910   (   0.016   -0.016    1.118)    1.118
   1.931   (  -1.838    1.838   -1.527)    3.015
   1.962   (   3.104   -3.104   -4.089)    5.999
   1.978   (  -1.760    1.760    3.410)    4.222
   2.075   (   4.135   -4.135    1.536)    6.046
   2.089   (   2.167   -2.167    1.985)    3.651
   2.194   (  -3.136    3.136   -0.238)    4.441
   2.322   (  -2.221    2.221    0.865)    3.258
   2.331   (  -1.746    1.746    1.754)    3.028
   2.334   (  -3.667    3.667    4.175)    6.658
   2.417   (  -0.092    0.092   -0.433)    0.452
   2.470   (   3.211   -3.211   -1.501)    4.782
   2.522   (  -1.224    1.224   -3.137)    3.582
   2.546   (   0.874   -0.874    0.303)    1.272
   2.702   (  -0.577    0.577   -0.893)    1.210
   2.837   (   3.582   -3.582   -0.313)    5.075
   3.115   (  -1.407    1.407   -0.238)    2.004
   3.176   (   1.233   -1.233   -0.144)    1.750
   3.229   (  -6.330    6.330    1.541)    9.084
   3.559   (  -2.686    2.686    4.899)    6.199
   3.568   (   3.107   -3.107   -1.103)    4.530
   3.624   (  -2.453    2.453   -1.960)    3.985
   3.654   (  -2.904    2.904   -1.076)    4.245
   3.733   (   0.333   -0.333    2.880)    2.918
   4.016   (   0.206   -0.206    0.968)    1.011
   4.056   (  -3.644    3.644    5.176)    7.304
   4.094   (   0.005   -0.005   -1.572)    1.572
   4.533   (   7.348   -7.348   -5.900)   11.950
   4.757   (  -0.764    0.764    5.277)    5.386
   4.937   (   1.263   -1.263   -2.770)    3.296
======================= Grid point 9 (8/32) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 28
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.371   (   0.000    0.000    0.000)    0.000
   0.475   (   0.000    0.000    0.000)    0.000
   0.982   (   0.000    0.000    0.000)    0.000
   0.993   (   0.000    0.000    0.000)    0.000
   1.079   (   0.000    0.000    0.000)    0.000
   1.103   (   0.000    0.000    0.000)    0.000
   1.136   (   0.000    0.000    0.000)    0.000
   1.169   (   0.000    0.000    0.000)    0.000
   1.201   (   0.000    0.000    0.000)    0.000
   1.338   (   0.000    0.000    0.000)    0.000
   1.392   (   0.000    0.000    0.000)    0.000
   1.456   (  -0.000   -0.000    0.000)    0.000
   1.509   (   0.000    0.000    0.000)    0.000
   1.602   (   0.000    0.000    0.000)    0.000
   1.633   (   0.000    0.000    0.000)    0.000
   1.758   (  -0.000   -0.000    0.000)    0.000
   1.839   (   0.000    0.000    0.000)    0.000
   1.880   (  -0.000   -0.000    0.000)    0.000
   1.928   (  -0.000   -0.000    0.000)    0.000
   1.930   (   0.000    0.000    0.000)    0.000
   1.939   (   0.000    0.000    0.000)    0.000
   2.006   (   0.000    0.000    0.000)    0.000
   2.035   (   0.000    0.000    0.000)    0.000
   2.078   (   0.000    0.000    0.000)    0.000
   2.281   (   0.000    0.000    0.000)    0.000
   2.337   (   0.000    0.000    0.000)    0.000
   2.366   (   0.000    0.000    0.000)    0.000
   2.377   (   0.000    0.000    0.000)    0.000
   2.411   (   0.000    0.000    0.000)    0.000
   2.445   (   0.000    0.000    0.000)    0.000
   2.485   (   0.000    0.000    0.000)    0.000
   2.589   (   0.000    0.000    0.000)    0.000
   2.714   (   0.000    0.000    0.000)    0.000
   2.736   (   0.000    0.000    0.000)    0.000
   3.127   (   0.000    0.000    0.000)    0.000
   3.158   (   0.000    0.000    0.000)    0.000
   3.425   (   0.000    0.000    0.000)    0.000
   3.520   (   0.000    0.000    0.000)    0.000
   3.525   (  -0.000   -0.000    0.000)    0.000
   3.623   (   0.000    0.000    0.000)    0.000
   3.687   (   0.000    0.000    0.000)    0.000
   3.842   (   0.000    0.000    0.000)    0.000
   4.010   (   0.000    0.000    0.000)    0.000
   4.069   (   0.000    0.000    0.000)    0.000
   4.125   (   0.000    0.000    0.000)    0.000
   4.318   (   0.000    0.000    0.000)    0.000
   4.733   (   0.000    0.000    0.000)    0.000
   4.968   (   0.000    0.000    0.000)    0.000
======================= Grid point 19 (9/32) =======================
q-point: ( 0.08 -0.08 -0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 42
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.280   ( -13.662   -0.000   -0.000)   13.662
   0.361   ( -17.163   -0.000   -0.000)   17.163
   0.506   ( -23.954   -0.000   -0.000)   23.954
   0.845   ( -11.025   -0.000   -0.000)   11.025
   0.903   (   6.887    0.000    0.000)    6.887
   1.089   (   0.986    0.000    0.000)    0.986
   1.096   (   2.114    0.000    0.000)    2.114
   1.327   ( -10.388   -0.000   -0.000)   10.388
   1.349   (   1.961    0.000    0.000)    1.961
   1.428   (   3.395    0.000    0.000)    3.395
   1.453   (   2.318    0.000    0.000)    2.318
   1.469   (  -0.205   -0.000   -0.000)    0.205
   1.571   (   3.751    0.000    0.000)    3.751
   1.579   (  10.523    0.000    0.000)   10.523
   1.727   (   3.708    0.000    0.000)    3.708
   1.878   (  -5.728   -0.000   -0.000)    5.728
   1.898   (   4.971    0.000    0.000)    4.971
   1.990   (  -2.054   -0.000   -0.000)    2.054
   2.000   (   4.820    0.000    0.000)    4.820
   2.004   (   6.394    0.000    0.000)    6.394
   2.110   (   3.534    0.000    0.000)    3.534
   2.126   (   1.572    0.000    0.000)    1.572
   2.158   (   0.403    0.000    0.000)    0.403
   2.164   (  -3.191   -0.000   -0.000)    3.191
   2.252   (  -3.042   -0.000   -0.000)    3.042
   2.260   (   5.191    0.000    0.000)    5.191
   2.260   (   3.084    0.000    0.000)    3.084
   2.346   (   0.075    0.000    0.000)    0.075
   2.350   (  -3.556   -0.000   -0.000)    3.556
   2.406   (  -3.872   -0.000   -0.000)    3.872
   2.441   (  -7.078   -0.000   -0.000)    7.078
   2.571   (   4.206    0.000    0.000)    4.206
   2.633   (   2.395    0.000    0.000)    2.395
   2.960   (  -4.932   -0.000   -0.000)    4.932
   2.971   (   2.103    0.000    0.000)    2.103
   3.017   (  -4.494   -0.000   -0.000)    4.494
   3.242   (   2.314    0.000    0.000)    2.314
   3.430   (  -3.171   -0.000   -0.000)    3.171
   3.538   (   1.160    0.000    0.000)    1.160
   3.619   (  -4.528   -0.000   -0.000)    4.528
   3.778   (   2.794    0.000    0.000)    2.794
   3.782   (   1.908    0.000    0.000)    1.908
   3.942   ( -11.707   -0.000   -0.000)   11.707
   4.083   (  -0.229   -0.000   -0.000)    0.229
   4.166   (   1.365    0.000    0.000)    1.365
   4.624   (   5.799    0.000    0.000)    5.799
   4.923   (   0.142    0.000    0.000)    0.142
   4.975   (   1.918    0.000    0.000)    1.918
======================= Grid point 20 (10/32) =======================
q-point: ( 0.08  0.25 -0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.341   (   6.692    0.000    7.840)   10.308
   0.450   (  12.713    0.000    6.213)   14.150
   0.664   (  17.907    0.000    9.542)   20.291
   0.857   (  10.886    0.000    1.441)   10.981
   0.943   (  -6.443   -0.000    3.418)    7.293
   1.094   (  -1.160   -0.000   -0.485)    1.258
   1.119   (   0.014    0.000    2.691)    2.691
   1.286   (   5.484    0.000   -0.782)    5.539
   1.349   (   1.239    0.000   -0.153)    1.248
   1.436   (  -2.080   -0.000   -1.035)    2.323
   1.471   (  -0.688   -0.000    0.913)    1.144
   1.536   (   1.793    0.000    5.136)    5.440
   1.547   (  -5.740   -0.000   -2.129)    6.123
   1.610   (  -3.221   -0.000    3.454)    4.723
   1.747   (  -3.471   -0.000    1.362)    3.729
   1.815   (   2.645    0.000   -4.184)    4.950
   1.865   (   1.834    0.000   -1.236)    2.212
   1.934   (  -2.881   -0.000   -5.633)    6.327
   1.960   (  -0.707   -0.000   -5.196)    5.244
   1.991   (  -3.453   -0.000   -2.231)    4.112
   2.015   (   0.074   -0.000   -3.544)    3.544
   2.021   (   2.448    0.000   -3.984)    4.676
   2.211   (  -4.008   -0.000   -1.159)    4.172
   2.218   (  -1.293   -0.000    1.006)    1.638
   2.232   (  -3.148   -0.000   -0.846)    3.260
   2.243   (  -1.678   -0.000    0.655)    1.801
   2.277   (   1.143    0.000   -3.339)    3.529
   2.341   (   2.667    0.000    0.044)    2.667
   2.374   (   5.098    0.000   -0.222)    5.102
   2.390   (  -0.272    0.000    2.341)    2.357
   2.452   (   6.729    0.000    1.689)    6.938
   2.553   (  -3.116   -0.000   -1.313)    3.382
   2.623   (   3.278    0.000    0.628)    3.338
   2.949   (   5.341    0.000   -1.314)    5.501
   2.973   (  -2.194   -0.000    0.244)    2.208
   3.018   (   4.327    0.000    0.085)    4.328
   3.239   (  18.616    0.000  -14.379)   23.522
   3.247   (  -2.452   -0.000    0.397)    2.484
   3.529   (  -0.692   -0.000   -0.751)    1.021
   3.600   (   1.739    0.000   -4.225)    4.569
   3.626   (   4.690    0.000    0.592)    4.727
   3.770   (  -1.371   -0.000   -2.076)    2.488
   3.902   (  10.427    0.000   -3.132)   10.887
   4.061   (  -0.482   -0.000   -0.832)    0.961
   4.176   (  -0.999   -0.000    0.659)    1.197
   4.365   (  12.628    0.000   -5.620)   13.822
   4.988   (  -1.883   -0.000    1.026)    2.145
   5.077   ( -12.693   -0.000    4.290)   13.399
======================= Grid point 21 (11/32) =======================
q-point: ( 0.08 -0.42 -0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.676   (   6.653    0.000   -1.445)    6.808
   0.739   ( -10.115   -0.000    3.739)   10.784
   0.973   (  -9.542   -0.000   -1.760)    9.703
   1.012   (   8.658    0.000    1.403)    8.771
   1.041   (   2.369    0.000   -1.585)    2.850
   1.103   (  -3.689   -0.000    3.853)    5.334
   1.288   (  -0.349   -0.000   -2.281)    2.308
   1.295   (  -2.769   -0.000   -1.638)    3.217
   1.438   (  -6.964   -0.000   -1.490)    7.122
   1.466   (   4.704    0.000    0.721)    4.759
   1.474   (   0.109   -0.000   -1.595)    1.598
   1.500   (   1.908    0.000    0.247)    1.924
   1.561   (  -6.341   -0.000   -2.589)    6.849
   1.611   (  -1.359   -0.000    1.272)    1.862
   1.666   (  -4.357   -0.000    0.343)    4.370
   1.675   (   2.118    0.000    0.214)    2.129
   1.723   (   2.879    0.000    2.606)    3.884
   1.764   (   4.392    0.000    3.814)    5.817
   1.834   (  -2.364   -0.000   -2.339)    3.326
   1.886   (   2.337    0.000   -1.607)    2.837
   1.967   (   2.599    0.000    1.909)    3.225
   2.018   (   0.667    0.000    3.829)    3.887
   2.094   (   1.935    0.000   -3.620)    4.105
   2.140   (  -1.514   -0.000    0.108)    1.518
   2.276   (  -1.197   -0.000   -2.863)    3.103
   2.283   (   1.128    0.000   -1.073)    1.556
   2.296   (   3.772    0.000    0.633)    3.824
   2.321   (   3.746    0.000    1.045)    3.889
   2.401   (   1.022    0.000   -0.052)    1.023
   2.421   (  -4.199   -0.000    1.913)    4.614
   2.516   (  -2.191   -0.000   -1.485)    2.647
   2.532   (   0.152    0.000   -0.062)    0.164
   2.802   (  -6.400   -0.000   -0.244)    6.404
   2.819   (   5.002    0.000    1.144)    5.132
   3.066   (   1.254    0.000   -1.276)    1.789
   3.093   (  -0.940   -0.000    0.906)    1.305
   3.379   (   6.232    0.000   -0.606)    6.262
   3.396   (  -6.589   -0.000    0.818)    6.640
   3.564   (   0.855    0.000   -0.844)    1.201
   3.595   (   2.878    0.000    1.363)    3.185
   3.645   (   3.275    0.000   -3.788)    5.008
   3.736   (  -1.722   -0.000    3.920)    4.282
   4.077   (   2.077    0.000   -0.259)    2.093
   4.090   (  -1.133   -0.000    1.065)    1.555
   4.334   (  11.962    0.000   -7.432)   14.083
   4.495   ( -13.728   -0.000    6.153)   15.044
   4.807   (   6.268    0.000   -3.603)    7.230
   4.895   (  -3.289   -0.000    3.695)    4.947
======================= Grid point 23 (12/32) =======================
q-point: ( 0.08  0.25  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.585   (   9.327    5.293   -0.526)   10.738
   0.671   (  11.683   -3.846    2.941)   12.646
   0.874   (  -5.835    4.558    2.623)    7.856
   0.984   (   7.947   -3.528    0.937)    8.745
   1.051   (   3.143    0.925    0.480)    3.311
   1.094   (  10.458   -1.185    3.928)   11.234
   1.139   (  -4.697    1.142   -0.359)    4.847
   1.242   (  -3.862   -6.684    0.427)    7.732
   1.285   (  -2.905   -3.329    0.839)    4.497
   1.430   (   2.405   -1.791    0.295)    3.013
   1.442   (   3.044   -1.672    0.230)    3.481
   1.488   (   4.543    5.868    4.760)    8.817
   1.559   (   0.197    0.568    2.934)    2.995
   1.611   (   1.515   -1.286    4.215)    4.660
   1.628   (  -4.612   -1.163    0.463)    4.779
   1.691   (  -4.627   -6.183   -1.378)    7.844
   1.813   (  -5.964    1.275   -3.923)    7.252
   1.874   (  -3.582    2.055   -1.363)    4.348
   1.931   (  -2.539    2.159   -4.717)    5.775
   1.945   (  -3.987   -1.619   -0.622)    4.348
   1.990   (  -2.861   -3.023   -2.126)    4.673
   2.056   (   0.051   -1.682   -3.094)    3.522
   2.108   (  -1.445    0.322    1.784)    2.318
   2.136   (  -2.516   -0.953    0.502)    2.737
   2.217   (   4.451   -1.485   -1.988)    5.095
   2.307   (   0.332    4.067   -0.754)    4.149
   2.333   (   1.413    2.728   -1.476)    3.408
   2.353   (  -0.319   -0.997   -2.030)    2.284
   2.394   (  -0.756    1.935    3.066)    3.704
   2.426   (   1.877   -1.020   -0.497)    2.193
   2.479   (  -2.321   -2.301    1.929)    3.795
   2.549   (   0.146    0.264    0.966)    1.012
   2.673   (   4.748    0.345    1.072)    4.880
   2.890   (  -3.785   -1.988   -0.030)    4.275
   3.076   (   1.947    1.597    0.144)    2.523
   3.114   (   1.208    4.349   -0.349)    4.527
   3.240   (   5.489   -1.729   -0.382)    5.768
   3.514   (  -3.672    3.824    0.213)    5.305
   3.549   (   1.974    1.489    0.387)    2.503
   3.638   (   2.852    1.860   -5.272)    6.276
   3.676   (   0.142   -0.112   -1.797)    1.806
   3.769   (  -2.188    2.443   -2.041)    3.862
   4.057   (   0.361   -0.740   -1.571)    1.774
   4.068   (  -0.182   -2.257   -0.733)    2.380
   4.177   (  11.739   -2.237   -2.868)   12.290
   4.490   (  -5.279   -2.535    5.167)    7.809
   4.859   (  -1.796   -1.280   -2.042)    3.005
   4.922   (  -3.462   -0.436   -0.060)    3.490
======================= Grid point 25 (13/32) =======================
q-point: ( 0.08 -0.08  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.615   (   5.937   11.420   -0.000)   12.870
   0.616   (  -4.479   13.279   -0.000)   14.014
   0.843   (   5.569   -6.059    0.000)    8.229
   0.983   (  -7.513   11.481   -0.000)   13.721
   1.031   (  -0.654    5.173   -0.000)    5.214
   1.037   (   5.818    9.143   -0.000)   10.837
   1.175   (   5.434   -3.370    0.000)    6.394
   1.229   (  -8.324   -5.990   -0.000)   10.256
   1.270   (  -5.394   -6.027    0.000)    8.088
   1.426   (   1.337    3.076   -0.000)    3.354
   1.429   (  -7.262    6.328   -0.000)    9.632
   1.438   (   4.776    5.599   -0.000)    7.359
   1.515   (   1.116    0.065    0.000)    1.118
   1.561   (  -0.547    3.650   -0.000)    3.690
   1.617   (  -1.814   -5.516    0.000)    5.806
   1.717   (  -4.351   -3.347   -0.000)    5.489
   1.875   (   1.679   -8.981    0.000)    9.137
   1.889   (   0.955   -3.115    0.000)    3.258
   1.914   (  -2.579   -8.915    0.000)    9.280
   2.034   (   4.443   -2.395    0.000)    5.048
   2.041   (  -2.118   -1.900   -0.000)    2.846
   2.074   (  -5.074   -1.017   -0.000)    5.175
   2.089   (   1.131   -1.615    0.000)    1.972
   2.102   (  -2.343   -0.987   -0.000)    2.542
   2.256   (  -1.100    5.339   -0.000)    5.451
   2.304   (   5.247   -1.147    0.000)    5.371
   2.328   (   0.982    1.990   -0.000)    2.219
   2.367   (   3.094    1.542    0.000)    3.457
   2.377   (  -1.776   -1.336   -0.000)    2.222
   2.430   (  -1.987    2.221   -0.000)    2.980
   2.468   (  -2.089   -4.461    0.000)    4.926
   2.526   (   1.086   -0.095    0.000)    1.090
   2.665   (   0.728    5.595   -0.000)    5.643
   2.892   (  -2.053   -3.451    0.000)    4.015
   3.072   (   1.123    1.565   -0.000)    1.926
   3.119   (   4.671    0.946    0.000)    4.766
   3.245   (  -1.470    6.154   -0.000)    6.327
   3.510   (   3.247   -4.506    0.000)    5.554
   3.544   (   1.110    1.649   -0.000)    1.988
   3.649   (  -2.195    0.613   -0.000)    2.280
   3.760   (   4.854   -4.660    0.000)    6.729
   3.792   (   1.343   -0.210    0.000)    1.359
   4.070   (  -3.281   -1.561   -0.000)    3.634
   4.093   (  -0.826    0.461   -0.000)    0.946
   4.213   (  -1.662   14.433   -0.000)   14.528
   4.410   (  -4.526   -9.702    0.000)   10.706
   4.864   (  -2.317   -3.031    0.000)    3.815
   4.954   (   1.052   -1.249    0.000)    1.633
======================= Grid point 26 (14/32) =======================
q-point: ( 0.08  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.634   (   3.585  -15.938    2.222)   16.487
   0.688   (  -2.626  -16.627   -2.066)   16.960
   0.865   (  -8.166    7.902    1.854)   11.514
   0.908   (   7.166    5.838   -1.629)    9.386
   1.116   (  -1.784    0.180   -0.069)    1.794
   1.130   (   2.685   -0.097   -1.078)    2.895
   1.159   (   6.423   -1.001    1.380)    6.645
   1.174   (  -0.862    0.802   -0.185)    1.192
   1.241   (   5.502   -2.117    2.414)    6.370
   1.321   (  -2.525   -5.868   -4.111)    7.597
   1.371   (  -4.494   -6.764   -0.561)    8.140
   1.422   (   0.688   -6.211   -3.975)    7.406
   1.630   (  -0.278    2.526    4.985)    5.596
   1.646   (   4.492    3.749    3.678)    6.911
   1.694   (  -2.151    3.902   -1.538)    4.713
   1.746   (  -1.751   -0.692    1.530)    2.426
   1.780   (  -3.083   -1.069    0.820)    3.365
   1.823   (   6.735    4.383    1.614)    8.196
   1.884   (  -3.374    1.309   -1.888)    4.082
   1.923   (   3.111    0.183   -2.460)    3.970
   1.945   (   1.067    0.212   -1.401)    1.774
   2.011   (  -4.626    0.189   -2.617)    5.319
   2.053   (   0.904   -2.327    1.491)    2.908
   2.063   (  -2.082   -1.087    0.717)    2.456
   2.281   (  -3.008   -1.650    0.082)    3.432
   2.302   (   1.817   -1.154   -1.484)    2.614
   2.360   (   1.416    2.261    0.221)    2.677
   2.373   (  -1.364    3.508   -0.602)    3.812
   2.412   (  -2.076    2.640    3.989)    5.215
   2.483   (  -0.024   -0.687    1.100)    1.297
   2.498   (   2.319    2.842   -1.968)    4.162
   2.513   (   0.939   -1.603   -0.478)    1.918
   2.755   (  -2.234   -2.442   -0.164)    3.313
   2.756   (   2.473   -2.072   -0.243)    3.235
   3.130   (   0.806    0.871    0.809)    1.436
   3.148   (  -0.469    0.590   -0.655)    0.999
   3.393   (  -8.770    2.480    0.953)    9.163
   3.419   (   8.129    3.259   -1.118)    8.830
   3.606   (  -0.192   -2.171    1.977)    2.943
   3.665   (  -0.088   -1.994   -2.011)    2.834
   3.695   (   1.594    1.118    1.472)    2.440
   3.745   (  -4.578    3.472   -1.913)    6.056
   4.028   (  -1.197    0.641    1.507)    2.028
   4.066   (   0.879    1.282   -1.747)    2.339
   4.210   (  -8.571   -5.008    5.531)   11.364
   4.336   (  10.307   -6.252   -4.978)   13.042
   4.804   (  -2.658   -0.330    4.579)    5.305
   4.911   (   1.042    0.293   -4.219)    4.356
======================= Grid point 28 (15/32) =======================
q-point: ( 0.08 -0.08  0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.606   ( -16.011    4.272   -0.000)   16.571
   0.713   ( -17.367   -2.257   -0.000)   17.513
   0.842   (   8.592   -8.490    0.000)   12.079
   0.928   (   4.521    6.641   -0.000)    8.034
   1.120   (  -1.630    0.801   -0.000)    1.816
   1.134   (  -2.655    5.883   -0.000)    6.454
   1.150   (   0.522    4.169   -0.000)    4.202
   1.188   (  -0.141    4.075   -0.000)    4.077
   1.197   (  -0.269   -0.474    0.000)    0.545
   1.375   (  -3.930   -0.386   -0.000)    3.949
   1.430   (  -5.697    4.168   -0.000)    7.058
   1.503   (   1.446   -2.344    0.000)    2.754
   1.527   (  -4.838   -4.435   -0.000)    6.564
   1.590   (   3.054    2.093    0.000)    3.703
   1.645   (  -3.065   -1.720   -0.000)    3.515
   1.727   (   2.758   -4.242    0.000)    5.060
   1.783   (   4.451    5.424   -0.000)    7.017
   1.843   (   3.219    5.225   -0.000)    6.137
   1.884   (   2.207    0.365    0.000)    2.237
   1.954   (  -1.982   -0.693   -0.000)    2.099
   1.978   (  -3.092   -0.080   -0.000)    3.093
   2.023   (   3.120   -4.519    0.000)    5.491
   2.054   (  -0.488   -4.981    0.000)    5.005
   2.058   (  -2.246    0.578   -0.000)    2.319
   2.285   (  -0.421   -2.335    0.000)    2.373
   2.321   (   2.207   -0.193    0.000)    2.215
   2.336   (   2.846   -2.567    0.000)    3.833
   2.381   (   4.073   -2.143    0.000)    4.603
   2.382   (  -3.242    3.901   -0.000)    5.073
   2.450   (  -1.475   -1.061   -0.000)    1.817
   2.515   (  -1.418    0.962   -0.000)    1.713
   2.528   (   4.564    3.526    0.000)    5.767
   2.758   (  -2.297    0.785   -0.000)    2.427
   2.760   (  -2.604   -0.285   -0.000)    2.619
   3.120   (   1.086    0.988    0.000)    1.468
   3.156   (   0.502   -0.313    0.000)    0.591
   3.382   (   2.072   -9.167    0.000)    9.398
   3.433   (   3.592    7.631   -0.000)    8.434
   3.587   (  -2.382   -0.967   -0.000)    2.571
   3.650   (   2.965    3.271   -0.000)    4.415
   3.717   (  -5.489   -3.028   -0.000)    6.269
   3.767   (   5.821   -5.571    0.000)    8.057
   4.008   (  -0.343   -0.601    0.000)    0.692
   4.111   (   0.929    0.607    0.000)    1.110
   4.115   (  -3.754   -8.068    0.000)    8.898
   4.400   (  -6.461   11.311   -0.000)   13.027
   4.744   (  -0.825   -3.925    0.000)    4.011
   4.963   (   0.431    0.283    0.000)    0.516
======================= Grid point 33 (16/32) =======================
q-point: ( 0.08  0.58 -0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.608   (   2.676   -3.498    2.327)    4.981
   0.725   (  -5.559  -10.033    1.451)   11.561
   0.945   (   1.817    4.948    3.251)    6.193
   1.034   (   1.683   -1.248   -1.233)    2.431
   1.041   (   1.622   -1.157    3.324)    3.875
   1.128   (  -3.139   -1.242    0.892)    3.492
   1.225   (  -3.543   -3.804    3.958)    6.534
   1.253   (  -5.051    2.303    1.117)    5.663
   1.323   (  -0.683    0.536    2.849)    2.979
   1.426   (  -2.146   -0.004   -0.878)    2.319
   1.436   (   0.527   -2.644   -0.413)    2.727
   1.573   (  -2.148    7.038    0.136)    7.360
   1.616   (   0.961    2.677    0.653)    2.918
   1.632   (  -0.415    0.628    2.623)    2.729
   1.699   (  -2.431    4.007   -1.662)    4.973
   1.704   (   0.826    0.074   -0.070)    0.833
   1.733   (  -0.111   -1.952    3.712)    4.195
   1.807   (  -0.602    1.908   -3.768)    4.266
   1.871   (   2.284    0.751   -0.287)    2.421
   1.899   (   0.629    0.658   -0.084)    0.914
   1.963   (  -0.194   -2.007   -3.417)    3.968
   2.039   (  -1.976    4.951    1.945)    5.674
   2.092   (  -0.031    4.591   -2.619)    5.286
   2.134   (   0.420   -4.125    0.027)    4.146
   2.182   (  -2.074   -0.490   -1.283)    2.488
   2.257   (   2.879   -3.426   -3.174)    5.487
   2.274   (   2.429   -2.375   -3.310)    4.743
   2.295   (   0.307   -0.575   -2.828)    2.902
   2.441   (  -1.573   -1.332    1.374)    2.477
   2.458   (   3.532    1.539    2.167)    4.420
   2.524   (  -1.076    1.367    1.675)    2.415
   2.562   (  -0.858    0.644    0.318)    1.119
   2.716   (  -0.999   -2.597    2.369)    3.655
   2.894   (  -1.901    4.258    0.339)    4.675
   3.064   (   3.865   -2.391   -1.237)    4.710
   3.111   (   2.600   -2.650    0.064)    3.713
   3.230   (  -1.985   -4.241   -0.394)    4.699
   3.452   (   6.918  -11.618   -6.068)   14.820
   3.540   (   1.894    6.200    0.501)    6.502
   3.560   (   2.429   -3.126    0.443)    3.984
   3.687   (   0.143    0.116   -0.403)    0.443
   3.722   (  -0.463   -0.290   -0.658)    0.856
   4.017   (  -0.387   -7.526   -3.353)    8.248
   4.034   (  -0.904   -0.537   -0.936)    1.408
   4.130   (  -2.761   -2.039   -0.727)    3.508
   4.571   (   2.150   -5.204    1.453)    5.815
   4.817   (  -4.934    9.860   -0.291)   11.029
   4.957   (  -0.639    2.586    1.904)    3.274
======================= Grid point 37 (17/32) =======================
q-point: ( 0.17 -0.17 -0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 50
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.551   ( -12.884   -0.000   -0.000)   12.884
   0.690   ( -15.037   -0.000   -0.000)   15.037
   0.784   (   6.255    0.000    0.000)    6.255
   0.934   ( -16.303   -0.000   -0.000)   16.303
   1.059   (   1.812    0.000    0.000)    1.812
   1.088   (  -2.843   -0.000   -0.000)    2.843
   1.101   ( -12.263   -0.000   -0.000)   12.263
   1.300   (   4.357    0.000    0.000)    4.357
   1.346   (   3.405    0.000    0.000)    3.405
   1.376   (   0.915    0.000    0.000)    0.915
   1.472   (  -0.087   -0.000   -0.000)    0.087
   1.474   (   0.003    0.000    0.000)    0.003
   1.521   (  -7.883   -0.000   -0.000)    7.883
   1.567   (  -4.150   -0.000   -0.000)    4.150
   1.626   (   6.139    0.000    0.000)    6.139
   1.780   (   5.491    0.000    0.000)    5.491
   1.839   (   8.634    0.000    0.000)    8.634
   1.907   (   4.776    0.000    0.000)    4.776
   1.933   (   0.422    0.000    0.000)    0.422
   1.936   (   7.123    0.000    0.000)    7.123
   2.062   (   1.611    0.000    0.000)    1.611
   2.085   (   2.532    0.000    0.000)    2.532
   2.129   (   1.293    0.000    0.000)    1.293
   2.171   (   5.353    0.000    0.000)    5.353
   2.210   (  -0.252   -0.000   -0.000)    0.252
   2.268   (  -4.849   -0.000   -0.000)    4.849
   2.314   (  -2.758   -0.000   -0.000)    2.758
   2.338   (  -1.080   -0.000   -0.000)    1.080
   2.400   (  -0.410   -0.000   -0.000)    0.410
   2.483   (  -3.139   -0.000   -0.000)    3.139
   2.485   (   3.176    0.000    0.000)    3.176
   2.489   (   4.493    0.000    0.000)    4.493
   2.657   (  -6.571   -0.000   -0.000)    6.571
   2.913   (   3.448    0.000    0.000)    3.448
   3.056   (  -3.817   -0.000   -0.000)    3.817
   3.068   (  -0.001   -0.000   -0.000)    0.001
   3.271   (  -5.386   -0.000   -0.000)    5.386
   3.462   (   1.495    0.000    0.000)    1.495
   3.535   (  -0.946   -0.000   -0.000)    0.946
   3.674   (   0.192    0.000    0.000)    0.192
   3.703   (   4.396    0.000    0.000)    4.396
   3.762   (  -0.740   -0.000   -0.000)    0.740
   4.093   (  -0.885   -0.000   -0.000)    0.885
   4.127   (   2.434    0.000    0.000)    2.434
   4.235   ( -16.321   -0.000   -0.000)   16.321
   4.457   (  10.338    0.000    0.000)   10.338
   4.906   (   5.155    0.000    0.000)    5.155
   4.923   (  -0.345   -0.000   -0.000)    0.345
======================= Grid point 38 (18/32) =======================
q-point: ( 0.17  0.17 -0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 26
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.288   (  -0.000    0.000   11.080)   11.080
   0.288   (  -0.000    0.000   11.080)   11.080
   0.415   (  -0.000    0.000   17.419)   17.419
   0.728   (  -0.000    0.000    1.266)    1.266
   0.728   (  -0.000    0.000    1.266)    1.266
   1.101   (  -0.000    0.000    1.748)    1.748
   1.112   (   0.000   -0.000   -2.066)    2.066
   1.132   (  -0.000    0.000    3.439)    3.439
   1.132   (  -0.000    0.000    3.439)    3.439
   1.480   (  -0.000    0.000    0.035)    0.035
   1.480   (  -0.000    0.000    0.035)    0.035
   1.481   (   0.000   -0.000   -0.046)    0.046
   1.518   (  -0.000    0.000    4.363)    4.363
   1.689   (  -0.000    0.000    2.433)    2.433
   1.747   (   0.000   -0.000   -0.462)    0.462
   1.794   (  -0.000    0.000    2.196)    2.196
   1.794   (  -0.000    0.000    2.196)    2.196
   1.947   (   0.000   -0.000   -7.472)    7.472
   1.947   (   0.000   -0.000   -7.472)    7.472
   2.030   (   0.000   -0.000   -6.255)    6.255
   2.030   (   0.000   -0.000   -6.255)    6.255
   2.080   (   0.000   -0.000   -0.941)    0.941
   2.173   (   0.000   -0.000   -5.619)    5.619
   2.212   (   0.000   -0.000   -0.231)    0.231
   2.257   (   0.000   -0.000   -0.416)    0.416
   2.274   (   0.000   -0.000   -3.307)    3.307
   2.296   (  -0.000    0.000    5.219)    5.219
   2.296   (  -0.000    0.000    5.219)    5.219
   2.332   (  -0.000    0.000    0.208)    0.208
   2.389   (  -0.000    0.000    2.812)    2.812
   2.389   (  -0.000    0.000    2.812)    2.812
   2.448   (  -0.000    0.000    4.149)    4.149
   2.594   (   0.000   -0.000   -2.433)    2.433
   2.880   (   0.000   -0.000   -0.895)    0.895
   2.961   (  -0.000    0.000    0.493)    0.493
   2.961   (  -0.000    0.000    0.493)    0.493
   2.996   (  -0.000    0.000    0.289)    0.289
   3.294   (  -0.000    0.000    0.517)    0.517
   3.545   (   0.000   -0.000   -1.303)    1.303
   3.558   (  -0.000    0.000    0.450)    0.450
   3.558   (  -0.000    0.000    0.450)    0.450
   3.782   (   0.000   -0.000   -2.060)    2.060
   3.782   (   0.000   -0.000   -2.060)    2.060
   4.081   (  -0.000    0.000    0.053)    0.053
   4.188   (  -0.000    0.000    0.549)    0.549
   4.188   (  -0.000    0.000    0.549)    0.549
   5.006   (  -0.000    0.000    1.032)    1.032
   5.254   (  -0.000    0.000    0.193)    0.193
======================= Grid point 39 (19/32) =======================
q-point: ( 0.17  0.50 -0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 67
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.584   (   2.014    0.000    4.843)    5.245
   0.789   (  11.683    0.000    2.324)   11.912
   0.918   (  -6.726   -0.000    4.494)    8.089
   1.011   (   0.321    0.000    0.209)    0.383
   1.014   (   6.503    0.000   -2.251)    6.881
   1.189   (  -3.797   -0.000    5.035)    6.306
   1.244   (  12.995    0.000    1.886)   13.131
   1.328   (  -2.196   -0.000    0.244)    2.210
   1.345   (   0.185    0.000    2.259)    2.266
   1.406   (  -2.698   -0.000   -0.966)    2.866
   1.466   (   0.196   -0.000   -1.657)    1.668
   1.543   (   4.005    0.000    6.292)    7.458
   1.611   (  -1.380   -0.000    0.452)    1.452
   1.659   (  -7.060   -0.000    0.391)    7.070
   1.692   (   1.014    0.000    1.841)    2.102
   1.730   (  -1.450   -0.000   -0.885)    1.699
   1.750   (  -4.500   -0.000   -3.130)    5.481
   1.792   (  -5.709   -0.000   -3.568)    6.733
   1.849   (   3.103    0.000    2.524)    4.000
   1.901   (  -3.246   -0.000    0.388)    3.269
   1.990   (   3.795    0.000   -3.325)    5.046
   2.011   (  -7.391   -0.000   -2.861)    7.925
   2.111   (  -4.476   -0.000    0.107)    4.478
   2.127   (   6.931    0.000   -0.107)    6.932
   2.210   (  -1.616   -0.000   -1.569)    2.252
   2.228   (   5.661    0.000   -3.837)    6.839
   2.244   (   1.986    0.000   -1.905)    2.752
   2.300   (  -0.821   -0.000   -3.559)    3.652
   2.416   (  -1.267   -0.000    1.387)    1.879
   2.459   (   0.287    0.000    1.901)    1.923
   2.529   (   0.100    0.000    0.846)    0.852
   2.573   (  -0.934   -0.000    1.001)    1.369
   2.726   (   2.980    0.000    3.130)    4.322
   2.913   (  -4.480   -0.000    0.306)    4.490
   3.015   (   4.432    0.000   -1.962)    4.847
   3.087   (   1.754    0.000    0.361)    1.791
   3.262   (   3.875    0.000   -0.185)    3.879
   3.354   (  15.331    0.000   -5.388)   16.251
   3.523   (  -4.024   -0.000   -0.467)    4.051
   3.536   (   2.331    0.000    0.051)    2.331
   3.670   (  -2.834   -0.000   -2.558)    3.818
   3.731   (   2.516    0.000    2.130)    3.297
   4.029   (  10.036    0.000   -3.943)   10.782
   4.045   (   0.875    0.000   -0.606)    1.064
   4.170   (   0.615    0.000   -1.120)    1.278
   4.535   (   9.752    0.000    0.126)    9.753
   4.875   ( -13.060   -0.000    1.334)   13.128
   4.963   (  -3.092   -0.000    1.651)    3.505
======================= Grid point 44 (20/32) =======================
q-point: ( 0.17  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 64
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.576   (  -2.291   -2.291    2.333)    3.992
   0.744   (   8.282    8.282   -1.861)   11.860
   0.881   (  -7.013   -7.013   -0.293)    9.922
   0.982   (  -1.756   -1.756    1.328)    2.817
   1.048   (   3.037    3.037    1.515)    4.554
   1.054   (  -0.647   -0.647    3.267)    3.392
   1.121   (   2.244    2.244    0.367)    3.194
   1.181   (  -1.961   -1.961   -0.214)    2.782
   1.242   (  -1.955   -1.955    1.247)    3.033
   1.319   (  -2.615   -2.615    0.011)    3.699
   1.489   (  -6.307   -6.307    1.240)    9.005
   1.531   (   5.940    5.940    1.581)    8.548
   1.539   (  -7.228   -7.228   -5.496)   11.606
   1.608   (   1.479    1.479    5.641)    6.016
   1.634   (   6.775    6.775    2.521)    9.907
   1.679   (   5.160    5.160    2.311)    7.655
   1.833   (  -1.700   -1.700   -3.351)    4.124
   1.837   (  -4.055   -4.055    2.887)    6.421
   1.910   (  -1.741   -1.741    0.372)    2.491
   1.913   (  -1.307   -1.307   -4.549)    4.910
   1.971   (  -0.836   -0.836   -4.440)    4.595
   2.034   (   0.053    0.053    0.075)    0.106
   2.095   (  -2.320   -2.320   -1.129)    3.470
   2.107   (  -2.015   -2.015    0.751)    2.947
   2.214   (   2.363    2.363   -1.899)    3.843
   2.336   (   0.132    0.132   -0.659)    0.685
   2.378   (  -1.712   -1.712   -1.964)    3.118
   2.405   (   2.350    2.350   -0.991)    3.468
   2.416   (   2.243    2.243    0.555)    3.220
   2.424   (  -0.314   -0.314    2.226)    2.269
   2.435   (   1.539    1.539    2.383)    3.227
   2.543   (  -1.073   -1.073    0.760)    1.697
   2.690   (   1.890    1.890    0.075)    2.674
   2.833   (  -3.487   -3.487    0.005)    4.931
   3.105   (   1.316    1.316    0.616)    1.961
   3.174   (  -0.925   -0.925   -0.027)    1.308
   3.265   (   6.238    6.238   -1.367)    8.927
   3.551   (  -3.822   -3.822    0.381)    5.419
   3.586   (   0.684    0.684    1.137)    1.493
   3.608   (   1.752    1.752    1.304)    2.800
   3.706   (   0.327    0.327   -5.531)    5.551
   3.799   (  -0.052   -0.052   -2.662)    2.663
   4.050   (   0.277    0.277   -2.157)    2.192
   4.050   (  -0.549   -0.549    2.008)    2.153
   4.173   (   3.960    3.960   -4.363)    7.099
   4.384   (  -7.612   -7.612    6.862)   12.766
   4.840   (  -0.885   -0.885    0.903)    1.544
   4.907   (  -0.667   -0.667   -3.046)    3.189
======================= Grid point 46 (21/32) =======================
q-point: ( 0.17 -0.17  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.706   (   4.375   -6.020    0.000)    7.442
   0.836   (   3.807    9.791   -0.000)   10.505
   0.864   (  -8.495    7.673   -0.000)   11.448
   0.980   (  -6.308   -5.645   -0.000)    8.465
   1.082   (   4.576   -0.203    0.000)    4.581
   1.145   (  -3.486   -2.862   -0.000)    4.510
   1.196   (  -0.018    4.262   -0.000)    4.262
   1.214   (  -4.143    5.369   -0.000)    6.782
   1.313   (  -9.784   12.804   -0.000)   16.114
   1.459   (  -3.351   -4.221    0.000)    5.390
   1.467   (   1.000   -3.801    0.000)    3.930
   1.503   (   4.946    0.933    0.000)    5.034
   1.570   (  -6.217    4.650   -0.000)    7.764
   1.619   (   5.225    4.838    0.000)    7.121
   1.620   (  -3.941   -3.786   -0.000)    5.465
   1.723   (   0.715    5.863   -0.000)    5.906
   1.730   (   0.895   -6.248    0.000)    6.312
   1.777   (  -0.719    3.782   -0.000)    3.850
   1.831   (   2.419   -1.217    0.000)    2.708
   1.960   (   4.241   -4.883    0.000)    6.467
   1.997   (  -1.972   -2.746    0.000)    3.381
   2.026   (   1.401   -4.226    0.000)    4.452
   2.049   (  -3.018   -0.321   -0.000)    3.035
   2.107   (  -2.009    1.223   -0.000)    2.352
   2.272   (   3.358   -1.582    0.000)    3.712
   2.274   (   1.595   -2.372    0.000)    2.858
   2.312   (  -0.367   -7.250    0.000)    7.259
   2.327   (   0.232   -3.532    0.000)    3.539
   2.428   (   1.317    4.030   -0.000)    4.240
   2.452   (   0.033    5.739   -0.000)    5.739
   2.469   (  -0.590    0.852   -0.000)    1.036
   2.530   (  -0.074    0.540   -0.000)    0.545
   2.790   (  -1.119    6.068   -0.000)    6.171
   2.815   (  -2.042   -3.923    0.000)    4.423
   3.078   (   2.453   -0.744    0.000)    2.564
   3.128   (   1.935    0.408    0.000)    1.978
   3.368   (  -0.229   -7.851    0.000)    7.855
   3.401   (  -0.069    7.491   -0.000)    7.491
   3.595   (   2.054    3.151   -0.000)    3.761
   3.596   (   1.049   -6.347    0.000)    6.433
   3.657   (   4.267   -5.095    0.000)    6.646
   3.786   (  -1.125   -1.081   -0.000)    1.560
   4.040   (  -2.515   -1.269   -0.000)    2.817
   4.108   (  -0.713    0.894   -0.000)    1.143
   4.189   (  -2.937  -11.531    0.000)   11.900
   4.524   (  -4.849   15.486   -0.000)   16.228
   4.756   (  -0.281   -7.475    0.000)    7.481
   4.951   (   0.682    0.938   -0.000)    1.160
======================= Grid point 47 (22/32) =======================
q-point: ( 0.17  0.17  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 39
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.403   (  -0.000   -0.000    2.375)    2.375
   0.454   (  -0.000   -0.000   -1.793)    1.793
   0.987   (  -0.000   -0.000    0.343)    0.343
   0.993   (  -0.000   -0.000   -0.104)    0.104
   1.130   (   0.000    0.000    4.463)    4.463
   1.132   (  -0.000   -0.000   -0.157)    0.157
   1.139   (  -0.680   -0.680   -0.663)    1.168
   1.139   (   0.680    0.680    3.026)    3.175
   1.249   (  -0.000   -0.000   -5.844)    5.844
   1.262   (   0.000    0.000    3.578)    3.578
   1.273   (  -0.000   -0.000   -8.078)    8.078
   1.343   (   0.000    0.000   -3.433)    3.433
   1.665   (  -0.000   -0.000    4.613)    4.613
   1.686   (  -0.000   -0.000    2.965)    2.965
   1.714   (  -0.000   -0.000   10.088)   10.088
   1.738   (   0.000    0.000   -1.659)    1.659
   1.778   (  -0.000   -0.000   -4.906)    4.906
   1.892   (  -0.000   -0.000    0.379)    0.379
   1.914   (   0.000    0.000   -1.179)    1.179
   1.925   (  -0.000   -0.000   -8.171)    8.171
   1.927   (   0.000    0.000    3.617)    3.617
   1.995   (  -0.000   -0.000   -1.579)    1.579
   2.021   (  -0.000   -0.000    3.136)    3.136
   2.068   (  -0.000   -0.000   -1.069)    1.069
   2.261   (  -0.000   -0.000   -0.909)    0.909
   2.276   (  -0.000   -0.000   -2.298)    2.298
   2.418   (   0.000    0.000    0.595)    0.595
   2.431   (  -0.000   -0.000    2.883)    2.883
   2.435   (  -0.000   -0.000   -0.792)    0.792
   2.450   (  -0.000   -0.000    4.886)    4.886
   2.476   (  -0.000   -0.000   -1.050)    1.050
   2.549   (  -0.000   -0.000   -3.408)    3.408
   2.719   (  -0.000   -0.000    0.473)    0.473
   2.731   (  -0.000   -0.000   -0.445)    0.445
   3.135   (  -0.000   -0.000    0.671)    0.671
   3.151   (  -0.000   -0.000   -0.588)    0.588
   3.444   (   0.000    0.000    1.630)    1.630
   3.490   (  -0.000   -0.000   -2.131)    2.131
   3.547   (  -0.000   -0.000    1.798)    1.798
   3.594   (  -0.000   -0.000   -2.113)    2.113
   3.726   (  -0.000   -0.000    2.918)    2.918
   3.798   (  -0.000   -0.000   -3.130)    3.130
   4.048   (  -0.000   -0.000    2.595)    2.595
   4.102   (  -0.000   -0.000   -1.919)    1.919
   4.132   (  -0.000   -0.000    4.915)    4.915
   4.253   (  -0.000   -0.000   -5.006)    5.006
   4.799   (  -0.000   -0.000    4.996)    4.996
   4.914   (  -0.000   -0.000   -4.396)    4.396
======================= Grid point 50 (23/32) =======================
q-point: ( 0.17  0.83 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 44
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.620   (  -3.649    3.649   -0.000)    5.160
   0.692   (   5.228   -5.228    0.000)    7.394
   0.912   ( -13.068   13.068   -0.000)   18.481
   1.008   (  -0.912    0.912   -0.000)    1.290
   1.074   (   2.000   -2.000    0.000)    2.829
   1.168   (   2.299   -2.299    0.000)    3.251
   1.209   (  -1.780    1.780   -0.000)    2.517
   1.216   (  -0.489    0.489   -0.000)    0.691
   1.369   ( -11.637   11.637   -0.000)   16.457
   1.378   (  -1.399    1.399   -0.000)    1.978
   1.456   (  -0.111    0.111   -0.000)    0.156
   1.481   (   0.365   -0.365    0.000)    0.517
   1.512   (   5.795   -5.795    0.000)    8.196
   1.590   (  -7.714    7.714   -0.000)   10.910
   1.654   (   0.658   -0.658    0.000)    0.931
   1.685   (  -1.661    1.661   -0.000)    2.349
   1.822   (   0.965   -0.965    0.000)    1.365
   1.876   (   1.009   -1.009    0.000)    1.427
   1.913   (   5.710   -5.710    0.000)    8.075
   1.914   (   0.718   -0.718    0.000)    1.015
   1.958   (   2.019   -2.019    0.000)    2.855
   1.999   (  -3.198    3.198   -0.000)    4.522
   2.046   (  -5.290    5.290   -0.000)    7.481
   2.069   (  -0.978    0.978   -0.000)    1.383
   2.208   (   3.741   -3.741    0.000)    5.291
   2.270   (   3.392   -3.392    0.000)    4.797
   2.293   (   3.199   -3.199    0.000)    4.524
   2.320   (   4.230   -4.230    0.000)    5.982
   2.424   (  -0.551    0.551   -0.000)    0.779
   2.486   (  -4.654    4.654   -0.000)    6.582
   2.542   (  -0.913    0.913   -0.000)    1.291
   2.559   (   1.374   -1.374    0.000)    1.944
   2.716   (  -0.294    0.294   -0.000)    0.415
   2.842   (  -3.648    3.648   -0.000)    5.159
   3.117   (   1.309   -1.309    0.000)    1.852
   3.178   (  -1.263    1.263   -0.000)    1.786
   3.209   (   6.661   -6.661    0.000)    9.420
   3.498   (   4.496   -4.496    0.000)    6.358
   3.588   (  -2.919    2.919   -0.000)    4.128
   3.650   (   5.304   -5.304    0.000)    7.501
   3.671   (   0.401   -0.401    0.000)    0.568
   3.692   (  -0.200    0.200   -0.000)    0.283
   3.987   (   3.726   -3.726    0.000)    5.270
   4.005   (  -0.506    0.506   -0.000)    0.715
   4.112   (   0.190   -0.190    0.000)    0.268
   4.635   ( -11.073   11.073   -0.000)   15.659
   4.660   (   4.141   -4.141    0.000)    5.856
   4.973   (  -0.528    0.528   -0.000)    0.747
======================= Grid point 56 (24/32) =======================
q-point: ( 0.17  0.50 -0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 66
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.566   (   3.863   -3.863    2.391)    5.963
   0.587   (  -0.068    0.068    8.348)    8.348
   0.912   (   8.112   -8.112   -0.486)   11.483
   0.925   (   3.537   -3.537    6.609)    8.288
   1.051   (  -0.034    0.034    4.348)    4.349
   1.056   (   0.116   -0.116    0.758)    0.775
   1.238   (  -1.064    1.064    4.320)    4.575
   1.302   (   1.732   -1.732    2.893)    3.791
   1.353   (  -5.568    5.568    0.473)    7.889
   1.446   (  -1.068    1.068   -2.695)    3.089
   1.452   (  -1.890    1.890   -0.118)    2.676
   1.532   (   4.445   -4.445    5.641)    8.446
   1.665   (  -1.350    1.350    4.022)    4.452
   1.678   (   1.941   -1.941    5.371)    6.031
   1.720   (  -1.744    1.744   -8.912)    9.247
   1.768   (  -1.957    1.957   -1.135)    2.991
   1.770   (  -5.573    5.573    0.073)    7.882
   1.827   (   4.748   -4.748   -4.050)    7.842
   1.878   (  -1.339    1.339   -3.570)    4.041
   1.891   (   0.353   -0.353   -2.272)    2.327
   1.933   (  -1.176    1.176   -1.878)    2.508
   1.997   (   6.001   -6.001   -1.774)    8.671
   2.140   (  -2.960    2.960   -0.669)    4.239
   2.170   (   4.081   -4.081   -6.772)    8.898
   2.179   (  -4.050    4.050   -1.131)    5.839
   2.224   (  -0.693    0.693    0.114)    0.987
   2.232   (   1.197   -1.197   -2.360)    2.904
   2.321   (  -2.852    2.852   -2.017)    4.509
   2.415   (   0.212   -0.212    2.205)    2.225
   2.466   (   1.197   -1.197    2.195)    2.771
   2.520   (   2.486   -2.486    1.702)    3.906
   2.537   (   0.258   -0.258    0.555)    0.664
   2.705   (   0.691   -0.691    3.814)    3.937
   2.957   (   7.428   -7.428   -4.646)   11.486
   2.960   (  -2.199    2.199    0.485)    3.148
   3.040   (   2.230   -2.230    0.244)    3.162
   3.190   (  11.305  -11.305   -6.011)   17.080
   3.226   (  -1.773    1.773    0.273)    2.523
   3.518   (   1.048   -1.048   -0.340)    1.520
   3.619   (   1.191   -1.191    0.503)    1.758
   3.632   (   0.670   -0.670    0.866)    1.283
   3.739   (   0.295   -0.295   -2.597)    2.630
   3.861   (   5.106   -5.106   -3.968)    8.240
   4.041   (  -0.705    0.705   -0.526)    1.127
   4.144   (  -2.221    2.221    0.701)    3.218
   4.417   (   7.622   -7.622   -1.010)   10.826
   4.998   (  -1.539    1.539    1.158)    2.465
   5.035   (  -9.060    9.060    2.028)   12.972
======================= Grid point 57 (25/32) =======================
q-point: ( 0.25 -0.25 -0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 42
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.664   (   4.041    0.000    0.000)    4.041
   0.794   ( -10.757   -0.000   -0.000)   10.757
   0.949   (  -9.094   -0.000   -0.000)    9.094
   1.026   (   1.088    0.000    0.000)    1.088
   1.029   (   7.621    0.000    0.000)    7.621
   1.195   (  -6.211   -0.000   -0.000)    6.211
   1.213   (  -0.612   -0.000   -0.000)    0.612
   1.257   (  -4.848   -0.000   -0.000)    4.848
   1.418   (  -5.548   -0.000   -0.000)    5.548
   1.453   (  -1.472   -0.000   -0.000)    1.472
   1.472   (   2.160    0.000    0.000)    2.160
   1.495   (   6.173    0.000    0.000)    6.173
   1.537   ( -10.222   -0.000   -0.000)   10.222
   1.623   (  -4.374   -0.000   -0.000)    4.374
   1.682   (   4.491    0.000    0.000)    4.491
   1.703   (   3.607    0.000    0.000)    3.607
   1.732   (  -8.300   -0.000   -0.000)    8.300
   1.808   (   4.637    0.000    0.000)    4.637
   1.809   (   0.738    0.000    0.000)    0.738
   1.885   (   3.187    0.000    0.000)    3.187
   2.018   (   4.100    0.000    0.000)    4.100
   2.021   (   4.575    0.000    0.000)    4.575
   2.078   (   1.076    0.000    0.000)    1.076
   2.128   (  -1.291   -0.000   -0.000)    1.291
   2.229   (  -1.214   -0.000   -0.000)    1.214
   2.258   (   2.379    0.000    0.000)    2.379
   2.309   (   6.702    0.000    0.000)    6.702
   2.333   (   5.887    0.000    0.000)    5.887
   2.402   (  -0.248   -0.000   -0.000)    0.248
   2.467   (  -7.831   -0.000   -0.000)    7.831
   2.476   (  -0.855   -0.000   -0.000)    0.855
   2.527   (  -1.281   -0.000   -0.000)    1.281
   2.802   (  -6.731   -0.000   -0.000)    6.731
   2.836   (   3.989    0.000    0.000)    3.989
   3.049   (   1.585    0.000    0.000)    1.585
   3.103   (  -0.867   -0.000   -0.000)    0.867
   3.372   (   6.207    0.000    0.000)    6.207
   3.407   (  -6.917   -0.000   -0.000)    6.917
   3.574   (  -2.712   -0.000   -0.000)    2.712
   3.580   (   9.589    0.000    0.000)    9.589
   3.611   (   4.316    0.000    0.000)    4.316
   3.790   (  -0.746   -0.000   -0.000)    0.746
   4.075   (   2.460    0.000    0.000)    2.460
   4.120   (  -1.683   -0.000   -0.000)    1.683
   4.221   (  12.468    0.000    0.000)   12.468
   4.575   ( -16.627   -0.000   -0.000)   16.627
   4.758   (   9.477    0.000    0.000)    9.477
   4.943   (  -1.590   -0.000   -0.000)    1.590
======================= Grid point 59 (26/32) =======================
q-point: ( 0.25  0.42 -0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.564   (   9.471    0.000    3.149)    9.980
   0.570   (   1.358    0.000   10.624)   10.710
   0.815   (   8.983    0.000    2.656)    9.367
   0.930   (   9.172    0.000    4.888)   10.393
   1.035   (  -4.300   -0.000    3.597)    5.606
   1.075   (  -0.191   -0.000   -0.669)    0.696
   1.198   (   4.481    0.000    4.423)    6.297
   1.321   (   2.505    0.000    3.278)    4.125
   1.391   (   3.762    0.000    4.119)    5.579
   1.434   (   2.955    0.000    0.977)    3.112
   1.455   (   1.078    0.000   -1.718)    2.029
   1.484   (   0.254   -0.000   -3.570)    3.579
   1.649   (   0.352    0.000    4.717)    4.730
   1.675   (  -2.293   -0.000   -1.126)    2.554
   1.713   (   5.074    0.000   -3.401)    6.108
   1.791   (   1.588    0.000   -3.201)    3.573
   1.795   (  -5.043   -0.000    0.571)    5.075
   1.831   (   4.683    0.000   -5.689)    7.369
   1.849   (  -4.106   -0.000   -0.621)    4.153
   1.886   (  -0.539   -0.000   -4.623)    4.654
   1.963   (   6.910    0.000   -1.276)    7.026
   1.965   (  -2.671   -0.000   -1.714)    3.173
   2.123   (   1.928    0.000   -5.354)    5.690
   2.151   (  -4.287   -0.000   -3.339)    5.434
   2.210   (   0.050   -0.000   -0.672)    0.674
   2.241   (   0.177   -0.000   -0.785)    0.805
   2.244   (  -5.961   -0.000    0.495)    5.981
   2.343   (  -3.686   -0.000   -0.391)    3.707
   2.430   (   2.276    0.000    2.649)    3.492
   2.434   (  -0.390    0.000    1.421)    1.474
   2.504   (   5.483    0.000    2.115)    5.877
   2.530   (  -0.110   -0.000   -0.501)    0.513
   2.667   (   7.940    0.000    2.809)    8.422
   2.880   (   5.124    0.000   -4.305)    6.692
   2.982   (  -2.292   -0.000    0.484)    2.342
   3.021   (   3.835    0.000    0.144)    3.838
   3.064   (   9.194    0.000   -2.268)    9.470
   3.256   (  -2.653   -0.000    0.272)    2.667
   3.510   (   0.260    0.000   -0.764)    0.807
   3.572   (  -0.672   -0.000    0.154)    0.690
   3.645   (   5.236    0.000    1.099)    5.350
   3.716   (  -1.442   -0.000   -2.373)    2.777
   3.824   (   8.407    0.000   -3.346)    9.048
   4.050   (  -1.443   -0.000   -0.321)    1.479
   4.188   (  -1.083   -0.000    0.345)    1.137
   4.318   (   9.660    0.000   -0.438)    9.670
   5.013   (  -1.730   -0.000    1.050)    2.024
   5.136   ( -10.930   -0.000    1.582)   11.044
======================= Grid point 69 (27/32) =======================
q-point: (-0.75 -0.25  0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.646   (  -0.393    0.738   -0.000)    0.836
   0.738   (  -6.207   -9.772    0.000)   11.576
   0.989   (   1.720    4.806   -0.000)    5.105
   1.012   (   2.894   -2.836    0.000)    4.052
   1.087   (   4.565    0.983    0.000)    4.670
   1.151   (  -8.427    1.718   -0.000)    8.600
   1.258   (  -3.711   -1.022   -0.000)    3.849
   1.258   (  -2.915    1.180   -0.000)    3.144
   1.389   (   2.808    0.109    0.000)    2.810
   1.406   (  -0.430   -0.752    0.000)    0.866
   1.480   (   2.637    0.368    0.000)    2.663
   1.512   (  -4.641   -4.149   -0.000)    6.225
   1.575   (  -6.971   11.586   -0.000)   13.521
   1.664   (  -1.302    1.671   -0.000)    2.119
   1.701   (  -2.019    5.368   -0.000)    5.736
   1.706   (  -0.626    3.126   -0.000)    3.188
   1.749   (   4.055   -1.143    0.000)    4.213
   1.833   (   3.519   -5.539    0.000)    6.562
   1.868   (  -0.842   -1.775    0.000)    1.964
   1.900   (  -1.700    2.716   -0.000)    3.204
   1.932   (   1.393   -3.972    0.000)    4.209
   2.022   (   3.414    6.568   -0.000)    7.402
   2.083   (  -3.448    3.697   -0.000)    5.055
   2.152   (   2.070   -7.335    0.000)    7.622
   2.153   (  -3.117    2.321   -0.000)    3.886
   2.210   (   3.764   -5.790    0.000)    6.906
   2.222   (   0.968   -2.258    0.000)    2.457
   2.253   (   1.300   -1.470    0.000)    1.962
   2.460   (  -2.089   -1.273   -0.000)    2.447
   2.483   (   3.961    1.813    0.000)    4.356
   2.546   (   0.395    0.616   -0.000)    0.732
   2.563   (  -2.135    2.701   -0.000)    3.443
   2.750   (  -1.949   -1.214   -0.000)    2.296
   2.899   (  -1.897    4.414   -0.000)    4.804
   3.042   (   5.124   -3.568    0.000)    6.244
   3.112   (   1.950   -3.419    0.000)    3.936
   3.225   (  -1.863   -3.791    0.000)    4.224
   3.380   (   6.158  -13.069    0.000)   14.448
   3.554   (   2.818    5.862   -0.000)    6.504
   3.566   (   2.863   -3.664    0.000)    4.650
   3.661   (   2.229    3.038   -0.000)    3.768
   3.735   (  -3.167   -3.874    0.000)    5.004
   3.964   (  -0.474  -10.034    0.000)   10.046
   4.025   (  -1.108   -0.153   -0.000)    1.118
   4.126   (  -1.845    0.889   -0.000)    2.048
   4.579   (   3.659   -9.282    0.000)    9.978
   4.823   (  -6.904   13.662   -0.000)   15.307
   4.981   (  -0.295    1.777   -0.000)    1.801
======================= Grid point 77 (28/32) =======================
q-point: ( 0.33  0.67 -0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 50
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.663   (  -4.169   -0.000   -0.000)    4.169
   0.816   (  11.538    0.000    0.000)   11.538
   0.966   (   7.450    0.000    0.000)    7.450
   0.968   (  -5.817   -0.000   -0.000)    5.817
   1.030   (  -1.795   -0.000   -0.000)    1.795
   1.269   (  11.204    0.000    0.000)   11.204
   1.297   (  -6.499   -0.000   -0.000)    6.499
   1.297   (  -3.260   -0.000   -0.000)    3.260
   1.384   (  -0.395   -0.000   -0.000)    0.395
   1.397   (   3.264    0.000    0.000)    3.264
   1.417   (  -5.254   -0.000   -0.000)    5.254
   1.592   (   4.133    0.000    0.000)    4.133
   1.637   (  -8.270   -0.000   -0.000)    8.270
   1.676   (   2.018    0.000    0.000)    2.018
   1.687   (  -1.320   -0.000   -0.000)    1.320
   1.705   (  -7.352   -0.000   -0.000)    7.352
   1.734   (   2.305    0.000    0.000)    2.305
   1.773   (  -0.345   -0.000   -0.000)    0.345
   1.907   (   5.577    0.000    0.000)    5.577
   1.917   (   5.865    0.000    0.000)    5.865
   1.939   (  -6.712   -0.000   -0.000)    6.712
   1.955   ( -10.053   -0.000   -0.000)   10.053
   2.113   (  -0.572   -0.000   -0.000)    0.572
   2.126   (   9.642    0.000    0.000)    9.642
   2.171   (   7.920    0.000    0.000)    7.920
   2.189   (  -4.545   -0.000   -0.000)    4.545
   2.218   (   1.359    0.000    0.000)    1.359
   2.249   (  -0.847   -0.000   -0.000)    0.847
   2.433   (  -1.465   -0.000   -0.000)    1.465
   2.485   (   0.256    0.000    0.000)    0.256
   2.539   (   0.049    0.000    0.000)    0.049
   2.585   (  -2.048   -0.000   -0.000)    2.048
   2.772   (   0.989    0.000    0.000)    0.989
   2.918   (  -4.654   -0.000   -0.000)    4.654
   2.982   (   6.561    0.000    0.000)    6.561
   3.091   (   2.229    0.000    0.000)    2.229
   3.260   (   3.451    0.000    0.000)    3.451
   3.290   (  17.107    0.000    0.000)   17.107
   3.526   (  -4.382   -0.000   -0.000)    4.382
   3.537   (   2.847    0.000    0.000)    2.847
   3.635   (  -3.015   -0.000   -0.000)    3.015
   3.764   (   3.406    0.000    0.000)    3.406
   3.972   (  11.406    0.000    0.000)   11.406
   4.038   (   0.912    0.000    0.000)    0.912
   4.158   (  -1.934   -0.000   -0.000)    1.934
   4.534   (  11.825    0.000    0.000)   11.825
   4.896   ( -14.710   -0.000   -0.000)   14.710
   4.984   (  -2.316   -0.000   -0.000)    2.316
======================= Grid point 80 (29/32) =======================
q-point: ( 0.33  0.33 -0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 32
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.456   (  -0.000    0.000    3.664)    3.664
   0.456   (  -0.000    0.000    3.664)    3.664
   0.753   (  -0.000    0.000   11.844)   11.844
   0.811   (  -0.000    0.000    5.894)    5.894
   0.811   (  -0.000    0.000    5.894)    5.894
   1.068   (   0.000   -0.000   -1.544)    1.544
   1.139   (  -0.000    0.000    1.395)    1.395
   1.256   (  -0.000    0.000    6.647)    6.647
   1.256   (  -0.000    0.000    6.647)    6.647
   1.441   (   0.000   -0.000   -2.599)    2.599
   1.441   (   0.000   -0.000   -2.599)    2.599
   1.480   (   0.000   -0.000   -0.046)    0.046
   1.630   (  -0.000    0.000    4.861)    4.861
   1.732   (   0.000   -0.000   -0.855)    0.855
   1.755   (   0.000   -0.000   -9.173)    9.173
   1.755   (   0.000   -0.000   -9.173)    9.173
   1.758   (  -0.000    0.000    3.261)    3.261
   1.853   (  -0.000    0.000    1.204)    1.204
   1.853   (  -0.000    0.000    1.204)    1.204
   1.918   (   0.000   -0.000   -2.173)    2.173
   1.918   (   0.000   -0.000   -2.173)    2.173
   1.985   (   0.000   -0.000   -5.453)    5.453
   2.129   (  -0.000    0.000    0.226)    0.226
   2.181   (   0.000   -0.000   -4.373)    4.373
   2.206   (   0.000   -0.000   -0.214)    0.214
   2.242   (   0.000   -0.000   -0.861)    0.861
   2.337   (  -0.000    0.000    0.190)    0.190
   2.389   (  -0.000    0.000    2.918)    2.918
   2.389   (  -0.000    0.000    2.918)    2.918
   2.437   (  -0.000    0.000    1.318)    1.318
   2.437   (  -0.000    0.000    1.318)    1.318
   2.520   (  -0.000    0.000    2.199)    2.199
   2.531   (   0.000   -0.000   -2.597)    2.597
   2.857   (   0.000   -0.000   -0.980)    0.980
   2.975   (  -0.000    0.000    0.566)    0.566
   2.975   (  -0.000    0.000    0.566)    0.566
   3.006   (  -0.000    0.000    0.486)    0.486
   3.305   (  -0.000    0.000    0.360)    0.360
   3.513   (   0.000   -0.000   -1.295)    1.295
   3.572   (  -0.000    0.000    0.719)    0.719
   3.572   (  -0.000    0.000    0.719)    0.719
   3.730   (   0.000   -0.000   -2.254)    2.254
   3.730   (   0.000   -0.000   -2.254)    2.254
   4.083   (  -0.000    0.000    0.051)    0.051
   4.200   (  -0.000    0.000    0.470)    0.470
   4.200   (  -0.000    0.000    0.470)    0.470
   5.031   (  -0.000    0.000    1.007)    1.007
   5.259   (  -0.000    0.000    0.193)    0.193
======================= Grid point 95 (30/32) =======================
q-point: ( 0.33  0.67 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 42
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.589   (   3.462   -3.462    0.000)    4.896
   0.723   (  -5.582    5.582   -0.000)    7.895
   0.868   (   9.096   -9.096    0.000)   12.864
   1.037   (  -0.521    0.521   -0.000)    0.737
   1.054   (   5.114   -5.114    0.000)    7.232
   1.103   (   0.006   -0.006    0.000)    0.009
   1.285   (  -1.506    1.506   -0.000)    2.130
   1.358   (  -6.041    6.041   -0.000)    8.543
   1.394   (   1.528   -1.528    0.000)    2.161
   1.399   (  -0.495    0.495   -0.000)    0.700
   1.471   (  -0.636    0.636   -0.000)    0.899
   1.536   (  -1.850    1.850   -0.000)    2.616
   1.674   (   5.532   -5.532    0.000)    7.823
   1.692   (  -1.227    1.227   -0.000)    1.735
   1.728   (   5.629   -5.629    0.000)    7.961
   1.756   (   2.121   -2.121    0.000)    2.999
   1.759   (  -0.071    0.071   -0.000)    0.100
   1.773   (  -6.519    6.519   -0.000)    9.219
   1.814   (   0.174   -0.174    0.000)    0.247
   1.872   (   1.048   -1.048    0.000)    1.483
   1.935   (  -0.332    0.332   -0.000)    0.470
   1.996   (   7.533   -7.533    0.000)   10.654
   2.064   (   0.996   -0.996    0.000)    1.408
   2.107   (  -1.238    1.238   -0.000)    1.751
   2.187   (  -5.174    5.174   -0.000)    7.317
   2.195   (   0.217   -0.217    0.000)    0.307
   2.219   (   0.060   -0.060    0.000)    0.085
   2.295   (  -4.651    4.651   -0.000)    6.577
   2.441   (  -0.106    0.106   -0.000)    0.150
   2.492   (   1.675   -1.675    0.000)    2.369
   2.543   (   2.093   -2.093    0.000)    2.960
   2.544   (   0.858   -0.858    0.000)    1.213
   2.781   (  -2.622    2.622   -0.000)    3.708
   2.863   (  11.999  -11.999    0.000)   16.969
   2.967   (  -2.214    2.214   -0.000)    3.131
   3.043   (   2.139   -2.139    0.000)    3.025
   3.135   (   8.227   -8.227    0.000)   11.635
   3.229   (  -1.763    1.763   -0.000)    2.493
   3.513   (   1.480   -1.480    0.000)    2.092
   3.631   (   1.506   -1.506    0.000)    2.129
   3.651   (   0.893   -0.893    0.000)    1.263
   3.699   (  -0.085    0.085   -0.000)    0.120
   3.810   (   4.370   -4.370    0.000)    6.180
   4.034   (  -0.878    0.878   -0.000)    1.241
   4.152   (  -2.020    2.020   -0.000)    2.857
   4.407   (   7.228   -7.228    0.000)   10.222
   5.014   (  -1.221    1.221   -0.000)    1.726
   5.059   (  -8.850    8.850   -0.000)   12.516
======================= Grid point 99 (31/32) =======================
q-point: ( 0.42  0.58 -0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 42
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.596   (   9.289    0.000    0.000)    9.289
   0.738   (  12.983    0.000    0.000)   12.983
   0.800   (  -7.730   -0.000   -0.000)    7.730
   1.023   (  10.547    0.000    0.000)   10.547
   1.038   (  -1.563   -0.000   -0.000)    1.563
   1.092   (  -4.073   -0.000   -0.000)    4.073
   1.322   (   0.370    0.000    0.000)    0.370
   1.341   (  15.492    0.000    0.000)   15.492
   1.346   (   2.339    0.000    0.000)    2.339
   1.421   (  -5.079   -0.000   -0.000)    5.079
   1.522   (  -3.881   -0.000   -0.000)    3.881
   1.529   (   2.210    0.000    0.000)    2.210
   1.568   (   5.884    0.000    0.000)    5.884
   1.639   (  13.086    0.000    0.000)   13.086
   1.705   (  -0.967   -0.000   -0.000)    0.967
   1.716   (   0.246    0.000    0.000)    0.246
   1.798   (   4.074    0.000    0.000)    4.074
   1.825   (  -2.546   -0.000   -0.000)    2.546
   1.831   (  -4.421   -0.000   -0.000)    4.421
   1.845   (  -6.012   -0.000   -0.000)    6.012
   1.945   (   3.118    0.000    0.000)    3.118
   1.951   (   7.176    0.000    0.000)    7.176
   2.064   (  -5.505   -0.000   -0.000)    5.505
   2.102   (  -3.684   -0.000   -0.000)    3.684
   2.202   (   0.266    0.000    0.000)    0.266
   2.227   (   0.354    0.000    0.000)    0.354
   2.250   (  -6.432   -0.000   -0.000)    6.432
   2.334   (  -6.524   -0.000   -0.000)    6.524
   2.450   (  -0.251   -0.000   -0.000)    0.251
   2.461   (   2.071    0.000    0.000)    2.071
   2.526   (   1.205    0.000    0.000)    1.205
   2.529   (   5.493    0.000    0.000)    5.493
   2.718   (  15.391    0.000    0.000)   15.391
   2.811   (  -2.959   -0.000   -0.000)    2.959
   2.989   (  -2.308   -0.000   -0.000)    2.308
   3.023   (   3.536    0.000    0.000)    3.536
   3.047   (   6.306    0.000    0.000)    6.306
   3.259   (  -2.728   -0.000   -0.000)    2.728
   3.500   (   0.743    0.000    0.000)    0.743
   3.578   (  -0.949   -0.000   -0.000)    0.949
   3.663   (   5.832    0.000    0.000)    5.832
   3.684   (  -1.617   -0.000   -0.000)    1.617
   3.778   (   7.225    0.000    0.000)    7.225
   4.046   (  -1.735   -0.000   -0.000)    1.735
   4.192   (  -1.301   -0.000   -0.000)    1.301
   4.315   (   8.941    0.000    0.000)    8.941
   5.027   (  -1.614   -0.000   -0.000)    1.614
   5.154   (  -9.887   -0.000   -0.000)    9.887
======================= Grid point 121 (32/32) =======================
q-point: (-0.50 -0.50  0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 
Number of triplets: 20
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.489   (  -0.000    0.000    0.000)    0.000
   0.489   (  -0.000    0.000    0.000)    0.000
   0.898   (  -0.000    0.000    0.000)    0.000
   0.905   (  -0.000    0.000    0.000)    0.000
   0.905   (  -0.000    0.000    0.000)    0.000
   1.050   (  -0.000    0.000    0.000)    0.000
   1.156   (  -0.000    0.000    0.000)    0.000
   1.311   (  -0.000    0.000    0.000)    0.000
   1.311   (  -0.000    0.000    0.000)    0.000
   1.479   (  -0.000    0.000    0.000)    0.000
   1.496   (  -0.000    0.000    0.000)    0.000
   1.496   (  -0.000    0.000    0.000)    0.000
   1.562   (   0.000   -0.000   -0.000)    0.000
   1.562   (   0.000   -0.000   -0.000)    0.000
   1.695   (  -0.000    0.000    0.000)    0.000
   1.718   (  -0.000    0.000    0.000)    0.000
   1.803   (  -0.000    0.000    0.000)    0.000
   1.851   (  -0.000    0.000    0.000)    0.000
   1.851   (  -0.000    0.000    0.000)    0.000
   1.916   (  -0.000    0.000    0.000)    0.000
   1.918   (  -0.000    0.000    0.000)    0.000
   1.918   (  -0.000    0.000    0.000)    0.000
   2.120   (   0.000   -0.000   -0.000)    0.000
   2.134   (  -0.000    0.000    0.000)    0.000
   2.204   (  -0.000    0.000    0.000)    0.000
   2.229   (  -0.000    0.000    0.000)    0.000
   2.339   (  -0.000    0.000    0.000)    0.000
   2.423   (  -0.000    0.000    0.000)    0.000
   2.423   (  -0.000    0.000    0.000)    0.000
   2.450   (  -0.000    0.000    0.000)    0.000
   2.450   (  -0.000    0.000    0.000)    0.000
   2.499   (  -0.000    0.000    0.000)    0.000
   2.545   (  -0.000    0.000    0.000)    0.000
   2.845   (  -0.000    0.000    0.000)    0.000
   2.982   (  -0.000    0.000    0.000)    0.000
   2.982   (  -0.000    0.000    0.000)    0.000
   3.012   (  -0.000    0.000    0.000)    0.000
   3.309   (  -0.000    0.000    0.000)    0.000
   3.497   (  -0.000    0.000    0.000)    0.000
   3.583   (  -0.000    0.000    0.000)    0.000
   3.583   (  -0.000    0.000    0.000)    0.000
   3.700   (  -0.000    0.000    0.000)    0.000
   3.700   (  -0.000    0.000    0.000)    0.000
   4.083   (  -0.000    0.000    0.000)    0.000
   4.206   (  -0.000    0.000    0.000)    0.000
   4.206   (  -0.000    0.000    0.000)    0.000
   5.044   (  -0.000    0.000    0.000)    0.000
   5.261   (  -0.000    0.000    0.000)    0.000
=================== 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/10368
   10.0     12.201     12.201      5.650     -0.000     -0.000      0.000 3/10368
   20.0      3.200      3.200      1.349     -0.000     -0.000      0.000 3/10368
   30.0      1.803      1.803      0.747     -0.000     -0.000      0.000 3/10368
   40.0      1.290      1.290      0.533     -0.000     -0.000      0.000 3/10368
   50.0      1.022      1.022      0.420     -0.000     -0.000      0.000 3/10368
   60.0      0.855      0.855      0.350     -0.000     -0.000      0.000 3/10368
   70.0      0.738      0.738      0.301     -0.000     -0.000      0.000 3/10368
   80.0      0.650      0.650      0.264     -0.000     -0.000      0.000 3/10368
   90.0      0.582      0.582      0.236     -0.000     -0.000      0.000 3/10368
  100.0      0.527      0.527      0.213     -0.000     -0.000      0.000 3/10368
  110.0      0.481      0.481      0.194     -0.000     -0.000      0.000 3/10368
  120.0      0.443      0.443      0.178     -0.000     -0.000      0.000 3/10368
  130.0      0.410      0.410      0.165     -0.000     -0.000      0.000 3/10368
  140.0      0.382      0.382      0.153     -0.000     -0.000      0.000 3/10368
  150.0      0.357      0.357      0.143     -0.000     -0.000      0.000 3/10368
  160.0      0.336      0.336      0.135     -0.000     -0.000      0.000 3/10368
  170.0      0.316      0.316      0.127     -0.000     -0.000      0.000 3/10368
  180.0      0.299      0.299      0.120     -0.000     -0.000      0.000 3/10368
  190.0      0.284      0.284      0.114     -0.000     -0.000      0.000 3/10368
  200.0      0.270      0.270      0.108     -0.000     -0.000      0.000 3/10368
  210.0      0.258      0.258      0.103     -0.000     -0.000      0.000 3/10368
  220.0      0.246      0.246      0.098     -0.000     -0.000      0.000 3/10368
  230.0      0.236      0.236      0.094     -0.000     -0.000      0.000 3/10368
  240.0      0.226      0.226      0.090     -0.000     -0.000      0.000 3/10368
  250.0      0.217      0.217      0.087     -0.000     -0.000      0.000 3/10368
  260.0      0.209      0.209      0.083     -0.000     -0.000      0.000 3/10368
  270.0      0.201      0.201      0.080     -0.000     -0.000      0.000 3/10368
  280.0      0.194      0.194      0.077     -0.000     -0.000      0.000 3/10368
  290.0      0.187      0.187      0.075     -0.000     -0.000      0.000 3/10368
  300.0      0.181      0.181      0.072     -0.000     -0.000      0.000 3/10368
  310.0      0.175      0.175      0.070     -0.000     -0.000      0.000 3/10368
  320.0      0.170      0.170      0.068     -0.000     -0.000      0.000 3/10368
  330.0      0.165      0.165      0.066     -0.000     -0.000      0.000 3/10368
  340.0      0.160      0.160      0.064     -0.000     -0.000      0.000 3/10368
  350.0      0.156      0.156      0.062     -0.000     -0.000      0.000 3/10368
  360.0      0.151      0.151      0.060     -0.000     -0.000      0.000 3/10368
  370.0      0.147      0.147      0.059     -0.000     -0.000      0.000 3/10368
  380.0      0.143      0.143      0.057     -0.000     -0.000      0.000 3/10368
  390.0      0.140      0.140      0.056     -0.000     -0.000      0.000 3/10368
  400.0      0.136      0.136      0.054     -0.000     -0.000      0.000 3/10368
  410.0      0.133      0.133      0.053     -0.000     -0.000      0.000 3/10368
  420.0      0.130      0.130      0.052     -0.000     -0.000      0.000 3/10368
  430.0      0.127      0.127      0.051     -0.000     -0.000      0.000 3/10368
  440.0      0.124      0.124      0.049     -0.000     -0.000      0.000 3/10368
  450.0      0.121      0.121      0.048     -0.000     -0.000      0.000 3/10368
  460.0      0.119      0.119      0.047     -0.000     -0.000      0.000 3/10368
  470.0      0.116      0.116      0.046     -0.000     -0.000      0.000 3/10368
  480.0      0.114      0.114      0.045     -0.000     -0.000      0.000 3/10368
  490.0      0.111      0.111      0.044     -0.000     -0.000      0.000 3/10368
  500.0      0.109      0.109      0.044     -0.000     -0.000      0.000 3/10368
  510.0      0.107      0.107      0.043     -0.000     -0.000      0.000 3/10368
  520.0      0.105      0.105      0.042     -0.000     -0.000      0.000 3/10368
  530.0      0.103      0.103      0.041     -0.000     -0.000      0.000 3/10368
  540.0      0.101      0.101      0.040     -0.000     -0.000      0.000 3/10368
  550.0      0.099      0.099      0.040     -0.000     -0.000      0.000 3/10368
  560.0      0.098      0.098      0.039     -0.000     -0.000      0.000 3/10368
  570.0      0.096      0.096      0.038     -0.000     -0.000      0.000 3/10368
  580.0      0.094      0.094      0.038     -0.000     -0.000      0.000 3/10368
  590.0      0.093      0.093      0.037     -0.000     -0.000      0.000 3/10368
  600.0      0.091      0.091      0.036     -0.000     -0.000      0.000 3/10368
  610.0      0.090      0.090      0.036     -0.000     -0.000      0.000 3/10368
  620.0      0.088      0.088      0.035     -0.000     -0.000      0.000 3/10368
  630.0      0.087      0.087      0.035     -0.000     -0.000      0.000 3/10368
  640.0      0.085      0.085      0.034     -0.000     -0.000      0.000 3/10368
  650.0      0.084      0.084      0.033     -0.000     -0.000      0.000 3/10368
  660.0      0.083      0.083      0.033     -0.000     -0.000      0.000 3/10368
  670.0      0.082      0.082      0.032     -0.000     -0.000      0.000 3/10368
  680.0      0.080      0.080      0.032     -0.000     -0.000      0.000 3/10368
  690.0      0.079      0.079      0.032     -0.000     -0.000      0.000 3/10368
  700.0      0.078      0.078      0.031     -0.000     -0.000      0.000 3/10368
  710.0      0.077      0.077      0.031     -0.000     -0.000      0.000 3/10368
  720.0      0.076      0.076      0.030     -0.000     -0.000      0.000 3/10368
  730.0      0.075      0.075      0.030     -0.000     -0.000      0.000 3/10368
  740.0      0.074      0.074      0.029     -0.000     -0.000      0.000 3/10368
  750.0      0.073      0.073      0.029     -0.000     -0.000      0.000 3/10368
  760.0      0.072      0.072      0.029     -0.000     -0.000      0.000 3/10368
  770.0      0.071      0.071      0.028     -0.000     -0.000      0.000 3/10368
  780.0      0.070      0.070      0.028     -0.000     -0.000      0.000 3/10368
  790.0      0.069      0.069      0.028     -0.000     -0.000      0.000 3/10368
  800.0      0.068      0.068      0.027     -0.000     -0.000      0.000 3/10368
  810.0      0.067      0.067      0.027     -0.000     -0.000      0.000 3/10368
  820.0      0.067      0.067      0.027     -0.000     -0.000      0.000 3/10368
  830.0      0.066      0.066      0.026     -0.000     -0.000      0.000 3/10368
  840.0      0.065      0.065      0.026     -0.000     -0.000      0.000 3/10368
  850.0      0.064      0.064      0.026     -0.000     -0.000      0.000 3/10368
  860.0      0.064      0.064      0.025     -0.000     -0.000      0.000 3/10368
  870.0      0.063      0.063      0.025     -0.000     -0.000      0.000 3/10368
  880.0      0.062      0.062      0.025     -0.000     -0.000      0.000 3/10368
  890.0      0.061      0.061      0.024     -0.000     -0.000      0.000 3/10368
  900.0      0.061      0.061      0.024     -0.000     -0.000      0.000 3/10368
  910.0      0.060      0.060      0.024     -0.000     -0.000      0.000 3/10368
  920.0      0.059      0.059      0.024     -0.000     -0.000      0.000 3/10368
  930.0      0.059      0.059      0.023     -0.000     -0.000      0.000 3/10368
  940.0      0.058      0.058      0.023     -0.000     -0.000      0.000 3/10368
  950.0      0.058      0.058      0.023     -0.000     -0.000      0.000 3/10368
  960.0      0.057      0.057      0.023     -0.000     -0.000      0.000 3/10368
  970.0      0.056      0.056      0.022     -0.000     -0.000      0.000 3/10368
  980.0      0.056      0.056      0.022     -0.000     -0.000      0.000 3/10368
  990.0      0.055      0.055      0.022     -0.000     -0.000      0.000 3/10368
 1000.0      0.055      0.055      0.022     -0.000     -0.000      0.000 3/10368

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

