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

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

-------------------------[time 2026-01-07 18:30:19]-------------------------
Compiled with OpenMP support (max 32 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: [4 4 2]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
Number of symmetry operations in supercell: 768
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.827787792793303    0.000000000000000    0.000000000000000
  b   -1.913893896396651    3.314961468854965    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.322216750000000
Atomic positions (fractional):
   *1 Si  0.33333333333333  0.66666666666667  0.93710958016861  28.085
    2 Si  0.66666666666667  0.33333333333333  0.43710958016861  28.085
    3 Si  0.33333333333333  0.66666666666667  0.56289041983139  28.085
    4 Si  0.66666666666667  0.33333333333333  0.06289041983139  28.085
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.827787792793303    0.000000000000000    0.000000000000000
  b   -1.913893896396651    3.314961468854965    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.322216750000000
Atomic positions (fractional):
   *1 Si  0.33333333333333  0.66666666666667  0.93710958016861  28.085 > 1
    2 Si  0.66666666666667  0.33333333333333  0.43710958016861  28.085 > 2
    3 Si  0.33333333333333  0.66666666666667  0.56289041983139  28.085 > 3
    4 Si  0.66666666666667  0.33333333333333  0.06289041983139  28.085 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   15.311151171173211    0.000000000000000    0.000000000000000
  b   -7.655575585586604   13.259845875419860    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.644433500000000
Atomic positions (fractional):
   *1 Si  0.08333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    2 Si  0.33333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    3 Si  0.58333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    4 Si  0.83333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    5 Si  0.08333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    6 Si  0.33333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    7 Si  0.58333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    8 Si  0.83333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    9 Si  0.08333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   10 Si  0.33333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   11 Si  0.58333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   12 Si  0.83333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   13 Si  0.08333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   14 Si  0.33333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   15 Si  0.58333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   16 Si  0.83333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   17 Si  0.08333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   18 Si  0.33333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   19 Si  0.58333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   20 Si  0.83333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   21 Si  0.08333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   22 Si  0.33333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   23 Si  0.58333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   24 Si  0.83333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   25 Si  0.08333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   26 Si  0.33333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   27 Si  0.58333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   28 Si  0.83333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   29 Si  0.08333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   30 Si  0.33333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   31 Si  0.58333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   32 Si  0.83333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   33 Si  0.16666666666667  0.08333333333333  0.21855479008430  28.085 > 2
   34 Si  0.41666666666667  0.08333333333333  0.21855479008430  28.085 > 2
   35 Si  0.66666666666667  0.08333333333333  0.21855479008430  28.085 > 2
   36 Si  0.91666666666667  0.08333333333333  0.21855479008430  28.085 > 2
   37 Si  0.16666666666667  0.33333333333333  0.21855479008430  28.085 > 2
   38 Si  0.41666666666667  0.33333333333333  0.21855479008430  28.085 > 2
   39 Si  0.66666666666667  0.33333333333333  0.21855479008430  28.085 > 2
   40 Si  0.91666666666667  0.33333333333333  0.21855479008430  28.085 > 2
   41 Si  0.16666666666667  0.58333333333333  0.21855479008430  28.085 > 2
   42 Si  0.41666666666667  0.58333333333333  0.21855479008430  28.085 > 2
   43 Si  0.66666666666667  0.58333333333333  0.21855479008430  28.085 > 2
   44 Si  0.91666666666667  0.58333333333333  0.21855479008430  28.085 > 2
   45 Si  0.16666666666667  0.83333333333333  0.21855479008430  28.085 > 2
   46 Si  0.41666666666667  0.83333333333333  0.21855479008430  28.085 > 2
   47 Si  0.66666666666667  0.83333333333333  0.21855479008430  28.085 > 2
   48 Si  0.91666666666667  0.83333333333333  0.21855479008430  28.085 > 2
   49 Si  0.16666666666667  0.08333333333333  0.71855479008430  28.085 > 2
   50 Si  0.41666666666667  0.08333333333333  0.71855479008430  28.085 > 2
   51 Si  0.66666666666667  0.08333333333333  0.71855479008430  28.085 > 2
   52 Si  0.91666666666667  0.08333333333333  0.71855479008430  28.085 > 2
   53 Si  0.16666666666667  0.33333333333333  0.71855479008430  28.085 > 2
   54 Si  0.41666666666667  0.33333333333333  0.71855479008430  28.085 > 2
   55 Si  0.66666666666667  0.33333333333333  0.71855479008430  28.085 > 2
   56 Si  0.91666666666667  0.33333333333333  0.71855479008430  28.085 > 2
   57 Si  0.16666666666667  0.58333333333333  0.71855479008430  28.085 > 2
   58 Si  0.41666666666667  0.58333333333333  0.71855479008430  28.085 > 2
   59 Si  0.66666666666667  0.58333333333333  0.71855479008430  28.085 > 2
   60 Si  0.91666666666667  0.58333333333333  0.71855479008430  28.085 > 2
   61 Si  0.16666666666667  0.83333333333333  0.71855479008430  28.085 > 2
   62 Si  0.41666666666667  0.83333333333333  0.71855479008430  28.085 > 2
   63 Si  0.66666666666667  0.83333333333333  0.71855479008430  28.085 > 2
   64 Si  0.91666666666667  0.83333333333333  0.71855479008430  28.085 > 2
   65 Si  0.08333333333333  0.16666666666667  0.28144520991570  28.085 > 3
   66 Si  0.33333333333333  0.16666666666667  0.28144520991570  28.085 > 3
   67 Si  0.58333333333333  0.16666666666667  0.28144520991570  28.085 > 3
   68 Si  0.83333333333333  0.16666666666667  0.28144520991570  28.085 > 3
   69 Si  0.08333333333333  0.41666666666667  0.28144520991570  28.085 > 3
   70 Si  0.33333333333333  0.41666666666667  0.28144520991570  28.085 > 3
   71 Si  0.58333333333333  0.41666666666667  0.28144520991570  28.085 > 3
   72 Si  0.83333333333333  0.41666666666667  0.28144520991570  28.085 > 3
   73 Si  0.08333333333333  0.66666666666667  0.28144520991570  28.085 > 3
   74 Si  0.33333333333333  0.66666666666667  0.28144520991570  28.085 > 3
   75 Si  0.58333333333333  0.66666666666667  0.28144520991570  28.085 > 3
   76 Si  0.83333333333333  0.66666666666667  0.28144520991570  28.085 > 3
   77 Si  0.08333333333333  0.91666666666667  0.28144520991570  28.085 > 3
   78 Si  0.33333333333333  0.91666666666667  0.28144520991570  28.085 > 3
   79 Si  0.58333333333333  0.91666666666667  0.28144520991570  28.085 > 3
   80 Si  0.83333333333333  0.91666666666667  0.28144520991570  28.085 > 3
   81 Si  0.08333333333333  0.16666666666667  0.78144520991570  28.085 > 3
   82 Si  0.33333333333333  0.16666666666667  0.78144520991570  28.085 > 3
   83 Si  0.58333333333333  0.16666666666667  0.78144520991570  28.085 > 3
   84 Si  0.83333333333333  0.16666666666667  0.78144520991570  28.085 > 3
   85 Si  0.08333333333333  0.41666666666667  0.78144520991570  28.085 > 3
   86 Si  0.33333333333333  0.41666666666667  0.78144520991570  28.085 > 3
   87 Si  0.58333333333333  0.41666666666667  0.78144520991570  28.085 > 3
   88 Si  0.83333333333333  0.41666666666667  0.78144520991570  28.085 > 3
   89 Si  0.08333333333333  0.66666666666667  0.78144520991570  28.085 > 3
   90 Si  0.33333333333333  0.66666666666667  0.78144520991570  28.085 > 3
   91 Si  0.58333333333333  0.66666666666667  0.78144520991570  28.085 > 3
   92 Si  0.83333333333333  0.66666666666667  0.78144520991570  28.085 > 3
   93 Si  0.08333333333333  0.91666666666667  0.78144520991570  28.085 > 3
   94 Si  0.33333333333333  0.91666666666667  0.78144520991570  28.085 > 3
   95 Si  0.58333333333333  0.91666666666667  0.78144520991570  28.085 > 3
   96 Si  0.83333333333333  0.91666666666667  0.78144520991570  28.085 > 3
   97 Si  0.16666666666667  0.08333333333333  0.03144520991570  28.085 > 4
   98 Si  0.41666666666667  0.08333333333333  0.03144520991570  28.085 > 4
   99 Si  0.66666666666667  0.08333333333333  0.03144520991570  28.085 > 4
  100 Si  0.91666666666667  0.08333333333333  0.03144520991570  28.085 > 4
  101 Si  0.16666666666667  0.33333333333333  0.03144520991570  28.085 > 4
  102 Si  0.41666666666667  0.33333333333333  0.03144520991570  28.085 > 4
  103 Si  0.66666666666667  0.33333333333333  0.03144520991570  28.085 > 4
  104 Si  0.91666666666667  0.33333333333333  0.03144520991570  28.085 > 4
  105 Si  0.16666666666667  0.58333333333333  0.03144520991570  28.085 > 4
  106 Si  0.41666666666667  0.58333333333333  0.03144520991570  28.085 > 4
  107 Si  0.66666666666667  0.58333333333333  0.03144520991570  28.085 > 4
  108 Si  0.91666666666667  0.58333333333333  0.03144520991570  28.085 > 4
  109 Si  0.16666666666667  0.83333333333333  0.03144520991570  28.085 > 4
  110 Si  0.41666666666667  0.83333333333333  0.03144520991570  28.085 > 4
  111 Si  0.66666666666667  0.83333333333333  0.03144520991570  28.085 > 4
  112 Si  0.91666666666667  0.83333333333333  0.03144520991570  28.085 > 4
  113 Si  0.16666666666667  0.08333333333333  0.53144520991570  28.085 > 4
  114 Si  0.41666666666667  0.08333333333333  0.53144520991570  28.085 > 4
  115 Si  0.66666666666667  0.08333333333333  0.53144520991570  28.085 > 4
  116 Si  0.91666666666667  0.08333333333333  0.53144520991570  28.085 > 4
  117 Si  0.16666666666667  0.33333333333333  0.53144520991570  28.085 > 4
  118 Si  0.41666666666667  0.33333333333333  0.53144520991570  28.085 > 4
  119 Si  0.66666666666667  0.33333333333333  0.53144520991570  28.085 > 4
  120 Si  0.91666666666667  0.33333333333333  0.53144520991570  28.085 > 4
  121 Si  0.16666666666667  0.58333333333333  0.53144520991570  28.085 > 4
  122 Si  0.41666666666667  0.58333333333333  0.53144520991570  28.085 > 4
  123 Si  0.66666666666667  0.58333333333333  0.53144520991570  28.085 > 4
  124 Si  0.91666666666667  0.58333333333333  0.53144520991570  28.085 > 4
  125 Si  0.16666666666667  0.83333333333333  0.53144520991570  28.085 > 4
  126 Si  0.41666666666667  0.83333333333333  0.53144520991570  28.085 > 4
  127 Si  0.66666666666667  0.83333333333333  0.53144520991570  28.085 > 4
  128 Si  0.91666666666667  0.83333333333333  0.53144520991570  28.085 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           16.0455516    0.0000000    0.0000000
            0.0000000   16.0455516    0.0000000
            0.0000000    0.0000000   15.4688150
-------------------------- Born effective charges --------------------------
    1 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    2 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    3 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    4 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
----------------------------------------------------------------------------
Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
-------------------------------- Symfc start -------------------------------
Symfc version 1.5.4 (https://github.com/symfc/symfc)
Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024)
Computing [2] order force constants.
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 134
Number of blocks in projector: 102
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: 63
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 40
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 31
Use standard eigh solver.
Tree of FC basis block matrices:
- (134, 132), data: False
|-- (31, 30), data: True
|-- (40, 40), data: True
|-- (63, 62), data: True
-----
Solver_atoms: 1 -- 128 / 128
Time (Solver_compr_matrix_reshape): 0.012
Solver_block: 80 / 80
 - Time: 0.071
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.083
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 134
Number of blocks in projector: 102
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: 63
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 40
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 31
Use standard eigh solver.
Tree of FC basis block matrices:
- (134, 132), data: False
|-- (31, 30), data: True
|-- (40, 40), data: True
|-- (63, 62), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) 
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-07 18:30:23]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 18:30:24]-------------------------
Compiled with OpenMP support (max 32 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): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.827787792793303    0.000000000000000    0.000000000000000
  b   -1.913893896396651    3.314961468854965    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.322216750000000
Atomic positions (fractional):
    1 Si  0.33333333333333  0.66666666666667  0.93710958016861  28.085
    2 Si  0.66666666666667  0.33333333333333  0.43710958016861  28.085
    3 Si  0.33333333333333  0.66666666666667  0.56289041983139  28.085
    4 Si  0.66666666666667  0.33333333333333  0.06289041983139  28.085
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.311151171173211    0.000000000000000    0.000000000000000
  b   -7.655575585586604   13.259845875419860    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.644433500000000
Atomic positions (fractional):
    1 Si  0.08333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    2 Si  0.33333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    3 Si  0.58333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    4 Si  0.83333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    5 Si  0.08333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    6 Si  0.33333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    7 Si  0.58333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    8 Si  0.83333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    9 Si  0.08333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   10 Si  0.33333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   11 Si  0.58333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   12 Si  0.83333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   13 Si  0.08333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   14 Si  0.33333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   15 Si  0.58333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   16 Si  0.83333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   17 Si  0.08333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   18 Si  0.33333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   19 Si  0.58333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   20 Si  0.83333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   21 Si  0.08333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   22 Si  0.33333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   23 Si  0.58333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   24 Si  0.83333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   25 Si  0.08333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   26 Si  0.33333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   27 Si  0.58333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   28 Si  0.83333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   29 Si  0.08333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   30 Si  0.33333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   31 Si  0.58333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   32 Si  0.83333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   33 Si  0.16666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   34 Si  0.41666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   35 Si  0.66666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   36 Si  0.91666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   37 Si  0.16666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   38 Si  0.41666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   39 Si  0.66666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   40 Si  0.91666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   41 Si  0.16666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   42 Si  0.41666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   43 Si  0.66666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   44 Si  0.91666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   45 Si  0.16666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   46 Si  0.41666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   47 Si  0.66666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   48 Si  0.91666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   49 Si  0.16666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   50 Si  0.41666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   51 Si  0.66666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   52 Si  0.91666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   53 Si  0.16666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   54 Si  0.41666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   55 Si  0.66666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   56 Si  0.91666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   57 Si  0.16666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   58 Si  0.41666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   59 Si  0.66666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   60 Si  0.91666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   61 Si  0.16666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   62 Si  0.41666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   63 Si  0.66666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   64 Si  0.91666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   65 Si  0.08333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   66 Si  0.33333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   67 Si  0.58333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   68 Si  0.83333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   69 Si  0.08333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   70 Si  0.33333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   71 Si  0.58333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   72 Si  0.83333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   73 Si  0.08333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   74 Si  0.33333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   75 Si  0.58333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   76 Si  0.83333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   77 Si  0.08333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   78 Si  0.33333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   79 Si  0.58333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   80 Si  0.83333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   81 Si  0.08333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   82 Si  0.33333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   83 Si  0.58333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   84 Si  0.83333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   85 Si  0.08333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   86 Si  0.33333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   87 Si  0.58333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   88 Si  0.83333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   89 Si  0.08333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   90 Si  0.33333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   91 Si  0.58333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   92 Si  0.83333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   93 Si  0.08333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   94 Si  0.33333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   95 Si  0.58333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   96 Si  0.83333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   97 Si  0.16666666666667  0.08333333333333  0.03144520991570  28.085 > 97
   98 Si  0.41666666666667  0.08333333333333  0.03144520991570  28.085 > 97
   99 Si  0.66666666666667  0.08333333333333  0.03144520991570  28.085 > 97
  100 Si  0.91666666666667  0.08333333333333  0.03144520991570  28.085 > 97
  101 Si  0.16666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  102 Si  0.41666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  103 Si  0.66666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  104 Si  0.91666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  105 Si  0.16666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  106 Si  0.41666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  107 Si  0.66666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  108 Si  0.91666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  109 Si  0.16666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  110 Si  0.41666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  111 Si  0.66666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  112 Si  0.91666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  113 Si  0.16666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  114 Si  0.41666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  115 Si  0.66666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  116 Si  0.91666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  117 Si  0.16666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  118 Si  0.41666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  119 Si  0.66666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  120 Si  0.91666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  121 Si  0.16666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  122 Si  0.41666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  123 Si  0.66666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  124 Si  0.91666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  125 Si  0.16666666666667  0.83333333333333  0.53144520991570  28.085 > 97
  126 Si  0.41666666666667  0.83333333333333  0.53144520991570  28.085 > 97
  127 Si  0.66666666666667  0.83333333333333  0.53144520991570  28.085 > 97
  128 Si  0.91666666666667  0.83333333333333  0.53144520991570  28.085 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           16.0455516    0.0000000    0.0000000
            0.0000000   16.0455516    0.0000000
            0.0000000    0.0000000   15.4688150
-------------------------- Born effective charges --------------------------
    1 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    2 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    3 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    4 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
----------------------------------------------------------------------------
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]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000058 (zzz) 0.00000058 (zzz) 0.00000058 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) 
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-07 18:30:27]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 18:30:28]-------------------------
Compiled with OpenMP support (max 32 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): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.827787792793303    0.000000000000000    0.000000000000000
  b   -1.913893896396651    3.314961468854965    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.322216750000000
Atomic positions (fractional):
    1 Si  0.33333333333333  0.66666666666667  0.93710958016861  28.085
    2 Si  0.66666666666667  0.33333333333333  0.43710958016861  28.085
    3 Si  0.33333333333333  0.66666666666667  0.56289041983139  28.085
    4 Si  0.66666666666667  0.33333333333333  0.06289041983139  28.085
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.311151171173211    0.000000000000000    0.000000000000000
  b   -7.655575585586604   13.259845875419860    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.644433500000000
Atomic positions (fractional):
    1 Si  0.08333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    2 Si  0.33333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    3 Si  0.58333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    4 Si  0.83333333333333  0.16666666666667  0.46855479008430  28.085 > 1
    5 Si  0.08333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    6 Si  0.33333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    7 Si  0.58333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    8 Si  0.83333333333333  0.41666666666667  0.46855479008430  28.085 > 1
    9 Si  0.08333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   10 Si  0.33333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   11 Si  0.58333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   12 Si  0.83333333333333  0.66666666666667  0.46855479008430  28.085 > 1
   13 Si  0.08333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   14 Si  0.33333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   15 Si  0.58333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   16 Si  0.83333333333333  0.91666666666667  0.46855479008430  28.085 > 1
   17 Si  0.08333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   18 Si  0.33333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   19 Si  0.58333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   20 Si  0.83333333333333  0.16666666666667  0.96855479008430  28.085 > 1
   21 Si  0.08333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   22 Si  0.33333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   23 Si  0.58333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   24 Si  0.83333333333333  0.41666666666667  0.96855479008430  28.085 > 1
   25 Si  0.08333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   26 Si  0.33333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   27 Si  0.58333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   28 Si  0.83333333333333  0.66666666666667  0.96855479008430  28.085 > 1
   29 Si  0.08333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   30 Si  0.33333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   31 Si  0.58333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   32 Si  0.83333333333333  0.91666666666667  0.96855479008430  28.085 > 1
   33 Si  0.16666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   34 Si  0.41666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   35 Si  0.66666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   36 Si  0.91666666666667  0.08333333333333  0.21855479008430  28.085 > 33
   37 Si  0.16666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   38 Si  0.41666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   39 Si  0.66666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   40 Si  0.91666666666667  0.33333333333333  0.21855479008430  28.085 > 33
   41 Si  0.16666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   42 Si  0.41666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   43 Si  0.66666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   44 Si  0.91666666666667  0.58333333333333  0.21855479008430  28.085 > 33
   45 Si  0.16666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   46 Si  0.41666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   47 Si  0.66666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   48 Si  0.91666666666667  0.83333333333333  0.21855479008430  28.085 > 33
   49 Si  0.16666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   50 Si  0.41666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   51 Si  0.66666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   52 Si  0.91666666666667  0.08333333333333  0.71855479008430  28.085 > 33
   53 Si  0.16666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   54 Si  0.41666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   55 Si  0.66666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   56 Si  0.91666666666667  0.33333333333333  0.71855479008430  28.085 > 33
   57 Si  0.16666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   58 Si  0.41666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   59 Si  0.66666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   60 Si  0.91666666666667  0.58333333333333  0.71855479008430  28.085 > 33
   61 Si  0.16666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   62 Si  0.41666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   63 Si  0.66666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   64 Si  0.91666666666667  0.83333333333333  0.71855479008430  28.085 > 33
   65 Si  0.08333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   66 Si  0.33333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   67 Si  0.58333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   68 Si  0.83333333333333  0.16666666666667  0.28144520991570  28.085 > 65
   69 Si  0.08333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   70 Si  0.33333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   71 Si  0.58333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   72 Si  0.83333333333333  0.41666666666667  0.28144520991570  28.085 > 65
   73 Si  0.08333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   74 Si  0.33333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   75 Si  0.58333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   76 Si  0.83333333333333  0.66666666666667  0.28144520991570  28.085 > 65
   77 Si  0.08333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   78 Si  0.33333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   79 Si  0.58333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   80 Si  0.83333333333333  0.91666666666667  0.28144520991570  28.085 > 65
   81 Si  0.08333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   82 Si  0.33333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   83 Si  0.58333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   84 Si  0.83333333333333  0.16666666666667  0.78144520991570  28.085 > 65
   85 Si  0.08333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   86 Si  0.33333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   87 Si  0.58333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   88 Si  0.83333333333333  0.41666666666667  0.78144520991570  28.085 > 65
   89 Si  0.08333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   90 Si  0.33333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   91 Si  0.58333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   92 Si  0.83333333333333  0.66666666666667  0.78144520991570  28.085 > 65
   93 Si  0.08333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   94 Si  0.33333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   95 Si  0.58333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   96 Si  0.83333333333333  0.91666666666667  0.78144520991570  28.085 > 65
   97 Si  0.16666666666667  0.08333333333333  0.03144520991570  28.085 > 97
   98 Si  0.41666666666667  0.08333333333333  0.03144520991570  28.085 > 97
   99 Si  0.66666666666667  0.08333333333333  0.03144520991570  28.085 > 97
  100 Si  0.91666666666667  0.08333333333333  0.03144520991570  28.085 > 97
  101 Si  0.16666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  102 Si  0.41666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  103 Si  0.66666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  104 Si  0.91666666666667  0.33333333333333  0.03144520991570  28.085 > 97
  105 Si  0.16666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  106 Si  0.41666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  107 Si  0.66666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  108 Si  0.91666666666667  0.58333333333333  0.03144520991570  28.085 > 97
  109 Si  0.16666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  110 Si  0.41666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  111 Si  0.66666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  112 Si  0.91666666666667  0.83333333333333  0.03144520991570  28.085 > 97
  113 Si  0.16666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  114 Si  0.41666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  115 Si  0.66666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  116 Si  0.91666666666667  0.08333333333333  0.53144520991570  28.085 > 97
  117 Si  0.16666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  118 Si  0.41666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  119 Si  0.66666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  120 Si  0.91666666666667  0.33333333333333  0.53144520991570  28.085 > 97
  121 Si  0.16666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  122 Si  0.41666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  123 Si  0.66666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  124 Si  0.91666666666667  0.58333333333333  0.53144520991570  28.085 > 97
  125 Si  0.16666666666667  0.83333333333333  0.53144520991570  28.085 > 97
  126 Si  0.41666666666667  0.83333333333333  0.53144520991570  28.085 > 97
  127 Si  0.66666666666667  0.83333333333333  0.53144520991570  28.085 > 97
  128 Si  0.91666666666667  0.83333333333333  0.53144520991570  28.085 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           16.0455516    0.0000000    0.0000000
            0.0000000   16.0455516    0.0000000
            0.0000000    0.0000000   15.4688150
-------------------------- Born effective charges --------------------------
    1 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    2 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    3 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    4 Si    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000058 (zzz) 0.00000058 (zzz) 0.00000058 (zzz)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- 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 ... [ 15 15 8 ]
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.96, Number of G-points: 313, Lambda: 0.40
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/135) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 135
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
   2.781   (   0.000    0.000    0.000)    0.000
   2.781   (   0.000    0.000    0.000)    0.000
  11.015   (   0.000    0.000    0.000)    0.000
  12.359   (   0.000    0.000    0.000)    0.000
  14.664   (   0.000    0.000    0.000)    0.000
  14.664   (   0.000    0.000    0.000)    0.000
  15.079   (   0.000    0.000    0.000)    0.000
  15.098   (   0.000    0.000    0.000)    0.000
  15.098   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/135) =======================
q-point: ( 0.07  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.849   (  36.451   21.045    0.000)   42.090
   0.871   (  36.748   21.216    0.000)   42.433
   1.737   (  74.001   42.724    0.000)   85.449
   2.899   (   9.744    5.626    0.000)   11.252
   3.250   (  36.761   21.224    0.000)   42.448
  10.856   ( -13.230   -7.638    0.000)   15.276
  12.293   (  -5.520   -3.187    0.000)    6.374
  14.597   (  -5.628   -3.249    0.000)    6.499
  14.654   (  -1.102   -0.636    0.000)    1.272
  15.025   (  -4.415   -2.549    0.000)    5.098
  15.038   (  -5.000   -2.887    0.000)    5.773
  15.039   (  -4.989   -2.880    0.000)    5.760
======================= Grid point 2 (3/135) =======================
q-point: ( 0.13  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.685   (  32.825   18.951    0.000)   37.903
   1.693   (  36.206   20.903    0.000)   41.807
   3.196   (  14.753    8.518    0.000)   17.035
   3.418   (  70.364   40.625    0.000)   81.249
   4.282   (  48.188   27.822    0.000)   55.643
  10.445   ( -20.546  -11.862    0.000)   23.724
  12.112   (  -9.906   -5.719    0.000)   11.438
  14.417   (  -9.460   -5.461    0.000)   10.923
  14.589   (  -5.209   -3.007    0.000)    6.014
  14.869   (  -9.576   -5.529    0.000)   11.058
  14.878   (  -8.339   -4.814    0.000)    9.629
  14.904   (  -5.068   -2.926    0.000)    5.852
======================= Grid point 3 (4/135) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.367   (  25.240   14.572    0.000)   29.144
   2.526   (  35.310   20.386    0.000)   40.772
   3.530   (  12.901    7.448    0.000)   14.897
   4.984   (  64.021   36.963    0.000)   73.925
   5.335   (  39.530   22.823    0.000)   45.646
  10.003   ( -14.718   -8.497    0.000)   16.995
  11.843   ( -13.095   -7.560    0.000)   15.121
  14.182   ( -10.121   -5.844    0.000)   11.687
  14.406   (  -9.846   -5.685    0.000)   11.369
  14.598   ( -13.752   -7.939    0.000)   15.879
  14.674   (  -8.676   -5.009    0.000)   10.018
  14.812   (  -3.040   -1.755    0.000)    3.511
======================= Grid point 4 (5/135) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.833   (  14.633    8.449    0.000)   16.897
   3.320   (  32.692   18.875    0.000)   37.750
   3.749   (   5.236    3.023    0.000)    6.046
   5.996   (  15.950    9.209    0.000)   18.418
   6.370   (  55.007   31.758    0.000)   63.517
   9.898   (   7.650    4.417    0.000)    8.834
  11.509   ( -15.588   -9.000    0.000)   17.999
  13.972   (  -7.527   -4.346    0.000)    8.692
  14.172   (  -9.581   -5.532    0.000)   11.063
  14.233   ( -17.573  -10.146    0.000)   20.292
  14.500   (  -5.870   -3.389    0.000)    6.778
  14.748   (  -2.647   -1.529    0.000)    3.057
======================= Grid point 5 (6/135) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.048   (   4.243    2.450    0.000)    4.900
   3.762   (  -3.833   -2.213    0.000)    4.426
   4.025   (  27.576   15.921    0.000)   31.842
   6.115   (  -2.991   -1.727    0.000)    3.454
   7.522   (  43.937   25.367    0.000)   50.734
  10.340   (  27.939   16.130    0.000)   32.261
  11.125   ( -17.297   -9.986    0.000)   19.973
  13.785   ( -20.977  -12.111    0.000)   24.223
  13.844   (  -3.502   -2.022    0.000)    4.044
  13.985   (  -6.257   -3.613    0.000)    7.225
  14.412   (  -1.799   -1.039    0.000)    2.078
  14.687   (  -2.528   -1.459    0.000)    2.919
======================= Grid point 6 (7/135) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.066   (  -1.684   -0.973    0.000)    1.945
   3.609   (  -8.120   -4.688    0.000)    9.376
   4.578   (  19.573   11.300    0.000)   22.601
   5.980   (  -6.829   -3.943    0.000)    7.886
   8.401   (  31.415   18.138    0.000)   36.275
  10.721   ( -16.863   -9.736    0.000)   19.471
  11.071   (  32.866   18.975    0.000)   37.951
  13.264   ( -23.723  -13.696    0.000)   27.392
  13.800   (  -0.609   -0.351    0.000)    0.703
  13.883   (  -2.638   -1.523    0.000)    3.046
  14.404   (   0.747    0.431    0.000)    0.862
  14.634   (  -1.945   -1.123    0.000)    2.246
======================= Grid point 7 (8/135) =======================
q-point: ( 0.47  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.018   (  -1.526   -0.881    0.000)    1.762
   3.451   (  -4.160   -2.402    0.000)    4.804
   4.904   (   7.703    4.447    0.000)    8.895
   5.852   (  -3.385   -1.955    0.000)    3.909
   8.945   (  13.744    7.935    0.000)   15.870
  10.394   (  -9.259   -5.346    0.000)   10.691
  11.777   (  24.823   14.332    0.000)   28.663
  12.714   ( -21.145  -12.208    0.000)   24.416
  13.799   (   0.134    0.078    0.000)    0.155
  13.850   (  -0.545   -0.315    0.000)    0.630
  14.425   (   0.659    0.381    0.000)    0.761
  14.602   (  -0.757   -0.437    0.000)    0.874
======================= Grid point 16 (9/135) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.463   (  20.742   35.926    0.000)   41.484
   1.567   (  23.457   40.630    0.000)   46.915
   2.929   (  39.388   68.222    0.000)   78.776
   3.153   (  10.555   18.281    0.000)   21.109
   3.949   (  25.714   44.538    0.000)   51.428
  10.571   ( -11.280  -19.538    0.000)   22.560
  12.170   (  -5.132   -8.889    0.000)   10.264
  14.465   (  -5.376   -9.311    0.000)   10.751
  14.632   (  -1.206   -2.090    0.000)    2.413
  14.917   (  -5.258   -9.108    0.000)   10.517
  14.929   (  -3.758   -6.509    0.000)    7.516
  14.935   (  -4.098   -7.099    0.000)    8.197
======================= Grid point 17 (10/135) =======================
q-point: ( 0.13  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.222   (  27.235   30.251    0.000)   40.705
   2.352   (  16.652   42.498    0.000)   45.644
   3.558   (   8.347   25.967    0.000)   27.275
   4.340   (  53.134   48.638    0.000)   72.034
   4.902   (  36.128   31.400    0.000)   47.866
  10.141   ( -14.354  -16.767    0.000)   22.072
  11.945   (  -9.142  -10.703    0.000)   14.076
  14.249   (  -7.373  -10.206    0.000)   12.591
  14.538   (  -7.809   -2.407    0.000)    8.171
  14.678   (  -9.258  -12.983    0.000)   15.946
  14.768   (  -6.869   -6.309    0.000)    9.326
  14.798   (  -1.953   -7.561    0.000)    7.809
======================= Grid point 18 (11/135) =======================
q-point: ( 0.20  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.958   (   4.200   40.582    0.000)   40.798
   2.997   (  28.843   25.677    0.000)   38.617
   3.906   (  -1.131   28.418    0.000)   28.441
   5.711   (  55.456   34.181    0.000)   65.143
   5.712   (  27.644   13.429    0.000)   30.733
   9.852   (  -2.460   -5.253    0.000)    5.801
  11.643   ( -12.160  -12.294    0.000)   17.292
  14.031   (  -6.449   -8.067    0.000)   10.328
  14.325   ( -12.170   -2.971    0.000)   12.528
  14.341   ( -12.026  -16.882    0.000)   20.727
  14.590   (  -7.096   -3.733    0.000)    8.018
  14.715   (   0.654   -7.829    0.000)    7.856
======================= Grid point 19 (12/135) =======================
q-point: ( 0.27  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.302   (  -8.296   35.166    0.000)   36.131
   3.722   (  26.711   20.825    0.000)   33.870
   4.065   ( -12.455   26.802    0.000)   29.554
   6.074   (   7.331   -0.581    0.000)    7.354
   6.930   (  50.641   23.445    0.000)   55.804
  10.023   (  20.386    7.410    0.000)   21.691
  11.278   ( -14.360  -13.944    0.000)   20.016
  13.879   (  -3.730   -4.083    0.000)    5.530
  13.923   ( -14.550  -20.053    0.000)   24.775
  14.094   ( -10.205   -2.869    0.000)   10.601
  14.465   (  -4.490   -0.510    0.000)    4.518
  14.651   (   1.035   -7.919    0.000)    7.986
======================= Grid point 20 (13/135) =======================
q-point: ( 0.33  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.399   ( -17.041   29.909    0.000)   34.423
   3.999   ( -21.063   24.597    0.000)   32.383
   4.354   (  21.163   16.627    0.000)   26.913
   6.047   (  -4.620   -4.127    0.000)    6.195
   7.928   (  41.349   15.097    0.000)   44.019
  10.621   (  34.394   10.936    0.000)   36.090
  10.872   ( -15.053  -15.261    0.000)   21.435
  13.434   ( -17.127  -22.282    0.000)   28.103
  13.815   (  -1.308   -0.571    0.000)    1.427
  13.930   (  -5.466   -2.305    0.000)    5.932
  14.424   (  -1.366    1.921    0.000)    2.357
  14.592   (   1.394   -7.735    0.000)    7.859
======================= Grid point 21 (14/135) =======================
q-point: ( 0.40  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.358   ( -19.435   27.819    0.000)   33.935
   3.851   ( -23.856   27.244    0.000)   36.212
   4.797   (  12.912   10.454    0.000)   16.614
   5.909   (  -5.686   -2.873    0.000)    6.370
   8.655   (  27.991    7.397    0.000)   28.952
  10.475   ( -11.670  -14.726    0.000)   18.789
  11.345   (  35.579    7.734    0.000)   36.410
  12.888   ( -18.582  -22.935    0.000)   29.518
  13.801   (  -0.664    0.069    0.000)    0.668
  13.859   (  -1.075   -0.838    0.000)    1.363
  14.438   (  -0.253    2.650    0.000)    2.662
  14.548   (   2.271   -7.198    0.000)    7.547
======================= Grid point 22 (15/135) =======================
q-point: ( 0.47  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.323   ( -16.761   29.031    0.000)   33.522
   3.772   ( -18.439   31.938    0.000)   36.879
   4.972   (  -1.331    2.305    0.000)    2.662
   5.836   (  -0.247    0.428    0.000)    0.494
   8.966   (   3.434   -5.948    0.000)    6.868
  10.268   (   3.943   -6.829    0.000)    7.885
  11.827   (   6.810  -11.796    0.000)   13.621
  12.482   (   4.540   -7.863    0.000)    9.079
  13.793   (   0.296   -0.512    0.000)    0.592
  13.851   (  -0.202    0.350    0.000)    0.404
  14.451   (  -1.201    2.081    0.000)    2.403
  14.531   (   3.677   -6.368    0.000)    7.353
======================= Grid point 33 (16/135) =======================
q-point: ( 0.13  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.869   (  19.211   33.274    0.000)   38.422
   3.205   (  22.409   38.813    0.000)   44.818
   4.126   (  15.853   27.459    0.000)   31.707
   5.336   (  29.018   50.261    0.000)   58.036
   5.479   (  14.386   24.917    0.000)   28.771
   9.852   (  -5.615   -9.726    0.000)   11.231
  11.687   (  -8.596  -14.889    0.000)   17.192
  14.045   (  -5.409   -9.369    0.000)   10.818
  14.356   ( -10.802  -18.709    0.000)   21.604
  14.462   (  -3.659   -6.337    0.000)    7.317
  14.626   (  -4.654   -8.061    0.000)    9.308
  14.636   (  -3.707   -6.420    0.000)    7.413
======================= Grid point 34 (17/135) =======================
q-point: ( 0.20  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.534   (  21.620   27.234    0.000)   34.772
   3.826   (   1.337   41.808    0.000)   41.830
   4.566   (  -2.464   34.110    0.000)   34.198
   5.872   (  13.699    3.380    0.000)   14.110
   6.352   (  40.492   30.186    0.000)   50.505
   9.824   (   7.612    2.712    0.000)    8.080
  11.350   ( -11.689  -16.813    0.000)   20.477
  13.893   (  -3.326   -5.229    0.000)    6.197
  13.937   ( -13.011  -22.846    0.000)   26.292
  14.265   ( -11.673   -3.407    0.000)   12.160
  14.514   (   2.157  -11.370    0.000)   11.573
  14.525   (  -3.489   -2.547    0.000)    4.320
======================= Grid point 35 (18/135) =======================
q-point: ( 0.27  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.098   ( -16.527   40.550    0.000)   43.789
   4.150   (  20.306   21.478    0.000)   29.558
   4.725   ( -18.693   36.060    0.000)   40.618
   6.019   (   3.992   -3.878    0.000)    5.566
   7.322   (  43.094   16.421    0.000)   46.116
  10.168   (  25.155    6.599    0.000)   26.006
  10.947   ( -13.606  -18.995    0.000)   23.365
  13.451   ( -14.291  -26.354    0.000)   29.979
  13.825   (  -0.896   -1.218    0.000)    1.512
  14.039   (  -9.844   -2.781    0.000)   10.229
  14.447   (   3.310  -11.693    0.000)   12.153
  14.472   (  -1.786    1.168    0.000)    2.134
======================= Grid point 36 (19/135) =======================
q-point: ( 0.33  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.119   ( -25.027   38.574    0.000)   45.982
   4.610   (  -4.851   26.950    0.000)   27.383
   4.734   (  -6.546   27.359    0.000)   28.131
   5.970   (  -3.999   -3.392    0.000)    5.244
   8.151   (  37.017    7.968    0.000)   37.864
  10.510   ( -12.945  -20.694    0.000)   24.409
  10.763   (  35.225    3.332    0.000)   35.382
  12.915   ( -15.243  -28.605    0.000)   32.413
  13.820   (  -0.383    0.818    0.000)    0.903
  13.884   (  -4.534   -2.177    0.000)    5.030
  14.394   (   4.013  -11.350    0.000)   12.038
  14.478   (  -0.545    3.225    0.000)    3.271
======================= Grid point 37 (20/135) =======================
q-point: ( 0.40  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.062   ( -24.978   39.092    0.000)   46.391
   4.583   ( -27.057   42.041    0.000)   49.995
   5.002   (   4.570   10.041    0.000)   11.033
   5.865   (  -3.486   -1.690    0.000)    3.874
   8.715   (  20.518   -0.816    0.000)   20.535
  10.137   (  -3.961  -18.449    0.000)   18.869
  11.370   (  30.844   -5.283    0.000)   31.293
  12.386   ( -10.393  -25.434    0.000)   27.476
  13.792   (  -1.291   -0.795    0.000)    1.516
  13.856   (   0.184    0.414    0.000)    0.453
  14.363   (   5.166  -10.638    0.000)   11.826
  14.501   (  -1.036    3.346    0.000)    3.503
======================= Grid point 49 (21/135) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.076   (  15.198   26.323    0.000)   30.396
   4.516   (  14.411   24.961    0.000)   28.823
   5.156   (  12.670   21.945    0.000)   25.340
   5.905   (   0.799    1.383    0.000)    1.597
   6.967   (  18.484   32.015    0.000)   36.968
   9.967   (   6.752   11.696    0.000)   13.505
  10.969   ( -12.064  -20.895    0.000)   24.128
  13.438   ( -15.290  -26.483    0.000)   30.580
  13.826   (  -0.841   -1.457    0.000)    1.683
  14.162   (  -4.492   -7.780    0.000)    8.983
  14.307   (  -4.362   -7.554    0.000)    8.723
  14.495   (  -0.185   -0.320    0.000)    0.369
======================= Grid point 50 (22/135) =======================
q-point: ( 0.27  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.580   (  14.811   21.254    0.000)   25.906
   4.818   (  -9.422   28.918    0.000)   30.414
   5.408   ( -11.828   29.564    0.000)   31.843
   5.944   (   3.354   -3.420    0.000)    4.790
   7.645   (  29.946   17.210    0.000)   34.539
  10.300   (  20.293    7.481    0.000)   21.628
  10.518   ( -13.719  -23.516    0.000)   27.225
  12.883   ( -15.716  -29.713    0.000)   33.614
  13.824   (   0.543    0.830    0.000)    0.991
  13.978   (  -8.313   -3.368    0.000)    8.970
  14.210   (   2.740  -11.068    0.000)   11.402
  14.510   (   0.383    2.477    0.000)    2.506
======================= Grid point 51 (23/135) =======================
q-point: ( 0.33  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.847   ( -20.002   31.838    0.000)   37.599
   4.997   (   9.724   17.280    0.000)   19.828
   5.428   ( -20.236   33.220    0.000)   38.898
   5.907   (  -3.385   -3.057    0.000)    4.561
   8.301   (  28.290    8.443    0.000)   29.523
  10.045   ( -10.883  -25.255    0.000)   27.500
  10.792   (  29.190    1.072    0.000)   29.210
  12.310   ( -13.893  -30.857    0.000)   33.840
  13.830   (  -3.709   -0.593    0.000)    3.756
  13.855   (  -0.204    0.050    0.000)    0.210
  14.162   (   4.504  -11.089    0.000)   11.969
  14.545   (  -0.043    3.294    0.000)    3.294
======================= Grid point 52 (24/135) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.824   ( -19.813   34.317    0.000)   39.626
   5.195   (  -5.201    9.008    0.000)   10.402
   5.410   ( -21.194   36.710    0.000)   42.389
   5.838   (   0.488   -0.845    0.000)    0.976
   8.647   (   3.017   -5.226    0.000)    6.035
   9.767   (   9.861  -17.079    0.000)   19.722
  11.152   (   8.671  -15.018    0.000)   17.342
  11.938   (  10.178  -17.628    0.000)   20.356
  13.785   (  -0.292    0.506    0.000)    0.585
  13.865   (  -0.174    0.301    0.000)    0.348
  14.146   (   5.987  -10.369    0.000)   11.973
  14.563   (  -1.541    2.668    0.000)    3.081
======================= Grid point 66 (25/135) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.996   (  11.388   19.724    0.000)   22.775
   5.221   (   6.688   11.583    0.000)   13.375
   5.847   (   7.422   12.856    0.000)   14.844
   5.878   (  -1.566   -2.713    0.000)    3.133
   8.058   (  13.842   23.974    0.000)   27.683
  10.012   ( -15.185  -26.302    0.000)   30.370
  10.516   (   8.566   14.837    0.000)   17.132
  12.272   ( -17.625  -30.527    0.000)   35.250
  13.848   (   0.731    1.266    0.000)    1.462
  13.898   (  -2.659   -4.605    0.000)    5.318
  14.035   (  -3.333   -5.772    0.000)    6.665
  14.563   (   1.464    2.535    0.000)    2.927
======================= Grid point 67 (26/135) =======================
q-point: ( 0.33  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   5.335   (   4.504   16.010    0.000)   16.631
   5.341   (  -7.243   16.485    0.000)   18.006
   5.787   (  -2.918   -4.346    0.000)    5.234
   5.992   (  -6.920   16.522    0.000)   17.913
   8.545   (  13.802   16.733    0.000)   21.691
   9.510   (  -7.566  -27.433    0.000)   28.457
  10.878   (  15.509    9.053    0.000)   17.958
  11.701   (  -9.903  -28.917    0.000)   30.566
  13.819   (  -2.873    0.246    0.000)    2.884
  13.868   (   0.297    0.602    0.000)    0.672
  13.969   (   2.908   -7.797    0.000)    8.322
  14.601   (   0.034    2.067    0.000)    2.067
======================= Grid point 82 (27/135) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   5.505   (   0.000    0.000    0.000)    0.000
   5.618   (  -5.785  -10.020    0.000)   11.570
   5.618   (   5.785   10.020    0.000)   11.570
   6.169   (   0.000    0.000    0.000)    0.000
   8.973   ( -12.470  -21.599    0.000)   24.941
   8.973   (  12.470   21.599    0.000)   24.941
  11.185   ( -10.644  -18.436    0.000)   21.288
  11.185   (  10.644   18.436    0.000)   21.288
  13.855   (  -1.782   -3.087    0.000)    3.565
  13.855   (   1.782    3.087    0.000)    3.565
  13.874   (  -0.000   -0.000    0.000)    0.000
  14.623   (  -0.000   -0.000    0.000)    0.000
======================= Grid point 241 (28/135) =======================
q-point: ( 0.00  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.813   (  -0.000   -0.000   39.847)   39.847
   0.813   (   0.000   -0.000   39.847)   39.847
   1.870   (  -0.000    0.000   93.640)   93.640
   2.802   (   0.000    0.000    1.828)    1.828
   2.802   (   0.000    0.000    1.828)    1.828
  10.530   (   0.000    0.000  -41.160)   41.160
  12.712   (   0.000    0.000   27.807)   27.807
  14.666   (   0.000   -0.000    0.227)    0.227
  14.666   (   0.000    0.000    0.227)    0.227
  15.028   (  -0.000   -0.000   -5.167)    5.167
  15.068   (  -0.000    0.000   -2.893)    2.893
  15.068   (   0.000    0.000   -2.893)    2.893
======================= Grid point 242 (29/135) =======================
q-point: ( 0.07  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.188   (  26.629   15.374   26.913)   40.863
   1.532   (  32.846   18.963   36.106)   52.365
   2.328   (  44.547   25.719   62.597)   81.020
   2.919   (   9.659    5.577    1.735)   11.288
   3.267   (  36.481   21.062    1.489)   42.151
  10.394   ( -11.373   -6.566  -39.756)   41.869
  12.628   (  -6.942   -4.008   26.912)   28.080
  14.598   (  -5.641   -3.257    0.213)    6.517
  14.653   (  -1.344   -0.776   -0.017)    1.552
  14.971   (  -4.570   -2.638   -5.384)    7.539
  15.008   (  -5.056   -2.919   -2.959)    6.545
  15.012   (  -4.803   -2.773   -2.701)    6.168
======================= Grid point 243 (30/135) =======================
q-point: ( 0.13  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.862   (  29.339   16.939   16.511)   37.687
   2.136   (  24.276   14.016   37.299)   46.658
   3.214   (  14.667    8.468    1.554)   17.008
   3.696   (  65.046   37.555   27.971)   80.148
   4.291   (  47.642   27.506    0.667)   55.017
  10.039   ( -17.842  -10.301  -35.971)   41.454
  12.408   ( -11.549   -6.668   24.840)   28.193
  14.418   (  -9.477   -5.472    0.181)   10.945
  14.575   (  -6.168   -3.561   -1.490)    7.277
  14.841   (  -7.945   -4.587   -3.493)    9.816
  14.846   (  -8.434   -4.869   -3.134)   10.230
  14.859   (  -5.917   -3.416   -3.621)    7.732
======================= Grid point 244 (31/135) =======================
q-point: ( 0.20  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.488   (  23.659   13.659   11.635)   29.693
   2.727   (  26.472   15.284   19.693)   36.361
   3.548   (  12.919    7.459    1.489)   14.992
   5.190   (  61.671   35.606   19.009)   73.705
   5.324   (  38.788   22.394   -0.538)   44.792
   9.664   ( -11.570   -6.680  -30.336)   33.148
  12.111   ( -13.759   -7.944   23.202)   28.120
  14.183   ( -10.135   -5.851    0.150)   11.704
  14.370   ( -10.632   -6.138   -3.806)   12.853
  14.585   ( -13.102   -7.564   -1.161)   15.173
  14.640   (  -8.789   -5.075   -3.375)   10.696
  14.773   (  -2.776   -1.603   -3.740)    4.926
======================= Grid point 245 (32/135) =======================
q-point: ( 0.27  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.929   (  13.864    8.004    9.199)   18.463
   3.348   (  26.535   15.320    4.019)   30.902
   3.769   (   5.447    3.145    1.756)    6.530
   5.916   (  12.505    7.220   -8.340)   16.675
   6.569   (  54.973   31.738   18.864)   66.221
   9.639   (  11.183    6.457  -22.272)   25.744
  11.778   ( -14.705   -8.490   23.450)   28.952
  13.973   (  -7.534   -4.349    0.134)    8.700
  14.115   ( -11.314   -6.532   -5.354)   14.119
  14.247   ( -15.306   -8.837    0.977)   17.700
  14.463   (  -5.990   -3.458   -3.643)    7.818
  14.711   (  -2.663   -1.537   -3.603)    4.737
======================= Grid point 246 (33/135) =======================
q-point: ( 0.33  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.131   (   3.912    2.258    8.025)    9.209
   3.789   (  -3.407   -1.967    2.490)    4.656
   3.934   (  23.362   13.488   -7.472)   27.992
   5.957   (  -5.872   -3.390  -15.673)   17.076
   7.731   (  44.718   25.818   19.603)   55.232
  10.132   (  27.810   16.056  -15.261)   35.554
  11.442   ( -13.465   -7.774   25.887)   30.197
  13.762   ( -19.510  -11.264   -2.815)   22.703
  13.845   (  -3.504   -2.023    0.132)    4.048
  13.980   (  -6.931   -4.002    0.055)    8.004
  14.372   (  -1.917   -1.107   -3.919)    4.501
  14.650   (  -2.582   -1.491   -3.766)    4.803
======================= Grid point 247 (34/135) =======================
q-point: ( 0.40  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.144   (  -1.776   -1.025    7.572)    7.845
   3.647   (  -7.620   -4.399    3.580)    9.499
   4.408   (  16.915    9.766  -15.158)   24.724
   5.771   (  -8.324   -4.806  -20.044)   22.230
   8.632   (  32.499   18.763   21.565)   43.282
  10.659   (   7.741    4.469   -3.929)    9.764
  11.309   (  10.230    5.906   20.491)   23.652
  13.260   ( -23.240  -13.417   -0.692)   26.844
  13.801   (  -0.608   -0.351    0.137)    0.716
  13.878   (  -2.445   -1.411   -0.319)    2.841
  14.362   (   0.654    0.378   -4.166)    4.234
  14.596   (  -1.966   -1.135   -3.885)    4.499
======================= Grid point 248 (35/135) =======================
q-point: ( 0.47  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.095   (  -1.532   -0.884    7.474)    7.680
   3.499   (  -3.920   -2.263    4.475)    6.365
   4.691   (   6.695    3.865  -19.292)   20.783
   5.622   (  -3.778   -2.181  -21.778)   22.211
   9.199   (  14.449    8.342   23.591)   28.895
  10.553   (  -7.114   -4.107   14.333)   16.520
  11.821   (  23.084   13.328    4.424)   27.020
  12.718   ( -20.983  -12.115    0.204)   24.230
  13.799   (   0.134    0.078    0.142)    0.210
  13.848   (  -0.440   -0.254   -0.006)    0.508
  14.381   (   0.623    0.360   -4.318)    4.377
  14.563   (  -0.758   -0.438   -3.915)    4.012
======================= Grid point 257 (36/135) =======================
q-point: ( 0.07  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.758   (  20.709   35.870   17.721)   45.051
   1.971   (  13.520   23.417   39.755)   48.079
   3.154   (  13.537   23.447    4.695)   27.478
   3.278   (  31.975   55.383   30.897)   71.023
   3.959   (  25.435   44.055    0.802)   50.876
  10.147   (  -9.812  -16.995  -37.189)   42.048
  12.477   (  -6.111  -10.584   25.446)   28.229
  14.464   (  -5.483   -9.497   -0.064)   10.967
  14.628   (  -1.351   -2.340   -0.321)    2.721
  14.871   (  -4.064   -7.038   -5.500)    9.813
  14.896   (  -4.853   -8.406   -2.442)   10.009
  14.904   (  -4.135   -7.163   -3.049)    8.815
======================= Grid point 258 (37/135) =======================
q-point: ( 0.13  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.448   (  15.448   31.972   16.188)   39.025
   2.518   (  20.409   28.868   20.842)   41.040
   3.572   (   8.528   25.682    1.223)   27.088
   4.580   (  50.433   47.691   23.237)   73.198
   4.898   (  35.438   30.408   -0.604)   46.700
   9.775   ( -12.151  -14.399  -32.709)   37.747
  12.222   (  -9.979  -11.531   23.720)   28.200
  14.245   (  -7.318  -10.534   -0.426)   12.834
  14.520   (  -9.157   -2.384   -1.857)    9.643
  14.658   (  -8.730  -12.619   -2.037)   15.479
  14.736   (  -6.883   -6.407   -3.133)    9.911
  14.751   (  -1.272   -7.676   -4.463)    8.970
======================= Grid point 259 (38/135) =======================
q-point: ( 0.20  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.996   (  11.465   29.104    6.325)   31.913
   3.134   (  15.894   29.350   11.815)   35.407
   3.920   (  -0.906   28.118    1.239)   28.159
   5.663   (  25.360   11.676   -5.027)   28.368
   5.922   (  54.230   34.456   19.950)   67.276
   9.556   (   0.633   -2.929  -26.168)   26.339
  11.910   ( -12.068  -11.897   23.284)   28.798
  14.025   (  -6.263   -8.385   -0.558)   10.481
  14.279   ( -13.407   -5.807   -4.089)   15.172
  14.345   ( -10.869  -13.479   -0.146)   17.316
  14.556   (  -7.221   -3.831   -3.353)    8.835
  14.676   (   0.902   -7.878   -3.715)    8.756
======================= Grid point 260 (39/135) =======================
q-point: ( 0.27  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.358   (  -5.492   31.406    4.999)   32.272
   3.696   (  20.330   19.602   -0.762)   28.252
   4.080   ( -12.308   26.315    1.429)   29.086
   5.951   (   4.037   -2.728  -12.417)   13.339
   7.138   (  50.593   24.354   19.594)   59.470
   9.800   (  22.770    8.762  -18.192)   30.433
  11.569   ( -12.700  -12.334   24.841)   30.504
  13.870   (  -4.894   -6.270   -1.306)    8.060
  13.911   ( -13.207  -16.135   -1.156)   20.883
  14.079   (  -9.671   -4.157   -1.118)   10.586
  14.427   (  -4.643   -0.621   -3.719)    5.981
  14.614   (   1.025   -7.937   -3.685)    8.811
======================= Grid point 261 (40/135) =======================
q-point: ( 0.33  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.461   ( -15.824   28.233    5.878)   32.894
   4.026   ( -20.082   24.156    2.393)   31.504
   4.220   (  17.217   14.953  -11.415)   25.501
   5.858   (  -6.741   -5.548  -18.359)   20.329
   8.152   (  41.973   16.109   20.886)   49.573
  10.386   (  28.256    7.959  -12.682)   31.978
  11.260   (  -6.299  -10.191   27.180)   29.704
  13.429   ( -16.942  -21.576   -0.798)   27.444
  13.815   (  -1.191   -0.622    0.046)    1.345
  13.917   (  -4.988   -2.497   -0.989)    5.665
  14.382   (  -1.474    1.782   -4.082)    4.692
  14.555   (   1.346   -7.679   -3.777)    8.663
======================= Grid point 262 (41/135) =======================
q-point: ( 0.40  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.424   ( -18.685   26.767    6.296)   33.245
   3.887   ( -23.094   26.680    3.279)   35.439
   4.600   (  11.017    9.372  -17.590)   22.773
   5.684   (  -6.475   -3.671  -21.393)   22.651
   8.901   (  28.772    8.249   22.885)   37.678
  10.584   (  -3.687  -10.041    9.060)   14.017
  11.436   (  28.444    4.231    9.510)   30.289
  12.892   ( -18.450  -22.602    0.084)   29.176
  13.801   (  -0.650    0.051    0.146)    0.669
  13.853   (  -0.712   -0.865   -0.411)    1.193
  14.394   (  -0.280    2.480   -4.359)    5.023
  14.510   (   2.213   -7.051   -3.786)    8.304
======================= Grid point 263 (42/135) =======================
q-point: ( 0.47  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.390   ( -16.214   28.084    6.439)   33.062
   3.812   ( -17.953   31.095    3.686)   36.094
   4.752   (  -1.323    2.291  -19.792)   19.968
   5.601   (  -0.070    0.122  -22.191)   22.192
   9.226   (   3.305   -5.725   24.155)   25.043
  10.442   (   3.708   -6.423   15.822)   17.474
  11.864   (   6.712  -11.626    3.634)   13.908
  12.488   (   4.483   -7.765    0.404)    8.976
  13.794   (   0.313   -0.542    0.167)    0.648
  13.848   (  -0.143    0.248   -0.136)    0.317
  14.406   (  -1.099    1.903   -4.478)    4.988
  14.494   (   3.561   -6.167   -3.756)    8.051
======================= Grid point 274 (43/135) =======================
q-point: ( 0.13  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.959   (  14.721   25.498    9.957)   31.081
   3.293   (  21.653   37.503    8.527)   44.137
   4.136   (  15.784   27.339    0.885)   31.581
   5.444   (  13.098   22.686   -3.912)   26.486
   5.567   (  28.881   50.023   21.776)   61.731
   9.539   (  -3.925   -6.798  -27.759)   28.848
  11.956   (  -8.490  -14.705   23.446)   28.949
  14.036   (  -5.462   -9.461   -0.879)   10.960
  14.343   ( -10.558  -18.287   -1.449)   21.165
  14.446   (  -3.833   -6.639   -1.575)    7.826
  14.576   (  -4.554   -7.888   -4.806)   10.298
  14.602   (  -3.766   -6.522   -3.308)    8.225
======================= Grid point 275 (44/135) =======================
q-point: ( 0.20  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.480   (  16.652   22.638   -3.195)   28.284
   3.907   (   3.224   39.729    7.473)   40.554
   4.574   (  -2.429   33.860    0.769)   33.956
   5.777   (  11.265    0.616   -9.773)   14.926
   6.580   (  39.833   31.393   21.286)   55.003
   9.574   (   9.934    4.979  -21.099)   23.846
  11.635   ( -10.819  -15.394   24.713)   31.060
  13.882   (  -3.833   -6.206   -1.240)    7.399
  13.934   ( -12.250  -21.357   -0.292)   24.622
  14.240   ( -11.860   -3.320   -2.413)   12.550
  14.473   (   2.011  -10.651   -4.025)   11.562
  14.489   (  -3.050   -3.527   -3.498)    5.828
======================= Grid point 276 (45/135) =======================
q-point: ( 0.27  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.953   (   8.766   23.037  -11.497)   27.198
   4.226   (  -6.485   34.332    6.292)   35.501
   4.729   ( -18.662   35.419    0.426)   40.037
   5.855   (   1.509   -5.792  -16.041)   17.121
   7.556   (  42.790   17.946   21.625)   51.192
   9.962   (  25.361    6.972  -14.876)   30.218
  11.281   ( -10.728  -16.366   27.155)   33.471
  13.451   ( -14.195  -25.825   -0.211)   29.470
  13.821   (  -0.713   -1.245   -0.341)    1.475
  14.018   (  -9.464   -2.707   -1.997)   10.044
  14.411   (   3.340  -11.434   -3.636)   12.455
  14.431   (  -1.894    0.808   -4.053)    4.546
======================= Grid point 277 (46/135) =======================
q-point: ( 0.33  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.131   ( -19.558   34.618   -0.767)   39.768
   4.514   (   7.666   18.157   -9.376)   21.825
   4.693   ( -24.671   36.066   -0.611)   43.702
   5.755   (  -5.595   -4.464  -20.335)   21.558
   8.398   (  37.248    9.269   22.840)   44.665
  10.419   (  15.140   -5.061   -4.821)   16.676
  11.048   (   8.718  -10.185   23.092)   26.701
  12.921   ( -15.210  -28.235    0.358)   32.073
  13.819   (  -0.196    0.844    0.033)    0.867
  13.871   (  -4.139   -2.154   -1.241)    4.828
  14.359   (   3.902  -11.085   -3.522)   12.268
  14.433   (  -0.547    2.977   -4.432)    5.366
======================= Grid point 278 (47/135) =======================
q-point: ( 0.40  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.100   ( -23.397   37.274    3.072)   44.116
   4.601   ( -26.026   40.803    1.445)   48.418
   4.792   (   2.813   10.142  -18.032)   20.879
   5.625   (  -3.816   -2.268  -22.371)   22.807
   8.979   (  20.765    0.031   24.315)   31.975
  10.312   (  -1.979  -16.045   15.544)   22.427
  11.420   (  28.777   -6.344    5.190)   29.922
  12.395   ( -10.470  -25.248    0.662)   27.341
  13.790   (  -1.129   -0.988   -0.142)    1.507
  13.852   (   0.310    0.532   -0.289)    0.680
  14.328   (   4.999  -10.319   -3.486)   11.984
  14.454   (  -0.953    3.115   -4.601)    5.638
======================= Grid point 290 (48/135) =======================
q-point: ( 0.20  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.934   (  12.881   22.311  -11.915)   28.384
   4.572   (  14.164   24.533    5.420)   28.842
   5.162   (  12.609   21.839    0.652)   25.226
   5.749   (  -0.980   -1.697  -15.426)   15.550
   7.219   (  19.037   32.974   23.069)   44.518
   9.757   (   7.621   13.201  -16.025)   22.117
  11.293   ( -10.660  -18.464   27.045)   34.438
  13.440   ( -15.067  -26.097    0.016)   30.135
  13.820   (  -0.786   -1.361   -0.527)    1.657
  14.141   (  -4.459   -7.723   -2.127)    9.168
  14.268   (  -4.162   -7.208   -3.781)    9.142
  14.453   (  -0.330   -0.572   -4.105)    4.158
======================= Grid point 291 (49/135) =======================
q-point: ( 0.27  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.364   (  11.830   19.369  -18.473)   29.263
   4.882   (  -7.081   27.279    5.470)   28.709
   5.395   ( -13.006   28.285   -2.052)   31.199
   5.746   (   2.908   -4.549  -17.871)   18.669
   7.910   (  29.625   18.663   24.121)   42.518
  10.085   (  16.869    5.447  -10.671)   20.691
  10.921   (  -8.255  -18.547   28.369)   34.884
  12.891   ( -15.610  -29.409    0.632)   33.301
  13.821   (   0.660    0.893   -0.275)    1.144
  13.960   (  -8.052   -3.195   -1.798)    8.848
  14.179   (   2.823  -10.738   -3.121)   11.534
  14.463   (   0.291    2.241   -4.572)    5.100
======================= Grid point 292 (50/135) =======================
q-point: ( 0.33  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.694   (   0.923   20.310  -19.719)   28.323
   4.963   ( -11.110   25.978    3.220)   28.437
   5.393   ( -17.567   27.957   -6.031)   33.564
   5.694   (  -6.421    0.754  -16.970)   18.160
   8.574   (  28.087    9.605   24.872)   38.726
  10.187   (  -2.967  -17.422   11.418)   21.041
  10.887   (  21.933   -4.685   10.512)   24.769
  12.322   ( -13.914  -30.643    0.957)   33.667
  13.824   (  -4.107   -1.221   -0.740)    4.348
  13.849   (   0.605    0.734   -0.381)    1.025
  14.134   (   4.380  -10.692   -2.850)   11.900
  14.496   (  -0.045    3.089   -4.877)    5.773
======================= Grid point 293 (51/135) =======================
q-point: ( 0.40  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.782   ( -14.577   25.248  -11.187)   31.227
   5.032   (  -9.587   16.605   -6.774)   20.336
   5.378   ( -16.431   28.460   -7.363)   33.677
   5.616   (  -3.433    5.946  -16.256)   17.647
   8.927   (   2.646   -4.583   25.439)   25.984
   9.983   (   8.797  -15.237   19.438)   26.218
  11.188   (   8.931  -15.469    3.683)   18.238
  11.950   (  10.076  -17.452    1.013)   20.177
  13.779   (  -0.163    0.282   -0.588)    0.672
  13.862   (  -0.210    0.365   -0.107)    0.434
  14.118   (   5.767   -9.989   -2.774)   11.864
  14.513   (  -1.448    2.508   -4.982)    5.763
======================= Grid point 307 (52/135) =======================
q-point: ( 0.27  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.742   (  10.314   17.864  -22.589)   30.590
   5.268   (   6.598   11.428    4.097)   13.818
   5.624   (  -2.490   -4.313  -23.327)   23.852
   5.850   (   7.459   12.919    0.422)   14.924
   8.344   (  14.141   24.494   25.689)   38.208
  10.105   (  -5.196   -9.000    6.362)   12.185
  10.665   (  -0.080   -0.139   16.000)   16.001
  12.285   ( -17.530  -30.362    1.094)   35.076
  13.845   (   0.764    1.324   -0.134)    1.534
  13.883   (  -2.554   -4.423   -1.465)    5.313
  14.011   (  -3.143   -5.443   -2.405)    6.729
  14.512   (   1.368    2.369   -4.980)    5.682
======================= Grid point 308 (53/135) =======================
q-point: ( 0.33  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   5.045   (   2.852   15.699  -25.115)   29.755
   5.397   (  -4.637   14.464    2.889)   15.462
   5.511   (  -3.783   -5.212  -23.657)   24.518
   5.997   (  -6.851   17.271    0.532)   18.588
   8.835   (  13.621   17.030   26.033)   33.959
   9.755   (  -6.823  -24.360   22.041)   33.553
  10.908   (  14.874    7.635    3.057)   16.997
  11.716   (  -9.932  -28.792    1.291)   30.484
  13.808   (  -2.671    0.143   -1.089)    2.888
  13.866   (   0.305    0.654   -0.030)    0.722
  13.949   (   2.782   -7.414   -2.025)    8.173
  14.548   (   0.035    1.945   -5.208)    5.560
======================= Grid point 323 (54/135) =======================
q-point: ( 0.33  0.33  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   5.323   (  -5.772   -9.997  -26.568)   28.967
   5.323   (   5.772    9.997  -26.568)   28.967
   5.558   (  -0.000   -0.000    5.265)    5.265
   6.175   (   0.000    0.000    0.477)    0.477
   9.256   ( -12.005  -20.793   25.484)   35.013
   9.256   (  12.005   20.793   25.484)   35.013
  11.203   ( -10.547  -18.267    1.563)   21.151
  11.203   (  10.547   18.267    1.563)   21.151
  13.841   (  -1.669   -2.890   -1.405)    3.621
  13.841   (   1.669    2.890   -1.405)    3.621
  13.873   (  -0.000   -0.000    0.018)    0.018
  14.569   (  -0.000   -0.000   -5.333)    5.333
======================= Grid point 482 (55/135) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.552   (   0.000    0.000   34.299)   34.299
   1.552   (   0.000    0.000   34.299)   34.299
   2.829   (   0.000    0.000    0.079)    0.079
   2.829   (  -0.000    0.000    0.079)    0.079
   3.685   (   0.000    0.000   89.559)   89.559
   9.560   (  -0.000    0.000  -55.528)   55.528
  13.284   (  -0.000    0.000   28.555)   28.555
  14.677   (  -0.000   -0.000    1.090)    1.090
  14.677   (   0.000    0.000    1.090)    1.090
  14.875   (   0.000   -0.000  -10.337)   10.337
  14.989   (   0.000    0.000   -4.892)    4.892
  14.989   (  -0.000   -0.000   -4.892)    4.892
======================= Grid point 483 (56/135) =======================
q-point: ( 0.07  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.770   (  17.323   10.001   29.687)   35.797
   2.138   (  39.074   22.560   28.154)   53.182
   2.945   (   9.554    5.516    0.068)   11.032
   3.284   (  35.506   20.500   -0.136)   40.999
   3.831   (  14.706    8.490   81.132)   82.890
   9.450   (  -9.143   -5.279  -54.199)   55.218
  13.190   (  -7.767   -4.484   28.363)   29.747
  14.610   (  -5.661   -3.268    1.084)    6.626
  14.656   (  -2.425   -1.400    0.286)    2.815
  14.814   (  -4.532   -2.616  -10.346)   11.594
  14.927   (  -5.203   -3.004   -4.999)    7.816
  14.938   (  -4.327   -2.498   -4.534)    6.747
======================= Grid point 484 (57/135) =======================
q-point: ( 0.13  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.264   (  23.345   13.478   22.376)   35.034
   2.897   (  22.142   12.784   36.008)   44.162
   3.237   (  14.598    8.428    0.090)   16.856
   4.269   (  45.677   26.372    0.204)   52.744
   4.519   (  44.933   25.942   50.042)   72.084
   9.172   ( -13.306   -7.682  -49.880)   52.193
  12.946   ( -12.652   -7.304   27.857)   31.455
  14.429   (  -9.500   -5.485    1.076)   11.022
  14.514   (  -9.872   -5.699   -5.734)   12.760
  14.725   (  -2.847   -1.644   -6.215)    7.031
  14.761   (  -8.685   -5.015   -5.288)   11.338
  14.794   (  -7.798   -4.502   -3.497)    9.660
======================= Grid point 485 (58/135) =======================
q-point: ( 0.20  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.785   (  20.318   11.731   17.157)   29.065
   3.258   (  12.100    6.986   32.396)   35.280
   3.571   (  13.090    7.558    0.305)   15.119
   5.229   (  32.686   18.871   -7.048)   38.395
   5.744   (  56.789   32.787   32.825)   73.331
   8.940   (  -3.389   -1.956  -40.747)   40.934
  12.630   ( -14.059   -8.117   27.552)   31.979
  14.194   ( -10.142   -5.855    1.086)   11.761
  14.246   ( -12.359   -7.135   -9.067)   16.907
  14.548   (  -9.092   -5.249   -5.688)   11.941
  14.557   ( -11.236   -6.487   -1.508)   13.062
  14.680   (  -2.631   -1.519   -5.110)    5.945
======================= Grid point 486 (59/135) =======================
q-point: ( 0.27  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.167   (  12.067    6.967   14.002)   19.753
   3.534   (  12.210    7.049   15.724)   21.119
   3.801   (   6.061    3.499    0.953)    7.063
   5.653   (   3.670    2.119  -18.580)   19.057
   7.047   (  54.114   31.243   27.732)   68.363
   9.147   (  21.743   12.553  -25.720)   35.943
  12.308   ( -13.603   -7.854   28.207)   32.285
  13.962   ( -12.114   -6.994  -10.252)   17.343
  13.983   (  -7.519   -4.341    1.129)    8.755
  14.270   ( -12.878   -7.435    1.464)   14.942
  14.364   (  -6.314   -3.645   -6.139)    9.531
  14.614   (  -3.041   -1.756   -5.901)    6.867
======================= Grid point 487 (60/135) =======================
q-point: ( 0.33  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.339   (   3.056    1.765   12.259)   12.757
   3.825   (  12.521    7.229   -1.253)   14.512
   3.842   (  -2.299   -1.327    2.218)    3.459
   5.523   ( -11.199   -6.466  -27.156)   30.078
   8.216   (  45.932   26.519   27.432)   59.712
   9.841   (  34.366   19.841  -13.981)   42.073
  12.005   ( -12.189   -7.037   29.413)   32.607
  13.650   ( -15.843   -9.147   -8.802)   20.301
  13.856   (  -3.476   -2.007    1.199)    4.189
  14.009   (  -8.420   -4.861    2.990)   10.172
  14.265   (  -2.237   -1.292   -6.607)    7.094
  14.543   (  -2.931   -1.692   -6.780)    7.578
======================= Grid point 488 (61/135) =======================
q-point: ( 0.40  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.340   (  -2.077   -1.199   11.433)   11.682
   3.728   (  -6.370   -3.678    3.922)    8.336
   4.094   (  10.138    5.853  -13.663)   17.992
   5.254   ( -10.495   -6.059  -30.356)   32.685
   9.157   (  34.552   19.949   29.238)   49.464
  10.581   (  25.453   14.695   -4.655)   29.756
  11.788   (  -3.975   -2.295   26.423)   26.818
  13.214   ( -21.224  -12.254   -4.678)   24.950
  13.813   (  -0.584   -0.337    1.272)    1.439
  13.889   (  -2.528   -1.459    1.973)    3.523
  14.248   (   0.403    0.233   -7.030)    7.046
  14.484   (  -2.113   -1.220   -7.310)    7.707
======================= Grid point 489 (62/135) =======================
q-point: ( 0.47  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.288   (  -1.586   -0.916   11.145)   11.294
   3.603   (  -3.334   -1.925    5.262)    6.520
   4.268   (   4.224    2.438  -20.253)   20.832
   5.077   (  -4.197   -2.423  -31.009)   31.385
   9.769   (  15.946    9.206   31.532)   36.514
  10.855   (   0.290    0.167   13.337)   13.341
  11.956   (  15.651    9.036    9.459)   20.398
  12.708   ( -19.858  -11.465   -1.685)   22.992
  13.811   (   0.144    0.083    1.318)    1.329
  13.861   (  -0.330   -0.191    1.724)    1.766
  14.263   (   0.526    0.304   -7.292)    7.317
  14.449   (  -0.780   -0.451   -7.471)    7.526
======================= Grid point 498 (63/135) =======================
q-point: ( 0.07  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.185   (  16.259   28.161   23.486)   40.112
   2.709   (  15.351   26.588   32.713)   44.864
   3.189   (  10.304   17.847    0.201)   20.608
   3.961   (  24.668   42.726   -1.290)   49.353
   4.245   (  20.229   35.038   61.434)   73.560
   9.254   (  -7.608  -13.177  -51.409)   53.614
  13.023   (  -6.733  -11.662   28.022)   31.089
  14.463   (  -6.009  -10.408   -0.375)   12.024
  14.624   (  -1.777   -3.079    0.049)    3.555
  14.712   (  -3.839   -6.650  -10.049)   12.647
  14.821   (  -4.217   -7.305   -5.116)    9.865
  14.832   (  -4.541   -7.866   -3.876)    9.875
======================= Grid point 499 (64/135) =======================
q-point: ( 0.13  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.775   (  12.323   34.394   17.898)   40.684
   3.116   (  13.163    9.438   34.777)   38.364
   3.591   (   8.641   25.370    0.116)   26.801
   4.849   (  32.320   27.860   -4.282)   42.884
   5.213   (  42.645   41.293   38.641)   70.830
   8.987   (  -7.220   -9.087  -44.969)   46.443
  12.748   ( -10.507  -11.910   27.689)   31.920
  14.227   (  -7.285  -11.806   -1.916)   14.005
  14.455   ( -13.090   -2.401   -5.015)   14.222
  14.595   (  -4.808  -10.711   -4.095)   12.434
  14.649   (  -3.602   -6.712   -5.014)    9.119
  14.657   (  -4.843   -8.419   -4.213)   10.587
======================= Grid point 500 (65/135) =======================
q-point: ( 0.20  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.221   (   4.805   23.982   15.362)   28.883
   3.454   (  10.197   18.007   20.813)   29.350
   3.940   (  -0.436   27.752    0.299)   27.757
   5.490   (  18.167    6.002  -13.102)   23.189
   6.432   (  51.146   34.258   29.866)   68.422
   8.953   (   9.680    4.107  -32.696)   34.345
  12.436   ( -11.748  -10.980   28.070)   32.349
  13.998   (  -6.250   -9.614   -3.346)   11.945
  14.174   ( -14.398   -4.670   -5.523)   16.113
  14.335   (  -9.021  -12.534   -0.908)   15.470
  14.468   (  -7.529   -4.309   -5.046)   10.036
  14.581   (   0.892   -8.065   -5.433)    9.765
======================= Grid point 501 (66/135) =======================
q-point: ( 0.27  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.477   (   0.956   17.949    7.716)   19.561
   3.775   (   3.720   20.747    9.321)   23.047
   4.108   ( -11.296   25.571    0.864)   27.968
   5.590   (  -3.032   -7.050  -23.651)   24.865
   7.625   (  50.156   25.569   27.656)   62.723
   9.424   (  31.008   13.983  -18.776)   38.854
  12.123   ( -11.612  -10.605   29.321)   33.272
  13.780   ( -11.061  -11.108   -8.375)   17.773
  13.905   (  -6.576   -5.889    0.811)    8.865
  14.073   (  -7.873   -8.410    0.775)   11.546
  14.327   (  -5.161   -0.934   -5.958)    7.938
  14.512   (   0.650   -8.059   -6.283)   10.239
======================= Grid point 502 (67/135) =======================
q-point: ( 0.33  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.596   (  -9.902   20.714    6.566)   23.879
   4.023   (   5.472   13.702   -4.674)   15.477
   4.079   ( -19.207   23.519    1.748)   30.415
   5.374   ( -10.390   -7.446  -28.941)   31.639
   8.662   (  43.163   17.726   28.430)   54.640
  10.172   (  35.471   12.437   -9.735)   38.828
  11.832   (  -8.510  -10.461   29.608)   32.534
  13.377   ( -16.108  -18.657   -5.197)   25.190
  13.824   (  -0.960   -0.836    1.164)    1.724
  13.910   (  -3.749   -3.998    0.811)    5.541
  14.272   (  -1.834    1.429   -6.770)    7.158
  14.447   (   1.034   -7.562   -6.935)   10.313
======================= Grid point 503 (68/135) =======================
q-point: ( 0.40  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.579   ( -15.712   22.847    8.338)   28.955
   3.959   ( -21.042   25.034    3.241)   32.863
   4.233   (   4.837    8.368  -16.132)   18.806
   5.147   (  -7.699   -4.327  -30.620)   31.868
   9.454   (  30.490    9.667   30.492)   44.191
  10.733   (  15.475    0.431    4.062)   16.005
  11.728   (  10.280   -4.234   19.421)   22.378
  12.874   ( -17.705  -21.207   -2.458)   27.735
  13.814   (  -0.553    0.040    1.340)    1.451
  13.855   (   0.155   -1.382    1.015)    1.721
  14.275   (  -0.390    2.060   -7.321)    7.616
  14.401   (   1.988   -6.673   -7.175)    9.998
======================= Grid point 504 (69/135) =======================
q-point: ( 0.47  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.551   ( -14.337   24.832    8.887)   30.020
   3.894   ( -16.614   28.776    3.884)   33.454
   4.324   (  -2.070    3.585  -20.018)   20.441
   5.052   (  -0.039    0.067  -30.816)   30.816
   9.808   (   3.212   -5.563   32.030)   32.667
  10.788   (   3.448   -5.973   16.295)   17.695
  11.967   (   6.233  -10.796    6.790)   14.196
  12.488   (   4.336   -7.510   -0.765)    8.705
  13.807   (   0.369   -0.640    1.452)    1.629
  13.855   (   0.098   -0.170    1.165)    1.181
  14.284   (  -0.876    1.517   -7.541)    7.742
  14.384   (   3.273   -5.668   -7.183)    9.718
======================= Grid point 515 (70/135) =======================
q-point: ( 0.13  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.290   (   5.772    9.998   23.479)   26.164
   3.516   (  19.893   34.455   13.205)   41.920
   4.151   (  15.768   27.311    0.165)   31.536
   5.303   (   9.281   16.075  -10.968)   21.560
   6.118   (  27.513   47.654   32.167)   63.739
   8.891   (   1.049    1.818  -35.673)   35.735
  12.485   (  -8.040  -13.925   28.219)   32.479
  14.003   (  -5.715   -9.899   -3.064)   11.834
  14.292   (  -9.446  -16.361   -3.867)   19.284
  14.403   (  -4.413   -7.644   -2.755)    9.247
  14.461   (  -4.621   -8.003   -5.709)   10.863
  14.515   (  -3.876   -6.714   -5.120)    9.291
======================= Grid point 516 (71/135) =======================
q-point: ( 0.20  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.504   (   7.357   10.440    7.328)   14.724
   4.096   (   4.240   36.229   11.058)   38.115
   4.587   (  -2.713   33.578    0.079)   33.687
   5.480   (   5.773   -5.756  -20.163)   21.749
   7.102   (  38.316   32.599   29.407)   58.272
   9.120   (  17.076   12.376  -23.280)   31.412
  12.191   (  -9.895  -13.442   29.634)   34.011
  13.823   (  -6.888  -10.392   -5.638)   13.683
  13.934   (  -8.900  -15.064    0.678)   17.510
  14.179   ( -11.236   -3.973   -3.322)   12.372
  14.367   (   1.350   -9.627   -6.169)   11.513
  14.397   (  -2.420   -5.239   -5.379)    7.890
======================= Grid point 517 (72/135) =======================
q-point: ( 0.27  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.750   (   7.762   11.469   -6.385)   15.250
   4.386   ( -10.937   34.591    9.307)   37.454
   4.725   ( -18.099   32.411   -1.714)   37.161
   5.430   (  -3.044   -7.118  -25.096)   26.263
   8.077   (  42.332   20.072   28.868)   55.030
   9.684   (  31.244   11.727  -13.121)   35.859
  11.871   ( -10.345  -14.754   30.855)   35.732
  13.423   ( -14.124  -23.455   -3.458)   27.597
  13.819   (  -0.239   -1.709    0.383)    1.767
  13.973   (  -7.882   -3.152   -2.056)    8.734
  14.308   (   2.366   -8.619   -6.526)   11.067
  14.323   (  -1.539   -2.283   -6.333)    6.906
======================= Grid point 518 (73/135) =======================
q-point: ( 0.33  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.975   (   1.777   14.639  -12.516)   19.342
   4.466   ( -15.076   29.821    4.894)   33.772
   4.675   ( -22.730   31.234   -3.486)   38.786
   5.263   (  -9.912   -1.686  -25.888)   27.772
   8.944   (  38.000   11.147   29.873)   49.605
  10.335   (  30.170    4.146   -4.496)   30.783
  11.572   (  -3.744  -15.478   28.329)   32.498
  12.912   ( -15.140  -26.803   -1.892)   30.841
  13.827   (   0.815    0.961    0.869)    1.531
  13.847   (  -3.280   -2.597   -0.727)    4.247
  14.259   (   3.389  -10.295   -6.404)   12.589
  14.312   (  -0.560    2.280   -7.426)    7.788
======================= Grid point 519 (74/135) =======================
q-point: ( 0.40  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.081   (  -9.617   20.700   -8.691)   24.424
   4.507   ( -10.739   22.525   -4.242)   25.312
   4.605   ( -21.993   33.924   -3.569)   40.587
   5.101   (  -7.081    1.934  -26.151)   27.162
   9.558   (  21.705    1.145   31.530)   38.296
  10.627   (   6.384   -9.749   13.125)   17.552
  11.588   (  19.596  -10.360   12.146)   25.276
  12.399   ( -10.364  -24.490   -0.689)   26.602
  13.792   (  -0.408   -1.736    0.583)    1.876
  13.854   (   0.510    0.838    0.877)    1.316
  14.228   (   4.415   -9.325   -6.549)   12.220
  14.328   (  -0.701    2.441   -7.735)    8.141
======================= Grid point 531 (75/135) =======================
q-point: ( 0.20  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.729   (   7.069   12.245   -5.955)   15.342
   4.711   (  13.184   22.836    7.766)   27.489
   5.177   (  12.552   21.741    0.464)   25.109
   5.326   (  -4.012   -6.948  -25.911)   27.125
   7.766   (  19.665   34.061   29.891)   49.400
   9.446   (  11.572   20.043  -14.883)   27.517
  11.888   (  -9.708  -16.814   31.310)   36.841
  13.423   ( -14.291  -24.753   -2.391)   28.682
  13.811   (  -0.752   -1.303   -0.160)    1.513
  14.081   (  -4.352   -7.539   -3.873)    9.527
  14.171   (  -3.712   -6.430   -5.360)    9.158
  14.343   (  -0.771   -1.336   -6.697)    6.873
======================= Grid point 532 (76/135) =======================
q-point: ( 0.27  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.991   (   8.019   12.700  -16.178)   22.075
   4.981   (  -5.669   18.260   -3.010)   19.355
   5.169   ( -10.116    8.497  -16.134)   20.853
   5.497   (  -2.826   19.152   -3.213)   19.624
   8.474   (  29.376   20.520   30.355)   46.963
   9.914   (  22.520   12.026   -7.541)   26.620
  11.523   ( -10.272  -19.662   31.206)   38.287
  12.893   ( -15.423  -28.459   -1.017)   32.386
  13.817   (   0.973    1.080    0.022)    1.454
  13.914   (  -6.757   -3.025   -2.546)    7.829
  14.095   (   2.595   -9.665   -5.009)   11.191
  14.339   (   0.015    1.568   -7.631)    7.791
======================= Grid point 533 (77/135) =======================
q-point: ( 0.33  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.237   (   4.589   12.603  -22.346)   26.062
   4.959   (  -5.779    3.289  -25.917)   26.756
   5.116   ( -16.481   24.591    3.515)   29.811
   5.532   ( -17.375   26.750   -1.565)   31.936
   9.154   (  28.310   11.005   31.069)   43.450
  10.361   (  16.166   -1.363    4.392)   16.807
  11.228   (   3.515  -17.307   22.804)   28.843
  12.333   ( -13.924  -29.961   -0.340)   33.040
  13.801   (  -2.983   -1.828   -1.299)    3.732
  13.853   (   0.922    1.320    0.959)    1.875
  14.054   (   3.797   -9.497   -5.028)   11.397
  14.363   (  -0.060    2.538   -8.199)    8.584
======================= Grid point 534 (78/135) =======================
q-point: ( 0.40  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.353   (  -5.323    9.220  -24.264)   26.497
   4.916   (  -0.183    0.318  -31.543)   31.545
   5.091   ( -17.195   29.783    7.647)   35.230
   5.504   ( -18.405   31.878   -0.449)   36.812
   9.520   (   2.428   -4.206   31.805)   32.173
  10.399   (   6.544  -11.335   18.960)   23.039
  11.305   (   9.740  -16.870    8.626)   21.304
  11.965   (   9.836  -17.037    0.094)   19.673
  13.766   (   0.288   -0.499   -0.540)    0.790
  13.870   (  -0.312    0.541    1.166)    1.323
  14.039   (   5.102   -8.837   -5.114)   11.414
  14.377   (  -1.204    2.085   -8.383)    8.722
======================= Grid point 548 (79/135) =======================
q-point: ( 0.27  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.257   (   7.685   13.310  -23.113)   27.756
   5.004   (  -3.105   -5.378  -32.193)   32.787
   5.430   (   6.131   10.620    7.394)   14.320
   5.858   (   7.364   12.755   -0.011)   14.728
   8.934   (  14.538   25.180   31.185)   42.637
  10.179   (   7.395   12.809    0.711)   14.808
  11.117   ( -10.534  -18.245   27.157)   34.371
  12.301   ( -17.277  -29.925   -0.090)   34.555
  13.843   (  -2.240   -3.879   -2.423)    5.092
  13.847   (   0.857    1.485    0.207)    1.727
  13.949   (  -2.625   -4.546   -3.309)    6.205
  14.377   (   1.107    1.918   -8.365)    8.654
======================= Grid point 549 (80/135) =======================
q-point: ( 0.33  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.491   (   2.470   12.401  -27.534)   30.299
   4.888   (  -2.024   -6.472  -32.808)   33.502
   5.548   (  -6.202   15.147    8.500)   18.443
   6.004   (  -6.582   17.120   -0.188)   18.343
   9.431   (  13.696   17.137   31.389)   38.295
  10.229   (  -2.734  -14.023   21.765)   26.035
  11.009   (  10.666   -0.231    8.023)   13.348
  11.738   (  -9.912  -28.424    0.358)   30.105
  13.779   (  -2.033   -0.172   -1.679)    2.642
  13.872   (   1.119   -0.014    0.236)    1.144
  13.895   (   1.577   -5.513   -2.543)    6.273
  14.406   (   0.037    1.622   -8.752)    8.901
======================= Grid point 564 (81/135) =======================
q-point: ( 0.33  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 192
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.721   (  -5.030   -8.712  -31.038)   32.627
   4.721   (   5.030    8.712  -31.038)   32.627
   5.701   (  -0.000   -0.000    8.683)    8.683
   6.179   (   0.000    0.000   -0.479)    0.479
   9.835   ( -10.982  -19.021   30.136)   37.290
   9.835   (  10.982   19.021   30.136)   37.290
  11.234   ( -10.063  -17.430    1.203)   20.163
  11.234   (  10.063   17.430    1.203)   20.163
  13.802   (  -1.322   -2.290   -2.358)    3.543
  13.802   (   1.322    2.290   -2.358)    3.543
  13.883   (  -0.000   -0.000    1.345)    1.345
  14.423   (  -0.000   -0.000   -8.955)    8.955
======================= Grid point 723 (82/135) =======================
q-point: ( 0.00  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.150   (  -0.000    0.000   25.783)   25.783
   2.150   (   0.000    0.000   25.783)   25.783
   2.775   (   0.000    0.000   -6.256)    6.256
   2.775   (   0.000   -0.000   -6.256)    6.256
   5.396   (   0.000    0.000   83.186)   83.186
   8.358   (   0.000   -0.000  -65.825)   65.825
  13.816   (  -0.000    0.000   24.955)   24.955
  14.620   (   0.000    0.000  -15.458)   15.458
  14.714   (  -0.000    0.000    2.699)    2.699
  14.714   (   0.000    0.000    2.699)    2.699
  14.885   (   0.000    0.000   -5.410)    5.410
  14.885   (   0.000    0.000   -5.410)    5.410
======================= Grid point 724 (83/135) =======================
q-point: ( 0.07  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.307   (  12.814    7.398   23.759)   27.989
   2.649   (  36.747   21.216   22.884)   48.209
   2.893   (   9.760    5.635   -5.953)   12.745
   3.226   (  35.147   20.292   -6.175)   41.051
   5.426   (   3.539    2.043   78.942)   79.047
   8.278   (  -6.477   -3.740  -64.039)   64.474
  13.722   (  -7.779   -4.491   25.126)   26.683
  14.537   (  -7.185   -4.148  -15.289)   17.395
  14.646   (  -5.643   -3.258    2.729)    7.064
  14.696   (  -1.252   -0.723    1.544)    2.115
  14.820   (  -5.390   -3.112   -5.521)    8.319
  14.842   (  -3.587   -2.071   -4.919)    6.431
======================= Grid point 725 (84/135) =======================
q-point: ( 0.13  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.690   (  18.853   10.885   19.692)   29.355
   3.192   (  14.961    8.638   -5.216)   18.046
   3.514   (  31.107   17.960   26.406)   44.582
   4.185   (  42.756   24.685   -7.987)   50.012
   5.667   (  20.643   11.918   62.314)   66.718
   8.103   (  -6.697   -3.867  -57.690)   58.206
  13.477   ( -12.731   -7.350   25.393)   29.341
  14.304   ( -12.202   -7.045  -15.221)   20.741
  14.466   (  -9.456   -5.459    2.818)   11.276
  14.648   (  -9.006   -5.200   -5.823)   11.919
  14.655   (  -2.486   -1.435   -0.935)    3.019
  14.723   (  -6.472   -3.737   -3.406)    8.213
======================= Grid point 726 (85/135) =======================
q-point: ( 0.20  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.122   (  17.174    9.915   16.005)   25.484
   3.537   (  13.678    7.897   -4.266)   16.360
   3.937   (   6.707    3.872   34.516)   35.374
   5.000   (  23.171   13.378  -16.586)   31.479
   6.459   (  45.362   26.190   38.764)   65.164
   8.103   (  10.247    5.916  -43.135)   44.728
  13.162   ( -13.808   -7.972   25.667)   30.216
  14.001   ( -13.304   -7.681  -15.703)   21.967
  14.232   ( -10.063   -5.810    2.963)   11.992
  14.427   (  -9.484   -5.476   -6.246)   12.607
  14.534   (  -8.720   -5.035   -0.573)   10.086
  14.591   (  -3.988   -2.303   -3.631)    5.865
======================= Grid point 727 (86/135) =======================
q-point: ( 0.27  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.446   (  10.091    5.826   13.327)   17.703
   3.785   (   7.044    4.067   -3.075)    8.696
   3.970   (  -1.315   -0.759   27.587)   27.629
   5.185   (  -5.077   -2.931  -28.215)   28.818
   7.611   (  51.127   29.518   28.393)   65.509
   8.648   (  34.702   20.035  -24.840)   47.145
  12.852   ( -12.748   -7.360   26.265)   30.108
  13.703   ( -12.316   -7.111  -15.987)   21.397
  14.025   (  -7.406   -4.276    3.159)    9.117
  14.233   (  -6.741   -3.892   -6.735)   10.293
  14.311   ( -10.021   -5.785    2.804)   11.906
  14.492   (  -4.455   -2.572   -6.088)    7.971
======================= Grid point 728 (87/135) =======================
q-point: ( 0.33  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.582   (   1.907    1.101   11.511)   11.720
   3.854   (  -0.869   -0.502   -1.516)    1.818
   3.942   (  -0.439   -0.253   13.636)   13.645
   4.933   ( -13.542   -7.818  -31.184)   34.885
   8.754   (  46.191   26.668   25.887)   59.287
   9.560   (  40.858   23.590  -15.212)   49.571
  12.567   ( -11.909   -6.876   26.836)   30.154
  13.412   ( -13.440   -7.760  -15.228)   21.743
  13.900   (  -3.349   -1.933    3.390)    5.142
  14.096   (  -7.770   -4.486    5.775)   10.670
  14.124   (  -2.664   -1.538   -7.253)    7.878
  14.392   (  -3.977   -2.296   -8.220)    9.415
======================= Grid point 729 (88/135) =======================
q-point: ( 0.40  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.563   (  -2.658   -1.535   10.345)   10.791
   3.775   (  -4.891   -2.824    0.256)    5.653
   3.956   (   1.428    0.824    1.124)    1.995
   4.645   ( -10.243   -5.914  -29.454)   31.741
   9.718   (  36.080   20.831   26.307)   49.272
  10.443   (  33.293   19.222  -10.015)   39.726
  12.305   (  -9.970   -5.756   24.899)   27.431
  13.053   ( -17.576  -10.148  -11.831)   23.492
  13.859   (  -0.480   -0.277    3.609)    3.651
  13.970   (  -3.171   -1.831    6.264)    7.255
  14.099   (   0.066    0.038   -7.725)    7.725
  14.315   (  -2.540   -1.466   -9.427)    9.872
======================= Grid point 730 (89/135) =======================
q-point: ( 0.47  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.502   (  -1.759   -1.016    9.752)    9.961
   3.676   (  -2.666   -1.539    1.550)    3.446
   3.991   (   1.064    0.614   -6.040)    6.163
   4.480   (  -3.697   -2.134  -27.376)   27.706
  10.367   (  17.268    9.970   27.497)   33.966
  10.989   (  11.194    6.463   -1.030)   12.966
  12.191   (   2.213    1.278   13.630)   13.868
  12.632   ( -16.212   -9.360   -6.338)   19.764
  13.860   (   0.185    0.107    3.745)    3.751
  13.932   (  -0.552   -0.319    5.669)    5.705
  14.108   (   0.396    0.229   -8.018)    8.031
  14.276   (  -0.839   -0.485   -9.733)    9.781
======================= Grid point 739 (90/135) =======================
q-point: ( 0.07  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.629   (  13.107   22.702   20.385)   33.207
   3.144   (  10.482   18.156   -5.303)   21.624
   3.268   (  19.585   33.923   24.109)   45.995
   3.883   (  23.520   40.738   -7.176)   47.585
   5.555   (   7.654   13.258   68.590)   70.278
   8.148   (  -4.598   -7.964  -60.032)   60.733
  13.554   (  -6.758  -11.705   25.335)   28.715
  14.365   (  -6.894  -11.940  -12.994)   18.945
  14.566   (  -3.500   -6.062   -1.018)    7.074
  14.635   (  -2.502   -4.334    1.205)    5.148
  14.713   (  -4.209   -7.291   -5.441)   10.024
  14.751   (  -3.882   -6.723   -4.037)    8.750
======================= Grid point 740 (91/135) =======================
q-point: ( 0.13  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.124   (  10.685   30.264   16.500)   36.088
   3.551   (   8.882   25.810   -4.384)   27.645
   3.796   (  13.726    8.048   32.081)   35.810
   4.690   (  27.510   22.226  -12.304)   37.446
   6.072   (  28.744   28.606   47.117)   62.165
   8.042   (   0.940   -0.301  -49.823)   49.833
  13.280   ( -10.436  -11.611   25.632)   30.012
  14.097   (  -9.241  -12.039  -13.016)   19.994
  14.377   ( -11.000   -4.532   -0.599)   11.912
  14.512   (  -6.939   -8.190   -4.749)   11.738
  14.568   (  -2.394   -6.757   -2.252)    7.514
  14.592   (  -4.409   -9.101   -2.520)   10.422
======================= Grid point 741 (92/135) =======================
q-point: ( 0.20  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.534   (   1.335   27.299   14.863)   31.111
   3.909   (   0.082   26.522   -3.033)   26.696
   3.964   (   3.862    2.516   29.224)   29.586
   5.129   (   8.786   -1.013  -23.500)   25.109
   7.051   (  45.314   31.186   31.920)   63.599
   8.307   (  22.880   14.462  -32.250)   42.104
  12.979   ( -11.349   -9.926   26.228)   30.253
  13.820   ( -10.043   -9.333  -14.602)   20.029
  14.132   ( -10.411   -5.299    1.295)   11.753
  14.290   (  -7.674   -8.619   -4.248)   12.298
  14.395   (  -6.803   -7.485   -1.242)   10.191
  14.479   (  -0.822   -8.436   -4.513)    9.603
======================= Grid point 742 (93/135) =======================
q-point: ( 0.27  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.698   (  -3.164   11.242   14.047)   18.268
   4.029   (  -4.669   14.769   13.530)   20.567
   4.098   (  -9.410   23.841    0.225)   25.632
   5.050   (  -8.119   -8.673  -29.753)   32.037
   8.172   (  48.723   25.496   26.711)   61.135
   9.058   (  39.659   19.691  -18.762)   48.089
  12.687   ( -11.234   -9.318   27.027)   30.716
  13.544   ( -11.782   -9.290  -15.435)   21.526
  13.943   (  -5.036   -4.355    3.092)    7.341
  14.092   (  -4.534   -9.065   -0.287)   10.139
  14.222   (  -6.285   -2.776   -2.985)    7.491
  14.379   (  -0.863   -8.283   -6.848)   10.782
======================= Grid point 743 (94/135) =======================
q-point: ( 0.33  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.721   (  -4.063    7.539    7.950)   11.686
   4.053   ( -11.861   17.671    1.887)   21.366
   4.098   ( -14.172   19.933    3.582)   24.718
   4.786   ( -12.272   -5.087  -28.613)   31.546
   9.210   (  43.925   18.497   25.890)   54.238
   9.963   (  40.896   16.024  -12.347)   45.626
  12.396   ( -10.415  -10.700   26.725)   30.613
  13.200   ( -14.253  -14.360  -12.926)   24.009
  13.869   (  -0.927   -0.974    3.439)    3.693
  13.959   (  -1.513   -6.511    3.708)    7.644
  14.134   (  -2.832    0.802   -6.215)    6.877
  14.291   (   0.009   -7.348   -8.494)   11.231
======================= Grid point 744 (95/135) =======================
q-point: ( 0.40  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.713   (  -8.114   12.156    5.364)   15.569
   3.979   ( -16.804   20.565   -2.129)   26.643
   4.072   (  -8.138   13.748    0.664)   15.990
   4.564   (  -8.908   -0.984  -26.040)   27.539
  10.033   (  32.086   10.456   26.763)   43.070
  10.712   (  28.178    7.337   -6.422)   29.818
  12.147   (  -3.365  -10.333   21.317)   23.927
  12.769   ( -15.357  -18.188   -8.629)   25.320
  13.862   (  -0.225    0.322    3.725)    3.745
  13.909   (   1.082   -3.148    4.685)    5.747
  14.120   (  -0.751    1.514   -7.803)    7.984
  14.236   (   1.436   -6.035   -9.168)   11.070
======================= Grid point 745 (96/135) =======================
q-point: ( 0.47  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.699   (  -9.101   15.764    5.606)   19.046
   3.930   ( -13.789   23.882   -1.269)   27.606
   4.075   (  -6.321   10.948   -3.616)   13.148
   4.477   (  -1.701    2.946  -24.897)   25.128
  10.412   (   3.361   -5.821   27.555)   28.363
  10.985   (   3.597   -6.230    1.945)    7.452
  12.132   (   5.185   -8.980    9.624)   14.147
  12.443   (   4.200   -7.275   -4.110)    9.352
  13.861   (   0.247   -0.429    4.124)    4.154
  13.909   (   0.815   -1.412    4.555)    4.838
  14.123   (  -0.597    1.034   -8.307)    8.392
  14.219   (   2.796   -4.843   -9.242)   10.802
======================= Grid point 756 (97/135) =======================
q-point: ( 0.13  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.781   (  18.027   31.224   12.844)   38.274
   3.848   (  -0.394   -0.682   31.272)   31.282
   4.123   (  16.119   27.919   -3.226)   32.399
   4.988   (   3.814    6.606  -21.161)   22.494
   6.787   (  23.803   41.227   34.651)   58.881
   8.177   (   8.847   15.323  -36.022)   40.133
  13.031   (  -7.418  -12.848   26.379)   30.265
  13.860   (  -6.453  -11.178  -12.806)   18.182
  14.223   (  -6.869  -11.898   -1.270)   13.797
  14.335   (  -5.304   -9.186   -4.154)   11.392
  14.388   (  -5.380   -9.318   -1.726)   10.897
  14.422   (  -4.077   -7.061   -3.610)    8.917
======================= Grid point 757 (98/135) =======================
q-point: ( 0.20  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.792   (  -2.451   -1.003   21.170)   21.335
   4.317   (   4.437   31.011   10.432)   33.018
   4.554   (  -4.475   32.406   -3.971)   32.954
   5.004   (   2.932   -8.919  -26.470)   28.086
   7.683   (  35.899   31.636   28.480)   55.684
   8.670   (  25.442   21.214  -22.495)   40.042
  12.763   (  -9.330  -11.803   27.519)   31.363
  13.627   (  -9.110  -11.172  -14.468)   20.424
  13.969   (  -6.301  -10.413    2.915)   12.515
  14.122   (  -8.206   -6.545   -2.789)   10.861
  14.254   (  -1.210   -8.399   -4.360)    9.540
  14.295   (  -2.316   -7.154   -4.592)    8.811
======================= Grid point 758 (99/135) =======================
q-point: ( 0.27  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.771   (   0.422    0.933    9.160)    9.217
   4.499   ( -11.709   15.989   -9.155)   21.830
   4.645   ( -14.523   28.531    1.501)   32.050
   4.990   (  -7.893   11.473  -15.402)   20.764
   8.632   (  41.663   20.698   26.172)   53.377
   9.420   (  36.780   16.445  -14.335)   42.763
  12.458   ( -10.895  -13.946   27.813)   32.966
  13.277   ( -13.541  -18.342  -11.964)   25.747
  13.850   (   0.075   -3.661    2.991)    4.728
  13.952   (  -4.644   -4.657   -0.241)    6.581
  14.172   (   0.008   -5.475   -6.381)    8.408
  14.197   (  -1.409   -5.660   -6.009)    8.374
======================= Grid point 759 (100/135) =======================
q-point: ( 0.33  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.819   (   1.639    4.064   -0.671)    4.433
   4.391   (  -9.499    5.278  -21.135)   23.765
   4.637   ( -21.526   32.104    7.357)   39.347
   4.913   ( -20.318   23.341   -8.609)   32.121
   9.508   (  38.813   11.855   26.008)   48.202
  10.206   (  36.546    8.733   -9.333)   38.717
  12.116   (  -9.163  -17.434   25.601)   32.300
  12.815   ( -14.547  -23.643   -8.659)   29.079
  13.850   (   1.245   -3.235    2.079)    4.042
  13.871   (  -0.492   -0.171    2.820)    2.868
  14.120   (   2.010   -8.867   -7.095)   11.532
  14.156   (  -0.836    1.315   -7.824)    7.977
======================= Grid point 760 (101/135) =======================
q-point: ( 0.40  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.882   (  -0.796    5.501   -5.868)    8.083
   4.312   (  -5.256    5.006  -22.558)   23.697
   4.581   ( -21.480   33.221    6.870)   40.152
   4.804   ( -21.183   29.541   -5.467)   36.760
  10.148   (  22.954    1.309   26.700)   35.235
  10.748   (  18.040   -3.197   -1.525)   18.385
  11.884   (   5.300  -15.670   16.355)   23.262
  12.346   (  -8.864  -22.626   -5.268)   24.864
  13.825   (   1.094   -3.245    3.135)    4.643
  13.897   (   0.709    1.131    3.591)    3.831
  14.080   (   3.170   -7.719   -7.923)   11.507
  14.164   (  -0.413    1.559   -8.476)    8.628
======================= Grid point 772 (102/135) =======================
q-point: ( 0.20  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.769   (   0.136    0.235   10.575)   10.578
   4.666   (  -3.277   -5.675  -24.309)   25.176
   4.967   (  10.006   17.331    2.422)   20.158
   5.163   (  12.566   21.765   -2.149)   25.224
   8.336   (  19.468   33.720   26.770)   47.251
   9.153   (  15.531   26.900  -15.505)   34.717
  12.485   (  -9.357  -16.207   28.306)   33.933
  13.303   ( -12.215  -21.157  -10.865)   26.737
  13.827   (  -1.613   -2.793    2.120)    3.860
  13.993   (  -4.019   -6.961   -4.887)    9.407
  14.091   (  -3.332   -5.771   -1.976)    6.950
  14.207   (  -1.557   -2.697   -6.532)    7.237
======================= Grid point 773 (103/135) =======================
q-point: ( 0.27  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.817   (   2.655    4.050    0.012)    4.842
   4.497   (  -5.666   -5.873  -29.313)   30.428
   5.228   (  -6.826   25.274    8.146)   27.418
   5.456   (  -5.928   24.727   -2.599)   25.560
   9.040   (  29.187   20.958   25.872)   44.277
   9.746   (  26.434   16.880  -10.434)   33.054
  12.110   ( -11.937  -20.713   27.158)   36.181
  12.818   ( -15.039  -26.341   -7.589)   31.267
  13.818   (   0.437   -0.116    0.042)    0.454
  13.875   (  -2.075   -3.486   -1.456)    4.310
  14.009   (   0.374   -6.580   -2.319)    6.986
  14.179   (  -0.465    0.471   -8.051)    8.078
======================= Grid point 774 (104/135) =======================
q-point: ( 0.33  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.921   (   2.387    6.170   -7.813)   10.238
   4.373   (  -3.325   -3.178  -27.404)   27.788
   5.273   ( -17.293   28.661    8.500)   34.536
   5.502   ( -17.887   29.243   -2.450)   34.367
   9.727   (  28.911   11.135   25.716)   40.263
  10.350   (  26.019    6.582   -5.443)   27.385
  11.701   (  -6.898  -23.170   22.935)   33.324
  12.281   ( -13.470  -28.543   -5.736)   32.079
  13.787   (  -0.641   -2.551    0.294)    2.647
  13.881   (   2.599   -0.335    0.389)    2.649
  13.961   (   0.994   -5.766   -2.292)    6.284
  14.189   (  -0.122    1.783   -8.912)    9.089
======================= Grid point 775 (105/135) =======================
q-point: ( 0.40  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.982   (  -2.708    4.690  -11.074)   12.327
   4.324   (   0.626   -1.084  -25.947)   25.977
   5.258   ( -17.974   31.132    8.496)   36.939
   5.484   ( -18.810   32.580   -2.211)   37.685
  10.104   (   2.856   -4.946   25.832)   26.456
  10.621   (   4.636   -8.029    2.368)    9.569
  11.536   (  10.229  -17.718   13.912)   24.740
  11.932   (   9.667  -16.744   -3.966)   19.737
  13.769   (   1.183   -2.050    1.148)    2.630
  13.919   (  -0.442    0.766    4.045)    4.140
  13.926   (   3.986   -6.903   -5.942)    9.943
  14.199   (  -0.895    1.550   -9.157)    9.330
======================= Grid point 789 (106/135) =======================
q-point: ( 0.27  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.931   (   3.930    6.808   -8.073)   11.268
   4.385   (  -2.828   -4.898  -27.916)   28.483
   5.598   (   6.640   11.500    8.590)   15.816
   5.835   (   7.133   12.355   -2.625)   14.505
   9.504   (  14.517   25.143   25.377)   38.561
  10.126   (  11.986   20.760   -6.358)   24.801
  11.647   ( -13.953  -24.168   24.699)   37.267
  12.252   ( -16.782  -29.067   -5.615)   34.030
  13.797   (  -1.717   -2.974   -1.885)    3.918
  13.825   (  -0.373   -0.647   -2.718)    2.818
  13.934   (  -0.693   -1.200    2.771)    3.098
  14.199   (   0.742    1.286   -9.084)    9.204
======================= Grid point 790 (107/135) =======================
q-point: ( 0.33  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.060   (   1.378    7.312  -14.119)   15.959
   4.300   (  -1.331   -4.160  -24.593)   24.977
   5.726   (  -6.125   15.491    8.930)   18.900
   5.976   (  -6.345   16.547   -2.936)   17.964
   9.998   (  13.903   16.295   24.712)   32.703
  10.483   (   7.793    4.817    2.870)    9.601
  11.254   (  -0.627  -17.491   15.791)   23.573
  11.705   (  -9.375  -27.631   -4.322)   29.497
  13.751   (  -0.905   -0.891   -0.816)    1.509
  13.814   (   1.948   -4.122   -3.646)    5.838
  13.931   (  -0.032    0.487    4.170)    4.198
  14.221   (   0.033    1.211   -9.534)    9.611
======================= Grid point 805 (108/135) =======================
q-point: ( 0.33  0.33  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.198   (  -3.035   -5.256  -19.836)   20.744
   4.198   (   3.035    5.256  -19.836)   20.744
   5.882   (  -0.000   -0.000    9.145)    9.145
   6.144   (  -0.000   -0.000   -3.284)    3.284
  10.358   (  -8.433  -14.607   20.845)   26.814
  10.358   (   8.433   14.607   20.845)   26.814
  11.236   (  -7.975  -13.814   -1.125)   15.990
  11.236   (   7.975   13.814   -1.125)   15.990
  13.756   (  -0.744   -1.288   -1.997)    2.490
  13.756   (   0.744    1.288   -1.997)    2.490
  13.937   (  -0.000   -0.000    4.296)    4.296
  14.233   (  -0.000   -0.000   -9.746)    9.746
======================= Grid point 964 (109/135) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 81
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.561   (   0.000   -0.000  -15.700)   15.700
   2.561   (   0.000   -0.000  -15.700)   15.700
   2.561   (   0.000   -0.000   15.700)   15.700
   2.561   (   0.000    0.000   15.700)   15.700
   6.963   (  -0.000   -0.000  -75.081)   75.081
   6.963   (   0.000   -0.000   75.081)   75.081
  14.265   (   0.000    0.000  -20.402)   20.402
  14.265   (   0.000    0.000   20.402)   20.402
  14.785   (   0.000    0.000   -4.456)    4.456
  14.785   (  -0.000    0.000   -4.456)    4.456
  14.785   (  -0.000    0.000    4.456)    4.456
  14.785   (  -0.000    0.000    4.456)    4.456
======================= Grid point 966 (110/135) =======================
q-point: ( 0.07  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.690   (  10.664    6.157  -14.800)   19.253
   2.690   (  10.664    6.157   14.800)   19.253
   3.023   (  35.356   20.413  -14.640)   43.371
   3.023   (  35.356   20.413   14.640)   43.371
   6.926   (  -2.578   -1.488  -72.417)   72.478
   6.926   (  -2.578   -1.488   72.417)   72.478
  14.177   (  -7.398   -4.271  -20.750)   22.439
  14.177   (  -7.398   -4.271   20.750)   22.439
  14.719   (  -5.551   -3.205   -4.535)    7.852
  14.719   (  -5.551   -3.205    4.535)    7.852
  14.753   (  -2.630   -1.519   -3.870)    4.920
  14.753   (  -2.630   -1.519    3.870)    4.920
======================= Grid point 968 (111/135) =======================
q-point: ( 0.13  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.015   (  16.218    9.363  -12.789)   22.677
   3.015   (  16.218    9.363   12.789)   22.677
   3.942   (  38.196   22.052  -16.843)   47.211
   3.942   (  38.196   22.052   16.843)   47.211
   6.911   (   3.837    2.215  -62.398)   62.555
   6.911   (   3.837    2.215   62.398)   62.555
  13.940   ( -12.437   -7.180  -21.215)   25.618
  13.940   ( -12.437   -7.180   21.215)   25.618
  14.541   (  -9.289   -5.363   -4.752)   11.731
  14.541   (  -9.289   -5.363    4.752)   11.731
  14.668   (  -4.616   -2.665   -1.961)    5.679
  14.668   (  -4.616   -2.665    1.961)    5.679
======================= Grid point 970 (112/135) =======================
q-point: ( 0.20  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.390   (  14.963    8.639  -10.662)   20.303
   3.390   (  14.963    8.639   10.662)   20.303
   4.562   (  12.350    7.130  -27.581)   31.050
   4.562   (  12.350    7.130   27.581)   31.050
   7.259   (  28.225   16.296  -41.852)   53.046
   7.259   (  28.225   16.296   41.852)   53.046
  13.631   ( -13.526   -7.809  -21.501)   26.575
  13.631   ( -13.526   -7.809   21.501)   26.575
  14.312   (  -9.839   -5.681   -5.067)   12.440
  14.312   (  -9.839   -5.681    5.067)   12.440
  14.541   (  -6.263   -3.616   -1.367)    7.360
  14.541   (  -6.263   -3.616    1.367)    7.360
======================= Grid point 972 (113/135) =======================
q-point: ( 0.27  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.668   (   8.385    4.841   -8.767)   13.061
   3.668   (   8.385    4.841    8.767)   13.061
   4.575   (  -7.577   -4.374  -32.152)   33.321
   4.575   (  -7.577   -4.374   32.152)   33.321
   8.150   (  44.944   25.948  -25.993)   58.042
   8.150   (  44.944   25.948   25.993)   58.042
  13.329   ( -12.350   -7.130  -21.729)   25.991
  13.329   ( -12.350   -7.130   21.729)   25.991
  14.110   (  -7.139   -4.122   -5.447)    9.881
  14.110   (  -7.139   -4.122    5.447)    9.881
  14.384   (  -6.995   -4.039   -4.650)    9.320
  14.384   (  -6.995   -4.039    4.650)    9.320
======================= Grid point 974 (114/135) =======================
q-point: ( 0.33  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.770   (   0.580    0.335   -7.050)    7.082
   3.770   (   0.580    0.335    7.050)    7.082
   4.349   ( -10.060   -5.808  -26.402)   28.845
   4.349   ( -10.060   -5.808   26.402)   28.845
   9.213   (  44.724   25.821  -20.287)   55.485
   9.213   (  44.724   25.821   20.287)   55.485
  13.048   ( -12.265   -7.081  -21.536)   25.776
  13.048   ( -12.265   -7.081   21.536)   25.776
  13.992   (  -3.069   -1.772   -5.863)    6.851
  13.992   (  -3.069   -1.772    5.863)    6.851
  14.231   (  -5.816   -3.358   -7.710)   10.225
  14.231   (  -5.816   -3.358    7.710)   10.225
======================= Grid point 976 (115/135) =======================
q-point: ( 0.40  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.724   (  -3.594   -2.075   -5.536)    6.919
   3.724   (  -3.594   -2.075    5.536)    6.919
   4.156   (  -6.180   -3.568  -18.535)   19.862
   4.156   (  -6.180   -3.568   18.535)   19.862
  10.165   (  36.081   20.831  -18.421)   45.553
  10.165   (  36.081   20.831   18.421)   45.553
  12.745   ( -13.804   -7.970  -19.195)   24.950
  12.745   ( -13.804   -7.970   19.195)   24.950
  13.957   (  -0.256   -0.148   -6.248)    6.255
  13.957   (  -0.256   -0.148    6.248)    6.255
  14.127   (  -3.121   -1.802   -9.186)    9.867
  14.127   (  -3.121   -1.802    9.186)    9.867
======================= Grid point 978 (116/135) =======================
q-point: ( 0.47  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.648   (  -2.115   -1.221   -4.564)    5.176
   3.648   (  -2.115   -1.221    4.564)    5.176
   4.064   (  -1.876   -1.083  -13.293)   13.468
   4.064   (  -1.876   -1.083   13.293)   13.468
  10.811   (  16.816    9.709  -16.441)   25.443
  10.811   (  16.816    9.709   16.441)   25.443
  12.452   (  -9.060   -5.231  -11.719)   15.709
  12.452   (  -9.060   -5.231   11.719)   15.709
  13.962   (   0.271    0.157   -6.486)    6.494
  13.962   (   0.271    0.157    6.486)    6.494
  14.083   (  -0.823   -0.475   -9.229)    9.278
  14.083   (  -0.823   -0.475    9.229)    9.278
======================= Grid point 996 (117/135) =======================
q-point: ( 0.07  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.964   (  11.318   19.604  -13.131)   26.169
   2.964   (  11.318   19.604   13.131)   26.169
   3.661   (  21.958   38.032  -15.520)   46.578
   3.661   (  21.958   38.032   15.520)   46.578
   6.896   (   0.085    0.148  -66.178)   66.178
   6.896   (   0.085    0.148   66.178)   66.178
  14.014   (  -6.585  -11.406  -21.023)   24.808
  14.014   (  -6.585  -11.406   21.023)   24.808
  14.611   (  -4.333   -7.506   -4.028)    9.557
  14.611   (  -4.333   -7.506    4.028)    9.557
  14.681   (  -3.061   -5.301   -2.807)    6.734
  14.681   (  -3.061   -5.301    2.807)    6.734
======================= Grid point 998 (118/135) =======================
q-point: ( 0.13  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.400   (   9.492   27.341  -10.984)   30.956
   3.400   (   9.492   27.341   10.984)   30.956
   4.345   (  20.493   14.660  -22.800)   33.981
   4.345   (  20.493   14.660   22.800)   33.981
   7.044   (  13.337   12.864  -50.520)   53.811
   7.044   (  13.337   12.864   50.520)   53.811
  13.747   ( -10.169  -11.295  -21.255)   26.130
  13.747   ( -10.169  -11.295   21.255)   26.130
  14.419   (  -8.860   -6.815   -4.000)   11.872
  14.419   (  -8.860   -6.815    4.000)   11.872
  14.558   (  -3.214   -8.035   -0.806)    8.692
  14.558   (  -3.214   -8.035    0.806)    8.692
======================= Grid point 1000 (119/135) =======================
q-point: ( 0.20  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.779   (   0.631   28.041   -9.619)   29.652
   3.779   (   0.631   28.041    9.619)   29.652
   4.581   (   2.560   -4.399  -30.913)   31.329
   4.581   (   2.560   -4.399   30.913)   31.329
   7.678   (  35.866   24.462  -31.559)   53.672
   7.678   (  35.866   24.462   31.559)   53.672
  13.456   ( -11.073   -9.190  -21.694)   26.033
  13.456   ( -11.073   -9.190   21.694)   26.033
  14.194   (  -8.481   -6.307   -4.542)   11.504
  14.194   (  -8.481   -6.307    4.542)   11.504
  14.411   (  -3.795   -8.833   -2.333)    9.893
  14.411   (  -3.795   -8.833    2.333)    9.893
======================= Grid point 1002 (120/135) =======================
q-point: ( 0.27  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.955   (  -7.816   19.159  -10.118)   23.033
   3.955   (  -7.816   19.159   10.118)   23.033
   4.475   (  -7.978   -1.677  -26.653)   27.872
   4.475   (  -7.978   -1.677   26.653)   27.872
   8.658   (  45.531   23.620  -22.217)   55.898
   8.658   (  45.531   23.620   22.217)   55.898
  13.173   ( -11.404   -8.736  -21.904)   26.194
  13.173   ( -11.404   -8.736   21.904)   26.194
  14.024   (  -4.272   -4.575   -4.851)    7.919
  14.024   (  -4.272   -4.575    4.851)    7.919
  14.257   (  -3.813   -7.858   -5.000)   10.064
  14.257   (  -3.813   -7.858    5.000)   10.064
======================= Grid point 1004 (121/135) =======================
q-point: ( 0.33  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.930   (  -9.782    9.531  -10.990)   17.530
   3.930   (  -9.782    9.531   10.990)   17.530
   4.322   ( -13.162    8.391  -17.004)   23.082
   4.322   ( -13.162    8.391   17.004)   23.082
   9.657   (  43.458   17.991  -18.986)   50.723
   9.657   (  43.458   17.991   18.986)   50.723
  12.867   ( -12.222  -11.680  -20.597)   26.647
  12.867   ( -12.222  -11.680   20.597)   26.647
  13.955   (  -0.280   -1.838   -4.654)    5.011
  13.955   (  -0.280   -1.838    4.654)    5.011
  14.133   (  -2.338   -5.789   -6.539)    9.041
  14.133   (  -2.338   -5.789    6.539)    9.041
======================= Grid point 1006 (122/135) =======================
q-point: ( 0.40  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.856   (  -8.943    9.555   -9.078)   15.927
   3.856   (  -8.943    9.555    9.078)   15.927
   4.195   ( -12.580   11.893  -11.185)   20.610
   4.195   ( -12.580   11.893   11.185)   20.610
  10.475   (  32.183    9.993  -17.446)   37.947
  10.475   (  32.183    9.993   17.446)   37.947
  12.524   ( -10.797  -14.355  -16.122)   24.136
  12.524   ( -10.797  -14.355   16.122)   24.136
  13.958   (   0.939    0.101   -5.068)    5.155
  13.958   (   0.939    0.101    5.068)    5.155
  14.062   (  -0.162   -4.204   -7.117)    8.267
  14.062   (  -0.162   -4.204    7.117)    8.267
======================= Grid point 1008 (123/135) =======================
q-point: ( 0.47  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.824   (  -7.164   12.409   -7.860)   16.343
   3.824   (  -7.164   12.409    7.860)   16.343
   4.143   (  -7.878   13.645   -9.191)   18.241
   4.143   (  -7.878   13.645    9.191)   18.241
  10.847   (   3.616   -6.263  -15.387)   17.002
  10.847   (   3.616   -6.263   15.387)   17.002
  12.318   (   4.359   -7.550   -8.379)   12.092
  12.318   (   4.359   -7.550    8.379)   12.092
  13.970   (  -0.228    0.394   -6.828)    6.843
  13.970   (  -0.228    0.394    6.828)    6.843
  14.039   (   1.975   -3.421   -8.344)    9.232
  14.039   (   1.975   -3.421    8.344)    9.232
======================= Grid point 1030 (124/135) =======================
q-point: ( 0.13  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.998   (  16.502   28.582   -8.660)   34.121
   3.998   (  16.502   28.582    8.660)   34.121
   4.476   (  -0.329   -0.569  -29.868)   29.875
   4.476   (  -0.329   -0.569   29.868)   29.875
   7.475   (  17.306   29.975  -35.019)   49.238
   7.475   (  17.306   29.975   35.019)   49.238
  13.510   (  -6.905  -11.959  -21.627)   25.659
  13.510   (  -6.905  -11.959   21.627)   25.659
  14.253   (  -5.808  -10.060   -3.515)   12.136
  14.253   (  -5.808  -10.060    3.515)   12.136
  14.380   (  -4.910   -8.504   -0.820)    9.853
  14.380   (  -4.910   -8.504    0.820)    9.853
======================= Grid point 1032 (125/135) =======================
q-point: ( 0.20  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.263   (  -9.241    3.764  -21.978)   24.137
   4.263   (  -9.241    3.764   21.978)   24.137
   4.579   (   6.612   18.890  -12.502)   23.598
   4.579   (   6.612   18.890   12.502)   23.598
   8.209   (  31.864   27.930  -24.667)   49.029
   8.209   (  31.864   27.930   24.667)   49.029
  13.260   (  -9.239  -10.927  -22.270)   26.471
  13.260   (  -9.239  -10.927   22.270)   26.471
  14.046   (  -5.609   -8.393   -4.521)   11.060
  14.046   (  -5.609   -8.393    4.521)   11.060
  14.229   (  -3.086   -7.982   -1.635)    8.712
  14.229   (  -3.086   -7.982    1.635)    8.712
======================= Grid point 1034 (126/135) =======================
q-point: ( 0.27  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.095   (  -6.832   -3.927  -22.023)   23.390
   4.095   (  -6.832   -3.927   22.023)   23.390
   4.788   ( -11.971   29.499   -8.199)   32.874
   4.788   ( -11.971   29.499    8.199)   32.874
   9.088   (  40.092   19.487  -19.760)   48.761
   9.088   (  40.092   19.487   19.760)   48.761
  12.947   ( -12.085  -14.699  -21.116)   28.425
  12.947   ( -12.085  -14.699   21.116)   28.425
  13.924   (  -0.483   -5.132   -3.449)    6.203
  13.924   (  -0.483   -5.132    3.449)    6.203
  14.112   (  -2.353   -5.593   -1.445)    6.238
  14.112   (  -2.353   -5.593    1.445)    6.238
======================= Grid point 1036 (127/135) =======================
q-point: ( 0.33  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.980   (  -3.410   -1.163  -16.345)   16.737
   3.980   (  -3.410   -1.163   16.345)   16.737
   4.783   ( -22.284   31.686   -6.577)   39.292
   4.783   ( -22.284   31.686    6.577)   39.292
   9.943   (  38.755   11.083  -17.630)   43.996
   9.943   (  38.755   11.083   17.630)   43.996
  12.552   ( -12.515  -19.894  -17.959)   29.579
  12.552   ( -12.515  -19.894   17.959)   29.579
  13.907   (   2.607   -3.495   -1.673)    4.670
  13.907   (   2.607   -3.495    1.673)    4.670
  14.039   (  -1.328   -2.502   -1.842)    3.379
  14.039   (  -1.328   -2.502    1.842)    3.379
======================= Grid point 1038 (128/135) =======================
q-point: ( 0.40  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.941   (  -1.537    1.257  -12.045)   12.208
   3.941   (  -1.537    1.257   12.045)   12.208
   4.710   ( -22.716   34.228   -5.509)   41.448
   4.710   ( -22.716   34.228    5.509)   41.448
  10.574   (  22.897    0.122  -15.686)   27.755
  10.574   (  22.897    0.122   15.686)   27.755
  12.175   (  -4.184  -19.551  -12.081)   23.360
  12.175   (  -4.184  -19.551   12.081)   23.360
  13.920   (   3.194   -4.643   -5.863)    8.132
  13.920   (   3.194   -4.643    5.863)    8.132
  14.012   (  -0.735    0.313   -6.035)    6.088
  14.012   (  -0.735    0.313    6.035)    6.088
======================= Grid point 1062 (129/135) =======================
q-point: ( 0.20  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.131   (  -4.799   -8.312  -24.771)   26.565
   4.131   (  -4.799   -8.312   24.771)   26.565
   5.085   (  12.622   21.863   -5.748)   25.891
   5.085   (  12.622   21.863    5.748)   25.891
   8.804   (  18.161   31.456  -20.412)   41.665
   8.804   (  18.161   31.456   20.412)   41.665
  12.981   ( -10.031  -17.374  -21.246)   29.221
  12.981   ( -10.031  -17.374   21.246)   29.221
  13.896   (  -3.098   -5.365   -4.559)    7.692
  13.896   (  -3.098   -5.365    4.559)    7.692
  14.105   (  -2.621   -4.539   -3.276)    6.181
  14.105   (  -2.621   -4.539    3.276)    6.181
======================= Grid point 1064 (130/135) =======================
q-point: ( 0.27  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.006   (  -2.563   -3.341  -18.490)   18.963
   4.006   (  -2.563   -3.341   18.490)   18.963
   5.372   (  -6.435   25.297   -5.963)   26.775
   5.372   (  -6.435   25.297    5.963)   26.775
   9.473   (  28.430   19.761  -17.670)   38.872
   9.473   (  28.430   19.761   17.670)   38.872
  12.565   ( -13.707  -23.054  -18.201)   32.414
  12.565   ( -13.707  -23.054   18.201)   32.414
  13.836   (   0.555   -3.425   -1.965)    3.987
  13.836   (   0.555   -3.425    1.965)    3.987
  14.042   (  -1.102   -1.567   -5.139)    5.484
  14.042   (  -1.102   -1.567    5.139)    5.484
======================= Grid point 1066 (131/135) =======================
q-point: ( 0.33  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.965   (  -0.785    0.140  -12.306)   12.331
   3.965   (  -0.785    0.140   12.306)   12.331
   5.420   ( -17.542   29.045   -5.973)   34.453
   5.420   ( -17.542   29.045    5.973)   34.453
  10.141   (  28.779    9.799  -15.794)   34.259
  10.141   (  28.779    9.799   15.794)   34.259
  12.081   ( -11.474  -26.177  -14.828)   32.199
  12.081   ( -11.474  -26.177   14.828)   32.199
  13.824   (   2.036   -4.073   -3.386)    5.674
  13.824   (   2.036   -4.073    3.386)    5.674
  14.029   (  -0.322    0.737   -6.794)    6.842
  14.029   (  -0.322    0.737    6.794)    6.842
======================= Grid point 1068 (132/135) =======================
q-point: ( 0.40  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.962   (  -0.623    1.079   -9.285)    9.368
   3.962   (  -0.623    1.079    9.285)    9.368
   5.404   ( -18.383   31.841   -5.901)   37.237
   5.404   ( -18.383   31.841    5.901)   37.237
  10.502   (   3.664   -6.346  -13.731)   15.564
  10.502   (   3.664   -6.346   13.731)   15.564
  11.791   (   9.841  -17.045  -10.507)   22.311
  11.791   (   9.841  -17.045   10.507)   22.311
  13.821   (   2.540   -4.399   -4.189)    6.585
  13.821   (   2.540   -4.399    4.189)    6.585
  14.032   (  -0.625    1.082   -7.305)    7.411
  14.032   (  -0.625    1.082    7.305)    7.411
======================= Grid point 1096 (133/135) =======================
q-point: ( 0.27  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.973   (  -0.096   -0.166  -12.342)   12.344
   3.973   (  -0.096   -0.166   12.342)   12.344
   5.748   (   6.877   11.912   -6.204)   15.089
   5.748   (   6.877   11.912    6.204)   15.089
   9.912   (  13.805   23.911  -15.650)   31.737
   9.912   (  13.805   23.911   15.650)   31.737
  12.049   ( -15.696  -27.186  -15.395)   34.964
  12.049   ( -15.696  -27.186   15.395)   34.964
  13.784   (  -1.095   -1.896   -0.762)    2.318
  13.784   (  -1.095   -1.896    0.762)    2.318
  14.036   (   0.288    0.498   -7.066)    7.089
  14.036   (   0.288    0.498    7.066)    7.089
======================= Grid point 1098 (134/135) =======================
q-point: ( 0.33  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.981   (  -0.065    1.063   -6.511)    6.598
   3.981   (  -0.065    1.063    6.511)    6.598
   5.882   (  -6.182   15.951   -6.535)   18.313
   5.882   (  -6.182   15.951    6.535)   18.313
  10.375   (  13.023   13.579  -12.868)   22.794
  10.375   (  13.023   13.579   12.868)   22.794
  11.546   (  -7.195  -25.396  -12.014)   29.001
  11.546   (  -7.195  -25.396   12.014)   29.001
  13.758   (   0.596   -2.298   -1.651)    2.892
  13.758   (   0.596   -2.298    1.651)    2.892
  14.048   (   0.016    0.834   -7.576)    7.622
  14.048   (   0.016    0.834    7.576)    7.622
======================= Grid point 1128 (135/135) =======================
q-point: ( 0.33  0.33 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 2.01e-04 2.01e-04 
Number of triplets: 119
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.994   (  -0.000   -0.000    0.000)    0.000
   3.994   (  -0.000   -0.000    0.000)    0.000
   6.044   (  -0.000   -0.000   -6.838)    6.838
   6.044   (  -0.000   -0.000    6.838)    6.838
  10.582   (  -0.000   -0.000   -0.000)    0.000
  10.582   (  -0.000   -0.000    0.000)    0.000
  11.213   (  -0.000   -0.000   -0.000)    0.000
  11.213   (  -0.000   -0.000    0.000)    0.000
  13.735   (  -0.000   -0.000   -0.000)    0.000
  13.735   (  -0.000   -0.000   -0.000)    0.000
  14.057   (  -0.000   -0.000   -7.738)    7.738
  14.057   (  -0.000   -0.000    7.738)    7.738
=================== 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/21600
   10.0   3452.091   3452.091   4369.727      0.000     -0.000     -0.000 3/21600
   20.0   2859.121   2859.121   3130.548      0.000     -0.000     -0.000 3/21600
   30.0   1588.172   1588.172   1480.379      0.000     -0.000     -0.000 3/21600
   40.0   1157.174   1157.174   1085.061     -0.000     -0.000     -0.000 3/21600
   50.0    917.723    917.723    860.692     -0.000     -0.000     -0.000 3/21600
   60.0    752.655    752.655    699.136     -0.000     -0.000     -0.000 3/21600
   70.0    625.564    625.564    573.075     -0.000     -0.000     -0.000 3/21600
   80.0    524.019    524.019    472.808     -0.000     -0.000     -0.000 3/21600
   90.0    442.620    442.620    393.606     -0.000     -0.000     -0.000 3/21600
  100.0    377.648    377.648    331.615     -0.000     -0.000     -0.000 3/21600
  110.0    325.872    325.872    283.223     -0.000     -0.000     -0.000 3/21600
  120.0    284.485    284.485    245.286     -0.000     -0.000     -0.000 3/21600
  130.0    251.174    251.174    215.274     -0.000     -0.000     -0.000 3/21600
  140.0    224.115    224.115    191.249     -0.000     -0.000     -0.000 3/21600
  150.0    201.901    201.901    171.766     -0.000     -0.000     -0.000 3/21600
  160.0    183.465    183.465    155.757     -0.000     -0.000     -0.000 3/21600
  170.0    167.997    167.997    142.435     -0.000     -0.000     -0.000 3/21600
  180.0    154.881    154.881    131.212     -0.000     -0.000     -0.000 3/21600
  190.0    143.651    143.651    121.653     -0.000     -0.000     -0.000 3/21600
  200.0    133.946    133.946    113.425     -0.000     -0.000     -0.000 3/21600
  210.0    125.487    125.487    106.276     -0.000     -0.000     -0.000 3/21600
  220.0    118.055    118.055    100.011     -0.000     -0.000     -0.000 3/21600
  230.0    111.478    111.478     94.477     -0.000     -0.000     -0.000 3/21600
  240.0    105.620    105.620     89.553     -0.000     -0.000     -0.000 3/21600
  250.0    100.369    100.369     85.144     -0.000     -0.000     -0.000 3/21600
  260.0     95.636     95.636     81.173     -0.000     -0.000     -0.000 3/21600
  270.0     91.349     91.349     77.576     -0.000     -0.000     -0.000 3/21600
  280.0     87.447     87.447     74.303     -0.000     -0.000     -0.000 3/21600
  290.0     83.880     83.880     71.310     -0.000     -0.000     -0.000 3/21600
  300.0     80.606     80.606     68.563     -0.000     -0.000     -0.000 3/21600
  310.0     77.590     77.590     66.032     -0.000     -0.000     -0.000 3/21600
  320.0     74.803     74.803     63.691     -0.000     -0.000     -0.000 3/21600
  330.0     72.218     72.218     61.520     -0.000     -0.000     -0.000 3/21600
  340.0     69.814     69.814     59.499     -0.000     -0.000     -0.000 3/21600
  350.0     67.573     67.573     57.615     -0.000     -0.000     -0.000 3/21600
  360.0     65.478     65.478     55.852     -0.000     -0.000     -0.000 3/21600
  370.0     63.514     63.514     54.199     -0.000     -0.000     -0.000 3/21600
  380.0     61.671     61.671     52.645     -0.000     -0.000     -0.000 3/21600
  390.0     59.936     59.936     51.183     -0.000     -0.000     -0.000 3/21600
  400.0     58.300     58.300     49.803     -0.000     -0.000     -0.000 3/21600
  410.0     56.755     56.755     48.499     -0.000     -0.000     -0.000 3/21600
  420.0     55.293     55.293     47.265     -0.000     -0.000     -0.000 3/21600
  430.0     53.907     53.907     46.094     -0.000     -0.000     -0.000 3/21600
  440.0     52.592     52.592     44.982     -0.000     -0.000     -0.000 3/21600
  450.0     51.342     51.342     43.925     -0.000     -0.000     -0.000 3/21600
  460.0     50.152     50.152     42.918     -0.000     -0.000     -0.000 3/21600
  470.0     49.018     49.018     41.958     -0.000     -0.000     -0.000 3/21600
  480.0     47.936     47.936     41.042     -0.000     -0.000     -0.000 3/21600
  490.0     46.902     46.902     40.166     -0.000     -0.000     -0.000 3/21600
  500.0     45.914     45.914     39.328     -0.000     -0.000     -0.000 3/21600
  510.0     44.968     44.968     38.525     -0.000     -0.000     -0.000 3/21600
  520.0     44.061     44.061     37.755     -0.000     -0.000     -0.000 3/21600
  530.0     43.191     43.191     37.017     -0.000     -0.000     -0.000 3/21600
  540.0     42.356     42.356     36.308     -0.000     -0.000     -0.000 3/21600
  550.0     41.553     41.553     35.626     -0.000     -0.000     -0.000 3/21600
  560.0     40.781     40.781     34.970     -0.000     -0.000     -0.000 3/21600
  570.0     40.038     40.038     34.338     -0.000     -0.000     -0.000 3/21600
  580.0     39.323     39.323     33.729     -0.000     -0.000     -0.000 3/21600
  590.0     38.633     38.633     33.142     -0.000     -0.000     -0.000 3/21600
  600.0     37.968     37.968     32.576     -0.000     -0.000     -0.000 3/21600
  610.0     37.326     37.326     32.029     -0.000     -0.000     -0.000 3/21600
  620.0     36.705     36.705     31.501     -0.000     -0.000     -0.000 3/21600
  630.0     36.106     36.106     30.990     -0.000     -0.000     -0.000 3/21600
  640.0     35.526     35.526     30.496     -0.000     -0.000     -0.000 3/21600
  650.0     34.965     34.965     30.017     -0.000     -0.000     -0.000 3/21600
  660.0     34.421     34.421     29.554     -0.000     -0.000     -0.000 3/21600
  670.0     33.895     33.895     29.105     -0.000     -0.000     -0.000 3/21600
  680.0     33.385     33.385     28.670     -0.000     -0.000     -0.000 3/21600
  690.0     32.890     32.890     28.248     -0.000     -0.000     -0.000 3/21600
  700.0     32.410     32.410     27.838     -0.000     -0.000     -0.000 3/21600
  710.0     31.945     31.945     27.441     -0.000     -0.000     -0.000 3/21600
  720.0     31.492     31.492     27.054     -0.000     -0.000     -0.000 3/21600
  730.0     31.053     31.053     26.679     -0.000     -0.000     -0.000 3/21600
  740.0     30.626     30.626     26.314     -0.000     -0.000     -0.000 3/21600
  750.0     30.210     30.210     25.959     -0.000     -0.000     -0.000 3/21600
  760.0     29.806     29.806     25.614     -0.000     -0.000     -0.000 3/21600
  770.0     29.413     29.413     25.278     -0.000     -0.000     -0.000 3/21600
  780.0     29.030     29.030     24.950     -0.000     -0.000     -0.000 3/21600
  790.0     28.657     28.657     24.632     -0.000     -0.000     -0.000 3/21600
  800.0     28.294     28.294     24.321     -0.000     -0.000     -0.000 3/21600
  810.0     27.940     27.940     24.018     -0.000     -0.000     -0.000 3/21600
  820.0     27.595     27.595     23.723     -0.000     -0.000     -0.000 3/21600
  830.0     27.258     27.258     23.435     -0.000     -0.000     -0.000 3/21600
  840.0     26.930     26.930     23.154     -0.000     -0.000     -0.000 3/21600
  850.0     26.610     26.610     22.880     -0.000     -0.000     -0.000 3/21600
  860.0     26.297     26.297     22.612     -0.000     -0.000     -0.000 3/21600
  870.0     25.992     25.992     22.351     -0.000     -0.000     -0.000 3/21600
  880.0     25.693     25.693     22.095     -0.000     -0.000     -0.000 3/21600
  890.0     25.402     25.402     21.846     -0.000     -0.000     -0.000 3/21600
  900.0     25.117     25.117     21.602     -0.000     -0.000     -0.000 3/21600
  910.0     24.839     24.839     21.363     -0.000     -0.000     -0.000 3/21600
  920.0     24.566     24.566     21.130     -0.000     -0.000     -0.000 3/21600
  930.0     24.300     24.300     20.902     -0.000     -0.000     -0.000 3/21600
  940.0     24.040     24.040     20.679     -0.000     -0.000     -0.000 3/21600
  950.0     23.785     23.785     20.460     -0.000     -0.000     -0.000 3/21600
  960.0     23.535     23.535     20.246     -0.000     -0.000     -0.000 3/21600
  970.0     23.291     23.291     20.037     -0.000     -0.000     -0.000 3/21600
  980.0     23.052     23.052     19.832     -0.000     -0.000     -0.000 3/21600
  990.0     22.818     22.818     19.631     -0.000     -0.000     -0.000 3/21600
 1000.0     22.588     22.588     19.434     -0.000     -0.000     -0.000 3/21600

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

