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

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

-------------------------[time 2026-01-09 04:38:03]-------------------------
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: P6_3mc (186)
Number of symmetry operations in supercell: 384
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.134469640253094    0.000000000000000    0.000000000000000
  b   -2.067234820126547    3.580555739634689    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.609144710000000
Atomic positions (fractional):
   *1 Li  0.33333333333333  0.66666666666667  0.37807293312149   6.941
    2 Li  0.66666666666667  0.33333333333333  0.87807293312149   6.941
   *3 Br  0.33333333333333  0.66666666666667  0.99692506687852  79.904
    4 Br  0.66666666666667  0.33333333333333  0.49692506687852  79.904
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.134469640253094    0.000000000000000    0.000000000000000
  b   -2.067234820126547    3.580555739634689    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.609144710000000
Atomic positions (fractional):
   *1 Li  0.33333333333333  0.66666666666667  0.37807293312149   6.941 > 1
    2 Li  0.66666666666667  0.33333333333333  0.87807293312149   6.941 > 2
   *3 Br  0.33333333333333  0.66666666666667  0.99692506687852  79.904 > 3
    4 Br  0.66666666666667  0.33333333333333  0.49692506687852  79.904 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   16.537878561012377    0.000000000000000    0.000000000000000
  b   -8.268939280506189   14.322222958538756    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.218289420000000
Atomic positions (fractional):
   *1 Li  0.08333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    3 Li  0.58333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    4 Li  0.83333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    5 Li  0.08333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    6 Li  0.33333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    7 Li  0.58333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    8 Li  0.83333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    9 Li  0.08333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   10 Li  0.33333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   11 Li  0.58333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   12 Li  0.83333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   13 Li  0.08333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   14 Li  0.33333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   15 Li  0.58333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   16 Li  0.83333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   17 Li  0.08333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   18 Li  0.33333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   19 Li  0.58333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   20 Li  0.83333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   21 Li  0.08333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   22 Li  0.33333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   23 Li  0.58333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   24 Li  0.83333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   25 Li  0.08333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   26 Li  0.33333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   27 Li  0.58333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   28 Li  0.83333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   29 Li  0.08333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   30 Li  0.33333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   31 Li  0.58333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   32 Li  0.83333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   33 Li  0.16666666666667  0.08333333333333  0.43903646656074   6.941 > 2
   34 Li  0.41666666666667  0.08333333333333  0.43903646656074   6.941 > 2
   35 Li  0.66666666666667  0.08333333333333  0.43903646656074   6.941 > 2
   36 Li  0.91666666666667  0.08333333333333  0.43903646656074   6.941 > 2
   37 Li  0.16666666666667  0.33333333333333  0.43903646656074   6.941 > 2
   38 Li  0.41666666666667  0.33333333333333  0.43903646656074   6.941 > 2
   39 Li  0.66666666666667  0.33333333333333  0.43903646656074   6.941 > 2
   40 Li  0.91666666666667  0.33333333333333  0.43903646656074   6.941 > 2
   41 Li  0.16666666666667  0.58333333333333  0.43903646656074   6.941 > 2
   42 Li  0.41666666666667  0.58333333333333  0.43903646656074   6.941 > 2
   43 Li  0.66666666666667  0.58333333333333  0.43903646656074   6.941 > 2
   44 Li  0.91666666666667  0.58333333333333  0.43903646656074   6.941 > 2
   45 Li  0.16666666666667  0.83333333333333  0.43903646656074   6.941 > 2
   46 Li  0.41666666666667  0.83333333333333  0.43903646656074   6.941 > 2
   47 Li  0.66666666666667  0.83333333333333  0.43903646656074   6.941 > 2
   48 Li  0.91666666666667  0.83333333333333  0.43903646656074   6.941 > 2
   49 Li  0.16666666666667  0.08333333333333  0.93903646656074   6.941 > 2
   50 Li  0.41666666666667  0.08333333333333  0.93903646656074   6.941 > 2
   51 Li  0.66666666666667  0.08333333333333  0.93903646656074   6.941 > 2
   52 Li  0.91666666666667  0.08333333333333  0.93903646656074   6.941 > 2
   53 Li  0.16666666666667  0.33333333333333  0.93903646656074   6.941 > 2
   54 Li  0.41666666666667  0.33333333333333  0.93903646656074   6.941 > 2
   55 Li  0.66666666666667  0.33333333333333  0.93903646656074   6.941 > 2
   56 Li  0.91666666666667  0.33333333333333  0.93903646656074   6.941 > 2
   57 Li  0.16666666666667  0.58333333333333  0.93903646656074   6.941 > 2
   58 Li  0.41666666666667  0.58333333333333  0.93903646656074   6.941 > 2
   59 Li  0.66666666666667  0.58333333333333  0.93903646656074   6.941 > 2
   60 Li  0.91666666666667  0.58333333333333  0.93903646656074   6.941 > 2
   61 Li  0.16666666666667  0.83333333333333  0.93903646656074   6.941 > 2
   62 Li  0.41666666666667  0.83333333333333  0.93903646656074   6.941 > 2
   63 Li  0.66666666666667  0.83333333333333  0.93903646656074   6.941 > 2
   64 Li  0.91666666666667  0.83333333333333  0.93903646656074   6.941 > 2
  *65 Br  0.08333333333333  0.16666666666667  0.49846253343926  79.904 > 3
   66 Br  0.33333333333333  0.16666666666667  0.49846253343926  79.904 > 3
   67 Br  0.58333333333333  0.16666666666667  0.49846253343926  79.904 > 3
   68 Br  0.83333333333333  0.16666666666667  0.49846253343926  79.904 > 3
   69 Br  0.08333333333333  0.41666666666667  0.49846253343926  79.904 > 3
   70 Br  0.33333333333333  0.41666666666667  0.49846253343926  79.904 > 3
   71 Br  0.58333333333333  0.41666666666667  0.49846253343926  79.904 > 3
   72 Br  0.83333333333333  0.41666666666667  0.49846253343926  79.904 > 3
   73 Br  0.08333333333333  0.66666666666667  0.49846253343926  79.904 > 3
   74 Br  0.33333333333333  0.66666666666667  0.49846253343926  79.904 > 3
   75 Br  0.58333333333333  0.66666666666667  0.49846253343926  79.904 > 3
   76 Br  0.83333333333333  0.66666666666667  0.49846253343926  79.904 > 3
   77 Br  0.08333333333333  0.91666666666667  0.49846253343926  79.904 > 3
   78 Br  0.33333333333333  0.91666666666667  0.49846253343926  79.904 > 3
   79 Br  0.58333333333333  0.91666666666667  0.49846253343926  79.904 > 3
   80 Br  0.83333333333333  0.91666666666667  0.49846253343926  79.904 > 3
   81 Br  0.08333333333333  0.16666666666667  0.99846253343926  79.904 > 3
   82 Br  0.33333333333333  0.16666666666667  0.99846253343926  79.904 > 3
   83 Br  0.58333333333333  0.16666666666667  0.99846253343926  79.904 > 3
   84 Br  0.83333333333333  0.16666666666667  0.99846253343926  79.904 > 3
   85 Br  0.08333333333333  0.41666666666667  0.99846253343926  79.904 > 3
   86 Br  0.33333333333333  0.41666666666667  0.99846253343926  79.904 > 3
   87 Br  0.58333333333333  0.41666666666667  0.99846253343926  79.904 > 3
   88 Br  0.83333333333333  0.41666666666667  0.99846253343926  79.904 > 3
   89 Br  0.08333333333333  0.66666666666667  0.99846253343926  79.904 > 3
   90 Br  0.33333333333333  0.66666666666667  0.99846253343926  79.904 > 3
   91 Br  0.58333333333333  0.66666666666667  0.99846253343926  79.904 > 3
   92 Br  0.83333333333333  0.66666666666667  0.99846253343926  79.904 > 3
   93 Br  0.08333333333333  0.91666666666667  0.99846253343926  79.904 > 3
   94 Br  0.33333333333333  0.91666666666667  0.99846253343926  79.904 > 3
   95 Br  0.58333333333333  0.91666666666667  0.99846253343926  79.904 > 3
   96 Br  0.83333333333333  0.91666666666667  0.99846253343926  79.904 > 3
   97 Br  0.16666666666667  0.08333333333333  0.24846253343926  79.904 > 4
   98 Br  0.41666666666667  0.08333333333333  0.24846253343926  79.904 > 4
   99 Br  0.66666666666667  0.08333333333333  0.24846253343926  79.904 > 4
  100 Br  0.91666666666667  0.08333333333333  0.24846253343926  79.904 > 4
  101 Br  0.16666666666667  0.33333333333333  0.24846253343926  79.904 > 4
  102 Br  0.41666666666667  0.33333333333333  0.24846253343926  79.904 > 4
  103 Br  0.66666666666667  0.33333333333333  0.24846253343926  79.904 > 4
  104 Br  0.91666666666667  0.33333333333333  0.24846253343926  79.904 > 4
  105 Br  0.16666666666667  0.58333333333333  0.24846253343926  79.904 > 4
  106 Br  0.41666666666667  0.58333333333333  0.24846253343926  79.904 > 4
  107 Br  0.66666666666667  0.58333333333333  0.24846253343926  79.904 > 4
  108 Br  0.91666666666667  0.58333333333333  0.24846253343926  79.904 > 4
  109 Br  0.16666666666667  0.83333333333333  0.24846253343926  79.904 > 4
  110 Br  0.41666666666667  0.83333333333333  0.24846253343926  79.904 > 4
  111 Br  0.66666666666667  0.83333333333333  0.24846253343926  79.904 > 4
  112 Br  0.91666666666667  0.83333333333333  0.24846253343926  79.904 > 4
  113 Br  0.16666666666667  0.08333333333333  0.74846253343926  79.904 > 4
  114 Br  0.41666666666667  0.08333333333333  0.74846253343926  79.904 > 4
  115 Br  0.66666666666667  0.08333333333333  0.74846253343926  79.904 > 4
  116 Br  0.91666666666667  0.08333333333333  0.74846253343926  79.904 > 4
  117 Br  0.16666666666667  0.33333333333333  0.74846253343926  79.904 > 4
  118 Br  0.41666666666667  0.33333333333333  0.74846253343926  79.904 > 4
  119 Br  0.66666666666667  0.33333333333333  0.74846253343926  79.904 > 4
  120 Br  0.91666666666667  0.33333333333333  0.74846253343926  79.904 > 4
  121 Br  0.16666666666667  0.58333333333333  0.74846253343926  79.904 > 4
  122 Br  0.41666666666667  0.58333333333333  0.74846253343926  79.904 > 4
  123 Br  0.66666666666667  0.58333333333333  0.74846253343926  79.904 > 4
  124 Br  0.91666666666667  0.58333333333333  0.74846253343926  79.904 > 4
  125 Br  0.16666666666667  0.83333333333333  0.74846253343926  79.904 > 4
  126 Br  0.41666666666667  0.83333333333333  0.74846253343926  79.904 > 4
  127 Br  0.66666666666667  0.83333333333333  0.74846253343926  79.904 > 4
  128 Br  0.91666666666667  0.83333333333333  0.74846253343926  79.904 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.8334873    0.0000000    0.0000000
            0.0000000    2.8334873    0.0000000
            0.0000000    0.0000000    2.8782661
-------------------------- Born effective charges --------------------------
    1 Li    1.0116389    0.0000000    0.0000000
            0.0000000    1.0116389    0.0000000
            0.0000000    0.0000000    1.0353697
    2 Li    1.0116389    0.0000000    0.0000000
            0.0000000    1.0116389    0.0000000
            0.0000000    0.0000000    1.0353697
    3 Br   -1.0116389   -0.0000000    0.0000000
            0.0000000   -1.0116389    0.0000000
            0.0000000    0.0000000   -1.0353697
    4 Br   -1.0116389   -0.0000000    0.0000000
            0.0000000   -1.0116389    0.0000000
            0.0000000    0.0000000   -1.0353697
----------------------------------------------------------------------------
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: 224
Number of blocks in projector: 172
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: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
-----
Solver_atoms: 1 -- 128 / 128
Time (Solver_compr_matrix_reshape): 0.010
Solver_block: 80 / 80
 - Time: 0.159
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.171
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) 
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: 224
Number of blocks in projector: 172
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: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) 
Force constants are written into "force_constants.hdf5".

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

Summary of calculation was written in "phonopy.yaml".
-------------------------[time 2026-01-09 04:38:07]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-09 04:38:08]-------------------------
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: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.134469640253094    0.000000000000000    0.000000000000000
  b   -2.067234820126547    3.580555739634689    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.609144710000000
Atomic positions (fractional):
    1 Li  0.33333333333333  0.66666666666667  0.37807293312149   6.941
    2 Li  0.66666666666667  0.33333333333333  0.87807293312149   6.941
    3 Br  0.33333333333333  0.66666666666667  0.99692506687852  79.904
    4 Br  0.66666666666667  0.33333333333333  0.49692506687852  79.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.537878561012377    0.000000000000000    0.000000000000000
  b   -8.268939280506189   14.322222958538756    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.218289420000000
Atomic positions (fractional):
    1 Li  0.08333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    3 Li  0.58333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    4 Li  0.83333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    5 Li  0.08333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    6 Li  0.33333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    7 Li  0.58333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    8 Li  0.83333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    9 Li  0.08333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   10 Li  0.33333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   11 Li  0.58333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   12 Li  0.83333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   13 Li  0.08333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   14 Li  0.33333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   15 Li  0.58333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   16 Li  0.83333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   17 Li  0.08333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   18 Li  0.33333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   19 Li  0.58333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   20 Li  0.83333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   21 Li  0.08333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   22 Li  0.33333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   23 Li  0.58333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   24 Li  0.83333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   25 Li  0.08333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   26 Li  0.33333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   27 Li  0.58333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   28 Li  0.83333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   29 Li  0.08333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   30 Li  0.33333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   31 Li  0.58333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   32 Li  0.83333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   33 Li  0.16666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   34 Li  0.41666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   35 Li  0.66666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   36 Li  0.91666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   37 Li  0.16666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   38 Li  0.41666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   39 Li  0.66666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   40 Li  0.91666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   41 Li  0.16666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   42 Li  0.41666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   43 Li  0.66666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   44 Li  0.91666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   45 Li  0.16666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   46 Li  0.41666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   47 Li  0.66666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   48 Li  0.91666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   49 Li  0.16666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   50 Li  0.41666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   51 Li  0.66666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   52 Li  0.91666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   53 Li  0.16666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   54 Li  0.41666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   55 Li  0.66666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   56 Li  0.91666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   57 Li  0.16666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   58 Li  0.41666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   59 Li  0.66666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   60 Li  0.91666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   61 Li  0.16666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   62 Li  0.41666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   63 Li  0.66666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   64 Li  0.91666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   65 Br  0.08333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   66 Br  0.33333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   67 Br  0.58333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   68 Br  0.83333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   69 Br  0.08333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   70 Br  0.33333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   71 Br  0.58333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   72 Br  0.83333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   73 Br  0.08333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   74 Br  0.33333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   75 Br  0.58333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   76 Br  0.83333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   77 Br  0.08333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   78 Br  0.33333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   79 Br  0.58333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   80 Br  0.83333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   81 Br  0.08333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   82 Br  0.33333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   83 Br  0.58333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   84 Br  0.83333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   85 Br  0.08333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   86 Br  0.33333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   87 Br  0.58333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   88 Br  0.83333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   89 Br  0.08333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   90 Br  0.33333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   91 Br  0.58333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   92 Br  0.83333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   93 Br  0.08333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   94 Br  0.33333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   95 Br  0.58333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   96 Br  0.83333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   97 Br  0.16666666666667  0.08333333333333  0.24846253343926  79.904 > 97
   98 Br  0.41666666666667  0.08333333333333  0.24846253343926  79.904 > 97
   99 Br  0.66666666666667  0.08333333333333  0.24846253343926  79.904 > 97
  100 Br  0.91666666666667  0.08333333333333  0.24846253343926  79.904 > 97
  101 Br  0.16666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  102 Br  0.41666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  103 Br  0.66666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  104 Br  0.91666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  105 Br  0.16666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  106 Br  0.41666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  107 Br  0.66666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  108 Br  0.91666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  109 Br  0.16666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  110 Br  0.41666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  111 Br  0.66666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  112 Br  0.91666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  113 Br  0.16666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  114 Br  0.41666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  115 Br  0.66666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  116 Br  0.91666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  117 Br  0.16666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  118 Br  0.41666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  119 Br  0.66666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  120 Br  0.91666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  121 Br  0.16666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  122 Br  0.41666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  123 Br  0.66666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  124 Br  0.91666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  125 Br  0.16666666666667  0.83333333333333  0.74846253343926  79.904 > 97
  126 Br  0.41666666666667  0.83333333333333  0.74846253343926  79.904 > 97
  127 Br  0.66666666666667  0.83333333333333  0.74846253343926  79.904 > 97
  128 Br  0.91666666666667  0.83333333333333  0.74846253343926  79.904 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.8334873    0.0000000    0.0000000
            0.0000000    2.8334873    0.0000000
            0.0000000    0.0000000    2.8782661
-------------------------- Born effective charges --------------------------
    1 Li    1.0116389    0.0000000    0.0000000
            0.0000000    1.0116389    0.0000000
            0.0000000    0.0000000    1.0353697
    2 Li    1.0116389    0.0000000    0.0000000
            0.0000000    1.0116389    0.0000000
            0.0000000    0.0000000    1.0353697
    3 Br   -1.0116389   -0.0000000    0.0000000
            0.0000000   -1.0116389    0.0000000
            0.0000000    0.0000000   -1.0353697
    4 Br   -1.0116389   -0.0000000    0.0000000
            0.0000000   -1.0116389    0.0000000
            0.0000000    0.0000000   -1.0353697
----------------------------------------------------------------------------
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.00000011 (zzz) 0.00000011 (zzz) 0.00000011 (zzz)
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-09 04:38:13]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-09 04:38:13]-------------------------
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: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.134469640253094    0.000000000000000    0.000000000000000
  b   -2.067234820126547    3.580555739634689    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.609144710000000
Atomic positions (fractional):
    1 Li  0.33333333333333  0.66666666666667  0.37807293312149   6.941
    2 Li  0.66666666666667  0.33333333333333  0.87807293312149   6.941
    3 Br  0.33333333333333  0.66666666666667  0.99692506687852  79.904
    4 Br  0.66666666666667  0.33333333333333  0.49692506687852  79.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.537878561012377    0.000000000000000    0.000000000000000
  b   -8.268939280506189   14.322222958538756    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.218289420000000
Atomic positions (fractional):
    1 Li  0.08333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    3 Li  0.58333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    4 Li  0.83333333333333  0.16666666666667  0.18903646656074   6.941 > 1
    5 Li  0.08333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    6 Li  0.33333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    7 Li  0.58333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    8 Li  0.83333333333333  0.41666666666667  0.18903646656074   6.941 > 1
    9 Li  0.08333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   10 Li  0.33333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   11 Li  0.58333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   12 Li  0.83333333333333  0.66666666666667  0.18903646656074   6.941 > 1
   13 Li  0.08333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   14 Li  0.33333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   15 Li  0.58333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   16 Li  0.83333333333333  0.91666666666667  0.18903646656074   6.941 > 1
   17 Li  0.08333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   18 Li  0.33333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   19 Li  0.58333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   20 Li  0.83333333333333  0.16666666666667  0.68903646656074   6.941 > 1
   21 Li  0.08333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   22 Li  0.33333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   23 Li  0.58333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   24 Li  0.83333333333333  0.41666666666667  0.68903646656074   6.941 > 1
   25 Li  0.08333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   26 Li  0.33333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   27 Li  0.58333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   28 Li  0.83333333333333  0.66666666666667  0.68903646656074   6.941 > 1
   29 Li  0.08333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   30 Li  0.33333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   31 Li  0.58333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   32 Li  0.83333333333333  0.91666666666667  0.68903646656074   6.941 > 1
   33 Li  0.16666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   34 Li  0.41666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   35 Li  0.66666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   36 Li  0.91666666666667  0.08333333333333  0.43903646656074   6.941 > 33
   37 Li  0.16666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   38 Li  0.41666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   39 Li  0.66666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   40 Li  0.91666666666667  0.33333333333333  0.43903646656074   6.941 > 33
   41 Li  0.16666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   42 Li  0.41666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   43 Li  0.66666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   44 Li  0.91666666666667  0.58333333333333  0.43903646656074   6.941 > 33
   45 Li  0.16666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   46 Li  0.41666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   47 Li  0.66666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   48 Li  0.91666666666667  0.83333333333333  0.43903646656074   6.941 > 33
   49 Li  0.16666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   50 Li  0.41666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   51 Li  0.66666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   52 Li  0.91666666666667  0.08333333333333  0.93903646656074   6.941 > 33
   53 Li  0.16666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   54 Li  0.41666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   55 Li  0.66666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   56 Li  0.91666666666667  0.33333333333333  0.93903646656074   6.941 > 33
   57 Li  0.16666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   58 Li  0.41666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   59 Li  0.66666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   60 Li  0.91666666666667  0.58333333333333  0.93903646656074   6.941 > 33
   61 Li  0.16666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   62 Li  0.41666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   63 Li  0.66666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   64 Li  0.91666666666667  0.83333333333333  0.93903646656074   6.941 > 33
   65 Br  0.08333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   66 Br  0.33333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   67 Br  0.58333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   68 Br  0.83333333333333  0.16666666666667  0.49846253343926  79.904 > 65
   69 Br  0.08333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   70 Br  0.33333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   71 Br  0.58333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   72 Br  0.83333333333333  0.41666666666667  0.49846253343926  79.904 > 65
   73 Br  0.08333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   74 Br  0.33333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   75 Br  0.58333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   76 Br  0.83333333333333  0.66666666666667  0.49846253343926  79.904 > 65
   77 Br  0.08333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   78 Br  0.33333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   79 Br  0.58333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   80 Br  0.83333333333333  0.91666666666667  0.49846253343926  79.904 > 65
   81 Br  0.08333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   82 Br  0.33333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   83 Br  0.58333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   84 Br  0.83333333333333  0.16666666666667  0.99846253343926  79.904 > 65
   85 Br  0.08333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   86 Br  0.33333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   87 Br  0.58333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   88 Br  0.83333333333333  0.41666666666667  0.99846253343926  79.904 > 65
   89 Br  0.08333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   90 Br  0.33333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   91 Br  0.58333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   92 Br  0.83333333333333  0.66666666666667  0.99846253343926  79.904 > 65
   93 Br  0.08333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   94 Br  0.33333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   95 Br  0.58333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   96 Br  0.83333333333333  0.91666666666667  0.99846253343926  79.904 > 65
   97 Br  0.16666666666667  0.08333333333333  0.24846253343926  79.904 > 97
   98 Br  0.41666666666667  0.08333333333333  0.24846253343926  79.904 > 97
   99 Br  0.66666666666667  0.08333333333333  0.24846253343926  79.904 > 97
  100 Br  0.91666666666667  0.08333333333333  0.24846253343926  79.904 > 97
  101 Br  0.16666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  102 Br  0.41666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  103 Br  0.66666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  104 Br  0.91666666666667  0.33333333333333  0.24846253343926  79.904 > 97
  105 Br  0.16666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  106 Br  0.41666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  107 Br  0.66666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  108 Br  0.91666666666667  0.58333333333333  0.24846253343926  79.904 > 97
  109 Br  0.16666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  110 Br  0.41666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  111 Br  0.66666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  112 Br  0.91666666666667  0.83333333333333  0.24846253343926  79.904 > 97
  113 Br  0.16666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  114 Br  0.41666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  115 Br  0.66666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  116 Br  0.91666666666667  0.08333333333333  0.74846253343926  79.904 > 97
  117 Br  0.16666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  118 Br  0.41666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  119 Br  0.66666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  120 Br  0.91666666666667  0.33333333333333  0.74846253343926  79.904 > 97
  121 Br  0.16666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  122 Br  0.41666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  123 Br  0.66666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  124 Br  0.91666666666667  0.58333333333333  0.74846253343926  79.904 > 97
  125 Br  0.16666666666667  0.83333333333333  0.74846253343926  79.904 > 97
  126 Br  0.41666666666667  0.83333333333333  0.74846253343926  79.904 > 97
  127 Br  0.66666666666667  0.83333333333333  0.74846253343926  79.904 > 97
  128 Br  0.91666666666667  0.83333333333333  0.74846253343926  79.904 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.8334873    0.0000000    0.0000000
            0.0000000    2.8334873    0.0000000
            0.0000000    0.0000000    2.8782661
-------------------------- Born effective charges --------------------------
    1 Li    1.0116389    0.0000000    0.0000000
            0.0000000    1.0116389    0.0000000
            0.0000000    0.0000000    1.0353697
    2 Li    1.0116389    0.0000000    0.0000000
            0.0000000    1.0116389    0.0000000
            0.0000000    0.0000000    1.0353697
    3 Br   -1.0116389   -0.0000000    0.0000000
            0.0000000   -1.0116389    0.0000000
            0.0000000    0.0000000   -1.0353697
    4 Br   -1.0116389   -0.0000000    0.0000000
            0.0000000   -1.0116389    0.0000000
            0.0000000    0.0000000   -1.0353697
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000011 (zzz) 0.00000011 (zzz) 0.00000011 (zzz)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 14 14 8 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.90, Number of G-points: 311, Lambda: 0.16
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/120) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 120
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.279   (   0.000    0.000   -0.000)    0.000
   1.279   (   0.000    0.000    0.000)    0.000
   3.359   (   0.000    0.000    0.000)    0.000
   8.458   (   0.000   -0.000    0.000)    0.000
   8.834   (  -0.000    0.000   -0.000)    0.000
   8.834   (   0.000   -0.000    0.000)    0.000
   8.938   (   0.000   -0.000    0.000)    0.000
   8.938   (   0.000   -0.000    0.000)    0.000
  10.117   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/120) =======================
q-point: ( 0.07  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.264   (  11.317    6.534    0.000)   13.068
   0.288   (  12.638    7.297    0.000)   14.593
   0.630   (  26.830   15.490    0.000)   30.981
   1.309   (   2.599    1.500    0.000)    3.001
   1.443   (  12.934    7.467    0.000)   14.935
   3.340   (  -1.648   -0.951    0.000)    1.903
   8.454   (  -0.273   -0.157    0.000)    0.315
   8.790   (  -2.365   -1.365    0.000)    2.731
   8.866   (   2.388    1.379    0.000)    2.758
   8.969   (   1.278    0.738    0.000)    1.476
  10.085   (  -2.475   -1.429    0.000)    2.858
  11.416   (  -2.513   -1.451    0.000)    2.902
======================= Grid point 2 (3/120) =======================
q-point: ( 0.14  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.517   (  10.488    6.055    0.000)   12.110
   0.579   (  12.680    7.321    0.000)   14.642
   1.228   (  24.922   14.389    0.000)   28.778
   1.399   (   5.109    2.950    0.000)    5.899
   1.803   (  16.931    9.775    0.000)   19.550
   3.286   (  -2.970   -1.715    0.000)    3.430
   8.447   (  -0.333   -0.192    0.000)    0.384
   8.757   (  -0.419   -0.242    0.000)    0.484
   8.913   (   0.904    0.522    0.000)    1.044
   8.967   (  -1.496   -0.864    0.000)    1.727
  10.028   (  -1.789   -1.033    0.000)    2.066
  11.338   (  -4.051   -2.339    0.000)    4.678
======================= Grid point 3 (4/120) =======================
q-point: ( 0.21  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.741   (   8.839    5.103    0.000)   10.206
   0.872   (  12.691    7.327    0.000)   14.655
   1.534   (   6.260    3.614    0.000)    7.228
   1.770   (  21.923   12.657    0.000)   25.315
   2.174   (  14.588    8.422    0.000)   16.845
   3.210   (  -3.332   -1.923    0.000)    3.847
   8.442   (  -0.056   -0.032    0.000)    0.065
   8.763   (   0.515    0.297    0.000)    0.594
   8.885   (  -3.282   -1.895    0.000)    3.790
   8.912   (  -2.795   -1.614    0.000)    3.227
  10.025   (   1.568    0.905    0.000)    1.810
  11.236   (  -4.731   -2.731    0.000)    5.462
======================= Grid point 4 (5/120) =======================
q-point: ( 0.29  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.919   (   6.581    3.799    0.000)    7.599
   1.161   (  12.360    7.136    0.000)   14.272
   1.671   (   5.327    3.075    0.000)    6.151
   2.229   (  17.730   10.236    0.000)   20.473
   2.446   (   8.448    4.877    0.000)    9.755
   3.151   (  -1.233   -0.712    0.000)    1.423
   8.447   (   0.559    0.323    0.000)    0.646
   8.756   (  -1.273   -0.735    0.000)    1.470
   8.785   (  -4.784   -2.762    0.000)    5.524
   8.868   (  -0.733   -0.423    0.000)    0.847
  10.089   (   3.781    2.183    0.000)    4.365
  11.118   (  -5.719   -3.302    0.000)    6.604
======================= Grid point 5 (6/120) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.043   (   4.234    2.445    0.000)    4.889
   1.435   (  11.147    6.435    0.000)   12.871
   1.769   (   3.087    1.782    0.000)    3.565
   2.501   (   1.762    1.017    0.000)    2.035
   2.625   (  11.383    6.572    0.000)   13.144
   3.170   (   2.607    1.505    0.000)    3.011
   8.469   (   1.317    0.760    0.000)    1.521
   8.684   (  -3.758   -2.170    0.000)    4.339
   8.725   (  -0.987   -0.570    0.000)    1.140
   8.861   (  -0.261   -0.151    0.000)    0.302
  10.195   (   5.463    3.154    0.000)    6.308
  10.965   (  -7.612   -4.395    0.000)    8.790
======================= Grid point 6 (7/120) =======================
q-point: ( 0.43  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.116   (   2.096    1.210    0.000)    2.421
   1.659   (   7.770    4.486    0.000)    8.972
   1.815   (   1.051    0.607    0.000)    1.213
   2.457   (  -3.359   -1.939    0.000)    3.878
   2.858   (   7.464    4.310    0.000)    8.619
   3.237   (   2.450    1.415    0.000)    2.829
   8.505   (   1.802    1.041    0.000)    2.081
   8.612   (  -2.576   -1.487    0.000)    2.974
   8.712   (  -0.276   -0.159    0.000)    0.318
   8.854   (  -0.225   -0.130    0.000)    0.260
  10.344   (   7.461    4.307    0.000)    8.615
  10.770   (  -9.093   -5.250    0.000)   10.500
======================= Grid point 7 (8/120) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 160
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.140   (   0.000    0.000    0.000)    0.000
   1.757   (   0.000    0.000    0.000)    0.000
   1.825   (   0.000    0.000    0.000)    0.000
   2.407   (   0.000    0.000    0.000)    0.000
   2.947   (   0.000    0.000    0.000)    0.000
   3.266   (   0.000    0.000    0.000)    0.000
   8.540   (   0.000    0.000    0.000)    0.000
   8.570   (  -0.000   -0.000    0.000)    0.000
   8.709   (   0.000    0.000    0.000)    0.000
   8.852   (   0.000    0.000    0.000)    0.000
  10.500   (  -0.000   -0.000    0.000)    0.000
  10.594   (   0.000    0.000    0.000)    0.000
======================= Grid point 16 (9/120) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.480   (   7.398   12.814    0.000)   14.796
   0.500   (   7.260   12.575    0.000)   14.520
   1.056   (  13.939   24.143    0.000)   27.878
   1.424   (   4.184    7.247    0.000)    8.368
   1.654   (   8.599   14.893    0.000)   17.197
   3.303   (  -1.546   -2.679    0.000)    3.093
   8.449   (  -0.199   -0.345    0.000)    0.398
   8.776   (  -1.330   -2.304    0.000)    2.660
   8.819   (  -0.282   -0.489    0.000)    0.565
   9.056   (   2.279    3.947    0.000)    4.558
  10.041   (  -1.416   -2.453    0.000)    2.832
  11.352   (  -2.700   -4.676    0.000)    5.400
======================= Grid point 17 (10/120) =======================
q-point: ( 0.14  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.730   (   5.706   13.780    0.000)   14.915
   0.765   (   9.641   10.824    0.000)   14.495
   1.492   (  10.302   11.108    0.000)   15.150
   1.620   (   9.650   17.540    0.000)   20.020
   1.995   (  14.706    9.931    0.000)   17.745
   3.237   (  -2.574   -3.019    0.000)    3.967
   8.443   (  -0.153   -0.171    0.000)    0.230
   8.737   (   0.525   -2.982    0.000)    3.028
   8.818   (   0.264   -1.867    0.000)    1.885
   9.085   (  -4.336    5.325    0.000)    6.867
  10.010   (   0.123   -0.529    0.000)    0.543
  11.243   (  -2.947   -6.972    0.000)    7.569
======================= Grid point 18 (11/120) =======================
q-point: ( 0.21  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.943   (   2.713   13.914    0.000)   14.176
   1.043   (  10.607    9.629    0.000)   14.326
   1.677   (   2.744   10.040    0.000)   10.409
   2.043   (  16.449   14.093    0.000)   21.660
   2.314   (  12.242    5.582    0.000)   13.454
   3.165   (  -2.103   -2.404    0.000)    3.195
   8.443   (   0.224    0.172    0.000)    0.282
   8.719   (   1.503   -5.213    0.000)    5.425
   8.792   (  -1.982    0.569    0.000)    2.062
   9.027   (  -7.872    5.545    0.000)    9.629
  10.039   (   3.000    0.393    0.000)    3.026
  11.118   (  -2.351   -8.841    0.000)    9.148
======================= Grid point 19 (12/120) =======================
q-point: ( 0.29  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.099   (  -0.274   13.233    0.000)   13.236
   1.318   (  10.586    8.491    0.000)   13.571
   1.804   (   0.332    9.788    0.000)    9.794
   2.414   (  11.662    6.917    0.000)   13.559
   2.514   (   7.732    3.227    0.000)    8.379
   3.139   (   1.261   -0.481    0.000)    1.350
   8.456   (   0.895    0.566    0.000)    1.060
   8.679   (  -0.343   -4.319    0.000)    4.332
   8.766   (  -2.480    1.664    0.000)    2.986
   8.952   (  -6.389    4.920    0.000)    8.064
  10.114   (   5.309    0.179    0.000)    5.312
  10.980   (  -2.943  -10.146    0.000)   10.564
======================= Grid point 20 (13/120) =======================
q-point: ( 0.36  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.199   (  -2.820   12.209    0.000)   12.531
   1.568   (   8.989    6.933    0.000)   11.352
   1.880   (  -2.348    8.890    0.000)    9.195
   2.478   (  -1.396   -2.630    0.000)    2.977
   2.751   (  10.180    5.568    0.000)   11.603
   3.189   (   3.678    0.232    0.000)    3.685
   8.485   (   1.575    0.791    0.000)    1.762
   8.629   (  -1.505   -2.368    0.000)    2.806
   8.743   (  -2.129    1.991    0.000)    2.915
   8.906   (  -3.823    3.522    0.000)    5.198
  10.229   (   7.919    0.202    0.000)    7.921
  10.809   (  -4.719  -10.755    0.000)   11.745
======================= Grid point 21 (14/120) =======================
q-point: ( 0.43  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.247   (  -4.992   11.065    0.000)   12.139
   1.743   (   3.840    3.945    0.000)    5.505
   1.905   (  -4.136    7.893    0.000)    8.911
   2.413   (  -1.749   -2.226    0.000)    2.831
   2.916   (   4.252    1.471    0.000)    4.499
   3.240   (   2.243   -1.095    0.000)    2.496
   8.522   (   1.695    0.454    0.000)    1.755
   8.584   (  -1.565   -1.008    0.000)    1.862
   8.732   (  -1.467    1.945    0.000)    2.436
   8.887   (  -2.201    2.974    0.000)    3.699
  10.389   (   9.870   -0.043    0.000)    9.870
  10.612   (  -5.321   -9.900    0.000)   11.240
======================= Grid point 31 (15/120) =======================
q-point: ( 0.14  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.998   (   7.089   12.278    0.000)   14.178
   1.010   (   7.684   13.309    0.000)   15.368
   1.710   (   6.083   10.536    0.000)   12.166
   1.983   (  10.383   17.983    0.000)   20.765
   2.205   (   6.178   10.701    0.000)   12.356
   3.172   (  -1.882   -3.260    0.000)    3.765
   8.442   (   0.071    0.123    0.000)    0.142
   8.676   (  -1.299   -2.250    0.000)    2.599
   8.817   (   0.283    0.490    0.000)    0.565
   9.140   (   0.087    0.151    0.000)    0.174
  10.012   (   0.420    0.727    0.000)    0.839
  11.081   (  -5.041   -8.732    0.000)   10.082
======================= Grid point 32 (16/120) =======================
q-point: ( 0.21  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.239   (   3.581   14.065    0.000)   14.514
   1.251   (   7.453   11.625    0.000)   13.809
   1.901   (   2.839   11.725    0.000)   12.064
   2.316   (   9.104   12.658    0.000)   15.592
   2.400   (   8.426    3.398    0.000)    9.086
   3.122   (  -0.542   -1.786    0.000)    1.866
   8.449   (   0.478    0.501    0.000)    0.693
   8.645   (  -0.028   -2.212    0.000)    2.213
   8.821   (  -1.172    1.713    0.000)    2.075
   9.104   (  -4.879    2.010    0.000)    5.277
  10.039   (   2.516   -0.231    0.000)    2.526
  10.908   (  -3.136  -11.442    0.000)   11.864
======================= Grid point 33 (17/120) =======================
q-point: ( 0.29  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 500
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.396   (  -2.001   15.366    0.000)   15.495
   1.497   (   8.373    9.337    0.000)   12.541
   2.034   (  -1.636   12.481    0.000)   12.588
   2.446   (   2.044   -2.135    0.000)    2.956
   2.636   (  10.100    7.495    0.000)   12.577
   3.128   (   2.589   -0.783    0.000)    2.705
   8.469   (   1.072    0.788    0.000)    1.331
   8.624   (  -0.411   -1.541    0.000)    1.595
   8.810   (  -2.393    2.473    0.000)    3.442
   9.032   (  -5.739    2.831    0.000)    6.399
  10.099   (   5.553   -1.418    0.000)    5.732
  10.741   (  -2.352  -12.965    0.000)   13.176
======================= Grid point 34 (18/120) =======================
q-point: ( 0.36  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.482   (  -5.475   15.020    0.000)   15.986
   1.711   (   6.087    7.291    0.000)    9.498
   2.091   (  -5.217   11.119    0.000)   12.282
   2.415   (  -1.254   -3.105    0.000)    3.348
   2.840   (   6.729    3.356    0.000)    7.519
   3.173   (   3.445   -1.985    0.000)    3.976
   8.500   (   1.437    0.816    0.000)    1.652
   8.598   (  -1.142   -0.898    0.000)    1.453
   8.793   (  -2.510    2.754    0.000)    3.727
   8.975   (  -3.995    2.915    0.000)    4.945
  10.207   (   8.912   -2.171    0.000)    9.173
  10.566   (  -2.927  -12.800    0.000)   13.130
======================= Grid point 35 (19/120) =======================
q-point: ( 0.43  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.509   (  -8.072   13.981    0.000)   16.144
   1.812   (  -1.676    2.903    0.000)    3.352
   2.093   (  -5.547    9.608    0.000)   11.094
   2.382   (   0.163   -0.283    0.000)    0.327
   2.918   (   0.609   -1.056    0.000)    1.219
   3.194   (   2.017   -3.494    0.000)    4.034
   8.523   (   0.110   -0.190    0.000)    0.220
   8.577   (  -0.140    0.242    0.000)    0.279
   8.784   (  -1.795    3.110    0.000)    3.591
   8.955   (  -1.978    3.426    0.000)    3.956
  10.325   (   4.157   -7.200    0.000)    8.314
  10.424   (   4.032   -6.984    0.000)    8.064
======================= Grid point 47 (20/120) =======================
q-point: ( 0.21  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.468   (   6.413   11.108    0.000)   12.826
   1.506   (   6.453   11.177    0.000)   12.907
   2.124   (   5.871   10.168    0.000)   11.741
   2.406   (  -0.432   -0.748    0.000)    0.864
   2.578   (   6.562   11.365    0.000)   13.124
   3.100   (  -0.224   -0.388    0.000)    0.448
   8.464   (   0.568    0.985    0.000)    1.137
   8.619   (  -0.431   -0.747    0.000)    0.863
   8.848   (   0.491    0.850    0.000)    0.981
   9.107   (  -0.892   -1.545    0.000)    1.785
  10.039   (   0.288    0.499    0.000)    0.576
  10.686   (  -6.087  -10.543    0.000)   12.174
======================= Grid point 48 (21/120) =======================
q-point: ( 0.29  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.673   (   2.054   11.426    0.000)   11.609
   1.691   (   4.460   10.330    0.000)   11.252
   2.270   (  -2.829    9.416    0.000)    9.832
   2.380   (   2.100   -2.818    0.000)    3.515
   2.775   (   5.723    6.261    0.000)    8.482
   3.102   (   1.626   -1.923    0.000)    2.518
   8.489   (   1.031    1.246    0.000)    1.617
   8.605   (  -0.570   -0.642    0.000)    0.858
   8.853   (  -1.266    1.669    0.000)    2.094
   9.058   (  -3.149   -0.173    0.000)    3.153
  10.072   (   3.996   -0.897    0.000)    4.095
  10.492   (  -3.701  -11.488    0.000)   12.070
======================= Grid point 49 (22/120) =======================
q-point: ( 0.36  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.759   (  -4.620   12.014    0.000)   12.871
   1.856   (   2.448    7.171    0.000)    7.578
   2.222   (  -3.510   -0.815    0.000)    3.603
   2.430   (  -1.577    6.812    0.000)    6.992
   2.889   (   2.510    1.727    0.000)    3.047
   3.110   (   2.690   -4.221    0.000)    5.005
   8.520   (   0.771    1.164    0.000)    1.396
   8.585   (  -0.780   -0.520    0.000)    0.937
   8.845   (  -1.737    2.313    0.000)    2.892
   9.014   (  -2.049    0.986    0.000)    2.274
  10.157   (   7.097   -2.479    0.000)    7.518
  10.330   (  -1.969  -10.108    0.000)   10.298
======================= Grid point 62 (23/120) =======================
q-point: ( 0.29  0.29  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.867   (   5.007    8.673    0.000)   10.014
   1.871   (   3.761    6.514    0.000)    7.522
   2.240   (  -3.577   -6.196    0.000)    7.155
   2.492   (   4.331    7.502    0.000)    8.662
   2.879   (   2.304    3.990    0.000)    4.607
   3.053   (  -1.595   -2.762    0.000)    3.189
   8.519   (   0.965    1.672    0.000)    1.930
   8.590   (  -0.547   -0.948    0.000)    1.095
   8.875   (   0.286    0.496    0.000)    0.573
   9.038   (  -0.793   -1.373    0.000)    1.586
  10.077   (   0.962    1.666    0.000)    1.923
  10.294   (  -4.739   -8.208    0.000)    9.478
======================= Grid point 63 (24/120) =======================
q-point: ( 0.36  0.29  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.935   (  -2.836    4.912    0.000)    5.672
   1.988   (  -3.246    5.622    0.000)    6.492
   2.144   (   2.819   -4.883    0.000)    5.638
   2.574   (  -2.818    4.881    0.000)    5.636
   2.924   (  -1.263    2.187    0.000)    2.526
   3.016   (   2.633   -4.560    0.000)    5.266
   8.544   (  -0.646    1.119    0.000)    1.293
   8.573   (   0.411   -0.713    0.000)    0.823
   8.878   (  -0.510    0.884    0.000)    1.021
   9.021   (  -0.015    0.026    0.000)    0.030
  10.122   (   0.425   -0.736    0.000)    0.850
  10.184   (   2.377   -4.117    0.000)    4.754
======================= Grid point 211 (25/120) =======================
q-point: ( 0.00  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 96
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.269   (   0.000    0.000   13.965)   13.965
   0.269   (   0.000    0.000   13.965)   13.965
   0.615   (   0.000   -0.000   32.198)   32.198
   1.259   (  -0.000   -0.000   -2.178)    2.178
   1.259   (   0.000   -0.000   -2.178)    2.178
   3.284   (  -0.000    0.000   -7.902)    7.902
   8.838   (  -0.000    0.000    0.462)    0.462
   8.838   (  -0.000    0.000    0.462)    0.462
   8.935   (   0.000   -0.000   -0.365)    0.365
   8.935   (   0.000    0.000   -0.365)    0.365
  10.163   (   0.000    0.000    4.776)    4.776
  11.113   (  -0.000    0.000   -3.546)    3.546
======================= Grid point 212 (26/120) =======================
q-point: ( 0.07  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.378   (   8.017    4.629   10.063)   13.674
   0.381   (   0.831    0.480   20.142)   20.165
   0.829   (  22.108   12.764   16.259)   30.266
   1.290   (   2.685    1.550   -2.052)    3.718
   1.413   (  12.573    7.259   -2.630)   14.754
   3.265   (  -1.599   -0.923   -7.845)    8.059
   8.657   ( -10.124   -5.845   11.196)   16.187
   8.798   (  -2.233   -1.289    0.749)    2.685
   8.877   (   1.606    0.927    1.099)    2.156
   8.963   (   1.135    0.655   -0.665)    1.470
  10.140   (  -1.696   -0.979    5.670)    5.998
  11.260   (   4.274    2.467  -10.807)   11.881
======================= Grid point 213 (27/120) =======================
q-point: ( 0.14  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.540   (  11.349    6.552    3.977)   13.694
   0.586   (   9.423    5.440    6.799)   12.830
   1.364   (  23.034   13.299   11.371)   28.926
   1.381   (   5.142    2.969   -1.888)    6.231
   1.768   (  16.732    9.660   -3.388)   19.615
   3.214   (  -2.699   -1.558   -7.463)    8.087
   8.522   (  -3.060   -1.767    6.548)    7.440
   8.765   (  -0.469   -0.271    0.858)    1.015
   8.905   (   0.037    0.021   -0.742)    0.743
   8.960   (  -1.464   -0.845   -0.809)    1.874
  10.107   (  -0.638   -0.368    7.935)    7.969
  11.275   (  -1.787   -1.032   -6.134)    6.472
======================= Grid point 214 (28/120) =======================
q-point: ( 0.21  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.794   (   8.389    4.843    5.363)   11.073
   0.825   (  12.613    7.282   -1.644)   14.657
   1.516   (   6.234    3.599   -1.886)    7.441
   1.870   (  20.697   11.950    9.027)   25.547
   2.132   (  14.049    8.111   -4.286)   16.779
   3.151   (  -2.449   -1.414   -6.091)    6.716
   8.480   (  -0.969   -0.560    3.578)    3.749
   8.769   (   0.383    0.221    0.594)    0.741
   8.861   (  -3.672   -2.120   -2.070)    4.718
   8.906   (  -2.677   -1.545   -0.583)    3.145
  10.126   (   2.221    1.282    9.973)   10.298
  11.210   (  -3.555   -2.052   -2.776)    4.955
======================= Grid point 215 (29/120) =======================
q-point: ( 0.29  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.965   (   6.361    3.672    4.672)    8.705
   1.108   (  11.892    6.866   -3.854)   14.262
   1.652   (   5.313    3.067   -1.941)    6.434
   2.295   (  14.318    8.267    3.513)   16.903
   2.393   (   9.647    5.570   -3.019)   11.541
   3.123   (   0.619    0.357   -2.772)    2.862
   8.470   (   0.096    0.056    2.152)    2.154
   8.758   (  -4.725   -2.728   -2.292)    5.917
   8.760   (  -1.259   -0.727    0.436)    1.517
   8.864   (  -0.760   -0.439   -0.452)    0.987
  10.198   (   3.795    2.191   10.647)   11.513
  11.114   (  -4.840   -2.794   -0.463)    5.608
======================= Grid point 216 (30/120) =======================
q-point: ( 0.36  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.086   (   4.131    2.385    4.330)    6.442
   1.368   (  10.483    6.053   -5.855)   13.447
   1.751   (   3.116    1.799   -1.924)    4.080
   2.406   (  -1.137   -0.656  -10.133)   10.218
   2.687   (  12.776    7.376    6.476)   16.111
   3.183   (   3.948    2.279    0.973)    4.661
   8.482   (   0.923    0.533    1.175)    1.587
   8.661   (  -3.536   -2.041   -1.945)    4.522
   8.730   (  -0.970   -0.560    0.508)    1.230
   8.855   (  -0.299   -0.172   -0.561)    0.659
  10.298   (   5.005    2.890   10.133)   11.665
  10.981   (  -6.724   -3.882    1.582)    7.924
======================= Grid point 217 (31/120) =======================
q-point: ( 0.43  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.157   (   2.051    1.184    4.168)    4.794
   1.578   (   7.335    4.235   -7.479)   11.299
   1.798   (   1.097    0.633   -1.823)    2.220
   2.336   (  -3.793   -2.190  -11.842)   12.626
   2.922   (   7.309    4.220    6.176)   10.459
   3.269   (   2.923    1.687    2.733)    4.342
   8.511   (   1.475    0.852    0.281)    1.726
   8.595   (  -2.335   -1.348   -1.385)    3.031
   8.717   (  -0.284   -0.164    0.508)    0.605
   8.848   (  -0.237   -0.137   -0.612)    0.671
  10.432   (   6.699    3.868    8.703)   11.644
  10.807   (  -8.163   -4.713    3.739)   10.140
======================= Grid point 218 (32/120) =======================
q-point: (-0.50  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 228
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.181   (   0.000    0.000    4.114)    4.114
   1.672   (   0.000    0.000   -8.079)    8.079
   1.809   (   0.000    0.000   -1.758)    1.758
   2.282   (  -0.000   -0.000  -12.224)   12.224
   3.009   (   0.000    0.000    5.927)    5.927
   3.303   (   0.000    0.000    3.228)    3.228
   8.547   (   0.000    0.000    0.637)    0.637
   8.550   (  -0.000   -0.000   -1.861)    1.861
   8.714   (   0.000    0.000    0.499)    0.499
   8.846   (   0.000    0.000   -0.628)    0.628
  10.578   (  -0.000   -0.000    7.636)    7.636
  10.643   (   0.000    0.000    5.016)    5.016
======================= Grid point 227 (33/120) =======================
q-point: ( 0.07  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.469   (   5.495    9.518    7.194)   13.136
   0.554   (   6.540   11.327    7.151)   14.907
   1.204   (  12.354   21.397   11.725)   27.348
   1.397   (   4.616    7.996   -2.158)    9.482
   1.630   (   8.478   14.684   -2.452)   17.132
   3.230   (  -1.444   -2.501   -7.630)    8.159
   8.547   (  -2.341   -4.054    8.151)    9.399
   8.774   (  -1.552   -2.689   -0.150)    3.108
   8.830   (  -0.434   -0.752    1.033)    1.350
   9.040   (   2.002    3.468   -1.585)    4.307
  10.113   (  -0.754   -1.306    7.261)    7.416
  11.274   (  -0.968   -1.676   -7.203)    7.458
======================= Grid point 228 (34/120) =======================
q-point: ( 0.14  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.711   (   9.372   10.506   -0.832)   14.103
   0.787   (   5.828   12.979    5.413)   15.222
   1.506   (   4.529    9.280   -0.863)   10.363
   1.698   (  14.641   18.186    8.654)   24.899
   1.962   (  14.014    9.841   -3.459)   17.470
   3.171   (  -2.124   -2.543   -6.794)    7.559
   8.492   (  -1.322   -1.171    4.688)    5.009
   8.729   (   0.692   -4.135   -0.729)    4.255
   8.821   (  -0.305   -1.275    0.378)    1.365
   9.064   (  -4.104    4.913   -2.013)    6.711
  10.105   (   0.796    0.285    9.407)    9.445
  11.208   (  -1.962   -5.542   -3.552)    6.869
======================= Grid point 229 (35/120) =======================
q-point: ( 0.21  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.965   (   6.407   10.739   -0.356)   12.510
   1.013   (   6.803   11.802    1.569)   13.713
   1.661   (   2.432   10.230   -1.817)   10.671
   2.126   (  15.539   12.900    7.133)   21.419
   2.265   (  11.579    5.358   -4.546)   13.544
   3.121   (  -0.836   -1.381   -4.482)    4.764
   8.472   (  -0.355   -0.244    2.736)    2.770
   8.700   (   1.264   -5.835   -1.702)    6.208
   8.796   (  -1.998    0.873    0.450)    2.226
   9.010   (  -7.363    5.447   -1.536)    9.287
  10.146   (   3.177    0.598   10.547)   11.031
  11.109   (  -1.709   -7.874   -1.055)    8.126
======================= Grid point 230 (36/120) =======================
q-point: ( 0.29  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.133   (   0.468   12.193    3.048)   12.577
   1.261   (   9.494    8.434   -4.190)   13.373
   1.786   (   0.244    9.846   -1.919)   10.035
   2.380   (   4.010    0.377   -7.598)    8.599
   2.536   (  13.431    8.563    5.765)   16.940
   3.135   (   2.979    0.706   -0.564)    3.113
   8.473   (   0.525    0.254    1.595)    1.698
   8.655   (  -0.324   -4.427   -2.112)    4.916
   8.771   (  -2.496    1.722    0.530)    3.079
   8.943   (  -6.037    4.852   -0.882)    7.795
  10.220   (   5.067   -0.019   10.464)   11.627
  10.989   (  -2.330   -9.363    0.898)    9.691
======================= Grid point 231 (37/120) =======================
q-point: ( 0.36  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.234   (  -2.446   11.581    3.446)   12.328
   1.494   (   8.339    6.567   -6.637)   12.519
   1.861   (  -2.333    8.760   -2.012)    9.286
   2.366   (  -2.340   -2.975  -11.105)   11.732
   2.818   (  10.152    5.568    6.389)   13.224
   3.215   (   4.436    0.799    2.153)    4.996
   8.493   (   1.251    0.520    0.704)    1.527
   8.607   (  -1.242   -2.431   -1.909)    3.331
   8.748   (  -2.137    2.043    0.566)    3.010
   8.900   (  -3.740    3.521   -0.616)    5.173
  10.325   (   7.319   -0.235    9.432)   11.941
  10.838   (  -3.962  -10.065    2.839)   11.183
======================= Grid point 232 (38/120) =======================
q-point: ( 0.43  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.282   (  -4.702   10.557    3.510)   12.078
   1.658   (   3.647    3.719   -7.934)    9.491
   1.885   (  -3.956    7.490   -2.169)    8.744
   2.291   (  -2.093   -1.925  -11.745)   12.085
   2.980   (   4.124    1.463    6.023)    7.445
   3.277   (   2.417   -0.882    3.180)    4.091
   8.524   (   1.677    0.410    0.049)    1.727
   8.565   (  -1.445   -1.270   -1.543)    2.466
   8.737   (  -1.499    2.018    0.583)    2.581
   8.882   (  -2.218    3.019   -0.561)    3.788
  10.468   (   9.270   -0.424    7.809)   12.128
  10.660   (  -4.686   -9.446    4.822)   11.595
======================= Grid point 242 (39/120) =======================
q-point: ( 0.14  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.939   (   6.850   11.865   -3.629)   14.173
   1.054   (   7.492   12.977    4.507)   15.647
   1.706   (   5.841   10.117   -1.122)   11.735
   2.063   (  10.297   17.836    7.837)   22.036
   2.159   (   5.557    9.625   -4.862)   12.131
   3.122   (  -1.215   -2.105   -5.095)    5.645
   8.475   (  -0.347   -0.601    3.185)    3.260
   8.655   (  -1.526   -2.643   -1.910)    3.600
   8.822   (   0.278    0.481    0.598)    0.817
   9.117   (   0.198    0.343   -2.166)    2.202
  10.119   (   0.611    1.058   10.489)   10.560
  11.071   (  -4.410   -7.638   -1.052)    8.882
======================= Grid point 243 (40/120) =======================
q-point: ( 0.21  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.174   (   7.352   10.313   -5.634)   13.861
   1.286   (   3.393   14.204    4.019)   15.146
   1.891   (   2.310   12.068   -1.368)   12.363
   2.306   (   4.657    2.621   -6.614)    8.503
   2.412   (  11.696   11.593    6.113)   17.566
   3.103   (   0.824   -0.260   -1.884)    2.072
   8.470   (   0.117    0.054    1.938)    1.942
   8.619   (  -0.053   -2.402   -2.343)    3.356
   8.827   (  -1.182    1.749    0.621)    2.201
   9.088   (  -4.500    2.092   -1.585)    5.210
  10.148   (   2.380   -0.356   10.590)   10.860
  10.918   (  -2.695  -10.534    0.887)   10.909
======================= Grid point 244 (41/120) =======================
q-point: ( 0.29  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.383   (   3.660   10.863   -5.071)   12.535
   1.463   (   2.495   12.551    1.279)   12.861
   2.016   (  -2.122   12.346   -2.037)   12.692
   2.348   (   0.995   -3.065   -9.926)   10.436
   2.701   (  10.240    7.765    6.473)   14.389
   3.146   (   3.666    0.212    1.416)    3.936
   8.480   (   0.782    0.405    0.982)    1.319
   8.598   (  -0.291   -1.580   -2.344)    2.841
   8.816   (  -2.395    2.504    0.632)    3.522
   9.022   (  -5.426    2.821   -0.982)    6.194
  10.200   (   5.122   -1.806    9.864)   11.260
  10.766   (  -1.881  -12.179    2.515)   12.578
======================= Grid point 245 (42/120) =======================
q-point: ( 0.36  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.497   (  -3.912   13.376    0.125)   13.937
   1.634   (   4.598    7.834   -5.546)   10.643
   2.059   (  -5.238    9.092   -3.912)   11.199
   2.310   (  -1.750   -1.329   -9.684)    9.931
   2.906   (   6.636    3.344    6.206)    9.681
   3.209   (   3.773   -1.516    3.080)    5.101
   8.503   (   1.247    0.470    0.198)    1.348
   8.574   (  -0.955   -0.909   -2.096)    2.476
   8.799   (  -2.493    2.793    0.666)    3.803
   8.969   (  -3.860    2.920   -0.628)    4.881
  10.293   (   8.267   -2.653    8.494)   12.146
  10.608   (  -2.319  -12.166    4.193)   13.075
======================= Grid point 246 (43/120) =======================
q-point: ( 0.43  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.529   (  -7.355   12.740    1.185)   14.759
   1.727   (  -1.931    3.345   -7.165)    8.140
   2.046   (  -3.820    6.616   -5.715)    9.540
   2.286   (  -1.524    2.640   -8.330)    8.870
   2.983   (   0.553   -0.958    6.036)    6.136
   3.235   (   1.906   -3.302    3.547)    5.207
   8.522   (   0.270   -0.468   -0.070)    0.545
   8.555   (  -0.131    0.228   -2.018)    2.035
   8.791   (  -1.811    3.137    0.701)    3.690
   8.950   (  -1.978    3.425   -0.521)    3.989
  10.401   (   4.201   -7.276    7.474)   11.245
  10.477   (   3.949   -6.840    5.270)    9.495
======================= Grid point 258 (44/120) =======================
q-point: ( 0.21  0.21  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.379   (   5.863   10.156   -7.931)   14.157
   1.542   (   6.372   11.037    3.673)   13.263
   2.119   (   5.878   10.181   -0.946)   11.794
   2.304   (  -1.144   -1.981  -10.246)   10.499
   2.646   (   6.436   11.147    6.556)   14.445
   3.112   (   0.611    1.058    0.961)    1.555
   8.475   (   0.308    0.533    1.074)    1.238
   8.592   (  -0.413   -0.715   -2.483)    2.617
   8.854   (   0.490    0.849    0.653)    1.178
   9.093   (  -0.778   -1.347   -1.332)    2.048
  10.141   (   0.041    0.071   10.034)   10.034
  10.712   (  -5.660   -9.803    2.548)   11.602
======================= Grid point 259 (45/120) =======================
q-point: ( 0.29  0.21  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.573   (   5.246    8.943   -9.766)   14.244
   1.719   (   1.059   11.675    3.384)   12.201
   2.205   (  -4.913    0.713   -9.144)   10.405
   2.317   (   3.448    5.382   -3.771)    7.421
   2.841   (   5.742    6.152    6.185)   10.444
   3.138   (   2.121   -1.116    3.007)    3.845
   8.491   (   0.788    0.809    0.223)    1.151
   8.579   (  -0.461   -0.547   -2.325)    2.433
   8.859   (  -1.244    1.659    0.644)    2.171
   9.049   (  -2.941   -0.063   -0.900)    3.077
  10.163   (   3.516   -1.428    8.953)    9.724
  10.532   (  -3.288  -10.832    3.954)   11.991
======================= Grid point 260 (46/120) =======================
q-point: ( 0.36  0.21  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.714   (  -0.020    8.238   -9.135)   12.301
   1.821   (  -1.803   10.064    1.664)   10.358
   2.107   (  -2.786   -3.457  -11.545)   12.369
   2.403   (  -2.652    9.147   -2.437)    9.831
   2.955   (   2.547    1.771    6.089)    6.834
   3.155   (   2.660   -3.732    3.950)    6.051
   8.515   (   0.716    0.772   -0.476)    1.156
   8.563   (  -0.698   -0.412   -2.053)    2.207
   8.852   (  -1.698    2.293    0.681)    2.933
   9.008   (  -1.966    1.044   -0.563)    2.296
  10.234   (   6.656   -2.885    7.603)   10.509
  10.383   (  -1.609   -9.644    5.275)   11.110
======================= Grid point 273 (47/120) =======================
q-point: ( 0.29  0.29  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.745   (   4.667    8.083  -11.515)   14.822
   1.902   (   3.700    6.409    3.035)    7.998
   2.112   (  -3.688   -6.389  -12.161)   14.223
   2.481   (   4.151    7.190   -1.361)    8.413
   2.943   (   2.318    4.015    5.913)    7.514
   3.103   (  -1.244   -2.155    4.319)    4.985
   8.513   (   0.720    1.247   -0.603)    1.561
   8.567   (  -0.438   -0.758   -2.087)    2.264
   8.880   (   0.289    0.501    0.625)    0.852
   9.032   (  -0.703   -1.218   -0.609)    1.533
  10.157   (   0.606    1.049    7.825)    7.918
  10.346   (  -4.405   -7.630    5.155)   10.208
======================= Grid point 274 (48/120) =======================
q-point: ( 0.36  0.29  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.855   (  -3.422    5.927  -11.995)   13.810
   1.974   (  -2.393    4.145    3.272)    5.797
   2.011   (   2.877   -4.983  -12.750)   13.988
   2.558   (  -2.765    4.790   -1.764)    5.805
   2.990   (  -1.177    2.038    5.902)    6.354
   3.072   (   2.311   -4.002    4.940)    6.765
   8.532   (  -0.471    0.815   -1.150)    1.487
   8.554   (   0.305   -0.529   -1.794)    1.895
   8.884   (  -0.502    0.869    0.636)    1.188
   9.017   (  -0.038    0.067   -0.439)    0.445
  10.193   (   0.594   -1.028    6.918)    7.019
  10.245   (   2.201   -3.812    6.014)    7.453
======================= Grid point 422 (49/120) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 120
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.524   (   0.000    0.000   12.842)   12.842
   0.524   (   0.000    0.000   12.842)   12.842
   1.197   (   0.000   -0.000   -4.420)    4.420
   1.197   (  -0.000    0.000   -4.420)    4.420
   1.212   (   0.000   -0.000   30.864)   30.864
   3.065   (  -0.000    0.000  -14.983)   14.983
   8.851   (   0.000   -0.000    0.832)    0.832
   8.851   (   0.000    0.000    0.832)    0.832
   8.924   (   0.000    0.000   -0.695)    0.695
   8.924   (   0.000   -0.000   -0.695)    0.695
  10.290   (   0.000    0.000    8.457)    8.457
  11.015   (   0.000    0.000   -6.746)    6.746
======================= Grid point 423 (50/120) =======================
q-point: ( 0.07  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.590   (   5.364    3.097   11.580)   13.133
   0.708   (   8.006    4.622   13.163)   16.085
   1.231   (   2.913    1.682   -4.218)    5.395
   1.251   (   7.731    4.463   24.974)   26.521
   1.352   (  12.273    7.086   -2.101)   14.326
   3.049   (  -1.329   -0.767  -14.747)   14.827
   8.784   (  -6.271   -3.620    3.643)    8.106
   8.818   (  -1.853   -1.070    1.369)    2.540
   8.897   (  -0.082   -0.047    0.701)    0.707
   8.944   (   0.729    0.421   -1.249)    1.507
  10.288   (  -0.017   -0.010    9.544)    9.544
  11.074   (   3.556    2.053   -9.340)   10.203
======================= Grid point 424 (51/120) =======================
q-point: ( 0.14  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.747   (   7.766    4.484    9.559)   13.107
   0.796   (   3.079    1.778   18.578)   18.916
   1.327   (   5.257    3.035   -3.916)    7.224
   1.600   (  18.733   10.815   10.506)   24.047
   1.698   (  16.629    9.601   -0.662)   19.213
   3.012   (  -1.613   -0.931  -13.560)   13.688
   8.644   (  -5.046   -2.913    5.764)    8.196
   8.788   (  -0.607   -0.351    1.579)    1.728
   8.889   (  -1.443   -0.833   -1.011)    1.949
   8.937   (  -1.369   -0.790   -1.509)    2.185
  10.303   (   1.621    0.936   12.109)   12.253
  11.131   (   1.233    0.712   -8.571)    8.689
======================= Grid point 425 (52/120) =======================
q-point: ( 0.21  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.916   (   7.003    4.043   11.536)   14.088
   0.927   (   7.530    4.348    8.194)   11.948
   1.462   (   6.192    3.575   -3.842)    8.117
   1.979   (  12.567    7.255   -5.681)   15.583
   2.101   (  18.211   10.514    8.676)   22.747
   2.993   (   0.472    0.273  -10.139)   10.153
   8.557   (  -2.686   -1.551    4.217)    5.235
   8.785   (   0.009    0.005    1.099)    1.099
   8.822   (  -4.147   -2.394   -2.033)    5.202
   8.890   (  -2.336   -1.349   -1.083)    2.907
  10.364   (   3.421    1.975   14.185)   14.725
  11.134   (  -0.814   -0.470   -5.088)    5.174
======================= Grid point 426 (53/120) =======================
q-point: ( 0.29  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.083   (   5.915    3.415    7.441)   10.100
   1.094   (   8.065    4.656    3.793)   10.055
   1.597   (   5.282    3.050   -3.894)    7.237
   2.159   (   3.112    1.797  -12.963)   13.452
   2.504   (  16.278    9.398    9.473)   21.049
   3.054   (   4.857    2.804   -4.401)    7.129
   8.515   (  -1.097   -0.633    2.363)    2.681
   8.717   (  -4.463   -2.577   -1.889)    5.489
   8.772   (  -1.215   -0.702    0.811)    1.621
   8.851   (  -0.832   -0.480   -0.834)    1.272
  10.449   (   3.758    2.170   14.729)   15.355
  11.096   (  -2.549   -1.472   -1.612)    3.356
======================= Grid point 427 (54/120) =======================
q-point: ( 0.36  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.197   (   3.902    2.253    7.037)    8.356
   1.275   (   7.475    4.316   -2.407)    8.961
   1.696   (   3.197    1.846   -3.853)    5.336
   2.149   (  -3.150   -1.818  -16.171)   16.575
   2.836   (  12.300    7.101    8.475)   16.539
   3.192   (   6.313    3.645   -0.228)    7.294
   8.502   (  -0.063   -0.036    0.874)    0.877
   8.629   (  -3.052   -1.762   -1.302)    3.757
   8.744   (  -0.917   -0.530    0.954)    1.425
   8.840   (  -0.402   -0.232   -1.029)    1.129
  10.536   (   3.936    2.272   13.858)   14.585
  11.015   (  -4.447   -2.568    1.615)    5.383
======================= Grid point 428 (55/120) =======================
q-point: ( 0.43  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.264   (   1.935    1.117    6.804)    7.162
   1.429   (   5.534    3.195   -6.877)    9.388
   1.746   (   1.223    0.706   -3.675)    3.937
   2.053   (  -4.289   -2.476  -16.965)   17.673
   3.059   (   6.809    3.931    7.564)   10.910
   3.314   (   3.848    2.222    1.622)    4.730
   8.508   (   0.470    0.272   -0.388)    0.667
   8.576   (  -1.612   -0.931   -0.590)    1.953
   8.731   (  -0.302   -0.174    0.965)    1.027
   8.832   (  -0.266   -0.154   -1.112)    1.154
  10.636   (   4.832    2.790   11.723)   12.983
  10.896   (  -5.767   -3.329    4.978)    8.314
======================= Grid point 429 (56/120) =======================
q-point: (-0.50  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 235
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.286   (   0.000    0.000    6.714)    6.714
   1.501   (   0.000    0.000   -8.608)    8.608
   1.759   (   0.000    0.000   -3.561)    3.561
   1.994   (  -0.000   -0.000  -17.014)   17.014
   3.139   (   0.000    0.000    7.155)    7.155
   3.359   (   0.000    0.000    2.190)    2.190
   8.516   (   0.000    0.000   -1.308)    1.308
   8.556   (  -0.000   -0.000    0.123)    0.123
   8.728   (   0.000    0.000    0.953)    0.953
   8.829   (   0.000    0.000   -1.136)    1.136
  10.752   (  -0.000   -0.000    9.615)    9.615
  10.767   (   0.000    0.000    7.468)    7.468
======================= Grid point 438 (57/120) =======================
q-point: ( 0.07  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.721   (   5.289    9.161    9.860)   14.461
   0.771   (   1.125    1.949   18.665)   18.800
   1.308   (   3.924    6.796   -3.360)    8.536
   1.507   (  10.319   17.872   16.906)   26.678
   1.561   (   8.100   14.030   -4.873)   16.918
   3.022   (  -1.002   -1.736  -14.041)   14.184
   8.690   (  -3.030   -5.247    6.352)    8.778
   8.774   (  -2.000   -3.463    0.349)    4.014
   8.851   (  -0.717   -1.242    1.003)    1.750
   9.000   (   1.440    2.495   -2.503)    3.816
  10.295   (   0.561    0.972   11.346)   11.401
  11.116   (   0.736    1.275   -9.112)    9.230
======================= Grid point 439 (58/120) =======================
q-point: ( 0.14  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.844   (   4.364    4.421   13.846)   15.176
   0.924   (   5.391   11.366    8.666)   15.275
   1.456   (   3.517    9.306   -3.900)   10.686
   1.848   (  10.950   10.977   -3.939)   15.997
   1.922   (  15.769   14.303    9.608)   23.357
   2.992   (  -0.523   -0.815  -11.774)   11.814
   8.592   (  -3.433   -2.421    5.135)    6.634
   8.718   (   0.841   -6.073   -0.224)    6.135
   8.832   (  -1.023   -0.312    0.762)    1.313
   9.017   (  -3.409    4.089   -2.881)    6.053
  10.331   (   2.076    1.722   13.545)   13.811
  11.115   (   0.101   -2.536   -6.086)    6.594
======================= Grid point 440 (59/120) =======================
q-point: ( 0.21  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.993   (   6.676    5.137    6.207)   10.464
   1.119   (   3.073   12.522    7.547)   14.940
   1.607   (   2.264   10.385   -3.847)   11.304
   2.086   (   7.455    3.115  -10.481)   13.234
   2.329   (  15.590   12.021    9.913)   22.041
   3.009   (   2.664    1.476   -6.970)    7.606
   8.529   (  -1.754   -1.193    2.980)    3.658
   8.670   (   0.844   -6.592   -1.194)    6.752
   8.810   (  -2.064    1.347    1.031)    2.671
   8.975   (  -6.174    4.998   -2.184)    8.239
  10.395   (   3.532    0.860   14.603)   15.048
  11.075   (  -0.023   -5.439   -2.582)    6.021
======================= Grid point 441 (60/120) =======================
q-point: ( 0.29  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.158   (   6.272    4.831    0.712)    7.949
   1.276   (   1.098   12.021    5.851)   13.414
   1.731   (   0.146    9.763   -3.981)   10.544
   2.157   (   0.422   -1.762  -14.801)   14.912
   2.687   (  13.569    8.717    8.901)   18.421
   3.114   (   6.177    2.964   -1.783)    7.079
   8.504   (  -0.428   -0.425    1.392)    1.517
   8.619   (  -0.169   -4.626   -1.472)    4.857
   8.786   (  -2.580    1.945    1.120)    3.419
   8.920   (  -5.341    4.540   -1.477)    7.163
  10.466   (   4.553   -0.573   14.293)   15.012
  11.005   (  -0.712   -7.377    0.510)    7.429
======================= Grid point 442 (61/120) =======================
q-point: ( 0.36  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.291   (   2.546    5.796    0.239)    6.335
   1.404   (   1.637    9.360    0.764)    9.532
   1.801   (  -2.168    7.797   -4.454)    9.238
   2.099   (  -3.525   -2.428  -15.956)   16.520
   2.960   (   9.640    5.170    7.910)   13.499
   3.249   (   5.818    1.850    0.997)    6.186
   8.502   (   0.355    0.036    0.176)    0.398
   8.576   (  -0.570   -2.662   -1.161)    2.960
   8.764   (  -2.221    2.186    1.114)    3.309
   8.883   (  -3.553    3.451   -1.200)    5.096
  10.546   (   5.954   -1.314   12.770)   14.150
  10.903   (  -2.017   -8.350    3.543)    9.293
======================= Grid point 443 (62/120) =======================
q-point: ( 0.43  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.356   (  -2.091    6.744    2.582)    7.519
   1.509   (   0.860    5.585   -4.905)    7.483
   1.815   (  -2.955    4.599   -5.881)    8.030
   2.024   (  -3.261    0.902  -14.689)   15.074
   3.111   (   3.855    1.265    7.204)    8.268
   3.332   (   2.774   -0.448    2.141)    3.533
   8.512   (   0.536    0.194   -1.028)    1.175
   8.547   (  -0.172   -1.665   -0.397)    1.720
   8.753   (  -1.591    2.188    1.092)    2.918
   8.866   (  -2.226    3.062   -1.125)    3.949
  10.649   (   7.545   -1.510   10.252)   12.818
  10.776   (  -2.769   -8.205    6.731)   10.968
======================= Grid point 453 (63/120) =======================
q-point: ( 0.14  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.960   (   4.033    6.985    7.098)   10.745
   1.169   (   7.071   12.247    7.207)   15.872
   1.660   (   5.946   10.299   -3.621)   12.432
   2.020   (   3.804    6.589   -9.871)   12.463
   2.244   (   9.679   16.764   10.594)   22.067
   2.993   (   0.765    1.325   -8.065)    8.209
   8.545   (  -1.278   -2.214    3.752)    4.540
   8.620   (  -1.940   -3.360   -1.531)    4.171
   8.840   (   0.350    0.606    1.310)    1.486
   9.066   (   0.402    0.696   -3.118)    3.221
  10.365   (   0.914    1.584   14.437)   14.552
  11.038   (  -2.812   -4.871   -2.634)    6.210
======================= Grid point 454 (64/120) =======================
q-point: ( 0.21  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.109   (   5.249    6.318    0.460)    8.227
   1.389   (   3.127   13.686    6.593)   15.510
   1.843   (   1.572   12.499   -3.648)   13.115
   2.105   (   2.847   -1.098  -13.537)   13.876
   2.569   (  10.989   11.793    9.166)   18.543
   3.054   (   3.700    2.987   -3.226)    5.746
   8.510   (  -0.916   -0.781    1.949)    2.291
   8.576   (   0.100   -2.914   -1.803)    3.428
   8.846   (  -1.244    1.880    1.438)    2.673
   9.049   (  -3.708    2.176   -2.470)    4.959
  10.395   (   2.107   -0.785   14.302)   14.477
  10.933   (  -1.510   -8.220    0.407)    8.368
======================= Grid point 455 (65/120) =======================
q-point: ( 0.29  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 788
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.261   (   5.377    5.410   -4.927)    9.081
   1.546   (  -0.920   14.062    5.970)   15.305
   1.941   (  -3.343    8.749   -6.906)   11.637
   2.116   (   0.397   -0.305  -12.661)   12.670
   2.847   (   9.866    7.228    8.101)   14.670
   3.165   (   5.509    1.995    0.287)    5.866
   8.498   (  -0.003   -0.178    0.730)    0.751
   8.555   (   0.057   -1.883   -1.787)    2.597
   8.835   (  -2.472    2.574    1.353)    3.817
   8.996   (  -4.774    2.772   -1.734)    5.787
  10.429   (   4.228   -2.799   13.230)   14.169
  10.823   (  -0.651  -10.246    2.964)   10.686
======================= Grid point 456 (66/120) =======================
q-point: ( 0.36  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.394   (   3.636    4.717   -8.207)   10.140
   1.643   (  -3.434   13.319    4.778)   14.561
   1.905   (  -4.212   -0.485  -12.786)   13.470
   2.133   (  -3.470    8.199   -7.398)   11.576
   3.041   (   6.315    2.945    7.270)   10.070
   3.262   (   4.428   -0.584    2.039)    4.910
   8.503   (   0.594    0.042   -0.242)    0.643
   8.537   (  -0.390   -1.231   -1.529)    2.001
   8.818   (  -2.517    2.875    1.280)    4.029
   8.951   (  -3.577    2.929   -1.279)    4.796
  10.491   (   6.780   -3.847   11.299)   13.727
  10.706   (  -0.719  -10.646    5.551)   12.028
======================= Grid point 457 (67/120) =======================
q-point: ( 0.43  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.455   (  -1.601    2.772   -8.733)    9.301
   1.687   (  -6.611   11.450    3.928)   13.793
   1.849   (   0.488   -0.846  -14.472)   14.505
   2.145   (  -5.794   10.036   -6.073)   13.083
   3.112   (   0.607   -1.051    6.895)    7.001
   3.298   (   1.676   -2.903    2.602)    4.244
   8.517   (  -0.061    0.106   -1.719)    1.723
   8.519   (   0.531   -0.920   -0.325)    1.111
   8.810   (  -1.856    3.216    1.276)    3.926
   8.934   (  -1.975    3.420   -1.137)    4.110
  10.572   (   4.269   -7.395    9.523)   12.790
  10.604   (   3.850   -6.668    7.509)   10.755
======================= Grid point 469 (68/120) =======================
q-point: ( 0.21  0.21  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.240   (   3.863    6.691   -5.040)    9.225
   1.636   (   6.101   10.568    5.847)   13.531
   2.049   (  -2.139   -3.704  -15.775)   16.345
   2.081   (   5.923   10.260   -3.050)   12.233
   2.792   (   6.014   10.417    8.104)   14.503
   3.122   (   2.092    3.623   -0.118)    4.185
   8.498   (  -0.232   -0.402    1.111)    1.204
   8.544   (  -0.406   -0.703   -2.173)    2.320
   8.874   (   0.464    0.803    1.456)    1.727
   9.060   (  -0.530   -0.917   -2.160)    2.405
  10.374   (  -0.540   -0.935   13.333)   13.377
  10.768   (  -4.591   -7.952    2.955)    9.646
======================= Grid point 470 (69/120) =======================
q-point: ( 0.29  0.21  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.374   (   4.216    5.892   -9.638)   12.057
   1.799   (   0.591    9.930    3.420)   10.519
   1.953   (  -3.778   -3.519  -14.221)   15.129
   2.256   (   1.662   10.555   -3.508)   11.246
   2.974   (   5.536    5.524    7.175)   10.614
   3.189   (   3.021    0.390    1.987)    3.637
   8.496   (   0.232   -0.005    0.168)    0.287
   8.534   (  -0.248   -0.499   -2.104)    2.176
   8.878   (  -1.250    1.602    1.323)    2.425
   9.025   (  -2.523    0.211   -1.574)    2.981
  10.370   (   2.492   -2.666   11.825)   12.375
  10.623   (  -2.244   -9.273    5.067)   10.802
======================= Grid point 471 (70/120) =======================
q-point: ( 0.36  0.21  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.487   (   1.802    4.792  -12.432)   13.445
   1.817   (  -0.959   -3.533  -15.199)   15.634
   1.910   (  -4.561    9.457    4.131)   11.283
   2.345   (  -3.623   10.688   -3.756)   11.894
   3.082   (   2.564    1.371    6.573)    7.188
   3.227   (   2.669   -2.722    3.104)    4.916
   8.504   (   0.546    0.090   -0.619)    0.830
   8.523   (  -0.502   -0.350   -1.891)    1.987
   8.870   (  -1.641    2.237    1.271)    3.052
   8.992   (  -1.821    1.221   -1.135)    2.468
  10.410   (   5.559   -3.844    9.939)   12.019
  10.507   (  -0.572   -8.604    7.081)   11.157
======================= Grid point 484 (71/120) =======================
q-point: ( 0.29  0.29  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.495   (   3.543    6.136  -13.310)   15.079
   1.806   (  -3.290   -5.699  -16.487)   17.752
   2.005   (   3.203    5.547    5.161)    8.226
   2.437   (   3.923    6.794   -3.346)    8.529
   3.066   (   2.052    3.553    6.275)    7.497
   3.182   (  -0.544   -0.943    3.493)    3.658
   8.500   (   0.239    0.415   -0.654)    0.810
   8.526   (  -0.273   -0.473   -2.030)    2.102
   8.898   (   0.272    0.472    1.209)    1.326
   9.015   (  -0.498   -0.863   -1.132)    1.508
  10.338   (  -0.181   -0.314   10.275)   10.281
  10.466   (  -3.632   -6.291    6.763)    9.925
======================= Grid point 485 (72/120) =======================
q-point: ( 0.36  0.29  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.584   (  -2.613    4.526  -15.075)   15.955
   1.719   (   2.774   -4.805  -16.778)   17.671
   2.063   (  -2.483    4.300    5.860)    7.681
   2.507   (  -2.696    4.670   -3.594)    6.480
   3.107   (  -0.815    1.411    5.760)    5.986
   3.166   (   1.686   -2.919    4.325)    5.484
   8.508   (  -0.150    0.259   -1.223)    1.259
   8.517   (   0.138   -0.240   -1.793)    1.815
   8.902   (  -0.475    0.822    1.170)    1.507
   9.004   (  -0.108    0.188   -0.890)    0.916
  10.352   (   0.931   -1.612    9.037)    9.227
  10.385   (   1.861   -3.224    7.981)    8.806
======================= Grid point 633 (73/120) =======================
q-point: ( 0.00  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 96
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.751   (   0.000    0.000   11.128)   11.128
   0.751   (   0.000    0.000   11.128)   11.128
   1.091   (   0.000    0.000   -6.730)    6.730
   1.091   (   0.000    0.000   -6.730)    6.730
   1.776   (  -0.000    0.000   28.591)   28.591
   2.725   (  -0.000   -0.000  -20.794)   20.794
   8.869   (  -0.000    0.000    1.046)    1.046
   8.869   (   0.000    0.000    1.046)    1.046
   8.909   (  -0.000    0.000   -0.949)    0.949
   8.909   (   0.000    0.000   -0.949)    0.949
  10.472   (   0.000    0.000   10.418)   10.418
  10.863   (   0.000   -0.000   -9.228)    9.228
======================= Grid point 634 (74/120) =======================
q-point: ( 0.07  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.801   (   4.193    2.421   10.598)   11.652
   0.921   (  11.298    6.523   10.075)   16.484
   1.129   (   3.221    1.860   -6.498)    7.487
   1.238   (  11.407    6.586   -6.488)   14.684
   1.791   (   2.636    1.522   27.621)   27.788
   2.715   (  -0.642   -0.371  -20.326)   20.340
   8.835   (  -3.218   -1.858    2.161)    4.298
   8.848   (  -1.275   -0.736    1.757)    2.292
   8.897   (  -0.769   -0.444   -0.774)    1.178
   8.916   (   0.125    0.072   -1.672)    1.678
  10.486   (   1.380    0.797   11.071)   11.185
  10.893   (   2.402    1.387   -9.946)   10.326
======================= Grid point 635 (75/120) =======================
q-point: ( 0.14  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.930   (   6.615    3.819    9.505)   12.194
   1.125   (   5.360    3.094   15.125)   16.343
   1.232   (   5.494    3.172   -6.062)    8.774
   1.541   (  13.406    7.740   -7.715)   17.296
   1.965   (  12.602    7.276   19.492)   24.325
   2.712   (   0.880    0.508  -17.868)   17.897
   8.737   (  -4.721   -2.726    4.234)    6.902
   8.823   (  -0.804   -0.464    2.048)    2.249
   8.861   (  -2.603   -1.503   -2.048)    3.637
   8.904   (  -1.216   -0.702   -1.998)    2.442
  10.541   (   3.328    1.921   12.489)   13.067
  10.959   (   2.997    1.730   -9.638)   10.240
======================= Grid point 636 (76/120) =======================
q-point: ( 0.21  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.089   (   6.826    3.941    8.603)   11.668
   1.204   (   2.531    1.462   17.220)   17.467
   1.370   (   6.217    3.590   -5.814)    9.238
   1.793   (   7.447    4.300  -12.382)   15.075
   2.315   (  16.417    9.478   12.203)   22.544
   2.782   (   5.526    3.190  -11.929)   13.529
   8.635   (  -3.907   -2.256    4.129)    6.116
   8.778   (  -4.316   -2.492   -2.799)    5.716
   8.809   (  -0.545   -0.314    1.435)    1.566
   8.866   (  -1.816   -1.048   -1.423)    2.534
  10.630   (   4.132    2.385   13.358)   14.184
  11.018   (   1.980    1.143   -7.258)    7.609
======================= Grid point 637 (77/120) =======================
q-point: ( 0.29  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.233   (   5.542    3.200    8.105)   10.326
   1.260   (   2.381    1.375   13.044)   13.331
   1.506   (   5.272    3.044   -5.774)    8.390
   1.874   (   0.051    0.029  -16.602)   16.602
   2.682   (  14.941    8.626    8.960)   19.440
   2.959   (   9.189    5.305   -5.679)   12.034
   8.562   (  -2.387   -1.378    2.860)    3.972
   8.679   (  -3.988   -2.302   -2.447)    5.214
   8.790   (  -1.142   -0.659    1.069)    1.697
   8.833   (  -0.934   -0.539   -1.085)    1.529
  10.718   (   3.351    1.935   13.050)   13.611
  11.044   (   0.142    0.082   -4.236)    4.240
======================= Grid point 638 (78/120) =======================
q-point: ( 0.36  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.314   (   2.265    1.307    6.769)    7.256
   1.340   (   3.668    2.118    7.811)    8.886
   1.605   (   3.318    1.916   -5.688)    6.858
   1.825   (  -3.682   -2.126  -17.353)   17.866
   2.987   (  11.289    6.518    7.178)   14.881
   3.169   (   8.348    4.820   -2.343)    9.920
   8.524   (  -1.066   -0.615    1.604)    2.022
   8.603   (  -2.534   -1.463   -1.741)    3.404
   8.765   (  -0.824   -0.476    1.273)    1.589
   8.817   (  -0.544   -0.314   -1.326)    1.467
  10.784   (   2.503    1.445   11.630)   11.984
  11.025   (  -1.725   -0.996   -1.139)    2.294
======================= Grid point 639 (79/120) =======================
q-point: ( 0.43  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.364   (   2.025    1.169    0.860)    2.492
   1.402   (   1.779    1.027    7.567)    7.840
   1.659   (   1.402    0.809   -5.465)    5.700
   1.731   (  -3.857   -2.227  -16.219)   16.820
   3.190   (   6.149    3.550    6.034)    9.317
   3.322   (   4.705    2.717   -0.828)    5.496
   8.509   (  -0.260   -0.150    0.717)    0.777
   8.561   (  -1.158   -0.669   -1.185)    1.787
   8.753   (  -0.318   -0.184    1.310)    1.361
   8.807   (  -0.298   -0.172   -1.414)    1.455
  10.840   (   2.483    1.434    9.131)    9.571
  10.970   (  -2.829   -1.634    2.245)    3.964
======================= Grid point 640 (80/120) =======================
q-point: (-0.50  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 228
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.392   (   0.000    0.000   -1.868)    1.868
   1.423   (   0.000    0.000    7.450)    7.450
   1.674   (   0.000    0.000   -5.328)    5.328
   1.679   (  -0.000   -0.000  -15.286)   15.286
   3.262   (   0.000    0.000    5.589)    5.589
   3.377   (   0.000    0.000   -0.349)    0.349
   8.507   (   0.000    0.000    0.354)    0.354
   8.548   (   0.000    0.000   -0.945)    0.945
   8.749   (   0.000    0.000    1.305)    1.305
   8.804   (   0.000    0.000   -1.434)    1.434
  10.899   (   0.000    0.000    5.710)    5.710
  10.906   (   0.000    0.000    5.944)    5.944
======================= Grid point 649 (81/120) =======================
q-point: ( 0.07  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.909   (   4.474    7.750    9.680)   13.183
   1.073   (   3.717    6.438   13.072)   15.038
   1.216   (   3.804    6.589   -5.908)    9.632
   1.447   (   7.406   12.828   -7.199)   16.470
   1.892   (   5.582    9.667   22.449)   25.071
   2.709   (   0.053    0.092  -18.860)   18.860
   8.784   (  -2.510   -4.347    3.538)    6.142
   8.789   (  -2.336   -4.046    1.145)    4.810
   8.863   (  -0.504   -0.873    0.482)    1.118
   8.948   (   0.860    1.490   -2.945)    3.410
  10.521   (   1.543    2.673   12.042)   12.431
  10.936   (   1.475    2.554   -9.876)   10.307
======================= Grid point 650 (82/120) =======================
q-point: ( 0.14  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.078   (   3.671   10.166    9.059)   14.102
   1.166   (   3.614    1.624   16.611)   17.077
   1.365   (   3.572    9.446   -5.667)   11.581
   1.695   (   8.597    6.699  -10.304)   14.998
   2.170   (  12.244   12.702   14.463)   22.813
   2.740   (   2.675    2.637  -14.524)   15.001
   8.670   (  -4.170   -3.763    3.041)    6.388
   8.722   (  -0.074   -6.929    0.307)    6.936
   8.852   (  -1.264    1.057    1.550)    2.262
   8.958   (  -2.580    3.220   -3.294)    5.280
  10.588   (   2.994    2.497   12.992)   13.564
  10.982   (   1.983    0.268   -7.971)    8.218
======================= Grid point 651 (83/120) =======================
q-point: ( 0.21  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.190   (   2.326    2.615   12.723)   13.195
   1.282   (   2.882    9.438   10.176)   14.175
   1.516   (   2.132   10.431   -5.737)   12.094
   1.841   (   3.887    0.529  -14.894)   15.402
   2.522   (  14.014   10.823   10.068)   20.370
   2.866   (   7.165    5.154   -8.077)   11.964
   8.580   (  -2.373   -2.866    2.421)    4.440
   8.652   (  -0.184   -6.270   -1.043)    6.359
   8.837   (  -2.394    2.162    1.806)    3.697
   8.928   (  -4.802    4.238   -2.641)    6.928
  10.663   (   3.686    0.597   13.078)   13.601
  11.005   (   1.916   -2.670   -5.093)    6.061
======================= Grid point 652 (84/120) =======================
q-point: ( 0.29  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.247   (   2.934    0.130    8.907)    9.379
   1.413   (  -0.162   11.373    7.943)   13.874
   1.634   (   0.056    8.668   -6.404)   10.777
   1.854   (  -1.164   -1.133  -16.285)   16.366
   2.848   (  12.559    7.824    7.827)   16.739
   3.064   (   8.949    4.971   -3.541)   10.833
   8.529   (  -0.995   -1.583    1.539)    2.422
   8.595   (  -0.364   -4.195   -1.352)    4.422
   8.813   (  -2.912    2.450    1.679)    4.159
   8.887   (  -4.508    3.970   -1.976)    6.324
  10.723   (   3.861   -1.542   12.169)   12.860
  10.993   (   1.235   -5.033   -2.210)    5.634
======================= Grid point 653 (85/120) =======================
q-point: ( 0.36  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.297   (   3.125   -0.559    3.350)    4.616
   1.501   (  -2.491   10.972    6.720)   13.106
   1.674   (  -1.792    1.585   -9.902)   10.187
   1.817   (  -3.950    4.532  -11.993)   13.415
   3.098   (   8.906    4.480    6.413)   11.854
   3.246   (   7.043    2.815   -1.333)    7.701
   8.508   (  -0.226   -0.721    0.748)    1.063
   8.556   (  -0.278   -2.457   -1.145)    2.725
   8.790   (  -2.451    2.456    1.566)    3.807
   8.856   (  -3.234    3.222   -1.673)    4.862
  10.771   (   4.313   -2.749   10.273)   11.475
  10.949   (   0.337   -6.321    0.784)    6.379
======================= Grid point 654 (86/120) =======================
q-point: ( 0.43  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.339   (   2.280   -1.131   -0.681)    2.634
   1.549   (  -4.209   10.381    5.717)   12.576
   1.628   (  -1.120   -2.953  -12.556)   12.947
   1.812   (  -5.144    9.203   -7.715)   13.065
   3.235   (   3.594    0.871    5.588)    6.701
   3.350   (   3.098   -0.007   -0.367)    3.119
   8.503   (   0.143   -0.396    0.181)    0.458
   8.536   (   0.256   -1.515   -0.862)    1.762
   8.778   (  -1.749    2.432    1.502)    3.351
   8.840   (  -2.138    2.969   -1.566)    3.979
  10.823   (   5.113   -3.087    7.440)    9.540
  10.886   (  -0.065   -6.529    4.111)    7.716
======================= Grid point 664 (87/120) =======================
q-point: ( 0.14  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.186   (   0.836    1.448   15.282)   15.373
   1.313   (   6.640   11.500    7.775)   15.388
   1.574   (   6.144   10.642   -5.373)   13.411
   1.790   (   1.626    2.815  -14.117)   14.487
   2.448   (   8.445   14.627   10.693)   19.990
   2.828   (   3.593    6.223   -9.287)   11.742
   8.599   (  -2.136   -3.699    1.658)    4.582
   8.604   (  -2.273   -3.936   -0.076)    4.546
   8.875   (   0.559    0.968    2.483)    2.724
   9.002   (   0.559    0.969   -3.600)    3.770
  10.630   (   0.908    1.573   12.851)   12.978
  10.966   (  -1.010   -1.749   -5.101)    5.486
======================= Grid point 665 (88/120) =======================
q-point: ( 0.21  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.216   (   1.375    0.847   10.389)   10.514
   1.522   (   3.001   12.593    6.996)   14.715
   1.742   (  -1.079    9.488   -8.418)   12.730
   1.832   (   3.740    0.780  -13.108)   13.653
   2.736   (  10.020   10.499    8.267)   16.703
   2.980   (   6.430    6.112   -4.646)   10.015
   8.534   (  -1.600   -1.581    0.727)    2.364
   8.557   (   0.015   -3.279   -0.418)    3.306
   8.883   (  -1.511    2.082    2.452)    3.554
   8.996   (  -2.881    2.217   -3.053)    4.747
  10.651   (   1.718   -1.688   12.057)   12.295
  10.919   (  -0.079   -5.502   -2.292)    5.961
======================= Grid point 666 (89/120) =======================
q-point: ( 0.29  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.251   (   2.074    0.300    4.396)    4.870
   1.647   (  -1.586    7.432   -0.764)    7.638
   1.738   (  -2.769    3.776   -7.666)    8.982
   1.940   (  -0.382   10.881   -6.753)   12.812
   2.989   (   9.182    6.215    6.646)   12.926
   3.151   (   7.058    3.592   -1.854)    8.133
   8.507   (  -0.487   -0.702    0.395)    0.941
   8.533   (   0.134   -2.078   -0.787)    2.226
   8.868   (  -2.780    2.660    2.075)    4.371
   8.957   (  -4.060    2.744   -2.318)    5.421
  10.662   (   3.249   -4.202   10.655)   11.905
  10.858   (   0.848   -7.994    0.198)    8.041
======================= Grid point 667 (90/120) =======================
q-point: ( 0.36  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.293   (   2.409    0.110   -1.137)    2.666
   1.607   (  -2.098   -3.801  -14.253)   14.900
   1.774   (  -5.030   12.973    5.783)   15.068
   2.014   (  -4.554   12.569   -5.934)   14.626
   3.164   (   5.926    2.165    5.544)    8.399
   3.278   (   5.030    0.355   -0.397)    5.058
   8.498   (  -0.021   -0.355   -0.057)    0.360
   8.517   (   0.084   -1.427   -0.764)    1.621
   8.847   (  -2.665    2.950    1.797)    4.363
   8.921   (  -3.256    2.957   -1.830)    4.764
  10.685   (   5.039   -5.396    8.528)   11.280
  10.791   (   1.297   -8.863    2.832)    9.394
======================= Grid point 668 (91/120) =======================
q-point: ( 0.43  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.320   (   0.504   -0.873   -3.708)    3.842
   1.561   (   1.692   -2.930  -14.465)   14.856
   1.795   (  -7.402   12.820    6.583)   16.201
   2.037   (  -7.081   12.265   -5.756)   15.288
   3.227   (   0.839   -1.453    5.101)    5.369
   3.324   (   1.415   -2.450    0.100)    2.831
   8.497   (   0.129   -0.223   -0.392)    0.469
   8.510   (   0.583   -1.010   -0.600)    1.312
   8.838   (  -1.913    3.313    1.720)    4.194
   8.907   (  -1.971    3.413   -1.676)    4.283
  10.727   (   4.248   -7.357    6.120)   10.470
  10.733   (   3.915   -6.781    5.390)    9.506
======================= Grid point 680 (92/120) =======================
q-point: ( 0.21  0.21  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.234   (   0.632    1.094    4.897)    5.058
   1.689   (  -0.437   -0.758  -10.380)   10.416
   1.800   (   3.817    6.610    0.043)    7.633
   2.005   (   5.971   10.342   -4.856)   12.891
   2.935   (   5.326    9.224    6.691)   12.578
   3.102   (   3.382    5.858   -2.123)    7.089
   8.512   (  -0.597   -1.034    0.318)    1.235
   8.516   (  -0.482   -0.835   -0.711)    1.198
   8.910   (   0.342    0.592    2.331)    2.429
   9.013   (  -0.227   -0.393   -2.709)    2.747
  10.608   (  -1.298   -2.248   10.652)   10.964
  10.803   (  -3.337   -5.779    0.179)    6.676
======================= Grid point 681 (93/120) =======================
q-point: ( 0.29  0.21  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.263   (   1.563    1.047   -1.008)    2.134
   1.606   (  -3.016   -5.259  -15.912)   17.028
   1.946   (   0.473   11.122    6.307)   12.794
   2.175   (   0.977   11.181   -5.020)   12.295
   3.095   (   5.086    4.444    5.388)    8.640
   3.205   (   3.836    1.840   -0.366)    4.270
   8.497   (  -0.256   -0.432   -0.041)    0.504
   8.506   (  -0.077   -0.715   -0.845)    1.109
   8.909   (  -1.391    1.413    1.939)    2.774
   8.991   (  -2.091    0.608   -2.066)    3.002
  10.574   (   1.372   -4.169    8.980)    9.996
  10.698   (  -0.975   -7.485    2.290)    7.888
======================= Grid point 682 (94/120) =======================
q-point: ( 0.36  0.21  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.304   (   1.396    1.099   -5.657)    5.930
   1.513   (  -0.874   -4.919  -14.765)   15.587
   2.025   (  -4.154   11.089    6.671)   13.591
   2.260   (  -3.987   11.081   -5.203)   12.874
   3.189   (   2.596    0.488    4.547)    5.259
   3.263   (   2.680   -1.639    0.641)    3.206
   8.493   (   0.155   -0.199   -0.474)    0.537
   8.498   (  -0.121   -0.578   -0.761)    0.963
   8.899   (  -1.626    2.097    1.700)    3.152
   8.966   (  -1.716    1.509   -1.615)    2.798
  10.576   (   4.122   -5.015    6.977)    9.530
  10.622   (   0.921   -7.436    4.422)    8.701
======================= Grid point 695 (95/120) =======================
q-point: ( 0.29  0.29  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.297   (   1.369    2.371   -6.324)    6.891
   1.496   (  -3.077   -5.329  -14.855)   16.079
   2.123   (   3.576    6.194    6.698)    9.798
   2.356   (   3.772    6.533   -5.079)    9.095
   3.166   (   1.507    2.610    4.158)    5.135
   3.225   (   0.178    0.308    1.031)    1.091
   8.490   (  -0.102   -0.176   -0.359)    0.413
   8.497   (  -0.214   -0.370   -0.970)    1.060
   8.925   (   0.181    0.313    1.631)    1.670
   8.989   (  -0.250   -0.432   -1.552)    1.630
  10.512   (  -1.053   -1.824    7.367)    7.662
  10.574   (  -2.757   -4.775    3.990)    6.806
======================= Grid point 696 (96/120) =======================
q-point: ( 0.36  0.29  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.338   (  -1.250    2.165   -9.532)    9.855
   1.423   (   2.119   -3.671  -13.177)   13.842
   2.185   (  -2.525    4.373    6.795)    8.466
   2.423   (  -2.635    4.564   -5.232)    7.426
   3.195   (  -0.255    0.441    3.426)    3.463
   3.225   (   1.024   -1.774    1.855)    2.764
   8.490   (   0.048   -0.083   -0.660)    0.667
   8.492   (   0.074   -0.127   -0.875)    0.887
   8.927   (  -0.419    0.726    1.497)    1.715
   8.984   (  -0.214    0.370   -1.296)    1.365
  10.502   (   1.217   -2.109    6.100)    6.568
  10.517   (   1.593   -2.759    5.250)    6.140
======================= Grid point 844 (97/120) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.942   (   0.000    0.000   -9.022)    9.022
   0.942   (   0.000   -0.000   -9.022)    9.022
   0.942   (  -0.000    0.000    9.022)    9.022
   0.942   (  -0.000    0.000    9.022)    9.022
   2.287   (   0.000    0.000  -25.284)   25.284
   2.287   (  -0.000    0.000   25.284)   25.284
   8.889   (   0.000   -0.000   -1.080)    1.080
   8.889   (   0.000    0.000   -1.080)    1.080
   8.889   (   0.000   -0.000    1.080)    1.080
   8.889   (  -0.000    0.000    1.080)    1.080
  10.673   (  -0.000    0.000  -10.585)   10.585
  10.673   (  -0.000    0.000   10.585)   10.585
======================= Grid point 846 (98/120) =======================
q-point: ( 0.07  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.985   (   3.604    2.081   -8.727)    9.669
   0.985   (   3.604    2.081    8.727)    9.669
   1.096   (  11.482    6.629   -8.415)   15.704
   1.096   (  11.482    6.629    8.415)   15.704
   2.288   (   0.549    0.317  -24.652)   24.660
   2.288   (   0.549    0.317   24.652)   24.660
   8.872   (  -1.651   -0.953   -1.727)    2.572
   8.872   (  -1.651   -0.953    1.727)    2.572
   8.882   (  -0.580   -0.335   -1.857)    1.975
   8.882   (  -0.580   -0.335    1.857)    1.975
  10.696   (   2.004    1.157  -10.872)   11.116
  10.696   (   2.004    1.157   10.872)   11.116
======================= Grid point 848 (99/120) =======================
q-point: ( 0.14  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.098   (   5.913    3.414   -8.080)   10.578
   1.098   (   5.913    3.414    8.080)   10.578
   1.367   (  10.178    5.876  -10.797)   15.959
   1.367   (  10.178    5.876   10.797)   15.959
   2.348   (   5.445    3.143  -20.280)   21.232
   2.348   (   5.445    3.143   20.280)   21.232
   8.809   (  -3.701   -2.137   -3.365)    5.439
   8.809   (  -3.701   -2.137    3.365)    5.439
   8.864   (  -1.019   -0.588   -2.193)    2.489
   8.864   (  -1.019   -0.588    2.193)    2.489
  10.764   (   3.764    2.173  -11.033)   11.858
  10.764   (   3.764    2.173   11.033)   11.858
======================= Grid point 850 (100/120) =======================
q-point: ( 0.21  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.243   (   6.403    3.697   -7.559)   10.573
   1.243   (   6.403    3.697    7.559)   10.573
   1.523   (   3.490    2.015  -15.902)   16.405
   1.523   (   3.490    2.015   15.902)   16.405
   2.550   (  11.497    6.638  -12.466)   18.211
   2.550   (  11.497    6.638   12.466)   18.211
   8.713   (  -4.327   -2.498   -3.966)    6.379
   8.713   (  -4.327   -2.498    3.966)    6.379
   8.838   (  -1.187   -0.685   -1.548)    2.068
   8.838   (  -1.187   -0.685    1.548)    2.068
  10.854   (   3.737    2.157  -10.223)   11.097
  10.854   (   3.737    2.157   10.223)   11.097
======================= Grid point 852 (101/120) =======================
q-point: ( 0.29  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.381   (   5.338    3.082   -7.329)    9.576
   1.381   (   5.338    3.082    7.329)    9.576
   1.552   (  -0.496   -0.286  -16.913)   16.922
   1.552   (  -0.496   -0.286   16.913)   16.922
   2.837   (  12.585    7.266   -7.334)   16.278
   2.837   (  12.585    7.266    7.334)   16.278
   8.623   (  -3.335   -1.925   -3.321)    5.085
   8.623   (  -3.335   -1.925    3.321)    5.085
   8.811   (  -1.044   -0.603   -1.166)    1.677
   8.811   (  -1.044   -0.603    1.166)    1.677
  10.925   (   2.215    1.279   -8.574)    8.947
  10.925   (   2.215    1.279    8.574)    8.947
======================= Grid point 854 (102/120) =======================
q-point: ( 0.36  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.483   (   3.474    2.005   -7.160)    8.207
   1.483   (   3.474    2.005    7.160)    8.207
   1.519   (  -2.025   -1.169  -14.263)   14.453
   1.519   (  -2.025   -1.169   14.263)   14.453
   3.101   (   9.971    5.757   -4.815)   12.480
   3.101   (   9.971    5.757    4.815)   12.480
   8.563   (  -1.886   -1.089   -2.345)    3.201
   8.563   (  -1.886   -1.089    2.345)    3.201
   8.791   (  -0.695   -0.401   -1.408)    1.620
   8.791   (  -0.695   -0.401    1.408)    1.620
  10.956   (   0.656    0.378   -6.333)    6.378
  10.956   (   0.656    0.378    6.333)    6.378
======================= Grid point 856 (103/120) =======================
q-point: ( 0.43  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.472   (  -1.720   -0.993  -10.320)   10.510
   1.472   (  -1.720   -0.993   10.320)   10.510
   1.541   (   1.596    0.921   -6.919)    7.161
   1.541   (   1.596    0.921    6.919)    7.161
   3.281   (   5.445    3.144   -3.513)    7.203
   3.281   (   5.445    3.144    3.513)    7.203
   8.533   (  -0.757   -0.437   -1.596)    1.820
   8.533   (  -0.757   -0.437    1.596)    1.820
   8.780   (  -0.318   -0.184   -1.476)    1.521
   8.780   (  -0.318   -0.184    1.476)    1.521
  10.962   (  -0.052   -0.030   -3.432)    3.432
  10.962   (  -0.052   -0.030    3.432)    3.432
======================= Grid point 858 (104/120) =======================
q-point: (-0.50  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.451   (   0.000    0.000   -8.053)    8.053
   1.451   (   0.000    0.000    8.053)    8.053
   1.559   (   0.000    0.000   -6.785)    6.785
   1.559   (   0.000    0.000    6.785)    6.785
   3.345   (   0.000    0.000   -3.072)    3.072
   3.345   (   0.000    0.000    3.072)    3.072
   8.525   (   0.000    0.000   -1.301)    1.301
   8.525   (   0.000    0.000    1.301)    1.301
   8.776   (   0.000    0.000   -1.485)    1.485
   8.776   (   0.000    0.000    1.485)    1.485
  10.960   (   0.000    0.000   -0.342)    0.342
  10.960   (   0.000    0.000    0.342)    0.342
======================= Grid point 876 (105/120) =======================
q-point: ( 0.07  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.079   (   3.969    6.875   -8.079)   11.327
   1.079   (   3.969    6.875    8.079)   11.327
   1.289   (   6.192   10.724   -9.539)   15.632
   1.289   (   6.192   10.724    9.539)   15.632
   2.319   (   2.040    3.533  -22.055)   22.429
   2.319   (   2.040    3.533   22.055)   22.429
   8.810   (  -2.568   -4.449   -0.630)    5.175
   8.810   (  -2.568   -4.449    0.630)    5.175
   8.893   (   0.353    0.612   -2.708)    2.799
   8.893   (   0.353    0.612    2.708)    2.799
  10.740   (   1.797    3.113  -10.990)   11.563
  10.740   (   1.797    3.113   10.990)   11.563
======================= Grid point 878 (106/120) =======================
q-point: ( 0.14  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.235   (   3.839    9.794   -7.479)   12.907
   1.235   (   3.839    9.794    7.479)   12.907
   1.465   (   5.290    3.342  -13.933)   15.273
   1.465   (   5.290    3.342   13.933)   15.273
   2.455   (   7.353    7.627  -15.381)   18.676
   2.455   (   7.353    7.627   15.381)   18.676
   8.712   (  -2.420   -5.791   -1.457)    6.443
   8.712   (  -2.420   -5.791    1.457)    6.443
   8.896   (  -1.776    2.347   -3.032)    4.226
   8.896   (  -1.776    2.347    3.032)    4.226
  10.810   (   3.021    2.039  -10.377)   10.998
  10.810   (   3.021    2.039   10.377)   10.998
======================= Grid point 880 (107/120) =======================
q-point: ( 0.21  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.385   (   1.872    9.299   -7.819)   12.292
   1.385   (   1.872    9.299    7.819)   12.292
   1.543   (   2.153    0.706  -15.894)   16.055
   1.543   (   2.153    0.706   15.894)   16.055
   2.704   (  11.180    8.453   -9.112)   16.718
   2.704   (  11.180    8.453    9.112)   16.718
   8.623   (  -1.574   -4.841   -2.056)    5.490
   8.623   (  -1.574   -4.841    2.056)    5.490
   8.878   (  -3.418    3.255   -2.539)    5.359
   8.878   (  -3.418    3.255    2.539)    5.359
  10.873   (   3.239   -0.547   -8.984)    9.565
  10.873   (   3.239   -0.547    8.984)    9.565
======================= Grid point 882 (108/120) =======================
q-point: ( 0.29  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.458   (   0.300    0.899  -11.113)   11.153
   1.458   (   0.300    0.899   11.113)   11.153
   1.589   (  -0.816    7.414  -10.211)   12.645
   1.589   (  -0.816    7.414   10.211)   12.645
   2.977   (  11.042    6.564   -5.714)   14.059
   2.977   (  11.042    6.564    5.714)   14.059
   8.563   (  -0.833   -3.022   -1.903)    3.667
   8.563   (  -0.833   -3.022    1.903)    3.667
   8.848   (  -3.613    3.210   -2.035)    5.244
   8.848   (  -3.613    3.210    2.035)    5.244
  10.908   (   2.819   -3.019   -7.110)    8.222
  10.908   (   2.819   -3.019    7.110)    8.222
======================= Grid point 884 (109/120) =======================
q-point: ( 0.36  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.439   (  -0.124   -3.332  -10.992)   11.487
   1.439   (  -0.124   -3.332   10.992)   11.487
   1.648   (  -3.355   10.654   -7.362)   13.378
   1.648   (  -3.355   10.654    7.362)   13.378
   3.197   (   8.045    3.671   -3.936)    9.680
   3.197   (   8.045    3.671    3.936)    9.680
   8.531   (  -0.318   -1.677   -1.441)    2.234
   8.531   (  -0.318   -1.677    1.441)    2.234
   8.822   (  -2.819    2.846   -1.796)    4.390
   8.822   (  -2.819    2.846    1.796)    4.390
  10.915   (   2.481   -4.424   -4.695)    6.912
  10.915   (   2.481   -4.424    4.695)    6.912
======================= Grid point 886 (110/120) =======================
q-point: ( 0.43  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.409   (   0.826   -3.342   -8.218)    8.910
   1.409   (   0.826   -3.342    8.218)    8.910
   1.680   (  -5.169   10.524   -6.852)   13.580
   1.680   (  -5.169   10.524    6.852)   13.580
   3.317   (   3.355    0.424   -3.067)    4.566
   3.317   (   3.355    0.424    3.067)    4.566
   8.516   (   0.197   -1.024   -1.057)    1.484
   8.516   (   0.197   -1.024    1.057)    1.484
   8.809   (  -1.953    2.725   -1.699)    3.758
   8.809   (  -1.953    2.725    1.699)    3.758
  10.911   (   2.587   -4.797   -1.673)    5.701
  10.911   (   2.587   -4.797    1.673)    5.701
======================= Grid point 906 (111/120) =======================
q-point: ( 0.14  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.447   (   5.245    9.084   -3.433)   11.037
   1.447   (   5.245    9.084    3.433)   11.037
   1.504   (   1.320    2.287  -12.983)   13.249
   1.504   (   1.320    2.287   12.983)   13.249
   2.645   (   6.363   11.022  -10.081)   16.236
   2.645   (   6.363   11.022   10.081)   16.236
   8.609   (  -2.425   -4.201   -0.251)    4.857
   8.609   (  -2.425   -4.201    0.251)    4.857
   8.933   (   0.639    1.106   -3.508)    3.734
   8.933   (   0.639    1.106    3.508)    3.734
  10.836   (   0.306    0.530   -8.837)    8.859
  10.836   (   0.306    0.530    8.837)    8.859
======================= Grid point 908 (112/120) =======================
q-point: ( 0.21  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.467   (  -1.269   -1.452  -14.899)   15.023
   1.467   (  -1.269   -1.452   14.899)   15.023
   1.670   (   3.594   11.742   -6.463)   13.877
   1.670   (   3.594   11.742    6.463)   13.877
   2.876   (   8.531    8.605   -6.454)   13.729
   2.876   (   8.531    8.605    6.454)   13.729
   8.547   (  -0.921   -2.648   -0.682)    2.885
   8.547   (  -0.921   -2.648    0.682)    2.885
   8.936   (  -2.112    2.210   -3.118)    4.367
   8.936   (  -2.112    2.210    3.118)    4.367
  10.833   (   1.046   -3.236   -7.044)    7.822
  10.833   (   1.046   -3.236    7.044)    7.822
======================= Grid point 910 (113/120) =======================
q-point: ( 0.29  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 492
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.420   (  -1.239   -3.226  -12.715)   13.177
   1.420   (  -1.239   -3.226   12.715)   13.177
   1.817   (  -1.099   13.143   -6.590)   14.744
   1.817   (  -1.099   13.143    6.590)   14.744
   3.093   (   8.254    4.971   -4.286)   10.546
   3.093   (   8.254    4.971    4.286)   10.546
   8.518   (  -0.248   -1.518   -0.769)    1.720
   8.518   (  -0.248   -1.518    0.769)    1.720
   8.911   (  -3.359    2.725   -2.469)    4.980
   8.911   (  -3.359    2.725    2.469)    4.980
  10.814   (   2.172   -5.939   -5.130)    8.143
  10.814   (   2.172   -5.939    5.130)    8.143
======================= Grid point 912 (114/120) =======================
q-point: ( 0.36  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.370   (  -0.316   -3.478   -9.118)    9.764
   1.370   (  -0.316   -3.478    9.118)    9.764
   1.897   (  -4.846   13.421   -6.539)   15.696
   1.897   (  -4.846   13.421    6.539)   15.696
   3.246   (   5.511    1.259   -3.028)    6.413
   3.246   (   5.511    1.259    3.028)    6.413
   8.504   (   0.000   -0.993   -0.609)    1.165
   8.504   (   0.000   -0.993    0.609)    1.165
   8.884   (  -2.935    2.980   -2.028)    4.648
   8.884   (  -2.935    2.980    2.028)    4.648
  10.795   (   3.239   -7.097   -2.787)    8.284
  10.795   (   3.239   -7.097    2.787)    8.284
======================= Grid point 914 (115/120) =======================
q-point: ( 0.43  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.349   (   1.519   -2.631   -6.899)    7.538
   1.349   (   1.519   -2.631    6.899)    7.538
   1.921   (  -7.364   12.755   -6.498)   16.098
   1.921   (  -7.364   12.755    6.498)   16.098
   3.301   (   1.133   -1.962   -2.571)    3.427
   3.301   (   1.133   -1.962    2.571)    3.427
   8.499   (   0.396   -0.685   -0.479)    0.925
   8.499   (   0.396   -0.685    0.479)    0.925
   8.873   (  -1.954    3.385   -1.894)    4.343
   8.873   (  -1.954    3.385    1.894)    4.343
  10.787   (   4.100   -7.101   -0.053)    8.199
  10.787   (   4.100   -7.101    0.053)    8.199
======================= Grid point 938 (116/120) =======================
q-point: ( 0.21  0.21 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.415   (  -1.871   -3.240  -13.437)   13.948
   1.415   (  -1.871   -3.240   13.437)   13.948
   1.900   (   5.991   10.377   -6.094)   13.443
   1.900   (   5.991   10.377    6.094)   13.443
   3.040   (   4.439    7.689   -4.427)    9.921
   3.040   (   4.439    7.689    4.427)    9.921
   8.512   (  -0.614   -1.064   -0.073)    1.230
   8.512   (  -0.614   -1.064    0.073)    1.230
   8.960   (   0.089    0.154   -2.827)    2.833
   8.960   (   0.089    0.154    2.827)    2.833
  10.759   (  -2.218   -3.841   -5.107)    6.764
  10.759   (  -2.218   -3.841    5.107)    6.764
======================= Grid point 940 (117/120) =======================
q-point: ( 0.29  0.21 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.350   (  -1.413   -3.501   -9.902)   10.598
   1.350   (  -1.413   -3.501    9.902)   10.598
   2.068   (   0.693   11.317   -6.229)   12.937
   2.068   (   0.693   11.317    6.229)   12.937
   3.174   (   4.511    3.176   -2.908)    6.236
   3.174   (   4.511    3.176    2.908)    6.236
   8.497   (  -0.224   -0.685   -0.210)    0.751
   8.497   (  -0.224   -0.685    0.210)    0.751
   8.949   (  -1.695    1.053   -2.228)    2.991
   8.949   (  -1.695    1.053    2.228)    2.991
  10.692   (   0.239   -5.777   -3.246)    6.631
  10.692   (   0.239   -5.777    3.246)    6.631
======================= Grid point 942 (118/120) =======================
q-point: ( 0.36  0.21 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.306   (   0.191   -2.791   -6.072)    6.685
   1.306   (   0.191   -2.791    6.072)    6.685
   2.150   (  -4.114   11.150   -6.345)   13.472
   2.150   (  -4.114   11.150    6.345)   13.472
   3.251   (   2.649   -0.566   -1.987)    3.360
   3.251   (   2.649   -0.566    1.987)    3.360
   8.490   (   0.013   -0.481   -0.152)    0.504
   8.490   (   0.013   -0.481    0.152)    0.504
   8.933   (  -1.655    1.837   -1.835)    3.079
   8.933   (  -1.655    1.837    1.835)    3.079
  10.656   (   2.543   -6.238   -1.192)    6.841
  10.656   (   2.543   -6.238    1.192)    6.841
======================= Grid point 968 (119/120) =======================
q-point: ( 0.29  0.29 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.290   (  -1.328   -2.300   -5.815)    6.393
   1.290   (  -1.328   -2.300    5.815)    6.393
   2.248   (   3.665    6.348   -6.328)    9.684
   2.248   (   3.665    6.348    6.328)    9.684
   3.220   (   0.860    1.490   -1.582)    2.337
   3.220   (   0.860    1.490    1.582)    2.337
   8.487   (  -0.210   -0.363   -0.122)    0.437
   8.487   (  -0.210   -0.363    0.122)    0.437
   8.958   (  -0.005   -0.009   -1.754)    1.754
   8.958   (  -0.005   -0.009    1.754)    1.754
  10.599   (  -1.905   -3.299   -1.618)    4.138
  10.599   (  -1.905   -3.299    1.618)    4.138
======================= Grid point 970 (120/120) =======================
q-point: ( 0.36  0.29 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.266   (   0.627   -1.086   -2.496)    2.793
   1.266   (   0.627   -1.086    2.496)    2.793
   2.312   (  -2.576    4.462   -6.433)    8.242
   2.312   (  -2.576    4.462    6.433)    8.242
   3.235   (   0.375   -0.649   -0.793)    1.091
   3.235   (   0.375   -0.649    0.793)    1.091
   8.483   (   0.085   -0.148   -0.024)    0.172
   8.483   (   0.085   -0.148    0.024)    0.172
   8.956   (  -0.327    0.567   -1.534)    1.668
   8.956   (  -0.327    0.567    1.534)    1.668
  10.565   (   1.411   -2.444   -0.368)    2.846
  10.565   (   1.411   -2.444    0.368)    2.846
=================== 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/18816
   10.0    112.506    112.506    169.531     -0.000     -0.000      0.000 3/18816
   20.0     77.347     77.347     94.581     -0.000     -0.000      0.000 3/18816
   30.0     59.472     59.472     72.714     -0.000      0.000      0.000 3/18816
   40.0     48.048     48.048     58.931     -0.000      0.000      0.000 3/18816
   50.0     40.108     40.108     49.232     -0.000      0.000      0.000 3/18816
   60.0     34.121     34.121     41.869     -0.000      0.000      0.000 3/18816
   70.0     29.359     29.359     36.000     -0.000      0.000      0.000 3/18816
   80.0     25.492     25.492     31.235     -0.000      0.000      0.000 3/18816
   90.0     22.341     22.341     27.354     -0.000      0.000      0.000 3/18816
  100.0     19.771     19.771     24.189     -0.000      0.000      0.000 3/18816
  110.0     17.666     17.666     21.599     -0.000      0.000      0.000 3/18816
  120.0     15.930     15.930     19.464     -0.000      0.000      0.000 3/18816
  130.0     14.485     14.485     17.690     -0.000      0.000      0.000 3/18816
  140.0     13.270     13.270     16.199     -0.000      0.000      0.000 3/18816
  150.0     12.239     12.239     14.934     -0.000      0.000      0.000 3/18816
  160.0     11.354     11.354     13.851     -0.000      0.000      0.000 3/18816
  170.0     10.588     10.588     12.913     -0.000      0.000      0.000 3/18816
  180.0      9.920      9.920     12.095     -0.000      0.000      0.000 3/18816
  190.0      9.331      9.331     11.376     -0.000      0.000      0.000 3/18816
  200.0      8.810      8.810     10.738     -0.000      0.000      0.000 3/18816
  210.0      8.344      8.344     10.170     -0.000      0.000      0.000 3/18816
  220.0      7.926      7.926      9.659     -0.000      0.000      0.000 3/18816
  230.0      7.549      7.549      9.199     -0.000      0.000      0.000 3/18816
  240.0      7.207      7.207      8.782     -0.000      0.000      0.000 3/18816
  250.0      6.895      6.895      8.401     -0.000      0.000      0.000 3/18816
  260.0      6.610      6.610      8.053     -0.000      0.000      0.000 3/18816
  270.0      6.348      6.348      7.733     -0.000      0.000      0.000 3/18816
  280.0      6.106      6.106      7.439     -0.000      0.000      0.000 3/18816
  290.0      5.882      5.882      7.166     -0.000      0.000      0.000 3/18816
  300.0      5.675      5.675      6.913     -0.000      0.000      0.000 3/18816
  310.0      5.482      5.482      6.678     -0.000      0.000      0.000 3/18816
  320.0      5.302      5.302      6.459     -0.000      0.000      0.000 3/18816
  330.0      5.134      5.134      6.254     -0.000      0.000      0.000 3/18816
  340.0      4.976      4.976      6.061     -0.000      0.000      0.000 3/18816
  350.0      4.828      4.828      5.881     -0.000      0.000      0.000 3/18816
  360.0      4.689      4.689      5.711     -0.000      0.000      0.000 3/18816
  370.0      4.557      4.557      5.551     -0.000      0.000      0.000 3/18816
  380.0      4.433      4.433      5.400     -0.000      0.000      0.000 3/18816
  390.0      4.316      4.316      5.257     -0.000      0.000      0.000 3/18816
  400.0      4.204      4.204      5.121     -0.000      0.000      0.000 3/18816
  410.0      4.099      4.099      4.992     -0.000      0.000      0.000 3/18816
  420.0      3.998      3.998      4.870     -0.000      0.000      0.000 3/18816
  430.0      3.903      3.903      4.754     -0.000      0.000      0.000 3/18816
  440.0      3.812      3.812      4.643     -0.000      0.000      0.000 3/18816
  450.0      3.725      3.725      4.537     -0.000      0.000      0.000 3/18816
  460.0      3.643      3.643      4.436     -0.000      0.000      0.000 3/18816
  470.0      3.563      3.563      4.340     -0.000      0.000      0.000 3/18816
  480.0      3.488      3.488      4.248     -0.000      0.000      0.000 3/18816
  490.0      3.415      3.415      4.159     -0.000      0.000      0.000 3/18816
  500.0      3.346      3.346      4.075     -0.000      0.000      0.000 3/18816
  510.0      3.279      3.279      3.993     -0.000      0.000      0.000 3/18816
  520.0      3.215      3.215      3.915     -0.000      0.000      0.000 3/18816
  530.0      3.153      3.153      3.840     -0.000      0.000      0.000 3/18816
  540.0      3.094      3.094      3.768     -0.000      0.000      0.000 3/18816
  550.0      3.037      3.037      3.698     -0.000      0.000      0.000 3/18816
  560.0      2.982      2.982      3.631     -0.000      0.000      0.000 3/18816
  570.0      2.929      2.929      3.567     -0.000      0.000      0.000 3/18816
  580.0      2.878      2.878      3.505     -0.000      0.000      0.000 3/18816
  590.0      2.828      2.828      3.445     -0.000      0.000      0.000 3/18816
  600.0      2.781      2.781      3.386     -0.000      0.000      0.000 3/18816
  610.0      2.735      2.735      3.330     -0.000      0.000      0.000 3/18816
  620.0      2.690      2.690      3.276     -0.000      0.000      0.000 3/18816
  630.0      2.647      2.647      3.224     -0.000      0.000      0.000 3/18816
  640.0      2.605      2.605      3.173     -0.000      0.000      0.000 3/18816
  650.0      2.565      2.565      3.123     -0.000      0.000      0.000 3/18816
  660.0      2.526      2.526      3.076     -0.000      0.000      0.000 3/18816
  670.0      2.488      2.488      3.029     -0.000      0.000      0.000 3/18816
  680.0      2.451      2.451      2.984     -0.000      0.000      0.000 3/18816
  690.0      2.415      2.415      2.941     -0.000      0.000      0.000 3/18816
  700.0      2.380      2.380      2.899     -0.000      0.000      0.000 3/18816
  710.0      2.346      2.346      2.857     -0.000      0.000      0.000 3/18816
  720.0      2.314      2.314      2.817     -0.000      0.000      0.000 3/18816
  730.0      2.282      2.282      2.779     -0.000      0.000      0.000 3/18816
  740.0      2.251      2.251      2.741     -0.000      0.000      0.000 3/18816
  750.0      2.220      2.220      2.704     -0.000      0.000      0.000 3/18816
  760.0      2.191      2.191      2.668     -0.000      0.000      0.000 3/18816
  770.0      2.162      2.162      2.633     -0.000      0.000      0.000 3/18816
  780.0      2.135      2.135      2.599     -0.000      0.000      0.000 3/18816
  790.0      2.107      2.107      2.566     -0.000      0.000      0.000 3/18816
  800.0      2.081      2.081      2.534     -0.000      0.000      0.000 3/18816
  810.0      2.055      2.055      2.503     -0.000      0.000      0.000 3/18816
  820.0      2.030      2.030      2.472     -0.000      0.000      0.000 3/18816
  830.0      2.005      2.005      2.442     -0.000      0.000      0.000 3/18816
  840.0      1.981      1.981      2.413     -0.000      0.000      0.000 3/18816
  850.0      1.958      1.958      2.384     -0.000      0.000      0.000 3/18816
  860.0      1.935      1.935      2.357     -0.000      0.000      0.000 3/18816
  870.0      1.913      1.913      2.329     -0.000      0.000      0.000 3/18816
  880.0      1.891      1.891      2.303     -0.000      0.000      0.000 3/18816
  890.0      1.870      1.870      2.277     -0.000      0.000      0.000 3/18816
  900.0      1.849      1.849      2.252     -0.000      0.000      0.000 3/18816
  910.0      1.829      1.829      2.227     -0.000      0.000      0.000 3/18816
  920.0      1.809      1.809      2.202     -0.000      0.000      0.000 3/18816
  930.0      1.789      1.789      2.179     -0.000      0.000      0.000 3/18816
  940.0      1.770      1.770      2.155     -0.000      0.000      0.000 3/18816
  950.0      1.751      1.751      2.133     -0.000      0.000      0.000 3/18816
  960.0      1.733      1.733      2.110     -0.000      0.000      0.000 3/18816
  970.0      1.715      1.715      2.089     -0.000      0.000      0.000 3/18816
  980.0      1.698      1.698      2.067     -0.000      0.000      0.000 3/18816
  990.0      1.681      1.681      2.046     -0.000      0.000      0.000 3/18816
 1000.0      1.664      1.664      2.026     -0.000      0.000      0.000 3/18816

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

