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

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

-------------------------[time 2026-01-08 04:05: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: [1 1 1]
  Primitive matrix:
    [-0.5  0.5  0.5]
    [ 0.5 -0.5  0.5]
    [ 0.5  0.5 -0.5]
Spacegroup: I-43m (217)
Number of symmetry operations in supercell: 48
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   -5.542292829999998    5.542292829999998    5.542292829999998
  b    5.542292829999998   -5.542292829999998    5.542292829999998
  c    5.542292829999998    5.542292829999998   -5.542292829999998
Atomic positions (fractional):
   *1 Na  0.50000000000000  0.25000000000000  0.75000000000000  22.990
    2 Na  0.75000000000000  0.25000000000000  0.50000000000000  22.990
    3 Na  0.25000000000000  0.75000000000000  0.50000000000000  22.990
    4 Na  0.50000000000000  0.75000000000000  0.25000000000000  22.990
   *5 Na  0.61379061177228  0.00000000000000  0.00000000000000  22.990
    6 Na  0.38620938822772  0.38620938822772  0.38620938822772  22.990
    7 Na  0.00000000000000  0.61379061177228  0.00000000000000  22.990
    8 Na  0.00000000000000  0.00000000000000  0.61379061177228  22.990
    9 Na  0.75000000000000  0.50000000000000  0.25000000000000  22.990
   10 Na  0.25000000000000  0.50000000000000  0.75000000000000  22.990
  *11 Sn  0.00000000000000  0.18827691529716  0.37044940180797 118.710
   12 Sn  0.81172308470284  0.18217248651082  0.81172308470284 118.710
   13 Sn  0.62955059819203  0.62955059819203  0.81782751348918 118.710
   14 Sn  0.37044940180797  0.18827691529716  0.00000000000000 118.710
   15 Sn  0.18827691529716  0.37044940180797  0.00000000000000 118.710
   16 Sn  0.37044940180797  0.00000000000000  0.18827691529716 118.710
   17 Sn  0.62955059819203  0.81782751348918  0.62955059819203 118.710
   18 Sn  0.81782751348918  0.62955059819203  0.62955059819203 118.710
   19 Sn  0.00000000000000  0.37044940180797  0.18827691529716 118.710
   20 Sn  0.81172308470284  0.81172308470284  0.18217248651082 118.710
   21 Sn  0.18217248651082  0.81172308470284  0.81172308470284 118.710
   22 Sn  0.18827691529716  0.00000000000000  0.37044940180797 118.710
  *23 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a   11.084585659999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.084585659999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.084585659999997
Atomic positions (fractional):
   *1 Na  0.25000000000000  0.50000000000000  0.00000000000000  22.990 > 1
    2 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 2
    3 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 3
    4 Na  0.25000000000000  0.00000000000000  0.50000000000000  22.990 > 4
   *5 Na  0.19310469411386  0.80689530588614  0.80689530588614  22.990 > 5
    6 Na  0.19310469411386  0.19310469411386  0.19310469411386  22.990 > 6
    7 Na  0.80689530588614  0.19310469411386  0.80689530588614  22.990 > 7
    8 Na  0.30689530588614  0.30689530588614  0.69310469411386  22.990 > 8
    9 Na  0.00000000000000  0.25000000000000  0.50000000000000  22.990 > 9
   10 Na  0.50000000000000  0.25000000000000  0.00000000000000  22.990 > 10
   11 Na  0.75000000000000  0.00000000000000  0.50000000000000  22.990 > 1
   12 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 2
   13 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 3
   14 Na  0.75000000000000  0.50000000000000  0.00000000000000  22.990 > 4
   15 Na  0.69310469411386  0.30689530588614  0.30689530588614  22.990 > 5
   16 Na  0.69310469411386  0.69310469411386  0.69310469411386  22.990 > 6
   17 Na  0.30689530588614  0.69310469411386  0.30689530588614  22.990 > 7
   18 Na  0.80689530588614  0.80689530588614  0.19310469411386  22.990 > 8
   19 Na  0.50000000000000  0.75000000000000  0.00000000000000  22.990 > 9
   20 Na  0.00000000000000  0.75000000000000  0.50000000000000  22.990 > 10
  *21 Sn  0.27936315855256  0.09108624325541  0.90891375674459 118.710 > 11
   22 Sn  0.09108624325541  0.72063684144744  0.09108624325541 118.710 > 12
   23 Sn  0.40891375674459  0.40891375674459  0.22063684144744 118.710 > 13
   24 Sn  0.90891375674459  0.09108624325541  0.27936315855256 118.710 > 14
   25 Sn  0.09108624325541  0.90891375674459  0.27936315855256 118.710 > 15
   26 Sn  0.40891375674459  0.77936315855256  0.59108624325541 118.710 > 16
   27 Sn  0.90891375674459  0.72063684144744  0.90891375674459 118.710 > 17
   28 Sn  0.72063684144744  0.90891375674459  0.90891375674459 118.710 > 18
   29 Sn  0.77936315855256  0.40891375674459  0.59108624325541 118.710 > 19
   30 Sn  0.09108624325541  0.09108624325541  0.72063684144744 118.710 > 20
   31 Sn  0.72063684144744  0.09108624325541  0.09108624325541 118.710 > 21
   32 Sn  0.09108624325541  0.27936315855256  0.90891375674459 118.710 > 22
   33 Sn  0.77936315855256  0.59108624325541  0.40891375674459 118.710 > 11
   34 Sn  0.59108624325541  0.22063684144744  0.59108624325541 118.710 > 12
   35 Sn  0.90891375674459  0.90891375674459  0.72063684144744 118.710 > 13
   36 Sn  0.40891375674459  0.59108624325541  0.77936315855256 118.710 > 14
   37 Sn  0.59108624325541  0.40891375674459  0.77936315855256 118.710 > 15
   38 Sn  0.90891375674459  0.27936315855256  0.09108624325541 118.710 > 16
   39 Sn  0.40891375674459  0.22063684144744  0.40891375674459 118.710 > 17
   40 Sn  0.22063684144744  0.40891375674459  0.40891375674459 118.710 > 18
   41 Sn  0.27936315855256  0.90891375674459  0.09108624325541 118.710 > 19
   42 Sn  0.59108624325541  0.59108624325541  0.22063684144744 118.710 > 20
   43 Sn  0.22063684144744  0.59108624325541  0.59108624325541 118.710 > 21
   44 Sn  0.59108624325541  0.77936315855256  0.40891375674459 118.710 > 22
  *45 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054 > 23
   46 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 23
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   11.084585659999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.084585659999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.084585659999997
Atomic positions (fractional):
   *1 Na  0.25000000000000  0.50000000000000  0.00000000000000  22.990 > 1
    2 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 2
    3 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 3
    4 Na  0.25000000000000  0.00000000000000  0.50000000000000  22.990 > 4
   *5 Na  0.19310469411386  0.80689530588614  0.80689530588614  22.990 > 5
    6 Na  0.19310469411386  0.19310469411386  0.19310469411386  22.990 > 6
    7 Na  0.80689530588614  0.19310469411386  0.80689530588614  22.990 > 7
    8 Na  0.30689530588614  0.30689530588614  0.69310469411386  22.990 > 8
    9 Na  0.00000000000000  0.25000000000000  0.50000000000000  22.990 > 9
   10 Na  0.50000000000000  0.25000000000000  0.00000000000000  22.990 > 10
   11 Na  0.75000000000000  0.00000000000000  0.50000000000000  22.990 > 1
   12 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 2
   13 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 3
   14 Na  0.75000000000000  0.50000000000000  0.00000000000000  22.990 > 4
   15 Na  0.69310469411386  0.30689530588614  0.30689530588614  22.990 > 5
   16 Na  0.69310469411386  0.69310469411386  0.69310469411386  22.990 > 6
   17 Na  0.30689530588614  0.69310469411386  0.30689530588614  22.990 > 7
   18 Na  0.80689530588614  0.80689530588614  0.19310469411386  22.990 > 8
   19 Na  0.50000000000000  0.75000000000000  0.00000000000000  22.990 > 9
   20 Na  0.00000000000000  0.75000000000000  0.50000000000000  22.990 > 10
  *21 Sn  0.27936315855256  0.09108624325541  0.90891375674459 118.710 > 11
   22 Sn  0.09108624325541  0.72063684144744  0.09108624325541 118.710 > 12
   23 Sn  0.40891375674459  0.40891375674459  0.22063684144744 118.710 > 13
   24 Sn  0.90891375674459  0.09108624325541  0.27936315855256 118.710 > 14
   25 Sn  0.09108624325541  0.90891375674459  0.27936315855256 118.710 > 15
   26 Sn  0.40891375674459  0.77936315855256  0.59108624325541 118.710 > 16
   27 Sn  0.90891375674459  0.72063684144744  0.90891375674459 118.710 > 17
   28 Sn  0.72063684144744  0.90891375674459  0.90891375674459 118.710 > 18
   29 Sn  0.77936315855256  0.40891375674459  0.59108624325541 118.710 > 19
   30 Sn  0.09108624325541  0.09108624325541  0.72063684144744 118.710 > 20
   31 Sn  0.72063684144744  0.09108624325541  0.09108624325541 118.710 > 21
   32 Sn  0.09108624325541  0.27936315855256  0.90891375674459 118.710 > 22
   33 Sn  0.77936315855256  0.59108624325541  0.40891375674459 118.710 > 11
   34 Sn  0.59108624325541  0.22063684144744  0.59108624325541 118.710 > 12
   35 Sn  0.90891375674459  0.90891375674459  0.72063684144744 118.710 > 13
   36 Sn  0.40891375674459  0.59108624325541  0.77936315855256 118.710 > 14
   37 Sn  0.59108624325541  0.40891375674459  0.77936315855256 118.710 > 15
   38 Sn  0.90891375674459  0.27936315855256  0.09108624325541 118.710 > 16
   39 Sn  0.40891375674459  0.22063684144744  0.40891375674459 118.710 > 17
   40 Sn  0.22063684144744  0.40891375674459  0.40891375674459 118.710 > 18
   41 Sn  0.27936315855256  0.90891375674459  0.09108624325541 118.710 > 19
   42 Sn  0.59108624325541  0.59108624325541  0.22063684144744 118.710 > 20
   43 Sn  0.22063684144744  0.59108624325541  0.59108624325541 118.710 > 21
   44 Sn  0.59108624325541  0.77936315855256  0.40891375674459 118.710 > 22
  *45 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054 > 23
   46 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 23
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           14.6569546    0.0000000    0.0000000
            0.0000000   14.6569546    0.0000000
            0.0000000    0.0000000   14.6569546
-------------------------- Born effective charges --------------------------
    1 Na    1.2139859    0.0000000    0.0000000
           -0.0000000    1.7042895    0.3452216
            0.0000000   -0.3452216    1.7042895
    2 Na    1.7042895   -0.3452216    0.0000000
            0.3452216    1.7042895   -0.0000000
            0.0000000   -0.0000000    1.2139859
    3 Na    1.7042895    0.3452216    0.0000000
           -0.3452216    1.7042895   -0.0000000
            0.0000000   -0.0000000    1.2139859
    4 Na    1.2139859    0.0000000    0.0000000
           -0.0000000    1.7042895   -0.3452216
            0.0000000    0.3452216    1.7042895
    5 Na    1.0793004   -0.6671811   -0.6671811
           -0.6671811    1.0793004    0.6671811
           -0.6671811    0.6671811    1.0793004
    6 Na    1.0793004    0.6671811    0.6671811
            0.6671811    1.0793004    0.6671811
            0.6671811    0.6671811    1.0793004
    7 Na    1.0793004   -0.6671811    0.6671811
           -0.6671811    1.0793004   -0.6671811
            0.6671811   -0.6671811    1.0793004
    8 Na    1.0793004    0.6671811   -0.6671811
            0.6671811    1.0793004   -0.6671811
           -0.6671811   -0.6671811    1.0793004
    9 Na    1.7042895    0.0000000   -0.3452216
           -0.0000000    1.2139859   -0.0000000
            0.3452216   -0.0000000    1.7042895
   10 Na    1.7042895    0.0000000    0.3452216
           -0.0000000    1.2139859   -0.0000000
           -0.3452216   -0.0000000    1.7042895
   11 Sn   -1.8122793    0.1231054   -0.1231054
           -0.3039164   -0.7036430    0.0687038
            0.3039164    0.0687038   -0.7036430
   12 Sn   -0.7036430    0.3039164   -0.0687038
           -0.1231054   -1.8122793   -0.1231054
           -0.0687038    0.3039164   -0.7036430
   13 Sn   -0.7036430   -0.0687038   -0.3039164
           -0.0687038   -0.7036430   -0.3039164
            0.1231054    0.1231054   -1.8122793
   14 Sn   -0.7036430    0.0687038    0.3039164
            0.0687038   -0.7036430   -0.3039164
           -0.1231054    0.1231054   -1.8122793
   15 Sn   -0.7036430    0.0687038   -0.3039164
            0.0687038   -0.7036430    0.3039164
            0.1231054   -0.1231054   -1.8122793
   16 Sn   -0.7036430    0.3039164    0.0687038
           -0.1231054   -1.8122793    0.1231054
            0.0687038   -0.3039164   -0.7036430
   17 Sn   -0.7036430   -0.3039164   -0.0687038
            0.1231054   -1.8122793    0.1231054
           -0.0687038   -0.3039164   -0.7036430
   18 Sn   -1.8122793    0.1231054    0.1231054
           -0.3039164   -0.7036430   -0.0687038
           -0.3039164   -0.0687038   -0.7036430
   19 Sn   -1.8122793   -0.1231054    0.1231054
            0.3039164   -0.7036430    0.0687038
           -0.3039164    0.0687038   -0.7036430
   20 Sn   -0.7036430   -0.0687038    0.3039164
           -0.0687038   -0.7036430    0.3039164
           -0.1231054   -0.1231054   -1.8122793
   21 Sn   -1.8122793   -0.1231054   -0.1231054
            0.3039164   -0.7036430   -0.0687038
            0.3039164   -0.0687038   -0.7036430
   22 Sn   -0.7036430   -0.3039164    0.0687038
            0.1231054   -1.8122793   -0.1231054
            0.0687038    0.3039164   -0.7036430
   23 Yb   -0.6840696    0.0000000    0.0000000
           -0.0000000   -0.6840696   -0.0000000
            0.0000000   -0.0000000   -0.6840696
----------------------------------------------------------------------------
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: 138/138
Permutation basis: 5829/5829
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 235
Number of blocks in projector: 235
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 140
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 95
Use standard eigh solver.
Tree of FC basis block matrices:
- (235, 224), data: False
|-- (95, 89), data: True
|-- (140, 135), data: True
-----
Solver_atoms: 1 -- 46 / 46
Time (Solver_compr_matrix_reshape): 0.001
Solver_block: 100 / 160
 - Time: 0.097
Solver_block: 160 / 160
 - Time: 0.044
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.143
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 138/138
Permutation basis: 5829/5829
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 235
Number of blocks in projector: 235
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 140
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 95
Use standard eigh solver.
Tree of FC basis block matrices:
- (235, 224), data: False
|-- (95, 89), data: True
|-- (140, 135), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) 
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 04:05:52]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 04:05:52]-------------------------
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): [1 1 1]
Primitive matrix:
  [-0.5  0.5  0.5]
  [ 0.5 -0.5  0.5]
  [ 0.5  0.5 -0.5]
Spacegroup: I-43m (217)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   -5.542292829999998    5.542292829999998    5.542292829999998
  b    5.542292829999998   -5.542292829999998    5.542292829999998
  c    5.542292829999998    5.542292829999998   -5.542292829999998
Atomic positions (fractional):
    1 Na  0.50000000000000  0.25000000000000  0.75000000000000  22.990
    2 Na  0.75000000000000  0.25000000000000  0.50000000000000  22.990
    3 Na  0.25000000000000  0.75000000000000  0.50000000000000  22.990
    4 Na  0.50000000000000  0.75000000000000  0.25000000000000  22.990
    5 Na  0.61379061177228  0.00000000000000  0.00000000000000  22.990
    6 Na  0.38620938822772  0.38620938822772  0.38620938822772  22.990
    7 Na  0.00000000000000  0.61379061177228  0.00000000000000  22.990
    8 Na  0.00000000000000  0.00000000000000  0.61379061177228  22.990
    9 Na  0.75000000000000  0.50000000000000  0.25000000000000  22.990
   10 Na  0.25000000000000  0.50000000000000  0.75000000000000  22.990
   11 Sn  0.00000000000000  0.18827691529716  0.37044940180797 118.710
   12 Sn  0.81172308470284  0.18217248651082  0.81172308470284 118.710
   13 Sn  0.62955059819203  0.62955059819203  0.81782751348918 118.710
   14 Sn  0.37044940180797  0.18827691529716  0.00000000000000 118.710
   15 Sn  0.18827691529716  0.37044940180797  0.00000000000000 118.710
   16 Sn  0.37044940180797  0.00000000000000  0.18827691529716 118.710
   17 Sn  0.62955059819203  0.81782751348918  0.62955059819203 118.710
   18 Sn  0.81782751348918  0.62955059819203  0.62955059819203 118.710
   19 Sn  0.00000000000000  0.37044940180797  0.18827691529716 118.710
   20 Sn  0.81172308470284  0.81172308470284  0.18217248651082 118.710
   21 Sn  0.18217248651082  0.81172308470284  0.81172308470284 118.710
   22 Sn  0.18827691529716  0.00000000000000  0.37044940180797 118.710
   23 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.084585659999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.084585659999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.084585659999997
Atomic positions (fractional):
    1 Na  0.25000000000000  0.50000000000000  0.00000000000000  22.990 > 1
    2 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 2
    3 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 3
    4 Na  0.25000000000000  0.00000000000000  0.50000000000000  22.990 > 4
    5 Na  0.19310469411386  0.80689530588614  0.80689530588614  22.990 > 5
    6 Na  0.19310469411386  0.19310469411386  0.19310469411386  22.990 > 6
    7 Na  0.80689530588614  0.19310469411386  0.80689530588614  22.990 > 7
    8 Na  0.30689530588614  0.30689530588614  0.69310469411386  22.990 > 8
    9 Na  0.00000000000000  0.25000000000000  0.50000000000000  22.990 > 9
   10 Na  0.50000000000000  0.25000000000000  0.00000000000000  22.990 > 10
   11 Na  0.75000000000000  0.00000000000000  0.50000000000000  22.990 > 1
   12 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 2
   13 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 3
   14 Na  0.75000000000000  0.50000000000000  0.00000000000000  22.990 > 4
   15 Na  0.69310469411386  0.30689530588614  0.30689530588614  22.990 > 5
   16 Na  0.69310469411386  0.69310469411386  0.69310469411386  22.990 > 6
   17 Na  0.30689530588614  0.69310469411386  0.30689530588614  22.990 > 7
   18 Na  0.80689530588614  0.80689530588614  0.19310469411386  22.990 > 8
   19 Na  0.50000000000000  0.75000000000000  0.00000000000000  22.990 > 9
   20 Na  0.00000000000000  0.75000000000000  0.50000000000000  22.990 > 10
   21 Sn  0.27936315855256  0.09108624325541  0.90891375674459 118.710 > 21
   22 Sn  0.09108624325541  0.72063684144744  0.09108624325541 118.710 > 22
   23 Sn  0.40891375674459  0.40891375674459  0.22063684144744 118.710 > 23
   24 Sn  0.90891375674459  0.09108624325541  0.27936315855256 118.710 > 24
   25 Sn  0.09108624325541  0.90891375674459  0.27936315855256 118.710 > 25
   26 Sn  0.40891375674459  0.77936315855256  0.59108624325541 118.710 > 26
   27 Sn  0.90891375674459  0.72063684144744  0.90891375674459 118.710 > 27
   28 Sn  0.72063684144744  0.90891375674459  0.90891375674459 118.710 > 28
   29 Sn  0.77936315855256  0.40891375674459  0.59108624325541 118.710 > 29
   30 Sn  0.09108624325541  0.09108624325541  0.72063684144744 118.710 > 30
   31 Sn  0.72063684144744  0.09108624325541  0.09108624325541 118.710 > 31
   32 Sn  0.09108624325541  0.27936315855256  0.90891375674459 118.710 > 32
   33 Sn  0.77936315855256  0.59108624325541  0.40891375674459 118.710 > 21
   34 Sn  0.59108624325541  0.22063684144744  0.59108624325541 118.710 > 22
   35 Sn  0.90891375674459  0.90891375674459  0.72063684144744 118.710 > 23
   36 Sn  0.40891375674459  0.59108624325541  0.77936315855256 118.710 > 24
   37 Sn  0.59108624325541  0.40891375674459  0.77936315855256 118.710 > 25
   38 Sn  0.90891375674459  0.27936315855256  0.09108624325541 118.710 > 26
   39 Sn  0.40891375674459  0.22063684144744  0.40891375674459 118.710 > 27
   40 Sn  0.22063684144744  0.40891375674459  0.40891375674459 118.710 > 28
   41 Sn  0.27936315855256  0.90891375674459  0.09108624325541 118.710 > 29
   42 Sn  0.59108624325541  0.59108624325541  0.22063684144744 118.710 > 30
   43 Sn  0.22063684144744  0.59108624325541  0.59108624325541 118.710 > 31
   44 Sn  0.59108624325541  0.77936315855256  0.40891375674459 118.710 > 32
   45 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054 > 45
   46 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 45
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           14.6569546    0.0000000    0.0000000
            0.0000000   14.6569546    0.0000000
            0.0000000    0.0000000   14.6569546
-------------------------- Born effective charges --------------------------
    1 Na    1.2139859    0.0000000    0.0000000
           -0.0000000    1.7042895    0.3452216
            0.0000000   -0.3452216    1.7042895
    2 Na    1.7042895   -0.3452216    0.0000000
            0.3452216    1.7042895   -0.0000000
            0.0000000   -0.0000000    1.2139859
    3 Na    1.7042895    0.3452216    0.0000000
           -0.3452216    1.7042895   -0.0000000
            0.0000000   -0.0000000    1.2139859
    4 Na    1.2139859    0.0000000    0.0000000
           -0.0000000    1.7042895   -0.3452216
            0.0000000    0.3452216    1.7042895
    5 Na    1.0793004   -0.6671811   -0.6671811
           -0.6671811    1.0793004    0.6671811
           -0.6671811    0.6671811    1.0793004
    6 Na    1.0793004    0.6671811    0.6671811
            0.6671811    1.0793004    0.6671811
            0.6671811    0.6671811    1.0793004
    7 Na    1.0793004   -0.6671811    0.6671811
           -0.6671811    1.0793004   -0.6671811
            0.6671811   -0.6671811    1.0793004
    8 Na    1.0793004    0.6671811   -0.6671811
            0.6671811    1.0793004   -0.6671811
           -0.6671811   -0.6671811    1.0793004
    9 Na    1.7042895    0.0000000   -0.3452216
           -0.0000000    1.2139859   -0.0000000
            0.3452216   -0.0000000    1.7042895
   10 Na    1.7042895    0.0000000    0.3452216
           -0.0000000    1.2139859   -0.0000000
           -0.3452216   -0.0000000    1.7042895
   11 Sn   -1.8122793    0.1231054   -0.1231054
           -0.3039164   -0.7036430    0.0687038
            0.3039164    0.0687038   -0.7036430
   12 Sn   -0.7036430    0.3039164   -0.0687038
           -0.1231054   -1.8122793   -0.1231054
           -0.0687038    0.3039164   -0.7036430
   13 Sn   -0.7036430   -0.0687038   -0.3039164
           -0.0687038   -0.7036430   -0.3039164
            0.1231054    0.1231054   -1.8122793
   14 Sn   -0.7036430    0.0687038    0.3039164
            0.0687038   -0.7036430   -0.3039164
           -0.1231054    0.1231054   -1.8122793
   15 Sn   -0.7036430    0.0687038   -0.3039164
            0.0687038   -0.7036430    0.3039164
            0.1231054   -0.1231054   -1.8122793
   16 Sn   -0.7036430    0.3039164    0.0687038
           -0.1231054   -1.8122793    0.1231054
            0.0687038   -0.3039164   -0.7036430
   17 Sn   -0.7036430   -0.3039164   -0.0687038
            0.1231054   -1.8122793    0.1231054
           -0.0687038   -0.3039164   -0.7036430
   18 Sn   -1.8122793    0.1231054    0.1231054
           -0.3039164   -0.7036430   -0.0687038
           -0.3039164   -0.0687038   -0.7036430
   19 Sn   -1.8122793   -0.1231054    0.1231054
            0.3039164   -0.7036430    0.0687038
           -0.3039164    0.0687038   -0.7036430
   20 Sn   -0.7036430   -0.0687038    0.3039164
           -0.0687038   -0.7036430    0.3039164
           -0.1231054   -0.1231054   -1.8122793
   21 Sn   -1.8122793   -0.1231054   -0.1231054
            0.3039164   -0.7036430   -0.0687038
            0.3039164   -0.0687038   -0.7036430
   22 Sn   -0.7036430   -0.3039164    0.0687038
            0.1231054   -1.8122793   -0.1231054
            0.0687038    0.3039164   -0.7036430
   23 Yb   -0.6840696    0.0000000    0.0000000
           -0.0000000   -0.6840696   -0.0000000
            0.0000000   -0.0000000   -0.6840696
----------------------------------------------------------------------------
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.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  0.0000]
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 5, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 21, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  0.0000]
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 45, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000001 (xxx) 0.00000001 (xxx) 0.00000001 (xxx)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (xz) -0.00000000 (xz) 
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 04:05:55]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 04:05:56]-------------------------
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): [1 1 1]
Primitive matrix:
  [-0.5  0.5  0.5]
  [ 0.5 -0.5  0.5]
  [ 0.5  0.5 -0.5]
Spacegroup: I-43m (217)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a   -5.542292829999998    5.542292829999998    5.542292829999998
  b    5.542292829999998   -5.542292829999998    5.542292829999998
  c    5.542292829999998    5.542292829999998   -5.542292829999998
Atomic positions (fractional):
    1 Na  0.50000000000000  0.25000000000000  0.75000000000000  22.990
    2 Na  0.75000000000000  0.25000000000000  0.50000000000000  22.990
    3 Na  0.25000000000000  0.75000000000000  0.50000000000000  22.990
    4 Na  0.50000000000000  0.75000000000000  0.25000000000000  22.990
    5 Na  0.61379061177228  0.00000000000000  0.00000000000000  22.990
    6 Na  0.38620938822772  0.38620938822772  0.38620938822772  22.990
    7 Na  0.00000000000000  0.61379061177228  0.00000000000000  22.990
    8 Na  0.00000000000000  0.00000000000000  0.61379061177228  22.990
    9 Na  0.75000000000000  0.50000000000000  0.25000000000000  22.990
   10 Na  0.25000000000000  0.50000000000000  0.75000000000000  22.990
   11 Sn  0.00000000000000  0.18827691529716  0.37044940180797 118.710
   12 Sn  0.81172308470284  0.18217248651082  0.81172308470284 118.710
   13 Sn  0.62955059819203  0.62955059819203  0.81782751348918 118.710
   14 Sn  0.37044940180797  0.18827691529716  0.00000000000000 118.710
   15 Sn  0.18827691529716  0.37044940180797  0.00000000000000 118.710
   16 Sn  0.37044940180797  0.00000000000000  0.18827691529716 118.710
   17 Sn  0.62955059819203  0.81782751348918  0.62955059819203 118.710
   18 Sn  0.81782751348918  0.62955059819203  0.62955059819203 118.710
   19 Sn  0.00000000000000  0.37044940180797  0.18827691529716 118.710
   20 Sn  0.81172308470284  0.81172308470284  0.18217248651082 118.710
   21 Sn  0.18217248651082  0.81172308470284  0.81172308470284 118.710
   22 Sn  0.18827691529716  0.00000000000000  0.37044940180797 118.710
   23 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.084585659999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.084585659999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.084585659999997
Atomic positions (fractional):
    1 Na  0.25000000000000  0.50000000000000  0.00000000000000  22.990 > 1
    2 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 2
    3 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 3
    4 Na  0.25000000000000  0.00000000000000  0.50000000000000  22.990 > 4
    5 Na  0.19310469411386  0.80689530588614  0.80689530588614  22.990 > 5
    6 Na  0.19310469411386  0.19310469411386  0.19310469411386  22.990 > 6
    7 Na  0.80689530588614  0.19310469411386  0.80689530588614  22.990 > 7
    8 Na  0.30689530588614  0.30689530588614  0.69310469411386  22.990 > 8
    9 Na  0.00000000000000  0.25000000000000  0.50000000000000  22.990 > 9
   10 Na  0.50000000000000  0.25000000000000  0.00000000000000  22.990 > 10
   11 Na  0.75000000000000  0.00000000000000  0.50000000000000  22.990 > 1
   12 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 2
   13 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 3
   14 Na  0.75000000000000  0.50000000000000  0.00000000000000  22.990 > 4
   15 Na  0.69310469411386  0.30689530588614  0.30689530588614  22.990 > 5
   16 Na  0.69310469411386  0.69310469411386  0.69310469411386  22.990 > 6
   17 Na  0.30689530588614  0.69310469411386  0.30689530588614  22.990 > 7
   18 Na  0.80689530588614  0.80689530588614  0.19310469411386  22.990 > 8
   19 Na  0.50000000000000  0.75000000000000  0.00000000000000  22.990 > 9
   20 Na  0.00000000000000  0.75000000000000  0.50000000000000  22.990 > 10
   21 Sn  0.27936315855256  0.09108624325541  0.90891375674459 118.710 > 21
   22 Sn  0.09108624325541  0.72063684144744  0.09108624325541 118.710 > 22
   23 Sn  0.40891375674459  0.40891375674459  0.22063684144744 118.710 > 23
   24 Sn  0.90891375674459  0.09108624325541  0.27936315855256 118.710 > 24
   25 Sn  0.09108624325541  0.90891375674459  0.27936315855256 118.710 > 25
   26 Sn  0.40891375674459  0.77936315855256  0.59108624325541 118.710 > 26
   27 Sn  0.90891375674459  0.72063684144744  0.90891375674459 118.710 > 27
   28 Sn  0.72063684144744  0.90891375674459  0.90891375674459 118.710 > 28
   29 Sn  0.77936315855256  0.40891375674459  0.59108624325541 118.710 > 29
   30 Sn  0.09108624325541  0.09108624325541  0.72063684144744 118.710 > 30
   31 Sn  0.72063684144744  0.09108624325541  0.09108624325541 118.710 > 31
   32 Sn  0.09108624325541  0.27936315855256  0.90891375674459 118.710 > 32
   33 Sn  0.77936315855256  0.59108624325541  0.40891375674459 118.710 > 21
   34 Sn  0.59108624325541  0.22063684144744  0.59108624325541 118.710 > 22
   35 Sn  0.90891375674459  0.90891375674459  0.72063684144744 118.710 > 23
   36 Sn  0.40891375674459  0.59108624325541  0.77936315855256 118.710 > 24
   37 Sn  0.59108624325541  0.40891375674459  0.77936315855256 118.710 > 25
   38 Sn  0.90891375674459  0.27936315855256  0.09108624325541 118.710 > 26
   39 Sn  0.40891375674459  0.22063684144744  0.40891375674459 118.710 > 27
   40 Sn  0.22063684144744  0.40891375674459  0.40891375674459 118.710 > 28
   41 Sn  0.27936315855256  0.90891375674459  0.09108624325541 118.710 > 29
   42 Sn  0.59108624325541  0.59108624325541  0.22063684144744 118.710 > 30
   43 Sn  0.22063684144744  0.59108624325541  0.59108624325541 118.710 > 31
   44 Sn  0.59108624325541  0.77936315855256  0.40891375674459 118.710 > 32
   45 Yb  0.00000000000000  0.00000000000000  0.00000000000000 173.054 > 45
   46 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 45
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           14.6569546    0.0000000    0.0000000
            0.0000000   14.6569546    0.0000000
            0.0000000    0.0000000   14.6569546
-------------------------- Born effective charges --------------------------
    1 Na    1.2139859    0.0000000    0.0000000
           -0.0000000    1.7042895    0.3452216
            0.0000000   -0.3452216    1.7042895
    2 Na    1.7042895   -0.3452216    0.0000000
            0.3452216    1.7042895   -0.0000000
            0.0000000   -0.0000000    1.2139859
    3 Na    1.7042895    0.3452216    0.0000000
           -0.3452216    1.7042895   -0.0000000
            0.0000000   -0.0000000    1.2139859
    4 Na    1.2139859    0.0000000    0.0000000
           -0.0000000    1.7042895   -0.3452216
            0.0000000    0.3452216    1.7042895
    5 Na    1.0793004   -0.6671811   -0.6671811
           -0.6671811    1.0793004    0.6671811
           -0.6671811    0.6671811    1.0793004
    6 Na    1.0793004    0.6671811    0.6671811
            0.6671811    1.0793004    0.6671811
            0.6671811    0.6671811    1.0793004
    7 Na    1.0793004   -0.6671811    0.6671811
           -0.6671811    1.0793004   -0.6671811
            0.6671811   -0.6671811    1.0793004
    8 Na    1.0793004    0.6671811   -0.6671811
            0.6671811    1.0793004   -0.6671811
           -0.6671811   -0.6671811    1.0793004
    9 Na    1.7042895    0.0000000   -0.3452216
           -0.0000000    1.2139859   -0.0000000
            0.3452216   -0.0000000    1.7042895
   10 Na    1.7042895    0.0000000    0.3452216
           -0.0000000    1.2139859   -0.0000000
           -0.3452216   -0.0000000    1.7042895
   11 Sn   -1.8122793    0.1231054   -0.1231054
           -0.3039164   -0.7036430    0.0687038
            0.3039164    0.0687038   -0.7036430
   12 Sn   -0.7036430    0.3039164   -0.0687038
           -0.1231054   -1.8122793   -0.1231054
           -0.0687038    0.3039164   -0.7036430
   13 Sn   -0.7036430   -0.0687038   -0.3039164
           -0.0687038   -0.7036430   -0.3039164
            0.1231054    0.1231054   -1.8122793
   14 Sn   -0.7036430    0.0687038    0.3039164
            0.0687038   -0.7036430   -0.3039164
           -0.1231054    0.1231054   -1.8122793
   15 Sn   -0.7036430    0.0687038   -0.3039164
            0.0687038   -0.7036430    0.3039164
            0.1231054   -0.1231054   -1.8122793
   16 Sn   -0.7036430    0.3039164    0.0687038
           -0.1231054   -1.8122793    0.1231054
            0.0687038   -0.3039164   -0.7036430
   17 Sn   -0.7036430   -0.3039164   -0.0687038
            0.1231054   -1.8122793    0.1231054
           -0.0687038   -0.3039164   -0.7036430
   18 Sn   -1.8122793    0.1231054    0.1231054
           -0.3039164   -0.7036430   -0.0687038
           -0.3039164   -0.0687038   -0.7036430
   19 Sn   -1.8122793   -0.1231054    0.1231054
            0.3039164   -0.7036430    0.0687038
           -0.3039164    0.0687038   -0.7036430
   20 Sn   -0.7036430   -0.0687038    0.3039164
           -0.0687038   -0.7036430    0.3039164
           -0.1231054   -0.1231054   -1.8122793
   21 Sn   -1.8122793   -0.1231054   -0.1231054
            0.3039164   -0.7036430   -0.0687038
            0.3039164   -0.0687038   -0.7036430
   22 Sn   -0.7036430   -0.3039164    0.0687038
            0.1231054   -1.8122793   -0.1231054
            0.0687038    0.3039164   -0.7036430
   23 Yb   -0.6840696    0.0000000    0.0000000
           -0.0000000   -0.6840696   -0.0000000
            0.0000000   -0.0000000   -0.6840696
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000001 (xxx) 0.00000001 (xxx) 0.00000001 (xxx)
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
--------------------------- 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 ... [ 6 6 6 ]
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.47, Number of G-points: 321, Lambda: 0.19
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/16) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 16
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.718   (   0.000   -0.000    0.000)    0.000
   0.718   (  -0.000    0.000    0.000)    0.000
   0.718   (   0.000    0.000    0.000)    0.000
   1.026   (  -0.000   -0.000    0.000)    0.000
   1.026   (  -0.000   -0.000    0.000)    0.000
   1.026   (   0.000   -0.000    0.000)    0.000
   1.370   (   0.000    0.000    0.000)    0.000
   1.433   (  -0.000    0.000    0.000)    0.000
   1.433   (   0.000    0.000    0.000)    0.000
   1.844   (   0.000    0.000    0.000)    0.000
   1.844   (   0.000    0.000    0.000)    0.000
   1.844   (   0.000   -0.000    0.000)    0.000
   2.263   (   0.000    0.000    0.000)    0.000
   2.263   (   0.000    0.000    0.000)    0.000
   2.263   (   0.000    0.000    0.000)    0.000
   2.472   (   0.000    0.000    0.000)    0.000
   2.472   (   0.000    0.000    0.000)    0.000
   2.472   (   0.000   -0.000    0.000)    0.000
   2.628   (   0.000   -0.000    0.000)    0.000
   2.628   (  -0.000    0.000    0.000)    0.000
   2.628   (  -0.000    0.000    0.000)    0.000
   2.660   (   0.000    0.000    0.000)    0.000
   2.660   (   0.000    0.000    0.000)    0.000
   2.755   (   0.000    0.000    0.000)    0.000
   2.755   (   0.000   -0.000    0.000)    0.000
   2.755   (   0.000    0.000    0.000)    0.000
   3.040   (  -0.000   -0.000    0.000)    0.000
   3.040   (  -0.000    0.000    0.000)    0.000
   3.040   (  -0.000    0.000    0.000)    0.000
   3.260   (   0.000    0.000    0.000)    0.000
   3.260   (   0.000   -0.000    0.000)    0.000
   3.260   (   0.000    0.000    0.000)    0.000
   3.374   (  -0.000    0.000    0.000)    0.000
   3.413   (  -0.000    0.000    0.000)    0.000
   3.413   (  -0.000    0.000    0.000)    0.000
   3.413   (  -0.000   -0.000    0.000)    0.000
   3.413   (   0.000   -0.000    0.000)    0.000
   3.475   (  -0.000   -0.000    0.000)    0.000
   3.475   (   0.000   -0.000    0.000)    0.000
   4.115   (   0.000    0.000    0.000)    0.000
   4.115   (   0.000    0.000    0.000)    0.000
   4.115   (  -0.000   -0.000    0.000)    0.000
   4.160   (   0.000    0.000    0.000)    0.000
   4.160   (   0.000    0.000    0.000)    0.000
   4.160   (   0.000   -0.000    0.000)    0.000
   4.176   (  -0.000   -0.000    0.000)    0.000
   4.176   (   0.000    0.000    0.000)    0.000
   4.209   (   0.000    0.000    0.000)    0.000
   4.209   (   0.000   -0.000    0.000)    0.000
   4.209   (  -0.000   -0.000    0.000)    0.000
   4.298   (   0.000   -0.000    0.000)    0.000
   4.298   (   0.000    0.000    0.000)    0.000
   4.298   (   0.000    0.000    0.000)    0.000
   4.341   (   0.000    0.000    0.000)    0.000
   4.479   (   0.000    0.000    0.000)    0.000
   4.479   (  -0.000   -0.000    0.000)    0.000
   4.479   (   0.000    0.000    0.000)    0.000
   4.620   (   0.000    0.000    0.000)    0.000
   4.620   (   0.000    0.000    0.000)    0.000
   4.620   (  -0.000   -0.000    0.000)    0.000
   4.700   (   0.000    0.000    0.000)    0.000
   4.700   (   0.000    0.000    0.000)    0.000
   5.194   (   0.000    0.000    0.000)    0.000
   5.194   (   0.000   -0.000    0.000)    0.000
   5.194   (   0.000    0.000    0.000)    0.000
   5.361   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/16) =======================
q-point: ( 0.17  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 42
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.198   (   0.000    6.677    6.677)    9.443
   0.362   (   0.000   11.026   11.026)   15.593
   0.632   (   0.000   19.570   19.570)   27.677
   0.704   (  -0.000   -0.997   -0.997)    1.410
   0.716   (   0.000    0.083    0.083)    0.117
   0.770   (   0.000    1.964    1.964)    2.778
   0.935   (  -0.000   -4.412   -4.412)    6.240
   0.958   (  -0.000   -3.825   -3.825)    5.410
   0.994   (  -0.000   -0.191   -0.191)    0.270
   1.408   (   0.000    0.539    0.539)    0.762
   1.455   (   0.000    1.008    1.008)    1.425
   1.463   (   0.000    2.437    2.437)    3.447
   1.796   (  -0.000   -3.079   -3.079)    4.354
   1.830   (  -0.000   -0.907   -0.907)    1.283
   1.875   (   0.000    1.685    1.685)    2.384
   2.258   (  -0.000   -0.280   -0.280)    0.396
   2.265   (   0.000    0.129    0.129)    0.183
   2.266   (   0.000    0.092    0.092)    0.130
   2.415   (  -0.000   -2.270   -2.270)    3.211
   2.453   (  -0.000   -2.373   -2.373)    3.356
   2.488   (   0.000    1.031    1.031)    1.458
   2.574   (  -0.000   -2.306   -2.306)    3.261
   2.578   (  -0.000   -1.498   -1.498)    2.118
   2.616   (  -0.000   -0.162   -0.162)    0.230
   2.648   (  -0.000   -0.030   -0.030)    0.043
   2.721   (   0.000    1.993    1.993)    2.818
   2.793   (   0.000    2.010    2.010)    2.843
   2.800   (   0.000    2.255    2.255)    3.188
   2.835   (   0.000    3.189    3.189)    4.510
   3.014   (  -0.000   -1.233   -1.233)    1.744
   3.041   (  -0.000   -0.061   -0.061)    0.087
   3.121   (  -0.000   -2.956   -2.956)    4.181
   3.251   (  -0.000   -0.681   -0.681)    0.963
   3.251   (  -0.000   -0.742   -0.742)    1.049
   3.285   (   0.000    1.667    1.667)    2.357
   3.329   (  -0.000   -2.046   -2.046)    2.893
   3.391   (  -0.000   -1.044   -1.044)    1.476
   3.433   (   0.000    1.309    1.309)    1.851
   3.448   (   0.000    4.225    4.225)    5.974
   3.478   (  -0.000   -0.067   -0.067)    0.094
   3.501   (   0.000    1.368    1.368)    1.935
   3.595   (   0.000    6.945    6.945)    9.822
   3.969   (  -0.000   -6.420   -6.420)    9.079
   4.020   (  -0.000   -4.222   -4.222)    5.971
   4.048   (  -0.000   -5.628   -5.628)    7.959
   4.075   (  -0.000   -3.396   -3.396)    4.802
   4.122   (  -0.000   -2.285   -2.285)    3.232
   4.134   (  -0.000   -1.241   -1.241)    1.754
   4.155   (  -0.000   -1.814   -1.814)    2.566
   4.179   (  -0.000   -0.759   -0.759)    1.074
   4.230   (   0.000    1.218    1.218)    1.722
   4.238   (   0.000    0.397    0.397)    0.562
   4.259   (  -0.000   -0.086   -0.086)    0.121
   4.339   (   0.000    1.990    1.990)    2.815
   4.363   (   0.000    2.959    2.959)    4.185
   4.374   (   0.000    1.147    1.147)    1.621
   4.416   (   0.000    3.342    3.342)    4.727
   4.483   (   0.000    0.018    0.018)    0.025
   4.502   (   0.000    1.105    1.105)    1.563
   4.513   (   0.000    0.402    0.402)    0.569
   4.590   (  -0.000   -1.175   -1.175)    1.662
   4.616   (  -0.000   -0.285   -0.285)    0.402
   4.695   (  -0.000   -0.443   -0.443)    0.626
   4.703   (   0.000    0.144    0.144)    0.204
   4.891   (  -0.000   -0.563   -0.563)    0.797
   5.190   (  -0.000   -0.194   -0.194)    0.275
   5.209   (   0.000    0.861    0.861)    1.218
   5.219   (   0.000    0.463    0.463)    0.654
   5.347   (  -0.000   -0.886   -0.886)    1.253
======================= Grid point 2 (3/16) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 46
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.390   (   0.000    5.740    5.740)    8.118
   0.584   (   0.000    1.777    1.777)    2.513
   0.650   (  -0.000   -2.849   -2.849)    4.029
   0.723   (   0.000    0.257    0.257)    0.363
   0.808   (   0.000    3.471    3.471)    4.909
   0.817   (  -0.000   -3.122   -3.122)    4.415
   0.871   (  -0.000   -0.918   -0.918)    1.299
   1.078   (   0.000    7.283    7.283)   10.300
   1.150   (   0.000   10.573   10.573)   14.952
   1.399   (  -0.000   -0.818   -0.818)    1.157
   1.467   (  -0.000   -0.404   -0.404)    0.571
   1.527   (   0.000    1.361    1.361)    1.924
   1.687   (  -0.000   -3.339   -3.339)    4.723
   1.787   (  -0.000   -2.020   -2.020)    2.857
   1.933   (   0.000    2.144    2.144)    3.032
   2.254   (   0.000    0.095    0.095)    0.135
   2.265   (  -0.000   -0.386   -0.386)    0.546
   2.277   (   0.000    0.705    0.705)    0.997
   2.333   (  -0.000   -4.321   -4.321)    6.111
   2.390   (   0.000    1.075    1.075)    1.520
   2.511   (  -0.000   -1.697   -1.697)    2.400
   2.519   (   0.000    0.897    0.897)    1.268
   2.565   (   0.000    0.014    0.014)    0.019
   2.617   (  -0.000   -0.080   -0.080)    0.113
   2.673   (   0.000    1.561    1.561)    2.207
   2.756   (   0.000    0.095    0.095)    0.135
   2.815   (  -0.000   -1.163   -1.163)    1.645
   2.874   (   0.000    2.337    2.337)    3.305
   2.895   (   0.000    0.757    0.757)    1.070
   3.004   (   0.000    0.721    0.721)    1.019
   3.020   (  -0.000   -1.328   -1.328)    1.878
   3.081   (   0.000    0.707    0.707)    1.000
   3.206   (  -0.000   -2.171   -2.171)    3.070
   3.216   (  -0.000   -1.482   -1.482)    2.096
   3.278   (  -0.000   -1.209   -1.209)    1.709
   3.331   (   0.000    0.924    0.924)    1.306
   3.376   (   0.000    0.116    0.116)    0.164
   3.394   (  -0.000   -4.475   -4.475)    6.329
   3.525   (   0.000    3.928    3.928)    5.555
   3.546   (   0.000    2.294    2.294)    3.244
   3.591   (   0.000    4.132    4.132)    5.843
   3.797   (   0.000    5.509    5.509)    7.791
   3.866   (  -0.000   -6.121   -6.121)    8.656
   3.878   (   0.000    1.335    1.335)    1.888
   3.895   (  -0.000   -3.734   -3.734)    5.281
   3.990   (  -0.000   -2.651   -2.651)    3.750
   4.049   (  -0.000   -1.749   -1.749)    2.473
   4.076   (  -0.000   -2.975   -2.975)    4.207
   4.106   (  -0.000   -1.039   -1.039)    1.470
   4.136   (  -0.000   -1.624   -1.624)    2.296
   4.242   (  -0.000   -1.026   -1.026)    1.450
   4.252   (   0.000    0.754    0.754)    1.066
   4.277   (   0.000    1.488    1.488)    2.105
   4.340   (  -0.000   -1.945   -1.945)    2.750
   4.360   (  -0.000   -1.153   -1.153)    1.631
   4.428   (   0.000    1.443    1.443)    2.040
   4.471   (  -0.000   -0.395   -0.395)    0.558
   4.477   (   0.000    0.178    0.178)    0.252
   4.525   (   0.000    1.114    1.114)    1.576
   4.565   (   0.000    1.725    1.725)    2.439
   4.597   (   0.000    1.446    1.446)    2.046
   4.610   (  -0.000   -0.034   -0.034)    0.048
   4.676   (  -0.000   -0.750   -0.750)    1.061
   4.705   (  -0.000   -0.005   -0.005)    0.007
   4.875   (  -0.000   -0.394   -0.394)    0.557
   5.186   (  -0.000   -0.002   -0.002)    0.002
   5.235   (   0.000    0.710    0.710)    1.004
   5.237   (   0.000    0.682    0.682)    0.964
   5.313   (  -0.000   -1.276   -1.276)    1.804
======================= Grid point 3 (4/16) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 26
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.521   (   0.000    0.000    0.000)    0.000
   0.555   (   0.000    0.000    0.000)    0.000
   0.576   (   0.000    0.000    0.000)    0.000
   0.734   (   0.000    0.000    0.000)    0.000
   0.762   (   0.000    0.000    0.000)    0.000
   0.873   (   0.000    0.000    0.000)    0.000
   0.901   (   0.000    0.000    0.000)    0.000
   1.198   (   0.000    0.000    0.000)    0.000
   1.326   (   0.000    0.000    0.000)    0.000
   1.381   (   0.000    0.000    0.000)    0.000
   1.445   (   0.000    0.000    0.000)    0.000
   1.546   (   0.000    0.000    0.000)    0.000
   1.642   (   0.000    0.000    0.000)    0.000
   1.737   (   0.000    0.000    0.000)    0.000
   1.979   (  -0.000   -0.000   -0.000)    0.000
   2.225   (   0.000    0.000    0.000)    0.000
   2.256   (   0.000    0.000    0.000)    0.000
   2.281   (  -0.000   -0.000   -0.000)    0.000
   2.292   (   0.000    0.000    0.000)    0.000
   2.454   (   0.000    0.000    0.000)    0.000
   2.481   (   0.000    0.000    0.000)    0.000
   2.541   (   0.000    0.000    0.000)    0.000
   2.556   (  -0.000   -0.000   -0.000)    0.000
   2.618   (   0.000    0.000    0.000)    0.000
   2.704   (   0.000    0.000    0.000)    0.000
   2.709   (   0.000    0.000    0.000)    0.000
   2.794   (   0.000    0.000    0.000)    0.000
   2.904   (   0.000    0.000    0.000)    0.000
   2.918   (   0.000    0.000    0.000)    0.000
   2.992   (   0.000    0.000    0.000)    0.000
   3.023   (   0.000    0.000    0.000)    0.000
   3.119   (   0.000    0.000    0.000)    0.000
   3.149   (   0.000    0.000    0.000)    0.000
   3.188   (   0.000    0.000    0.000)    0.000
   3.272   (   0.000    0.000    0.000)    0.000
   3.291   (   0.000    0.000    0.000)    0.000
   3.349   (   0.000    0.000    0.000)    0.000
   3.373   (   0.000    0.000    0.000)    0.000
   3.619   (   0.000    0.000    0.000)    0.000
   3.639   (  -0.000   -0.000   -0.000)    0.000
   3.667   (   0.000    0.000    0.000)    0.000
   3.688   (   0.000    0.000    0.000)    0.000
   3.839   (   0.000    0.000    0.000)    0.000
   3.893   (   0.000    0.000    0.000)    0.000
   3.925   (   0.000    0.000    0.000)    0.000
   3.959   (   0.000    0.000    0.000)    0.000
   4.004   (   0.000    0.000    0.000)    0.000
   4.034   (   0.000    0.000    0.000)    0.000
   4.076   (   0.000    0.000    0.000)    0.000
   4.105   (   0.000    0.000    0.000)    0.000
   4.210   (   0.000    0.000    0.000)    0.000
   4.278   (  -0.000   -0.000   -0.000)    0.000
   4.308   (  -0.000   -0.000   -0.000)    0.000
   4.309   (   0.000    0.000    0.000)    0.000
   4.316   (   0.000    0.000    0.000)    0.000
   4.456   (  -0.000   -0.000   -0.000)    0.000
   4.461   (   0.000    0.000    0.000)    0.000
   4.468   (   0.000    0.000    0.000)    0.000
   4.545   (   0.000    0.000    0.000)    0.000
   4.597   (   0.000    0.000    0.000)    0.000
   4.606   (   0.000    0.000    0.000)    0.000
   4.633   (   0.000    0.000    0.000)    0.000
   4.656   (  -0.000   -0.000   -0.000)    0.000
   4.705   (   0.000    0.000    0.000)    0.000
   4.869   (   0.000    0.000    0.000)    0.000
   5.188   (   0.000    0.000    0.000)    0.000
   5.246   (   0.000    0.000    0.000)    0.000
   5.253   (   0.000    0.000    0.000)    0.000
   5.286   (   0.000    0.000    0.000)    0.000
======================= Grid point 8 (5/16) =======================
q-point: ( 0.17  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 66
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.422   (   1.219    1.219   10.658)   10.797
   0.490   (   0.982    0.982    8.035)    8.154
   0.665   (  -2.800   -2.800    1.832)    4.363
   0.724   (  -0.082   -0.082   -1.703)    1.707
   0.764   (  -1.766   -1.766   -0.157)    2.502
   0.858   (  -1.351   -1.351   -2.052)    2.804
   0.917   (   2.454    2.454   -3.049)    4.619
   1.017   (   2.295    2.295    4.115)    5.241
   1.025   (   9.924    9.924   13.026)   19.148
   1.418   (   1.293    1.293   -1.225)    2.201
   1.452   (  -0.847   -0.847    1.285)    1.757
   1.534   (  -2.644   -2.644    6.581)    7.569
   1.718   (  -3.103   -3.103   -1.864)    4.768
   1.835   (   1.410    1.410   -1.926)    2.773
   1.897   (   0.667    0.667    1.927)    2.145
   2.254   (   0.235    0.235   -0.768)    0.837
   2.262   (   0.358    0.358   -0.069)    0.511
   2.274   (  -0.050   -0.050    0.686)    0.690
   2.355   (  -1.060   -1.060   -2.732)    3.116
   2.382   (  -4.152   -4.152   -2.494)    6.380
   2.501   (   0.157    0.157   -2.442)    2.452
   2.539   (  -0.599   -0.599   -1.419)    1.652
   2.542   (   0.005    0.005   -1.686)    1.686
   2.615   (  -0.014   -0.014    3.245)    3.245
   2.694   (   2.780    2.780    0.291)    3.943
   2.755   (  -0.279   -0.279    1.785)    1.828
   2.791   (   0.968    0.968    0.251)    1.391
   2.881   (   2.271    2.271    1.804)    3.683
   2.887   (   1.674    1.674    1.087)    2.605
   3.002   (   0.633    0.633   -1.693)    1.915
   3.036   (   0.016    0.016    0.872)    0.872
   3.076   (  -0.201   -0.201   -1.509)    1.535
   3.208   (  -0.121   -0.121   -3.525)    3.529
   3.257   (  -0.909   -0.909   -1.584)    2.040
   3.290   (  -0.392   -0.392   -1.727)    1.813
   3.293   (  -0.872   -0.872    1.639)    2.051
   3.400   (  -2.363   -2.363    1.502)    3.664
   3.401   (  -2.184   -2.184   -0.431)    3.118
   3.504   (   0.954    0.954    0.715)    1.527
   3.524   (   0.885    0.885    1.541)    1.985
   3.560   (   3.474    3.474    5.639)    7.480
   3.752   (   5.126    5.126    7.037)   10.103
   3.872   (  -1.877   -1.877   -2.353)    3.547
   3.890   (  -3.253   -3.253   -4.338)    6.323
   3.927   (  -6.041   -6.041   -3.507)    9.235
   4.028   (  -1.903   -1.903   -4.075)    4.884
   4.073   (   0.401    0.401   -2.098)    2.173
   4.086   (   1.032    1.032   -4.757)    4.976
   4.137   (   0.735    0.735   -2.644)    2.841
   4.175   (   1.290    1.290   -2.500)    3.095
   4.218   (   1.740    1.740   -1.826)    3.065
   4.226   (   0.478    0.478   -2.226)    2.327
   4.299   (  -1.110   -1.110    2.534)    2.981
   4.356   (  -1.459   -1.459    4.388)    4.849
   4.367   (  -0.255   -0.255   -2.315)    2.343
   4.408   (  -0.780   -0.780    0.890)    1.417
   4.446   (  -1.064   -1.064    1.653)    2.235
   4.488   (   0.942    0.942   -0.800)    1.554
   4.521   (   0.070    0.070    0.656)    0.663
   4.545   (   1.552    1.552    2.197)    3.105
   4.598   (  -0.189   -0.189    1.725)    1.745
   4.602   (  -0.482   -0.482    0.851)    1.090
   4.684   (   0.077    0.077   -0.962)    0.968
   4.703   (  -0.434   -0.434    0.220)    0.652
   4.877   (  -0.240   -0.240   -0.874)    0.938
   5.202   (   0.897    0.897   -0.629)    1.416
   5.203   (   0.988    0.988   -0.306)    1.431
   5.243   (  -0.457   -0.457    2.017)    2.118
   5.324   (  -0.606   -0.606   -1.444)    1.679
======================= Grid point 9 (6/16) =======================
q-point: ( 0.33  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.499   (  -0.794   -0.465    0.309)    0.971
   0.555   (  -1.636   -0.291    1.931)    2.548
   0.640   (  -0.957   -4.461    2.693)    5.298
   0.687   (  -2.170   -2.208    0.305)    3.111
   0.831   (   4.633    1.346    1.491)    5.050
   0.871   (   0.828    0.764    4.904)    5.032
   0.919   (   3.287    1.143    1.665)    3.858
   1.116   (  -1.046    5.286   -1.911)    5.717
   1.215   (  -4.915    3.019   -0.608)    5.800
   1.416   (  -0.461   -2.495    2.622)    3.648
   1.447   (  -2.406   -2.025    0.741)    3.231
   1.568   (  -3.119   -4.719    5.833)    8.126
   1.720   (   5.884    0.565    3.847)    7.053
   1.801   (   4.340   -2.232   -1.057)    4.994
   1.961   (   0.298    2.809    1.782)    3.340
   2.229   (  -1.357   -3.148   -2.095)    4.017
   2.268   (   0.144   -0.642    0.696)    0.958
   2.277   (  -0.492    0.931   -0.507)    1.169
   2.295   (   0.294   -0.536    0.175)    0.636
   2.351   (  -6.089    2.670    0.376)    6.660
   2.454   (  -1.669    0.188   -4.225)    4.546
   2.513   (  -2.150   -1.188    0.917)    2.622
   2.527   (  -2.847   -1.001   -0.272)    3.030
   2.657   (   3.302    1.112   -0.011)    3.484
   2.723   (   2.982    1.381   -0.388)    3.309
   2.784   (   2.895    1.539    0.862)    3.390
   2.809   (   3.766    0.262   -0.420)    3.799
   2.901   (   0.767    0.905   -0.428)    1.261
   2.910   (   0.300    0.140   -0.354)    0.485
   3.002   (   1.787    1.446   -1.113)    2.554
   3.043   (   3.048    0.102    0.241)    3.059
   3.097   (   0.255    3.068    0.585)    3.134
   3.149   (   0.244    0.404   -3.691)    3.722
   3.206   (   0.933   -1.402   -2.749)    3.224
   3.248   (  -1.745   -0.262   -1.487)    2.307
   3.275   (  -2.951   -0.112   -3.912)    4.902
   3.353   (  -2.095   -2.769    2.155)    4.087
   3.373   (  -2.993   -3.009    1.161)    4.401
   3.568   (  -2.146    2.289    2.857)    4.244
   3.593   (  -1.319    2.458    5.259)    5.953
   3.673   (   4.175    0.426    3.146)    5.245
   3.747   (  -3.494   -6.566   -4.278)    8.580
   3.831   (  -1.950    0.013   -0.636)    2.051
   3.872   (  -0.121    0.909    0.734)    1.175
   3.933   (  -1.572   -2.242   -1.921)    3.345
   3.944   (   2.922    4.498    4.718)    7.144
   4.040   (   1.827    1.164   -2.390)    3.226
   4.056   (   1.295   -0.297   -1.723)    2.175
   4.082   (   2.276   -1.950   -2.567)    3.946
   4.157   (   4.273    0.194   -1.237)    4.453
   4.195   (  -1.518    1.117   -2.824)    3.395
   4.233   (   0.397    3.509   -1.840)    3.982
   4.275   (   0.729    1.631   -4.792)    5.114
   4.337   (   0.332   -2.600    2.461)    3.596
   4.373   (  -0.535   -3.378    1.207)    3.627
   4.428   (  -2.629    0.767    2.979)    4.047
   4.445   (  -2.462   -0.233   -0.039)    2.474
   4.477   (  -0.051   -1.803   -0.191)    1.814
   4.520   (  -2.225    0.159   -0.212)    2.240
   4.597   (   1.044    1.074    1.359)    2.023
   4.613   (   0.481    0.001    0.364)    0.603
   4.630   (  -0.020    0.020    1.777)    1.777
   4.669   (   0.845   -0.242   -0.732)    1.144
   4.699   (  -0.414    0.151   -0.379)    0.581
   4.864   (  -1.034   -0.591    0.154)    1.201
   5.199   (   1.641    0.245   -0.259)    1.680
   5.228   (   0.093    2.053   -0.310)    2.078
   5.253   (  -0.206   -0.385   -0.189)    0.476
   5.297   (  -0.228   -0.929   -0.065)    0.959
======================= Grid point 15 (7/16) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 32
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.458   (  -0.349   -0.349    0.349)    0.605
   0.589   (  -3.023   -3.023    0.087)    4.277
   0.589   (  -1.066   -1.066    4.002)    4.277
   0.734   (  -2.616   -2.616    2.616)    4.531
   0.929   (   1.648    1.648   -3.649)    4.330
   0.929   (   2.982    2.982   -0.981)    4.330
   1.048   (  -6.456   -6.456   -3.852)    9.910
   1.048   (   0.416    0.416    9.892)    9.910
   1.049   (   3.126    3.126   -3.126)    5.414
   1.433   (  -2.055   -2.055    1.864)    3.453
   1.433   (  -1.928   -1.928    2.119)    3.453
   1.595   (  -6.358   -6.358    6.358)   11.013
   1.835   (   1.010    1.010   -2.259)    2.673
   1.835   (   1.843    1.843   -0.594)    2.673
   1.983   (   0.959    0.959   -0.959)    1.661
   2.230   (  -3.591   -3.591   -0.253)    5.085
   2.230   (  -1.028   -1.028    4.872)    5.085
   2.286   (  -0.546   -0.546   -0.480)    0.909
   2.286   (   0.138    0.138    0.888)    0.909
   2.306   (   0.290    0.290   -0.290)    0.503
   2.391   (   1.928    1.928   -1.928)    3.340
   2.499   (  -1.793   -1.793    1.792)    3.105
   2.499   (  -1.792   -1.792    1.793)    3.105
   2.708   (   1.772    1.772   -4.535)    5.181
   2.708   (   3.614    3.614   -0.850)    5.181
   2.842   (  -0.264   -0.264    0.264)    0.457
   2.874   (   1.415    1.415   -2.009)    2.836
   2.874   (   1.811    1.811   -1.217)    2.836
   2.901   (   1.628    1.628   -1.628)    2.820
   3.030   (   1.786    1.786   -1.786)    3.094
   3.093   (   1.272    1.272   -3.335)    3.790
   3.093   (   2.647    2.647   -0.585)    3.790
   3.111   (   1.804    1.804   -1.804)    3.124
   3.185   (  -0.033   -0.033   -0.254)    0.258
   3.185   (   0.158    0.158    0.129)    0.258
   3.216   (   0.847    0.847   -0.847)    1.467
   3.353   (  -3.169   -3.169    2.938)    5.359
   3.353   (  -3.015   -3.015    3.247)    5.359
   3.615   (  -2.887   -2.887   -0.725)    4.147
   3.615   (  -0.479   -0.479    4.091)    4.147
   3.658   (  -4.710   -4.710    4.710)    8.158
   3.781   (  -0.082   -0.082    0.082)    0.142
   3.810   (   0.814    0.814   -2.364)    2.630
   3.810   (   1.847    1.847   -0.298)    2.630
   3.917   (  -1.635   -1.635    1.635)    2.832
   4.038   (   1.201    1.201   -1.796)    2.472
   4.038   (   1.598    1.598   -1.003)    2.472
   4.085   (   1.047    1.047   -1.047)    1.814
   4.087   (   0.987    0.987   -3.327)    3.608
   4.087   (   2.547    2.547   -0.207)    3.608
   4.224   (   2.351    2.351   -3.433)    4.779
   4.224   (   3.072    3.072   -1.991)    4.779
   4.226   (   1.209    1.209   -1.209)    2.094
   4.353   (  -1.964   -1.964    0.777)    2.885
   4.353   (  -1.173   -1.173    2.360)    2.885
   4.437   (  -0.924   -0.924    0.924)    1.601
   4.465   (  -1.817   -1.817    0.467)    2.612
   4.465   (  -0.917   -0.917    2.267)    2.612
   4.473   (  -1.359   -1.359    1.359)    2.353
   4.619   (   0.196    0.196   -0.490)    0.563
   4.619   (   0.392    0.392   -0.097)    0.563
   4.643   (  -0.958   -0.958    0.958)    1.659
   4.683   (   0.090    0.090   -1.213)    1.220
   4.683   (   0.839    0.839    0.285)    1.220
   4.856   (  -1.194   -1.194    1.194)    2.067
   5.216   (   0.985    0.985   -0.985)    1.705
   5.233   (   0.256    0.256   -1.334)    1.382
   5.233   (   0.975    0.975    0.103)    1.382
   5.299   (  -0.754   -0.754    0.754)    1.307
======================= Grid point 50 (8/16) =======================
q-point: ( 0.17  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 28
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.515   (   0.444    0.444    0.444)    0.769
   0.515   (   0.444    0.444    0.444)    0.769
   0.582   (  -2.139   -2.139   -2.139)    3.706
   0.719   (   1.402    1.402    1.402)    2.429
   0.719   (   1.402    1.402    1.402)    2.429
   0.877   (   2.762    2.762    2.762)    4.784
   1.053   (   2.530    2.530    2.530)    4.382
   1.053   (   2.530    2.530    2.530)    4.382
   1.166   (  -2.273   -2.273   -2.273)    3.936
   1.433   (  -1.035   -1.035   -1.035)    1.792
   1.433   (  -1.035   -1.035   -1.035)    1.792
   1.466   (  -0.775   -0.775   -0.775)    1.342
   1.733   (   4.863    4.863    4.863)    8.422
   1.898   (   1.362    1.362    1.362)    2.359
   1.898   (   1.362    1.362    1.362)    2.359
   2.247   (  -2.314   -2.314   -2.314)    4.008
   2.247   (  -2.314   -2.314   -2.314)    4.008
   2.279   (   0.639    0.639    0.639)    1.107
   2.300   (  -0.723   -0.723   -0.723)    1.252
   2.300   (  -0.723   -0.723   -0.723)    1.252
   2.460   (  -1.416   -1.416   -1.416)    2.453
   2.498   (  -2.089   -2.089   -2.089)    3.618
   2.498   (  -2.089   -2.089   -2.089)    3.618
   2.683   (   3.060    3.060    3.060)    5.299
   2.774   (   1.678    1.678    1.678)    2.906
   2.774   (   1.678    1.678    1.678)    2.906
   2.846   (   2.296    2.296    2.296)    3.976
   2.918   (   0.475    0.475    0.475)    0.823
   2.918   (   0.475    0.475    0.475)    0.823
   3.061   (   0.069    0.069    0.069)    0.120
   3.069   (   2.224    2.224    2.224)    3.852
   3.069   (   2.224    2.224    2.224)    3.852
   3.206   (  -1.169   -1.169   -1.169)    2.025
   3.206   (  -1.169   -1.169   -1.169)    2.025
   3.260   (  -0.856   -0.856   -0.856)    1.483
   3.262   (  -1.812   -1.812   -1.812)    3.139
   3.326   (  -2.788   -2.788   -2.788)    4.830
   3.326   (  -2.788   -2.788   -2.788)    4.830
   3.528   (   0.473    0.473    0.473)    0.819
   3.528   (   0.473    0.473    0.473)    0.819
   3.686   (   4.154    4.154    4.154)    7.194
   3.751   (  -5.437   -5.437   -5.437)    9.418
   3.837   (  -0.024   -0.024   -0.024)    0.042
   3.837   (  -0.024   -0.024   -0.024)    0.042
   3.931   (   6.087    6.087    6.087)   10.543
   3.943   (  -2.989   -2.989   -2.989)    5.177
   4.071   (  -0.938   -0.938   -0.938)    1.625
   4.071   (  -0.938   -0.938   -0.938)    1.625
   4.179   (   0.303    0.303    0.303)    0.524
   4.179   (   0.303    0.303    0.303)    0.524
   4.209   (  -1.451   -1.451   -1.451)    2.514
   4.305   (   1.653    1.653    1.653)    2.863
   4.305   (   1.653    1.653    1.653)    2.863
   4.321   (   0.739    0.739    0.739)    1.280
   4.354   (  -1.904   -1.904   -1.904)    3.297
   4.354   (  -1.904   -1.904   -1.904)    3.297
   4.402   (  -0.577   -0.577   -0.577)    0.999
   4.491   (  -1.284   -1.284   -1.284)    2.225
   4.491   (  -1.284   -1.284   -1.284)    2.225
   4.594   (   1.500    1.500    1.500)    2.598
   4.611   (   0.707    0.707    0.707)    1.224
   4.611   (   0.707    0.707    0.707)    1.224
   4.691   (  -0.036   -0.036   -0.036)    0.063
   4.691   (  -0.036   -0.036   -0.036)    0.063
   4.855   (  -0.965   -0.965   -0.965)    1.672
   5.229   (   0.481    0.481    0.481)    0.833
   5.233   (   0.515    0.515    0.515)    0.892
   5.233   (   0.515    0.515    0.515)    0.892
   5.300   (  -0.862   -0.862   -0.862)    1.493
======================= Grid point 51 (9/16) =======================
q-point: ( 0.33  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 39
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.491   (  -2.598   -0.000   -0.000)    2.598
   0.493   (  -0.635   -0.000   -0.000)    0.635
   0.559   (  -6.135   -0.000   -0.000)    6.135
   0.716   (   8.684   -0.000   -0.000)    8.684
   0.875   (  -1.050   -0.000   -0.000)    1.050
   0.984   (   9.976   -0.000   -0.000)    9.976
   0.994   (   5.358    0.000    0.000)    5.358
   1.033   ( -10.692   -0.000   -0.000)   10.692
   1.159   (  -5.352   -0.000   -0.000)    5.352
   1.396   (   0.552   -0.000   -0.000)    0.552
   1.408   (  -1.283   -0.000   -0.000)    1.283
   1.463   (  -5.074   -0.000   -0.000)    5.074
   1.862   (   7.713    0.000    0.000)    7.713
   1.872   (   7.223   -0.000   -0.000)    7.223
   1.989   (  -0.920    0.000    0.000)    0.920
   2.141   (  -2.996   -0.000   -0.000)    2.996
   2.216   (  -9.502   -0.000   -0.000)    9.502
   2.287   (   1.859   -0.000   -0.000)    1.859
   2.301   (   1.172   -0.000   -0.000)    1.172
   2.308   (  -0.211   -0.000   -0.000)    0.211
   2.440   (  -0.981   -0.000   -0.000)    0.981
   2.454   (  -2.880   -0.000   -0.000)    2.880
   2.456   (  -2.793   -0.000   -0.000)    2.793
   2.744   (   4.200   -0.000   -0.000)    4.200
   2.778   (   2.265   -0.000   -0.000)    2.265
   2.821   (  -0.334   -0.000   -0.000)    0.334
   2.897   (   5.311    0.000    0.000)    5.311
   2.912   (   1.593   -0.000   -0.000)    1.593
   2.934   (   1.202   -0.000   -0.000)    1.202
   3.060   (   4.595   -0.000   -0.000)    4.595
   3.108   (   4.770   -0.000   -0.000)    4.770
   3.142   (   1.699   -0.000   -0.000)    1.699
   3.177   (  -0.038   -0.000   -0.000)    0.038
   3.177   (   0.134   -0.000   -0.000)    0.134
   3.222   (  -2.770    0.000    0.000)    2.770
   3.223   (  -2.913   -0.000   -0.000)    2.913
   3.266   (  -5.120   -0.000   -0.000)    5.120
   3.280   (  -5.664   -0.000   -0.000)    5.664
   3.543   (  -0.915   -0.000   -0.000)    0.915
   3.591   (  -3.600   -0.000   -0.000)    3.600
   3.610   (  -3.264    0.000    0.000)    3.264
   3.757   (   5.504   -0.000   -0.000)    5.504
   3.850   (   0.381   -0.000   -0.000)    0.381
   3.866   (  -1.647   -0.000   -0.000)    1.647
   3.869   (  -1.175   -0.000   -0.000)    1.175
   4.045   (  -1.205   -0.000   -0.000)    1.205
   4.049   (  -0.996   -0.000   -0.000)    0.996
   4.075   (   4.833   -0.000   -0.000)    4.833
   4.123   (   5.881   -0.000   -0.000)    5.881
   4.186   (  -0.060   -0.000   -0.000)    0.060
   4.227   (   1.563   -0.000   -0.000)    1.563
   4.270   (   0.009    0.000    0.000)    0.009
   4.311   (   3.063   -0.000   -0.000)    3.063
   4.311   (  -2.435   -0.000   -0.000)    2.435
   4.348   (   2.404   -0.000   -0.000)    2.404
   4.400   (  -1.651   -0.000   -0.000)    1.651
   4.404   (  -1.961   -0.000   -0.000)    1.961
   4.435   (  -2.363    0.000    0.000)    2.363
   4.467   (  -4.410   -0.000   -0.000)    4.410
   4.620   (   0.745   -0.000   -0.000)    0.745
   4.623   (   0.591   -0.000   -0.000)    0.591
   4.625   (  -0.031   -0.000   -0.000)    0.031
   4.685   (   0.972   -0.000   -0.000)    0.972
   4.696   (  -0.366   -0.000   -0.000)    0.366
   4.837   (  -1.411   -0.000   -0.000)    1.411
   5.229   (   1.852   -0.000   -0.000)    1.852
   5.246   (   0.003   -0.000   -0.000)    0.003
   5.248   (  -0.177   -0.000   -0.000)    0.177
   5.282   (  -0.453   -0.000   -0.000)    0.453
======================= Grid point 55 (10/16) =======================
q-point: (-0.17  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 28
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.473   (  10.919    0.000    0.000)   10.919
   0.473   (  10.919    0.000    0.000)   10.919
   0.718   (   2.186    0.000    0.000)    2.186
   0.718   (   2.186    0.000    0.000)    2.186
   0.755   (  21.130    0.000    0.000)   21.130
   0.793   (  -1.573   -0.000   -0.000)    1.573
   0.917   (  -4.191   -0.000   -0.000)    4.191
   0.917   (  -4.191   -0.000   -0.000)    4.191
   0.947   (   1.899    0.000    0.000)    1.899
   1.392   (   0.911    0.000    0.000)    0.911
   1.472   (   3.767    0.000    0.000)    3.767
   1.586   (   8.902    0.000    0.000)    8.902
   1.805   (  -2.131   -0.000   -0.000)    2.131
   1.805   (  -2.131   -0.000   -0.000)    2.131
   1.862   (   0.976    0.000    0.000)    0.976
   2.255   (  -0.300   -0.000   -0.000)    0.300
   2.255   (  -0.300   -0.000   -0.000)    0.300
   2.274   (   0.512    0.000    0.000)    0.512
   2.365   (  -4.161   -0.000   -0.000)    4.161
   2.479   (  -0.203   -0.000   -0.000)    0.203
   2.479   (  -0.203   -0.000   -0.000)    0.203
   2.511   (  -6.490   -0.000   -0.000)    6.490
   2.544   (  -3.440   -0.000   -0.000)    3.440
   2.635   (   1.704    0.000    0.000)    1.704
   2.635   (   1.704    0.000    0.000)    1.704
   2.795   (  -0.297   -0.000   -0.000)    0.297
   2.795   (  -0.297   -0.000   -0.000)    0.297
   2.809   (   5.265    0.000    0.000)    5.265
   2.836   (   2.516    0.000    0.000)    2.516
   3.043   (  -0.244   -0.000   -0.000)    0.244
   3.043   (  -0.244   -0.000   -0.000)    0.244
   3.047   (  -5.782   -0.000   -0.000)    5.782
   3.217   (  -2.229   -0.000   -0.000)    2.229
   3.265   (  -4.885   -0.000   -0.000)    4.885
   3.296   (   1.838    0.000    0.000)    1.838
   3.296   (   1.838    0.000    0.000)    1.838
   3.460   (   3.498    0.000    0.000)    3.498
   3.460   (   3.498    0.000    0.000)    3.498
   3.460   (  -1.390   -0.000   -0.000)    1.390
   3.494   (   0.005    0.000    0.000)    0.005
   3.510   (   6.507    0.000    0.000)    6.507
   3.662   (   8.615    0.000    0.000)    8.615
   3.950   (  -5.803   -0.000   -0.000)    5.803
   3.950   (  -5.803   -0.000   -0.000)    5.803
   3.975   (  -7.017   -0.000   -0.000)    7.017
   4.042   (  -3.536   -0.000   -0.000)    3.536
   4.078   (  -3.998   -0.000   -0.000)    3.998
   4.078   (  -3.998   -0.000   -0.000)    3.998
   4.157   (  -3.229   -0.000   -0.000)    3.229
   4.159   (  -0.540   -0.000   -0.000)    0.540
   4.186   (  -2.741   -0.000   -0.000)    2.741
   4.186   (  -2.741   -0.000   -0.000)    2.741
   4.313   (   0.098    0.000    0.000)    0.098
   4.398   (   2.741    0.000    0.000)    2.741
   4.399   (   3.732    0.000    0.000)    3.732
   4.399   (   3.732    0.000    0.000)    3.732
   4.426   (   0.889    0.000    0.000)    0.889
   4.480   (   1.777    0.000    0.000)    1.777
   4.511   (   0.968    0.000    0.000)    0.968
   4.511   (   0.968    0.000    0.000)    0.968
   4.617   (   0.011    0.000    0.000)    0.011
   4.617   (   0.011    0.000    0.000)    0.011
   4.684   (  -0.668   -0.000   -0.000)    0.668
   4.712   (   0.356    0.000    0.000)    0.356
   4.873   (  -0.775   -0.000   -0.000)    0.775
   5.186   (  -0.427   -0.000   -0.000)    0.427
   5.186   (  -0.427   -0.000   -0.000)    0.427
   5.252   (   1.863    0.000    0.000)    1.863
   5.333   (  -1.665   -0.000   -0.000)    1.665
======================= Grid point 62 (11/16) =======================
q-point: (-0.33  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.528   (   2.808    2.354    0.000)    3.664
   0.562   (  -0.898   -1.329    0.000)    1.603
   0.715   (  -3.717    0.652   -0.000)    3.774
   0.747   (   8.195    6.236    0.000)   10.298
   0.801   (   3.276   -0.039    0.000)    3.276
   0.853   (   1.035    1.709   -0.000)    1.998
   0.904   (   5.549   -0.154    0.000)    5.551
   1.031   (   0.143   -4.733    0.000)    4.735
   1.116   (   7.558   -9.541    0.000)   12.171
   1.444   (   6.032    0.857    0.000)    6.093
   1.538   (   7.637    3.675    0.000)    8.475
   1.668   (   7.162    4.578    0.000)    8.501
   1.725   (  -1.007    5.158   -0.000)    5.255
   1.795   (   0.572    0.138    0.000)    0.588
   1.911   (   0.888   -3.298    0.000)    3.415
   2.248   (  -0.268    0.399   -0.000)    0.480
   2.258   (   0.085   -0.783    0.000)    0.787
   2.283   (   0.620    0.028    0.000)    0.620
   2.323   (  -2.315   -0.637   -0.000)    2.400
   2.375   (  -2.771    3.055   -0.000)    4.124
   2.469   (  -2.754   -0.096   -0.000)    2.756
   2.499   (  -1.451   -2.664    0.000)    3.034
   2.517   (  -2.906   -2.557   -0.000)    3.871
   2.666   (   2.392    0.692    0.000)    2.490
   2.681   (   2.501   -0.787    0.000)    2.622
   2.749   (  -1.222    1.131   -0.000)    1.664
   2.782   (  -2.261   -0.888   -0.000)    2.429
   2.874   (   0.984   -0.815    0.000)    1.278
   2.897   (   1.378   -1.449    0.000)    1.999
   2.966   (  -2.007   -0.717   -0.000)    2.131
   3.024   (  -1.609    0.360   -0.000)    1.649
   3.053   (  -1.622   -2.589    0.000)    3.055
   3.161   (  -3.316    1.744   -0.000)    3.747
   3.200   (  -4.507   -0.697   -0.000)    4.561
   3.282   (  -3.084    2.207   -0.000)    3.792
   3.307   (  -1.916    0.066   -0.000)    1.917
   3.396   (  -1.895    1.314   -0.000)    2.306
   3.444   (  -0.161    3.101   -0.000)    3.105
   3.547   (   5.137   -1.765    0.000)    5.432
   3.560   (   5.231   -2.599    0.000)    5.841
   3.624   (   5.998   -1.110    0.000)    6.100
   3.812   (   5.694   -3.041    0.000)    6.455
   3.852   (  -3.609    3.960   -0.000)    5.358
   3.877   (  -0.039   -0.349    0.000)    0.351
   3.926   (  -2.375    1.826   -0.000)    2.996
   3.970   (  -2.780    0.830   -0.000)    2.901
   4.021   (  -2.132   -0.458   -0.000)    2.181
   4.028   (  -1.965   -0.384   -0.000)    2.003
   4.086   (  -3.465    1.835   -0.000)    3.921
   4.124   (  -2.103    0.453   -0.000)    2.151
   4.167   (  -2.885   -2.571   -0.000)    3.865
   4.190   (  -2.958   -3.660    0.000)    4.706
   4.262   (  -4.742    0.712   -0.000)    4.795
   4.386   (   2.935    3.967   -0.000)    4.935
   4.417   (   0.934    1.852   -0.000)    2.074
   4.448   (   2.024   -0.517    0.000)    2.089
   4.471   (   0.453   -0.508    0.000)    0.681
   4.487   (   0.679    1.354   -0.000)    1.514
   4.522   (   0.188   -0.962    0.000)    0.980
   4.566   (   2.286   -2.075    0.000)    3.087
   4.611   (   0.473    0.588   -0.000)    0.755
   4.623   (   0.552   -0.417    0.000)    0.692
   4.670   (  -1.028    0.080   -0.000)    1.031
   4.706   (  -0.547    0.475   -0.000)    0.724
   4.877   (   1.078    0.242    0.000)    1.105
   5.181   (  -0.282   -0.188   -0.000)    0.339
   5.195   (  -0.379   -1.837    0.000)    1.876
   5.262   (  -0.596    0.444   -0.000)    0.744
   5.311   (  -0.081    0.562   -0.000)    0.568
======================= Grid point 67 (12/16) =======================
q-point: (-0.50  0.50  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 39
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.524   (   0.078    0.078   -0.000)    0.110
   0.552   (  -1.225   -1.225   -0.000)    1.732
   0.683   (  -3.029   -3.029   -0.000)    4.284
   0.764   (  10.697   10.697   -0.000)   15.128
   0.805   (   2.472    2.472   -0.000)    3.495
   0.873   (   1.075    1.075   -0.000)    1.520
   0.987   (   3.816    3.816   -0.000)    5.396
   1.057   (  -5.748   -5.748   -0.000)    8.128
   1.234   (  -2.894   -2.894   -0.000)    4.093
   1.501   (   4.042    4.042   -0.000)    5.716
   1.551   (   7.708    7.708   -0.000)   10.901
   1.668   (   3.463    3.463   -0.000)    4.897
   1.725   (   6.645    6.645    0.000)    9.397
   1.775   (   2.310    2.310   -0.000)    3.267
   1.955   (  -1.515   -1.515   -0.000)    2.142
   2.228   (   0.568    0.568   -0.000)    0.803
   2.272   (  -0.416   -0.416   -0.000)    0.589
   2.286   (   1.361    1.361   -0.000)    1.924
   2.294   (   0.059    0.059    0.000)    0.083
   2.393   (  -3.189   -3.189   -0.000)    4.510
   2.441   (  -1.288   -1.288   -0.000)    1.822
   2.481   (  -4.451   -4.451   -0.000)    6.295
   2.524   (  -2.442   -2.442   -0.000)    3.453
   2.673   (   3.438    3.438   -0.000)    4.862
   2.707   (   0.649    0.649   -0.000)    0.918
   2.748   (  -1.760   -1.760    0.000)    2.489
   2.758   (   1.832    1.832   -0.000)    2.591
   2.882   (  -1.170   -1.170   -0.000)    1.654
   2.907   (  -0.924   -0.924   -0.000)    1.306
   2.974   (  -1.164   -1.164   -0.000)    1.647
   3.012   (  -1.115   -1.115   -0.000)    1.577
   3.089   (  -1.719   -1.719    0.000)    2.431
   3.091   (  -3.399   -3.399   -0.000)    4.806
   3.152   (  -1.170   -1.170   -0.000)    1.655
   3.212   (  -3.170   -3.170   -0.000)    4.484
   3.286   (  -0.462   -0.462   -0.000)    0.654
   3.381   (  -0.379   -0.379   -0.000)    0.536
   3.414   (   0.519    0.519   -0.000)    0.733
   3.623   (   0.691    0.691   -0.000)    0.977
   3.674   (   2.960    2.960   -0.000)    4.186
   3.704   (   2.196    2.196    0.000)    3.106
   3.712   (   1.949    1.949   -0.000)    2.756
   3.855   (  -1.485   -1.485    0.000)    2.100
   3.858   (   0.937    0.937   -0.000)    1.325
   3.958   (   1.288    1.288   -0.000)    1.822
   3.968   (   0.285    0.285    0.000)    0.403
   3.983   (  -0.374   -0.374   -0.000)    0.530
   4.033   (  -2.179   -2.179   -0.000)    3.081
   4.035   (  -0.178   -0.178   -0.000)    0.252
   4.095   (  -0.991   -0.991   -0.000)    1.402
   4.167   (  -3.558   -3.558   -0.000)    5.031
   4.189   (  -3.592   -3.592   -0.000)    5.080
   4.212   (  -0.733   -0.733   -0.000)    1.036
   4.377   (   2.579    2.579    0.000)    3.647
   4.399   (   1.579    1.579   -0.000)    2.233
   4.466   (   0.350    0.350   -0.000)    0.496
   4.473   (   0.451    0.451   -0.000)    0.637
   4.482   (   0.877    0.877   -0.000)    1.240
   4.535   (  -0.724   -0.724   -0.000)    1.024
   4.595   (   0.041    0.041   -0.000)    0.058
   4.617   (   0.485    0.485   -0.000)    0.686
   4.647   (  -0.069   -0.069   -0.000)    0.098
   4.648   (  -0.326   -0.326    0.000)    0.461
   4.695   (  -0.633   -0.633   -0.000)    0.895
   4.885   (   1.355    1.355   -0.000)    1.916
   5.181   (  -0.275   -0.275   -0.000)    0.388
   5.219   (  -1.220   -1.220   -0.000)    1.726
   5.236   (  -1.033   -1.033   -0.000)    1.461
   5.310   (   1.076    1.076   -0.000)    1.522
======================= Grid point 69 (13/16) =======================
q-point: (-0.33  0.50  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.497   (  -1.203    1.777   -1.777)    2.787
   0.574   (  -0.057   -1.864    1.864)    2.637
   0.638   (  -4.197    0.044   -0.044)    4.197
   0.859   (   3.653   -0.911    0.911)    3.873
   0.862   (   6.262    2.183   -2.183)    6.982
   0.885   (   3.619    2.478   -2.478)    5.037
   0.989   (  -1.863   -0.623    0.623)    2.061
   1.036   (   6.693   -0.776    0.776)    6.782
   1.174   (  -1.293    0.074   -0.074)    1.298
   1.530   (   7.002    4.013   -4.013)    9.013
   1.571   (   4.651    7.166   -7.166)   11.151
   1.744   (   8.044    2.516   -2.516)    8.796
   1.788   (   5.301   -1.035    1.035)    5.500
   1.832   (   1.963    1.029   -1.029)    2.444
   1.932   (  -0.377   -1.961    1.961)    2.799
   2.248   (   0.681   -0.462    0.462)    0.944
   2.256   (  -0.013    0.504   -0.504)    0.713
   2.294   (   0.529   -0.234    0.234)    0.624
   2.308   (  -1.311   -1.332    1.332)    2.295
   2.323   (   0.596    2.208   -2.208)    3.178
   2.392   (  -4.845    2.199   -2.199)    5.757
   2.485   (  -3.254   -2.947    2.947)    5.287
   2.491   (  -3.270   -2.645    2.645)    4.969
   2.690   (   0.255    0.556   -0.556)    0.827
   2.693   (   2.285    3.439   -3.439)    5.373
   2.769   (   1.982   -3.469    3.469)    5.291
   2.802   (   1.819   -2.055    2.055)    3.429
   2.872   (  -1.425   -0.258    0.258)    1.471
   2.891   (  -3.421    0.861   -0.861)    3.631
   2.964   (  -0.042   -2.154    2.154)    3.046
   3.023   (  -1.219   -2.313    2.313)    3.490
   3.026   (  -2.653   -1.869    1.869)    3.745
   3.094   (  -4.767    0.317   -0.317)    4.788
   3.157   (  -1.922   -1.355    1.355)    2.714
   3.203   (  -3.968    1.102   -1.102)    4.263
   3.228   (  -4.994    0.907   -0.907)    5.156
   3.382   (  -1.266   -1.603    1.603)    2.596
   3.410   (  -1.990    0.506   -0.506)    2.114
   3.637   (   6.176    0.816   -0.816)    6.282
   3.656   (   6.417   -0.197    0.197)    6.423
   3.728   (   5.314   -1.501    1.501)    5.722
   3.786   (  -0.987    3.062   -3.062)    4.441
   3.841   (  -1.668    3.780   -3.780)    5.601
   3.842   (  -3.334    0.968   -0.968)    3.605
   3.951   (   3.427   -2.377    2.377)    4.801
   3.956   (   2.193   -0.462    0.462)    2.288
   3.987   (  -0.408    0.352   -0.352)    0.644
   4.032   (  -1.414   -0.866    0.866)    1.870
   4.071   (  -3.729   -1.746    1.746)    4.473
   4.074   (  -2.265   -0.078    0.078)    2.268
   4.145   (  -2.350   -2.034    2.034)    3.714
   4.149   (  -2.549   -1.748    1.748)    3.551
   4.212   (  -1.761    0.124   -0.124)    1.770
   4.404   (   1.079    1.910   -1.910)    2.909
   4.410   (   0.278    1.599   -1.599)    2.278
   4.459   (   0.385    0.702   -0.702)    1.065
   4.469   (   2.240   -0.117    0.117)    2.246
   4.503   (  -0.021    0.439   -0.439)    0.621
   4.509   (  -0.729    0.340   -0.340)    0.873
   4.601   (   0.770   -0.478    0.478)    1.025
   4.618   (  -0.254    0.123   -0.123)    0.308
   4.636   (   1.041   -1.431    1.431)    2.275
   4.661   (  -0.750   -0.571    0.571)    1.102
   4.686   (  -1.461    0.370   -0.370)    1.552
   4.897   (   2.157    1.239   -1.239)    2.779
   5.187   (  -0.303   -0.705    0.705)    1.042
   5.198   (  -0.005   -1.159    1.159)    1.640
   5.232   (  -2.462    0.077   -0.077)    2.464
   5.319   (   1.105    0.553   -0.553)    1.353
======================= Grid point 110 (14/16) =======================
q-point: (-0.33  0.33  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 30
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.533   (  -1.531   -0.000   -0.000)    1.531
   0.533   (  -1.531   -0.000   -0.000)    1.531
   0.625   (  -6.317   -0.000   -0.000)    6.317
   0.854   (   3.499    0.000    0.000)    3.499
   0.854   (   3.499    0.000    0.000)    3.499
   0.987   (  -2.198   -0.000   -0.000)    2.198
   1.019   (   8.404    0.000    0.000)    8.404
   1.019   (   8.404    0.000    0.000)    8.404
   1.115   (   3.584    0.000    0.000)    3.584
   1.613   (  13.008    0.000    0.000)   13.008
   1.726   (  10.128    0.000    0.000)   10.128
   1.770   (   1.333    0.000    0.000)    1.333
   1.770   (   1.333    0.000    0.000)    1.333
   1.847   (   7.823    0.000    0.000)    7.823
   1.894   (   0.891    0.000    0.000)    0.891
   2.261   (  -3.302   -0.000   -0.000)    3.302
   2.263   (   0.645    0.000    0.000)    0.645
   2.263   (   0.645    0.000    0.000)    0.645
   2.289   (   0.396    0.000    0.000)    0.396
   2.389   (   0.363    0.000    0.000)    0.363
   2.404   (  -6.525   -0.000   -0.000)    6.525
   2.404   (  -6.525   -0.000   -0.000)    6.525
   2.460   (  -1.304   -0.000   -0.000)    1.304
   2.696   (   1.664    0.000    0.000)    1.664
   2.696   (   1.664    0.000    0.000)    1.664
   2.777   (   2.708    0.000    0.000)    2.708
   2.777   (   2.708    0.000    0.000)    2.708
   2.869   (  -0.318   -0.000   -0.000)    0.318
   2.897   (  -1.010   -0.000   -0.000)    1.010
   2.939   (  -1.171   -0.000   -0.000)    1.171
   2.986   (  -3.270   -0.000   -0.000)    3.270
   2.986   (  -3.270   -0.000   -0.000)    3.270
   3.109   (  -5.372   -0.000   -0.000)    5.372
   3.135   (  -2.911   -0.000   -0.000)    2.911
   3.241   (  -6.463   -0.000   -0.000)    6.463
   3.241   (  -6.463   -0.000   -0.000)    6.463
   3.312   (  -5.145   -0.000   -0.000)    5.145
   3.410   (  -5.997   -0.000   -0.000)    5.997
   3.652   (   8.376    0.000    0.000)    8.376
   3.652   (   8.376    0.000    0.000)    8.376
   3.708   (   5.720    0.000    0.000)    5.720
   3.857   (  -1.520   -0.000   -0.000)    1.520
   3.857   (  -1.520   -0.000   -0.000)    1.520
   3.874   (   3.301    0.000    0.000)    3.301
   3.934   (   3.325    0.000    0.000)    3.325
   3.939   (  -2.200   -0.000   -0.000)    2.200
   4.006   (   0.111    0.000    0.000)    0.111
   4.006   (   0.111    0.000    0.000)    0.111
   4.066   (  -2.944   -0.000   -0.000)    2.944
   4.090   (  -2.952   -0.000   -0.000)    2.952
   4.090   (  -2.952   -0.000   -0.000)    2.952
   4.112   (  -2.564   -0.000   -0.000)    2.564
   4.221   (  -1.187   -0.000   -0.000)    1.187
   4.448   (  -0.813   -0.000   -0.000)    0.813
   4.448   (  -0.813   -0.000   -0.000)    0.813
   4.457   (   1.583    0.000    0.000)    1.583
   4.462   (   0.345    0.000    0.000)    0.345
   4.506   (  -0.499   -0.000   -0.000)    0.499
   4.506   (  -0.499   -0.000   -0.000)    0.499
   4.588   (   3.045    0.000    0.000)    3.045
   4.621   (  -0.075   -0.000   -0.000)    0.075
   4.621   (  -0.075   -0.000   -0.000)    0.075
   4.639   (  -1.463   -0.000   -0.000)    1.463
   4.696   (  -1.511   -0.000   -0.000)    1.511
   4.917   (   2.902    0.000    0.000)    2.902
   5.177   (  -0.061   -0.000   -0.000)    0.061
   5.177   (  -0.061   -0.000   -0.000)    0.061
   5.233   (  -2.787   -0.000   -0.000)    2.787
   5.328   (   1.226    0.000    0.000)    1.226
======================= Grid point 115 (15/16) =======================
q-point: (-0.50  0.50  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 39
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.517   (  -0.533   -0.533   -0.000)    0.753
   0.520   (  -0.834   -0.834   -0.000)    1.179
   0.572   (  -3.724   -3.724   -0.000)    5.266
   0.886   (   2.554    2.554   -0.000)    3.612
   0.910   (   1.343    1.343   -0.000)    1.899
   0.946   (  -1.507   -1.507    0.000)    2.131
   1.051   (   6.541    6.541   -0.000)    9.251
   1.095   (   3.026    3.026   -0.000)    4.280
   1.170   (  -1.494   -1.494   -0.000)    2.112
   1.696   (   8.400    8.400   -0.000)   11.879
   1.752   (   5.505    5.505   -0.000)    7.785
   1.815   (   5.497    5.497    0.000)    7.775
   1.849   (   2.075    2.075   -0.000)    2.934
   1.906   (  -1.349   -1.349   -0.000)    1.908
   1.909   (   4.977    4.977   -0.000)    7.039
   2.231   (  -2.032   -2.032   -0.000)    2.874
   2.262   (  -0.876   -0.876   -0.000)    1.239
   2.273   (  -1.441   -1.441   -0.000)    2.037
   2.295   (   0.046    0.046   -0.000)    0.065
   2.329   (  -3.541   -3.541   -0.000)    5.008
   2.333   (  -4.463   -4.463    0.000)    6.311
   2.425   (   0.306    0.306    0.000)    0.432
   2.437   (   1.051    1.051   -0.000)    1.486
   2.720   (   0.874    0.874   -0.000)    1.236
   2.775   (   1.904    1.904   -0.000)    2.693
   2.778   (  -0.529   -0.529    0.000)    0.748
   2.795   (   4.137    4.137   -0.000)    5.850
   2.850   (  -1.830   -1.830   -0.000)    2.588
   2.865   (   0.504    0.504   -0.000)    0.713
   2.934   (  -0.972   -0.972   -0.000)    1.374
   2.953   (  -2.180   -2.180   -0.000)    3.083
   2.981   (  -3.377   -3.377   -0.000)    4.775
   3.049   (  -0.723   -0.723   -0.000)    1.022
   3.123   (  -1.426   -1.426    0.000)    2.017
   3.133   (  -2.559   -2.559   -0.000)    3.619
   3.181   (  -5.637   -5.637   -0.000)    7.972
   3.302   (  -3.378   -3.378   -0.000)    4.778
   3.341   (  -4.747   -4.747   -0.000)    6.713
   3.725   (   4.822    4.822   -0.000)    6.820
   3.766   (   2.914    2.914    0.000)    4.122
   3.778   (   2.467    2.467    0.000)    3.489
   3.792   (   2.433    2.433   -0.000)    3.440
   3.833   (  -0.321   -0.321   -0.000)    0.454
   3.859   (  -1.288   -1.288   -0.000)    1.822
   3.969   (   0.494    0.494    0.000)    0.699
   3.975   (   0.802    0.802   -0.000)    1.134
   3.996   (   0.076    0.076   -0.000)    0.107
   4.016   (   0.600    0.600   -0.000)    0.848
   4.034   (   0.963    0.963   -0.000)    1.362
   4.045   (  -1.954   -1.954   -0.000)    2.764
   4.090   (  -2.772   -2.772   -0.000)    3.921
   4.092   (  -1.879   -1.879    0.000)    2.657
   4.218   (   1.242    1.242   -0.000)    1.756
   4.425   (   0.273    0.273   -0.000)    0.385
   4.427   (   0.258    0.258   -0.000)    0.365
   4.460   (  -1.360   -1.360   -0.000)    1.923
   4.484   (   0.308    0.308    0.000)    0.436
   4.490   (   0.051    0.051   -0.000)    0.073
   4.506   (  -0.949   -0.949   -0.000)    1.343
   4.607   (   0.586    0.586   -0.000)    0.828
   4.619   (  -0.109   -0.109   -0.000)    0.154
   4.628   (  -0.805   -0.805   -0.000)    1.138
   4.647   (   0.304    0.304   -0.000)    0.430
   4.670   (  -0.788   -0.788    0.000)    1.114
   4.941   (   1.787    1.787   -0.000)    2.528
   5.177   (   0.000    0.000   -0.000)    0.001
   5.189   (  -0.717   -0.717    0.000)    1.014
   5.199   (  -1.207   -1.207   -0.000)    1.707
   5.338   (   0.725    0.725   -0.000)    1.025
======================= Grid point 163 (16/16) =======================
q-point: (-0.50 -0.50  0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 8.55e-05 
Number of triplets: 10
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.506   (   0.000   -0.000    0.000)    0.000
   0.506   (   0.000   -0.000    0.000)    0.000
   0.506   (   0.000   -0.000    0.000)    0.000
   0.932   (   0.000   -0.000    0.000)    0.000
   0.932   (   0.000   -0.000    0.000)    0.000
   0.932   (   0.000   -0.000    0.000)    0.000
   1.145   (   0.000   -0.000    0.000)    0.000
   1.145   (   0.000   -0.000    0.000)    0.000
   1.145   (   0.000   -0.000    0.000)    0.000
   1.878   (   0.000   -0.000    0.000)    0.000
   1.878   (   0.000   -0.000    0.000)    0.000
   1.878   (   0.000   -0.000    0.000)    0.000
   1.897   (   0.000   -0.000    0.000)    0.000
   1.897   (   0.000   -0.000    0.000)    0.000
   1.997   (   0.000   -0.000    0.000)    0.000
   2.176   (   0.000   -0.000    0.000)    0.000
   2.176   (   0.000   -0.000    0.000)    0.000
   2.176   (   0.000   -0.000    0.000)    0.000
   2.296   (   0.000   -0.000    0.000)    0.000
   2.296   (   0.000   -0.000    0.000)    0.000
   2.296   (   0.000   -0.000    0.000)    0.000
   2.471   (   0.000   -0.000    0.000)    0.000
   2.471   (   0.000   -0.000    0.000)    0.000
   2.751   (   0.000   -0.000    0.000)    0.000
   2.751   (   0.000   -0.000    0.000)    0.000
   2.751   (   0.000   -0.000    0.000)    0.000
   2.875   (  -0.000    0.000   -0.000)    0.000
   2.875   (  -0.000    0.000   -0.000)    0.000
   2.875   (   0.000   -0.000    0.000)    0.000
   2.944   (  -0.000    0.000   -0.000)    0.000
   2.944   (  -0.000    0.000   -0.000)    0.000
   2.944   (   0.000   -0.000    0.000)    0.000
   3.041   (   0.000   -0.000    0.000)    0.000
   3.041   (   0.000   -0.000    0.000)    0.000
   3.041   (   0.000   -0.000    0.000)    0.000
   3.077   (   0.000   -0.000    0.000)    0.000
   3.250   (  -0.000    0.000   -0.000)    0.000
   3.250   (   0.000   -0.000    0.000)    0.000
   3.821   (   0.000   -0.000    0.000)    0.000
   3.821   (   0.000   -0.000    0.000)    0.000
   3.821   (   0.000   -0.000    0.000)    0.000
   3.828   (   0.000   -0.000    0.000)    0.000
   3.828   (   0.000   -0.000    0.000)    0.000
   3.828   (   0.000   -0.000    0.000)    0.000
   3.996   (  -0.000    0.000   -0.000)    0.000
   3.996   (   0.000   -0.000    0.000)    0.000
   3.999   (   0.000   -0.000    0.000)    0.000
   3.999   (   0.000   -0.000    0.000)    0.000
   3.999   (   0.000   -0.000    0.000)    0.000
   4.058   (   0.000   -0.000    0.000)    0.000
   4.058   (   0.000   -0.000    0.000)    0.000
   4.058   (   0.000   -0.000    0.000)    0.000
   4.255   (   0.000   -0.000    0.000)    0.000
   4.420   (   0.000   -0.000    0.000)    0.000
   4.420   (   0.000   -0.000    0.000)    0.000
   4.420   (   0.000   -0.000    0.000)    0.000
   4.490   (   0.000   -0.000    0.000)    0.000
   4.490   (   0.000   -0.000    0.000)    0.000
   4.490   (   0.000   -0.000    0.000)    0.000
   4.617   (   0.000   -0.000    0.000)    0.000
   4.617   (   0.000   -0.000    0.000)    0.000
   4.617   (   0.000   -0.000    0.000)    0.000
   4.656   (  -0.000    0.000   -0.000)    0.000
   4.656   (   0.000   -0.000    0.000)    0.000
   4.970   (   0.000   -0.000    0.000)    0.000
   5.178   (   0.000   -0.000    0.000)    0.000
   5.178   (   0.000   -0.000    0.000)    0.000
   5.178   (   0.000   -0.000    0.000)    0.000
   5.350   (   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/14904
   10.0      9.314      9.314      9.314     -0.000      0.000      0.000 3/14904
   20.0      4.357      4.357      4.357     -0.000      0.000      0.000 3/14904
   30.0      2.542      2.542      2.542     -0.000      0.000      0.000 3/14904
   40.0      1.692      1.692      1.692     -0.000      0.000      0.000 3/14904
   50.0      1.243      1.243      1.243     -0.000      0.000      0.000 3/14904
   60.0      0.978      0.978      0.978     -0.000      0.000      0.000 3/14904
   70.0      0.807      0.807      0.807     -0.000      0.000      0.000 3/14904
   80.0      0.688      0.688      0.688     -0.000      0.000      0.000 3/14904
   90.0      0.600      0.600      0.600     -0.000      0.000      0.000 3/14904
  100.0      0.532      0.532      0.532     -0.000      0.000      0.000 3/14904
  110.0      0.479      0.479      0.479     -0.000      0.000      0.000 3/14904
  120.0      0.436      0.436      0.436     -0.000      0.000      0.000 3/14904
  130.0      0.400      0.400      0.400     -0.000      0.000      0.000 3/14904
  140.0      0.369      0.369      0.369     -0.000      0.000      0.000 3/14904
  150.0      0.343      0.343      0.343     -0.000      0.000      0.000 3/14904
  160.0      0.321      0.321      0.321     -0.000      0.000      0.000 3/14904
  170.0      0.301      0.301      0.301     -0.000      0.000      0.000 3/14904
  180.0      0.284      0.284      0.284     -0.000      0.000      0.000 3/14904
  190.0      0.269      0.269      0.269     -0.000      0.000      0.000 3/14904
  200.0      0.255      0.255      0.255     -0.000      0.000      0.000 3/14904
  210.0      0.242      0.242      0.242     -0.000      0.000      0.000 3/14904
  220.0      0.231      0.231      0.231     -0.000      0.000      0.000 3/14904
  230.0      0.221      0.221      0.221     -0.000      0.000      0.000 3/14904
  240.0      0.211      0.211      0.211     -0.000      0.000      0.000 3/14904
  250.0      0.203      0.203      0.203     -0.000      0.000      0.000 3/14904
  260.0      0.195      0.195      0.195     -0.000      0.000      0.000 3/14904
  270.0      0.188      0.188      0.188     -0.000      0.000      0.000 3/14904
  280.0      0.181      0.181      0.181     -0.000      0.000      0.000 3/14904
  290.0      0.174      0.174      0.174     -0.000      0.000      0.000 3/14904
  300.0      0.169      0.169      0.169     -0.000      0.000      0.000 3/14904
  310.0      0.163      0.163      0.163     -0.000      0.000      0.000 3/14904
  320.0      0.158      0.158      0.158     -0.000      0.000      0.000 3/14904
  330.0      0.153      0.153      0.153     -0.000      0.000      0.000 3/14904
  340.0      0.149      0.149      0.149     -0.000      0.000      0.000 3/14904
  350.0      0.144      0.144      0.144     -0.000      0.000      0.000 3/14904
  360.0      0.140      0.140      0.140     -0.000      0.000      0.000 3/14904
  370.0      0.136      0.136      0.136     -0.000      0.000      0.000 3/14904
  380.0      0.133      0.133      0.133     -0.000      0.000      0.000 3/14904
  390.0      0.129      0.129      0.129     -0.000      0.000      0.000 3/14904
  400.0      0.126      0.126      0.126     -0.000      0.000      0.000 3/14904
  410.0      0.123      0.123      0.123     -0.000      0.000      0.000 3/14904
  420.0      0.120      0.120      0.120     -0.000      0.000      0.000 3/14904
  430.0      0.117      0.117      0.117     -0.000      0.000      0.000 3/14904
  440.0      0.115      0.115      0.115     -0.000      0.000      0.000 3/14904
  450.0      0.112      0.112      0.112     -0.000      0.000      0.000 3/14904
  460.0      0.110      0.110      0.110     -0.000      0.000      0.000 3/14904
  470.0      0.107      0.107      0.107     -0.000      0.000      0.000 3/14904
  480.0      0.105      0.105      0.105     -0.000      0.000      0.000 3/14904
  490.0      0.103      0.103      0.103     -0.000      0.000      0.000 3/14904
  500.0      0.101      0.101      0.101     -0.000      0.000      0.000 3/14904
  510.0      0.099      0.099      0.099     -0.000      0.000      0.000 3/14904
  520.0      0.097      0.097      0.097     -0.000      0.000      0.000 3/14904
  530.0      0.095      0.095      0.095     -0.000      0.000      0.000 3/14904
  540.0      0.093      0.093      0.093     -0.000      0.000      0.000 3/14904
  550.0      0.092      0.092      0.092     -0.000      0.000      0.000 3/14904
  560.0      0.090      0.090      0.090     -0.000      0.000      0.000 3/14904
  570.0      0.088      0.088      0.088     -0.000      0.000      0.000 3/14904
  580.0      0.087      0.087      0.087     -0.000      0.000      0.000 3/14904
  590.0      0.085      0.085      0.085     -0.000      0.000      0.000 3/14904
  600.0      0.084      0.084      0.084     -0.000      0.000      0.000 3/14904
  610.0      0.083      0.083      0.083     -0.000      0.000      0.000 3/14904
  620.0      0.081      0.081      0.081     -0.000      0.000      0.000 3/14904
  630.0      0.080      0.080      0.080     -0.000      0.000      0.000 3/14904
  640.0      0.079      0.079      0.079     -0.000      0.000      0.000 3/14904
  650.0      0.078      0.078      0.078     -0.000      0.000      0.000 3/14904
  660.0      0.076      0.076      0.076     -0.000      0.000      0.000 3/14904
  670.0      0.075      0.075      0.075     -0.000      0.000      0.000 3/14904
  680.0      0.074      0.074      0.074     -0.000      0.000      0.000 3/14904
  690.0      0.073      0.073      0.073     -0.000      0.000      0.000 3/14904
  700.0      0.072      0.072      0.072     -0.000      0.000      0.000 3/14904
  710.0      0.071      0.071      0.071     -0.000      0.000      0.000 3/14904
  720.0      0.070      0.070      0.070     -0.000      0.000      0.000 3/14904
  730.0      0.069      0.069      0.069     -0.000      0.000      0.000 3/14904
  740.0      0.068      0.068      0.068     -0.000      0.000      0.000 3/14904
  750.0      0.067      0.067      0.067     -0.000      0.000      0.000 3/14904
  760.0      0.066      0.066      0.066     -0.000      0.000      0.000 3/14904
  770.0      0.065      0.065      0.065     -0.000      0.000      0.000 3/14904
  780.0      0.065      0.065      0.065     -0.000      0.000      0.000 3/14904
  790.0      0.064      0.064      0.064     -0.000      0.000      0.000 3/14904
  800.0      0.063      0.063      0.063     -0.000      0.000      0.000 3/14904
  810.0      0.062      0.062      0.062     -0.000      0.000      0.000 3/14904
  820.0      0.061      0.061      0.061     -0.000      0.000      0.000 3/14904
  830.0      0.061      0.061      0.061     -0.000      0.000      0.000 3/14904
  840.0      0.060      0.060      0.060     -0.000      0.000      0.000 3/14904
  850.0      0.059      0.059      0.059     -0.000      0.000      0.000 3/14904
  860.0      0.059      0.059      0.059     -0.000      0.000      0.000 3/14904
  870.0      0.058      0.058      0.058     -0.000      0.000      0.000 3/14904
  880.0      0.057      0.057      0.057     -0.000      0.000      0.000 3/14904
  890.0      0.057      0.057      0.057     -0.000      0.000      0.000 3/14904
  900.0      0.056      0.056      0.056     -0.000      0.000      0.000 3/14904
  910.0      0.055      0.055      0.055     -0.000      0.000      0.000 3/14904
  920.0      0.055      0.055      0.055     -0.000      0.000      0.000 3/14904
  930.0      0.054      0.054      0.054     -0.000      0.000      0.000 3/14904
  940.0      0.054      0.054      0.054     -0.000      0.000      0.000 3/14904
  950.0      0.053      0.053      0.053     -0.000      0.000      0.000 3/14904
  960.0      0.052      0.052      0.052     -0.000      0.000      0.000 3/14904
  970.0      0.052      0.052      0.052     -0.000      0.000      0.000 3/14904
  980.0      0.051      0.051      0.051     -0.000      0.000      0.000 3/14904
  990.0      0.051      0.051      0.051     -0.000      0.000      0.000 3/14904
 1000.0      0.050      0.050      0.050     -0.000      0.000      0.000 3/14904

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

