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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [4 4 2]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P-3m1 (164)
Number of symmetry operations in supercell: 384
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.050840731741600    0.000000000000000    0.000000000000000
  b   -2.025420365870800    3.508130980372969    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.128222920000000
Atomic positions (fractional):
   *1 Cu  0.66666666666667  0.33333333333333  0.59368810756266  63.546
    2 Cu  0.33333333333333  0.66666666666667  0.40631189243734  63.546
   *3 I   0.33333333333333  0.66666666666667  0.76578333646945 126.904
    4 I   0.66666666666667  0.33333333333333  0.23421666353055 126.904
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.050840731741600    0.000000000000000    0.000000000000000
  b   -2.025420365870800    3.508130980372969    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.128222920000000
Atomic positions (fractional):
   *1 Cu  0.66666666666667  0.33333333333333  0.59368810756266  63.546 > 1
    2 Cu  0.33333333333333  0.66666666666667  0.40631189243734  63.546 > 2
   *3 I   0.33333333333333  0.66666666666667  0.76578333646945 126.904 > 3
    4 I   0.66666666666667  0.33333333333333  0.23421666353055 126.904 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   16.203362926966399    0.000000000000000    0.000000000000000
  b   -8.101681463483200   14.032523921491876    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.256445840000000
Atomic positions (fractional):
   *1 Cu  0.16666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    2 Cu  0.41666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    3 Cu  0.66666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    4 Cu  0.91666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    5 Cu  0.16666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    6 Cu  0.41666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    7 Cu  0.66666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    8 Cu  0.91666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    9 Cu  0.16666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   10 Cu  0.41666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   11 Cu  0.66666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   12 Cu  0.91666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   13 Cu  0.16666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   14 Cu  0.41666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   15 Cu  0.66666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   16 Cu  0.91666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   17 Cu  0.16666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   18 Cu  0.41666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   19 Cu  0.66666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   20 Cu  0.91666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   21 Cu  0.16666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   22 Cu  0.41666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   23 Cu  0.66666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   24 Cu  0.91666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   25 Cu  0.16666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   26 Cu  0.41666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   27 Cu  0.66666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   28 Cu  0.91666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   29 Cu  0.16666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   30 Cu  0.41666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   31 Cu  0.66666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   32 Cu  0.91666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   33 Cu  0.08333333333333  0.16666666666667  0.20315594621867  63.546 > 2
   34 Cu  0.33333333333333  0.16666666666667  0.20315594621867  63.546 > 2
   35 Cu  0.58333333333333  0.16666666666667  0.20315594621867  63.546 > 2
   36 Cu  0.83333333333333  0.16666666666667  0.20315594621867  63.546 > 2
   37 Cu  0.08333333333333  0.41666666666667  0.20315594621867  63.546 > 2
   38 Cu  0.33333333333333  0.41666666666667  0.20315594621867  63.546 > 2
   39 Cu  0.58333333333333  0.41666666666667  0.20315594621867  63.546 > 2
   40 Cu  0.83333333333333  0.41666666666667  0.20315594621867  63.546 > 2
   41 Cu  0.08333333333333  0.66666666666667  0.20315594621867  63.546 > 2
   42 Cu  0.33333333333333  0.66666666666667  0.20315594621867  63.546 > 2
   43 Cu  0.58333333333333  0.66666666666667  0.20315594621867  63.546 > 2
   44 Cu  0.83333333333333  0.66666666666667  0.20315594621867  63.546 > 2
   45 Cu  0.08333333333333  0.91666666666667  0.20315594621867  63.546 > 2
   46 Cu  0.33333333333333  0.91666666666667  0.20315594621867  63.546 > 2
   47 Cu  0.58333333333333  0.91666666666667  0.20315594621867  63.546 > 2
   48 Cu  0.83333333333333  0.91666666666667  0.20315594621867  63.546 > 2
   49 Cu  0.08333333333333  0.16666666666667  0.70315594621867  63.546 > 2
   50 Cu  0.33333333333333  0.16666666666667  0.70315594621867  63.546 > 2
   51 Cu  0.58333333333333  0.16666666666667  0.70315594621867  63.546 > 2
   52 Cu  0.83333333333333  0.16666666666667  0.70315594621867  63.546 > 2
   53 Cu  0.08333333333333  0.41666666666667  0.70315594621867  63.546 > 2
   54 Cu  0.33333333333333  0.41666666666667  0.70315594621867  63.546 > 2
   55 Cu  0.58333333333333  0.41666666666667  0.70315594621867  63.546 > 2
   56 Cu  0.83333333333333  0.41666666666667  0.70315594621867  63.546 > 2
   57 Cu  0.08333333333333  0.66666666666667  0.70315594621867  63.546 > 2
   58 Cu  0.33333333333333  0.66666666666667  0.70315594621867  63.546 > 2
   59 Cu  0.58333333333333  0.66666666666667  0.70315594621867  63.546 > 2
   60 Cu  0.83333333333333  0.66666666666667  0.70315594621867  63.546 > 2
   61 Cu  0.08333333333333  0.91666666666667  0.70315594621867  63.546 > 2
   62 Cu  0.33333333333333  0.91666666666667  0.70315594621867  63.546 > 2
   63 Cu  0.58333333333333  0.91666666666667  0.70315594621867  63.546 > 2
   64 Cu  0.83333333333333  0.91666666666667  0.70315594621867  63.546 > 2
  *65 I   0.08333333333333  0.16666666666667  0.38289166823473 126.904 > 3
   66 I   0.33333333333333  0.16666666666667  0.38289166823473 126.904 > 3
   67 I   0.58333333333333  0.16666666666667  0.38289166823473 126.904 > 3
   68 I   0.83333333333333  0.16666666666667  0.38289166823473 126.904 > 3
   69 I   0.08333333333333  0.41666666666667  0.38289166823473 126.904 > 3
   70 I   0.33333333333333  0.41666666666667  0.38289166823473 126.904 > 3
   71 I   0.58333333333333  0.41666666666667  0.38289166823473 126.904 > 3
   72 I   0.83333333333333  0.41666666666667  0.38289166823473 126.904 > 3
   73 I   0.08333333333333  0.66666666666667  0.38289166823473 126.904 > 3
   74 I   0.33333333333333  0.66666666666667  0.38289166823473 126.904 > 3
   75 I   0.58333333333333  0.66666666666667  0.38289166823473 126.904 > 3
   76 I   0.83333333333333  0.66666666666667  0.38289166823473 126.904 > 3
   77 I   0.08333333333333  0.91666666666667  0.38289166823473 126.904 > 3
   78 I   0.33333333333333  0.91666666666667  0.38289166823473 126.904 > 3
   79 I   0.58333333333333  0.91666666666667  0.38289166823473 126.904 > 3
   80 I   0.83333333333333  0.91666666666667  0.38289166823473 126.904 > 3
   81 I   0.08333333333333  0.16666666666667  0.88289166823473 126.904 > 3
   82 I   0.33333333333333  0.16666666666667  0.88289166823473 126.904 > 3
   83 I   0.58333333333333  0.16666666666667  0.88289166823473 126.904 > 3
   84 I   0.83333333333333  0.16666666666667  0.88289166823473 126.904 > 3
   85 I   0.08333333333333  0.41666666666667  0.88289166823473 126.904 > 3
   86 I   0.33333333333333  0.41666666666667  0.88289166823473 126.904 > 3
   87 I   0.58333333333333  0.41666666666667  0.88289166823473 126.904 > 3
   88 I   0.83333333333333  0.41666666666667  0.88289166823473 126.904 > 3
   89 I   0.08333333333333  0.66666666666667  0.88289166823473 126.904 > 3
   90 I   0.33333333333333  0.66666666666667  0.88289166823473 126.904 > 3
   91 I   0.58333333333333  0.66666666666667  0.88289166823473 126.904 > 3
   92 I   0.83333333333333  0.66666666666667  0.88289166823473 126.904 > 3
   93 I   0.08333333333333  0.91666666666667  0.88289166823473 126.904 > 3
   94 I   0.33333333333333  0.91666666666667  0.88289166823473 126.904 > 3
   95 I   0.58333333333333  0.91666666666667  0.88289166823473 126.904 > 3
   96 I   0.83333333333333  0.91666666666667  0.88289166823473 126.904 > 3
   97 I   0.16666666666667  0.08333333333333  0.11710833176527 126.904 > 4
   98 I   0.41666666666667  0.08333333333333  0.11710833176527 126.904 > 4
   99 I   0.66666666666667  0.08333333333333  0.11710833176527 126.904 > 4
  100 I   0.91666666666667  0.08333333333333  0.11710833176527 126.904 > 4
  101 I   0.16666666666667  0.33333333333333  0.11710833176527 126.904 > 4
  102 I   0.41666666666667  0.33333333333333  0.11710833176527 126.904 > 4
  103 I   0.66666666666667  0.33333333333333  0.11710833176527 126.904 > 4
  104 I   0.91666666666667  0.33333333333333  0.11710833176527 126.904 > 4
  105 I   0.16666666666667  0.58333333333333  0.11710833176527 126.904 > 4
  106 I   0.41666666666667  0.58333333333333  0.11710833176527 126.904 > 4
  107 I   0.66666666666667  0.58333333333333  0.11710833176527 126.904 > 4
  108 I   0.91666666666667  0.58333333333333  0.11710833176527 126.904 > 4
  109 I   0.16666666666667  0.83333333333333  0.11710833176527 126.904 > 4
  110 I   0.41666666666667  0.83333333333333  0.11710833176527 126.904 > 4
  111 I   0.66666666666667  0.83333333333333  0.11710833176527 126.904 > 4
  112 I   0.91666666666667  0.83333333333333  0.11710833176527 126.904 > 4
  113 I   0.16666666666667  0.08333333333333  0.61710833176527 126.904 > 4
  114 I   0.41666666666667  0.08333333333333  0.61710833176527 126.904 > 4
  115 I   0.66666666666667  0.08333333333333  0.61710833176527 126.904 > 4
  116 I   0.91666666666667  0.08333333333333  0.61710833176527 126.904 > 4
  117 I   0.16666666666667  0.33333333333333  0.61710833176527 126.904 > 4
  118 I   0.41666666666667  0.33333333333333  0.61710833176527 126.904 > 4
  119 I   0.66666666666667  0.33333333333333  0.61710833176527 126.904 > 4
  120 I   0.91666666666667  0.33333333333333  0.61710833176527 126.904 > 4
  121 I   0.16666666666667  0.58333333333333  0.61710833176527 126.904 > 4
  122 I   0.41666666666667  0.58333333333333  0.61710833176527 126.904 > 4
  123 I   0.66666666666667  0.58333333333333  0.61710833176527 126.904 > 4
  124 I   0.91666666666667  0.58333333333333  0.61710833176527 126.904 > 4
  125 I   0.16666666666667  0.83333333333333  0.61710833176527 126.904 > 4
  126 I   0.41666666666667  0.83333333333333  0.61710833176527 126.904 > 4
  127 I   0.66666666666667  0.83333333333333  0.61710833176527 126.904 > 4
  128 I   0.91666666666667  0.83333333333333  0.61710833176527 126.904 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            6.7156849    0.0000000    0.0000000
            0.0000000    6.7156849    0.0000000
            0.0000000    0.0000000    5.9496968
-------------------------- Born effective charges --------------------------
    1 Cu    1.0552415    0.0000000    0.0000000
            0.0000000    1.0552415    0.0000000
            0.0000000    0.0000000    0.7420857
    2 Cu    1.0552415    0.0000000    0.0000000
            0.0000000    1.0552415    0.0000000
            0.0000000    0.0000000    0.7420857
    3 I    -1.0552415    0.0000000    0.0000000
            0.0000000   -1.0552415    0.0000000
            0.0000000    0.0000000   -0.7420857
    4 I    -1.0552415    0.0000000    0.0000000
            0.0000000   -1.0552415    0.0000000
            0.0000000    0.0000000   -0.7420857
----------------------------------------------------------------------------
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: 252
Number of blocks in projector: 200
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: 116
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 56
Use standard eigh solver.
Tree of FC basis block matrices:
- (252, 248), data: False
|-- (56, 54), data: True
|-- (80, 80), data: True
|-- (116, 114), data: True
-----
Solver_atoms: 1 -- 128 / 128
Time (Solver_compr_matrix_reshape): 0.010
Solver_block: 80 / 80
 - Time: 0.198
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.208
--------------------------------- 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: 252
Number of blocks in projector: 200
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: 116
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 56
Use standard eigh solver.
Tree of FC basis block matrices:
- (252, 248), data: False
|-- (56, 54), data: True
|-- (80, 80), data: True
|-- (116, 114), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 18:48:09]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P-3m1 (164)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.050840731741600    0.000000000000000    0.000000000000000
  b   -2.025420365870800    3.508130980372969    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.128222920000000
Atomic positions (fractional):
    1 Cu  0.66666666666667  0.33333333333333  0.59368810756266  63.546
    2 Cu  0.33333333333333  0.66666666666667  0.40631189243734  63.546
    3 I   0.33333333333333  0.66666666666667  0.76578333646945 126.904
    4 I   0.66666666666667  0.33333333333333  0.23421666353055 126.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.203362926966399    0.000000000000000    0.000000000000000
  b   -8.101681463483200   14.032523921491876    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.256445840000000
Atomic positions (fractional):
    1 Cu  0.16666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    2 Cu  0.41666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    3 Cu  0.66666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    4 Cu  0.91666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    5 Cu  0.16666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    6 Cu  0.41666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    7 Cu  0.66666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    8 Cu  0.91666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    9 Cu  0.16666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   10 Cu  0.41666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   11 Cu  0.66666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   12 Cu  0.91666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   13 Cu  0.16666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   14 Cu  0.41666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   15 Cu  0.66666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   16 Cu  0.91666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   17 Cu  0.16666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   18 Cu  0.41666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   19 Cu  0.66666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   20 Cu  0.91666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   21 Cu  0.16666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   22 Cu  0.41666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   23 Cu  0.66666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   24 Cu  0.91666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   25 Cu  0.16666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   26 Cu  0.41666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   27 Cu  0.66666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   28 Cu  0.91666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   29 Cu  0.16666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   30 Cu  0.41666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   31 Cu  0.66666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   32 Cu  0.91666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   33 Cu  0.08333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   34 Cu  0.33333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   35 Cu  0.58333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   36 Cu  0.83333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   37 Cu  0.08333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   38 Cu  0.33333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   39 Cu  0.58333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   40 Cu  0.83333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   41 Cu  0.08333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   42 Cu  0.33333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   43 Cu  0.58333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   44 Cu  0.83333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   45 Cu  0.08333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   46 Cu  0.33333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   47 Cu  0.58333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   48 Cu  0.83333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   49 Cu  0.08333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   50 Cu  0.33333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   51 Cu  0.58333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   52 Cu  0.83333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   53 Cu  0.08333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   54 Cu  0.33333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   55 Cu  0.58333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   56 Cu  0.83333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   57 Cu  0.08333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   58 Cu  0.33333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   59 Cu  0.58333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   60 Cu  0.83333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   61 Cu  0.08333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   62 Cu  0.33333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   63 Cu  0.58333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   64 Cu  0.83333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   65 I   0.08333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   66 I   0.33333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   67 I   0.58333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   68 I   0.83333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   69 I   0.08333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   70 I   0.33333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   71 I   0.58333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   72 I   0.83333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   73 I   0.08333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   74 I   0.33333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   75 I   0.58333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   76 I   0.83333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   77 I   0.08333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   78 I   0.33333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   79 I   0.58333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   80 I   0.83333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   81 I   0.08333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   82 I   0.33333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   83 I   0.58333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   84 I   0.83333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   85 I   0.08333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   86 I   0.33333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   87 I   0.58333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   88 I   0.83333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   89 I   0.08333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   90 I   0.33333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   91 I   0.58333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   92 I   0.83333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   93 I   0.08333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   94 I   0.33333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   95 I   0.58333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   96 I   0.83333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   97 I   0.16666666666667  0.08333333333333  0.11710833176527 126.904 > 97
   98 I   0.41666666666667  0.08333333333333  0.11710833176527 126.904 > 97
   99 I   0.66666666666667  0.08333333333333  0.11710833176527 126.904 > 97
  100 I   0.91666666666667  0.08333333333333  0.11710833176527 126.904 > 97
  101 I   0.16666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  102 I   0.41666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  103 I   0.66666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  104 I   0.91666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  105 I   0.16666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  106 I   0.41666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  107 I   0.66666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  108 I   0.91666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  109 I   0.16666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  110 I   0.41666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  111 I   0.66666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  112 I   0.91666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  113 I   0.16666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  114 I   0.41666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  115 I   0.66666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  116 I   0.91666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  117 I   0.16666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  118 I   0.41666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  119 I   0.66666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  120 I   0.91666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  121 I   0.16666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  122 I   0.41666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  123 I   0.66666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  124 I   0.91666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  125 I   0.16666666666667  0.83333333333333  0.61710833176527 126.904 > 97
  126 I   0.41666666666667  0.83333333333333  0.61710833176527 126.904 > 97
  127 I   0.66666666666667  0.83333333333333  0.61710833176527 126.904 > 97
  128 I   0.91666666666667  0.83333333333333  0.61710833176527 126.904 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            6.7156849    0.0000000    0.0000000
            0.0000000    6.7156849    0.0000000
            0.0000000    0.0000000    5.9496968
-------------------------- Born effective charges --------------------------
    1 Cu    1.0552415    0.0000000    0.0000000
            0.0000000    1.0552415    0.0000000
            0.0000000    0.0000000    0.7420857
    2 Cu    1.0552415    0.0000000    0.0000000
            0.0000000    1.0552415    0.0000000
            0.0000000    0.0000000    0.7420857
    3 I    -1.0552415    0.0000000    0.0000000
            0.0000000   -1.0552415    0.0000000
            0.0000000    0.0000000   -0.7420857
    4 I    -1.0552415    0.0000000    0.0000000
            0.0000000   -1.0552415    0.0000000
            0.0000000    0.0000000   -0.7420857
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 65, 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.00000010 (xzx) 0.00000010 (xzx) 0.00000010 (xxz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 18:48:14]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 18:48:14]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P-3m1 (164)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.050840731741600    0.000000000000000    0.000000000000000
  b   -2.025420365870800    3.508130980372969    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.128222920000000
Atomic positions (fractional):
    1 Cu  0.66666666666667  0.33333333333333  0.59368810756266  63.546
    2 Cu  0.33333333333333  0.66666666666667  0.40631189243734  63.546
    3 I   0.33333333333333  0.66666666666667  0.76578333646945 126.904
    4 I   0.66666666666667  0.33333333333333  0.23421666353055 126.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.203362926966399    0.000000000000000    0.000000000000000
  b   -8.101681463483200   14.032523921491876    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.256445840000000
Atomic positions (fractional):
    1 Cu  0.16666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    2 Cu  0.41666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    3 Cu  0.66666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    4 Cu  0.91666666666667  0.08333333333333  0.29684405378133  63.546 > 1
    5 Cu  0.16666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    6 Cu  0.41666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    7 Cu  0.66666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    8 Cu  0.91666666666667  0.33333333333333  0.29684405378133  63.546 > 1
    9 Cu  0.16666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   10 Cu  0.41666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   11 Cu  0.66666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   12 Cu  0.91666666666667  0.58333333333333  0.29684405378133  63.546 > 1
   13 Cu  0.16666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   14 Cu  0.41666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   15 Cu  0.66666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   16 Cu  0.91666666666667  0.83333333333333  0.29684405378133  63.546 > 1
   17 Cu  0.16666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   18 Cu  0.41666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   19 Cu  0.66666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   20 Cu  0.91666666666667  0.08333333333333  0.79684405378133  63.546 > 1
   21 Cu  0.16666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   22 Cu  0.41666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   23 Cu  0.66666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   24 Cu  0.91666666666667  0.33333333333333  0.79684405378133  63.546 > 1
   25 Cu  0.16666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   26 Cu  0.41666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   27 Cu  0.66666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   28 Cu  0.91666666666667  0.58333333333333  0.79684405378133  63.546 > 1
   29 Cu  0.16666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   30 Cu  0.41666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   31 Cu  0.66666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   32 Cu  0.91666666666667  0.83333333333333  0.79684405378133  63.546 > 1
   33 Cu  0.08333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   34 Cu  0.33333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   35 Cu  0.58333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   36 Cu  0.83333333333333  0.16666666666667  0.20315594621867  63.546 > 33
   37 Cu  0.08333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   38 Cu  0.33333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   39 Cu  0.58333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   40 Cu  0.83333333333333  0.41666666666667  0.20315594621867  63.546 > 33
   41 Cu  0.08333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   42 Cu  0.33333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   43 Cu  0.58333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   44 Cu  0.83333333333333  0.66666666666667  0.20315594621867  63.546 > 33
   45 Cu  0.08333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   46 Cu  0.33333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   47 Cu  0.58333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   48 Cu  0.83333333333333  0.91666666666667  0.20315594621867  63.546 > 33
   49 Cu  0.08333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   50 Cu  0.33333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   51 Cu  0.58333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   52 Cu  0.83333333333333  0.16666666666667  0.70315594621867  63.546 > 33
   53 Cu  0.08333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   54 Cu  0.33333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   55 Cu  0.58333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   56 Cu  0.83333333333333  0.41666666666667  0.70315594621867  63.546 > 33
   57 Cu  0.08333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   58 Cu  0.33333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   59 Cu  0.58333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   60 Cu  0.83333333333333  0.66666666666667  0.70315594621867  63.546 > 33
   61 Cu  0.08333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   62 Cu  0.33333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   63 Cu  0.58333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   64 Cu  0.83333333333333  0.91666666666667  0.70315594621867  63.546 > 33
   65 I   0.08333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   66 I   0.33333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   67 I   0.58333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   68 I   0.83333333333333  0.16666666666667  0.38289166823473 126.904 > 65
   69 I   0.08333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   70 I   0.33333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   71 I   0.58333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   72 I   0.83333333333333  0.41666666666667  0.38289166823473 126.904 > 65
   73 I   0.08333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   74 I   0.33333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   75 I   0.58333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   76 I   0.83333333333333  0.66666666666667  0.38289166823473 126.904 > 65
   77 I   0.08333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   78 I   0.33333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   79 I   0.58333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   80 I   0.83333333333333  0.91666666666667  0.38289166823473 126.904 > 65
   81 I   0.08333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   82 I   0.33333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   83 I   0.58333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   84 I   0.83333333333333  0.16666666666667  0.88289166823473 126.904 > 65
   85 I   0.08333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   86 I   0.33333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   87 I   0.58333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   88 I   0.83333333333333  0.41666666666667  0.88289166823473 126.904 > 65
   89 I   0.08333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   90 I   0.33333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   91 I   0.58333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   92 I   0.83333333333333  0.66666666666667  0.88289166823473 126.904 > 65
   93 I   0.08333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   94 I   0.33333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   95 I   0.58333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   96 I   0.83333333333333  0.91666666666667  0.88289166823473 126.904 > 65
   97 I   0.16666666666667  0.08333333333333  0.11710833176527 126.904 > 97
   98 I   0.41666666666667  0.08333333333333  0.11710833176527 126.904 > 97
   99 I   0.66666666666667  0.08333333333333  0.11710833176527 126.904 > 97
  100 I   0.91666666666667  0.08333333333333  0.11710833176527 126.904 > 97
  101 I   0.16666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  102 I   0.41666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  103 I   0.66666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  104 I   0.91666666666667  0.33333333333333  0.11710833176527 126.904 > 97
  105 I   0.16666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  106 I   0.41666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  107 I   0.66666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  108 I   0.91666666666667  0.58333333333333  0.11710833176527 126.904 > 97
  109 I   0.16666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  110 I   0.41666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  111 I   0.66666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  112 I   0.91666666666667  0.83333333333333  0.11710833176527 126.904 > 97
  113 I   0.16666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  114 I   0.41666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  115 I   0.66666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  116 I   0.91666666666667  0.08333333333333  0.61710833176527 126.904 > 97
  117 I   0.16666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  118 I   0.41666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  119 I   0.66666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  120 I   0.91666666666667  0.33333333333333  0.61710833176527 126.904 > 97
  121 I   0.16666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  122 I   0.41666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  123 I   0.66666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  124 I   0.91666666666667  0.58333333333333  0.61710833176527 126.904 > 97
  125 I   0.16666666666667  0.83333333333333  0.61710833176527 126.904 > 97
  126 I   0.41666666666667  0.83333333333333  0.61710833176527 126.904 > 97
  127 I   0.66666666666667  0.83333333333333  0.61710833176527 126.904 > 97
  128 I   0.91666666666667  0.83333333333333  0.61710833176527 126.904 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            6.7156849    0.0000000    0.0000000
            0.0000000    6.7156849    0.0000000
            0.0000000    0.0000000    5.9496968
-------------------------- Born effective charges --------------------------
    1 Cu    1.0552415    0.0000000    0.0000000
            0.0000000    1.0552415    0.0000000
            0.0000000    0.0000000    0.7420857
    2 Cu    1.0552415    0.0000000    0.0000000
            0.0000000    1.0552415    0.0000000
            0.0000000    0.0000000    0.7420857
    3 I    -1.0552415    0.0000000    0.0000000
            0.0000000   -1.0552415    0.0000000
            0.0000000    0.0000000   -0.7420857
    4 I    -1.0552415    0.0000000    0.0000000
            0.0000000   -1.0552415    0.0000000
            0.0000000    0.0000000   -0.7420857
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000010 (xzx) 0.00000010 (xzx) 0.00000010 (xxz)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 14 14 7 ]
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.89, Number of G-points: 313, Lambda: 0.24
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/144) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 144
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   1.104   (   0.000    0.000    0.000)    0.000
   1.104   (   0.000    0.000    0.000)    0.000
   2.793   (   0.000    0.000    0.000)    0.000
   2.793   (   0.000    0.000    0.000)    0.000
   3.359   (   0.000    0.000    0.000)    0.000
   3.835   (   0.000    0.000    0.000)    0.000
   3.835   (   0.000    0.000    0.000)    0.000
   4.796   (   0.000    0.000    0.000)    0.000
   5.161   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/144) =======================
q-point: ( 0.07  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.144   (   9.412    5.434    0.887)   10.904
   0.314   (  12.887    7.440    0.515)   14.889
   0.652   (  26.726   15.430   -0.513)   30.865
   1.140   (   2.975    1.718    0.243)    3.444
   1.298   (  14.297    8.254    0.741)   16.526
   2.813   (   1.656    0.956    0.123)    1.916
   3.315   (  -2.219   -1.281   -0.140)    2.566
   3.369   (   0.276    0.159   -0.089)    0.331
   3.843   (   0.668    0.386   -0.204)    0.798
   3.849   (   1.174    0.678    0.208)    1.372
   4.814   (   1.424    0.822    0.001)    1.645
   5.158   (  -0.242   -0.140   -0.155)    0.319
======================= Grid point 2 (3/144) =======================
q-point: ( 0.14  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.396   (  11.373    6.566    0.435)   13.140
   0.596   (  10.953    6.324    0.217)   12.650
   1.233   (   4.638    2.678    0.365)    5.368
   1.233   (  22.251   12.847   -0.530)   25.699
   1.673   (  16.090    9.290    0.779)   18.596
   2.866   (   2.760    1.593    0.111)    3.189
   3.252   (  -3.033   -1.751   -0.263)    3.512
   3.378   (   0.710    0.410   -0.144)    0.832
   3.865   (   1.247    0.720   -0.232)    1.459
   3.889   (   2.182    1.260    0.432)    2.557
   4.860   (   2.449    1.414    0.004)    2.827
   5.152   (  -0.248   -0.143   -0.229)    0.366
======================= Grid point 3 (4/144) =======================
q-point: ( 0.21  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.665   (  11.336    6.545    0.237)   13.092
   0.823   (   8.304    4.794   -0.117)    9.589
   1.344   (   4.529    2.615    0.307)    5.238
   1.677   (  15.168    8.757   -0.571)   17.524
   2.010   (  12.134    7.005    0.485)   14.019
   2.937   (   3.160    1.825    0.003)    3.649
   3.179   (  -3.013   -1.739   -0.343)    3.496
   3.414   (   2.548    1.471   -0.038)    2.942
   3.900   (   1.629    0.940   -0.083)    1.883
   3.949   (   2.898    1.673    0.527)    3.387
   4.924   (   2.899    1.674    0.003)    3.347
   5.148   (  -0.079   -0.046   -0.204)    0.224
======================= Grid point 4 (5/144) =======================
q-point: ( 0.29  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.923   (  10.509    6.067    0.170)   12.136
   0.986   (   5.603    3.235   -0.333)    6.478
   1.436   (   3.132    1.808    0.136)    3.619
   1.940   (   7.238    4.179   -0.511)    8.374
   2.238   (   7.433    4.292    0.235)    8.587
   3.010   (   3.003    1.734   -0.077)    3.469
   3.115   (  -2.355   -1.360   -0.314)    2.738
   3.501   (   4.709    2.718    0.060)    5.437
   3.939   (   1.680    0.970    0.093)    1.942
   4.018   (   2.807    1.620    0.438)    3.270
   4.992   (   2.730    1.576    0.011)    3.152
   5.148   (   0.124    0.072   -0.144)    0.203
======================= Grid point 5 (6/144) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.090   (   3.345    1.931   -0.359)    3.879
   1.152   (   8.811    5.087    0.141)   10.175
   1.489   (   1.430    0.826   -0.015)    1.651
   2.033   (   1.169    0.675   -0.160)    1.359
   2.375   (   4.554    2.629   -0.047)    5.259
   3.066   (  -1.895   -1.094   -0.138)    2.192
   3.074   (   2.391    1.381   -0.080)    2.762
   3.623   (   5.391    3.112   -0.006)    6.224
   3.976   (   1.355    0.782    0.166)    1.573
   4.072   (   1.619    0.935    0.268)    1.889
   5.048   (   1.984    1.145    0.045)    2.291
   5.154   (   0.329    0.190   -0.115)    0.397
======================= Grid point 6 (7/144) =======================
q-point: ( 0.43  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.148   (   1.579    0.911   -0.222)    1.836
   1.325   (   5.599    3.233    0.081)    6.466
   1.508   (   0.316    0.183   -0.055)    0.370
   2.023   (  -1.323   -0.764    0.080)    1.530
   2.461   (   2.710    1.565   -0.175)    3.135
   3.026   (  -1.421   -0.821    0.016)    1.641
   3.119   (   1.346    0.777   -0.042)    1.555
   3.736   (   3.771    2.177   -0.056)    4.355
   4.001   (   0.747    0.431    0.117)    0.871
   4.093   (   0.248    0.143    0.118)    0.310
   5.083   (   0.959    0.554    0.057)    1.109
   5.162   (   0.347    0.201   -0.088)    0.410
======================= Grid point 7 (8/144) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 200
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.166   (   0.000    0.000    0.000)    0.000
   1.395   (   0.000    0.000   -0.000)    0.000
   1.510   (   0.000    0.000    0.000)    0.000
   2.002   (   0.000    0.000   -0.000)    0.000
   2.495   (   0.000    0.000    0.000)    0.000
   3.007   (   0.000    0.000   -0.000)    0.000
   3.135   (   0.000    0.000   -0.000)    0.000
   3.783   (   0.000    0.000   -0.000)    0.000
   4.010   (   0.000    0.000   -0.000)    0.000
   4.093   (   0.000    0.000    0.000)    0.000
   5.094   (   0.000    0.000    0.000)    0.000
   5.167   (   0.000    0.000    0.000)    0.000
======================= Grid point 16 (9/144) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.321   (   6.330   10.964    0.000)   12.660
   0.563   (   8.072   13.982    0.000)   16.145
   1.054   (  11.616   20.119    0.000)   23.232
   1.257   (   5.004    8.667    0.000)   10.008
   1.535   (   8.505   14.731    0.000)   17.010
   2.860   (   2.013    3.487    0.000)    4.027
   3.267   (  -1.822   -3.156    0.000)    3.645
   3.371   (   0.092    0.160    0.000)    0.185
   3.861   (   0.822    1.424    0.000)    1.644
   3.872   (   0.939    1.626    0.000)    1.877
   4.845   (   1.277    2.213    0.000)    2.555
   5.153   (  -0.165   -0.285    0.000)    0.329
======================= Grid point 17 (10/144) =======================
q-point: ( 0.14  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.561   (   8.769    9.282    0.194)   12.771
   0.828   (   5.015   14.441   -0.039)   15.287
   1.329   (   4.954    7.200    0.229)    8.743
   1.541   (  11.338   14.665   -0.378)   18.540
   1.850   (  12.633    8.233    0.357)   15.084
   2.942   (   1.710    5.543    0.056)    5.801
   3.195   (  -2.156   -3.657   -0.202)    4.250
   3.385   (   1.446    0.478   -0.046)    1.524
   3.896   (   0.788    2.253   -0.223)    2.397
   3.917   (   2.418    1.528    0.421)    2.891
   4.900   (   2.123    2.380    0.002)    3.189
   5.148   (  -0.109   -0.190   -0.119)    0.249
======================= Grid point 18 (11/144) =======================
q-point: ( 0.21  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.810   (   9.417    7.691    0.240)   12.160
   1.027   (   1.400   13.692   -0.327)   13.767
   1.459   (   1.401    8.044    0.298)    8.171
   1.839   (   9.177    6.651   -0.562)   11.348
   2.120   (  10.420    4.178    0.292)   11.230
   3.018   (   0.158    5.078   -0.016)    5.081
   3.130   (  -1.155   -2.472   -0.284)    2.743
   3.441   (   3.985    1.303    0.036)    4.193
   3.933   (   0.851    2.084   -0.026)    2.251
   3.985   (   3.020    1.962    0.416)    3.625
   4.964   (   2.435    2.172    0.002)    3.263
   5.146   (   0.067   -0.099   -0.123)    0.171
======================= Grid point 19 (12/144) =======================
q-point: ( 0.29  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.042   (   8.589    6.237    0.337)   10.620
   1.160   (  -0.977   12.062   -0.592)   12.116
   1.544   (  -1.190    8.218    0.219)    8.306
   1.992   (   4.037    1.023   -0.395)    4.184
   2.303   (   6.767    2.586    0.121)    7.245
   3.029   (  -0.231   -2.059   -0.177)    2.080
   3.127   (  -0.622    4.276   -0.108)    4.322
   3.548   (   5.880    1.956    0.026)    6.197
   3.965   (   0.822    1.100    0.158)    1.383
   4.055   (   2.283    2.283    0.265)    3.240
   5.025   (   2.143    1.679    0.022)    2.723
   5.149   (   0.291   -0.035   -0.102)    0.310
======================= Grid point 20 (13/144) =======================
q-point: ( 0.36  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.210   (   0.604    8.057    0.310)    8.086
   1.264   (   3.302    6.863   -0.585)    7.638
   1.588   (  -3.177    8.226    0.147)    8.819
   2.019   (   0.197   -1.908   -0.056)    1.919
   2.422   (   4.070    2.230   -0.104)    4.642
   2.997   (   0.367   -3.835   -0.064)    3.853
   3.163   (  -1.242    5.242   -0.040)    5.388
   3.673   (   5.630    1.848   -0.073)    5.926
   3.980   (   0.677   -0.712    0.210)    1.005
   4.106   (   0.545    2.574    0.114)    2.633
   5.070   (   1.329    1.033    0.053)    1.684
   5.156   (   0.480   -0.023   -0.095)    0.490
======================= Grid point 21 (14/144) =======================
q-point: ( 0.43  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.257   (  -4.081    8.954   -0.019)    9.840
   1.374   (   2.482    2.049   -0.119)    3.221
   1.606   (  -4.302    8.362    0.074)    9.404
   1.990   (  -0.087   -2.335    0.034)    2.337
   2.490   (   1.348    1.362   -0.090)    1.918
   2.970   (   1.125   -3.574    0.013)    3.747
   3.185   (  -1.978    4.583   -0.011)    4.992
   3.765   (   2.659    0.799   -0.067)    2.777
   3.976   (   1.007   -2.516    0.110)    2.712
   4.128   (  -1.073    2.848    0.026)    3.044
   5.093   (   0.369    0.419    0.033)    0.559
   5.163   (   0.358   -0.172   -0.046)    0.400
======================= Grid point 31 (15/144) =======================
q-point: ( 0.14  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.764   (   6.009   10.408    0.000)   12.018
   1.114   (   7.231   12.524    0.000)   14.462
   1.488   (   4.493    7.781    0.000)    8.985
   1.798   (   5.847   10.127    0.000)   11.693
   2.021   (   5.077    8.793    0.000)   10.153
   3.067   (   3.466    6.004    0.000)    6.933
   3.112   (  -2.424   -4.199    0.000)    4.849
   3.407   (   1.128    1.953    0.000)    2.256
   3.944   (   0.903    1.564    0.000)    1.806
   3.960   (   1.980    3.430    0.000)    3.961
   4.953   (   1.611    2.790    0.000)    3.222
   5.145   (  -0.062   -0.107    0.000)    0.123
======================= Grid point 32 (16/144) =======================
q-point: ( 0.21  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.977   (   6.998    8.676    0.120)   11.147
   1.316   (   0.728   13.441   -0.290)   13.464
   1.631   (   2.397    7.960    0.323)    8.319
   1.946   (   2.952    3.882   -0.393)    4.893
   2.203   (   7.360    4.506    0.152)    8.631
   3.026   (  -2.267   -3.604   -0.101)    4.259
   3.177   (   1.455    6.145   -0.044)    6.315
   3.473   (   4.078    1.929    0.023)    4.512
   3.967   (   0.297    0.889    0.065)    0.939
   4.040   (   2.658    3.572    0.116)    4.454
   5.010   (   1.790    2.249    0.005)    2.874
   5.144   (   0.101   -0.082   -0.042)    0.137
======================= Grid point 33 (17/144) =======================
q-point: ( 0.29  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.179   (   6.524    7.194    0.201)    9.713
   1.422   (  -3.363   12.605   -0.534)   13.057
   1.729   (  -0.957    8.765    0.539)    8.834
   1.990   (   0.566   -0.881   -0.387)    1.116
   2.361   (   5.850    3.354    0.056)    6.743
   2.960   (  -0.830   -3.929   -0.100)    4.017
   3.244   (  -0.924    6.399   -0.033)    6.465
   3.584   (   6.045    1.677   -0.054)    6.274
   3.969   (  -0.170   -0.743    0.189)    0.785
   4.117   (   1.733    3.641    0.060)    4.033
   5.059   (   1.440    1.567    0.027)    2.128
   5.147   (   0.329   -0.113   -0.060)    0.353
======================= Grid point 34 (18/144) =======================
q-point: ( 0.36  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.343   (   4.079    5.809    0.283)    7.103
   1.464   (  -5.244   10.826   -0.586)   12.043
   1.782   (  -3.516    9.662    0.538)   10.296
   1.963   (  -0.600   -3.289   -0.278)    3.355
   2.473   (   3.066    2.921   -0.046)    4.235
   2.915   (   0.411   -4.115   -0.018)    4.136
   3.279   (  -2.273    6.105   -0.011)    6.514
   3.704   (   5.325    1.247   -0.142)    5.471
   3.945   (  -0.496   -2.618    0.211)    2.673
   4.170   (   0.081    3.399    0.012)    3.400
   5.090   (   0.660    0.901    0.043)    1.118
   5.153   (   0.459   -0.247   -0.060)    0.525
======================= Grid point 35 (19/144) =======================
q-point: ( 0.43  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.419   (  -1.416    2.453    0.000)    2.833
   1.466   (  -5.961   10.324   -0.000)   11.921
   1.802   (  -5.753    9.964   -0.000)   11.506
   1.932   (   1.838   -3.184    0.000)    3.676
   2.517   (  -0.775    1.342    0.000)    1.549
   2.896   (   2.056   -3.561    0.000)    4.112
   3.290   (  -3.246    5.621    0.000)    6.491
   3.767   (   0.438   -0.759    0.000)    0.876
   3.918   (   1.603   -2.777   -0.000)    3.207
   4.188   (  -1.592    2.758   -0.000)    3.185
   5.100   (  -0.204    0.354   -0.000)    0.409
   5.156   (   0.265   -0.459   -0.000)    0.530
======================= Grid point 47 (20/144) =======================
q-point: ( 0.21  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.158   (   5.128    8.881    0.000)   10.255
   1.540   (   4.637    8.032    0.000)    9.275
   1.765   (   2.928    5.072    0.000)    5.856
   1.995   (   0.537    0.930    0.000)    1.074
   2.315   (   3.674    6.364    0.000)    7.349
   2.948   (  -2.144   -3.714    0.000)    4.288
   3.296   (   2.920    5.057    0.000)    5.839
   3.527   (   2.087    3.615    0.000)    4.174
   3.968   (  -0.437   -0.756    0.000)    0.873
   4.120   (   2.299    3.983    0.000)    4.599
   5.054   (   1.146    1.985    0.000)    2.292
   5.143   (  -0.000   -0.000    0.000)    0.000
======================= Grid point 48 (21/144) =======================
q-point: ( 0.29  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.331   (   4.772    7.619    0.068)    8.991
   1.641   (  -2.269    8.456   -0.405)    8.764
   1.839   (   2.240    0.568    0.705)    2.416
   1.993   (  -2.142    2.378   -0.405)    3.226
   2.446   (   4.140    5.063    0.021)    6.541
   2.879   (  -1.149   -3.926   -0.023)    4.091
   3.374   (  -0.523    6.079    0.004)    6.102
   3.621   (   4.959    2.121   -0.080)    5.394
   3.936   (  -1.255   -2.456    0.105)    2.760
   4.192   (   1.360    3.442    0.002)    3.701
   5.088   (   0.813    1.317    0.020)    1.548
   5.144   (   0.186   -0.168   -0.028)    0.252
======================= Grid point 49 (22/144) =======================
q-point: ( 0.36  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.464   (   1.450    6.072    0.076)    6.244
   1.651   (  -4.704    6.826   -0.307)    8.295
   1.858   (   2.129   -3.603    0.532)    4.219
   1.983   (  -3.941    6.408   -0.309)    7.529
   2.544   (   0.988    4.027   -0.002)    4.146
   2.831   (   0.652   -4.174   -0.002)    4.225
   3.407   (  -2.733    6.217    0.008)    6.792
   3.726   (   4.202    0.928   -0.119)    4.305
   3.879   (  -1.051   -3.717    0.123)    3.865
   4.233   (  -0.300    2.608   -0.002)    2.626
   5.107   (   0.093    0.788    0.020)    0.794
   5.146   (   0.305   -0.417   -0.024)    0.517
======================= Grid point 62 (23/144) =======================
q-point: ( 0.29  0.29  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.483   (   4.086    7.076    0.000)    8.171
   1.761   (   1.718    2.976    0.000)    3.437
   1.786   (  -2.214   -3.835    0.000)    4.428
   2.054   (   1.397    2.419    0.000)    2.793
   2.558   (   3.127    5.417    0.000)    6.255
   2.801   (  -2.174   -3.765    0.000)    4.347
   3.477   (   2.118    3.669    0.000)    4.236
   3.677   (   2.033    3.521    0.000)    4.066
   3.871   (  -2.194   -3.800    0.000)    4.387
   4.249   (   1.168    2.024    0.000)    2.337
   5.111   (   0.516    0.894    0.000)    1.033
   5.140   (  -0.138   -0.239    0.000)    0.276
======================= Grid point 63 (24/144) =======================
q-point: ( 0.36  0.29  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.583   (  -2.979    5.160    0.000)    5.959
   1.710   (   1.055   -1.828   -0.000)    2.110
   1.799   (   0.313   -0.543    0.000)    0.627
   2.079   (  -1.414    2.448    0.000)    2.827
   2.632   (  -2.406    4.168   -0.000)    4.812
   2.746   (   2.276   -3.942    0.000)    4.552
   3.515   (  -2.065    3.577   -0.000)    4.131
   3.741   (  -0.348    0.603    0.000)    0.696
   3.806   (   1.638   -2.837   -0.000)    3.276
   4.270   (  -0.557    0.964   -0.000)    1.113
   5.122   (  -0.350    0.606   -0.000)    0.700
   5.136   (   0.279   -0.484    0.000)    0.559
======================= Grid point 211 (25/144) =======================
q-point: ( 0.00  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 144
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.178   (   0.000    0.000    8.421)    8.421
   0.178   (   0.000   -0.000    8.421)    8.421
   0.405   (  -0.000    0.000   19.004)   19.004
   1.094   (   0.000    0.000   -0.923)    0.923
   1.094   (   0.000   -0.000   -0.923)    0.923
   2.790   (   0.000    0.000   -0.279)    0.279
   2.790   (   0.000   -0.000   -0.279)    0.279
   3.350   (   0.000    0.000   -0.902)    0.902
   3.836   (   0.000    0.000    0.088)    0.088
   3.836   (   0.000   -0.000    0.088)    0.088
   4.785   (  -0.000    0.000   -1.072)    1.072
   5.349   (   0.000   -0.000    0.590)    0.590
======================= Grid point 212 (26/144) =======================
q-point: ( 0.07  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.368   (  11.137    6.430    4.314)   13.564
   0.383   (   3.284    1.896   14.598)   15.083
   0.683   (  24.282   14.019    3.586)   28.267
   1.135   (   3.171    1.831   -0.705)    3.729
   1.301   (  14.926    8.618   -0.575)   17.245
   2.812   (   1.715    0.990   -0.194)    1.990
   3.070   (  10.308    5.951  -13.391)   17.917
   3.344   (  -0.300   -0.173   -1.125)    1.177
   3.840   (   0.556    0.321   -0.052)    0.644
   3.854   (   1.333    0.769    0.214)    1.554
   4.797   (   1.211    0.699   -1.197)    1.841
   5.259   (  -4.040   -2.333    5.182)    6.973
======================= Grid point 213 (27/144) =======================
q-point: ( 0.14  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.511   (   7.525    4.345    9.357)   12.769
   0.625   (  10.238    5.911    2.455)   12.074
   1.232   (   4.699    2.713   -0.520)    5.451
   1.235   (  21.441   12.379    0.596)   24.765
   1.682   (  16.046    9.264   -0.046)   18.528
   2.866   (   2.739    1.582   -0.149)    3.167
   3.177   (   0.504    0.291   -6.085)    6.113
   3.351   (   1.165    0.672   -1.788)    2.237
   3.862   (   1.277    0.737   -0.114)    1.479
   3.897   (   2.252    1.300    0.281)    2.616
   4.842   (   2.503    1.445   -1.437)    3.227
   5.198   (  -1.505   -0.869    3.822)    4.199
======================= Grid point 214 (28/144) =======================
q-point: ( 0.21  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.714   (   9.394    5.424    4.386)   11.701
   0.839   (   7.886    4.553    1.598)    9.245
   1.342   (   4.424    2.554   -0.514)    5.135
   1.662   (  14.500    8.372   -0.966)   16.771
   2.017   (  12.082    6.975    0.125)   13.952
   2.936   (   3.144    1.815   -0.097)    3.632
   3.151   (  -1.964   -1.134   -2.084)    3.080
   3.404   (   3.342    1.930   -0.860)    3.954
   3.898   (   1.732    1.000   -0.096)    2.002
   3.957   (   2.847    1.644    0.236)    3.296
   4.909   (   3.043    1.757   -1.313)    3.751
   5.176   (  -0.549   -0.317    2.598)    2.675
======================= Grid point 215 (29/144) =======================
q-point: ( 0.29  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.939   (   9.464    5.464    1.376)   11.014
   0.995   (   5.416    3.127    1.165)    6.361
   1.430   (   2.964    1.712   -0.666)    3.487
   1.912   (   6.973    4.026   -2.085)    8.317
   2.246   (   7.482    4.320    0.446)    8.651
   3.009   (   3.038    1.754    0.003)    3.508
   3.103   (  -1.948   -1.125   -0.790)    2.384
   3.505   (   5.123    2.958    0.387)    5.928
   3.940   (   1.745    1.008   -0.072)    2.017
   4.024   (   2.680    1.548    0.116)    3.097
   4.980   (   2.877    1.661   -1.053)    3.484
   5.168   (  -0.128   -0.074    1.913)    1.919
======================= Grid point 216 (30/144) =======================
q-point: ( 0.36  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.097   (   3.340    1.928    1.000)    3.984
   1.149   (   8.221    4.746   -0.355)    9.499
   1.480   (   1.324    0.765   -0.837)    1.743
   2.005   (   1.354    0.782   -2.397)    2.862
   2.381   (   4.415    2.549    0.620)    5.136
   3.061   (  -1.624   -0.938   -0.300)    1.899
   3.075   (   2.448    1.414    0.115)    2.830
   3.633   (   5.419    3.129    0.884)    6.319
   3.977   (   1.340    0.774   -0.088)    1.550
   4.075   (   1.525    0.880    0.023)    1.761
   5.039   (   2.110    1.218   -0.810)    2.568
   5.169   (   0.133    0.077    1.441)    1.449
======================= Grid point 217 (31/144) =======================
q-point: ( 0.43  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.156   (   1.689    0.975    0.993)    2.189
   1.312   (   5.280    3.048   -1.295)    6.232
   1.498   (   0.319    0.184   -0.912)    0.984
   2.000   (  -1.239   -0.716   -2.264)    2.678
   2.464   (   2.658    1.534    0.481)    3.107
   3.026   (  -1.277   -0.738    0.013)    1.475
   3.121   (   1.391    0.803    0.192)    1.618
   3.744   (   3.735    2.156    0.876)    4.401
   4.001   (   0.671    0.388   -0.145)    0.789
   4.095   (   0.207    0.120    0.030)    0.241
   5.076   (   1.012    0.585   -0.629)    1.327
   5.174   (   0.253    0.146    1.098)    1.136
======================= Grid point 218 (32/144) =======================
q-point: (-0.50  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.178   (   0.178    0.103    1.099)    1.118
   1.378   (  -0.058   -0.034   -1.638)    1.640
   1.501   (   0.052    0.030   -0.896)    0.898
   1.978   (  -0.119   -0.068   -2.255)    2.259
   2.500   (   0.183    0.106    0.470)    0.515
   3.009   (  -0.046   -0.026    0.124)    0.135
   3.138   (   0.031    0.018    0.232)    0.235
   3.792   (   0.055    0.032    0.831)    0.833
   4.007   (  -0.096   -0.055   -0.225)    0.250
   4.094   (  -0.081   -0.047    0.035)    0.100
   5.087   (  -0.043   -0.025   -0.568)    0.570
   5.178   (   0.065    0.037    0.983)    0.986
======================= Grid point 220 (33/144) =======================
q-point: (-0.43  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.163   (  -1.396   -0.806    1.288)    2.063
   1.309   (  -5.391   -3.112   -1.386)    6.377
   1.499   (  -0.297   -0.171   -0.864)    0.929
   1.997   (   1.314    0.758   -2.388)    2.829
   2.470   (  -2.590   -1.496    0.702)    3.072
   3.026   (   1.389    0.802    0.016)    1.604
   3.122   (  -1.327   -0.766    0.247)    1.552
   3.746   (  -3.729   -2.153    0.913)    4.402
   3.996   (  -0.812   -0.469   -0.293)    0.982
   4.091   (  -0.396   -0.229   -0.103)    0.469
   5.075   (  -1.037   -0.598   -0.646)    1.360
   5.176   (  -0.180   -0.104    1.153)    1.172
======================= Grid point 221 (34/144) =======================
q-point: (-0.36  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.110   (  -3.237   -1.869    1.481)    4.020
   1.144   (  -8.268   -4.773   -0.496)    9.560
   1.480   (  -1.482   -0.856   -0.858)    1.914
   2.010   (  -0.802   -0.463   -2.276)    2.457
   2.383   (  -4.786   -2.763    0.722)    5.574
   3.067   (   1.924    1.111   -0.098)    2.223
   3.078   (  -2.410   -1.391    0.217)    2.791
   3.633   (  -5.547   -3.202    0.866)    6.463
   3.971   (  -1.330   -0.768   -0.297)    1.564
   4.066   (  -1.780   -1.028   -0.300)    2.078
   5.038   (  -2.064   -1.191   -0.813)    2.518
   5.172   (  -0.107   -0.062    1.529)    1.534
======================= Grid point 222 (35/144) =======================
q-point: (-0.29  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.933   (  -9.496   -5.483    1.235)   11.034
   1.007   (  -5.598   -3.232    1.626)    6.665
   1.425   (  -3.238   -1.869   -0.887)    3.843
   1.929   (  -6.645   -3.836   -1.564)    7.831
   2.239   (  -7.815   -4.512    0.295)    9.029
   3.012   (  -3.097   -1.788    0.100)    3.578
   3.114   (   2.090    1.207   -0.421)    2.450
   3.503   (  -5.111   -2.951    0.322)    5.911
   3.936   (  -1.541   -0.890   -0.188)    1.790
   4.009   (  -2.909   -1.679   -0.432)    3.386
   4.980   (  -2.845   -1.643   -1.029)    3.443
   5.173   (   0.197    0.114    2.042)    2.055
======================= Grid point 223 (36/144) =======================
q-point: (-0.21  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.707   (  -9.528   -5.501    4.225)   11.785
   0.844   (  -8.322   -4.805    1.795)    9.776
   1.331   (  -4.622   -2.668   -0.947)    5.420
   1.681   ( -14.553   -8.402   -0.337)   16.807
   2.002   ( -12.438   -7.181   -0.189)   14.363
   2.936   (  -3.309   -1.911   -0.099)    3.823
   3.163   (   1.938    1.119   -1.698)    2.809
   3.405   (  -3.161   -1.825   -0.847)    3.747
   3.901   (  -1.448   -0.836    0.008)    1.672
   3.939   (  -2.848   -1.645   -0.416)    3.315
   4.909   (  -3.040   -1.755   -1.274)    3.734
   5.182   (   0.609    0.352    2.764)    2.852
======================= Grid point 224 (37/144) =======================
q-point: (-0.14  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.499   (  -7.719   -4.456    9.193)   12.804
   0.618   ( -10.724   -6.192    2.266)   12.589
   1.219   (  -4.646   -2.682   -1.010)    5.459
   1.253   ( -21.496  -12.411    1.157)   24.848
   1.660   ( -16.154   -9.327   -0.340)   18.657
   2.862   (  -2.858   -1.650   -0.284)    3.312
   3.188   (  -0.517   -0.299   -5.620)    5.651
   3.353   (  -1.255   -0.724   -1.828)    2.332
   3.870   (  -1.163   -0.671    0.174)    1.354
   3.882   (  -1.988   -1.148   -0.250)    2.309
   4.842   (  -2.507   -1.447   -1.403)    3.217
   5.205   (   1.427    0.824    3.931)    4.262
======================= Grid point 225 (38/144) =======================
q-point: (-0.07  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.353   ( -11.282   -6.514    3.952)   13.614
   0.371   (  -2.998   -1.731   14.655)   15.058
   0.700   ( -24.361  -14.065    4.066)   28.422
   1.127   (  -2.868   -1.656   -1.025)    3.467
   1.288   ( -14.019   -8.094   -0.242)   16.189
   2.808   (  -1.619   -0.935   -0.344)    1.901
   3.078   ( -10.613   -6.127  -13.163)   17.984
   3.344   (   0.262    0.151   -1.128)    1.168
   3.846   (  -1.019   -0.588   -0.085)    1.180
   3.847   (  -0.754   -0.435    0.200)    0.893
   4.797   (  -1.220   -0.705   -1.190)    1.844
   5.262   (   3.834    2.214    5.152)    6.793
======================= Grid point 227 (39/144) =======================
q-point: ( 0.07  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.456   (   4.133    6.029   10.872)   13.101
   0.590   (   8.066   13.159    2.551)   15.644
   1.068   (  10.437   19.096    1.224)   21.796
   1.246   (   5.431    9.005   -0.969)   10.560
   1.539   (   9.395   14.291    0.226)   17.104
   2.857   (   2.206    3.386   -0.267)    4.050
   3.169   (   0.624    2.074   -7.905)    8.196
   3.345   (   0.201    0.488   -1.652)    1.735
   3.860   (   0.420    1.593   -0.095)    1.650
   3.875   (   1.469    1.283    0.177)    1.958
   4.827   (   1.270    2.222   -1.396)    2.915
   5.212   (  -1.317   -1.796    4.316)    4.856
======================= Grid point 228 (40/144) =======================
q-point: ( 0.14  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.626   (   7.168    6.949    5.742)   11.517
   0.846   (   4.702   14.009    1.713)   14.876
   1.327   (   5.081    7.025   -0.495)    8.683
   1.532   (  10.780   14.279   -0.522)   17.899
   1.858   (  12.661    8.140    0.290)   15.054
   2.941   (   1.767    5.500   -0.188)    5.780
   3.157   (  -1.048   -2.163   -3.065)    3.895
   3.370   (   1.897    1.180   -1.273)    2.571
   3.891   (   0.892    2.178   -0.163)    2.359
   3.925   (   2.542    1.419    0.224)    2.920
   4.883   (   2.213    2.505   -1.370)    3.613
   5.182   (  -0.722   -0.660    2.984)    3.140
======================= Grid point 229 (41/144) =======================
q-point: ( 0.21  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.835   (   8.477    6.512    2.254)   10.924
   1.036   (   1.126   13.383    1.143)   13.479
   1.458   (   1.314    8.132   -0.468)    8.251
   1.815   (   8.704    6.310   -1.687)   10.882
   2.127   (  10.435    4.263    0.391)   11.279
   3.016   (   0.119    4.922   -0.165)    4.926
   3.115   (  -0.685   -1.732   -1.021)    2.124
   3.441   (   4.480    1.711   -0.072)    4.796
   3.931   (   1.021    2.114   -0.109)    2.350
   3.991   (   2.991    1.750    0.134)    3.468
   4.951   (   2.551    2.306   -1.140)    3.623
   5.169   (  -0.243   -0.311    2.154)    2.189
======================= Grid point 230 (42/144) =======================
q-point: ( 0.29  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.047   (   8.267    5.567    0.222)    9.969
   1.162   (  -1.249   11.820    0.766)   11.911
   1.540   (  -1.471    8.457   -0.557)    8.602
   1.963   (   4.087    0.994   -2.337)    4.812
   2.311   (   6.691    2.603    0.608)    7.205
   3.024   (  -0.133   -2.106   -0.227)    2.122
   3.125   (  -0.478    4.658   -0.121)    4.684
   3.555   (   6.051    2.000    0.698)    6.411
   3.966   (   0.881    1.196   -0.036)    1.486
   4.058   (   2.239    2.051   -0.007)    3.036
   5.015   (   2.256    1.796   -0.885)    3.016
   5.166   (   0.086   -0.196    1.611)    1.625
======================= Grid point 231 (43/144) =======================
q-point: ( 0.36  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.219   (   1.639    7.407    0.382)    7.596
   1.249   (   2.125    6.931   -0.681)    7.281
   1.582   (  -3.448    8.468   -0.716)    9.172
   1.994   (   0.420   -1.852   -2.317)    2.995
   2.427   (   3.897    2.208    0.577)    4.516
   2.996   (   0.516   -3.710   -0.045)    3.746
   3.164   (  -1.178    5.343    0.144)    5.473
   3.681   (   5.629    1.734    0.842)    5.949
   3.983   (   0.608   -0.579    0.011)    0.840
   4.106   (   0.552    2.385   -0.094)    2.450
   5.063   (   1.408    1.124   -0.671)    1.923
   5.168   (   0.341   -0.159    1.198)    1.256
======================= Grid point 232 (44/144) =======================
q-point: ( 0.43  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.266   (  -4.068    8.927    0.861)    9.848
   1.355   (   2.703    1.826   -1.629)    3.646
   1.597   (  -4.499    8.378   -0.897)    9.552
   1.967   (  -0.087   -2.314   -2.173)    3.176
   2.494   (   1.379    1.423    0.441)    2.030
   2.972   (   1.161   -3.533    0.133)    3.721
   3.187   (  -1.968    4.646    0.247)    5.052
   3.773   (   2.745    0.746    0.784)    2.950
   3.978   (   0.800   -2.434   -0.002)    2.562
   4.127   (  -1.038    2.711   -0.115)    2.905
   5.087   (   0.350    0.461   -0.540)    0.792
   5.173   (   0.370   -0.235    0.950)    1.047
======================= Grid point 233 (45/144) =======================
q-point: (-0.50  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.265   (  -5.411    7.196    0.708)    9.031
   1.360   (  -2.913   -0.321   -1.289)    3.202
   1.594   (  -5.040    7.700   -1.003)    9.257
   1.966   (   1.906   -0.988   -2.227)    3.094
   2.497   (  -1.660   -0.428    0.558)    1.803
   2.972   (   2.470   -2.746    0.125)    3.696
   3.188   (  -2.995    4.035    0.263)    5.032
   3.775   (  -1.934   -1.832    0.857)    2.798
   3.974   (   1.491   -2.177   -0.144)    2.642
   4.126   (  -1.946    2.264   -0.141)    2.989
   5.086   (  -0.640   -0.105   -0.552)    0.852
   5.174   (   0.128   -0.390    0.978)    1.061
======================= Grid point 234 (46/144) =======================
q-point: (-0.43  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.203   (  -6.820    0.121   -0.507)    6.840
   1.274   (  -7.494    4.025    0.618)    8.529
   1.576   (  -5.375    6.603   -0.958)    8.567
   1.995   (   1.607   -0.860   -2.290)    2.926
   2.431   (  -3.964   -2.521    0.726)    4.754
   2.998   (   3.161   -2.160    0.058)    3.829
   3.165   (  -3.962    3.690    0.201)    5.418
   3.684   (  -4.471   -3.982    0.908)    6.056
   3.975   (   0.232   -0.945   -0.257)    1.007
   4.102   (  -2.627    0.788   -0.224)    2.752
   5.062   (  -1.662   -0.633   -0.684)    1.906
   5.171   (   0.011   -0.395    1.262)    1.323
======================= Grid point 235 (47/144) =======================
q-point: (-0.36  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.036   (  -8.833   -4.140   -0.125)    9.756
   1.183   (  -9.913    6.785    1.548)   12.112
   1.532   (  -6.442    4.921   -0.921)    8.159
   1.976   (  -2.483   -2.635   -1.948)    4.111
   2.308   (  -5.904   -4.742    0.551)    7.593
   3.031   (   1.865   -0.899    0.014)    2.071
   3.128   (  -3.511    2.666    0.006)    4.409
   3.555   (  -4.903   -4.264    0.662)    6.532
   3.960   (  -1.274   -0.088   -0.237)    1.299
   4.049   (  -3.312   -0.778   -0.334)    3.419
   5.015   (  -2.645   -1.026   -0.877)    2.970
   5.169   (   0.186   -0.225    1.697)    1.722
======================= Grid point 236 (48/144) =======================
q-point: (-0.29  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.827   (  -9.868   -3.920    1.967)   10.799
   1.049   ( -12.548    5.425    1.690)   13.774
   1.447   (  -7.675    2.749   -0.937)    8.206
   1.834   (  -9.782   -4.511   -1.075)   10.826
   2.118   (  -9.140   -6.826    0.184)   11.409
   3.016   (  -4.519    2.164   -0.123)    5.012
   3.124   (   2.053   -0.450   -0.719)    2.222
   3.440   (  -3.623   -2.967   -0.108)    4.684
   3.932   (  -2.050    0.459   -0.077)    2.102
   3.976   (  -3.351   -1.503   -0.385)    3.693
   4.951   (  -3.260   -1.056   -1.116)    3.604
   5.173   (   0.498   -0.064    2.264)    2.319
======================= Grid point 237 (49/144) =======================
q-point: (-0.21  0.07  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.620   (  -9.717   -2.267    5.534)   11.410
   0.848   ( -14.880    2.746    1.889)   15.249
   1.318   (  -8.552   -0.214   -0.868)    8.598
   1.544   ( -17.870   -3.149   -0.103)   18.145
   1.847   ( -13.717   -6.146    0.089)   15.031
   2.939   (  -5.708    1.280   -0.256)    5.855
   3.165   (   2.432   -0.712   -2.810)    3.784
   3.371   (  -1.880   -1.144   -1.263)    2.537
   3.899   (  -2.350    0.448    0.093)    2.394
   3.910   (  -2.376   -0.994   -0.274)    2.590
   4.883   (  -3.275   -0.680   -1.348)    3.606
   5.186   (   1.007   -0.010    3.070)    3.231
======================= Grid point 242 (50/144) =======================
q-point: ( 0.14  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.789   (   5.399    8.830    2.428)   10.630
   1.128   (   6.765   12.543    1.330)   14.313
   1.481   (   4.974    7.485   -0.725)    9.017
   1.785   (   4.887   10.091   -1.279)   11.285
   2.027   (   5.464    8.680    0.495)   10.269
   3.065   (   3.758    5.896   -0.149)    6.994
   3.098   (  -2.494   -3.135   -1.152)    4.168
   3.404   (   1.437    2.428   -0.281)    2.835
   3.942   (   0.607    1.937   -0.162)    2.036
   3.960   (   2.571    2.809   -0.018)    3.808
   4.940   (   1.696    2.951   -1.151)    3.593
   5.171   (  -0.385   -0.434    2.297)    2.369
======================= Grid point 243 (51/144) =======================
q-point: ( 0.21  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.981   (   6.609    7.732    0.317)   10.177
   1.322   (   0.273   13.430    0.803)   13.456
   1.630   (   2.790    7.802   -0.475)    8.300
   1.919   (   2.441    3.943   -2.116)    5.097
   2.212   (   7.418    4.522    0.646)    8.712
   3.021   (  -2.195   -3.303   -0.330)    3.979
   3.176   (   1.453    6.357   -0.000)    6.521
   3.478   (   4.295    2.079    0.453)    4.793
   3.967   (   0.440    0.947   -0.048)    1.045
   4.040   (   2.722    3.320   -0.101)    4.295
   5.000   (   1.880    2.390   -0.902)    3.172
   5.163   (  -0.120   -0.298    1.714)    1.744
======================= Grid point 244 (52/144) =======================
q-point: ( 0.29  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.171   (   6.387    6.574   -0.957)    9.215
   1.421   (  -3.576   12.528    0.435)   13.035
   1.731   (  -1.056    9.008   -0.368)    9.077
   1.960   (   0.677   -0.957   -2.395)    2.667
   2.369   (   5.715    3.357    0.682)    6.663
   2.959   (  -0.770   -3.753   -0.050)    3.832
   3.246   (  -0.895    6.557    0.189)    6.620
   3.591   (   6.051    1.601    0.714)    6.300
   3.972   (  -0.062   -0.703    0.087)    0.711
   4.116   (   1.731    3.491   -0.201)    3.902
   5.051   (   1.526    1.690   -0.672)    2.374
   5.161   (   0.161   -0.301    1.254)    1.299
======================= Grid point 245 (53/144) =======================
q-point: ( 0.36  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.328   (   4.206    5.203   -1.618)    6.884
   1.461   (  -5.241   10.954    0.300)   12.147
   1.781   (  -4.121   10.060   -0.583)   10.887
   1.936   (  -0.150   -3.602   -2.223)    4.236
   2.479   (   2.924    2.964    0.548)    4.200
   2.917   (   0.492   -4.015    0.125)    4.047
   3.282   (  -2.273    6.220    0.323)    6.630
   3.709   (   5.332    1.144    0.671)    5.495
   3.949   (  -0.527   -2.568    0.192)    2.628
   4.167   (   0.100    3.307   -0.256)    3.318
   5.085   (   0.691    1.002   -0.486)    1.310
   5.162   (   0.392   -0.395    0.903)    1.061
======================= Grid point 246 (54/144) =======================
q-point: ( 0.43  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.396   (  -1.887    2.282   -2.026)    3.588
   1.474   (  -4.966   10.278    0.643)   11.433
   1.789   (  -6.167    9.216   -1.204)   11.154
   1.913   (   1.915   -2.519   -1.795)    3.638
   2.523   (  -0.721    1.410    0.486)    1.657
   2.899   (   2.024   -3.525    0.216)    4.071
   3.294   (  -3.280    5.693    0.370)    6.581
   3.775   (   0.686   -0.675    0.734)    1.210
   3.919   (   1.288   -2.866    0.090)    3.143
   4.185   (  -1.558    2.695   -0.265)    3.124
   5.096   (  -0.284    0.412   -0.413)    0.649
   5.165   (   0.373   -0.536    0.773)    1.012
======================= Grid point 248 (55/144) =======================
q-point: (-0.50  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.320   (  -6.864   -0.536   -1.812)    7.119
   1.481   (  -6.414    9.788    0.940)   11.740
   1.761   (  -6.159    7.034   -1.403)    9.454
   1.946   (   2.655   -0.364   -1.823)    3.241
   2.480   (  -3.982   -1.210    0.626)    4.209
   2.917   (   3.330   -2.375    0.156)    4.093
   3.283   (  -4.207    5.063    0.337)    6.591
   3.714   (  -3.722   -4.020    0.844)    5.543
   3.942   (   2.386   -0.833   -0.085)    2.528
   4.167   (  -2.982    1.592   -0.264)    3.391
   5.084   (  -1.234   -0.088   -0.500)    1.334
   5.164   (   0.202   -0.552    0.938)    1.107
======================= Grid point 249 (56/144) =======================
q-point: (-0.43  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.164   (  -8.984   -1.871   -1.177)    9.252
   1.441   (  -9.099    8.953    1.200)   12.822
   1.711   (  -6.708    4.878   -1.234)    8.386
   1.973   (   0.414   -0.590   -1.910)    2.041
   2.367   (  -5.934   -3.445    0.664)    6.893
   2.962   (   3.805   -1.271    0.090)    4.012
   3.247   (  -5.164    4.019    0.231)    6.548
   3.593   (  -4.596   -4.513    0.775)    6.488
   3.966   (   0.683   -0.088   -0.157)    0.706
   4.114   (  -4.074    0.346   -0.268)    4.098
   5.051   (  -2.202   -0.440   -0.676)    2.345
   5.162   (   0.202   -0.366    1.296)    1.362
======================= Grid point 250 (57/144) =======================
q-point: (-0.36  0.14  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.977   ( -10.050   -1.459    0.171)   10.156
   1.333   ( -11.918    5.635    1.274)   13.244
   1.618   (  -7.909    2.166   -1.006)    8.261
   1.932   (  -4.533   -0.967   -1.689)    4.933
   2.207   (  -7.759   -3.908    0.547)    8.705
   3.025   (   3.971   -0.056   -0.201)    3.976
   3.178   (  -6.166    1.786    0.049)    6.420
   3.477   (  -4.005   -2.627    0.429)    4.809
   3.965   (  -0.860    0.310   -0.130)    0.923
   4.036   (  -4.504   -0.333   -0.245)    4.523
   4.999   (  -2.996   -0.426   -0.895)    3.156
   5.165   (   0.354   -0.168    1.752)    1.795
======================= Grid point 258 (58/144) =======================
q-point: ( 0.21  0.21  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.145   (   4.961    8.277   -1.060)    9.708
   1.548   (   3.890    8.354    0.843)    9.254
   1.753   (   3.763    4.208   -1.158)    5.762
   1.974   (  -0.021    1.419   -1.893)    2.365
   2.323   (   3.809    6.325    0.744)    7.421
   2.948   (  -2.191   -3.551    0.029)    4.173
   3.299   (   2.953    5.186    0.262)    5.974
   3.534   (   2.075    3.688    0.646)    4.281
   3.968   (  -0.236   -0.798   -0.007)    0.832
   4.118   (   2.320    3.918   -0.252)    4.560
   5.046   (   1.232    2.104   -0.664)    2.527
   5.158   (  -0.178   -0.214    1.284)    1.313
======================= Grid point 259 (59/144) =======================
q-point: ( 0.29  0.21  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.312   (   4.715    7.233   -1.802)    8.820
   1.639   (  -2.771    8.469    0.198)    8.913
   1.835   (   3.100   -0.519   -1.117)    3.336
   1.970   (  -2.627    3.364   -1.686)    4.589
   2.455   (   4.055    5.062    0.728)    6.526
   2.881   (  -1.126   -3.811    0.146)    3.977
   3.379   (  -0.508    6.202    0.429)    6.238
   3.627   (   4.858    2.089    0.613)    5.324
   3.939   (  -1.087   -2.470    0.180)    2.705
   4.189   (   1.350    3.387   -0.319)    3.660
   5.083   (   0.894    1.437   -0.444)    1.750
   5.154   (   0.047   -0.357    0.874)    0.945
======================= Grid point 260 (60/144) =======================
q-point: ( 0.36  0.21  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.441   (   1.395    5.792   -2.185)    6.346
   1.646   (  -4.245    6.478   -0.245)    7.748
   1.848   (   1.476   -3.997   -1.393)    4.483
   1.968   (  -3.617    7.055   -1.110)    8.005
   2.550   (   0.893    4.064    0.593)    4.203
   2.834   (   0.677   -4.088    0.284)    4.153
   3.412   (  -2.773    6.303    0.503)    6.905
   3.730   (   4.281    0.937    0.549)    4.417
   3.884   (  -1.111   -3.753    0.289)    3.924
   4.229   (  -0.284    2.567   -0.359)    2.607
   5.104   (   0.081    0.902   -0.270)    0.945
   5.152   (   0.310   -0.568    0.592)    0.877
======================= Grid point 261 (61/144) =======================
q-point: (-0.57  0.21  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.439   (  -5.887    1.740   -2.243)    6.536
   1.657   (  -2.734    7.325    0.146)    7.820
   1.830   (   0.885   -3.144   -1.965)    3.812
   1.977   (  -3.018    6.071   -0.873)    6.836
   2.550   (  -3.948    1.226    0.608)    4.178
   2.834   (   3.219   -2.621    0.286)    4.161
   3.412   (  -4.078    5.525    0.494)    6.884
   3.734   (  -2.778   -3.230    0.704)    4.318
   3.879   (   3.582   -0.887    0.125)    3.693
   4.229   (  -2.076    1.532   -0.354)    2.604
   5.104   (  -0.861    0.376   -0.277)    0.979
   5.153   (   0.392   -0.548    0.604)    0.905
======================= Grid point 262 (62/144) =======================
q-point: (-0.50  0.21  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.310   (  -8.674   -0.238   -1.878)    8.878
   1.655   (  -6.075    6.008    0.933)    8.594
   1.809   (  -1.103   -1.227   -2.208)    2.756
   1.984   (  -1.321    2.563   -1.234)    3.137
   2.454   (  -6.460   -0.983    0.725)    6.575
   2.882   (   3.929   -0.990    0.178)    4.055
   3.378   (  -5.078    3.528    0.424)    6.198
   3.629   (  -4.311   -3.387    0.712)    5.529
   3.936   (   2.657    0.039    0.044)    2.657
   4.189   (  -3.631    0.531   -0.318)    3.684
   5.083   (  -1.687   -0.009   -0.449)    1.746
   5.155   (   0.295   -0.291    0.890)    0.982
======================= Grid point 273 (63/144) =======================
q-point: ( 0.29  0.29  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.459   (   4.058    6.881   -2.294)    8.312
   1.755   (  -4.281   -0.212   -1.977)    4.720
   1.774   (   3.837   -0.839    0.268)    3.936
   2.041   (   1.201    2.528   -1.191)    3.041
   2.565   (   3.104    5.360    0.700)    6.233
   2.804   (  -2.130   -3.674    0.311)    4.258
   3.484   (   2.183    3.707    0.582)    4.342
   3.683   (   1.849    3.552    0.544)    4.041
   3.874   (  -2.024   -3.811    0.254)    4.323
   4.245   (   1.149    1.999   -0.388)    2.338
   5.109   (   0.616    1.018   -0.206)    1.208
   5.146   (  -0.269   -0.407    0.520)    0.713
======================= Grid point 274 (64/144) =======================
q-point: ( 0.36  0.29  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.556   (  -2.987    5.051   -2.493)    6.376
   1.685   (   1.763   -2.623   -2.352)    3.940
   1.804   (  -0.391    0.265    0.602)    0.765
   2.066   (  -1.321    2.394   -1.223)    2.995
   2.638   (  -2.394    4.143    0.601)    4.822
   2.751   (   2.235   -3.869    0.441)    4.490
   3.521   (  -2.112    3.626    0.634)    4.244
   3.747   (  -0.162    0.600    0.489)    0.791
   3.809   (   1.478   -2.854    0.318)    3.230
   4.266   (  -0.546    0.949   -0.415)    1.171
   5.121   (  -0.437    0.727   -0.033)    0.849
   5.140   (   0.376   -0.620    0.302)    0.785
======================= Grid point 422 (65/144) =======================
q-point: ( 0.00  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 144
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.328   (   0.000    0.000    6.294)    6.294
   0.328   (   0.000   -0.000    6.294)    6.294
   0.738   (   0.000    0.000   13.566)   13.566
   1.070   (   0.000    0.000   -1.364)    1.364
   1.070   (   0.000   -0.000   -1.364)    1.364
   2.783   (  -0.000    0.000   -0.345)    0.345
   2.783   (   0.000   -0.000   -0.345)    0.345
   3.327   (   0.000    0.000   -1.245)    1.245
   3.838   (  -0.000    0.000    0.108)    0.108
   3.838   (   0.000   -0.000    0.108)    0.108
   4.759   (   0.000   -0.000   -1.296)    1.296
   5.363   (   0.000   -0.000    0.686)    0.686
======================= Grid point 423 (66/144) =======================
q-point: ( 0.07  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.461   (   8.887    5.131    4.374)   11.155
   0.625   (   9.368    5.409    8.909)   14.014
   0.785   (  12.789    7.384    5.886)   15.897
   1.114   (   3.330    1.923   -1.325)    4.067
   1.285   (  15.322    8.846   -0.685)   17.705
   2.806   (   1.751    1.011   -0.362)    2.054
   2.907   (   8.151    4.706   -4.141)   10.283
   3.318   (  -0.559   -0.323   -1.323)    1.473
   3.841   (   0.503    0.290    0.138)    0.597
   3.856   (   1.315    0.759   -0.006)    1.518
   4.771   (   1.104    0.637   -1.246)    1.782
   5.327   (  -2.505   -1.446    2.089)    3.568
======================= Grid point 424 (67/144) =======================
q-point: ( 0.14  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.682   (   9.260    5.346    2.898)   11.079
   0.711   (   3.400    1.963    9.561)   10.336
   1.213   (   4.805    2.774   -1.193)    5.675
   1.246   (  19.689   11.367    0.187)   22.736
   1.674   (  16.333    9.430   -0.538)   18.867
   2.861   (   2.814    1.625   -0.294)    3.262
   3.067   (   4.498    2.597   -4.052)    6.588
   3.317   (   1.042    0.602   -1.465)    1.896
   3.861   (   1.229    0.710    0.090)    1.422
   3.898   (   2.202    1.271   -0.146)    2.547
   4.813   (   2.434    1.406   -1.288)    3.092
   5.270   (  -2.063   -1.191    2.975)    3.811
======================= Grid point 425 (68/144) =======================
q-point: ( 0.21  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.820   (   5.829    3.366    5.596)    8.753
   0.880   (   7.409    4.278    2.230)    8.841
   1.325   (   4.395    2.538   -1.112)    5.195
   1.637   (  13.151    7.593   -1.501)   15.259
   2.016   (  12.461    7.195   -0.101)   14.390
   2.933   (   3.274    1.890   -0.120)    3.782
   3.112   (  -0.167   -0.096   -1.556)    1.568
   3.386   (   4.603    2.657   -0.804)    5.375
   3.896   (   1.678    0.969   -0.037)    1.938
   3.957   (   2.764    1.596   -0.258)    3.202
   4.881   (   3.191    1.843   -1.265)    3.896
   5.233   (  -1.170   -0.675    2.701)    3.020
======================= Grid point 426 (69/144) =======================
q-point: ( 0.29  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.975   (   7.108    4.104    2.127)    8.479
   1.028   (   5.200    3.002    1.938)    6.310
   1.411   (   2.848    1.644   -1.138)    3.480
   1.862   (   6.263    3.616   -2.640)    7.698
   2.255   (   7.878    4.548    0.430)    9.107
   3.010   (   3.175    1.833    0.084)    3.668
   3.089   (  -1.375   -0.794   -0.518)    1.670
   3.515   (   6.027    3.480    0.537)    6.980
   3.937   (   1.690    0.976   -0.183)    1.960
   4.022   (   2.611    1.507   -0.323)    3.032
   4.957   (   3.087    1.782   -1.095)    3.729
   5.213   (  -0.614   -0.354    2.187)    2.299
======================= Grid point 427 (70/144) =======================
q-point: ( 0.36  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.127   (   3.249    1.876    1.793)    4.158
   1.143   (   6.845    3.952   -0.197)    7.906
   1.458   (   1.222    0.706   -1.227)    1.870
   1.947   (   1.275    0.736   -3.101)    3.433
   2.397   (   4.499    2.597    0.817)    5.259
   3.056   (  -1.379   -0.796   -0.228)    1.609
   3.079   (   2.539    1.466    0.227)    2.941
   3.655   (   5.619    3.244    1.194)    6.598
   3.973   (   1.297    0.749   -0.278)    1.523
   4.073   (   1.549    0.894   -0.291)    1.812
   5.021   (   2.295    1.325   -0.896)    2.797
   5.202   (  -0.259   -0.150    1.704)    1.730
======================= Grid point 428 (71/144) =======================
q-point: ( 0.43  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.184   (   1.662    0.960    1.676)    2.548
   1.282   (   4.622    2.669   -1.471)    5.536
   1.474   (   0.283    0.163   -1.305)    1.345
   1.942   (  -1.122   -0.648   -3.079)    3.341
   2.478   (   2.489    1.437    0.775)    2.976
   3.026   (  -1.070   -0.618   -0.045)    1.236
   3.126   (   1.430    0.825    0.287)    1.676
   3.767   (   3.640    2.101    1.230)    4.379
   3.996   (   0.640    0.369   -0.302)    0.798
   4.093   (   0.301    0.174   -0.149)    0.378
   5.062   (   1.141    0.659   -0.723)    1.503
   5.200   (  -0.006   -0.004    1.315)    1.315
======================= Grid point 429 (72/144) =======================
q-point: (-0.50  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.206   (   0.233    0.135    1.597)    1.619
   1.340   (  -0.069   -0.040   -1.872)    1.874
   1.477   (   0.052    0.030   -1.351)    1.352
   1.922   (  -0.159   -0.092   -2.966)    2.971
   2.511   (   0.231    0.133    0.593)    0.650
   3.011   (  -0.048   -0.028    0.096)    0.111
   3.144   (   0.040    0.023    0.295)    0.299
   3.813   (   0.050    0.029    1.124)    1.126
   4.002   (  -0.121   -0.070   -0.284)    0.317
   4.094   (  -0.092   -0.053   -0.033)    0.111
   5.075   (  -0.034   -0.020   -0.630)    0.631
   5.201   (   0.061    0.035    1.132)    1.134
======================= Grid point 431 (73/144) =======================
q-point: (-0.43  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.194   (  -1.276   -0.736    1.606)    2.179
   1.279   (  -4.745   -2.739   -1.424)    5.661
   1.476   (  -0.280   -0.162   -1.361)    1.399
   1.938   (   1.189    0.686   -3.063)    3.357
   2.485   (  -2.402   -1.387    0.699)    2.860
   3.026   (   1.231    0.710   -0.023)    1.421
   3.128   (  -1.348   -0.778    0.271)    1.580
   3.768   (  -3.648   -2.106    1.180)    4.375
   3.991   (  -0.816   -0.471   -0.252)    0.975
   4.089   (  -0.532   -0.307   -0.085)    0.620
   5.061   (  -1.154   -0.666   -0.686)    1.499
   5.202   (   0.080    0.046    1.264)    1.267
======================= Grid point 432 (74/144) =======================
q-point: (-0.36  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.138   (  -6.877   -3.970   -0.105)    7.941
   1.144   (  -3.103   -1.792    1.686)    3.960
   1.457   (  -1.438   -0.830   -1.291)    2.103
   1.952   (  -0.645   -0.372   -3.221)    3.306
   2.400   (  -4.894   -2.825    0.815)    5.709
   3.063   (   1.735    1.002   -0.267)    2.021
   3.082   (  -2.492   -1.439    0.194)    2.884
   3.654   (  -5.751   -3.320    1.172)    6.743
   3.965   (  -1.283   -0.741   -0.206)    1.495
   4.061   (  -1.880   -1.086   -0.165)    2.177
   5.020   (  -2.248   -1.298   -0.863)    2.736
   5.206   (   0.300    0.173    1.639)    1.675
======================= Grid point 433 (75/144) =======================
q-point: (-0.29  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.970   (  -7.120   -4.111    2.248)    8.523
   1.044   (  -5.421   -3.130    1.854)    6.528
   1.403   (  -3.194   -1.844   -1.147)    3.863
   1.881   (  -5.830   -3.366   -2.927)    7.340
   2.249   (  -8.209   -4.739    0.615)    9.499
   3.014   (  -3.251   -1.877    0.051)    3.754
   3.102   (   1.493    0.862   -0.670)    1.850
   3.513   (  -5.994   -3.460    0.571)    6.945
   3.933   (  -1.433   -0.828   -0.143)    1.661
   4.003   (  -2.902   -1.675   -0.133)    3.353
   4.957   (  -3.051   -1.761   -1.081)    3.685
   5.217   (   0.689    0.398    2.106)    2.252
======================= Grid point 434 (76/144) =======================
q-point: (-0.21  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.814   (  -5.948   -3.434    5.747)    8.956
   0.887   (  -7.948   -4.589    2.243)    9.447
   1.310   (  -4.631   -2.674   -1.039)    5.447
   1.660   ( -13.201   -7.622   -1.775)   15.346
   2.004   ( -12.695   -7.329    0.300)   14.662
   2.933   (  -3.478   -2.008   -0.118)    4.018
   3.125   (   0.141    0.082   -1.768)    1.775
   3.387   (  -4.443   -2.565   -0.818)    5.195
   3.900   (  -1.325   -0.765   -0.075)    1.532
   3.934   (  -2.756   -1.591   -0.028)    3.182
   4.882   (  -3.188   -1.841   -1.257)    3.890
   5.239   (   1.204    0.695    2.579)    2.930
======================= Grid point 435 (77/144) =======================
q-point: (-0.14  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.675   (  -9.822   -5.670    3.066)   11.748
   0.700   (  -3.538   -2.043    9.748)   10.570
   1.196   (  -4.719   -2.724   -1.075)    5.554
   1.267   ( -19.743  -11.399   -0.033)   22.797
   1.659   ( -16.168   -9.335    0.148)   18.670
   2.856   (  -2.960   -1.709   -0.242)    3.426
   3.080   (  -4.550   -2.627   -4.292)    6.785
   3.319   (  -1.117   -0.645   -1.498)    1.977
   3.871   (  -1.088   -0.628   -0.015)    1.256
   3.879   (  -1.888   -1.090    0.041)    2.181
   4.814   (  -2.455   -1.417   -1.286)    3.113
   5.276   (   1.935    1.117    2.833)    3.608
======================= Grid point 436 (78/144) =======================
q-point: (-0.07  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.446   (  -8.878   -5.126    4.704)   11.279
   0.619   (  -7.947   -4.588    9.350)   13.100
   0.801   ( -14.147   -8.168    5.390)   17.201
   1.103   (  -2.931   -1.692   -1.238)    3.604
   1.283   ( -14.623   -8.443   -0.252)   16.887
   2.801   (  -1.633   -0.943   -0.304)    1.910
   2.915   (  -8.519   -4.918   -4.307)   10.738
   3.319   (   0.532    0.307   -1.325)    1.461
   3.846   (  -0.913   -0.527    0.091)    1.058
   3.850   (  -0.749   -0.433    0.045)    0.866
   4.772   (  -1.122   -0.648   -1.247)    1.798
   5.329   (   2.350    1.357    2.030)    3.389
======================= Grid point 438 (79/144) =======================
q-point: ( 0.07  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.648   (   7.014   10.449    3.611)   13.093
   0.683   (   2.198    3.727   10.035)   10.928
   1.088   (   8.687   15.839    0.332)   18.068
   1.225   (   6.070   10.115   -1.054)   11.844
   1.541   (   9.243   14.870    0.056)   17.509
   2.851   (   2.267    3.377   -0.324)    4.080
   3.034   (   2.695    5.823   -4.808)    8.018
   3.313   (   0.055    0.186   -1.433)    1.446
   3.859   (   0.372    1.534    0.091)    1.581
   3.877   (   1.498    1.198   -0.063)    1.919
   4.799   (   1.195    2.118   -1.269)    2.744
   5.286   (  -1.462   -2.143    2.837)    3.844
======================= Grid point 439 (80/144) =======================
q-point: ( 0.14  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.761   (   4.301    3.329    6.835)    8.735
   0.889   (   4.221   13.350    2.385)   14.203
   1.310   (   5.059    6.873   -1.170)    8.614
   1.517   (  10.081   13.146   -0.983)   16.596
   1.861   (  12.381    8.757    0.089)   15.165
   2.936   (   1.861    5.491   -0.285)    5.805
   3.098   (   1.095    0.058   -2.409)    2.647
   3.343   (   2.512    2.005   -1.187)    3.426
   3.890   (   0.884    2.082    0.018)    2.262
   3.925   (   2.496    1.361   -0.182)    2.849
   4.855   (   2.268    2.607   -1.284)    3.686
   5.245   (  -1.215   -1.237    2.823)    3.313
======================= Grid point 440 (81/144) =======================
q-point: ( 0.21  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.893   (   6.250    4.058    3.192)    8.107
   1.068   (   0.966   12.846    1.861)   13.016
   1.442   (   1.298    8.228   -1.054)    8.396
   1.775   (   7.821    5.654   -2.171)    9.892
   2.134   (  10.543    4.665    0.315)   11.534
   3.012   (   0.309    4.615   -0.168)    4.628
   3.096   (   0.137   -0.550   -0.694)    0.897
   3.440   (   5.450    2.520    0.005)    6.004
   3.929   (   1.013    2.076   -0.112)    2.312
   3.989   (   2.913    1.671   -0.272)    3.369
   4.927   (   2.705    2.509   -1.152)    3.866
   5.218   (  -0.688   -0.763    2.367)    2.581
======================= Grid point 441 (82/144) =======================
q-point: ( 0.29  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.054   (   6.977    4.038    0.464)    8.074
   1.187   (  -1.282   11.334    1.615)   11.520
   1.524   (  -1.712    8.633   -1.056)    8.864
   1.908   (   3.738    0.873   -2.873)    4.795
   2.325   (   6.887    2.764    0.666)    7.451
   3.021   (   0.011   -2.138   -0.084)    2.140
   3.123   (  -0.273    5.136   -0.004)    5.143
   3.573   (   6.535    2.200    0.934)    6.959
   3.963   (   0.847    1.264   -0.210)    1.535
   4.055   (   2.214    1.961   -0.310)    2.974
   4.996   (   2.406    2.000   -0.947)    3.269
   5.203   (  -0.275   -0.572    1.866)    1.971
======================= Grid point 442 (83/144) =======================
q-point: ( 0.36  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.209   (   5.784    3.575   -1.064)    6.882
   1.258   (  -2.536    9.431    1.355)    9.860
   1.562   (  -3.704    8.580   -1.173)    9.419
   1.936   (   0.448   -1.674   -3.025)    3.487
   2.442   (   3.810    2.139    0.820)    4.446
   2.995   (   0.604   -3.540   -0.021)    3.591
   3.168   (  -1.078    5.502    0.231)    5.611
   3.703   (   5.694    1.597    1.191)    6.033
   3.981   (   0.543   -0.390   -0.215)    0.702
   4.102   (   0.633    2.267   -0.271)    2.370
   5.048   (   1.520    1.307   -0.755)    2.142
   5.196   (   0.073   -0.478    1.425)    1.505
======================= Grid point 443 (84/144) =======================
q-point: ( 0.43  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.289   (  -3.919    8.318    1.326)    9.290
   1.317   (   2.698    1.628   -1.888)    3.673
   1.573   (  -4.647    8.337   -1.301)    9.633
   1.913   (  -0.110   -2.083   -2.891)    3.564
   2.506   (   1.279    1.367    0.651)    1.982
   2.975   (   1.225   -3.411    0.104)    3.626
   3.193   (  -1.974    4.757    0.324)    5.161
   3.793   (   2.754    0.619    1.121)    3.038
   3.976   (   0.665   -2.221   -0.155)    2.323
   4.124   (  -0.911    2.534   -0.206)    2.701
   5.074   (   0.373    0.598   -0.611)    0.933
   5.195   (   0.298   -0.443    1.121)    1.242
======================= Grid point 444 (85/144) =======================
q-point: (-0.50  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.279   (  -4.641    3.872    0.392)    6.056
   1.334   (  -2.939    2.598   -0.982)    4.044
   1.570   (  -4.974    7.656   -1.303)    9.223
   1.911   (   1.659   -0.833   -2.882)    3.428
   2.510   (  -1.494   -0.347    0.614)    1.652
   2.974   (   2.347   -2.730    0.115)    3.602
   3.194   (  -3.073    4.098    0.321)    5.133
   3.796   (  -1.838   -1.853    1.088)    2.827
   3.971   (   1.319   -2.039   -0.115)    2.431
   4.123   (  -1.873    2.073   -0.188)    2.800
   5.074   (  -0.755   -0.049   -0.589)    0.959
   5.196   (   0.339   -0.434    1.094)    1.225
======================= Grid point 445 (86/144) =======================
q-point: (-0.43  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.190   (  -5.888   -2.496   -0.617)    6.425
   1.291   (  -6.983    6.559    0.888)    9.621
   1.553   (  -5.289    6.624   -1.208)    8.563
   1.938   (   1.442   -0.689   -3.074)    3.465
   2.447   (  -3.867   -2.515    0.782)    4.679
   2.999   (   3.032   -2.119   -0.031)    3.699
   3.170   (  -4.127    3.679    0.217)    5.533
   3.705   (  -4.435   -4.081    1.143)    6.135
   3.971   (   0.118   -0.837   -0.132)    0.856
   4.098   (  -2.638    0.674   -0.209)    2.731
   5.047   (  -1.872   -0.641   -0.719)    2.106
   5.199   (   0.430   -0.319    1.370)    1.471
======================= Grid point 446 (87/144) =======================
q-point: (-0.36  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.041   (  -6.896   -3.865    0.621)    7.930
   1.215   (  -9.478    6.729    1.497)   11.720
   1.511   (  -6.401    4.957   -1.099)    8.170
   1.922   (  -2.151   -2.270   -3.095)    4.400
   2.322   (  -6.142   -4.854    0.766)    7.865
   3.029   (   1.833   -1.052   -0.146)    2.119
   3.128   (  -4.012    2.706   -0.058)    4.840
   3.572   (  -5.331   -4.569    0.943)    7.083
   3.956   (  -1.258   -0.008   -0.145)    1.266
   4.043   (  -3.338   -0.772   -0.188)    3.431
   4.996   (  -2.893   -1.061   -0.929)    3.219
   5.207   (   0.706   -0.100    1.810)    1.945
======================= Grid point 447 (88/144) =======================
q-point: (-0.29  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.882   (  -6.658   -3.210    3.243)    8.071
   1.088   ( -12.069    5.268    1.935)   13.310
   1.425   (  -7.687    2.690   -1.071)    8.214
   1.797   (  -8.734   -4.058   -2.448)    9.937
   2.126   (  -9.505   -6.729    0.541)   11.658
   3.013   (  -4.454    1.833   -0.156)    4.819
   3.107   (   0.685   -0.635   -0.861)    1.271
   3.439   (  -4.794   -3.396    0.030)    5.875
   3.930   (  -1.938    0.518   -0.125)    2.010
   3.971   (  -3.327   -1.418   -0.091)    3.618
   4.927   (  -3.504   -1.101   -1.145)    3.847
   5.222   (   1.113    0.085    2.297)    2.554
======================= Grid point 448 (89/144) =======================
q-point: (-0.21  0.07  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.751   (  -5.251   -1.423    6.745)    8.666
   0.896   ( -14.050    2.539    2.659)   14.523
   1.297   (  -8.393   -0.053   -1.179)    8.476
   1.532   ( -16.537   -3.286   -1.149)   16.899
   1.853   ( -13.962   -5.720    0.391)   15.094
   2.933   (  -5.769    1.194   -0.256)    5.897
   3.106   (  -0.535   -1.525   -2.542)    3.012
   3.344   (  -2.934   -1.250   -1.204)    3.409
   3.899   (  -2.311    0.351   -0.057)    2.338
   3.907   (  -2.225   -0.775   -0.021)    2.356
   4.856   (  -3.388   -0.704   -1.281)    3.690
   5.248   (   1.727    0.128    2.749)    3.249
======================= Grid point 453 (90/144) =======================
q-point: ( 0.14  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.853   (   3.699    5.717    3.630)    7.717
   1.161   (   6.245   12.253    1.851)   13.876
   1.460   (   5.228    7.291   -1.279)    9.063
   1.751   (   4.279    9.302   -2.009)   10.435
   2.038   (   5.525    8.880    0.579)   10.474
   3.062   (   4.170    5.644   -0.172)    7.019
   3.075   (  -2.106   -1.396   -1.025)    2.727
   3.399   (   2.035    3.462   -0.218)    4.022
   3.939   (   0.586    2.005   -0.094)    2.091
   3.959   (   2.602    2.567   -0.160)    3.659
   4.915   (   1.815    3.182   -1.157)    3.842
   5.222   (  -0.728   -1.011    2.393)    2.698
======================= Grid point 454 (91/144) =======================
q-point: ( 0.21  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.993   (   5.326    5.740    0.850)    7.876
   1.345   (  -0.095   13.165    1.379)   13.238
   1.613   (   3.201    7.441   -1.210)    8.190
   1.870   (   1.930    3.821   -2.484)    4.949
   2.225   (   7.413    4.769    0.649)    8.838
   3.015   (  -1.870   -2.871   -0.220)    3.434
   3.177   (   1.537    6.678    0.059)    6.853
   3.489   (   4.751    2.477    0.634)    5.395
   3.965   (   0.479    1.005   -0.144)    1.123
   4.036   (   2.683    3.186   -0.259)    4.173
   4.980   (   2.018    2.645   -0.943)    3.458
   5.202   (  -0.448   -0.783    1.919)    2.120
======================= Grid point 455 (92/144) =======================
q-point: ( 0.29  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.149   (   5.679    5.354   -1.027)    7.872
   1.437   (  -3.638   12.085    1.039)   12.663
   1.715   (  -1.010    8.841   -1.144)    8.971
   1.905   (   0.488   -0.730   -2.711)    2.849
   2.384   (   5.690    3.429    0.783)    6.690
   2.958   (  -0.703   -3.554    0.031)    3.623
   3.251   (  -0.838    6.835    0.283)    6.892
   3.609   (   6.218    1.529    1.002)    6.481
   3.972   (  -0.016   -0.561   -0.114)    0.572
   4.110   (   1.719    3.376   -0.324)    3.802
   5.036   (   1.635    1.937   -0.726)    2.637
   5.190   (  -0.108   -0.734    1.452)    1.631
======================= Grid point 456 (93/144) =======================
q-point: ( 0.36  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.288   (   3.967    4.436   -2.028)    6.287
   1.473   (  -4.999   10.368    0.760)   11.536
   1.761   (  -4.490   10.007   -1.247)   11.038
   1.884   (   0.026   -3.246   -2.592)    4.154
   2.492   (   2.750    2.943    0.736)    4.095
   2.919   (   0.567   -3.838    0.135)    3.882
   3.290   (  -2.284    6.437    0.424)    6.844
   3.727   (   5.368    0.950    1.021)    5.546
   3.951   (  -0.545   -2.377   -0.015)    2.438
   4.161   (   0.151    3.181   -0.350)    3.204
   5.074   (   0.741    1.236   -0.541)    1.539
   5.183   (   0.251   -0.759    1.067)    1.333
======================= Grid point 457 (94/144) =======================
q-point: ( 0.43  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.350   (  -1.499    1.866   -2.262)    3.293
   1.486   (  -4.599    9.634    0.493)   10.687
   1.761   (  -6.096    8.499   -1.351)   10.546
   1.869   (   1.488   -1.457   -2.383)    3.165
   2.535   (  -0.730    1.487    0.618)    1.768
   2.904   (   1.968   -3.427    0.223)    3.959
   3.303   (  -3.369    5.853    0.472)    6.770
   3.794   (   0.826   -0.749    0.986)    1.488
   3.921   (   1.106   -2.769    0.037)    2.982
   4.179   (  -1.475    2.562   -0.351)    2.977
   5.086   (  -0.374    0.583   -0.442)    0.822
   5.182   (   0.497   -0.758    0.886)    1.267
======================= Grid point 459 (95/144) =======================
q-point: (-0.50  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.280   (  -6.119   -0.884   -1.879)    6.462
   1.497   (  -5.871    9.633    0.574)   11.296
   1.732   (  -5.691    6.441   -1.241)    8.684
   1.899   (   1.921    0.362   -2.599)    3.252
   2.495   (  -3.857   -1.093    0.726)    4.074
   2.920   (   3.181   -2.340    0.133)    3.951
   3.291   (  -4.377    5.174    0.420)    6.790
   3.733   (  -3.610   -4.139    0.961)    5.576
   3.942   (   2.209   -0.718    0.061)    2.324
   4.161   (  -2.908    1.489   -0.336)    3.284
   5.073   (  -1.453   -0.015   -0.513)    1.542
   5.185   (   0.588   -0.611    1.031)    1.335
======================= Grid point 460 (96/144) =======================
q-point: (-0.43  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.141   (  -7.605   -1.843   -0.937)    7.881
   1.465   (  -8.697    8.863    1.000)   12.457
   1.684   (  -6.290    4.384   -1.204)    7.761
   1.923   (   0.113   -0.077   -2.815)    2.818
   2.383   (  -5.984   -3.415    0.832)    6.940
   2.963   (   3.642   -1.259   -0.002)    3.853
   3.253   (  -5.418    4.094    0.267)    6.796
   3.611   (  -4.660   -4.698    0.975)    6.689
   3.963   (   0.552    0.003   -0.042)    0.554
   4.108   (  -4.017    0.317   -0.289)    4.040
   5.036   (  -2.469   -0.422   -0.707)    2.603
   5.191   (   0.721   -0.347    1.418)    1.628
======================= Grid point 461 (97/144) =======================
q-point: (-0.36  0.14  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.987   (  -7.704   -1.247    0.894)    7.855
   1.362   ( -11.571    5.472    1.440)   12.880
   1.594   (  -7.653    1.966   -1.277)    8.004
   1.886   (  -4.200   -0.748   -2.660)    5.028
   2.221   (  -7.960   -3.803    0.758)    8.854
   3.020   (   3.441   -0.208   -0.262)    3.457
   3.179   (  -6.489    1.850    0.036)    6.747
   3.488   (  -4.583   -2.815    0.648)    5.417
   3.962   (  -0.880    0.364   -0.119)    0.959
   4.031   (  -4.441   -0.278   -0.206)    4.455
   4.980   (  -3.283   -0.432   -0.938)    3.441
   5.204   (   0.943   -0.136    1.896)    2.122
======================= Grid point 469 (98/144) =======================
q-point: ( 0.21  0.21  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.123   (   4.233    6.937   -0.960)    8.184
   1.571   (   3.260    8.502    1.268)    9.193
   1.720   (   3.994    3.089   -2.023)    5.440
   1.929   (   0.067    1.929   -2.227)    2.947
   2.340   (   3.840    6.412    0.855)    7.523
   2.949   (  -2.113   -3.340    0.037)    3.952
   3.305   (   3.091    5.439    0.337)    6.266
   3.550   (   2.115    3.782    0.867)    4.419
   3.967   (  -0.120   -0.707   -0.052)    0.719
   4.111   (   2.274    3.824   -0.318)    4.460
   5.031   (   1.370    2.355   -0.697)    2.812
   5.187   (  -0.453   -0.689    1.439)    1.658
======================= Grid point 470 (99/144) =======================
q-point: ( 0.29  0.21  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.270   (   4.355    6.423   -2.131)    8.047
   1.645   (  -3.500    7.763    0.270)    8.520
   1.801   (   3.727   -1.257   -1.943)    4.387
   1.935   (  -2.609    4.368   -1.641)    5.346
   2.472   (   3.919    5.109    0.856)    6.496
   2.884   (  -1.066   -3.637    0.179)    3.794
   3.389   (  -0.481    6.438    0.544)    6.479
   3.643   (   4.853    1.946    0.874)    5.301
   3.942   (  -0.977   -2.319    0.073)    2.517
   4.181   (   1.324    3.287   -0.410)    3.567
   5.074   (   0.998    1.702   -0.468)    2.028
   5.174   (  -0.166   -0.789    1.005)    1.288
======================= Grid point 471 (100/144) =======================
q-point: ( 0.36  0.21  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.388   (   1.360    5.289   -2.686)    6.086
   1.635   (  -3.591    4.315   -0.956)    5.695
   1.819   (   0.676   -2.041   -0.992)    2.368
   1.940   (  -3.315    6.874   -1.563)    7.790
   2.564   (   0.730    4.110    0.738)    4.240
   2.841   (   0.708   -3.897    0.321)    3.974
   3.424   (  -2.837    6.487    0.633)    7.108
   3.745   (   4.381    0.797    0.823)    4.529
   3.889   (  -1.172   -3.632    0.187)    3.821
   4.220   (  -0.253    2.470   -0.462)    2.526
   5.098   (   0.073    1.167   -0.283)    1.203
   5.166   (   0.293   -0.924    0.690)    1.190
======================= Grid point 472 (101/144) =======================
q-point: (-0.57  0.21  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.386   (  -5.457    1.505   -2.648)    6.249
   1.648   (  -0.774    5.595   -1.196)    5.774
   1.799   (  -0.854   -1.446   -0.673)    1.809
   1.950   (  -2.972    5.786   -1.691)    6.721
   2.565   (  -3.885    1.387    0.744)    4.191
   2.841   (   3.046   -2.552    0.319)    3.987
   3.424   (  -4.206    5.663    0.637)    7.083
   3.750   (  -2.672   -3.397    0.781)    4.392
   3.883   (   3.451   -0.752    0.227)    3.539
   4.221   (  -2.001    1.456   -0.461)    2.517
   5.098   (  -1.078    0.516   -0.271)    1.225
   5.167   (   0.704   -0.713    0.676)    1.209
======================= Grid point 473 (102/144) =======================
q-point: (-0.50  0.21  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.267   (  -7.811   -0.292   -2.102)    8.094
   1.673   (  -5.236    6.801    0.490)    8.597
   1.762   (  -0.837   -2.351   -2.015)    3.208
   1.951   (  -1.980    2.723   -1.848)    3.841
   2.471   (  -6.432   -0.868    0.877)    6.549
   2.885   (   3.772   -0.976    0.171)    3.900
   3.389   (  -5.286    3.620    0.545)    6.430
   3.646   (  -4.207   -3.509    0.839)    5.542
   3.937   (   2.464    0.116    0.113)    2.470
   4.181   (  -3.535    0.497   -0.406)    3.593
   5.073   (  -1.968    0.022   -0.455)    2.020
   5.175   (   0.778   -0.318    0.988)    1.297
======================= Grid point 484 (103/144) =======================
q-point: ( 0.29  0.29  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.404   (   3.821    6.456   -2.814)    8.013
   1.696   (  -2.900   -3.576   -3.217)    5.616
   1.793   (   2.402    2.379    1.265)    3.610
   2.010   (   1.110    2.425   -1.750)    3.190
   2.582   (   3.034    5.257    0.838)    6.128
   2.812   (  -2.016   -3.477    0.361)    4.035
   3.498   (   2.267    3.831    0.740)    4.513
   3.696   (   1.725    3.433    0.681)    3.902
   3.880   (  -1.899   -3.683    0.273)    4.153
   4.236   (   1.109    1.936   -0.502)    2.287
   5.104   (   0.764    1.285   -0.197)    1.508
   5.158   (  -0.476   -0.774    0.591)    1.084
======================= Grid point 485 (104/144) =======================
q-point: ( 0.36  0.29  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.496   (  -2.849    4.823   -3.101)    6.403
   1.625   (   2.162   -3.478   -3.208)    5.202
   1.824   (  -0.834    1.168    1.301)    1.937
   2.033   (  -1.218    2.230   -1.875)    3.158
   2.652   (  -2.353    4.083    0.725)    4.768
   2.761   (   2.136   -3.704    0.521)    4.307
   3.537   (  -2.193    3.756    0.805)    4.424
   3.758   (   0.008    0.413    0.605)    0.732
   3.817   (   1.352   -2.736    0.352)    3.072
   4.256   (  -0.523    0.911   -0.539)    1.181
   5.121   (  -0.581    0.982   -0.011)    1.140
   5.147   (   0.539   -0.906    0.354)    1.112
======================= Grid point 633 (105/144) =======================
q-point: ( 0.00  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 144
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.418   (   0.000    0.000    2.436)    2.436
   0.418   (   0.000   -0.000    2.436)    2.436
   0.928   (  -0.000    0.000    4.978)    4.978
   1.047   (   0.000    0.000   -0.718)    0.718
   1.047   (   0.000   -0.000   -0.718)    0.718
   2.777   (  -0.000    0.000   -0.152)    0.152
   2.777   (   0.000   -0.000   -0.152)    0.152
   3.307   (   0.000   -0.000   -0.605)    0.605
   3.839   (   0.000    0.000    0.047)    0.047
   3.839   (   0.000   -0.000    0.047)    0.047
   4.739   (   0.000    0.000   -0.565)    0.565
   5.374   (   0.000   -0.000    0.290)    0.290
======================= Grid point 634 (106/144) =======================
q-point: ( 0.07  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.525   (   7.699    4.445    1.673)    9.046
   0.743   (  19.636   11.337    3.037)   22.876
   0.876   (  -0.185   -0.107    2.473)    2.482
   1.089   (   3.327    1.921   -0.913)    3.949
   1.280   (  15.196    8.773    0.107)   17.547
   2.799   (   1.737    1.003   -0.255)    2.022
   2.863   (   6.382    3.685   -0.805)    7.413
   3.298   (  -0.595   -0.344   -0.598)    0.911
   3.845   (   0.549    0.317    0.222)    0.672
   3.854   (   1.133    0.654   -0.191)    1.322
   4.753   (   1.147    0.662   -0.533)    1.428
   5.354   (  -1.547   -0.893    0.638)    1.897
======================= Grid point 635 (107/144) =======================
q-point: ( 0.14  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.727   (   8.753    5.053    1.258)   10.185
   0.858   (   1.090    0.629    4.273)    4.454
   1.190   (   4.876    2.815   -0.966)    5.712
   1.241   (  18.829   10.871   -0.323)   21.744
   1.666   (  16.409    9.474   -0.195)   18.948
   2.855   (   2.926    1.689   -0.215)    3.385
   3.019   (   5.678    3.278   -0.952)    6.625
   3.296   (   1.091    0.630   -0.558)    1.378
   3.865   (   1.140    0.658    0.224)    1.335
   3.892   (   2.054    1.186   -0.419)    2.409
   4.795   (   2.423    1.399   -0.505)    2.843
   5.311   (  -1.827   -1.055    1.088)    2.374
======================= Grid point 636 (108/144) =======================
q-point: ( 0.21  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.910   (   3.338    1.927    2.646)    4.675
   0.917   (   7.238    4.179    1.228)    8.447
   1.303   (   4.458    2.574   -0.896)    5.225
   1.611   (  12.332    7.120   -0.679)   14.256
   2.014   (  12.791    7.385   -0.106)   14.770
   2.932   (   3.451    1.992   -0.053)    3.985
   3.093   (   0.787    0.455   -0.278)    0.951
   3.375   (   5.440    3.141   -0.275)    6.287
   3.896   (   1.509    0.871    0.050)    1.743
   3.949   (   2.700    1.559   -0.546)    3.165
   4.863   (   3.250    1.877   -0.486)    3.785
   5.273   (  -1.391   -0.803    1.151)    1.976
======================= Grid point 637 (109/144) =======================
q-point: ( 0.29  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.012   (   5.137    2.966    1.127)    6.038
   1.063   (   5.132    2.963    1.284)    6.063
   1.390   (   2.860    1.651   -0.783)    3.394
   1.821   (   5.761    3.326   -1.090)    6.741
   2.261   (   8.235    4.754    0.145)    9.510
   3.012   (   3.313    1.913    0.101)    3.826
   3.084   (  -1.099   -0.634    0.053)    1.270
   3.524   (   6.714    3.876    0.240)    7.757
   3.933   (   1.553    0.897   -0.157)    1.800
   4.013   (   2.646    1.527   -0.511)    3.097
   4.940   (   3.199    1.847   -0.438)    3.720
   5.246   (  -0.927   -0.535    0.983)    1.454
======================= Grid point 638 (110/144) =======================
q-point: ( 0.36  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.141   (   5.620    3.244    0.001)    6.489
   1.160   (   3.153    1.820    1.252)    3.850
   1.437   (   1.192    0.688   -0.710)    1.548
   1.897   (   1.097    0.633   -1.442)    1.919
   2.410   (   4.692    2.709    0.434)    5.435
   3.054   (  -1.358   -0.784    0.040)    1.569
   3.083   (   2.596    1.499    0.169)    3.002
   3.674   (   5.821    3.361    0.544)    6.743
   3.967   (   1.256    0.725   -0.260)    1.473
   4.065   (   1.673    0.966   -0.357)    1.965
   5.007   (   2.406    1.389   -0.380)    2.804
   5.229   (  -0.547   -0.316    0.777)    1.002
======================= Grid point 639 (111/144) =======================
q-point: ( 0.43  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.214   (   1.520    0.877    1.089)    2.065
   1.260   (   4.067    2.348   -0.612)    4.735
   1.452   (   0.237    0.137   -0.707)    0.758
   1.892   (  -1.021   -0.589   -1.537)    1.936
   2.491   (   2.330    1.345    0.489)    2.735
   3.025   (  -0.975   -0.563   -0.027)    1.127
   3.131   (   1.432    0.827    0.166)    1.662
   3.787   (   3.567    2.060    0.593)    4.162
   3.990   (   0.677    0.391   -0.232)    0.815
   4.090   (   0.452    0.261   -0.175)    0.550
   5.051   (   1.235    0.713   -0.319)    1.461
   5.220   (  -0.217   -0.125    0.598)    0.649
======================= Grid point 640 (112/144) =======================
q-point: (-0.50  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.233   (   0.111    0.064    0.855)    0.864
   1.312   (  -0.029   -0.017   -0.762)    0.763
   1.454   (   0.015    0.009   -0.748)    0.748
   1.874   (  -0.075   -0.043   -1.385)    1.388
   2.520   (   0.104    0.060    0.267)    0.293
   3.012   (  -0.019   -0.011    0.024)    0.032
   3.148   (   0.018    0.010    0.134)    0.135
   3.831   (   0.016    0.009    0.535)    0.535
   3.997   (  -0.054   -0.031   -0.128)    0.142
   4.093   (  -0.037   -0.022   -0.046)    0.063
   5.065   (  -0.011   -0.006   -0.269)    0.270
   5.218   (   0.023    0.013    0.487)    0.487
======================= Grid point 642 (113/144) =======================
q-point: (-0.43  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.219   (  -1.334   -0.770    0.642)    1.669
   1.259   (  -4.115   -2.376   -0.487)    4.776
   1.452   (  -0.250   -0.145   -0.755)    0.808
   1.890   (   1.041    0.601   -1.348)    1.806
   2.495   (  -2.291   -1.323    0.171)    2.651
   3.025   (   1.054    0.609   -0.037)    1.218
   3.132   (  -1.395   -0.805    0.086)    1.613
   3.787   (  -3.576   -2.065    0.549)    4.166
   3.988   (  -0.756   -0.436   -0.017)    0.873
   4.088   (  -0.553   -0.319    0.011)    0.639
   5.050   (  -1.238   -0.714   -0.283)    1.457
   5.221   (   0.247    0.142    0.507)    0.582
======================= Grid point 643 (114/144) =======================
q-point: (-0.36  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.139   (  -5.625   -3.248    0.192)    6.499
   1.168   (  -3.080   -1.778    0.518)    3.594
   1.436   (  -1.296   -0.748   -0.640)    1.628
   1.899   (  -0.830   -0.479   -1.604)    1.868
   2.412   (  -4.849   -2.800    0.287)    5.607
   3.057   (   1.511    0.873   -0.261)    1.765
   3.084   (  -2.575   -1.487    0.021)    2.974
   3.673   (  -5.870   -3.389    0.578)    6.802
   3.964   (  -1.249   -0.721    0.043)    1.443
   4.060   (  -1.826   -1.054    0.108)    2.111
   5.007   (  -2.387   -1.378   -0.365)    2.780
   5.230   (   0.566    0.327    0.639)    0.914
======================= Grid point 644 (115/144) =======================
q-point: (-0.29  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.010   (  -5.140   -2.968    1.316)    6.080
   1.071   (  -5.229   -3.019    0.574)    6.065
   1.386   (  -3.015   -1.741   -0.417)    3.506
   1.829   (  -5.557   -3.208   -1.835)    6.674
   2.259   (  -8.354   -4.823    0.347)    9.652
   3.014   (  -3.347   -1.932   -0.039)    3.865
   3.090   (   1.132    0.654   -0.479)    1.392
   3.523   (  -6.693   -3.864    0.332)    7.736
   3.931   (  -1.439   -0.831    0.011)    1.661
   4.005   (  -2.778   -1.604    0.284)    3.220
   4.941   (  -3.184   -1.838   -0.465)    3.706
   5.248   (   0.957    0.553    0.791)    1.359
======================= Grid point 645 (116/144) =======================
q-point: (-0.21  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.908   (  -3.398   -1.962    2.892)    4.874
   0.921   (  -7.479   -4.318    0.901)    8.683
   1.295   (  -4.557   -2.631   -0.229)    5.267
   1.621   ( -12.334   -7.121   -1.604)   14.333
   2.011   ( -12.834   -7.409    0.292)   14.822
   2.932   (  -3.541   -2.044   -0.051)    4.089
   3.099   (  -0.814   -0.470   -0.796)    1.231
   3.375   (  -5.376   -3.104   -0.312)    6.215
   3.898   (  -1.352   -0.781   -0.100)    1.565
   3.938   (  -2.695   -1.556    0.403)    3.137
   4.863   (  -3.250   -1.877   -0.531)    3.790
   5.276   (   1.398    0.807    0.904)    1.850
======================= Grid point 646 (117/144) =======================
q-point: (-0.14  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 372
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.725   (  -8.996   -5.194    1.496)   10.494
   0.854   (  -1.223   -0.706    4.676)    4.885
   1.182   (  -4.826   -2.786   -0.228)    5.578
   1.251   ( -18.779  -10.842   -1.192)   21.717
   1.663   ( -16.296   -9.409    0.192)   18.818
   2.853   (  -2.990   -1.726   -0.022)    3.453
   3.024   (  -5.706   -3.294   -1.450)    6.746
   3.297   (  -1.120   -0.646   -0.619)    1.434
   3.869   (  -1.078   -0.622   -0.191)    1.259
   3.883   (  -1.916   -1.106    0.357)    2.241
   4.795   (  -2.435   -1.406   -0.537)    2.863
   5.313   (   1.770    1.022    0.894)    2.231
======================= Grid point 647 (118/144) =======================
q-point: (-0.07  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 371
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.520   (  -7.679   -4.434    2.220)    9.141
   0.743   ( -18.527  -10.697    3.130)   21.621
   0.881   (  -0.923   -0.533    1.969)    2.239
   1.084   (  -3.140   -1.813   -0.437)    3.652
   1.281   ( -15.035   -8.680    0.059)   17.361
   2.797   (  -1.685   -0.973   -0.040)    1.947
   2.866   (  -6.514   -3.761   -1.102)    7.602
   3.298   (   0.582    0.336   -0.608)    0.906
   3.849   (  -0.659   -0.380   -0.141)    0.774
   3.849   (  -0.951   -0.549    0.226)    1.121
   4.753   (  -1.156   -0.667   -0.541)    1.440
   5.354   (   1.496    0.864    0.589)    1.825
======================= Grid point 649 (119/144) =======================
q-point: ( 0.07  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.700   (   6.507   10.813    1.457)   12.703
   0.847   (   0.556    0.675    5.115)    5.189
   1.077   (   8.040   14.174   -0.994)   16.326
   1.207   (   6.300   10.839   -0.609)   12.551
   1.543   (   8.959   15.312    0.091)   17.741
   2.846   (   2.145    3.467   -0.140)    4.079
   2.976   (   3.292    6.155   -1.304)    7.101
   3.292   (   0.070    0.163   -0.593)    0.619
   3.863   (   0.477    1.444    0.179)    1.531
   3.874   (   1.214    1.252   -0.173)    1.752
   4.780   (   1.199    2.101   -0.521)    2.474
   5.324   (  -1.123   -1.811    0.924)    2.323
======================= Grid point 650 (120/144) =======================
q-point: ( 0.14  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.864   (   2.703    2.106    2.611)    4.308
   0.934   (   3.580   11.959    1.859)   12.621
   1.285   (   5.147    6.806   -1.040)    8.596
   1.498   (   9.677   12.398   -0.552)   15.737
   1.862   (  12.139    9.287    0.005)   15.284
   2.931   (   1.934    5.523   -0.165)    5.854
   3.068   (   2.330    1.037   -0.613)    2.622
   3.327   (   2.918    2.489   -0.438)    3.860
   3.892   (   0.746    2.034    0.193)    2.175
   3.918   (   2.335    1.389   -0.433)    2.751
   4.836   (   2.297    2.637   -0.502)    3.533
   5.285   (  -1.233   -1.409    1.111)    2.177
======================= Grid point 651 (121/144) =======================
q-point: ( 0.21  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.945   (   4.434    2.234    1.451)    5.173
   1.103   (   0.972   12.413    1.400)   12.530
   1.420   (   1.462    8.175   -0.964)    8.360
   1.741   (   7.225    5.359   -0.839)    9.034
   2.139   (  10.642    4.969    0.070)   11.745
   3.010   (   0.586    4.471   -0.026)    4.509
   3.089   (   0.590    0.048   -0.048)    0.594
   3.441   (   6.162    3.108    0.029)    6.901
   3.927   (   0.830    1.998   -0.030)    2.164
   3.981   (   2.838    1.778   -0.467)    3.381
   4.910   (   2.786    2.618   -0.456)    3.850
   5.254   (  -0.900   -1.056    1.028)    1.727
======================= Grid point 652 (122/144) =======================
q-point: ( 0.29  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.062   (   5.730    2.744    0.188)    6.356
   1.219   (  -1.180   10.920    1.368)   11.068
   1.503   (  -1.691    8.595   -0.861)    8.802
   1.863   (   3.315    0.881   -1.187)    3.630
   2.335   (   7.152    2.905    0.278)    7.724
   3.021   (   0.068   -2.152    0.133)    2.157
   3.124   (  -0.148    5.357    0.089)    5.360
   3.588   (   6.942    2.388    0.423)    7.353
   3.959   (   0.744    1.241   -0.228)    1.465
   4.048   (   2.225    2.091   -0.361)    3.075
   4.982   (   2.488    2.134   -0.391)    3.301
   5.232   (  -0.512   -0.858    0.838)    1.304
======================= Grid point 653 (123/144) =======================
q-point: ( 0.36  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.191   (   5.338    2.663   -0.634)    5.999
   1.287   (  -2.647    9.119    1.267)    9.579
   1.540   (  -3.741    8.532   -0.808)    9.351
   1.888   (   0.301   -1.439   -1.421)    2.045
   2.456   (   3.856    2.074    0.463)    4.403
   2.996   (   0.550   -3.466    0.055)    3.509
   3.172   (  -1.020    5.595    0.144)    5.689
   3.722   (   5.775    1.542    0.593)    6.007
   3.975   (   0.534   -0.308   -0.282)    0.678
   4.097   (   0.720    2.324   -0.216)    2.442
   5.036   (   1.585    1.441   -0.328)    2.167
   5.218   (  -0.122   -0.730    0.648)    0.984
======================= Grid point 654 (124/144) =======================
q-point: ( 0.43  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.287   (   2.612    1.446   -0.968)    3.139
   1.313   (  -3.872    7.586    0.985)    8.574
   1.551   (  -4.652    8.384   -0.732)    9.616
   1.867   (  -0.125   -1.784   -1.386)    2.262
   2.517   (   1.117    1.239    0.370)    1.708
   2.976   (   1.258   -3.308    0.026)    3.540
   3.199   (  -1.995    4.831    0.158)    5.229
   3.811   (   2.699    0.532    0.576)    2.810
   3.972   (   0.712   -2.063   -0.182)    2.190
   4.120   (  -0.803    2.461   -0.121)    2.592
   5.065   (   0.397    0.712   -0.267)    0.858
   5.212   (   0.220   -0.621    0.499)    0.826
======================= Grid point 655 (125/144) =======================
q-point: (-0.50  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.279   (  -3.178   -0.628   -0.044)    3.240
   1.325   (  -3.796    6.457   -0.237)    7.494
   1.549   (  -4.996    7.949   -0.550)    9.405
   1.866   (   1.509   -0.731   -1.306)    2.125
   2.518   (  -1.483   -0.311    0.200)    1.528
   2.976   (   2.243   -2.726    0.037)    3.531
   3.199   (  -3.159    4.148    0.134)    5.216
   3.813   (  -1.787   -1.966    0.471)    2.698
   3.970   (   1.310   -1.816    0.030)    2.239
   4.120   (  -1.784    1.936   -0.087)    2.634
   5.065   (  -0.831    0.004   -0.244)    0.867
   5.213   (   0.467   -0.484    0.453)    0.810
======================= Grid point 656 (126/144) =======================
q-point: (-0.43  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.184   (  -4.974   -3.140    0.023)    5.883
   1.301   (  -6.553    6.906   -0.057)    9.521
   1.535   (  -5.391    7.060   -0.385)    8.891
   1.888   (   1.175   -0.731   -1.455)    2.008
   2.458   (  -3.768   -2.422    0.248)    4.486
   2.997   (   2.850   -2.132   -0.102)    3.561
   3.173   (  -4.291    3.678    0.062)    5.651
   3.723   (  -4.320   -4.199    0.502)    6.045
   3.971   (   0.028   -0.698    0.114)    0.708
   4.095   (  -2.533    0.577   -0.033)    2.598
   5.036   (  -2.032   -0.644   -0.298)    2.152
   5.219   (   0.716   -0.265    0.545)    0.938
======================= Grid point 657 (127/144) =======================
q-point: (-0.36  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.056   (  -5.224   -3.555    0.714)    6.359
   1.233   (  -8.986    6.537    0.113)   11.113
   1.495   (  -6.477    5.318   -0.236)    8.384
   1.869   (  -2.225   -2.174   -1.749)    3.569
   2.334   (  -6.205   -4.851    0.356)    7.885
   3.025   (   1.828   -1.131   -0.222)    2.161
   3.126   (  -4.448    2.763   -0.092)    5.237
   3.587   (  -5.602   -4.815    0.476)    7.402
   3.955   (  -1.327    0.016    0.066)    1.329
   4.043   (  -3.167   -0.809    0.110)    3.270
   4.982   (  -3.074   -1.080   -0.396)    3.282
   5.233   (   1.030   -0.006    0.707)    1.249
======================= Grid point 658 (128/144) =======================
q-point: (-0.29  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.939   (  -4.220   -2.646    1.928)    5.341
   1.114   ( -11.404    5.250    0.452)   12.562
   1.411   (  -7.738    2.632   -0.178)    8.175
   1.750   (  -8.209   -3.603   -1.736)    9.132
   2.136   (  -9.679   -6.710    0.344)   11.782
   3.011   (  -4.317    1.650   -0.096)    4.622
   3.093   (  -0.213   -0.605   -0.467)    0.794
   3.440   (  -5.724   -3.756    0.076)    6.846
   3.928   (  -1.981    0.440   -0.076)    2.031
   3.973   (  -3.149   -1.448    0.288)    3.478
   4.910   (  -3.650   -1.112   -0.485)    3.846
   5.256   (   1.405    0.198    0.865)    1.662
======================= Grid point 659 (129/144) =======================
q-point: (-0.21  0.07  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.857   (  -3.565   -0.804    3.116)    4.802
   0.941   ( -12.049    2.414    1.359)   12.363
   1.278   (  -8.403   -0.618   -0.432)    8.437
   1.505   ( -15.609   -2.692   -1.183)   15.883
   1.860   ( -14.135   -5.686    0.260)   15.238
   2.930   (  -5.791    1.108   -0.070)    5.897
   3.071   (  -2.044   -1.736   -0.937)    2.840
   3.327   (  -3.591   -1.316   -0.476)    3.854
   3.896   (  -2.122    0.497   -0.209)    2.190
   3.910   (  -2.323   -1.104    0.323)    2.592
   4.837   (  -3.431   -0.693   -0.526)    3.540
   5.287   (   1.849    0.245    0.979)    2.107
======================= Grid point 664 (130/144) =======================
q-point: ( 0.14  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.914   (   2.213    3.473    1.898)    4.535
   1.192   (   6.234   11.715    0.973)   13.306
   1.436   (   4.886    7.442   -0.825)    8.941
   1.716   (   4.474    8.578   -1.082)    9.735
   2.047   (   5.395    9.116    0.250)   10.596
   3.059   (   4.262    3.047   -0.076)    5.240
   3.061   (  -1.530    1.970   -0.327)    2.516
   3.396   (   2.471    4.258   -0.052)    4.924
   3.938   (   0.667    1.684    0.007)    1.811
   3.955   (   2.235    2.916   -0.130)    3.676
   4.898   (   1.895    3.304   -0.471)    3.838
   5.257   (  -0.832   -1.340    0.963)    1.848
======================= Grid point 665 (131/144) =======================
q-point: ( 0.21  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.009   (   4.058    4.127    0.523)    5.811
   1.371   (  -0.163   12.848    1.096)   12.896
   1.586   (   3.422    7.018   -1.188)    7.898
   1.833   (   1.781    3.856   -0.895)    4.340
   2.235   (   7.395    4.987    0.231)    8.922
   3.013   (  -1.574   -2.682    0.016)    3.110
   3.178   (   1.647    6.876    0.062)    7.070
   3.499   (   5.094    2.805    0.301)    5.823
   3.962   (   0.377    1.001   -0.127)    1.077
   4.031   (   2.573    3.286   -0.213)    4.179
   4.966   (   2.107    2.811   -0.387)    3.534
   5.231   (  -0.631   -1.136    0.824)    1.539
======================= Grid point 666 (132/144) =======================
q-point: ( 0.29  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.134   (   4.919    4.350   -0.452)    6.582
   1.459   (  -3.620   11.542    1.034)   12.140
   1.689   (  -0.632    8.200   -1.299)    8.326
   1.866   (   0.039   -0.022   -0.865)    0.866
   2.396   (   5.778    3.492    0.340)    6.760
   2.960   (  -0.688   -3.497    0.106)    3.566
   3.256   (  -0.792    7.022    0.157)    7.068
   3.625   (   6.407    1.507    0.505)    6.601
   3.968   (  -0.070   -0.432   -0.228)    0.494
   4.105   (   1.705    3.389   -0.189)    3.798
   5.025   (   1.698    2.116   -0.306)    2.731
   5.212   (  -0.285   -1.070    0.643)    1.281
======================= Grid point 667 (133/144) =======================
q-point: ( 0.36  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.257   (   3.737    3.874   -0.952)    5.466
   1.488   (  -4.918    9.417    0.711)   10.648
   1.736   (  -4.066    9.548   -1.163)   10.443
   1.846   (  -0.460   -2.215   -0.861)    2.421
   2.504   (   2.668    2.873    0.383)    3.940
   2.922   (   0.569   -3.728    0.070)    3.772
   3.297   (  -2.298    6.593    0.203)    6.985
   3.745   (   5.419    0.815    0.585)    5.511
   3.949   (  -0.546   -2.196   -0.211)    2.272
   4.155   (   0.189    3.116   -0.164)    3.126
   5.066   (   0.771    1.418   -0.235)    1.631
   5.200   (   0.145   -1.049    0.479)    1.163
======================= Grid point 668 (134/144) =======================
q-point: ( 0.43  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.316   (  -1.043    1.537   -0.942)    2.083
   1.492   (  -4.506    8.525    0.081)    9.643
   1.743   (  -5.929    9.174   -0.393)   10.930
   1.830   (   1.111   -1.290   -1.211)    2.089
   2.544   (  -0.807    1.503    0.279)    1.729
   2.907   (   1.925   -3.340    0.085)    3.856
   3.310   (  -3.449    5.981    0.214)    6.907
   3.809   (   0.738   -0.947    0.470)    1.289
   3.921   (   1.229   -2.545   -0.015)    2.826
   4.173   (  -1.415    2.459   -0.165)    2.842
   5.080   (  -0.431    0.725   -0.185)    0.864
   5.196   (   0.566   -0.941    0.381)    1.163
======================= Grid point 670 (135/144) =======================
q-point: (-0.50  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.253   (  -5.349   -1.195   -0.625)    5.516
   1.500   (  -5.375    9.096   -0.330)   10.571
   1.721   (  -5.833    6.937    0.189)    9.065
   1.854   (   1.619    0.419   -1.533)    2.269
   2.506   (  -3.791   -0.951    0.269)    3.918
   2.922   (   3.012   -2.326    0.022)    3.806
   3.297   (  -4.535    5.276    0.182)    6.959
   3.747   (  -3.470   -4.271    0.330)    5.513
   3.944   (   2.123   -0.623    0.202)    2.222
   4.155   (  -2.824    1.407   -0.157)    3.159
   5.065   (  -1.617    0.044   -0.207)    1.630
   5.200   (   0.858   -0.656    0.420)    1.159
======================= Grid point 671 (136/144) =======================
q-point: (-0.43  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.130   (  -6.287   -1.919   -0.121)    6.575
   1.473   (  -8.201    8.849   -0.227)   12.067
   1.672   (  -6.336    4.121    0.166)    7.561
   1.875   (  -0.198    0.314   -1.634)    1.676
   2.396   (  -5.977   -3.344    0.364)    6.859
   2.962   (   3.468   -1.199   -0.101)    3.671
   3.257   (  -5.649    4.173    0.097)    7.024
   3.626   (  -4.609   -4.826    0.417)    6.686
   3.964   (   0.438   -0.032    0.133)    0.458
   4.104   (  -3.892    0.267   -0.096)    3.902
   5.025   (  -2.671   -0.405   -0.294)    2.718
   5.213   (   1.083   -0.317    0.575)    1.267
======================= Grid point 672 (137/144) =======================
q-point: (-0.36  0.14  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.007   (  -5.644   -1.221    0.745)    5.822
   1.381   ( -11.173    5.976    0.274)   12.673
   1.576   (  -7.557    1.073   -0.267)    7.637
   1.840   (  -4.214   -0.063   -1.538)    4.487
   2.233   (  -8.058   -3.832    0.373)    8.931
   3.015   (   3.118   -0.156   -0.173)    3.127
   3.179   (  -6.754    1.943   -0.008)    7.028
   3.499   (  -5.001   -2.984    0.334)    5.833
   3.960   (  -0.950    0.300   -0.005)    0.997
   4.029   (  -4.287   -0.385   -0.006)    4.304
   4.966   (  -3.481   -0.426   -0.394)    3.529
   5.232   (   1.317   -0.076    0.768)    1.526
======================= Grid point 680 (138/144) =======================
q-point: ( 0.21  0.21  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.111   (   3.437    5.772   -0.270)    6.724
   1.593   (   3.509    8.325    0.796)    9.069
   1.683   (   3.109    2.675   -1.289)    4.299
   1.896   (   0.746    2.105   -0.946)    2.426
   2.353   (   3.813    6.523    0.364)    7.564
   2.949   (  -1.969   -3.262    0.017)    3.811
   3.311   (   3.228    5.628    0.154)    6.490
   3.564   (   2.171    3.814    0.407)    4.407
   3.966   (  -0.186   -0.550   -0.039)    0.582
   4.106   (   2.202    3.767   -0.143)    4.366
   5.021   (   1.468    2.538   -0.290)    2.946
   5.209   (  -0.625   -1.045    0.610)    1.362
======================= Grid point 681 (139/144) =======================
q-point: ( 0.29  0.21  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.237   (   3.941    5.725   -0.922)    7.011
   1.649   (  -4.039    5.795    0.245)    7.068
   1.769   (   3.943    0.449   -1.233)    4.156
   1.912   (  -2.283    4.237   -0.502)    4.839
   2.485   (   3.858    5.148    0.374)    6.444
   2.887   (  -1.018   -3.543    0.093)    3.688
   3.398   (  -0.462    6.611    0.243)    6.631
   3.657   (   4.944    1.792    0.464)    5.279
   3.942   (  -1.014   -2.113   -0.076)    2.345
   4.175   (   1.304    3.211   -0.184)    3.471
   5.067   (   1.065    1.903   -0.197)    2.190
   5.189   (  -0.307   -1.124    0.439)    1.245
======================= Grid point 682 (140/144) =======================
q-point: ( 0.36  0.21  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.347   (   1.379    4.871   -1.199)    5.202
   1.615   (  -3.000    1.662   -0.677)    3.496
   1.811   (   0.315    0.776   -0.106)    0.844
   1.913   (  -3.485    6.388   -0.805)    7.321
   2.576   (   0.620    4.117    0.338)    4.177
   2.845   (   0.719   -3.746    0.132)    3.816
   3.434   (  -2.875    6.630    0.281)    7.232
   3.759   (   4.453    0.607    0.481)    4.520
   3.890   (  -1.210   -3.442   -0.050)    3.648
   4.213   (  -0.233    2.386   -0.208)    2.407
   5.094   (   0.069    1.372   -0.120)    1.379
   5.177   (   0.275   -1.202    0.305)    1.270
======================= Grid point 683 (141/144) =======================
q-point: (-0.57  0.21  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.346   (  -4.994    1.255   -1.102)    5.265
   1.619   (   0.576    3.423   -1.138)    3.653
   1.803   (  -1.893    0.395    0.679)    2.050
   1.918   (  -3.158    5.927   -1.220)    6.825
   2.576   (  -3.852    1.507    0.316)    4.149
   2.845   (   2.896   -2.492    0.130)    3.823
   3.434   (  -4.307    5.788    0.295)    7.221
   3.762   (  -2.663   -3.561    0.248)    4.453
   3.887   (   3.460   -0.643    0.198)    3.525
   4.213   (  -1.942    1.396   -0.217)    2.402
   5.094   (  -1.233    0.625   -0.105)    1.386
   5.177   (   0.922   -0.839    0.283)    1.278
======================= Grid point 684 (142/144) =======================
q-point: (-0.50  0.21  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.236   (  -6.961   -0.431   -0.807)    7.021
   1.663   (  -2.558    6.469   -1.093)    7.042
   1.749   (  -2.946   -2.625    0.613)    3.993
   1.919   (  -2.309    3.320   -1.152)    4.205
   2.485   (  -6.406   -0.786    0.378)    6.465
   2.888   (   3.620   -0.929    0.045)    3.737
   3.397   (  -5.472    3.699    0.252)    6.610
   3.659   (  -4.068   -3.511    0.320)    5.383
   3.940   (   2.322    0.020    0.125)    2.326
   4.175   (  -3.444    0.474   -0.187)    3.481
   5.066   (  -2.178    0.041   -0.185)    2.186
   5.190   (   1.132   -0.322    0.412)    1.246
======================= Grid point 695 (143/144) =======================
q-point: ( 0.29  0.29  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 688
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.360   (   3.544    6.067   -1.230)    7.133
   1.645   (  -2.528   -4.056   -1.495)    5.008
   1.816   (   2.006    2.984    0.797)    3.683
   1.980   (   1.085    2.113   -0.995)    2.575
   2.595   (   2.986    5.180    0.363)    5.990
   2.817   (  -1.920   -3.320    0.152)    3.838
   3.510   (   2.305    3.949    0.335)    4.585
   3.707   (   1.754    3.241    0.305)    3.698
   3.884   (  -1.919   -3.505    0.106)    3.997
   4.228   (   1.080    1.878   -0.230)    2.178
   5.101   (   0.865    1.486   -0.078)    1.721
   5.167   (  -0.618   -1.052    0.255)    1.246
======================= Grid point 696 (144/144) =======================
q-point: ( 0.36  0.29  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 
Number of triplets: 686
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.448   (  -2.682    4.602   -1.380)    5.503
   1.575   (   2.204   -3.737   -1.482)    4.585
   1.848   (  -0.933    1.520    0.886)    1.992
   2.000   (  -1.166    2.077   -1.120)    2.632
   2.664   (  -2.320    4.025    0.316)    4.657
   2.769   (   2.057   -3.566    0.223)    4.123
   3.549   (  -2.248    3.874    0.365)    4.493
   3.768   (   0.028    0.155    0.269)    0.312
   3.822   (   1.363   -2.551    0.140)    2.895
   4.247   (  -0.505    0.878   -0.248)    1.043
   5.121   (  -0.681    1.172    0.001)    1.356
   5.152   (   0.653   -1.121    0.154)    1.306
=================== 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/16464
   10.0    113.866    113.866     24.937     -0.000      0.000     -0.000 3/16464
   20.0     38.270     38.270      6.123     -0.000      0.000     -0.000 3/16464
   30.0     19.737     19.737      3.071     -0.000      0.000     -0.000 3/16464
   40.0     12.502     12.502      2.013     -0.000      0.000     -0.000 3/16464
   50.0      9.009      9.009      1.496     -0.000      0.000     -0.000 3/16464
   60.0      7.027      7.027      1.193     -0.000      0.000     -0.000 3/16464
   70.0      5.766      5.766      0.993     -0.000      0.000     -0.000 3/16464
   80.0      4.896      4.896      0.852     -0.000      0.000     -0.000 3/16464
   90.0      4.260      4.260      0.746     -0.000      0.000     -0.000 3/16464
  100.0      3.775      3.775      0.665     -0.000      0.000     -0.000 3/16464
  110.0      3.391      3.391      0.600     -0.000      0.000     -0.000 3/16464
  120.0      3.080      3.080      0.546     -0.000      0.000     -0.000 3/16464
  130.0      2.823      2.823      0.502     -0.000      0.000     -0.000 3/16464
  140.0      2.607      2.607      0.464     -0.000      0.000     -0.000 3/16464
  150.0      2.422      2.422      0.432     -0.000      0.000     -0.000 3/16464
  160.0      2.262      2.262      0.404     -0.000      0.000     -0.000 3/16464
  170.0      2.122      2.122      0.379     -0.000      0.000     -0.000 3/16464
  180.0      1.999      1.999      0.358     -0.000      0.000     -0.000 3/16464
  190.0      1.889      1.889      0.338     -0.000      0.000     -0.000 3/16464
  200.0      1.791      1.791      0.321     -0.000      0.000     -0.000 3/16464
  210.0      1.703      1.703      0.305     -0.000      0.000     -0.000 3/16464
  220.0      1.623      1.623      0.291     -0.000      0.000     -0.000 3/16464
  230.0      1.551      1.551      0.278     -0.000      0.000     -0.000 3/16464
  240.0      1.485      1.485      0.266     -0.000      0.000     -0.000 3/16464
  250.0      1.424      1.424      0.256     -0.000      0.000     -0.000 3/16464
  260.0      1.368      1.368      0.246     -0.000      0.000     -0.000 3/16464
  270.0      1.316      1.316      0.236     -0.000      0.000     -0.000 3/16464
  280.0      1.268      1.268      0.228     -0.000      0.000     -0.000 3/16464
  290.0      1.224      1.224      0.220     -0.000      0.000     -0.000 3/16464
  300.0      1.182      1.182      0.212     -0.000      0.000     -0.000 3/16464
  310.0      1.144      1.144      0.206     -0.000      0.000     -0.000 3/16464
  320.0      1.107      1.107      0.199     -0.000      0.000     -0.000 3/16464
  330.0      1.073      1.073      0.193     -0.000      0.000     -0.000 3/16464
  340.0      1.041      1.041      0.187     -0.000      0.000     -0.000 3/16464
  350.0      1.011      1.011      0.182     -0.000      0.000     -0.000 3/16464
  360.0      0.983      0.983      0.177     -0.000      0.000     -0.000 3/16464
  370.0      0.956      0.956      0.172     -0.000      0.000     -0.000 3/16464
  380.0      0.931      0.931      0.167     -0.000      0.000     -0.000 3/16464
  390.0      0.907      0.907      0.163     -0.000      0.000     -0.000 3/16464
  400.0      0.884      0.884      0.159     -0.000      0.000     -0.000 3/16464
  410.0      0.862      0.862      0.155     -0.000      0.000     -0.000 3/16464
  420.0      0.841      0.841      0.151     -0.000      0.000     -0.000 3/16464
  430.0      0.822      0.822      0.148     -0.000      0.000     -0.000 3/16464
  440.0      0.803      0.803      0.144     -0.000      0.000     -0.000 3/16464
  450.0      0.785      0.785      0.141     -0.000      0.000     -0.000 3/16464
  460.0      0.768      0.768      0.138     -0.000      0.000     -0.000 3/16464
  470.0      0.751      0.751      0.135     -0.000      0.000     -0.000 3/16464
  480.0      0.735      0.735      0.132     -0.000      0.000     -0.000 3/16464
  490.0      0.720      0.720      0.130     -0.000      0.000     -0.000 3/16464
  500.0      0.706      0.706      0.127     -0.000      0.000     -0.000 3/16464
  510.0      0.692      0.692      0.125     -0.000      0.000     -0.000 3/16464
  520.0      0.679      0.679      0.122     -0.000      0.000     -0.000 3/16464
  530.0      0.666      0.666      0.120     -0.000      0.000     -0.000 3/16464
  540.0      0.653      0.653      0.118     -0.000      0.000     -0.000 3/16464
  550.0      0.641      0.641      0.115     -0.000      0.000     -0.000 3/16464
  560.0      0.630      0.630      0.113     -0.000      0.000     -0.000 3/16464
  570.0      0.619      0.619      0.111     -0.000      0.000     -0.000 3/16464
  580.0      0.608      0.608      0.109     -0.000      0.000     -0.000 3/16464
  590.0      0.598      0.598      0.108     -0.000      0.000     -0.000 3/16464
  600.0      0.588      0.588      0.106     -0.000      0.000     -0.000 3/16464
  610.0      0.578      0.578      0.104     -0.000      0.000     -0.000 3/16464
  620.0      0.569      0.569      0.102     -0.000      0.000     -0.000 3/16464
  630.0      0.560      0.560      0.101     -0.000      0.000     -0.000 3/16464
  640.0      0.551      0.551      0.099     -0.000      0.000     -0.000 3/16464
  650.0      0.542      0.542      0.098     -0.000      0.000     -0.000 3/16464
  660.0      0.534      0.534      0.096     -0.000      0.000     -0.000 3/16464
  670.0      0.526      0.526      0.095     -0.000      0.000     -0.000 3/16464
  680.0      0.518      0.518      0.093     -0.000      0.000     -0.000 3/16464
  690.0      0.511      0.511      0.092     -0.000      0.000     -0.000 3/16464
  700.0      0.504      0.504      0.091     -0.000      0.000     -0.000 3/16464
  710.0      0.496      0.496      0.089     -0.000      0.000     -0.000 3/16464
  720.0      0.489      0.489      0.088     -0.000      0.000     -0.000 3/16464
  730.0      0.483      0.483      0.087     -0.000      0.000     -0.000 3/16464
  740.0      0.476      0.476      0.086     -0.000      0.000     -0.000 3/16464
  750.0      0.470      0.470      0.085     -0.000      0.000     -0.000 3/16464
  760.0      0.464      0.464      0.083     -0.000      0.000     -0.000 3/16464
  770.0      0.458      0.458      0.082     -0.000      0.000     -0.000 3/16464
  780.0      0.452      0.452      0.081     -0.000      0.000     -0.000 3/16464
  790.0      0.446      0.446      0.080     -0.000      0.000     -0.000 3/16464
  800.0      0.440      0.440      0.079     -0.000      0.000     -0.000 3/16464
  810.0      0.435      0.435      0.078     -0.000      0.000     -0.000 3/16464
  820.0      0.430      0.430      0.077     -0.000      0.000     -0.000 3/16464
  830.0      0.424      0.424      0.076     -0.000      0.000     -0.000 3/16464
  840.0      0.419      0.419      0.076     -0.000      0.000     -0.000 3/16464
  850.0      0.414      0.414      0.075     -0.000      0.000     -0.000 3/16464
  860.0      0.410      0.410      0.074     -0.000      0.000     -0.000 3/16464
  870.0      0.405      0.405      0.073     -0.000      0.000     -0.000 3/16464
  880.0      0.400      0.400      0.072     -0.000      0.000     -0.000 3/16464
  890.0      0.396      0.396      0.071     -0.000      0.000     -0.000 3/16464
  900.0      0.391      0.391      0.070     -0.000      0.000     -0.000 3/16464
  910.0      0.387      0.387      0.070     -0.000      0.000     -0.000 3/16464
  920.0      0.383      0.383      0.069     -0.000      0.000     -0.000 3/16464
  930.0      0.379      0.379      0.068     -0.000      0.000     -0.000 3/16464
  940.0      0.375      0.375      0.067     -0.000      0.000     -0.000 3/16464
  950.0      0.371      0.371      0.067     -0.000      0.000     -0.000 3/16464
  960.0      0.367      0.367      0.066     -0.000      0.000     -0.000 3/16464
  970.0      0.363      0.363      0.065     -0.000      0.000     -0.000 3/16464
  980.0      0.359      0.359      0.065     -0.000      0.000     -0.000 3/16464
  990.0      0.356      0.356      0.064     -0.000      0.000     -0.000 3/16464
 1000.0      0.352      0.352      0.063     -0.000      0.000     -0.000 3/16464

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

