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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [2 2 1]
  Primitive matrix:
    [ 1. -0.  0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P-31c (163)
Number of symmetry operations in supercell: 48
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    6.893806931178462    0.000000000000000    0.000000000000000
  b   -3.446903465589231    5.970211931185789    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
   *1 Na  0.33333333333333  0.66666666666667  0.04001250876844  22.990
    2 Na  0.33333333333333  0.66666666666667  0.45998749123156  22.990
   *3 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990
    4 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990
    5 Na  0.66666666666667  0.33333333333333  0.95998749123156  22.990
    6 Na  0.66666666666667  0.33333333333333  0.54001250876844  22.990
   *7 Cl  0.36395380808473  0.04835563217144  0.12803350975519  35.453
    8 Cl  0.95164436782856  0.31559817591329  0.12803350975519  35.453
    9 Cl  0.68440182408671  0.63604619191527  0.12803350975519  35.453
   10 Cl  0.63604619191527  0.68440182408671  0.62803350975519  35.453
   11 Cl  0.31559817591329  0.95164436782856  0.62803350975519  35.453
   12 Cl  0.63604619191527  0.95164436782856  0.87196649024481  35.453
   13 Cl  0.31559817591329  0.36395380808473  0.87196649024481  35.453
   14 Cl  0.04835563217144  0.68440182408671  0.87196649024481  35.453
   15 Cl  0.36395380808473  0.31559817591329  0.37196649024481  35.453
   16 Cl  0.95164436782856  0.63604619191527  0.37196649024481  35.453
   17 Cl  0.04835563217144  0.36395380808473  0.62803350975519  35.453
   18 Cl  0.68440182408671  0.04835563217144  0.37196649024481  35.453
  *19 In  0.33333333333333  0.66666666666667  0.75000000000000 114.818
   20 In  0.66666666666667  0.33333333333333  0.25000000000000 114.818
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    6.893806931178462    0.000000000000000    0.000000000000000
  b   -3.446903465589231    5.970211931185789    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
   *1 Na  0.33333333333333  0.66666666666667  0.04001250876844  22.990 > 1
    2 Na  0.33333333333333  0.66666666666667  0.45998749123156  22.990 > 2
   *3 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990 > 3
    4 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990 > 4
    5 Na  0.66666666666667  0.33333333333333  0.95998749123156  22.990 > 5
    6 Na  0.66666666666667  0.33333333333333  0.54001250876844  22.990 > 6
   *7 Cl  0.36395380808473  0.04835563217144  0.12803350975519  35.453 > 7
    8 Cl  0.95164436782856  0.31559817591329  0.12803350975519  35.453 > 8
    9 Cl  0.68440182408671  0.63604619191527  0.12803350975519  35.453 > 9
   10 Cl  0.63604619191527  0.68440182408671  0.62803350975519  35.453 > 10
   11 Cl  0.31559817591329  0.95164436782856  0.62803350975519  35.453 > 11
   12 Cl  0.63604619191527  0.95164436782856  0.87196649024481  35.453 > 12
   13 Cl  0.31559817591329  0.36395380808473  0.87196649024481  35.453 > 13
   14 Cl  0.04835563217144  0.68440182408671  0.87196649024481  35.453 > 14
   15 Cl  0.36395380808473  0.31559817591329  0.37196649024481  35.453 > 15
   16 Cl  0.95164436782856  0.63604619191527  0.37196649024481  35.453 > 16
   17 Cl  0.04835563217144  0.36395380808473  0.62803350975519  35.453 > 17
   18 Cl  0.68440182408671  0.04835563217144  0.37196649024481  35.453 > 18
  *19 In  0.33333333333333  0.66666666666667  0.75000000000000 114.818 > 19
   20 In  0.66666666666667  0.33333333333333  0.25000000000000 114.818 > 20
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   13.787613862356924    0.000000000000000    0.000000000000000
  b   -6.893806931178462   11.940423862371578    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
   *1 Na  0.16666666666667  0.33333333333333  0.04001250876844  22.990 > 1
    2 Na  0.66666666666667  0.33333333333333  0.04001250876844  22.990 > 1
    3 Na  0.16666666666667  0.83333333333333  0.04001250876844  22.990 > 1
    4 Na  0.66666666666667  0.83333333333333  0.04001250876844  22.990 > 1
    5 Na  0.16666666666667  0.33333333333333  0.45998749123156  22.990 > 2
    6 Na  0.66666666666667  0.33333333333333  0.45998749123156  22.990 > 2
    7 Na  0.16666666666667  0.83333333333333  0.45998749123156  22.990 > 2
    8 Na  0.66666666666667  0.83333333333333  0.45998749123156  22.990 > 2
   *9 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990 > 3
   10 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 3
   11 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 3
   12 Na  0.50000000000000  0.50000000000000  0.75000000000000  22.990 > 3
   13 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990 > 4
   14 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 4
   15 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 4
   16 Na  0.50000000000000  0.50000000000000  0.25000000000000  22.990 > 4
   17 Na  0.33333333333333  0.16666666666667  0.95998749123156  22.990 > 5
   18 Na  0.83333333333333  0.16666666666667  0.95998749123156  22.990 > 5
   19 Na  0.33333333333333  0.66666666666667  0.95998749123156  22.990 > 5
   20 Na  0.83333333333333  0.66666666666667  0.95998749123156  22.990 > 5
   21 Na  0.33333333333333  0.16666666666667  0.54001250876844  22.990 > 6
   22 Na  0.83333333333333  0.16666666666667  0.54001250876844  22.990 > 6
   23 Na  0.33333333333333  0.66666666666667  0.54001250876844  22.990 > 6
   24 Na  0.83333333333333  0.66666666666667  0.54001250876844  22.990 > 6
  *25 Cl  0.18197690404237  0.02417781608572  0.12803350975519  35.453 > 7
   26 Cl  0.68197690404237  0.02417781608572  0.12803350975519  35.453 > 7
   27 Cl  0.18197690404237  0.52417781608572  0.12803350975519  35.453 > 7
   28 Cl  0.68197690404237  0.52417781608572  0.12803350975519  35.453 > 7
   29 Cl  0.47582218391428  0.15779908795665  0.12803350975519  35.453 > 8
   30 Cl  0.97582218391428  0.15779908795665  0.12803350975519  35.453 > 8
   31 Cl  0.47582218391428  0.65779908795665  0.12803350975519  35.453 > 8
   32 Cl  0.97582218391428  0.65779908795665  0.12803350975519  35.453 > 8
   33 Cl  0.34220091204335  0.31802309595763  0.12803350975519  35.453 > 9
   34 Cl  0.84220091204335  0.31802309595763  0.12803350975519  35.453 > 9
   35 Cl  0.34220091204335  0.81802309595763  0.12803350975519  35.453 > 9
   36 Cl  0.84220091204335  0.81802309595763  0.12803350975519  35.453 > 9
   37 Cl  0.31802309595763  0.34220091204335  0.62803350975519  35.453 > 10
   38 Cl  0.81802309595763  0.34220091204335  0.62803350975519  35.453 > 10
   39 Cl  0.31802309595763  0.84220091204335  0.62803350975519  35.453 > 10
   40 Cl  0.81802309595763  0.84220091204335  0.62803350975519  35.453 > 10
   41 Cl  0.15779908795665  0.47582218391428  0.62803350975519  35.453 > 11
   42 Cl  0.65779908795665  0.47582218391428  0.62803350975519  35.453 > 11
   43 Cl  0.15779908795665  0.97582218391428  0.62803350975519  35.453 > 11
   44 Cl  0.65779908795665  0.97582218391428  0.62803350975519  35.453 > 11
   45 Cl  0.31802309595763  0.47582218391428  0.87196649024481  35.453 > 12
   46 Cl  0.81802309595763  0.47582218391428  0.87196649024481  35.453 > 12
   47 Cl  0.31802309595763  0.97582218391428  0.87196649024481  35.453 > 12
   48 Cl  0.81802309595763  0.97582218391428  0.87196649024481  35.453 > 12
   49 Cl  0.15779908795665  0.18197690404237  0.87196649024481  35.453 > 13
   50 Cl  0.65779908795665  0.18197690404237  0.87196649024481  35.453 > 13
   51 Cl  0.15779908795665  0.68197690404237  0.87196649024481  35.453 > 13
   52 Cl  0.65779908795665  0.68197690404237  0.87196649024481  35.453 > 13
   53 Cl  0.02417781608572  0.34220091204335  0.87196649024481  35.453 > 14
   54 Cl  0.52417781608572  0.34220091204335  0.87196649024481  35.453 > 14
   55 Cl  0.02417781608572  0.84220091204335  0.87196649024481  35.453 > 14
   56 Cl  0.52417781608572  0.84220091204335  0.87196649024481  35.453 > 14
   57 Cl  0.18197690404237  0.15779908795665  0.37196649024481  35.453 > 15
   58 Cl  0.68197690404237  0.15779908795665  0.37196649024481  35.453 > 15
   59 Cl  0.18197690404237  0.65779908795665  0.37196649024481  35.453 > 15
   60 Cl  0.68197690404237  0.65779908795665  0.37196649024481  35.453 > 15
   61 Cl  0.47582218391428  0.31802309595763  0.37196649024481  35.453 > 16
   62 Cl  0.97582218391428  0.31802309595763  0.37196649024481  35.453 > 16
   63 Cl  0.47582218391428  0.81802309595763  0.37196649024481  35.453 > 16
   64 Cl  0.97582218391428  0.81802309595763  0.37196649024481  35.453 > 16
   65 Cl  0.02417781608572  0.18197690404237  0.62803350975519  35.453 > 17
   66 Cl  0.52417781608572  0.18197690404237  0.62803350975519  35.453 > 17
   67 Cl  0.02417781608572  0.68197690404237  0.62803350975519  35.453 > 17
   68 Cl  0.52417781608572  0.68197690404237  0.62803350975519  35.453 > 17
   69 Cl  0.34220091204335  0.02417781608572  0.37196649024481  35.453 > 18
   70 Cl  0.84220091204335  0.02417781608572  0.37196649024481  35.453 > 18
   71 Cl  0.34220091204335  0.52417781608572  0.37196649024481  35.453 > 18
   72 Cl  0.84220091204335  0.52417781608572  0.37196649024481  35.453 > 18
  *73 In  0.16666666666667  0.33333333333333  0.75000000000000 114.818 > 19
   74 In  0.66666666666667  0.33333333333333  0.75000000000000 114.818 > 19
   75 In  0.16666666666667  0.83333333333333  0.75000000000000 114.818 > 19
   76 In  0.66666666666667  0.83333333333333  0.75000000000000 114.818 > 19
   77 In  0.33333333333333  0.16666666666667  0.25000000000000 114.818 > 20
   78 In  0.83333333333333  0.16666666666667  0.25000000000000 114.818 > 20
   79 In  0.33333333333333  0.66666666666667  0.25000000000000 114.818 > 20
   80 In  0.83333333333333  0.66666666666667  0.25000000000000 114.818 > 20
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.7129661   -0.0000000    0.0000000
            0.0000000    2.7129661    0.0000000
            0.0000000    0.0000000    2.7894726
-------------------------- Born effective charges --------------------------
    1 Na    1.0949434    0.0922169   -0.0000000
           -0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    2 Na    1.0949434   -0.0922169   -0.0000000
            0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    3 Na    1.2201615    0.0000000   -0.0000000
           -0.0000000    1.2201615   -0.0000000
           -0.0000000   -0.0000000    0.8640476
    4 Na    1.2201615    0.0000000   -0.0000000
           -0.0000000    1.2201615   -0.0000000
           -0.0000000   -0.0000000    0.8640476
    5 Na    1.0949434    0.0922169   -0.0000000
           -0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    6 Na    1.0949434   -0.0922169   -0.0000000
            0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    7 Cl   -0.9181712   -0.2124562   -0.0588922
           -0.1928141   -1.1792767   -0.2944857
           -0.0341992   -0.3121885   -1.0695868
    8 Cl   -1.2894875   -0.0215655    0.2844782
           -0.0019234   -0.8079604    0.0962407
            0.2874628    0.1264769   -1.0695868
    9 Cl   -0.9385132    0.2045585   -0.2255860
            0.2242007   -1.1589348    0.1982450
           -0.2532636    0.1857116   -1.0695868
   10 Cl   -1.2894875    0.0215655    0.2844782
            0.0019234   -0.8079604   -0.0962407
            0.2874628   -0.1264769   -1.0695868
   11 Cl   -0.9181712    0.2124562   -0.0588922
            0.1928141   -1.1792767    0.2944857
           -0.0341992    0.3121885   -1.0695868
   12 Cl   -0.9181712   -0.2124562   -0.0588922
           -0.1928141   -1.1792767   -0.2944857
           -0.0341992   -0.3121885   -1.0695868
   13 Cl   -0.9385132    0.2045585   -0.2255860
            0.2242007   -1.1589348    0.1982450
           -0.2532636    0.1857116   -1.0695868
   14 Cl   -1.2894875   -0.0215655    0.2844782
           -0.0019234   -0.8079604    0.0962407
            0.2874628    0.1264769   -1.0695868
   15 Cl   -1.2894875    0.0215655    0.2844782
            0.0019234   -0.8079604   -0.0962407
            0.2874628   -0.1264769   -1.0695868
   16 Cl   -0.9385132   -0.2045585   -0.2255860
           -0.2242007   -1.1589348   -0.1982450
           -0.2532636   -0.1857116   -1.0695868
   17 Cl   -0.9385132   -0.2045585   -0.2255860
           -0.2242007   -1.1589348   -0.1982450
           -0.2532636   -0.1857116   -1.0695868
   18 Cl   -0.9181712    0.2124562   -0.0588922
            0.1928141   -1.1792767    0.2944857
           -0.0341992    0.3121885   -1.0695868
   19 In    2.8822955   -0.0000000   -0.0000000
            0.0000000    2.8822955   -0.0000000
           -0.0000000   -0.0000000    3.0514539
   20 In    2.8822955   -0.0000000   -0.0000000
            0.0000000    2.8822955   -0.0000000
           -0.0000000   -0.0000000    3.0514539
----------------------------------------------------------------------------
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: 240/240
Permutation basis: 7440/7440
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 654
Number of blocks in projector: 384
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 304
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 256
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 94
Use standard eigh solver.
Tree of FC basis block matrices:
- (654, 638), data: False
|-- (94, 90), data: True
|-- (256, 252), data: True
|-- (304, 296), data: True
-----
Solver_atoms: 1 -- 80 / 80
Time (Solver_compr_matrix_reshape): 0.005
Solver_block: 100 / 240
 - Time: 0.414
Solver_block: 200 / 240
 - Time: 0.422
Solver_block: 240 / 240
 - Time: 0.195
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.046
--------------------------------- Symfc end --------------------------------
Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) 
Permutation basis: 240/240
Permutation basis: 7440/7440
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 654
Number of blocks in projector: 384
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 304
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 256
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 94
Use standard eigh solver.
Tree of FC basis block matrices:
- (654, 638), data: False
|-- (94, 90), data: True
|-- (256, 252), data: True
|-- (304, 296), 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 05:28:28]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 05:28:28]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [2 2 1]
Primitive matrix:
  [ 1. -0.  0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P-31c (163)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    6.893806931178462    0.000000000000000    0.000000000000000
  b   -3.446903465589231    5.970211931185789    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
    1 Na  0.33333333333333  0.66666666666667  0.04001250876844  22.990
    2 Na  0.33333333333333  0.66666666666667  0.45998749123156  22.990
    3 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990
    4 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990
    5 Na  0.66666666666667  0.33333333333333  0.95998749123156  22.990
    6 Na  0.66666666666667  0.33333333333333  0.54001250876844  22.990
    7 Cl  0.36395380808473  0.04835563217144  0.12803350975519  35.453
    8 Cl  0.95164436782856  0.31559817591329  0.12803350975519  35.453
    9 Cl  0.68440182408671  0.63604619191527  0.12803350975519  35.453
   10 Cl  0.63604619191527  0.68440182408671  0.62803350975519  35.453
   11 Cl  0.31559817591329  0.95164436782856  0.62803350975519  35.453
   12 Cl  0.63604619191527  0.95164436782856  0.87196649024481  35.453
   13 Cl  0.31559817591329  0.36395380808473  0.87196649024481  35.453
   14 Cl  0.04835563217144  0.68440182408671  0.87196649024481  35.453
   15 Cl  0.36395380808473  0.31559817591329  0.37196649024481  35.453
   16 Cl  0.95164436782856  0.63604619191527  0.37196649024481  35.453
   17 Cl  0.04835563217144  0.36395380808473  0.62803350975519  35.453
   18 Cl  0.68440182408671  0.04835563217144  0.37196649024481  35.453
   19 In  0.33333333333333  0.66666666666667  0.75000000000000 114.818
   20 In  0.66666666666667  0.33333333333333  0.25000000000000 114.818
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.787613862356924    0.000000000000000    0.000000000000000
  b   -6.893806931178462   11.940423862371578    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
    1 Na  0.16666666666667  0.33333333333333  0.04001250876844  22.990 > 1
    2 Na  0.66666666666667  0.33333333333333  0.04001250876844  22.990 > 1
    3 Na  0.16666666666667  0.83333333333333  0.04001250876844  22.990 > 1
    4 Na  0.66666666666667  0.83333333333333  0.04001250876844  22.990 > 1
    5 Na  0.16666666666667  0.33333333333333  0.45998749123156  22.990 > 5
    6 Na  0.66666666666667  0.33333333333333  0.45998749123156  22.990 > 5
    7 Na  0.16666666666667  0.83333333333333  0.45998749123156  22.990 > 5
    8 Na  0.66666666666667  0.83333333333333  0.45998749123156  22.990 > 5
    9 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990 > 9
   10 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 9
   11 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 9
   12 Na  0.50000000000000  0.50000000000000  0.75000000000000  22.990 > 9
   13 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990 > 13
   14 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 13
   15 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 13
   16 Na  0.50000000000000  0.50000000000000  0.25000000000000  22.990 > 13
   17 Na  0.33333333333333  0.16666666666667  0.95998749123156  22.990 > 17
   18 Na  0.83333333333333  0.16666666666667  0.95998749123156  22.990 > 17
   19 Na  0.33333333333333  0.66666666666667  0.95998749123156  22.990 > 17
   20 Na  0.83333333333333  0.66666666666667  0.95998749123156  22.990 > 17
   21 Na  0.33333333333333  0.16666666666667  0.54001250876844  22.990 > 21
   22 Na  0.83333333333333  0.16666666666667  0.54001250876844  22.990 > 21
   23 Na  0.33333333333333  0.66666666666667  0.54001250876844  22.990 > 21
   24 Na  0.83333333333333  0.66666666666667  0.54001250876844  22.990 > 21
   25 Cl  0.18197690404237  0.02417781608572  0.12803350975519  35.453 > 25
   26 Cl  0.68197690404237  0.02417781608572  0.12803350975519  35.453 > 25
   27 Cl  0.18197690404237  0.52417781608572  0.12803350975519  35.453 > 25
   28 Cl  0.68197690404237  0.52417781608572  0.12803350975519  35.453 > 25
   29 Cl  0.47582218391428  0.15779908795665  0.12803350975519  35.453 > 29
   30 Cl  0.97582218391428  0.15779908795665  0.12803350975519  35.453 > 29
   31 Cl  0.47582218391428  0.65779908795665  0.12803350975519  35.453 > 29
   32 Cl  0.97582218391428  0.65779908795665  0.12803350975519  35.453 > 29
   33 Cl  0.34220091204335  0.31802309595763  0.12803350975519  35.453 > 33
   34 Cl  0.84220091204335  0.31802309595763  0.12803350975519  35.453 > 33
   35 Cl  0.34220091204335  0.81802309595763  0.12803350975519  35.453 > 33
   36 Cl  0.84220091204335  0.81802309595763  0.12803350975519  35.453 > 33
   37 Cl  0.31802309595763  0.34220091204335  0.62803350975519  35.453 > 37
   38 Cl  0.81802309595763  0.34220091204335  0.62803350975519  35.453 > 37
   39 Cl  0.31802309595763  0.84220091204335  0.62803350975519  35.453 > 37
   40 Cl  0.81802309595763  0.84220091204335  0.62803350975519  35.453 > 37
   41 Cl  0.15779908795665  0.47582218391428  0.62803350975519  35.453 > 41
   42 Cl  0.65779908795665  0.47582218391428  0.62803350975519  35.453 > 41
   43 Cl  0.15779908795665  0.97582218391428  0.62803350975519  35.453 > 41
   44 Cl  0.65779908795665  0.97582218391428  0.62803350975519  35.453 > 41
   45 Cl  0.31802309595763  0.47582218391428  0.87196649024481  35.453 > 45
   46 Cl  0.81802309595763  0.47582218391428  0.87196649024481  35.453 > 45
   47 Cl  0.31802309595763  0.97582218391428  0.87196649024481  35.453 > 45
   48 Cl  0.81802309595763  0.97582218391428  0.87196649024481  35.453 > 45
   49 Cl  0.15779908795665  0.18197690404237  0.87196649024481  35.453 > 49
   50 Cl  0.65779908795665  0.18197690404237  0.87196649024481  35.453 > 49
   51 Cl  0.15779908795665  0.68197690404237  0.87196649024481  35.453 > 49
   52 Cl  0.65779908795665  0.68197690404237  0.87196649024481  35.453 > 49
   53 Cl  0.02417781608572  0.34220091204335  0.87196649024481  35.453 > 53
   54 Cl  0.52417781608572  0.34220091204335  0.87196649024481  35.453 > 53
   55 Cl  0.02417781608572  0.84220091204335  0.87196649024481  35.453 > 53
   56 Cl  0.52417781608572  0.84220091204335  0.87196649024481  35.453 > 53
   57 Cl  0.18197690404237  0.15779908795665  0.37196649024481  35.453 > 57
   58 Cl  0.68197690404237  0.15779908795665  0.37196649024481  35.453 > 57
   59 Cl  0.18197690404237  0.65779908795665  0.37196649024481  35.453 > 57
   60 Cl  0.68197690404237  0.65779908795665  0.37196649024481  35.453 > 57
   61 Cl  0.47582218391428  0.31802309595763  0.37196649024481  35.453 > 61
   62 Cl  0.97582218391428  0.31802309595763  0.37196649024481  35.453 > 61
   63 Cl  0.47582218391428  0.81802309595763  0.37196649024481  35.453 > 61
   64 Cl  0.97582218391428  0.81802309595763  0.37196649024481  35.453 > 61
   65 Cl  0.02417781608572  0.18197690404237  0.62803350975519  35.453 > 65
   66 Cl  0.52417781608572  0.18197690404237  0.62803350975519  35.453 > 65
   67 Cl  0.02417781608572  0.68197690404237  0.62803350975519  35.453 > 65
   68 Cl  0.52417781608572  0.68197690404237  0.62803350975519  35.453 > 65
   69 Cl  0.34220091204335  0.02417781608572  0.37196649024481  35.453 > 69
   70 Cl  0.84220091204335  0.02417781608572  0.37196649024481  35.453 > 69
   71 Cl  0.34220091204335  0.52417781608572  0.37196649024481  35.453 > 69
   72 Cl  0.84220091204335  0.52417781608572  0.37196649024481  35.453 > 69
   73 In  0.16666666666667  0.33333333333333  0.75000000000000 114.818 > 73
   74 In  0.66666666666667  0.33333333333333  0.75000000000000 114.818 > 73
   75 In  0.16666666666667  0.83333333333333  0.75000000000000 114.818 > 73
   76 In  0.66666666666667  0.83333333333333  0.75000000000000 114.818 > 73
   77 In  0.33333333333333  0.16666666666667  0.25000000000000 114.818 > 77
   78 In  0.83333333333333  0.16666666666667  0.25000000000000 114.818 > 77
   79 In  0.33333333333333  0.66666666666667  0.25000000000000 114.818 > 77
   80 In  0.83333333333333  0.66666666666667  0.25000000000000 114.818 > 77
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.7129661   -0.0000000    0.0000000
            0.0000000    2.7129661    0.0000000
            0.0000000    0.0000000    2.7894726
-------------------------- Born effective charges --------------------------
    1 Na    1.0949434    0.0922169   -0.0000000
           -0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    2 Na    1.0949434   -0.0922169   -0.0000000
            0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    3 Na    1.2201615    0.0000000   -0.0000000
           -0.0000000    1.2201615   -0.0000000
           -0.0000000   -0.0000000    0.8640476
    4 Na    1.2201615    0.0000000   -0.0000000
           -0.0000000    1.2201615   -0.0000000
           -0.0000000   -0.0000000    0.8640476
    5 Na    1.0949434    0.0922169   -0.0000000
           -0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    6 Na    1.0949434   -0.0922169   -0.0000000
            0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    7 Cl   -0.9181712   -0.2124562   -0.0588922
           -0.1928141   -1.1792767   -0.2944857
           -0.0341992   -0.3121885   -1.0695868
    8 Cl   -1.2894875   -0.0215655    0.2844782
           -0.0019234   -0.8079604    0.0962407
            0.2874628    0.1264769   -1.0695868
    9 Cl   -0.9385132    0.2045585   -0.2255860
            0.2242007   -1.1589348    0.1982450
           -0.2532636    0.1857116   -1.0695868
   10 Cl   -1.2894875    0.0215655    0.2844782
            0.0019234   -0.8079604   -0.0962407
            0.2874628   -0.1264769   -1.0695868
   11 Cl   -0.9181712    0.2124562   -0.0588922
            0.1928141   -1.1792767    0.2944857
           -0.0341992    0.3121885   -1.0695868
   12 Cl   -0.9181712   -0.2124562   -0.0588922
           -0.1928141   -1.1792767   -0.2944857
           -0.0341992   -0.3121885   -1.0695868
   13 Cl   -0.9385132    0.2045585   -0.2255860
            0.2242007   -1.1589348    0.1982450
           -0.2532636    0.1857116   -1.0695868
   14 Cl   -1.2894875   -0.0215655    0.2844782
           -0.0019234   -0.8079604    0.0962407
            0.2874628    0.1264769   -1.0695868
   15 Cl   -1.2894875    0.0215655    0.2844782
            0.0019234   -0.8079604   -0.0962407
            0.2874628   -0.1264769   -1.0695868
   16 Cl   -0.9385132   -0.2045585   -0.2255860
           -0.2242007   -1.1589348   -0.1982450
           -0.2532636   -0.1857116   -1.0695868
   17 Cl   -0.9385132   -0.2045585   -0.2255860
           -0.2242007   -1.1589348   -0.1982450
           -0.2532636   -0.1857116   -1.0695868
   18 Cl   -0.9181712    0.2124562   -0.0588922
            0.1928141   -1.1792767    0.2944857
           -0.0341992    0.3121885   -1.0695868
   19 In    2.8822955   -0.0000000   -0.0000000
            0.0000000    2.8822955   -0.0000000
           -0.0000000   -0.0000000    3.0514539
   20 In    2.8822955   -0.0000000   -0.0000000
            0.0000000    2.8822955   -0.0000000
           -0.0000000   -0.0000000    3.0514539
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 9, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 25, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [-0.0050  0.0087  0.0000]
    [ 0.0050 -0.0087  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 73, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000045 (xxy) 0.00000045 (xxy) 0.00000045 (xyx)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000009 (zy) -0.00000009 (zy) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 05:28:36]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 05:28:36]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [2 2 1]
Primitive matrix:
  [ 1. -0.  0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P-31c (163)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    6.893806931178462    0.000000000000000    0.000000000000000
  b   -3.446903465589231    5.970211931185789    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
    1 Na  0.33333333333333  0.66666666666667  0.04001250876844  22.990
    2 Na  0.33333333333333  0.66666666666667  0.45998749123156  22.990
    3 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990
    4 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990
    5 Na  0.66666666666667  0.33333333333333  0.95998749123156  22.990
    6 Na  0.66666666666667  0.33333333333333  0.54001250876844  22.990
    7 Cl  0.36395380808473  0.04835563217144  0.12803350975519  35.453
    8 Cl  0.95164436782856  0.31559817591329  0.12803350975519  35.453
    9 Cl  0.68440182408671  0.63604619191527  0.12803350975519  35.453
   10 Cl  0.63604619191527  0.68440182408671  0.62803350975519  35.453
   11 Cl  0.31559817591329  0.95164436782856  0.62803350975519  35.453
   12 Cl  0.63604619191527  0.95164436782856  0.87196649024481  35.453
   13 Cl  0.31559817591329  0.36395380808473  0.87196649024481  35.453
   14 Cl  0.04835563217144  0.68440182408671  0.87196649024481  35.453
   15 Cl  0.36395380808473  0.31559817591329  0.37196649024481  35.453
   16 Cl  0.95164436782856  0.63604619191527  0.37196649024481  35.453
   17 Cl  0.04835563217144  0.36395380808473  0.62803350975519  35.453
   18 Cl  0.68440182408671  0.04835563217144  0.37196649024481  35.453
   19 In  0.33333333333333  0.66666666666667  0.75000000000000 114.818
   20 In  0.66666666666667  0.33333333333333  0.25000000000000 114.818
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.787613862356924    0.000000000000000    0.000000000000000
  b   -6.893806931178462   11.940423862371578    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.338465540000000
Atomic positions (fractional):
    1 Na  0.16666666666667  0.33333333333333  0.04001250876844  22.990 > 1
    2 Na  0.66666666666667  0.33333333333333  0.04001250876844  22.990 > 1
    3 Na  0.16666666666667  0.83333333333333  0.04001250876844  22.990 > 1
    4 Na  0.66666666666667  0.83333333333333  0.04001250876844  22.990 > 1
    5 Na  0.16666666666667  0.33333333333333  0.45998749123156  22.990 > 5
    6 Na  0.66666666666667  0.33333333333333  0.45998749123156  22.990 > 5
    7 Na  0.16666666666667  0.83333333333333  0.45998749123156  22.990 > 5
    8 Na  0.66666666666667  0.83333333333333  0.45998749123156  22.990 > 5
    9 Na  0.00000000000000  0.00000000000000  0.75000000000000  22.990 > 9
   10 Na  0.50000000000000  0.00000000000000  0.75000000000000  22.990 > 9
   11 Na  0.00000000000000  0.50000000000000  0.75000000000000  22.990 > 9
   12 Na  0.50000000000000  0.50000000000000  0.75000000000000  22.990 > 9
   13 Na  0.00000000000000  0.00000000000000  0.25000000000000  22.990 > 13
   14 Na  0.50000000000000  0.00000000000000  0.25000000000000  22.990 > 13
   15 Na  0.00000000000000  0.50000000000000  0.25000000000000  22.990 > 13
   16 Na  0.50000000000000  0.50000000000000  0.25000000000000  22.990 > 13
   17 Na  0.33333333333333  0.16666666666667  0.95998749123156  22.990 > 17
   18 Na  0.83333333333333  0.16666666666667  0.95998749123156  22.990 > 17
   19 Na  0.33333333333333  0.66666666666667  0.95998749123156  22.990 > 17
   20 Na  0.83333333333333  0.66666666666667  0.95998749123156  22.990 > 17
   21 Na  0.33333333333333  0.16666666666667  0.54001250876844  22.990 > 21
   22 Na  0.83333333333333  0.16666666666667  0.54001250876844  22.990 > 21
   23 Na  0.33333333333333  0.66666666666667  0.54001250876844  22.990 > 21
   24 Na  0.83333333333333  0.66666666666667  0.54001250876844  22.990 > 21
   25 Cl  0.18197690404237  0.02417781608572  0.12803350975519  35.453 > 25
   26 Cl  0.68197690404237  0.02417781608572  0.12803350975519  35.453 > 25
   27 Cl  0.18197690404237  0.52417781608572  0.12803350975519  35.453 > 25
   28 Cl  0.68197690404237  0.52417781608572  0.12803350975519  35.453 > 25
   29 Cl  0.47582218391428  0.15779908795665  0.12803350975519  35.453 > 29
   30 Cl  0.97582218391428  0.15779908795665  0.12803350975519  35.453 > 29
   31 Cl  0.47582218391428  0.65779908795665  0.12803350975519  35.453 > 29
   32 Cl  0.97582218391428  0.65779908795665  0.12803350975519  35.453 > 29
   33 Cl  0.34220091204335  0.31802309595763  0.12803350975519  35.453 > 33
   34 Cl  0.84220091204335  0.31802309595763  0.12803350975519  35.453 > 33
   35 Cl  0.34220091204335  0.81802309595763  0.12803350975519  35.453 > 33
   36 Cl  0.84220091204335  0.81802309595763  0.12803350975519  35.453 > 33
   37 Cl  0.31802309595763  0.34220091204335  0.62803350975519  35.453 > 37
   38 Cl  0.81802309595763  0.34220091204335  0.62803350975519  35.453 > 37
   39 Cl  0.31802309595763  0.84220091204335  0.62803350975519  35.453 > 37
   40 Cl  0.81802309595763  0.84220091204335  0.62803350975519  35.453 > 37
   41 Cl  0.15779908795665  0.47582218391428  0.62803350975519  35.453 > 41
   42 Cl  0.65779908795665  0.47582218391428  0.62803350975519  35.453 > 41
   43 Cl  0.15779908795665  0.97582218391428  0.62803350975519  35.453 > 41
   44 Cl  0.65779908795665  0.97582218391428  0.62803350975519  35.453 > 41
   45 Cl  0.31802309595763  0.47582218391428  0.87196649024481  35.453 > 45
   46 Cl  0.81802309595763  0.47582218391428  0.87196649024481  35.453 > 45
   47 Cl  0.31802309595763  0.97582218391428  0.87196649024481  35.453 > 45
   48 Cl  0.81802309595763  0.97582218391428  0.87196649024481  35.453 > 45
   49 Cl  0.15779908795665  0.18197690404237  0.87196649024481  35.453 > 49
   50 Cl  0.65779908795665  0.18197690404237  0.87196649024481  35.453 > 49
   51 Cl  0.15779908795665  0.68197690404237  0.87196649024481  35.453 > 49
   52 Cl  0.65779908795665  0.68197690404237  0.87196649024481  35.453 > 49
   53 Cl  0.02417781608572  0.34220091204335  0.87196649024481  35.453 > 53
   54 Cl  0.52417781608572  0.34220091204335  0.87196649024481  35.453 > 53
   55 Cl  0.02417781608572  0.84220091204335  0.87196649024481  35.453 > 53
   56 Cl  0.52417781608572  0.84220091204335  0.87196649024481  35.453 > 53
   57 Cl  0.18197690404237  0.15779908795665  0.37196649024481  35.453 > 57
   58 Cl  0.68197690404237  0.15779908795665  0.37196649024481  35.453 > 57
   59 Cl  0.18197690404237  0.65779908795665  0.37196649024481  35.453 > 57
   60 Cl  0.68197690404237  0.65779908795665  0.37196649024481  35.453 > 57
   61 Cl  0.47582218391428  0.31802309595763  0.37196649024481  35.453 > 61
   62 Cl  0.97582218391428  0.31802309595763  0.37196649024481  35.453 > 61
   63 Cl  0.47582218391428  0.81802309595763  0.37196649024481  35.453 > 61
   64 Cl  0.97582218391428  0.81802309595763  0.37196649024481  35.453 > 61
   65 Cl  0.02417781608572  0.18197690404237  0.62803350975519  35.453 > 65
   66 Cl  0.52417781608572  0.18197690404237  0.62803350975519  35.453 > 65
   67 Cl  0.02417781608572  0.68197690404237  0.62803350975519  35.453 > 65
   68 Cl  0.52417781608572  0.68197690404237  0.62803350975519  35.453 > 65
   69 Cl  0.34220091204335  0.02417781608572  0.37196649024481  35.453 > 69
   70 Cl  0.84220091204335  0.02417781608572  0.37196649024481  35.453 > 69
   71 Cl  0.34220091204335  0.52417781608572  0.37196649024481  35.453 > 69
   72 Cl  0.84220091204335  0.52417781608572  0.37196649024481  35.453 > 69
   73 In  0.16666666666667  0.33333333333333  0.75000000000000 114.818 > 73
   74 In  0.66666666666667  0.33333333333333  0.75000000000000 114.818 > 73
   75 In  0.16666666666667  0.83333333333333  0.75000000000000 114.818 > 73
   76 In  0.66666666666667  0.83333333333333  0.75000000000000 114.818 > 73
   77 In  0.33333333333333  0.16666666666667  0.25000000000000 114.818 > 77
   78 In  0.83333333333333  0.16666666666667  0.25000000000000 114.818 > 77
   79 In  0.33333333333333  0.66666666666667  0.25000000000000 114.818 > 77
   80 In  0.83333333333333  0.66666666666667  0.25000000000000 114.818 > 77
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.7129661   -0.0000000    0.0000000
            0.0000000    2.7129661    0.0000000
            0.0000000    0.0000000    2.7894726
-------------------------- Born effective charges --------------------------
    1 Na    1.0949434    0.0922169   -0.0000000
           -0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    2 Na    1.0949434   -0.0922169   -0.0000000
            0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    3 Na    1.2201615    0.0000000   -0.0000000
           -0.0000000    1.2201615   -0.0000000
           -0.0000000   -0.0000000    0.8640476
    4 Na    1.2201615    0.0000000   -0.0000000
           -0.0000000    1.2201615   -0.0000000
           -0.0000000   -0.0000000    0.8640476
    5 Na    1.0949434    0.0922169   -0.0000000
           -0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    6 Na    1.0949434   -0.0922169   -0.0000000
            0.0922169    1.0949434   -0.0000000
           -0.0000000   -0.0000000    1.2510097
    7 Cl   -0.9181712   -0.2124562   -0.0588922
           -0.1928141   -1.1792767   -0.2944857
           -0.0341992   -0.3121885   -1.0695868
    8 Cl   -1.2894875   -0.0215655    0.2844782
           -0.0019234   -0.8079604    0.0962407
            0.2874628    0.1264769   -1.0695868
    9 Cl   -0.9385132    0.2045585   -0.2255860
            0.2242007   -1.1589348    0.1982450
           -0.2532636    0.1857116   -1.0695868
   10 Cl   -1.2894875    0.0215655    0.2844782
            0.0019234   -0.8079604   -0.0962407
            0.2874628   -0.1264769   -1.0695868
   11 Cl   -0.9181712    0.2124562   -0.0588922
            0.1928141   -1.1792767    0.2944857
           -0.0341992    0.3121885   -1.0695868
   12 Cl   -0.9181712   -0.2124562   -0.0588922
           -0.1928141   -1.1792767   -0.2944857
           -0.0341992   -0.3121885   -1.0695868
   13 Cl   -0.9385132    0.2045585   -0.2255860
            0.2242007   -1.1589348    0.1982450
           -0.2532636    0.1857116   -1.0695868
   14 Cl   -1.2894875   -0.0215655    0.2844782
           -0.0019234   -0.8079604    0.0962407
            0.2874628    0.1264769   -1.0695868
   15 Cl   -1.2894875    0.0215655    0.2844782
            0.0019234   -0.8079604   -0.0962407
            0.2874628   -0.1264769   -1.0695868
   16 Cl   -0.9385132   -0.2045585   -0.2255860
           -0.2242007   -1.1589348   -0.1982450
           -0.2532636   -0.1857116   -1.0695868
   17 Cl   -0.9385132   -0.2045585   -0.2255860
           -0.2242007   -1.1589348   -0.1982450
           -0.2532636   -0.1857116   -1.0695868
   18 Cl   -0.9181712    0.2124562   -0.0588922
            0.1928141   -1.1792767    0.2944857
           -0.0341992    0.3121885   -1.0695868
   19 In    2.8822955   -0.0000000   -0.0000000
            0.0000000    2.8822955   -0.0000000
           -0.0000000   -0.0000000    3.0514539
   20 In    2.8822955   -0.0000000   -0.0000000
            0.0000000    2.8822955   -0.0000000
           -0.0000000   -0.0000000    3.0514539
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000045 (xxy) 0.00000045 (xxy) 0.00000045 (xyx)
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 ... [ 8 8 4 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.52, Number of G-points: 313, Lambda: 0.09
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/35) =======================
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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 35
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
   1.166   (   0.000    0.000    0.000)    0.000
   1.166   (   0.000    0.000    0.000)    0.000
   1.240   (   0.000    0.000    0.000)    0.000
   1.646   (   0.000    0.000    0.000)    0.000
   1.646   (   0.000    0.000    0.000)    0.000
   1.938   (   0.000    0.000    0.000)    0.000
   2.362   (   0.000    0.000    0.000)    0.000
   2.362   (   0.000    0.000    0.000)    0.000
   2.567   (   0.000    0.000    0.000)    0.000
   2.814   (   0.000    0.000    0.000)    0.000
   3.004   (   0.000    0.000    0.000)    0.000
   3.004   (   0.000    0.000    0.000)    0.000
   3.225   (   0.000    0.000    0.000)    0.000
   3.237   (   0.000    0.000    0.000)    0.000
   3.237   (   0.000    0.000    0.000)    0.000
   3.264   (   0.000    0.000    0.000)    0.000
   3.322   (   0.000    0.000    0.000)    0.000
   3.322   (   0.000    0.000    0.000)    0.000
   3.346   (   0.000    0.000    0.000)    0.000
   3.381   (   0.000    0.000    0.000)    0.000
   3.381   (   0.000    0.000    0.000)    0.000
   3.581   (   0.000    0.000    0.000)    0.000
   3.581   (   0.000    0.000    0.000)    0.000
   3.819   (   0.000    0.000    0.000)    0.000
   3.819   (   0.000    0.000    0.000)    0.000
   4.047   (   0.000    0.000    0.000)    0.000
   4.123   (   0.000    0.000    0.000)    0.000
   4.225   (   0.000    0.000    0.000)    0.000
   4.687   (   0.000    0.000    0.000)    0.000
   4.786   (   0.000    0.000    0.000)    0.000
   4.786   (   0.000    0.000    0.000)    0.000
   4.826   (   0.000    0.000    0.000)    0.000
   4.826   (   0.000    0.000    0.000)    0.000
   4.896   (   0.000    0.000    0.000)    0.000
   4.896   (   0.000    0.000    0.000)    0.000
   4.936   (   0.000    0.000    0.000)    0.000
   4.936   (   0.000    0.000    0.000)    0.000
   5.128   (   0.000    0.000    0.000)    0.000
   5.376   (   0.000    0.000    0.000)    0.000
   5.470   (   0.000    0.000    0.000)    0.000
   5.470   (   0.000    0.000    0.000)    0.000
   5.659   (   0.000    0.000    0.000)    0.000
   5.659   (   0.000    0.000    0.000)    0.000
   5.804   (   0.000    0.000    0.000)    0.000
   5.929   (   0.000    0.000    0.000)    0.000
   5.929   (   0.000    0.000    0.000)    0.000
   5.930   (   0.000    0.000    0.000)    0.000
   5.998   (   0.000    0.000    0.000)    0.000
   5.998   (   0.000    0.000    0.000)    0.000
   7.416   (   0.000    0.000    0.000)    0.000
   7.416   (   0.000    0.000    0.000)    0.000
   7.421   (   0.000    0.000    0.000)    0.000
   7.421   (   0.000    0.000    0.000)    0.000
   7.700   (   0.000    0.000    0.000)    0.000
   7.920   (   0.000    0.000    0.000)    0.000
   7.977   (   0.000    0.000    0.000)    0.000
   8.521   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/35) =======================
q-point: ( 0.12  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 76
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.313   (  12.578    7.262    0.000)   14.524
   0.377   (  15.082    8.708    0.000)   17.415
   0.703   (  27.246   15.730    0.000)   31.461
   1.183   (  -0.490   -0.283    0.000)    0.566
   1.208   (  -0.646   -0.373    0.000)    0.746
   1.252   (   6.470    3.735    0.000)    7.470
   1.550   (  -6.409   -3.700    0.000)    7.401
   1.699   (   1.306    0.754    0.000)    1.508
   1.922   (  -1.312   -0.758    0.000)    1.515
   2.364   (  -0.206   -0.119    0.000)    0.238
   2.366   (   0.222    0.128    0.000)    0.256
   2.591   (  -0.183   -0.105    0.000)    0.211
   2.687   (  -7.920   -4.572    0.000)    9.145
   2.924   (  -6.910   -3.989    0.000)    7.979
   2.981   (  -9.683   -5.590    0.000)   11.181
   3.001   (  -9.968   -5.755    0.000)   11.510
   3.128   (  -0.130   -0.075    0.000)    0.150
   3.135   (  -1.572   -0.908    0.000)    1.815
   3.260   (  -1.425   -0.823    0.000)    1.645
   3.348   (  -0.888   -0.513    0.000)    1.025
   3.389   (   1.479    0.854    0.000)    1.708
   3.423   (   5.533    3.194    0.000)    6.389
   3.501   (  -0.954   -0.551    0.000)    1.102
   3.542   (   4.736    2.734    0.000)    5.468
   3.598   (   4.780    2.760    0.000)    5.519
   3.603   (   4.641    2.679    0.000)    5.359
   3.808   (   0.884    0.511    0.000)    1.021
   3.845   (   2.115    1.221    0.000)    2.442
   4.025   (  -1.532   -0.884    0.000)    1.769
   4.159   (   3.376    1.949    0.000)    3.898
   4.231   (   0.699    0.404    0.000)    0.807
   4.704   (   0.588    0.339    0.000)    0.679
   4.718   (  -4.591   -2.650    0.000)    5.301
   4.767   (  -2.961   -1.710    0.000)    3.420
   4.812   (   0.446    0.257    0.000)    0.515
   4.877   (  -2.971   -1.716    0.000)    3.431
   4.948   (   3.715    2.145    0.000)    4.289
   4.949   (   1.176    0.679    0.000)    1.358
   4.965   (   2.240    1.294    0.000)    2.587
   5.146   (   3.998    2.308    0.000)    4.616
   5.256   (   7.540    4.353    0.000)    8.706
   5.372   (  -0.466   -0.269    0.000)    0.538
   5.470   (   0.003    0.002    0.000)    0.003
   5.579   (  -5.223   -3.015    0.000)    6.031
   5.639   (  -6.770   -3.909    0.000)    7.818
   5.666   (  -1.817   -1.049    0.000)    2.098
   5.796   (   1.331    0.769    0.000)    1.537
   5.914   (  -0.842   -0.486    0.000)    0.972
   5.951   (   1.849    1.067    0.000)    2.135
   5.957   (   1.546    0.893    0.000)    1.785
   6.149   (   8.192    4.730    0.000)    9.459
   6.431   (  -1.008   -0.582    0.000)    1.164
   7.410   (   0.105    0.061    0.000)    0.121
   7.453   (   1.907    1.101    0.000)    2.202
   7.456   (   2.754    1.590    0.000)    3.180
   7.685   (  -1.199   -0.692    0.000)    1.385
   7.921   (   0.088    0.051    0.000)    0.102
   7.984   (   0.558    0.322    0.000)    0.644
   8.406   (  -2.592   -1.496    0.000)    2.993
   8.426   (  -7.460   -4.307    0.000)    8.615
======================= Grid point 2 (3/35) =======================
q-point: ( 0.25  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 78
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.597   (  10.610    6.126    0.000)   12.251
   0.714   (  12.375    7.145    0.000)   14.289
   1.115   (  -2.961   -1.710    0.000)    3.420
   1.231   (   0.007    0.004    0.000)    0.008
   1.243   (  15.178    8.763    0.000)   17.526
   1.434   (  -1.453   -0.839    0.000)    1.678
   1.436   (   8.042    4.643    0.000)    9.286
   1.787   (   7.044    4.067    0.000)    8.134
   1.876   (  -2.339   -1.351    0.000)    2.701
   2.297   (  -6.785   -3.917    0.000)    7.834
   2.343   (  -3.626   -2.094    0.000)    4.187
   2.477   (  -8.361   -4.827    0.000)    9.654
   2.516   (  -4.462   -2.576    0.000)    5.152
   2.689   ( -12.581   -7.264    0.000)   14.528
   2.735   (  -5.022   -2.899    0.000)    5.799
   2.870   (  -1.715   -0.990    0.000)    1.981
   3.135   (   0.406    0.234    0.000)    0.468
   3.172   (   0.443    0.256    0.000)    0.511
   3.183   (  -3.534   -2.040    0.000)    4.080
   3.293   (  -3.266   -1.886    0.000)    3.771
   3.368   (  -1.896   -1.094    0.000)    2.189
   3.486   (   0.399    0.230    0.000)    0.460
   3.524   (   3.196    1.845    0.000)    3.690
   3.613   (   3.843    2.219    0.000)    4.438
   3.718   (   4.261    2.460    0.000)    4.920
   3.763   (   5.966    3.444    0.000)    6.889
   3.890   (   0.193    0.111    0.000)    0.222
   3.913   (   7.611    4.394    0.000)    8.789
   3.998   (   0.278    0.161    0.000)    0.321
   4.271   (   2.631    1.519    0.000)    3.039
   4.296   (   7.004    4.044    0.000)    8.087
   4.597   (  -5.569   -3.216    0.000)    6.431
   4.639   (  -5.013   -2.894    0.000)    5.789
   4.724   (  -0.404   -0.233    0.000)    0.467
   4.781   (  -4.543   -2.623    0.000)    5.246
   4.830   (   0.426    0.246    0.000)    0.492
   4.988   (   1.829    1.056    0.000)    2.112
   5.015   (   1.973    1.139    0.000)    2.278
   5.041   (   2.609    1.506    0.000)    3.012
   5.235   (   2.347    1.355    0.000)    2.711
   5.339   (  -2.668   -1.540    0.000)    3.081
   5.429   (   4.162    2.403    0.000)    4.806
   5.442   (  -6.737   -3.889    0.000)    7.779
   5.469   (  -0.952   -0.550    0.000)    1.100
   5.470   (  -0.464   -0.268    0.000)    0.536
   5.678   (   1.191    0.687    0.000)    1.375
   5.816   (  -0.884   -0.510    0.000)    1.021
   5.902   (   1.252    0.723    0.000)    1.445
   5.999   (   1.795    1.036    0.000)    2.073
   6.014   (   3.139    1.813    0.000)    3.625
   6.329   (   6.257    3.613    0.000)    7.225
   6.395   (  -1.906   -1.101    0.000)    2.201
   7.421   (   0.728    0.420    0.000)    0.841
   7.508   (   2.518    1.454    0.000)    2.908
   7.546   (   4.427    2.556    0.000)    5.111
   7.652   (  -1.265   -0.730    0.000)    1.461
   7.922   (  -0.053   -0.031    0.000)    0.062
   8.003   (   0.981    0.567    0.000)    1.133
   8.193   ( -10.988   -6.344    0.000)   12.688
   8.318   (  -4.445   -2.566    0.000)    5.133
======================= Grid point 3 (4/35) =======================
q-point: ( 0.38  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 76
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.815   (   7.094    4.096    0.000)    8.191
   0.944   (   5.706    3.294    0.000)    6.588
   1.082   (   1.474    0.851    0.000)    1.702
   1.220   (  -0.209   -0.121    0.000)    0.241
   1.421   (   2.333    1.347    0.000)    2.694
   1.442   (  -0.503   -0.290    0.000)    0.580
   1.610   (   5.261    3.037    0.000)    6.075
   1.843   (   1.023    0.590    0.000)    1.181
   1.988   (   5.792    3.344    0.000)    6.688
   2.081   (  -8.491   -4.903    0.000)    9.805
   2.194   (  -5.877   -3.393    0.000)    6.786
   2.247   (  -9.859   -5.692    0.000)   11.384
   2.475   (  -3.071   -1.773    0.000)    3.546
   2.489   (  -0.515   -0.298    0.000)    0.595
   2.711   (   1.554    0.897    0.000)    1.795
   2.827   (  -2.151   -1.242    0.000)    2.484
   3.040   (  -4.960   -2.863    0.000)    5.727
   3.139   (   0.199    0.115    0.000)    0.229
   3.220   (  -2.579   -1.489    0.000)    2.978
   3.240   (   2.164    1.249    0.000)    2.499
   3.355   (   1.765    1.019    0.000)    2.038
   3.483   (  -0.291   -0.168    0.000)    0.336
   3.604   (   3.080    1.778    0.000)    3.556
   3.707   (   3.133    1.809    0.000)    3.618
   3.780   (  -0.302   -0.175    0.000)    0.349
   3.836   (  -3.044   -1.758    0.000)    3.515
   3.891   (   5.343    3.085    0.000)    6.170
   4.054   (   3.500    2.021    0.000)    4.042
   4.133   (   7.731    4.463    0.000)    8.927
   4.332   (   0.716    0.413    0.000)    0.827
   4.371   (  -0.725   -0.419    0.000)    0.838
   4.475   (  -2.478   -1.431    0.000)    2.861
   4.625   (   4.405    2.543    0.000)    5.086
   4.668   (  -4.656   -2.688    0.000)    5.376
   4.685   (  -3.011   -1.739    0.000)    3.477
   4.832   (  -0.183   -0.105    0.000)    0.211
   5.022   (   0.705    0.407    0.000)    0.814
   5.028   (  -3.353   -1.936    0.000)    3.872
   5.060   (   1.557    0.899    0.000)    1.798
   5.232   (  -1.545   -0.892    0.000)    1.785
   5.256   (  -3.352   -1.935    0.000)    3.870
   5.335   (  -1.947   -1.124    0.000)    2.249
   5.459   (   2.812    1.623    0.000)    3.247
   5.465   (  -0.189   -0.109    0.000)    0.219
   5.558   (   2.932    1.693    0.000)    3.386
   5.709   (   1.104    0.637    0.000)    1.275
   5.751   (  -3.371   -1.946    0.000)    3.892
   5.975   (   3.436    1.984    0.000)    3.968
   6.037   (   1.272    0.734    0.000)    1.468
   6.087   (   2.505    1.446    0.000)    2.893
   6.344   (  -2.093   -1.208    0.000)    2.417
   6.443   (   3.143    1.814    0.000)    3.629
   7.440   (   0.688    0.397    0.000)    0.795
   7.564   (   1.903    1.099    0.000)    2.197
   7.637   (   0.230    0.133    0.000)    0.266
   7.655   (   4.213    2.432    0.000)    4.865
   7.915   (  -0.544   -0.314    0.000)    0.628
   7.940   (  -9.055   -5.228    0.000)   10.456
   8.029   (   1.136    0.656    0.000)    1.311
   8.206   (  -4.360   -2.517    0.000)    5.034
======================= Grid point 4 (5/35) =======================
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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 46
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.902   (   0.000    0.000    0.000)    0.000
   1.006   (   0.000    0.000    0.000)    0.000
   1.171   (  -0.000   -0.000    0.000)    0.000
   1.187   (   0.000    0.000    0.000)    0.000
   1.401   (   0.000    0.000    0.000)    0.000
   1.443   (   0.000    0.000    0.000)    0.000
   1.671   (   0.000    0.000    0.000)    0.000
   1.940   (  -0.000   -0.000    0.000)    0.000
   1.943   (   0.000    0.000    0.000)    0.000
   2.014   (   0.000    0.000    0.000)    0.000
   2.075   (   0.000    0.000    0.000)    0.000
   2.145   (   0.000    0.000    0.000)    0.000
   2.446   (   0.000    0.000    0.000)    0.000
   2.499   (   0.000    0.000    0.000)    0.000
   2.756   (   0.000    0.000    0.000)    0.000
   2.777   (  -0.000   -0.000    0.000)    0.000
   2.975   (   0.000    0.000    0.000)    0.000
   3.154   (   0.000    0.000    0.000)    0.000
   3.175   (   0.000    0.000    0.000)    0.000
   3.256   (   0.000    0.000    0.000)    0.000
   3.420   (   0.000    0.000    0.000)    0.000
   3.479   (   0.000    0.000    0.000)    0.000
   3.677   (   0.000    0.000    0.000)    0.000
   3.688   (  -0.000   -0.000    0.000)    0.000
   3.746   (   0.000    0.000    0.000)    0.000
   3.817   (   0.000    0.000    0.000)    0.000
   3.994   (   0.000    0.000    0.000)    0.000
   4.114   (   0.000    0.000    0.000)    0.000
   4.152   (   0.000    0.000    0.000)    0.000
   4.299   (  -0.000   -0.000    0.000)    0.000
   4.495   (   0.000    0.000    0.000)    0.000
   4.500   (   0.000    0.000    0.000)    0.000
   4.519   (   0.000    0.000    0.000)    0.000
   4.583   (   0.000    0.000    0.000)    0.000
   4.746   (   0.000    0.000    0.000)    0.000
   4.815   (  -0.000   -0.000    0.000)    0.000
   4.964   (   0.000    0.000    0.000)    0.000
   5.027   (   0.000    0.000    0.000)    0.000
   5.094   (   0.000    0.000    0.000)    0.000
   5.195   (   0.000    0.000    0.000)    0.000
   5.211   (   0.000    0.000    0.000)    0.000
   5.322   (   0.000    0.000    0.000)    0.000
   5.462   (   0.000    0.000    0.000)    0.000
   5.514   (   0.000    0.000    0.000)    0.000
   5.587   (   0.000    0.000    0.000)    0.000
   5.699   (   0.000    0.000    0.000)    0.000
   5.723   (   0.000    0.000    0.000)    0.000
   6.022   (   0.000    0.000    0.000)    0.000
   6.056   (   0.000    0.000    0.000)    0.000
   6.119   (   0.000    0.000    0.000)    0.000
   6.313   (   0.000    0.000    0.000)    0.000
   6.481   (   0.000    0.000    0.000)    0.000
   7.449   (   0.000    0.000    0.000)    0.000
   7.589   (   0.000    0.000    0.000)    0.000
   7.653   (   0.000    0.000    0.000)    0.000
   7.716   (   0.000    0.000    0.000)    0.000
   7.810   (   0.000    0.000    0.000)    0.000
   7.904   (   0.000    0.000    0.000)    0.000
   8.049   (   0.000    0.000    0.000)    0.000
   8.142   (   0.000    0.000    0.000)    0.000
======================= Grid point 10 (6/35) =======================
q-point: ( 0.12  0.12  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 76
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.573   (   7.536   13.052    2.050)   15.210
   0.600   (   7.370   12.766    2.833)   15.011
   1.121   (  10.954   18.973   -1.799)   21.982
   1.143   (  -1.747   -3.026   -1.226)    3.702
   1.240   (   1.956    3.389   -3.124)    5.007
   1.350   (   1.336    2.315   10.600)   10.931
   1.467   (  -0.230   -0.399  -11.022)   11.032
   1.744   (   2.606    4.513    4.641)    6.979
   1.894   (  -1.109   -1.921    2.026)    3.003
   2.350   (  -1.579   -2.734   -0.234)    3.166
   2.354   (  -0.684   -1.184   -5.827)    5.985
   2.521   (  -4.453   -7.713    2.489)    9.247
   2.538   (  -4.868   -8.432   -2.536)   10.061
   2.789   (  -4.343   -7.523    0.739)    8.718
   2.816   (  -4.196   -7.267   -0.502)    8.407
   2.841   (  -6.296  -10.904    4.322)   13.313
   3.071   (  -0.033   -0.057   -0.181)    0.192
   3.100   (  -0.854   -1.479    0.357)    1.745
   3.274   (  -0.351   -0.609    4.461)    4.516
   3.344   (   0.298    0.517   -3.526)    3.576
   3.452   (  -2.951   -5.111    3.162)    6.696
   3.478   (   1.152    1.996   -1.127)    2.566
   3.522   (   1.588    2.751   -4.995)    5.920
   3.574   (   3.200    5.543   -4.001)    7.548
   3.608   (   2.039    3.531   -3.083)    5.112
   3.752   (   4.300    7.448    2.964)    9.097
   3.837   (   2.430    4.209   -4.736)    6.786
   3.916   (   1.782    3.086    7.995)    8.753
   3.999   (  -0.548   -0.948    0.405)    1.168
   4.250   (   0.938    1.624   -0.473)    1.934
   4.262   (   4.685    8.114    1.607)    9.506
   4.588   (  -3.958   -6.855    0.018)    7.916
   4.732   (  -1.891   -3.274    3.076)    4.874
   4.753   (  -1.973   -3.417    1.088)    4.093
   4.814   (   0.589    1.020   -0.073)    1.180
   4.846   (   1.907    3.303    2.282)    4.445
   4.884   (  -0.836   -1.448   -3.964)    4.303
   4.971   (   0.870    1.506   -2.196)    2.801
   5.145   (   3.508    6.077   -2.130)    7.333
   5.150   (  -0.081   -0.140   -0.981)    0.995
   5.372   (   0.484    0.838   -0.806)    1.259
   5.397   (   3.388    5.868    4.123)    7.932
   5.441   (  -1.098   -1.902   -2.097)    3.037
   5.491   (  -4.814   -8.338   -0.069)    9.628
   5.558   (  -1.858   -3.218    3.306)    4.974
   5.568   (  -1.112   -1.926   -0.430)    2.265
   5.872   (  -1.587   -2.749   -2.782)    4.220
   5.891   (   3.287    5.694   -3.809)    7.598
   5.981   (   0.790    1.369    2.879)    3.285
   5.989   (   1.394    2.414   -2.215)    3.560
   6.280   (   3.729    6.458   -2.225)    7.782
   6.412   (  -0.684   -1.184    8.799)    8.904
   7.421   (   0.403    0.698    1.202)    1.447
   7.477   (   1.532    2.653   -0.930)    3.201
   7.530   (   2.337    4.048    0.822)    4.746
   7.660   (  -0.828   -1.435    0.393)    1.702
   7.922   (   0.004    0.007   -0.165)    0.166
   7.997   (   0.517    0.895   -0.036)    1.034
   8.263   (  -6.074  -10.521    0.690)   12.168
   8.345   (  -2.425   -4.201   -3.767)    6.142
======================= Grid point 11 (7/35) =======================
q-point: ( 0.25  0.12  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.801   (   3.915   11.590    2.262)   12.441
   0.874   (   8.764    7.883    1.648)   11.903
   1.089   (  -1.333   -0.808   -0.094)    1.562
   1.260   (  -2.264    3.095   -0.632)    3.886
   1.357   (   3.593   -0.061    0.963)    3.720
   1.421   (   5.799    3.247    0.872)    6.703
   1.556   (   3.538    4.415   -5.112)    7.625
   1.841   (  -0.807   -1.298    0.097)    1.532
   1.910   (   5.939    5.699    6.665)   10.591
   2.208   (  -8.799   -4.520   -2.174)   10.128
   2.266   (  -5.115   -6.435  -10.054)   12.987
   2.309   (  -7.847   -7.205    4.598)   11.603
   2.423   (   1.529   -7.536   -1.808)    7.899
   2.561   (  -8.698   -3.926    4.303)   10.469
   2.685   (  -1.105   -5.245    0.283)    5.368
   2.768   (   2.856   -5.480    1.653)    6.397
   3.045   (  -2.404   -3.843   -1.413)    4.748
   3.119   (   4.167   -0.670   -1.808)    4.591
   3.240   (  -3.965    4.164    5.511)    7.965
   3.320   (  -1.969   -1.191   -4.806)    5.329
   3.359   (  -1.973   -0.909    0.572)    2.247
   3.517   (   1.511    1.955    0.283)    2.487
   3.559   (  -0.340    3.485   -1.507)    3.812
   3.648   (   1.539    1.401    1.437)    2.529
   3.747   (   8.128    0.232   -9.509)   12.511
   3.860   (  -2.192    5.104    3.854)    6.761
   3.954   (  -0.128    6.071    3.030)    6.786
   4.010   (   1.555    1.691    0.044)    2.297
   4.033   (  10.342    7.611    1.648)   12.946
   4.304   (   2.859    1.803   -0.247)    3.389
   4.418   (  -2.645   -1.499    2.621)    4.014
   4.452   (  -0.534   -0.074   -0.798)    0.964
   4.621   (  -1.656    0.291    1.871)    2.516
   4.684   (  -1.343   -5.560    2.198)    6.128
   4.825   (  -6.831    5.529    0.891)    8.833
   4.842   (   0.162    1.154   -4.212)    4.370
   4.948   (   4.929    0.279    1.860)    5.275
   5.006   (   2.463   -3.374   -1.643)    4.489
   5.102   (  -3.715    0.480   -2.208)    4.348
   5.229   (   1.253   -1.113   -1.752)    2.424
   5.312   (  -4.699   -2.455    4.483)    6.943
   5.336   (  -2.885   -5.614   -1.173)    6.420
   5.442   (   1.679    1.927    0.615)    2.629
   5.478   (  -0.310    0.598    1.398)    1.551
   5.527   (   2.071    2.686   -1.955)    3.915
   5.585   (   4.180   -5.093   -0.485)    6.607
   5.794   (  -3.235   -2.923   -4.865)    6.533
   5.991   (   1.459    1.061   -2.423)    3.021
   6.024   (   1.258    6.372    5.883)    8.763
   6.065   (   3.624    3.135   -4.083)    6.295
   6.365   (  -1.574   -0.087    2.518)    2.970
   6.403   (   2.737    1.054    3.675)    4.701
   7.439   (   0.426    1.113    1.784)    2.146
   7.549   (   2.476    3.078   -0.783)    4.027
   7.618   (   2.714    3.458    1.106)    4.533
   7.638   (  -0.071   -0.230   -0.690)    0.731
   7.919   (  -0.280   -0.376   -0.478)    0.669
   8.003   (  -3.215   -4.748    1.048)    5.829
   8.040   (  -3.878   -3.416    0.307)    5.177
   8.240   (  -3.590   -4.787   -3.419)    6.891
======================= Grid point 12 (8/35) =======================
q-point: ( 0.38  0.12  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.937   (   0.154    5.896    2.398)    6.367
   1.026   (   0.401    3.824   -0.646)    3.899
   1.135   (   4.761    4.932    1.862)    7.103
   1.249   (  -3.454    2.586    0.168)    4.318
   1.407   (   2.157    0.368   -0.254)    2.203
   1.463   (  -2.996    0.254    0.312)    3.023
   1.638   (   2.834   -0.559   -1.584)    3.294
   1.859   (   4.381    0.557    1.728)    4.742
   1.959   (  -1.913   -3.612    3.704)    5.516
   2.050   (  -3.854   -0.725   -4.251)    5.783
   2.092   (  -4.507   -6.754   -4.683)    9.373
   2.189   (  -1.094    0.526    4.125)    4.300
   2.346   (   0.097   -6.824   -1.710)    7.036
   2.486   (   0.377   -0.278    1.241)    1.326
   2.647   (   4.507   -8.902   -3.393)   10.539
   2.777   (   2.437   -3.141    3.800)    5.499
   2.960   (  -2.097   -2.203   -2.056)    3.671
   3.120   (  -0.041   -2.659    2.920)    3.950
   3.269   (  -2.847    2.776   -7.250)    8.268
   3.303   (   0.964    3.004    3.014)    4.363
   3.370   (   0.307    5.742   -1.549)    5.955
   3.547   (  -1.099    3.504   -2.858)    4.654
   3.610   (   3.742   -0.700    7.397)    8.320
   3.673   (   1.314    0.001    1.337)    1.875
   3.776   (  -6.799   -0.399    5.557)    8.791
   3.871   (   2.838    1.262  -10.887)   11.321
   4.049   (   0.072    7.364   -3.088)    7.985
   4.084   (   3.393    5.513    3.681)    7.446
   4.211   (  -2.419   -0.618    3.373)    4.196
   4.288   (  -2.740   -2.282   -4.356)    5.629
   4.405   (   2.609    1.851    2.828)    4.270
   4.407   (   5.034   -1.263    0.147)    5.192
   4.600   (  -1.304   -1.564    0.925)    2.237
   4.706   ( -10.197    0.976   -0.379)   10.250
   4.755   (   4.410    6.720    4.260)    9.097
   4.852   (  -1.528    2.094   -2.426)    3.550
   4.976   (  -1.432   -1.514   -0.939)    2.286
   5.007   (   1.892   -1.068    0.133)    2.176
   5.058   (   0.678   -1.354   -1.547)    2.165
   5.187   (  -3.246   -5.215    5.988)    8.578
   5.225   (  -0.287    0.129    1.515)    1.548
   5.294   (   1.413   -2.026   -2.971)    3.863
   5.470   (  -0.347    0.516    1.326)    1.465
   5.512   (   1.960    1.119    0.395)    2.291
   5.576   (   3.447   -2.589   -2.635)    5.052
   5.601   (   1.903   -3.366   -2.328)    4.514
   5.711   (  -2.213   -1.896   -6.355)    6.992
   6.005   (   1.507   -2.217    0.200)    2.689
   6.112   (  -1.327    7.659    5.224)    9.366
   6.162   (   2.083    4.665   -2.912)    5.880
   6.302   (  -2.623   -2.505    1.320)    3.860
   6.452   (   2.289   -1.553    2.414)    3.671
   7.455   (  -0.018    0.969    2.124)    2.335
   7.610   (   0.045    3.319   -0.362)    3.340
   7.645   (   0.955    0.744    1.252)    1.742
   7.705   (   2.326    2.026   -2.120)    3.743
   7.823   (  -4.092   -5.282    1.327)    6.812
   7.904   (  -0.490   -0.955   -0.350)    1.128
   8.047   (   0.616    0.861   -0.568)    1.201
   8.138   (  -1.780   -3.716   -1.819)    4.504
======================= Grid point 19 (9/35) =======================
q-point: ( 0.25  0.25  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 78
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.979   (   2.733    4.734    3.801)    6.658
   1.047   (   2.055    3.560   -0.474)    4.138
   1.114   (   4.570    7.916    3.064)    9.640
   1.322   (   0.790    1.369   -2.893)    3.297
   1.377   (   1.998    3.462    2.949)    4.967
   1.476   (   0.388    0.672    0.086)    0.780
   1.607   (   0.264    0.458   -4.608)    4.638
   1.823   (   0.210    0.363    2.640)    2.673
   1.959   (  -0.741   -1.284    3.762)    4.043
   2.069   (  -5.250   -9.092    0.148)   10.500
   2.105   (  -2.825   -4.893   -3.543)    6.669
   2.190   (  -4.269   -7.394   -3.412)    9.194
   2.343   (   1.450    2.511   -2.341)    3.726
   2.503   (  -1.322   -2.290    4.171)    4.938
   2.524   (  -3.788   -6.562   -2.884)    8.107
   2.689   (  -2.044   -3.541    4.384)    5.995
   2.987   (  -1.041   -1.803   -2.585)    3.319
   3.102   (  -1.672   -2.896   -0.717)    3.420
   3.276   (  -0.950   -1.646   -6.753)    7.015
   3.310   (   0.823    1.426    7.237)    7.422
   3.434   (   6.191   10.722    0.973)   12.419
   3.565   (   1.199    2.077    2.245)    3.285
   3.607   (   0.052    0.090   -0.066)    0.123
   3.697   (   1.542    2.671   -8.135)    8.700
   3.765   (   1.393    2.413   -6.309)    6.897
   3.873   (  -1.447   -2.505    9.252)    9.694
   4.065   (   1.828    3.167    2.518)    4.440
   4.114   (   1.151    1.993    0.143)    2.306
   4.243   (  -2.055   -3.560   -3.767)    5.575
   4.248   (   6.809   11.793   -0.596)   13.631
   4.349   (   1.399    2.423    0.869)    2.930
   4.525   (  -2.269   -3.930    3.500)    5.731
   4.535   (  -2.325   -4.027    1.645)    4.933
   4.717   (   5.131    8.888    5.080)   11.451
   4.851   (  -4.634   -8.027   -0.028)    9.269
   4.878   (   1.234    2.137   -5.378)    5.917
   4.952   (   1.728    2.994    0.215)    3.463
   4.979   (   0.564    0.978   -0.529)    1.247
   5.080   (  -1.299   -2.249   -3.288)    4.190
   5.196   (  -0.629   -1.090    0.551)    1.374
   5.203   (  -3.993   -6.916    8.685)   11.799
   5.259   (  -0.954   -1.653   -2.173)    2.892
   5.458   (  -0.794   -1.375    0.621)    1.706
   5.482   (   0.002    0.003    1.907)    1.907
   5.539   (   1.345    2.330   -1.854)    3.267
   5.584   (   0.481    0.834   -4.317)    4.423
   5.725   (  -1.796   -3.112   -8.167)    8.923
   5.979   (  -0.433   -0.749   -0.593)    1.049
   6.144   (   2.829    4.900   -2.525)    6.196
   6.171   (   3.001    5.197    5.672)    8.258
   6.362   (  -0.759   -1.315    8.881)    9.010
   6.383   (  -1.099   -1.903   -3.893)    4.470
   7.462   (   0.566    0.980    3.341)    3.527
   7.625   (   1.676    2.903    2.273)    4.051
   7.652   (   1.588    2.751   -3.720)    4.892
   7.685   (   1.502    2.601   -0.198)    3.010
   7.835   (  -4.326   -7.493    2.634)    9.044
   7.902   (  -0.799   -1.384   -1.471)    2.172
   8.047   (   0.455    0.788   -0.517)    1.046
   8.133   (  -2.738   -4.743   -2.999)    6.244
======================= Grid point 20 (10/35) =======================
q-point: ( 0.38  0.25  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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 76
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.027   (  -1.485    2.573    3.682)    4.731
   1.060   (   0.020   -0.035    2.107)    2.108
   1.286   (  -3.684    6.382   -0.616)    7.395
   1.330   (  -3.297    5.711    0.647)    6.627
   1.408   (   1.044   -1.808    0.851)    2.255
   1.459   (   0.800   -1.386   -0.009)    1.600
   1.602   (   0.808   -1.399   -3.649)    3.990
   1.864   (   0.016   -0.028    2.077)    2.077
   1.904   (   0.349   -0.605    4.219)    4.276
   1.961   (   2.944   -5.099    1.647)    6.114
   2.049   (  -0.063    0.110   -5.886)    5.887
   2.095   (   2.521   -4.366   -2.099)    5.461
   2.396   (  -4.131    7.155    1.004)    8.322
   2.449   (   3.101   -5.371   -5.546)    8.320
   2.482   (   0.881   -1.526    1.982)    2.651
   2.637   (   4.895   -8.479    3.382)   10.358
   2.969   (  -1.365    2.364   -2.257)    3.543
   3.055   (   1.999   -3.462    1.151)    4.160
   3.257   (   0.382   -0.662   -8.454)    8.489
   3.330   (  -0.023    0.040    8.098)    8.098
   3.592   (  -0.569    0.986    0.643)    1.307
   3.598   (   0.108   -0.186    4.519)    4.524
   3.632   (  -7.829   13.560   -6.535)   16.967
   3.729   (  -1.797    3.113    5.746)    6.778
   3.775   (  -0.852    1.475    3.672)    4.048
   3.845   (   2.044   -3.540  -13.233)   13.850
   4.010   (   6.404  -11.092   -1.374)   12.881
   4.108   (  -1.109    1.921    2.067)    3.032
   4.340   (  -2.710    4.693    1.560)    5.639
   4.366   (  -6.418   11.117   -0.633)   12.853
   4.391   (   0.245   -0.424    1.892)    1.955
   4.402   (   0.539   -0.934    4.052)    4.193
   4.588   (   2.679   -4.639   -1.767)    5.641
   4.660   (   1.159   -2.007    0.663)    2.410
   4.871   (  -2.104    3.644    6.658)    7.876
   4.907   (  -1.699    2.943   -4.570)    5.695
   4.982   (   0.930   -1.611   -2.713)    3.290
   4.991   (   0.029   -0.050   -0.054)    0.079
   5.023   (   0.385   -0.667   -0.787)    1.101
   5.099   (  -0.061    0.105    5.687)    5.688
   5.217   (   3.410   -5.906   -3.578)    7.701
   5.227   (   0.154   -0.266    5.364)    5.373
   5.450   (   0.386   -0.668    0.522)    0.931
   5.476   (   0.156   -0.270    1.823)    1.849
   5.570   (  -0.255    0.442   -1.833)    1.903
   5.587   (   0.522   -0.903   -4.650)    4.766
   5.687   (   0.183   -0.317  -10.325)   10.332
   5.971   (   0.504   -0.873   -0.324)    1.059
   6.227   (  -1.689    2.926    7.009)    7.781
   6.254   (  -0.452    0.784   -1.778)    1.995
   6.273   (  -1.918    3.322    0.551)    3.875
   6.378   (   2.864   -4.960    2.489)    6.245
   7.474   (  -0.346    0.600    3.996)    4.055
   7.666   (  -0.738    1.278    2.239)    2.682
   7.692   (  -2.151    3.725   -1.137)    4.449
   7.731   (  -0.075    0.129   -2.684)    2.688
   7.734   (   1.668   -2.889    1.856)    3.817
   7.873   (   0.958   -1.659   -2.463)    3.120
   8.060   (  -0.111    0.191   -0.970)    0.995
   8.073   (   1.023   -1.771   -1.181)    2.361
======================= Grid point 73 (11/35) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 30
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.354   (   0.000    0.000   17.259)   17.259
   0.354   (  -0.000    0.000   17.259)   17.259
   0.684   (   0.000    0.000   30.894)   30.894
   0.999   (   0.000    0.000  -13.322)   13.322
   0.999   (   0.000    0.000  -13.322)   13.322
   1.259   (  -0.000   -0.000    0.598)    0.598
   1.796   (   0.000    0.000   11.675)   11.675
   1.796   (   0.000   -0.000   11.675)   11.675
   1.991   (   0.000    0.000   -0.025)    0.025
   2.265   (   0.000    0.000   -8.709)    8.709
   2.265   (   0.000    0.000   -8.709)    8.709
   2.383   (  -0.000   -0.000  -15.861)   15.861
   2.853   (   0.000    0.000    3.816)    3.816
   3.044   (   0.000    0.000    3.516)    3.516
   3.044   (   0.000   -0.000    3.516)    3.516
   3.124   (  -0.000   -0.000   -8.641)    8.641
   3.215   (   0.000   -0.000   -3.197)    3.197
   3.215   (  -0.000   -0.000   -3.197)    3.197
   3.351   (   0.000    0.000    4.191)    4.191
   3.351   (   0.000   -0.000    4.191)    4.191
   3.419   (   0.000    0.000    3.527)    3.527
   3.429   (   0.000    0.000    3.140)    3.140
   3.429   (   0.000   -0.000    3.140)    3.140
   3.496   (  -0.000    0.000   -0.533)    0.533
   3.576   (   0.000    0.000    0.070)    0.070
   3.576   (   0.000   -0.000    0.070)    0.070
   3.765   (   0.000    0.000   -5.117)    5.117
   3.765   (   0.000    0.000   -5.117)    5.117
   4.137   (   0.000    0.000    1.491)    1.491
   4.210   (  -0.000   -0.000   -1.422)    1.422
   4.232   (   0.000    0.000   15.217)   15.217
   4.783   (   0.000    0.000   -0.038)    0.038
   4.783   (  -0.000    0.000   -0.038)    0.038
   4.816   (   0.000    0.000   -0.916)    0.916
   4.816   (   0.000   -0.000   -0.916)    0.916
   4.900   (   0.000    0.000    0.226)    0.226
   4.900   (  -0.000   -0.000    0.226)    0.226
   4.903   (  -0.000   -0.000  -14.894)   14.894
   4.921   (   0.000    0.000   -1.057)    1.057
   4.921   (   0.000   -0.000   -1.057)    1.057
   5.427   (   0.000    0.000    9.531)    9.531
   5.507   (   0.000    0.000    3.375)    3.375
   5.507   (  -0.000    0.000    3.375)    3.375
   5.643   (   0.000    0.000   -1.738)    1.738
   5.643   (   0.000    0.000   -1.738)    1.738
   5.737   (   0.000    0.000   -5.497)    5.497
   5.936   (   0.000   -0.000    0.860)    0.860
   5.936   (   0.000    0.000    0.860)    0.860
   5.964   (  -0.000   -0.000    2.266)    2.266
   5.988   (   0.000    0.000   -0.921)    0.921
   5.988   (  -0.000   -0.000   -0.921)    0.921
   6.022   (   0.000    0.000   -0.590)    0.590
   7.417   (   0.000    0.000    0.068)    0.068
   7.417   (   0.000   -0.000    0.068)    0.068
   7.421   (   0.000    0.000   -0.066)    0.066
   7.421   (   0.000   -0.000   -0.066)    0.066
   7.928   (   0.000    0.000    0.805)    0.805
   7.969   (  -0.000   -0.000   -0.773)    0.773
   8.523   (   0.000    0.000    0.159)    0.159
   8.531   (  -0.000    0.000   -0.157)    0.157
======================= Grid point 74 (12/35) =======================
q-point: ( 0.12  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.478   (   7.992    8.607   11.096)   16.158
   0.493   (   7.717    6.182   12.273)   15.760
   0.943   (  18.113    7.863   17.858)   26.624
   1.021   (   2.200    1.142  -13.017)   13.250
   1.051   (   3.408   -0.124  -11.667)   12.156
   1.271   (   0.909    2.493    1.700)    3.152
   1.698   (  -3.293   -6.695   10.334)   12.746
   1.809   (  -1.975    3.063    8.255)    9.024
   1.926   (  -4.206    1.020    0.923)    4.426
   2.238   (   1.109   -4.395   -9.931)   10.917
   2.315   (   1.148    2.468   -5.347)    5.999
   2.475   (   5.762    2.161  -11.089)   12.682
   2.750   (  -6.408   -5.111    5.370)    9.799
   2.937   (  -8.386   -4.499   -0.467)    9.528
   2.993   ( -11.897   -8.193    0.398)   14.450
   3.033   ( -14.227   -4.415    0.300)   14.899
   3.118   (   1.139    3.797    1.522)    4.246
   3.132   (  -0.206    0.073    1.822)    1.835
   3.244   (  -2.524    4.733   -0.648)    5.403
   3.331   (   3.384   -4.773   -3.338)    6.736
   3.390   (   1.427    1.668   -0.557)    2.264
   3.453   (   2.764   -8.975    3.727)   10.103
   3.477   (  -1.700   -0.850   -2.527)    3.162
   3.518   (   2.016    1.598    0.118)    2.575
   3.592   (   2.705   -1.301    0.307)    3.017
   3.675   (  -0.089    5.782   -0.932)    5.857
   3.809   (   5.770    9.708    3.310)   11.768
   3.924   (   6.546   -1.013    2.874)    7.221
   4.151   (  -0.400   -0.003    4.702)    4.719
   4.196   (  -0.188    0.052    7.660)    7.663
   4.220   (   1.766    0.502   -0.084)    1.838
   4.640   (  -9.399   -2.530   -6.600)   11.760
   4.758   (  -2.547   -1.057   -0.043)    2.758
   4.769   (  -1.510   -0.761    1.043)    1.986
   4.779   (  -2.021   -0.266   -1.346)    2.443
   4.862   (   2.678    0.952   -3.395)    4.428
   4.935   (   0.768    0.202    1.770)    1.940
   4.955   (   2.598    3.580    0.647)    4.470
   5.007   (   7.136   -1.608    1.790)    7.531
   5.121   (   7.369    4.220   -4.320)    9.527
   5.300   (  -4.183   -1.073    7.951)    9.048
   5.444   (  -2.948   -0.841    4.292)    5.274
   5.507   (   0.811   -3.851    2.961)    4.924
   5.570   (  -4.739   -0.470   -0.686)    4.811
   5.622   (  -6.972   -0.634   -1.384)    7.136
   5.644   (  -3.292    0.389   -2.109)    3.929
   5.775   (   2.006   -4.156   -1.964)    5.016
   5.845   (   0.121   -2.972   -5.015)    5.830
   5.947   (   1.599   -0.873    1.160)    2.160
   5.988   (   1.676    1.904    0.688)    2.628
   6.204   (   9.892    2.654    4.282)   11.101
   6.380   (   1.529   11.780   -3.893)   12.500
   7.414   (  -0.165    0.914    0.415)    1.017
   7.444   (   2.101    0.154   -0.611)    2.194
   7.462   (   2.284    2.989    0.447)    3.788
   7.556   (   6.214    2.301   -6.285)    9.133
   7.943   (   0.805    0.510    1.085)    1.444
   7.978   (   0.704    0.397   -0.582)    0.996
   8.403   (  -5.577   -6.621    0.106)    8.658
   8.422   (  -8.048   -4.064   -0.438)    9.026
======================= Grid point 75 (13/35) =======================
q-point: ( 0.25  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.692   (   6.854    7.195    8.167)   12.862
   0.753   (   9.908    7.698    3.885)   13.135
   1.048   (   0.202   -0.273   -7.699)    7.707
   1.110   (   3.904   -0.284   -8.443)    9.306
   1.322   (   9.080    4.609    4.308)   11.057
   1.345   (   5.423    4.525    1.952)    7.327
   1.627   (   2.698   -3.196    9.178)   10.086
   1.767   (  -3.963    0.233   -3.597)    5.358
   1.942   (   1.026    7.374    7.024)   10.235
   2.199   (  -0.464   -5.450   -9.871)   11.285
   2.296   (  -3.641   -4.997    0.827)    6.237
   2.519   (  -4.626   -4.342    2.754)    6.916
   2.577   (  -4.698   -2.143   -0.650)    5.204
   2.679   ( -11.563   -4.455   -0.650)   12.408
   2.736   (  -6.533   -0.862    1.014)    6.667
   2.843   (  -2.899   -3.448   -1.722)    4.823
   3.121   (  -0.573   -2.402    1.055)    2.685
   3.164   (  -0.529    1.263    1.212)    1.829
   3.217   (  -6.729    6.748    2.252)    9.792
   3.337   (   3.178   -8.586   -1.108)    9.222
   3.400   (   1.277   -4.108    4.465)    6.200
   3.409   (  -1.176    2.746   -5.323)    6.104
   3.465   (   1.562    0.068   -4.550)    4.811
   3.580   (   1.861    3.703   -2.450)    4.814
   3.684   (   9.632   -2.916   -0.794)   10.095
   3.736   (   1.717    4.778   -3.866)    6.382
   3.967   (   3.832   -2.479    7.568)    8.837
   4.004   (   1.084    8.124    5.659)    9.960
   4.097   (  -1.484   -0.687    5.910)    6.133
   4.262   (   2.232    1.066   -0.866)    2.621
   4.325   (   5.294    6.409    2.544)    8.693
   4.475   (  -4.568   -2.910   -7.025)    8.871
   4.678   (  -4.114   -2.066    0.679)    4.653
   4.711   (  -3.311   -1.093   -1.531)    3.808
   4.767   (  -1.476    1.003    1.594)    2.393
   4.911   (  -0.858   -0.865   -0.927)    1.532
   4.941   (   2.304    1.354    2.063)    3.375
   5.042   (   2.280    1.763    0.146)    2.886
   5.072   (   3.385   -4.133    3.917)    6.625
   5.200   (   0.883   -2.261   -3.072)    3.915
   5.336   (  -0.366    3.844    0.220)    3.868
   5.404   (  -0.515   -1.668   -1.835)    2.533
   5.443   (  -2.375    0.928    0.822)    2.679
   5.479   (  -3.131   -1.349    0.520)    3.448
   5.508   (  -1.109   -0.071    2.066)    2.346
   5.644   (   0.469    0.114   -3.144)    3.181
   5.801   (   2.541   -3.129   -1.375)    4.259
   5.846   (   1.389    1.316   -3.749)    4.209
   5.984   (   2.730    0.112    1.065)    2.932
   6.049   (   3.222    1.298    0.032)    3.474
   6.353   (   3.829    0.022    3.648)    5.288
   6.447   (  -2.285    7.845    2.310)    8.492
   7.429   (   0.260    1.539    0.842)    1.773
   7.492   (   2.619    0.684   -1.367)    3.032
   7.549   (   2.935    3.003    0.587)    4.240
   7.629   (   1.729    0.592   -2.026)    2.729
   7.951   (  -0.245   -0.078    2.278)    2.293
   8.000   (   1.070    0.501   -0.312)    1.222
   8.184   ( -11.474   -4.889   -0.625)   12.488
   8.258   (  -2.963   -6.736   -4.549)    8.652
======================= Grid point 76 (14/35) =======================
q-point: ( 0.38  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.859   (   3.884    4.959    3.603)    7.257
   0.952   (   4.335    2.141    0.309)    4.845
   1.062   (   2.204    2.434   -1.574)    3.641
   1.168   (   2.018    0.745   -3.954)    4.502
   1.418   (   3.936    0.558    0.217)    3.981
   1.484   (   1.293    2.772    3.881)    4.942
   1.669   (   1.021   -0.874    3.318)    3.580
   1.738   (   2.213    1.443   -5.478)    6.081
   1.985   (  -5.907    0.898   -2.131)    6.344
   2.074   (  -6.514   -5.610   -7.907)   11.680
   2.173   (   0.443   -4.780    9.108)   10.296
   2.368   (  -3.382   -3.034    7.781)    9.010
   2.460   (  -4.607    1.584   -1.822)    5.201
   2.535   (  -3.335   -1.958    1.600)    4.185
   2.711   (  -0.625    2.630    0.353)    2.727
   2.791   (   0.474   -4.686   -3.208)    5.698
   3.083   (  -0.086   -2.050    3.029)    3.659
   3.121   (  -1.508   -1.557    0.759)    2.296
   3.180   (  -4.313    4.341   -0.272)    6.125
   3.298   (  -0.039   -3.123    1.752)    3.581
   3.400   (   3.833   -6.822   -1.472)    7.962
   3.452   (   2.607    3.699    2.390)    5.117
   3.512   (   0.090    5.072   -8.468)    9.871
   3.652   (   1.159    3.891   -1.991)    4.522
   3.791   (  -1.384    4.216   -0.092)    4.438
   3.843   (   9.198   -3.729   -0.163)    9.927
   3.956   (  -0.804   -2.806    4.544)    5.401
   4.082   (   3.925    2.488    2.967)    5.513
   4.156   (   1.049    4.842    0.053)    4.955
   4.249   (  -0.878   -0.495   -4.754)    4.860
   4.428   (  -0.027    2.856    3.179)    4.273
   4.451   (   2.049   -3.086    0.026)    3.704
   4.610   (  -4.280   -1.300   -2.686)    5.218
   4.649   (   0.492    3.043   -0.049)    3.082
   4.730   (  -3.269    0.685    2.800)    4.358
   4.866   (  -0.393   -4.292    1.397)    4.531
   4.981   (   0.352    0.148   -1.768)    1.809
   5.035   (  -2.751   -1.303    1.124)    3.244
   5.076   (   0.936   -1.197    0.326)    1.554
   5.162   (  -2.639    0.368   -4.464)    5.199
   5.301   (  -4.496    2.151    2.266)    5.475
   5.349   (   0.655   -3.094    1.032)    3.326
   5.471   (   1.508    2.528    1.383)    3.253
   5.490   (  -0.155    1.423    2.186)    2.613
   5.548   (   4.340   -0.475   -1.139)    4.512
   5.657   (   2.656   -2.931   -4.149)    5.732
   5.781   (   1.548   -6.140    2.018)    6.645
   5.903   (  -0.993    6.670   -4.149)    7.918
   6.038   (   2.369    0.508    1.085)    2.655
   6.133   (   4.749    0.982   -0.826)    4.919
   6.340   (  -4.023    0.587    4.078)    5.759
   6.499   (   0.812    2.819    1.821)    3.453
   7.450   (   0.050    1.698    1.207)    2.083
   7.544   (   2.019    0.673   -1.903)    2.855
   7.593   (   0.223    0.563   -0.588)    0.845
   7.687   (   3.435    1.363   -0.887)    3.801
   7.924   (  -4.053   -2.683    2.046)    5.274
   7.961   (  -5.588   -1.040    1.519)    5.884
   8.024   (   1.128    0.261   -0.432)    1.235
   8.152   (  -2.081   -4.545   -4.798)    6.929
======================= Grid point 77 (15/35) =======================
q-point: (-0.50  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.918   (  -1.062    1.840    1.297)    2.489
   0.984   (   0.519   -0.899   -1.870)    2.139
   1.154   (  -0.390    0.675   -1.128)    1.371
   1.166   (  -0.282    0.488   -1.618)    1.714
   1.463   (   0.531   -0.919    5.550)    5.651
   1.487   (  -0.867    1.501    3.699)    4.085
   1.675   (   0.555   -0.962   -0.405)    1.182
   1.853   (  -1.273    2.204   -7.371)    7.798
   1.856   (  -0.132    0.228   -2.890)    2.902
   1.935   (   0.352   -0.609   -6.769)    6.805
   2.155   (   2.573   -4.456    8.969)   10.340
   2.330   (  -0.614    1.064   10.171)   10.245
   2.453   (  -1.453    2.517   -1.146)    3.124
   2.491   (   0.001   -0.001   -0.365)    0.365
   2.720   (  -1.316    2.279   -3.288)    4.211
   2.762   (   2.979   -5.160   -0.342)    5.968
   3.062   (   0.354   -0.613    6.790)    6.827
   3.131   (  -3.531    6.115   -1.590)    7.238
   3.146   (   3.918   -6.786   -0.897)    7.887
   3.277   (  -3.100    5.369    1.787)    6.452
   3.399   (   6.093  -10.554   -2.995)   12.549
   3.529   (  -1.649    2.857    5.250)    6.200
   3.537   (  -2.827    4.897  -11.795)   13.080
   3.695   (  -0.737    1.277   -1.614)    2.186
   3.758   (  -2.909    5.038    5.548)    8.039
   3.928   (   6.025  -10.436    5.051)   13.066
   3.935   (  -0.433    0.751   -2.684)    2.820
   4.114   (  -1.817    3.147   -3.647)    5.149
   4.185   (   0.037   -0.064    3.889)    3.890
   4.274   (   1.547   -2.679   -1.218)    3.324
   4.449   (   1.537   -2.663   -3.234)    4.462
   4.450   (  -2.350    4.070   -4.079)    6.222
   4.562   (  -2.540    4.399    3.570)    6.208
   4.610   (   1.687   -2.921    2.600)    4.259
   4.786   (  -1.864    3.228    2.518)    4.498
   4.801   (   1.141   -1.977   -1.524)    2.745
   4.980   (   0.484   -0.839   -3.249)    3.391
   4.981   (   0.504   -0.873    1.417)    1.739
   5.091   (   0.075   -0.130   -0.650)    0.667
   5.126   (  -2.017    3.493   -4.097)    5.749
   5.274   (  -1.709    2.959    3.806)    5.115
   5.346   (   1.621   -2.808    1.757)    3.688
   5.483   (  -0.723    1.252    2.346)    2.756
   5.521   (  -0.885    1.532    0.699)    1.902
   5.581   (   1.771   -3.067   -0.778)    3.626
   5.675   (   2.688   -4.655   -4.675)    7.124
   5.754   (   3.306   -5.725    4.555)    8.028
   5.931   (  -3.859    6.684   -4.429)    8.899
   6.064   (   0.463   -0.803    0.981)    1.350
   6.218   (  -0.198    0.342   -6.443)    6.455
   6.270   (   0.007   -0.011    9.771)    9.771
   6.521   (  -1.126    1.951    0.703)    2.359
   7.460   (  -0.734    1.271    1.376)    2.011
   7.568   (   0.218   -0.378   -2.000)    2.047
   7.597   (  -0.192    0.332   -1.168)    1.230
   7.740   (  -0.101    0.174   -5.528)    5.532
   7.826   (   0.285   -0.494    6.624)    6.649
   7.933   (  -0.277    0.481    1.695)    1.784
   8.040   (   0.325   -0.562   -0.791)    1.023
   8.103   (   0.720   -1.247   -3.510)    3.794
======================= Grid point 83 (16/35) =======================
q-point: ( 0.12  0.12  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.695   (   6.370   11.033    5.401)   13.837
   0.696   (   6.218   10.771    9.070)   15.393
   1.043   (  -0.058   -0.100   -9.095)    9.096
   1.074   (   2.299    3.982  -11.165)   12.074
   1.211   (   8.190   14.186    8.781)   18.586
   1.338   (   1.404    2.431    0.179)    2.813
   1.585   (  -1.371   -2.374   10.248)   10.608
   1.818   (  -0.930   -1.611   -0.047)    1.860
   1.987   (   2.809    4.865    9.237)   10.811
   2.143   (  -2.481   -4.297  -11.802)   12.803
   2.319   (  -2.069   -3.583   -2.451)    4.809
   2.524   (   0.797    1.380   -5.609)    5.831
   2.572   (  -5.914  -10.243    4.625)   12.700
   2.797   (  -5.382   -9.321    0.936)   10.804
   2.807   (  -2.868   -4.968   -0.387)    5.749
   2.888   (  -6.516  -11.286    0.329)   13.036
   3.123   (  -1.908   -3.305    3.440)    5.138
   3.130   (  -1.101   -1.907    2.860)    3.609
   3.214   (   0.762    1.320   -4.965)    5.194
   3.335   (   0.857    1.484    0.705)    1.853
   3.407   (   0.304    0.526   -4.586)    4.626
   3.422   (   0.755    1.307   -2.485)    2.908
   3.499   (  -1.853   -3.210    1.463)    3.985
   3.509   (   2.474    4.285   -6.425)    8.109
   3.610   (   2.069    3.584    4.383)    6.028
   3.813   (   1.777    3.078    2.660)    4.439
   3.835   (   2.825    4.893    1.666)    5.890
   4.034   (   2.274    3.938    4.836)    6.638
   4.110   (  -1.417   -2.454    7.009)    7.560
   4.239   (   1.109    1.920   -0.615)    2.301
   4.319   (   6.351   11.000    3.746)   13.242
   4.505   (  -5.079   -8.797   -7.164)   12.430
   4.744   (  -0.992   -1.718   -3.098)    3.679
   4.748   (  -1.023   -1.773   -0.952)    2.257
   4.819   (  -0.399   -0.691    2.822)    2.933
   4.861   (   1.501    2.600   -2.518)    3.918
   4.918   (  -0.552   -0.956    4.300)    4.439
   4.980   (   1.036    1.795    4.515)    4.968
   5.096   (   2.545    4.409   -2.411)    5.633
   5.128   (  -0.903   -1.564   -0.183)    1.815
   5.354   (   2.807    4.863    0.242)    5.620
   5.408   (  -2.235   -3.872   -2.179)    4.974
   5.468   (   0.062    0.108    2.363)    2.366
   5.489   (  -3.261   -5.648    0.375)    6.533
   5.554   (  -0.906   -1.570   -0.955)    2.049
   5.627   (  -1.951   -3.379    2.573)    4.673
   5.765   (  -2.053   -3.556   -6.264)    7.490
   5.821   (   3.548    6.145   -1.849)    7.333
   5.958   (   0.969    1.678   -0.345)    1.968
   6.039   (   1.690    2.926    2.043)    3.949
   6.283   (   2.229    3.860    1.228)    4.624
   6.548   (   1.873    3.244    2.909)    4.743
   7.454   (   1.188    2.058   -1.010)    2.582
   7.455   (   0.974    1.687    1.736)    2.609
   7.545   (   2.424    4.198    0.881)    4.927
   7.608   (   0.990    1.715   -4.866)    5.253
   7.953   (  -0.000   -0.001    2.330)    2.330
   7.993   (   0.525    0.910   -0.408)    1.127
   8.234   (  -4.311   -7.467   -4.612)    9.779
   8.271   (  -5.516   -9.554    0.277)   11.035
======================= Grid point 84 (17/35) =======================
q-point: ( 0.25  0.12  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.885   (   2.920    9.581    4.679)   11.055
   0.918   (   5.277    7.163    1.883)    9.094
   1.058   (   1.222    2.738   -2.513)    3.912
   1.164   (   2.017    4.256   -7.131)    8.546
   1.371   (   1.894    1.287    1.008)    2.502
   1.443   (   5.917    3.921    2.174)    7.424
   1.606   (   3.076    0.237    6.214)    6.938
   1.771   (  -2.619    0.294   -3.514)    4.392
   2.046   (  -1.844   -3.695   -7.591)    8.641
   2.115   (   2.471    1.874    7.532)    8.145
   2.146   (  -7.949   -7.332   -3.817)   11.468
   2.391   (  -4.001   -4.851    6.604)    9.119
   2.492   (  -1.961  -10.182    3.990)   11.111
   2.589   (  -6.179   -1.612   -2.896)    7.012
   2.675   (  -3.493   -8.253   -0.474)    8.974
   2.790   (   0.771    0.652    0.287)    1.049
   3.062   (  -1.676   -3.635    3.674)    5.433
   3.132   (   4.346   -4.251    2.772)    6.681
   3.161   (  -3.888   -1.023   -7.545)    8.549
   3.343   (  -3.563    4.198    3.242)    6.390
   3.383   (  -1.608   -1.623   -1.590)    2.784
   3.454   (   1.368    3.568   -0.950)    3.938
   3.514   (   5.104   -0.526   -5.015)    7.175
   3.639   (   5.004    2.441   -3.271)    6.457
   3.689   (   2.310    3.905    3.065)    5.475
   3.872   (  -0.341    8.990   -0.078)    8.997
   3.985   (   2.005    4.437    0.969)    4.964
   4.071   (  -0.342    0.087    5.261)    5.273
   4.114   (   5.508    3.054    2.621)    6.822
   4.280   (  -1.503   -0.948   -3.836)    4.228
   4.346   (   1.192   -5.855   -2.547)    6.496
   4.533   (  -1.514   10.436    4.840)   11.603
   4.651   (  -3.191    0.586   -0.393)    3.268
   4.716   (  -1.531   -2.338   -0.922)    2.943
   4.785   (   0.263    0.010    1.274)    1.301
   4.853   (  -3.285   -1.543   -1.703)    4.009
   4.958   (   2.593    0.029    0.326)    2.614
   5.021   (   1.751   -1.236    3.348)    3.975
   5.062   (  -2.635   -1.632   -0.201)    3.106
   5.180   (   1.974    1.888   -2.432)    3.658
   5.290   (  -1.606   -6.837   -2.158)    7.348
   5.405   (  -3.158    1.347    2.981)    4.547
   5.466   (   1.467    1.402    0.848)    2.199
   5.485   (   1.760    0.529   -1.895)    2.640
   5.517   (  -1.168    0.328    2.109)    2.433
   5.574   (   1.040   -4.296   -0.914)    4.513
   5.699   (  -0.293   -5.117   -2.430)    5.672
   5.940   (   3.138    3.225   -2.439)    5.118
   6.012   (   2.151    2.862   -0.914)    3.695
   6.115   (   1.483    5.895    2.459)    6.557
   6.341   (   1.954   -0.862   -0.923)    2.327
   6.555   (  -2.522    1.396    5.460)    6.174
   7.482   (  -0.326    3.348    1.629)    3.738
   7.521   (   2.621    2.072   -1.074)    3.510
   7.605   (  -0.343    1.849   -0.947)    2.106
   7.641   (   2.719    0.691   -0.765)    2.908
   7.933   (  -1.117   -2.070    1.902)    3.025
   8.006   (  -0.283   -0.763   -0.608)    1.016
   8.064   (  -6.352   -5.189    2.229)    8.500
   8.120   (  -1.099   -5.283   -6.834)    8.708
======================= Grid point 85 (18/35) =======================
q-point: ( 0.38  0.12  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.986   (  -0.727    5.700    1.162)    5.862
   0.994   (  -0.446    2.787   -1.911)    3.409
   1.163   (   2.061    6.623    1.074)    7.019
   1.216   (  -1.252    3.535   -2.917)    4.751
   1.436   (   2.852    1.819    2.490)    4.200
   1.510   (  -0.427   -0.571    3.868)    3.933
   1.643   (   1.400   -2.181    1.523)    3.007
   1.791   (   4.238    3.436   -5.159)    7.509
   1.947   (  -3.167   -4.381   -6.555)    8.497
   1.962   (  -6.080   -2.743   -4.151)    7.856
   2.083   (  -1.608   -4.024    6.436)    7.759
   2.359   (   3.325   -7.078    6.386)   10.096
   2.365   (  -0.897   -1.475    5.325)    5.598
   2.489   (  -2.288   -2.732   -0.540)    3.605
   2.617   (   4.117   -7.479   -0.843)    8.579
   2.786   (  -1.463    0.481   -0.759)    1.717
   3.020   (   2.394   -5.648    6.107)    8.656
   3.082   (  -1.950   -0.135   -8.182)    8.412
   3.181   (   1.274   -0.439    2.986)    3.276
   3.302   (  -2.592    1.861   -2.706)    4.183
   3.391   (  -1.263    3.957    1.838)    4.542
   3.546   (   4.233    4.699   -1.409)    6.480
   3.597   (   1.669    1.772   -6.139)    6.604
   3.701   (   2.399   -2.823   -0.727)    3.775
   3.763   (   2.470    2.308    4.046)    5.272
   3.879   (  -7.843    0.972    3.075)    8.480
   3.992   (  -1.147    6.064   -2.148)    6.534
   4.152   (   4.055    4.971    2.791)    6.996
   4.190   (  -0.691   -3.864   -2.367)    4.584
   4.274   (  -0.411    5.629    0.510)    5.667
   4.392   (   1.851    0.338   -0.478)    1.941
   4.480   (   0.816    1.500    2.805)    3.284
   4.594   (  -3.894   -2.176   -0.629)    4.505
   4.687   (  -6.470    0.169   -1.376)    6.617
   4.786   (   0.005    1.508   -0.646)    1.640
   4.848   (   1.761    3.977   -0.627)    4.394
   4.965   (  -0.420   -0.173   -0.367)    0.584
   4.984   (  -2.423   -2.821    0.312)    3.732
   5.050   (   3.308   -1.255    0.286)    3.550
   5.199   (  -2.865    2.104   -3.244)    4.812
   5.248   (   3.234   -5.776   -1.542)    6.797
   5.334   (  -2.597    0.436    6.364)    6.887
   5.499   (   0.615   -0.842   -0.375)    1.108
   5.523   (   2.778   -0.637    0.101)    2.852
   5.532   (   2.272   -2.358   -0.398)    3.298
   5.567   (   0.445    0.079    0.859)    0.971
   5.643   (   0.846   -5.927    0.489)    6.006
   5.996   (   2.036   -1.160   -0.857)    2.495
   6.085   (  -0.237    6.420   -2.064)    6.748
   6.192   (   1.045    2.909   -1.717)    3.536
   6.328   (  -2.379    0.521    4.219)    4.871
   6.525   (   0.162   -1.805    2.054)    2.739
   7.504   (  -1.280    3.435    1.866)    4.113
   7.576   (   0.555    2.220   -2.176)    3.157
   7.607   (  -0.568    0.813   -1.134)    1.506
   7.706   (   3.592    0.562   -2.817)    4.600
   7.845   (  -4.210   -4.080    2.911)    6.546
   7.930   (  -0.400   -2.600    2.007)    3.309
   8.033   (   0.184    0.418   -0.486)    0.667
   8.073   (   0.156   -2.536   -3.654)    4.450
======================= Grid point 92 (19/35) =======================
q-point: ( 0.25  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.014   (   1.464    2.536   -2.801)    4.052
   1.031   (   2.323    4.023    0.410)    4.663
   1.201   (   5.282    9.148    5.281)   11.810
   1.251   (   2.172    3.763   -4.067)    5.951
   1.424   (   2.000    3.464    1.131)    4.157
   1.492   (   0.185    0.321    1.150)    1.208
   1.597   (  -0.763   -1.322    3.206)    3.551
   1.789   (   1.103    1.911   -3.344)    4.006
   1.945   (  -2.440   -4.226   -7.258)    8.746
   1.989   (  -2.877   -4.983   -1.980)    6.085
   2.116   (  -2.395   -4.147    3.789)    6.107
   2.251   (  -5.864  -10.156    5.048)   12.768
   2.388   (   2.027    3.511    6.417)    7.591
   2.486   (  -4.794   -8.304   -0.197)    9.590
   2.551   (  -1.140   -1.975   -1.384)    2.667
   2.749   (  -2.535   -4.391    1.375)    5.253
   2.997   (  -0.778   -1.347    4.633)    4.887
   3.088   (  -0.769   -1.331   -4.432)    4.691
   3.170   (  -0.987   -1.710    0.803)    2.131
   3.395   (   1.372    2.376   -4.870)    5.590
   3.408   (   2.887    5.000   -4.266)    7.178
   3.485   (   0.473    0.820   -0.489)    1.066
   3.558   (   2.579    4.468   -0.429)    5.177
   3.678   (   0.489    0.847    0.137)    0.988
   3.766   (   1.876    3.249    4.052)    5.522
   4.027   (   0.038    0.066    4.412)    4.413
   4.050   (   1.386    2.400   -1.338)    3.078
   4.077   (  -3.162   -5.476   -2.565)    6.824
   4.216   (   3.659    6.337    0.864)    7.368
   4.239   (   5.908   10.233   -0.923)   11.852
   4.370   (   2.177    3.771    0.917)    4.450
   4.582   (  -3.364   -5.827    1.966)    7.010
   4.595   (  -2.599   -4.501    2.696)    5.855
   4.780   (   0.441    0.763   -3.435)    3.546
   4.790   (  -3.377   -5.849   -4.899)    8.343
   4.801   (   4.600    7.967    1.264)    9.286
   4.984   (   0.810    1.403    2.944)    3.361
   4.988   (  -0.607   -1.052    2.104)    2.430
   5.013   (  -0.931   -1.613   -2.630)    3.223
   5.160   (  -2.206   -3.820   -2.876)    5.266
   5.254   (   1.995    3.456   -1.150)    4.152
   5.355   (  -3.559   -6.165    5.736)    9.142
   5.459   (  -1.056   -1.829   -1.665)    2.689
   5.485   (  -0.427   -0.739   -3.203)    3.315
   5.529   (   0.872    1.510    1.606)    2.371
   5.543   (   0.676    1.170    1.833)    2.277
   5.623   (  -1.344   -2.328    1.370)    3.017
   5.963   (   0.118    0.204   -0.868)    0.899
   6.098   (   2.882    4.992   -1.785)    6.034
   6.260   (   3.492    6.048    2.411)    7.388
   6.313   (  -0.836   -1.448   -3.545)    3.920
   6.515   (  -2.524   -4.371    5.379)    7.376
   7.552   (   1.542    2.671    4.434)    5.401
   7.573   (   1.527    2.644   -2.939)    4.238
   7.643   (   1.393    2.413   -1.976)    3.416
   7.650   (   0.073    0.127   -1.939)    1.945
   7.868   (  -2.175   -3.768   -1.644)    4.651
   7.907   (  -3.325   -5.759    3.622)    7.572
   8.037   (   0.163    0.282    0.073)    0.334
   8.056   (  -0.768   -1.331   -3.252)    3.597
======================= Grid point 93 (20/35) =======================
q-point: ( 0.38  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.047   (  -1.052    1.821   -2.932)    3.608
   1.063   (  -0.753    1.304   -0.970)    1.792
   1.297   (  -3.536    6.124    0.533)    7.092
   1.343   (  -3.911    6.774    1.394)    7.945
   1.455   (   2.188   -3.790    3.891)    5.856
   1.476   (  -0.354    0.612   -0.923)    1.163
   1.572   (   2.466   -4.271    2.323)    5.451
   1.842   (  -0.427    0.740   -3.432)    3.537
   1.879   (   0.471   -0.815   -4.569)    4.665
   1.950   (  -0.110    0.190   -3.312)    3.319
   2.036   (  -0.379    0.656    4.454)    4.518
   2.141   (   4.552   -7.885    4.740)   10.264
   2.384   (   2.860   -4.953   -0.421)    5.735
   2.453   (  -3.185    5.516    3.937)    7.488
   2.521   (   0.865   -1.498    0.733)    1.878
   2.681   (   4.848   -8.396    1.427)    9.800
   2.990   (  -1.782    3.086    4.992)    6.133
   3.088   (  -0.067    0.116   -8.255)    8.256
   3.115   (   2.929   -5.073    5.343)    7.928
   3.418   (  -3.914    6.780  -10.412)   13.027
   3.453   (  -1.087    1.883    2.012)    2.963
   3.543   (   2.514   -4.354  -10.346)   11.503
   3.657   (  -1.430    2.476    2.958)    4.114
   3.667   (  -0.276    0.478   -0.104)    0.562
   3.831   (  -0.194    0.336    1.539)    1.587
   3.893   (  -2.780    4.814    8.035)    9.770
   4.043   (   4.217   -7.304    4.605)    9.609
   4.089   (  -1.368    2.369   -3.448)    4.401
   4.355   (  -4.881    8.453   -0.909)    9.803
   4.405   (  -3.264    5.653    3.136)    7.242
   4.420   (  -1.556    2.695    0.663)    3.181
   4.479   (   0.947   -1.640    2.842)    3.415
   4.542   (   1.607   -2.784   -1.548)    3.568
   4.655   (   1.719   -2.977   -1.141)    3.622
   4.801   (  -0.419    0.726   -5.091)    5.159
   4.860   (   1.535   -2.659   -6.196)    6.915
   4.973   (  -1.247    2.160   -0.182)    2.501
   5.012   (  -1.101    1.906    1.970)    2.954
   5.028   (   0.760   -1.316    3.222)    3.563
   5.153   (   0.936   -1.622   -2.540)    3.155
   5.216   (   0.215   -0.373    2.251)    2.292
   5.311   (   1.171   -2.029    0.483)    2.392
   5.442   (   0.931   -1.613   -0.817)    2.034
   5.466   (   1.031   -1.785   -3.157)    3.770
   5.555   (   0.701   -1.215    1.633)    2.153
   5.563   (  -0.791    1.370    3.064)    3.448
   5.591   (   0.099   -0.172    3.513)    3.518
   5.964   (   0.558   -0.967   -0.310)    1.159
   6.188   (   0.046   -0.079   -3.855)    3.856
   6.248   (  -4.051    7.016   -1.186)    8.188
   6.375   (  -0.788    1.364    5.009)    5.250
   6.427   (   3.513   -6.085    0.780)    7.070
   7.580   (  -1.767    3.061    3.932)    5.288
   7.618   (  -0.283    0.490   -1.055)    1.197
   7.626   (  -0.734    1.272   -3.717)    3.997
   7.717   (  -0.095    0.165   -6.363)    6.366
   7.784   (   0.211   -0.366    4.144)    4.165
   7.849   (   2.592   -4.489    1.545)    5.409
   8.039   (  -0.133    0.231   -0.646)    0.699
   8.046   (   0.279   -0.484   -0.974)    1.122
======================= Grid point 98 (21/35) =======================
q-point: (-0.12  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.631   (  -6.233   10.795    7.329)   14.461
   0.641   (  -5.362    9.288    6.521)   12.552
   1.045   (  -0.353    0.611  -10.356)   10.380
   1.152   (  -3.433    5.945   -9.834)   11.994
   1.256   (  -7.857   13.609    9.211)   18.215
   1.269   (  -0.418    0.724    2.808)    2.929
   1.713   (   0.217   -0.377   11.277)   11.286
   1.740   (   0.453   -0.785    0.280)    0.948
   1.844   (   1.045   -1.810    0.888)    2.270
   2.265   (   0.829   -1.436   -9.489)    9.633
   2.340   (  -0.973    1.686   -1.141)    2.256
   2.530   (   1.138   -1.972   -3.283)    3.995
   2.644   (   3.859   -6.684    6.775)   10.270
   2.724   (   6.727  -11.652   -4.032)   14.046
   2.792   (   4.114   -7.126    0.263)    8.233
   2.841   (   6.304  -10.918    1.896)   12.749
   3.099   (   1.532   -2.654    1.098)    3.255
   3.138   (  -1.419    2.459    5.903)    6.550
   3.186   (   0.016   -0.027   -3.529)    3.529
   3.332   (   0.374   -0.648   -4.164)    4.231
   3.378   (   2.638   -4.569   -1.575)    5.506
   3.494   (  -2.112    3.657   -0.924)    4.323
   3.520   (  -1.166    2.020   -4.035)    4.661
   3.639   (  -1.311    2.271   -5.336)    5.946
   3.704   (  -2.448    4.240    7.115)    8.637
   3.723   (  -4.021    6.964    2.376)    8.385
   3.771   (  -1.469    2.545   -2.661)    3.964
   4.020   (  -1.091    1.889    9.028)    9.288
   4.104   (   1.616   -2.798    8.368)    8.970
   4.251   (  -3.051    5.285    0.535)    6.126
   4.264   (  -1.414    2.449    0.954)    2.984
   4.490   (   3.577   -6.195   -7.041)   10.037
   4.715   (   2.326   -4.029   -1.615)    4.925
   4.745   (   2.105   -3.646    1.480)    4.463
   4.790   (  -1.550    2.684   -0.734)    3.185
   4.822   (  -1.490    2.580   -0.383)    3.004
   4.984   (  -0.713    1.235    4.035)    4.280
   5.005   (  -0.953    1.651    0.711)    2.035
   5.177   (  -3.592    6.221    0.897)    7.239
   5.184   (   0.792   -1.373    1.812)    2.407
   5.282   (  -1.750    3.031   -4.427)    5.643
   5.396   (   0.926   -1.604    1.328)    2.279
   5.449   (   3.162   -5.477   -2.217)    6.701
   5.501   (   3.324   -5.757    0.528)    6.669
   5.526   (   0.298   -0.516   -0.008)    0.596
   5.590   (   0.138   -0.239    2.417)    2.433
   5.838   (  -1.268    2.197   -6.512)    6.988
   5.867   (  -0.168    0.291   -2.306)    2.330
   6.004   (  -1.869    3.238    1.041)    3.881
   6.033   (  -1.892    3.277    1.976)    4.269
   6.285   (   1.867   -3.233   -2.044)    4.256
   6.383   (  -4.404    7.627    6.980)   11.238
   7.415   (  -0.283    0.490    0.663)    0.871
   7.484   (  -1.473    2.551   -0.115)    2.948
   7.512   (  -2.093    3.624   -0.728)    4.248
   7.621   (  -0.918    1.590   -2.489)    3.093
   7.950   (  -0.035    0.060    1.870)    1.871
   7.994   (  -0.580    1.005   -0.365)    1.216
   8.239   (   6.590  -11.414   -1.699)   13.288
   8.342   (   2.445   -4.234   -2.877)    5.673
======================= Grid point 107 (22/35) =======================
q-point: (-0.25  0.38  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.812   (  -3.908    7.388    2.892)    8.844
   0.884   (  -8.501    8.048    3.016)   12.088
   1.046   (   0.102   -0.489   -4.310)    4.339
   1.221   (   2.710    6.759   -2.622)    7.740
   1.357   (  -7.593    1.550   -1.178)    7.839
   1.428   (  -2.787    1.503    3.972)    5.080
   1.707   (   2.396    0.235    1.264)    2.719
   1.724   (  -0.660    1.571   -1.151)    2.056
   1.904   (  -5.157   -0.462    2.197)    5.624
   2.157   (   6.190   -6.818   -7.886)   12.124
   2.299   (   6.730   -0.540    6.312)    9.243
   2.435   (   2.091   -3.558   -1.698)    4.462
   2.465   (   8.240   -4.438    3.836)   10.115
   2.530   (   2.843   -8.463    3.152)    9.468
   2.653   (   0.830   -3.386   -1.171)    3.677
   2.759   (  -4.050   -8.179    0.604)    9.146
   3.064   (   0.309   -2.007   -0.929)    2.233
   3.156   (   4.245    1.192   -1.086)    4.541
   3.200   (   0.134    2.913    5.182)    5.946
   3.288   (  -2.300   -3.792   -3.132)    5.429
   3.372   (  -3.670    0.384   -1.504)    3.984
   3.505   (   4.722    2.519   -0.160)    5.354
   3.573   (  -2.878   -0.234   -2.521)    3.833
   3.637   (   5.095    7.860   -2.902)    9.806
   3.731   (  -0.431    2.668   -2.396)    3.612
   3.873   (  -3.021    6.249   -0.948)    7.005
   3.949   ( -11.819   -1.186    4.020)   12.540
   4.041   (   0.443    2.982    5.922)    6.645
   4.095   (  -1.350    3.301    3.195)    4.789
   4.277   (   2.001   -0.311   -4.721)    5.137
   4.375   (  -2.572    2.980    0.758)    4.009
   4.430   (  -3.193   -1.222    0.086)    3.420
   4.626   (   2.165   -2.399   -0.106)    3.233
   4.667   (   2.986   -3.200    1.413)    4.599
   4.803   (   5.121    5.448   -1.094)    7.556
   4.913   (  -3.579    0.965   -1.501)    3.999
   4.969   (   1.551    0.687    6.079)    6.311
   5.023   (  -1.043   -3.094    0.512)    3.305
   5.120   (   3.619   -0.072   -0.428)    3.645
   5.192   (   4.508   -3.587   -4.687)    7.427
   5.243   (  -0.482   -5.752   -1.328)    5.923
   5.381   (   2.419   -2.721    2.450)    4.388
   5.442   (  -1.774    2.044    0.206)    2.714
   5.456   (  -2.747    1.055   -0.585)    3.000
   5.557   (  -2.020    4.422    0.391)    4.878
   5.602   (  -3.214   -3.600    0.958)    4.919
   5.858   (   1.451    0.843    0.669)    1.806
   5.879   (  -0.030    2.992   -6.488)    7.145
   6.056   (   0.944    4.606    2.980)    5.567
   6.123   (  -3.337    4.098   -1.919)    5.623
   6.244   (  -3.288   -6.912    0.089)    7.654
   6.487   (  -0.641    2.833    5.902)    6.577
   7.430   (  -0.682    0.477    1.183)    1.446
   7.541   (  -1.428    2.487   -1.091)    3.068
   7.587   (  -1.289    2.523   -1.237)    3.092
   7.654   (  -2.476    1.427   -0.033)    2.858
   7.945   (   0.454   -0.493    1.903)    2.018
   7.980   (   5.001   -8.139   -0.713)    9.580
   8.031   (   1.056   -0.506   -0.954)    1.510
   8.235   (   4.739   -3.122   -3.456)    6.644
======================= Grid point 116 (23/35) =======================
q-point: (-0.38  0.50  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.913   (  -0.930    3.247    0.171)    3.382
   1.014   (   0.699    1.158   -1.506)    2.024
   1.118   (  -5.527    6.880    0.732)    8.855
   1.230   (   4.094    4.201   -1.133)    5.974
   1.445   (  -1.344   -2.357    2.317)    3.568
   1.485   (  -1.583    2.754    2.907)    4.306
   1.682   (   0.299    0.477    1.233)    1.355
   1.774   (  -3.892    1.235   -4.053)    5.753
   1.879   (   3.222   -4.317   -2.858)    6.098
   1.977   (   4.920   -5.493   -7.092)   10.231
   2.218   (   0.775    0.570    9.164)    9.214
   2.314   (   1.746   -5.586    5.854)    8.278
   2.365   (   1.496   -2.563    0.731)    3.057
   2.474   (  -1.086   -1.680    0.643)    2.102
   2.654   (  -3.071   -3.690   -5.222)    7.094
   2.751   (  -4.262   -7.865    4.388)    9.964
   3.032   (   0.851   -1.259   -3.483)    3.800
   3.079   (   2.480   -1.120    6.375)    6.931
   3.210   (  -1.167   -0.487   -3.225)    3.464
   3.345   (   1.892    7.250   -2.430)    7.877
   3.441   (   1.471   -1.628   -6.364)    6.732
   3.459   (  -3.269    1.641    4.645)    5.913
   3.605   (  -2.112   -1.537   -3.164)    4.103
   3.653   (   2.608    8.348    3.103)    9.280
   3.719   (   2.868   -0.788    1.738)    3.445
   3.947   (  -0.158    3.092   -3.647)    4.784
   4.073   (  -1.354    1.269   -0.650)    1.967
   4.129   (  -6.392    2.118    5.445)    8.660
   4.159   (   4.859    6.499    0.542)    8.133
   4.282   (  -4.099    3.486   -3.772)    6.571
   4.391   (   2.996   -3.183    0.465)    4.396
   4.446   (  -3.118   -2.749    4.320)    5.995
   4.567   (   3.301    1.025   -1.923)    3.955
   4.679   (   0.162    1.176   -1.181)    1.674
   4.724   (   1.086    1.275    0.078)    1.676
   4.895   (   5.315    1.227    0.506)    5.478
   4.986   (  -0.909    0.659   -0.801)    1.379
   5.005   (   0.727   -0.208    1.595)    1.765
   5.058   (   0.190   -3.994   -3.309)    5.190
   5.086   (   1.403   -2.719   -0.915)    3.193
   5.210   (  -1.472   -3.673   -0.375)    3.975
   5.357   (  -0.672   -0.797    3.076)    3.247
   5.461   (  -0.508   -0.750    0.550)    1.059
   5.508   (  -1.570   -0.075   -0.431)    1.629
   5.617   (  -3.891   -4.044    0.951)    5.692
   5.640   (  -0.233    2.634    0.911)    2.797
   5.846   (   1.442    2.289    5.905)    6.496
   5.900   (   1.058    2.800   -8.857)    9.349
   6.070   (   0.626    1.173    1.819)    2.253
   6.192   (  -1.061    0.205   -6.933)    7.017
   6.278   (  -3.204    0.920    6.984)    7.739
   6.477   (   0.685   -2.705    2.408)    3.686
   7.445   (  -0.470    0.077    1.581)    1.652
   7.582   (  -0.929    1.416   -1.750)    2.435
   7.604   (   1.221    1.678   -1.920)    2.827
   7.723   (  -1.481    0.616   -2.900)    3.314
   7.834   (   1.057   -2.957    3.718)    4.867
   7.925   (   0.743   -1.313    1.097)    1.865
   8.048   (  -0.576    1.203   -0.547)    1.442
   8.128   (   2.855   -2.493   -2.703)    4.655
======================= Grid point 117 (24/35) =======================
q-point: (-0.25  0.50  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.927   (  -2.055    3.559   -0.736)    4.175
   1.024   (   0.460   -0.796   -1.905)    2.116
   1.106   (  -6.524   11.300    2.035)   13.206
   1.367   (  -2.433    4.214    2.442)    5.444
   1.371   (  -0.658    1.140    1.448)    1.957
   1.456   (  -1.116    1.933   -1.981)    2.984
   1.723   (   0.921   -1.596    5.583)    5.879
   1.724   (  -1.288    2.230   -5.335)    5.924
   1.884   (   1.643   -2.845   -2.612)    4.197
   1.980   (   4.624   -8.009   -5.329)   10.674
   2.235   (   2.860   -4.954    7.502)    9.434
   2.304   (   4.947   -8.569    1.620)   10.026
   2.414   (   0.969   -1.679    3.083)    3.642
   2.431   (  -0.270    0.468    0.351)    0.644
   2.585   (   1.683   -2.915   -5.267)    6.250
   2.622   (   1.751   -3.032    6.229)    7.146
   3.051   (   0.254   -0.440    1.439)    1.526
   3.130   (   2.453   -4.249    2.319)    5.427
   3.167   (   0.881   -1.525   -4.384)    4.725
   3.344   (  -2.555    4.426   -7.806)    9.330
   3.429   (  -1.844    3.193    8.393)    9.167
   3.540   (   1.199   -2.077   -6.170)    6.620
   3.553   (  -1.959    3.393    1.289)    4.124
   3.677   (  -0.559    0.968   -8.058)    8.135
   3.896   (  -2.544    4.406    8.255)    9.697
   3.925   (  -1.063    1.841    6.727)    7.054
   3.968   (  -1.909    3.307   -5.106)    6.376
   4.105   (  -1.950    3.378    0.160)    3.904
   4.251   (  -4.481    7.762   -0.691)    8.989
   4.259   (  -0.231    0.400   -2.996)    3.031
   4.389   (   0.375   -0.649    6.219)    6.264
   4.452   (  -0.633    1.096   -3.442)    3.667
   4.576   (   2.767   -4.792    4.594)    7.192
   4.628   (  -2.125    3.681   -2.813)    5.096
   4.832   (   3.748   -6.492   -1.922)    7.739
   4.962   (  -2.796    4.843    0.325)    5.602
   4.977   (  -0.857    1.484   -0.409)    1.761
   4.984   (   0.575   -0.997    4.077)    4.237
   5.034   (   4.156   -7.199   -5.545)    9.993
   5.127   (   0.716   -1.240    1.298)    1.932
   5.152   (   1.439   -2.493   -3.487)    4.521
   5.327   (   1.524   -2.640    3.715)    4.806
   5.451   (   0.159   -0.276   -1.126)    1.170
   5.454   (   1.271   -2.201    0.085)    2.543
   5.583   (  -1.524    2.641    1.582)    3.435
   5.654   (  -1.762    3.053    1.498)    3.830
   5.866   (  -0.009    0.016    5.226)    5.226
   5.925   (  -0.840    1.455   -7.511)    7.697
   6.112   (   1.041   -1.803    1.246)    2.426
   6.195   (  -1.954    3.384   -5.985)    7.148
   6.201   (  -2.504    4.337    3.147)    5.914
   6.480   (   1.998   -3.460    5.276)    6.618
   7.441   (  -0.279    0.483    1.659)    1.750
   7.578   (  -0.714    1.236   -1.611)    2.152
   7.640   (  -0.511    0.886   -3.753)    3.890
   7.706   (  -2.349    4.069    0.830)    4.771
   7.820   (   2.745   -4.755    0.675)    5.531
   7.927   (   0.883   -1.529    1.099)    2.080
   8.048   (  -0.535    0.926   -0.517)    1.188
   8.144   (   2.760   -4.781   -1.977)    5.863
======================= Grid point 124 (25/35) =======================
q-point: (-0.38  0.62  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.968   (   1.033    1.790   -1.300)    2.442
   1.007   (  -0.550   -0.953   -1.693)    2.019
   1.307   (   3.393    5.876   -1.076)    6.870
   1.326   (   4.054    7.022    2.725)    8.554
   1.404   (  -0.662   -1.146    0.313)    1.360
   1.500   (  -1.862   -3.225    2.961)    4.758
   1.689   (   0.575    0.997    4.771)    4.908
   1.787   (  -0.534   -0.925   -6.301)    6.391
   1.819   (  -0.128   -0.221   -1.390)    1.413
   1.877   (  -2.021   -3.500   -4.890)    6.344
   2.183   (  -1.871   -3.240    3.301)    4.989
   2.186   (  -1.896   -3.283    5.970)    7.072
   2.396   (   2.242    3.883    1.138)    4.626
   2.451   (  -0.505   -0.874   -0.574)    1.161
   2.548   (  -3.205   -5.551   -5.011)    8.136
   2.599   (  -2.408   -4.171    8.252)    9.555
   3.054   (   1.045    1.809   -0.475)    2.142
   3.061   (   0.002    0.003    6.878)    6.878
   3.155   (  -2.059   -3.566   -6.696)    7.861
   3.419   (   1.365    2.365  -11.443)   11.764
   3.420   (  -0.057   -0.099    5.539)    5.541
   3.517   (  -2.524   -4.372   -7.016)    8.643
   3.628   (   4.662    8.075   -0.499)    9.337
   3.694   (  -0.029   -0.050    0.601)    0.604
   3.869   (   4.838    8.380    7.976)   12.540
   4.011   (   1.309    2.267   -5.864)    6.421
   4.036   (  -2.119   -3.670    0.666)    4.290
   4.087   (  -2.437   -4.222    8.179)    9.522
   4.327   (   0.246    0.425   -2.470)    2.518
   4.353   (   1.439    2.493   -0.231)    2.888
   4.373   (   1.989    3.445    3.256)    5.140
   4.405   (   2.740    4.746    2.799)    6.154
   4.604   (  -3.509   -6.078    0.223)    7.022
   4.631   (   0.185    0.320   -2.307)    2.336
   4.754   (   1.195    2.070   -3.819)    4.505
   4.877   (  -2.132   -3.693   -6.688)    7.931
   4.994   (   0.683    1.184    2.381)    2.746
   5.008   (   0.353    0.612    0.077)    0.711
   5.049   (  -0.205   -0.354    1.894)    1.938
   5.078   (   2.922    5.061    1.028)    5.933
   5.146   (  -1.187   -2.056   -2.724)    3.613
   5.285   (  -4.316   -7.475    2.842)    9.087
   5.433   (  -0.703   -1.218   -0.808)    1.622
   5.446   (  -0.443   -0.767   -0.944)    1.295
   5.621   (   1.440    2.494    1.777)    3.384
   5.677   (  -0.288   -0.498    2.951)    3.007
   5.868   (   0.095    0.165    6.582)    6.585
   5.939   (   0.399    0.692   -6.235)    6.286
   6.089   (   0.527    0.912   -2.764)    2.958
   6.226   (   2.161    3.743   -2.672)    5.081
   6.304   (  -0.582   -1.009    3.307)    3.506
   6.395   (  -1.928   -3.339    2.913)    4.833
   7.447   (   0.037    0.064    1.767)    1.768
   7.599   (   0.377    0.652   -2.544)    2.653
   7.627   (  -0.232   -0.402   -4.105)    4.131
   7.753   (   1.983    3.435    0.584)    4.009
   7.812   (  -0.559   -0.969    3.658)    3.825
   7.887   (  -1.101   -1.907   -1.035)    2.434
   8.067   (   0.111    0.192   -0.437)    0.490
   8.076   (  -0.821   -1.423   -1.012)    1.929
======================= Grid point 146 (26/35) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 20
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.695   (   0.000    0.000  -16.185)   16.185
   0.695   (  -0.000    0.000  -16.185)   16.185
   0.695   (   0.000    0.000   16.185)   16.185
   0.695   (   0.000   -0.000   16.185)   16.185
   1.160   (  -0.000   -0.000  -13.630)   13.630
   1.160   (   0.000    0.000   13.630)   13.630
   2.046   (   0.000    0.000  -12.177)   12.177
   2.046   (   0.000    0.000  -12.177)   12.177
   2.046   (   0.000    0.000   12.177)   12.177
   2.046   (   0.000   -0.000   12.177)   12.177
   2.075   (  -0.000   -0.000  -10.550)   10.550
   2.075   (   0.000    0.000   10.550)   10.550
   2.963   (  -0.000   -0.000   -6.905)    6.905
   2.963   (   0.000    0.000    6.905)    6.905
   3.130   (   0.000    0.000   -4.641)    4.641
   3.130   (  -0.000   -0.000   -4.641)    4.641
   3.130   (   0.000    0.000    4.641)    4.641
   3.130   (   0.000   -0.000    4.641)    4.641
   3.437   (   0.000    0.000   -2.895)    2.895
   3.437   (   0.000   -0.000   -2.895)    2.895
   3.437   (   0.000    0.000    2.895)    2.895
   3.437   (   0.000   -0.000    2.895)    2.895
   3.473   (  -0.000   -0.000   -1.811)    1.811
   3.473   (   0.000    0.000    1.811)    1.811
   3.638   (   0.000    0.000   -6.158)    6.158
   3.638   (   0.000   -0.000   -6.158)    6.158
   3.638   (   0.000    0.000    6.158)    6.158
   3.638   (   0.000    0.000    6.158)    6.158
   4.173   (  -0.000   -0.000   -2.005)    2.005
   4.173   (   0.000    0.000    2.005)    2.005
   4.574   (  -0.000    0.000  -17.245)   17.245
   4.574   (   0.000    0.000   17.245)   17.245
   4.793   (   0.000    0.000   -1.050)    1.050
   4.793   (   0.000    0.000   -1.050)    1.050
   4.793   (   0.000    0.000    1.050)    1.050
   4.793   (  -0.000   -0.000    1.050)    1.050
   4.905   (   0.000    0.000   -0.469)    0.469
   4.905   (  -0.000    0.000   -0.469)    0.469
   4.905   (   0.000    0.000    0.469)    0.469
   4.905   (   0.000    0.000    0.469)    0.469
   5.585   (   0.000    0.000   -3.785)    3.785
   5.585   (   0.000   -0.000   -3.785)    3.785
   5.585   (   0.000   -0.000    3.785)    3.785
   5.585   (   0.000    0.000    3.785)    3.785
   5.602   (   0.000    0.000   -7.689)    7.689
   5.602   (  -0.000   -0.000    7.689)    7.689
   5.962   (   0.000    0.000   -1.502)    1.502
   5.962   (  -0.000   -0.000   -1.502)    1.502
   5.962   (   0.000    0.000    1.502)    1.502
   5.962   (   0.000    0.000    1.502)    1.502
   6.002   (   0.000    0.000   -1.436)    1.436
   6.002   (  -0.000   -0.000    1.436)    1.436
   7.419   (   0.000    0.000   -0.095)    0.095
   7.419   (   0.000    0.000    0.095)    0.095
   7.419   (   0.000    0.000   -0.095)    0.095
   7.419   (  -0.000    0.000    0.095)    0.095
   7.949   (  -0.000   -0.000   -1.115)    1.115
   7.949   (   0.000    0.000    1.115)    1.115
   8.527   (  -0.000   -0.000   -0.223)    0.223
   8.527   (   0.000    0.000    0.223)    0.223
======================= Grid point 148 (27/35) =======================
q-point: ( 0.12  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 68
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.736   (   4.502    2.048  -12.356)   13.309
   0.736   (   4.025    2.875   12.356)   13.309
   0.776   (   5.757    1.093  -12.545)   13.846
   0.776   (   3.825    4.439   12.545)   13.846
   1.242   (   5.404    4.784   -7.444)   10.369
   1.242   (   6.846    2.288    7.444)   10.369
   1.857   (  -7.763   -3.420   -2.559)    8.861
   1.857   (  -6.843   -5.013    2.559)    8.861
   2.047   (   3.158   -5.385   -7.460)    9.728
   2.047   (  -3.084    5.427    7.460)    9.728
   2.263   (   8.363    3.306   -5.084)   10.330
   2.263   (   7.045    5.589    5.084)   10.330
   2.869   (  -8.320   -2.427   -5.524)   10.278
   2.869   (  -6.262   -5.992    5.524)   10.278
   3.005   (  -8.010  -11.357   -0.595)   13.910
   3.005   ( -13.840   -1.258    0.595)   13.910
   3.183   (   2.052    3.381   -3.526)    5.298
   3.183   (   3.954    0.086    3.526)    5.298
   3.246   (  -2.517   -2.416   -4.432)    5.641
   3.246   (  -3.351   -0.972    4.432)    5.641
   3.408   (  -1.741   -1.088   -2.944)    3.589
   3.408   (  -1.813   -0.963    2.944)    3.589
   3.521   (   4.142   -7.830    1.109)    8.927
   3.521   (  -4.710    7.502   -1.109)    8.927
   3.613   (  -1.310    4.537   -2.667)    5.423
   3.613   (   3.274   -3.403    2.667)    5.423
   3.904   (  15.936   -2.480   -3.687)   16.544
   3.904   (   5.820   15.041    3.687)   16.544
   4.191   (   2.080    0.653   -1.489)    2.640
   4.191   (   1.605    1.475    1.489)    2.640
   4.446   (  -8.243   -2.726  -12.015)   14.823
   4.446   (  -6.483   -5.776   12.015)   14.823
   4.760   (  -2.141   -1.791   -0.581)    2.851
   4.760   (  -2.622   -0.959    0.581)    2.851
   4.795   (   1.459    0.540   -2.294)    2.772
   4.795   (   1.197    0.994    2.294)    2.772
   4.961   (   3.366    1.736   -0.366)    3.805
   4.961   (   3.186    2.047    0.366)    3.805
   5.047   (   5.829    7.174   -2.589)    9.599
   5.047   (   9.127    1.461    2.589)    9.599
   5.458   (  -6.577   -6.020   -5.437)   10.444
   5.458   (  -8.502   -2.686    5.437)   10.444
   5.560   (  -3.408   -0.438   -3.067)    4.605
   5.560   (  -2.083   -2.732    3.067)    4.605
   5.597   (   0.247   -6.060    0.367)    6.076
   5.597   (  -5.124    3.244   -0.367)    6.076
   5.765   (  -0.017   -2.784   -1.993)    3.424
   5.765   (  -2.420    1.377    1.993)    3.424
   5.978   (   1.781   -0.146   -0.613)    1.889
   5.978   (   0.764    1.615    0.613)    1.889
   6.295   (   6.526   10.767   -4.261)   13.292
   6.295   (  12.588    0.268    4.261)   13.292
   7.427   (   1.367   -0.320   -0.776)    1.604
   7.427   (   0.407    1.344    0.776)    1.604
   7.482   (   5.173    0.818   -1.779)    5.532
   7.482   (   3.295    4.071    1.779)    5.532
   7.963   (   1.017    0.599   -0.868)    1.465
   7.963   (   1.027    0.581    0.868)    1.465
   8.411   (  -8.131   -5.359    0.115)    9.739
   8.411   (  -8.706   -4.363   -0.115)    9.739
======================= Grid point 150 (28/35) =======================
q-point: ( 0.25  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 70
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.851   (   4.828    1.822   -3.456)    6.211
   0.851   (   3.992    3.271    3.456)    6.211
   0.922   (   8.433    1.676   -5.736)   10.336
   0.922   (   5.668    6.465    5.736)   10.336
   1.388   (   5.461    1.140   -2.356)    6.056
   1.388   (   3.717    4.160    2.356)    6.056
   1.733   (  -1.067   -4.342   -0.665)    4.520
   1.733   (  -4.293    1.247    0.665)    4.520
   2.045   (   1.944   -4.953   -4.004)    6.659
   2.045   (  -3.318    4.160    4.004)    6.659
   2.400   (   1.161    2.886   -9.101)    9.618
   2.400   (   3.080   -0.438    9.101)    9.618
   2.625   ( -10.132   -4.760   -3.561)   11.747
   2.625   (  -9.188   -6.395    3.561)   11.747
   2.799   (  -4.247   -4.831   -3.495)    7.321
   2.799   (  -6.307   -1.263    3.495)    7.321
   3.192   (  -2.081   -2.776   -3.559)    4.970
   3.192   (  -3.445   -0.414    3.559)    4.970
   3.249   (   3.860   -2.812   -2.179)    5.249
   3.249   (  -0.505    4.749    2.179)    5.249
   3.373   (   0.842   -3.473   -2.435)    4.325
   3.373   (  -2.587    2.466    2.435)    4.325
   3.511   (  -0.884    3.503   -3.534)    5.053
   3.511   (   2.591   -2.517    3.534)    5.053
   3.685   (   7.809   -2.761    0.984)    8.341
   3.685   (   1.514    8.143   -0.984)    8.341
   4.067   (   5.656   -2.177   -0.194)    6.064
   4.067   (   0.943    5.987    0.194)    6.064
   4.200   (  -3.922   -1.485   -5.209)    6.687
   4.200   (  -3.247   -2.654    5.209)    6.687
   4.381   (   5.862   -0.873   -2.501)    6.433
   4.381   (   2.175    5.513    2.501)    6.433
   4.687   (  -3.459   -2.617   -0.297)    4.347
   4.687   (  -3.996   -1.687    0.297)    4.347
   4.824   (   1.230   -1.342   -4.017)    4.410
   4.824   (  -0.547    1.737    4.017)    4.410
   5.027   (   1.012    2.687   -2.609)    3.880
   5.027   (   2.833   -0.467    2.609)    3.880
   5.140   (  -1.460   -0.386   -2.542)    2.957
   5.140   (  -1.064   -1.072    2.542)    2.957
   5.356   (   1.925   -2.802   -2.009)    3.949
   5.356   (  -1.464    3.068    2.009)    3.949
   5.474   (  -2.000   -1.789   -1.262)    2.966
   5.474   (  -2.549   -0.838    1.262)    2.966
   5.562   (   1.220   -3.413   -3.916)    5.336
   5.562   (  -2.345    2.763    3.916)    5.336
   5.794   (   3.923   -2.409   -0.703)    4.657
   5.794   (  -0.125    4.602    0.703)    4.657
   6.022   (   1.798    2.099   -1.829)    3.314
   6.022   (   2.717    0.507    1.829)    3.314
   6.437   (   5.466   -2.940   -0.859)    6.265
   6.437   (   0.187    6.203    0.859)    6.265
   7.457   (   1.935    0.095   -1.825)    2.662
   7.457   (   1.050    1.628    1.825)    2.662
   7.580   (   3.960    0.649   -2.371)    4.661
   7.580   (   2.542    3.105    2.371)    4.661
   7.988   (   0.966    0.195   -1.049)    1.439
   7.988   (   0.651    0.739    1.049)    1.439
   8.189   (  -6.776   -8.418   -1.622)   10.927
   8.189   ( -10.678   -1.659    1.622)   10.927
======================= Grid point 152 (29/35) =======================
q-point: ( 0.38  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 68
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.928   (   2.947   -0.961   -2.720)    4.124
   0.928   (   0.641    3.033    2.720)    4.124
   1.080   (   5.942    1.572   -3.838)    7.246
   1.080   (   4.333    4.360    3.838)    7.246
   1.476   (   3.677    2.588   -4.440)    6.319
   1.476   (   4.080    1.890    4.440)    6.319
   1.701   (   0.608   -2.907   -2.654)    3.982
   1.701   (  -2.213    1.980    2.654)    3.982
   1.966   (  -5.036   -5.275   -1.314)    7.410
   1.966   (  -7.086   -1.724    1.314)    7.410
   2.358   (  -4.665    3.189   -7.944)    9.749
   2.358   (   0.429   -5.634    7.944)    9.749
   2.506   (  -2.257   -1.482   -4.641)    5.370
   2.506   (  -2.412   -1.214    4.641)    5.370
   2.732   (   2.202   -6.252   -1.635)    6.827
   2.732   (  -4.313    5.033    1.635)    6.827
   3.130   (   1.199   -5.915   -1.578)    6.238
   3.130   (  -4.523    3.996    1.578)    6.238
   3.263   (   2.899   -4.685   -7.329)    9.169
   3.263   (  -2.608    4.853    7.329)    9.169
   3.358   (   2.391   -5.883   -2.835)    6.954
   3.358   (  -3.899    5.012    2.835)    6.954
   3.578   (   2.627    4.244   -5.522)    7.443
   3.578   (   4.989    0.153    5.522)    7.443
   3.810   (   7.644   -2.524   -1.759)    8.240
   3.810   (   1.636    7.882    1.759)    8.240
   4.062   (  -4.143    0.757   -5.208)    6.698
   4.062   (  -1.416   -3.967    5.208)    6.698
   4.177   (   3.246    0.682   -3.310)    4.685
   4.177   (   2.213    2.470    3.310)    4.685
   4.466   (   2.578   -5.208   -0.771)    5.862
   4.466   (  -3.221    4.837    0.771)    5.862
   4.609   (   0.240   -1.838   -2.695)    3.271
   4.609   (  -1.472    1.127    2.695)    3.271
   4.808   (   1.096   -4.730   -4.842)    6.857
   4.808   (  -3.548    3.314    4.842)    6.857
   5.035   (  -3.430   -0.570   -4.580)    5.750
   5.035   (  -2.209   -2.686    4.580)    5.750
   5.076   (  -0.810   -0.940   -2.559)    2.844
   5.076   (  -1.219   -0.232    2.559)    2.844
   5.347   (   2.447   -5.176   -1.686)    5.968
   5.347   (  -3.259    4.707    1.686)    5.968
   5.511   (   4.357   -0.862   -2.245)    4.977
   5.511   (   1.432    4.205    2.245)    4.977
   5.559   (   2.377   -3.575   -4.677)    6.349
   5.559   (  -1.908    3.846    4.677)    6.349
   5.833   (   2.123   -0.894    0.190)    2.312
   5.833   (   0.288    2.286   -0.190)    2.312
   6.079   (   1.451    2.498   -2.822)    4.039
   6.079   (   2.889    0.008    2.822)    4.039
   6.462   (   1.196   -2.666   -4.588)    5.439
   6.462   (  -1.711    2.369    4.588)    5.439
   7.494   (   1.729    0.097   -2.790)    3.284
   7.494   (   0.949    1.449    2.790)    3.284
   7.627   (   1.100   -0.279   -3.458)    3.640
   7.627   (   0.308    1.092    3.458)    3.640
   7.975   (  -3.025   -0.500   -1.289)    3.325
   7.975   (  -1.945   -2.369    1.289)    3.325
   8.048   (  -0.534   -3.486   -4.228)    5.506
   8.048   (  -3.286    1.280    4.228)    5.506
======================= Grid point 154 (30/35) =======================
q-point: (-0.50  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 42
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.947   (   1.010   -1.750   -1.569)    2.559
   0.947   (  -1.010    1.750    1.569)    2.559
   1.143   (  -0.039    0.067   -0.288)    0.298
   1.143   (   0.039   -0.067    0.288)    0.298
   1.551   (   1.635   -2.831   -0.504)    3.308
   1.551   (  -1.635    2.831    0.504)    3.308
   1.692   (  -0.348    0.602   -5.756)    5.798
   1.692   (   0.348   -0.602    5.756)    5.798
   1.851   (  -0.298    0.517   -1.884)    1.976
   1.851   (   0.298   -0.517    1.884)    1.976
   2.336   (  -2.422    4.196   -8.105)    9.443
   2.336   (   2.422   -4.196    8.105)    9.443
   2.479   (   0.565   -0.978   -2.129)    2.410
   2.479   (  -0.565    0.978    2.129)    2.410
   2.718   (   3.378   -5.851   -3.050)    7.413
   2.718   (  -3.378    5.851    3.050)    7.413
   3.113   (   5.882  -10.189    0.161)   11.766
   3.113   (  -5.882   10.189   -0.161)   11.766
   3.270   (  -1.469    2.545   -9.137)    9.598
   3.270   (   1.469   -2.545    9.137)    9.598
   3.342   (   4.698   -8.137   -0.437)    9.406
   3.342   (  -4.698    8.137    0.437)    9.406
   3.636   (  -0.021    0.037   -4.437)    4.437
   3.636   (   0.021   -0.037    4.437)    4.437
   3.864   (   1.425   -2.468   -4.185)    5.063
   3.864   (  -1.425    2.468    4.185)    5.063
   4.030   (  -3.852    6.672   -4.358)    8.852
   4.030   (   3.852   -6.672    4.358)    8.852
   4.243   (   0.124   -0.214   -2.254)    2.267
   4.243   (  -0.124    0.214    2.254)    2.267
   4.403   (   1.256   -2.175   -0.462)    2.554
   4.403   (  -1.256    2.175    0.462)    2.554
   4.636   (   1.432   -2.481   -1.533)    3.249
   4.636   (  -1.432    2.481    1.533)    3.249
   4.786   (  -0.852    1.475   -1.236)    2.105
   4.786   (   0.852   -1.475    1.236)    2.105
   4.976   (   0.016   -0.027   -2.453)    2.453
   4.976   (  -0.016    0.027    2.453)    2.453
   5.074   (   0.575   -0.995   -0.191)    1.165
   5.074   (  -0.575    0.995    0.191)    1.165
   5.347   (   2.885   -4.997   -2.343)    6.228
   5.347   (  -2.885    4.997    2.343)    6.228
   5.543   (   1.619   -2.804   -0.097)    3.239
   5.543   (  -1.619    2.804    0.097)    3.239
   5.569   (   1.264   -2.189   -3.444)    4.272
   5.569   (  -1.264    2.189    3.444)    4.272
   5.849   (  -4.574    7.922   -4.283)   10.100
   5.849   (   4.574   -7.922    4.283)   10.100
   6.109   (  -0.776    1.344   -3.811)    4.115
   6.109   (   0.776   -1.344    3.811)    4.115
   6.453   (  -0.626    1.083   -7.192)    7.300
   6.453   (   0.626   -1.083    7.192)    7.300
   7.512   (   0.298   -0.517   -3.294)    3.348
   7.512   (  -0.298    0.517    3.294)    3.348
   7.631   (   0.387   -0.671   -4.205)    4.276
   7.631   (  -0.387    0.671    4.205)    4.276
   7.935   (  -0.280    0.485   -2.827)    2.882
   7.935   (   0.280   -0.485    2.827)    2.882
   8.036   (   0.219   -0.379   -1.503)    1.565
   8.036   (  -0.219    0.379    1.503)    1.565
======================= Grid point 166 (31/35) =======================
q-point: ( 0.12  0.12 -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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 76
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.814   (   2.048    3.547   -8.324)    9.277
   0.819   (   3.392    5.876   -2.626)    7.275
   0.853   (   3.979    6.892   -0.100)    7.958
   0.929   (   5.009    8.675    9.544)   13.836
   1.323   (   2.586    4.478    0.045)    5.171
   1.374   (   3.339    5.784   -0.097)    6.680
   1.708   (  -2.900   -5.022    1.541)    6.000
   1.836   (  -1.362   -2.358    0.596)    2.788
   1.948   (  -1.261   -2.184   -6.612)    7.077
   2.121   (   0.840    1.455    3.174)    3.591
   2.364   (   0.760    1.317    8.915)    9.044
   2.369   (   3.211    5.561   -5.898)    8.719
   2.656   (  -6.031  -10.445    3.227)   12.486
   2.771   (  -4.491   -7.779   -4.555)   10.071
   2.805   (  -3.499   -6.061   -0.121)    7.000
   2.875   (  -6.744  -11.681   -1.248)   13.545
   3.118   (  -1.371   -2.375   -0.097)    2.744
   3.192   (  -1.467   -2.540   -3.432)    4.515
   3.214   (  -1.163   -2.014    5.013)    5.526
   3.288   (   3.464    6.000   -6.162)    9.272
   3.370   (  -2.446   -4.236   -0.857)    4.966
   3.393   (   1.140    1.974    5.867)    6.295
   3.465   (   1.888    3.270    1.422)    4.035
   3.546   (  -1.024   -1.774    3.089)    3.706
   3.680   (   2.950    5.110    2.868)    6.561
   3.815   (   1.669    2.891   -3.008)    4.493
   3.829   (   3.035    5.257   -2.628)    6.615
   4.098   (   1.517    2.627    0.687)    3.110
   4.221   (   0.935    1.620   -1.695)    2.524
   4.239   (  -1.403   -2.430    4.896)    5.643
   4.335   (  -3.284   -5.688   -8.466)   10.715
   4.387   (   3.783    6.553    2.897)    8.102
   4.707   (  -1.728   -2.994   -0.753)    3.538
   4.737   (  -0.561   -0.971    0.306)    1.163
   4.809   (  -0.284   -0.492   -1.730)    1.821
   4.828   (   1.207    2.090   -0.854)    2.560
   4.990   (  -0.164   -0.283    2.272)    2.296
   5.044   (   2.502    4.333   -1.761)    5.304
   5.109   (   1.970    3.411    6.084)    7.248
   5.166   (   0.758    1.313    2.301)    2.755
   5.325   (  -3.226   -5.587   -5.288)    8.342
   5.386   (  -0.081   -0.140    2.897)    2.902
   5.463   (  -2.611   -4.522   -4.056)    6.611
   5.496   (  -2.153   -3.730    0.702)    4.364
   5.536   (  -0.425   -0.736   -0.764)    1.143
   5.640   (  -0.709   -1.228   -1.388)    1.984
   5.719   (   0.030    0.052    3.349)    3.349
   5.821   (   2.113    3.660    1.581)    4.512
   5.978   (   0.962    1.665    2.051)    2.812
   6.058   (   2.073    3.591   -0.237)    4.153
   6.292   (   0.101    0.174   -0.672)    0.701
   6.517   (   4.246    7.354   -5.133)    9.923
   7.438   (   0.834    1.445   -0.977)    1.933
   7.478   (   1.311    2.271    0.547)    2.679
   7.521   (   1.813    3.141   -2.444)    4.373
   7.574   (   2.227    3.858    1.997)    4.882
   7.980   (   0.448    0.776   -0.947)    1.304
   7.984   (   0.748    1.295    0.716)    1.658
   8.205   (  -6.790  -11.761    1.081)   13.623
   8.288   (  -4.208   -7.289    1.749)    8.597
======================= Grid point 168 (32/35) =======================
q-point: ( 0.25  0.12 -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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.892   (   2.170    2.800   -4.375)    5.630
   0.943   (   2.897    4.710    1.489)    5.727
   1.016   (   5.539    6.270   -5.722)   10.135
   1.105   (   2.383   10.378    6.920)   12.699
   1.412   (   4.698    1.421   -4.465)    6.635
   1.445   (   1.303    1.074    3.940)    4.287
   1.656   (   0.613   -2.114   -1.481)    2.653
   1.780   (  -3.669    1.502    1.818)    4.362
   1.935   (   1.396   -4.268   -1.455)    4.721
   2.065   (  -5.945   -4.496   -0.765)    7.493
   2.305   (  -1.508   -4.565    7.874)    9.226
   2.453   (   0.186    1.463   -7.657)    7.798
   2.525   (  -2.826   -3.553    4.409)    6.328
   2.530   (  -7.628   -6.469   -0.948)   10.046
   2.666   (  -3.784   -7.735   -1.011)    8.670
   2.777   (   2.719   -2.035   -1.250)    3.619
   3.059   (  -2.581   -4.669   -2.004)    5.699
   3.160   (  -0.861   -1.121    1.037)    1.753
   3.188   (   1.965   -0.364    5.076)    5.456
   3.291   (  -0.029   -4.396   -3.888)    5.869
   3.410   (  -3.904    4.657    2.594)    6.607
   3.435   (   0.626    5.239   -2.729)    5.941
   3.503   (   4.199    2.510    3.782)    6.183
   3.580   (   5.372   -2.270    0.039)    5.832
   3.734   (  -2.458    8.223    0.860)    8.625
   3.865   (   1.020    7.749   -0.927)    7.870
   4.016   (   7.480   -0.157    0.178)    7.484
   4.107   (  -0.227   -0.343   -1.399)    1.458
   4.160   (  -2.958   -0.406   -0.447)    3.019
   4.188   (  -0.967   -0.194   -0.259)    1.020
   4.350   (   4.999   -0.149    1.528)    5.230
   4.545   (   2.291    6.447   -5.974)    9.083
   4.631   (  -3.856   -1.576   -0.510)    4.197
   4.688   (  -3.579   -1.199   -0.727)    3.844
   4.777   (  -3.278    0.580    0.122)    3.331
   4.869   (   1.319    0.727    2.791)    3.171
   4.989   (   2.643   -2.135    1.760)    3.826
   5.079   (  -4.329   -3.247    2.697)    6.046
   5.099   (  -0.577    2.395   -0.915)    2.628
   5.127   (  -3.402   -1.702   -0.370)    3.822
   5.244   (   0.892   -7.113   -1.917)    7.421
   5.426   (  -1.399    0.516   -1.531)    2.137
   5.447   (   2.226    0.346   -0.993)    2.462
   5.457   (   2.690    0.902   -0.370)    2.862
   5.542   (   1.073    1.956    1.165)    2.516
   5.600   (   0.013   -2.009    1.745)    2.661
   5.732   (   0.468   -0.591    5.272)    5.326
   5.896   (   1.690    3.552   -1.824)    4.336
   6.031   (   2.388    1.887    2.860)    4.176
   6.121   (  -0.914    5.665   -1.564)    5.948
   6.309   (   5.078   -5.143   -2.891)    7.784
   6.583   (  -2.620    3.111   -2.478)    4.763
   7.478   (   1.330    1.883   -2.463)    3.374
   7.516   (  -0.310    3.442    1.005)    3.599
   7.587   (   2.658    0.034   -1.307)    2.962
   7.632   (   1.064    1.824    1.327)    2.494
   7.973   (  -3.294   -5.933    0.581)    6.811
   7.987   (  -0.803   -2.390   -0.954)    2.695
   8.017   (   0.529   -0.792   -1.848)    2.078
   8.140   (  -5.343   -3.053    4.854)    7.837
======================= Grid point 170 (33/35) =======================
q-point: ( 0.38  0.12 -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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.940   (   0.364    1.947   -2.311)    3.044
   0.997   (  -0.783    2.475    0.004)    2.596
   1.154   (   2.140    5.594   -2.524)    6.499
   1.199   (  -2.525    7.061    1.929)    7.743
   1.476   (   2.492   -1.966   -0.004)    3.174
   1.545   (   5.664    2.720   -1.012)    6.364
   1.672   (   2.229    0.582    0.510)    2.359
   1.734   (  -2.776    1.523    0.435)    3.197
   1.862   (  -2.934   -3.976   -1.174)    5.078
   1.897   (  -5.819   -4.184   -0.436)    7.180
   2.234   (   0.682   -5.790    7.381)    9.405
   2.368   (  -0.863   -3.992   -3.508)    5.384
   2.463   (  -1.130   -2.095   -2.987)    3.820
   2.481   (   0.055   -1.149    1.311)    1.744
   2.594   (   3.252   -6.064    0.073)    6.881
   2.799   (   0.944   -1.270    0.739)    1.747
   2.970   (  -1.285   -3.145   -0.942)    3.525
   3.145   (   1.935   -3.154    4.030)    5.471
   3.240   (   2.484    0.381    3.159)    4.037
   3.278   (   0.371   -0.782   -9.110)    9.151
   3.327   (  -4.537    0.983    4.544)    6.496
   3.478   (  -2.295    4.634    0.330)    5.182
   3.637   (   3.412    4.178    2.981)    6.163
   3.652   (   3.721   -1.089   -4.477)    5.922
   3.793   (  -0.092    6.081   -1.945)    6.385
   3.894   (  -3.324    1.026    0.358)    3.498
   4.007   (  -2.189   -0.216    3.491)    4.127
   4.179   (   0.103    2.933    0.210)    2.942
   4.203   (   3.818    4.891    1.284)    6.337
   4.235   (   3.247    5.290   -2.540)    6.706
   4.389   (  -0.741   -1.141    0.026)    1.361
   4.520   (  -8.126    1.275    1.778)    8.415
   4.574   (   0.587   -3.219   -4.621)    5.663
   4.656   (   0.799    1.566    0.431)    1.810
   4.745   (   1.524   -0.549   -2.546)    3.018
   4.855   (  -4.615   -0.563    3.145)    5.613
   4.986   (  -1.134   -1.072    1.174)    1.953
   5.009   (  -3.416   -1.066    0.592)    3.627
   5.039   (   2.674   -0.140   -0.991)    2.855
   5.105   (  -2.227    1.425   -3.725)    4.568
   5.220   (   3.807   -6.647   -1.255)    7.762
   5.409   (  -0.380    0.266   -0.938)    1.047
   5.482   (   1.246   -1.620    1.465)    2.515
   5.486   (   2.424   -2.913   -1.771)    4.183
   5.599   (   1.318   -0.485    3.134)    3.434
   5.619   (   1.192    3.766    2.618)    4.739
   5.720   (  -0.141   -1.639    7.026)    7.216
   5.953   (   0.810    2.372   -3.766)    4.524
   6.086   (   0.830    2.211    2.702)    3.588
   6.128   (  -1.811    1.177   -3.506)    4.117
   6.386   (   6.335   -0.203   -1.150)    6.442
   6.511   (  -5.281   -1.673   -1.895)    5.855
   7.513   (   0.296    1.458   -4.035)    4.301
   7.538   (  -0.739    2.427    2.350)    3.458
   7.610   (   0.988   -0.961   -3.104)    3.397
   7.651   (  -0.887    1.371    3.134)    3.534
   7.894   (   1.916   -3.322   -0.422)    3.858
   7.942   (  -1.949   -2.868   -0.761)    3.550
   8.033   (   0.435    0.636    0.259)    0.813
   8.056   (  -1.894   -0.353    3.152)    3.694
======================= Grid point 184 (34/35) =======================
q-point: ( 0.25  0.25 -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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 78
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.950   (   1.369    2.371   -2.755)    3.884
   1.013   (   0.819    1.419   -0.952)    1.894
   1.169   (   4.617    7.996   -3.800)    9.985
   1.318   (   4.729    8.191    5.681)   11.033
   1.444   (  -0.488   -0.845    0.354)    1.038
   1.453   (   1.588    2.750   -4.790)    5.747
   1.650   (   0.808    1.400    1.852)    2.458
   1.783   (  -0.368   -0.637    0.908)    1.169
   1.861   (  -2.035   -3.525   -0.715)    4.133
   1.923   (  -3.613   -6.258   -1.696)    7.423
   2.224   (  -2.544   -4.407    6.026)    7.887
   2.328   (  -5.478   -9.487    0.826)   10.986
   2.469   (  -0.433   -0.749   -4.394)    4.478
   2.483   (  -0.410   -0.710    0.154)    0.834
   2.508   (  -3.262   -5.651    2.116)    6.860
   2.738   (  -0.672   -1.165   -3.392)    3.649
   2.996   (  -1.542   -2.670    0.187)    3.089
   3.157   (   0.052    0.090    9.655)    9.656
   3.223   (   1.471    2.548   -3.671)    4.704
   3.226   (  -0.488   -0.846   -7.387)    7.452
   3.365   (  -3.153   -5.461    6.767)    9.250
   3.535   (  -0.315   -0.546    1.039)    1.215
   3.585   (   3.282    5.684    2.464)    7.011
   3.615   (   3.747    6.490   -6.051)    9.632
   3.851   (   1.675    2.901    4.485)    5.598
   3.997   (   1.272    2.202   -1.512)    2.959
   4.035   (   0.535    0.926   -3.194)    3.368
   4.113   (   0.623    1.079    1.121)    1.676
   4.213   (   3.326    5.760   -0.044)    6.651
   4.232   (   3.563    6.171    0.681)    7.158
   4.395   (   2.264    3.921    2.006)    4.953
   4.555   (  -2.336   -4.045   -7.486)    8.824
   4.609   (  -0.520   -0.901   -0.448)    1.133
   4.653   (  -1.094   -1.894   -3.556)    4.175
   4.779   (  -2.271   -3.934    2.656)    5.262
   4.874   (  -1.011   -1.751    4.880)    5.282
   4.978   (   0.382    0.661   -0.739)    1.062
   5.011   (   0.183    0.318   -0.700)    0.790
   5.066   (  -0.908   -1.573    4.594)    4.940
   5.099   (  -2.152   -3.728   -3.959)    5.848
   5.166   (   0.360    0.624   -4.388)    4.447
   5.397   (  -2.277   -3.945   -2.695)    5.292
   5.437   (  -0.239   -0.414    0.108)    0.490
   5.446   (  -1.411   -2.443    0.907)    2.964
   5.583   (   1.281    2.219    1.426)    2.932
   5.618   (   2.449    4.242    4.719)    6.802
   5.742   (   0.035    0.061    8.782)    8.782
   5.937   (   0.546    0.945   -2.039)    2.313
   6.087   (   1.511    2.618    1.231)    3.264
   6.225   (  -0.329   -0.569   -5.019)    5.062
   6.260   (   3.029    5.246   -2.084)    6.406
   6.557   (  -2.905   -5.032   -1.321)    5.959
   7.522   (   1.073    1.859   -4.677)    5.146
   7.573   (   0.609    1.054   -0.930)    1.532
   7.583   (  -0.060   -0.103    0.250)    0.277
   7.678   (   1.672    2.896    3.177)    4.612
   7.841   (  -1.969   -3.410   -1.177)    4.110
   7.944   (  -1.893   -3.278   -0.164)    3.789
   8.032   (   0.458    0.793    0.526)    1.056
   8.076   (  -1.385   -2.400    3.678)    4.605
======================= Grid point 186 (35/35) =======================
q-point: ( 0.38  0.25 -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 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 1.25e-05 1.25e-05 
Number of triplets: 76
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.983   (  -0.950    1.646   -2.528)    3.163
   1.015   (   0.171   -0.296   -2.333)    2.357
   1.297   (  -3.555    6.158   -0.068)    7.111
   1.367   (  -4.844    8.391   -0.205)    9.691
   1.427   (   0.927   -1.606   -2.377)    3.014
   1.524   (   2.927   -5.070    1.554)    6.057
   1.665   (   1.570   -2.720    5.581)    6.404
   1.787   (  -1.559    2.701   -4.498)    5.473
   1.809   (   0.128   -0.221    1.000)    1.033
   1.847   (   0.322   -0.557   -4.743)    4.787
   2.152   (   0.538   -0.932    5.402)    5.508
   2.228   (   3.867   -6.697    2.371)    8.088
   2.417   (   1.124   -1.946    2.149)    3.109
   2.471   (  -0.090    0.156   -2.095)    2.103
   2.494   (   1.620   -2.806   -1.989)    3.802
   2.723   (   1.872   -3.243   -0.636)    3.798
   2.966   (  -0.221    0.383    0.254)    0.510
   3.156   (  -2.293    3.971    7.210)    8.545
   3.236   (   0.297   -0.515   -4.065)    4.108
   3.254   (   0.742   -1.285  -11.927)   12.019
   3.282   (   2.345   -4.062   10.497)   11.498
   3.524   (   0.497   -0.860   -0.451)    1.091
   3.656   (  -2.056    3.561   -1.701)    4.450
   3.699   (  -0.432    0.749    0.825)    1.196
   3.892   (  -0.692    1.199    4.867)    5.060
   3.990   (  -3.182    5.512   -3.515)    7.271
   4.025   (  -1.406    2.435    1.915)    3.402
   4.132   (   2.254   -3.904    2.778)    5.295
   4.320   (  -2.810    4.868   -1.892)    5.931
   4.397   (  -2.459    4.259   -2.885)    5.701
   4.418   (  -2.588    4.483   -0.815)    5.241
   4.488   (   1.618   -2.802   -3.870)    5.044
   4.565   (  -1.384    2.398    3.633)    4.568
   4.624   (   3.310   -5.732   -1.689)    6.831
   4.730   (   0.357   -0.619   -1.381)    1.555
   4.785   (   1.863   -3.227    0.885)    3.830
   4.989   (  -0.734    1.272    1.422)    2.045
   5.019   (  -0.680    1.178   -0.189)    1.373
   5.089   (   2.937   -5.088    1.344)    6.026
   5.116   (  -3.127    5.417   -0.727)    6.297
   5.129   (  -1.033    1.788   -6.116)    6.455
   5.324   (   4.802   -8.318   -0.491)    9.617
   5.424   (   0.520   -0.900   -0.375)    1.105
   5.435   (   1.007   -1.744    0.149)    2.019
   5.615   (  -1.591    2.756    2.180)    3.858
   5.672   (   0.211   -0.366    3.262)    3.289
   5.735   (  -0.731    1.266   11.432)   11.525
   5.948   (   0.463   -0.803   -1.784)    2.011
   6.114   (   0.153   -0.265   -3.199)    3.213
   6.214   (  -4.183    7.246   -1.472)    8.495
   6.390   (   2.922   -5.060   -4.051)    7.110
   6.430   (   0.493   -0.854    0.250)    1.017
   7.542   (  -0.601    1.042   -6.328)    6.442
   7.572   (  -0.440    0.762    0.797)    1.187
   7.588   (   0.610   -1.056   -2.320)    2.621
   7.715   (  -2.549    4.416    4.563)    6.842
   7.850   (   1.621   -2.807    2.327)    3.991
   7.862   (   1.509   -2.613   -1.012)    3.183
   8.046   (   0.022   -0.038    1.306)    1.307
   8.047   (  -0.096    0.167    1.019)    1.037
=================== 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/15360
   10.0    152.096    152.096    143.816     -0.000     -0.000     -0.000 3/15360
   20.0     35.521     35.521     34.041     -0.000     -0.000     -0.000 3/15360
   30.0     14.466     14.466     13.946     -0.000     -0.000     -0.000 3/15360
   40.0      7.545      7.545      7.372     -0.000     -0.000     -0.000 3/15360
   50.0      4.696      4.696      4.652     -0.000     -0.000     -0.000 3/15360
   60.0      3.305      3.305      3.308     -0.000     -0.000     -0.000 3/15360
   70.0      2.526      2.526      2.547     -0.000     -0.000     -0.000 3/15360
   80.0      2.041      2.041      2.069     -0.000     -0.000     -0.000 3/15360
   90.0      1.714      1.714      1.744     -0.000     -0.000     -0.000 3/15360
  100.0      1.481      1.481      1.510     -0.000     -0.000     -0.000 3/15360
  110.0      1.306      1.306      1.333     -0.000     -0.000     -0.000 3/15360
  120.0      1.169      1.169      1.195     -0.000     -0.000     -0.000 3/15360
  130.0      1.060      1.060      1.085     -0.000     -0.000     -0.000 3/15360
  140.0      0.970      0.970      0.994     -0.000     -0.000     -0.000 3/15360
  150.0      0.895      0.895      0.917     -0.000     -0.000     -0.000 3/15360
  160.0      0.832      0.832      0.852     -0.000     -0.000     -0.000 3/15360
  170.0      0.777      0.777      0.796     -0.000     -0.000     -0.000 3/15360
  180.0      0.729      0.729      0.747     -0.000     -0.000     -0.000 3/15360
  190.0      0.687      0.687      0.704     -0.000     -0.000     -0.000 3/15360
  200.0      0.650      0.650      0.666     -0.000     -0.000     -0.000 3/15360
  210.0      0.617      0.617      0.632     -0.000     -0.000     -0.000 3/15360
  220.0      0.587      0.587      0.601     -0.000     -0.000     -0.000 3/15360
  230.0      0.560      0.560      0.574     -0.000     -0.000     -0.000 3/15360
  240.0      0.535      0.535      0.548     -0.000     -0.000     -0.000 3/15360
  250.0      0.513      0.513      0.525     -0.000     -0.000     -0.000 3/15360
  260.0      0.492      0.492      0.504     -0.000     -0.000     -0.000 3/15360
  270.0      0.473      0.473      0.485     -0.000     -0.000     -0.000 3/15360
  280.0      0.455      0.455      0.467     -0.000     -0.000     -0.000 3/15360
  290.0      0.439      0.439      0.450     -0.000     -0.000     -0.000 3/15360
  300.0      0.424      0.424      0.435     -0.000     -0.000     -0.000 3/15360
  310.0      0.410      0.410      0.420     -0.000     -0.000     -0.000 3/15360
  320.0      0.397      0.397      0.407     -0.000     -0.000     -0.000 3/15360
  330.0      0.384      0.384      0.394     -0.000     -0.000     -0.000 3/15360
  340.0      0.373      0.373      0.382     -0.000     -0.000     -0.000 3/15360
  350.0      0.362      0.362      0.371     -0.000     -0.000     -0.000 3/15360
  360.0      0.352      0.352      0.360     -0.000     -0.000     -0.000 3/15360
  370.0      0.342      0.342      0.350     -0.000     -0.000     -0.000 3/15360
  380.0      0.333      0.333      0.341     -0.000     -0.000     -0.000 3/15360
  390.0      0.324      0.324      0.332     -0.000     -0.000     -0.000 3/15360
  400.0      0.316      0.316      0.324     -0.000     -0.000     -0.000 3/15360
  410.0      0.308      0.308      0.316     -0.000     -0.000     -0.000 3/15360
  420.0      0.300      0.300      0.308     -0.000     -0.000     -0.000 3/15360
  430.0      0.293      0.293      0.301     -0.000     -0.000     -0.000 3/15360
  440.0      0.287      0.287      0.294     -0.000     -0.000     -0.000 3/15360
  450.0      0.280      0.280      0.287     -0.000     -0.000     -0.000 3/15360
  460.0      0.274      0.274      0.281     -0.000     -0.000     -0.000 3/15360
  470.0      0.268      0.268      0.275     -0.000     -0.000     -0.000 3/15360
  480.0      0.262      0.262      0.269     -0.000     -0.000     -0.000 3/15360
  490.0      0.257      0.257      0.263     -0.000     -0.000     -0.000 3/15360
  500.0      0.252      0.252      0.258     -0.000     -0.000     -0.000 3/15360
  510.0      0.247      0.247      0.253     -0.000     -0.000     -0.000 3/15360
  520.0      0.242      0.242      0.248     -0.000     -0.000     -0.000 3/15360
  530.0      0.237      0.237      0.243     -0.000     -0.000     -0.000 3/15360
  540.0      0.233      0.233      0.239     -0.000     -0.000     -0.000 3/15360
  550.0      0.229      0.229      0.234     -0.000     -0.000     -0.000 3/15360
  560.0      0.225      0.225      0.230     -0.000     -0.000     -0.000 3/15360
  570.0      0.221      0.221      0.226     -0.000     -0.000     -0.000 3/15360
  580.0      0.217      0.217      0.222     -0.000     -0.000     -0.000 3/15360
  590.0      0.213      0.213      0.218     -0.000     -0.000     -0.000 3/15360
  600.0      0.210      0.210      0.215     -0.000     -0.000     -0.000 3/15360
  610.0      0.206      0.206      0.211     -0.000     -0.000     -0.000 3/15360
  620.0      0.203      0.203      0.208     -0.000     -0.000     -0.000 3/15360
  630.0      0.200      0.200      0.204     -0.000     -0.000     -0.000 3/15360
  640.0      0.196      0.196      0.201     -0.000     -0.000     -0.000 3/15360
  650.0      0.193      0.193      0.198     -0.000     -0.000     -0.000 3/15360
  660.0      0.190      0.190      0.195     -0.000     -0.000     -0.000 3/15360
  670.0      0.188      0.188      0.192     -0.000     -0.000     -0.000 3/15360
  680.0      0.185      0.185      0.189     -0.000     -0.000     -0.000 3/15360
  690.0      0.182      0.182      0.187     -0.000     -0.000     -0.000 3/15360
  700.0      0.179      0.179      0.184     -0.000     -0.000     -0.000 3/15360
  710.0      0.177      0.177      0.181     -0.000     -0.000     -0.000 3/15360
  720.0      0.174      0.174      0.179     -0.000     -0.000     -0.000 3/15360
  730.0      0.172      0.172      0.176     -0.000     -0.000     -0.000 3/15360
  740.0      0.170      0.170      0.174     -0.000     -0.000     -0.000 3/15360
  750.0      0.167      0.167      0.172     -0.000     -0.000     -0.000 3/15360
  760.0      0.165      0.165      0.169     -0.000     -0.000     -0.000 3/15360
  770.0      0.163      0.163      0.167     -0.000     -0.000     -0.000 3/15360
  780.0      0.161      0.161      0.165     -0.000     -0.000     -0.000 3/15360
  790.0      0.159      0.159      0.163     -0.000     -0.000     -0.000 3/15360
  800.0      0.157      0.157      0.161     -0.000     -0.000     -0.000 3/15360
  810.0      0.155      0.155      0.159     -0.000     -0.000     -0.000 3/15360
  820.0      0.153      0.153      0.157     -0.000     -0.000     -0.000 3/15360
  830.0      0.151      0.151      0.155     -0.000     -0.000     -0.000 3/15360
  840.0      0.149      0.149      0.153     -0.000     -0.000     -0.000 3/15360
  850.0      0.148      0.148      0.151     -0.000     -0.000     -0.000 3/15360
  860.0      0.146      0.146      0.150     -0.000     -0.000     -0.000 3/15360
  870.0      0.144      0.144      0.148     -0.000     -0.000     -0.000 3/15360
  880.0      0.143      0.143      0.146     -0.000     -0.000     -0.000 3/15360
  890.0      0.141      0.141      0.144     -0.000     -0.000     -0.000 3/15360
  900.0      0.139      0.139      0.143     -0.000     -0.000     -0.000 3/15360
  910.0      0.138      0.138      0.141     -0.000     -0.000     -0.000 3/15360
  920.0      0.136      0.136      0.140     -0.000     -0.000     -0.000 3/15360
  930.0      0.135      0.135      0.138     -0.000     -0.000     -0.000 3/15360
  940.0      0.134      0.134      0.137     -0.000     -0.000     -0.000 3/15360
  950.0      0.132      0.132      0.135     -0.000     -0.000     -0.000 3/15360
  960.0      0.131      0.131      0.134     -0.000     -0.000     -0.000 3/15360
  970.0      0.129      0.129      0.133     -0.000     -0.000     -0.000 3/15360
  980.0      0.128      0.128      0.131     -0.000     -0.000     -0.000 3/15360
  990.0      0.127      0.127      0.130     -0.000     -0.000     -0.000 3/15360
 1000.0      0.125      0.125      0.129     -0.000     -0.000     -0.000 3/15360

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

