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

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

-------------------------[time 2026-01-07 18:31:18]-------------------------
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    3.976815918746695    0.000000000000000    0.000000000000000
  b   -1.988407959373347    3.444023611808990    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.534142810000000
Atomic positions (fractional):
   *1 Zn  0.33333333333333  0.66666666666667  0.99946928815536  65.380
    2 Zn  0.66666666666667  0.33333333333333  0.49946928815536  65.380
   *3 Se  0.33333333333333  0.66666666666667  0.62553071184464  78.960
    4 Se  0.66666666666667  0.33333333333333  0.12553071184464  78.960
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.976815918746695    0.000000000000000    0.000000000000000
  b   -1.988407959373347    3.444023611808990    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.534142810000000
Atomic positions (fractional):
   *1 Zn  0.33333333333333  0.66666666666667  0.99946928815536  65.380 > 1
    2 Zn  0.66666666666667  0.33333333333333  0.49946928815536  65.380 > 2
   *3 Se  0.33333333333333  0.66666666666667  0.62553071184464  78.960 > 3
    4 Se  0.66666666666667  0.33333333333333  0.12553071184464  78.960 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   15.907263674986780    0.000000000000000    0.000000000000000
  b   -7.953631837493388   13.776094447235961    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.068285619999999
Atomic positions (fractional):
   *1 Zn  0.08333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    2 Zn  0.33333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    3 Zn  0.58333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    4 Zn  0.83333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    5 Zn  0.08333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    6 Zn  0.33333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    7 Zn  0.58333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    8 Zn  0.83333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    9 Zn  0.08333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   10 Zn  0.33333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   11 Zn  0.58333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   12 Zn  0.83333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   13 Zn  0.08333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   14 Zn  0.33333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   15 Zn  0.58333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   16 Zn  0.83333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   17 Zn  0.08333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   18 Zn  0.33333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   19 Zn  0.58333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   20 Zn  0.83333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   21 Zn  0.08333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   22 Zn  0.33333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   23 Zn  0.58333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   24 Zn  0.83333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   25 Zn  0.08333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   26 Zn  0.33333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   27 Zn  0.58333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   28 Zn  0.83333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   29 Zn  0.08333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   30 Zn  0.33333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   31 Zn  0.58333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   32 Zn  0.83333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   33 Zn  0.16666666666667  0.08333333333333  0.24973464407768  65.380 > 2
   34 Zn  0.41666666666667  0.08333333333333  0.24973464407768  65.380 > 2
   35 Zn  0.66666666666667  0.08333333333333  0.24973464407768  65.380 > 2
   36 Zn  0.91666666666667  0.08333333333333  0.24973464407768  65.380 > 2
   37 Zn  0.16666666666667  0.33333333333333  0.24973464407768  65.380 > 2
   38 Zn  0.41666666666667  0.33333333333333  0.24973464407768  65.380 > 2
   39 Zn  0.66666666666667  0.33333333333333  0.24973464407768  65.380 > 2
   40 Zn  0.91666666666667  0.33333333333333  0.24973464407768  65.380 > 2
   41 Zn  0.16666666666667  0.58333333333333  0.24973464407768  65.380 > 2
   42 Zn  0.41666666666667  0.58333333333333  0.24973464407768  65.380 > 2
   43 Zn  0.66666666666667  0.58333333333333  0.24973464407768  65.380 > 2
   44 Zn  0.91666666666667  0.58333333333333  0.24973464407768  65.380 > 2
   45 Zn  0.16666666666667  0.83333333333333  0.24973464407768  65.380 > 2
   46 Zn  0.41666666666667  0.83333333333333  0.24973464407768  65.380 > 2
   47 Zn  0.66666666666667  0.83333333333333  0.24973464407768  65.380 > 2
   48 Zn  0.91666666666667  0.83333333333333  0.24973464407768  65.380 > 2
   49 Zn  0.16666666666667  0.08333333333333  0.74973464407768  65.380 > 2
   50 Zn  0.41666666666667  0.08333333333333  0.74973464407768  65.380 > 2
   51 Zn  0.66666666666667  0.08333333333333  0.74973464407768  65.380 > 2
   52 Zn  0.91666666666667  0.08333333333333  0.74973464407768  65.380 > 2
   53 Zn  0.16666666666667  0.33333333333333  0.74973464407768  65.380 > 2
   54 Zn  0.41666666666667  0.33333333333333  0.74973464407768  65.380 > 2
   55 Zn  0.66666666666667  0.33333333333333  0.74973464407768  65.380 > 2
   56 Zn  0.91666666666667  0.33333333333333  0.74973464407768  65.380 > 2
   57 Zn  0.16666666666667  0.58333333333333  0.74973464407768  65.380 > 2
   58 Zn  0.41666666666667  0.58333333333333  0.74973464407768  65.380 > 2
   59 Zn  0.66666666666667  0.58333333333333  0.74973464407768  65.380 > 2
   60 Zn  0.91666666666667  0.58333333333333  0.74973464407768  65.380 > 2
   61 Zn  0.16666666666667  0.83333333333333  0.74973464407768  65.380 > 2
   62 Zn  0.41666666666667  0.83333333333333  0.74973464407768  65.380 > 2
   63 Zn  0.66666666666667  0.83333333333333  0.74973464407768  65.380 > 2
   64 Zn  0.91666666666667  0.83333333333333  0.74973464407768  65.380 > 2
  *65 Se  0.08333333333333  0.16666666666667  0.31276535592232  78.960 > 3
   66 Se  0.33333333333333  0.16666666666667  0.31276535592232  78.960 > 3
   67 Se  0.58333333333333  0.16666666666667  0.31276535592232  78.960 > 3
   68 Se  0.83333333333333  0.16666666666667  0.31276535592232  78.960 > 3
   69 Se  0.08333333333333  0.41666666666667  0.31276535592232  78.960 > 3
   70 Se  0.33333333333333  0.41666666666667  0.31276535592232  78.960 > 3
   71 Se  0.58333333333333  0.41666666666667  0.31276535592232  78.960 > 3
   72 Se  0.83333333333333  0.41666666666667  0.31276535592232  78.960 > 3
   73 Se  0.08333333333333  0.66666666666667  0.31276535592232  78.960 > 3
   74 Se  0.33333333333333  0.66666666666667  0.31276535592232  78.960 > 3
   75 Se  0.58333333333333  0.66666666666667  0.31276535592232  78.960 > 3
   76 Se  0.83333333333333  0.66666666666667  0.31276535592232  78.960 > 3
   77 Se  0.08333333333333  0.91666666666667  0.31276535592232  78.960 > 3
   78 Se  0.33333333333333  0.91666666666667  0.31276535592232  78.960 > 3
   79 Se  0.58333333333333  0.91666666666667  0.31276535592232  78.960 > 3
   80 Se  0.83333333333333  0.91666666666667  0.31276535592232  78.960 > 3
   81 Se  0.08333333333333  0.16666666666667  0.81276535592232  78.960 > 3
   82 Se  0.33333333333333  0.16666666666667  0.81276535592232  78.960 > 3
   83 Se  0.58333333333333  0.16666666666667  0.81276535592232  78.960 > 3
   84 Se  0.83333333333333  0.16666666666667  0.81276535592232  78.960 > 3
   85 Se  0.08333333333333  0.41666666666667  0.81276535592232  78.960 > 3
   86 Se  0.33333333333333  0.41666666666667  0.81276535592232  78.960 > 3
   87 Se  0.58333333333333  0.41666666666667  0.81276535592232  78.960 > 3
   88 Se  0.83333333333333  0.41666666666667  0.81276535592232  78.960 > 3
   89 Se  0.08333333333333  0.66666666666667  0.81276535592232  78.960 > 3
   90 Se  0.33333333333333  0.66666666666667  0.81276535592232  78.960 > 3
   91 Se  0.58333333333333  0.66666666666667  0.81276535592232  78.960 > 3
   92 Se  0.83333333333333  0.66666666666667  0.81276535592232  78.960 > 3
   93 Se  0.08333333333333  0.91666666666667  0.81276535592232  78.960 > 3
   94 Se  0.33333333333333  0.91666666666667  0.81276535592232  78.960 > 3
   95 Se  0.58333333333333  0.91666666666667  0.81276535592232  78.960 > 3
   96 Se  0.83333333333333  0.91666666666667  0.81276535592232  78.960 > 3
   97 Se  0.16666666666667  0.08333333333333  0.06276535592232  78.960 > 4
   98 Se  0.41666666666667  0.08333333333333  0.06276535592232  78.960 > 4
   99 Se  0.66666666666667  0.08333333333333  0.06276535592232  78.960 > 4
  100 Se  0.91666666666667  0.08333333333333  0.06276535592232  78.960 > 4
  101 Se  0.16666666666667  0.33333333333333  0.06276535592232  78.960 > 4
  102 Se  0.41666666666667  0.33333333333333  0.06276535592232  78.960 > 4
  103 Se  0.66666666666667  0.33333333333333  0.06276535592232  78.960 > 4
  104 Se  0.91666666666667  0.33333333333333  0.06276535592232  78.960 > 4
  105 Se  0.16666666666667  0.58333333333333  0.06276535592232  78.960 > 4
  106 Se  0.41666666666667  0.58333333333333  0.06276535592232  78.960 > 4
  107 Se  0.66666666666667  0.58333333333333  0.06276535592232  78.960 > 4
  108 Se  0.91666666666667  0.58333333333333  0.06276535592232  78.960 > 4
  109 Se  0.16666666666667  0.83333333333333  0.06276535592232  78.960 > 4
  110 Se  0.41666666666667  0.83333333333333  0.06276535592232  78.960 > 4
  111 Se  0.66666666666667  0.83333333333333  0.06276535592232  78.960 > 4
  112 Se  0.91666666666667  0.83333333333333  0.06276535592232  78.960 > 4
  113 Se  0.16666666666667  0.08333333333333  0.56276535592232  78.960 > 4
  114 Se  0.41666666666667  0.08333333333333  0.56276535592232  78.960 > 4
  115 Se  0.66666666666667  0.08333333333333  0.56276535592232  78.960 > 4
  116 Se  0.91666666666667  0.08333333333333  0.56276535592232  78.960 > 4
  117 Se  0.16666666666667  0.33333333333333  0.56276535592232  78.960 > 4
  118 Se  0.41666666666667  0.33333333333333  0.56276535592232  78.960 > 4
  119 Se  0.66666666666667  0.33333333333333  0.56276535592232  78.960 > 4
  120 Se  0.91666666666667  0.33333333333333  0.56276535592232  78.960 > 4
  121 Se  0.16666666666667  0.58333333333333  0.56276535592232  78.960 > 4
  122 Se  0.41666666666667  0.58333333333333  0.56276535592232  78.960 > 4
  123 Se  0.66666666666667  0.58333333333333  0.56276535592232  78.960 > 4
  124 Se  0.91666666666667  0.58333333333333  0.56276535592232  78.960 > 4
  125 Se  0.16666666666667  0.83333333333333  0.56276535592232  78.960 > 4
  126 Se  0.41666666666667  0.83333333333333  0.56276535592232  78.960 > 4
  127 Se  0.66666666666667  0.83333333333333  0.56276535592232  78.960 > 4
  128 Se  0.91666666666667  0.83333333333333  0.56276535592232  78.960 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            8.7723920    0.0000000    0.0000000
            0.0000000    8.7723920    0.0000000
            0.0000000    0.0000000    8.6626588
-------------------------- Born effective charges --------------------------
    1 Zn    2.0392631    0.0000000    0.0000000
            0.0000000    2.0392631    0.0000000
            0.0000000    0.0000000    2.1517870
    2 Zn    2.0392631    0.0000000    0.0000000
            0.0000000    2.0392631    0.0000000
            0.0000000    0.0000000    2.1517870
    3 Se   -2.0392631   -0.0000000    0.0000000
            0.0000000   -2.0392631    0.0000000
            0.0000000    0.0000000   -2.1517870
    4 Se   -2.0392631   -0.0000000    0.0000000
            0.0000000   -2.0392631    0.0000000
            0.0000000    0.0000000   -2.1517870
----------------------------------------------------------------------------
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: 188
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.156
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.168
--------------------------------- 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: 188
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-07 18:31:23]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 18:31:23]-------------------------
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    3.976815918746695    0.000000000000000    0.000000000000000
  b   -1.988407959373347    3.444023611808990    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.534142810000000
Atomic positions (fractional):
    1 Zn  0.33333333333333  0.66666666666667  0.99946928815536  65.380
    2 Zn  0.66666666666667  0.33333333333333  0.49946928815536  65.380
    3 Se  0.33333333333333  0.66666666666667  0.62553071184464  78.960
    4 Se  0.66666666666667  0.33333333333333  0.12553071184464  78.960
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.907263674986780    0.000000000000000    0.000000000000000
  b   -7.953631837493388   13.776094447235961    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.068285619999999
Atomic positions (fractional):
    1 Zn  0.08333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    2 Zn  0.33333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    3 Zn  0.58333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    4 Zn  0.83333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    5 Zn  0.08333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    6 Zn  0.33333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    7 Zn  0.58333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    8 Zn  0.83333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    9 Zn  0.08333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   10 Zn  0.33333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   11 Zn  0.58333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   12 Zn  0.83333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   13 Zn  0.08333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   14 Zn  0.33333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   15 Zn  0.58333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   16 Zn  0.83333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   17 Zn  0.08333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   18 Zn  0.33333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   19 Zn  0.58333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   20 Zn  0.83333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   21 Zn  0.08333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   22 Zn  0.33333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   23 Zn  0.58333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   24 Zn  0.83333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   25 Zn  0.08333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   26 Zn  0.33333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   27 Zn  0.58333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   28 Zn  0.83333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   29 Zn  0.08333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   30 Zn  0.33333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   31 Zn  0.58333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   32 Zn  0.83333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   33 Zn  0.16666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   34 Zn  0.41666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   35 Zn  0.66666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   36 Zn  0.91666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   37 Zn  0.16666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   38 Zn  0.41666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   39 Zn  0.66666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   40 Zn  0.91666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   41 Zn  0.16666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   42 Zn  0.41666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   43 Zn  0.66666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   44 Zn  0.91666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   45 Zn  0.16666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   46 Zn  0.41666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   47 Zn  0.66666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   48 Zn  0.91666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   49 Zn  0.16666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   50 Zn  0.41666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   51 Zn  0.66666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   52 Zn  0.91666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   53 Zn  0.16666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   54 Zn  0.41666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   55 Zn  0.66666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   56 Zn  0.91666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   57 Zn  0.16666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   58 Zn  0.41666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   59 Zn  0.66666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   60 Zn  0.91666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   61 Zn  0.16666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   62 Zn  0.41666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   63 Zn  0.66666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   64 Zn  0.91666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   65 Se  0.08333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   66 Se  0.33333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   67 Se  0.58333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   68 Se  0.83333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   69 Se  0.08333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   70 Se  0.33333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   71 Se  0.58333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   72 Se  0.83333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   73 Se  0.08333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   74 Se  0.33333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   75 Se  0.58333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   76 Se  0.83333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   77 Se  0.08333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   78 Se  0.33333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   79 Se  0.58333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   80 Se  0.83333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   81 Se  0.08333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   82 Se  0.33333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   83 Se  0.58333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   84 Se  0.83333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   85 Se  0.08333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   86 Se  0.33333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   87 Se  0.58333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   88 Se  0.83333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   89 Se  0.08333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   90 Se  0.33333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   91 Se  0.58333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   92 Se  0.83333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   93 Se  0.08333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   94 Se  0.33333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   95 Se  0.58333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   96 Se  0.83333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   97 Se  0.16666666666667  0.08333333333333  0.06276535592232  78.960 > 97
   98 Se  0.41666666666667  0.08333333333333  0.06276535592232  78.960 > 97
   99 Se  0.66666666666667  0.08333333333333  0.06276535592232  78.960 > 97
  100 Se  0.91666666666667  0.08333333333333  0.06276535592232  78.960 > 97
  101 Se  0.16666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  102 Se  0.41666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  103 Se  0.66666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  104 Se  0.91666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  105 Se  0.16666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  106 Se  0.41666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  107 Se  0.66666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  108 Se  0.91666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  109 Se  0.16666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  110 Se  0.41666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  111 Se  0.66666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  112 Se  0.91666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  113 Se  0.16666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  114 Se  0.41666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  115 Se  0.66666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  116 Se  0.91666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  117 Se  0.16666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  118 Se  0.41666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  119 Se  0.66666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  120 Se  0.91666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  121 Se  0.16666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  122 Se  0.41666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  123 Se  0.66666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  124 Se  0.91666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  125 Se  0.16666666666667  0.83333333333333  0.56276535592232  78.960 > 97
  126 Se  0.41666666666667  0.83333333333333  0.56276535592232  78.960 > 97
  127 Se  0.66666666666667  0.83333333333333  0.56276535592232  78.960 > 97
  128 Se  0.91666666666667  0.83333333333333  0.56276535592232  78.960 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            8.7723920    0.0000000    0.0000000
            0.0000000    8.7723920    0.0000000
            0.0000000    0.0000000    8.6626588
-------------------------- Born effective charges --------------------------
    1 Zn    2.0392631    0.0000000    0.0000000
            0.0000000    2.0392631    0.0000000
            0.0000000    0.0000000    2.1517870
    2 Zn    2.0392631    0.0000000    0.0000000
            0.0000000    2.0392631    0.0000000
            0.0000000    0.0000000    2.1517870
    3 Se   -2.0392631   -0.0000000    0.0000000
            0.0000000   -2.0392631    0.0000000
            0.0000000    0.0000000   -2.1517870
    4 Se   -2.0392631   -0.0000000    0.0000000
            0.0000000   -2.0392631    0.0000000
            0.0000000    0.0000000   -2.1517870
----------------------------------------------------------------------------
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.00000030 (xxz) 0.00000030 (xxz) 0.00000030 (xzx)
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-07 18:31:28]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 18:31:29]-------------------------
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    3.976815918746695    0.000000000000000    0.000000000000000
  b   -1.988407959373347    3.444023611808990    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.534142810000000
Atomic positions (fractional):
    1 Zn  0.33333333333333  0.66666666666667  0.99946928815536  65.380
    2 Zn  0.66666666666667  0.33333333333333  0.49946928815536  65.380
    3 Se  0.33333333333333  0.66666666666667  0.62553071184464  78.960
    4 Se  0.66666666666667  0.33333333333333  0.12553071184464  78.960
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.907263674986780    0.000000000000000    0.000000000000000
  b   -7.953631837493388   13.776094447235961    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.068285619999999
Atomic positions (fractional):
    1 Zn  0.08333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    2 Zn  0.33333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    3 Zn  0.58333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    4 Zn  0.83333333333333  0.16666666666667  0.49973464407768  65.380 > 1
    5 Zn  0.08333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    6 Zn  0.33333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    7 Zn  0.58333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    8 Zn  0.83333333333333  0.41666666666667  0.49973464407768  65.380 > 1
    9 Zn  0.08333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   10 Zn  0.33333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   11 Zn  0.58333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   12 Zn  0.83333333333333  0.66666666666667  0.49973464407768  65.380 > 1
   13 Zn  0.08333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   14 Zn  0.33333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   15 Zn  0.58333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   16 Zn  0.83333333333333  0.91666666666667  0.49973464407768  65.380 > 1
   17 Zn  0.08333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   18 Zn  0.33333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   19 Zn  0.58333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   20 Zn  0.83333333333333  0.16666666666667  0.99973464407768  65.380 > 1
   21 Zn  0.08333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   22 Zn  0.33333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   23 Zn  0.58333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   24 Zn  0.83333333333333  0.41666666666667  0.99973464407768  65.380 > 1
   25 Zn  0.08333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   26 Zn  0.33333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   27 Zn  0.58333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   28 Zn  0.83333333333333  0.66666666666667  0.99973464407768  65.380 > 1
   29 Zn  0.08333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   30 Zn  0.33333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   31 Zn  0.58333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   32 Zn  0.83333333333333  0.91666666666667  0.99973464407768  65.380 > 1
   33 Zn  0.16666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   34 Zn  0.41666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   35 Zn  0.66666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   36 Zn  0.91666666666667  0.08333333333333  0.24973464407768  65.380 > 33
   37 Zn  0.16666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   38 Zn  0.41666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   39 Zn  0.66666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   40 Zn  0.91666666666667  0.33333333333333  0.24973464407768  65.380 > 33
   41 Zn  0.16666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   42 Zn  0.41666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   43 Zn  0.66666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   44 Zn  0.91666666666667  0.58333333333333  0.24973464407768  65.380 > 33
   45 Zn  0.16666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   46 Zn  0.41666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   47 Zn  0.66666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   48 Zn  0.91666666666667  0.83333333333333  0.24973464407768  65.380 > 33
   49 Zn  0.16666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   50 Zn  0.41666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   51 Zn  0.66666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   52 Zn  0.91666666666667  0.08333333333333  0.74973464407768  65.380 > 33
   53 Zn  0.16666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   54 Zn  0.41666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   55 Zn  0.66666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   56 Zn  0.91666666666667  0.33333333333333  0.74973464407768  65.380 > 33
   57 Zn  0.16666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   58 Zn  0.41666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   59 Zn  0.66666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   60 Zn  0.91666666666667  0.58333333333333  0.74973464407768  65.380 > 33
   61 Zn  0.16666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   62 Zn  0.41666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   63 Zn  0.66666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   64 Zn  0.91666666666667  0.83333333333333  0.74973464407768  65.380 > 33
   65 Se  0.08333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   66 Se  0.33333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   67 Se  0.58333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   68 Se  0.83333333333333  0.16666666666667  0.31276535592232  78.960 > 65
   69 Se  0.08333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   70 Se  0.33333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   71 Se  0.58333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   72 Se  0.83333333333333  0.41666666666667  0.31276535592232  78.960 > 65
   73 Se  0.08333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   74 Se  0.33333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   75 Se  0.58333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   76 Se  0.83333333333333  0.66666666666667  0.31276535592232  78.960 > 65
   77 Se  0.08333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   78 Se  0.33333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   79 Se  0.58333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   80 Se  0.83333333333333  0.91666666666667  0.31276535592232  78.960 > 65
   81 Se  0.08333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   82 Se  0.33333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   83 Se  0.58333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   84 Se  0.83333333333333  0.16666666666667  0.81276535592232  78.960 > 65
   85 Se  0.08333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   86 Se  0.33333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   87 Se  0.58333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   88 Se  0.83333333333333  0.41666666666667  0.81276535592232  78.960 > 65
   89 Se  0.08333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   90 Se  0.33333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   91 Se  0.58333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   92 Se  0.83333333333333  0.66666666666667  0.81276535592232  78.960 > 65
   93 Se  0.08333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   94 Se  0.33333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   95 Se  0.58333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   96 Se  0.83333333333333  0.91666666666667  0.81276535592232  78.960 > 65
   97 Se  0.16666666666667  0.08333333333333  0.06276535592232  78.960 > 97
   98 Se  0.41666666666667  0.08333333333333  0.06276535592232  78.960 > 97
   99 Se  0.66666666666667  0.08333333333333  0.06276535592232  78.960 > 97
  100 Se  0.91666666666667  0.08333333333333  0.06276535592232  78.960 > 97
  101 Se  0.16666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  102 Se  0.41666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  103 Se  0.66666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  104 Se  0.91666666666667  0.33333333333333  0.06276535592232  78.960 > 97
  105 Se  0.16666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  106 Se  0.41666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  107 Se  0.66666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  108 Se  0.91666666666667  0.58333333333333  0.06276535592232  78.960 > 97
  109 Se  0.16666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  110 Se  0.41666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  111 Se  0.66666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  112 Se  0.91666666666667  0.83333333333333  0.06276535592232  78.960 > 97
  113 Se  0.16666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  114 Se  0.41666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  115 Se  0.66666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  116 Se  0.91666666666667  0.08333333333333  0.56276535592232  78.960 > 97
  117 Se  0.16666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  118 Se  0.41666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  119 Se  0.66666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  120 Se  0.91666666666667  0.33333333333333  0.56276535592232  78.960 > 97
  121 Se  0.16666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  122 Se  0.41666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  123 Se  0.66666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  124 Se  0.91666666666667  0.58333333333333  0.56276535592232  78.960 > 97
  125 Se  0.16666666666667  0.83333333333333  0.56276535592232  78.960 > 97
  126 Se  0.41666666666667  0.83333333333333  0.56276535592232  78.960 > 97
  127 Se  0.66666666666667  0.83333333333333  0.56276535592232  78.960 > 97
  128 Se  0.91666666666667  0.83333333333333  0.56276535592232  78.960 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            8.7723920    0.0000000    0.0000000
            0.0000000    8.7723920    0.0000000
            0.0000000    0.0000000    8.6626588
-------------------------- Born effective charges --------------------------
    1 Zn    2.0392631    0.0000000    0.0000000
            0.0000000    2.0392631    0.0000000
            0.0000000    0.0000000    2.1517870
    2 Zn    2.0392631    0.0000000    0.0000000
            0.0000000    2.0392631    0.0000000
            0.0000000    0.0000000    2.1517870
    3 Se   -2.0392631   -0.0000000    0.0000000
            0.0000000   -2.0392631    0.0000000
            0.0000000    0.0000000   -2.1517870
    4 Se   -2.0392631   -0.0000000    0.0000000
            0.0000000   -2.0392631    0.0000000
            0.0000000    0.0000000   -2.1517870
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000030 (xxz) 0.00000030 (xxz) 0.00000030 (xzx)
Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 15 15 8 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.93, Number of G-points: 313, Lambda: 0.29
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/135) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   1.519   (   0.000    0.000    0.000)    0.000
   1.519   (   0.000    0.000    0.000)    0.000
   5.103   (   0.000    0.000   -0.000)    0.000
   6.243   (  -0.000    0.000   -0.000)    0.000
   6.243   (   0.000    0.000    0.000)    0.000
   6.282   (   0.000    0.000    0.000)    0.000
   6.373   (  -0.000   -0.000    0.000)    0.000
   6.373   (   0.000    0.000    0.000)    0.000
   6.430   (   0.000   -0.000    0.000)    0.000
======================= Grid point 1 (2/135) =======================
q-point: ( 0.07  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.384   (  17.456   10.078    0.000)   20.157
   0.384   (  16.889    9.751    0.000)   19.501
   0.863   (  37.787   21.817    0.000)   43.633
   1.560   (   3.551    2.050    0.000)    4.101
   1.732   (  17.056    9.847    0.000)   19.695
   5.039   (  -5.205   -3.005    0.000)    6.011
   6.238   (  -0.331   -0.191    0.000)    0.383
   6.248   (  -0.451   -0.260    0.000)    0.521
   6.279   (  -0.251   -0.145    0.000)    0.290
   6.378   (   0.318    0.184    0.000)    0.367
   6.456   (   2.827    1.632    0.000)    3.264
   7.326   (  -0.271   -0.156    0.000)    0.313
======================= Grid point 2 (3/135) =======================
q-point: ( 0.13  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.747   (  15.307    8.837    0.000)   17.675
   0.768   (  16.926    9.772    0.000)   19.545
   1.670   (   6.033    3.483    0.000)    6.966
   1.679   (  34.759   20.068    0.000)   40.137
   2.174   (  20.396   11.776    0.000)   23.552
   4.905   (  -5.830   -3.366    0.000)    6.731
   6.201   (  -3.673   -2.121    0.000)    4.241
   6.236   (   0.243    0.140    0.000)    0.281
   6.271   (  -0.346   -0.200    0.000)    0.399
   6.379   (  -0.327   -0.189    0.000)    0.378
   6.558   (   5.626    3.248    0.000)    6.496
   7.314   (  -0.904   -0.522    0.000)    1.044
======================= Grid point 3 (4/135) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.058   (  12.305    7.104    0.000)   14.209
   1.138   (  16.160    9.330    0.000)   18.660
   1.814   (   6.463    3.731    0.000)    7.462
   2.375   (  24.617   14.213    0.000)   28.425
   2.599   (  19.318   11.153    0.000)   22.306
   4.825   (  -0.359   -0.207    0.000)    0.414
   6.094   (  -5.736   -3.312    0.000)    6.624
   6.250   (   1.009    0.583    0.000)    1.166
   6.268   (   0.117    0.067    0.000)    0.135
   6.364   (  -0.959   -0.554    0.000)    1.108
   6.678   (   4.605    2.659    0.000)    5.318
   7.280   (  -2.339   -1.351    0.000)    2.701
======================= Grid point 4 (5/135) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.288   (   8.235    4.754    0.000)    9.509
   1.487   (  14.923    8.616    0.000)   17.231
   1.942   (   4.620    2.668    0.000)    5.335
   2.656   (   4.110    2.373    0.000)    4.745
   3.140   (  25.894   14.950    0.000)   29.899
   4.911   (   7.958    4.594    0.000)    9.189
   5.948   (  -7.297   -4.213    0.000)    8.426
   6.279   (   1.455    0.840    0.000)    1.680
   6.280   (   0.965    0.557    0.000)    1.114
   6.343   (  -0.764   -0.441    0.000)    0.882
   6.751   (   1.787    1.032    0.000)    2.063
   7.201   (  -4.825   -2.786    0.000)    5.571
======================= Grid point 5 (6/135) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.427   (   4.290    2.477    0.000)    4.954
   1.800   (  12.986    7.497    0.000)   14.995
   2.011   (   1.466    0.846    0.000)    1.693
   2.672   (  -1.411   -0.815    0.000)    1.630
   3.682   (  22.273   12.859    0.000)   25.719
   5.156   (  13.145    7.589    0.000)   15.178
   5.768   (  -8.784   -5.071    0.000)   10.142
   6.310   (   1.670    0.964    0.000)    1.928
   6.311   (   1.368    0.790    0.000)    1.580
   6.335   (   0.196    0.113    0.000)    0.227
   6.758   (  -1.037   -0.598    0.000)    1.197
   7.061   (  -7.697   -4.444    0.000)    8.888
======================= Grid point 6 (7/135) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.491   (   1.699    0.981    0.000)    1.961
   2.012   (  -0.988   -0.570    0.000)    1.141
   2.059   (   9.873    5.700    0.000)   11.401
   2.626   (  -2.336   -1.349    0.000)    2.698
   4.127   (  17.265    9.968    0.000)   19.936
   5.468   (  14.089    8.134    0.000)   16.268
   5.561   (  -9.425   -5.441    0.000)   10.883
   6.336   (   0.929    0.536    0.000)    1.073
   6.345   (   1.137    0.656    0.000)    1.312
   6.350   (   0.858    0.496    0.000)    0.991
   6.707   (  -3.596   -2.076    0.000)    4.153
   6.876   (  -7.887   -4.553    0.000)    9.107
======================= Grid point 7 (8/135) =======================
q-point: ( 0.47  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.513   (   0.449    0.259    0.000)    0.519
   1.986   (  -0.888   -0.513    0.000)    1.026
   2.224   (   4.294    2.479    0.000)    4.958
   2.579   (  -1.499   -0.866    0.000)    1.731
   4.423   (   7.970    4.601    0.000)    9.202
   5.379   (  -5.446   -3.144    0.000)    6.288
   5.739   (   8.310    4.798    0.000)    9.596
   6.340   (  -1.555   -0.898    0.000)    1.796
   6.354   (   0.531    0.306    0.000)    0.613
   6.363   (   0.272    0.157    0.000)    0.314
   6.612   (  -3.221   -1.860    0.000)    3.719
   6.759   (  -2.336   -1.349    0.000)    2.698
======================= Grid point 16 (9/135) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.665   (   9.818   17.005    0.000)   19.636
   0.696   (  10.824   18.747    0.000)   21.647
   1.432   (  18.284   31.669    0.000)   36.568
   1.696   (   6.238   10.805    0.000)   12.476
   2.009   (  10.752   18.622    0.000)   21.503
   4.941   (  -3.633   -6.292    0.000)    7.265
   6.219   (  -1.620   -2.806    0.000)    3.240
   6.247   (   0.053    0.092    0.000)    0.107
   6.278   (   0.089    0.155    0.000)    0.179
   6.368   (   0.045    0.078    0.000)    0.091
   6.533   (   2.999    5.194    0.000)    5.998
   7.312   (  -0.728   -1.262    0.000)    1.457
======================= Grid point 17 (10/135) =======================
q-point: ( 0.13  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.006   (  12.458   14.098    0.000)   18.813
   1.044   (   7.944   19.390    0.000)   20.955
   1.806   (   5.218   10.427    0.000)   11.660
   2.150   (  21.425   25.981    0.000)   33.676
   2.390   (  16.411   10.307    0.000)   19.379
   4.836   (  -2.294   -2.983    0.000)    3.763
   6.133   (  -3.897   -4.924    0.000)    6.279
   6.248   (   0.105    0.328    0.000)    0.345
   6.283   (  -0.428    1.672    0.000)    1.725
   6.381   (   0.132    1.107    0.000)    1.115
   6.637   (   3.640    4.178    0.000)    5.541
   7.277   (  -0.957   -3.332    0.000)    3.467
======================= Grid point 18 (11/135) =======================
q-point: ( 0.20  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.323   (   3.551   18.216    0.000)   18.559
   1.348   (  12.775   12.175    0.000)   17.647
   1.967   (   1.752   11.154    0.000)   11.291
   2.582   (   8.414    6.190    0.000)   10.446
   2.851   (  25.234   15.476    0.000)   29.601
   4.843   (   4.110    2.328    0.000)    4.723
   5.999   (  -5.372   -6.455    0.000)    8.398
   6.259   (   0.865    0.282    0.000)    0.910
   6.299   (  -0.261    2.895    0.000)    2.906
   6.384   (  -1.859    2.701    0.000)    3.279
   6.722   (   2.745    1.586    0.000)    3.171
   7.212   (  -1.835   -5.572    0.000)    5.866
======================= Grid point 19 (12/135) =======================
q-point: ( 0.27  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.508   (  -1.305   16.454    0.000)   16.506
   1.668   (  12.408    9.994    0.000)   15.933
   2.073   (  -1.582   10.140    0.000)   10.263
   2.673   (   0.921   -0.103    0.000)    0.926
   3.413   (  24.140   13.203    0.000)   27.515
   5.015   (  10.938    5.943    0.000)   12.448
   5.829   (  -6.539   -7.983    0.000)   10.319
   6.283   (   1.613    0.337    0.000)    1.648
   6.328   (  -0.060    3.382    0.000)    3.382
   6.377   (  -2.446    3.627    0.000)    4.375
   6.755   (   0.613   -0.605    0.000)    0.861
   7.101   (  -4.028   -7.234    0.000)    8.280
======================= Grid point 20 (13/135) =======================
q-point: ( 0.33  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.604   (  -4.805   14.522    0.000)   15.296
   1.950   (  10.564    8.018    0.000)   13.262
   2.109   (  -4.799    8.643    0.000)    9.886
   2.653   (  -1.978   -0.939    0.000)    2.189
   3.902   (  20.595    9.958    0.000)   22.877
   5.296   (  14.404    6.719    0.000)   15.894
   5.628   (  -7.397   -9.280    0.000)   11.868
   6.315   (   1.835   -0.163    0.000)    1.843
   6.360   (  -0.119    3.444    0.000)    3.446
   6.376   (  -1.814    3.690    0.000)    4.112
   6.732   (  -1.768   -1.966    0.000)    2.644
   6.941   (  -6.242   -7.336    0.000)    9.632
======================= Grid point 21 (14/135) =======================
q-point: ( 0.40  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.641   (  -6.537   13.288    0.000)   14.809
   2.090   (  -6.091    8.070    0.000)   10.111
   2.169   (   6.640    5.657    0.000)    8.723
   2.604   (  -2.433   -0.918    0.000)    2.601
   4.282   (  14.565    6.076    0.000)   15.781
   5.417   (  -6.172   -9.156    0.000)   11.042
   5.596   (  13.374    4.887    0.000)   14.239
   6.324   (   0.270   -2.508    0.000)    2.523
   6.385   (  -0.910    2.640    0.000)    2.792
   6.388   (  -0.973    3.982    0.000)    4.099
   6.655   (  -4.178   -2.905    0.000)    5.088
   6.790   (  -3.942   -3.669    0.000)    5.385
======================= Grid point 22 (15/135) =======================
q-point: ( 0.47  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.650   (  -7.317   12.674    0.000)   14.634
   2.074   (  -5.146    8.913    0.000)   10.292
   2.261   (  -0.789    1.367    0.000)    1.578
   2.573   (  -0.132    0.228    0.000)    0.264
   4.451   (   1.006   -1.743    0.000)    2.013
   5.303   (   2.527   -4.377    0.000)    5.054
   5.759   (   1.598   -2.769    0.000)    3.197
   6.294   (   1.678   -2.906    0.000)    3.356
   6.391   (  -1.158    2.006    0.000)    2.317
   6.400   (  -2.281    3.951    0.000)    4.562
   6.596   (  -0.127    0.220    0.000)    0.254
   6.742   (   0.212   -0.367    0.000)    0.424
======================= Grid point 33 (16/135) =======================
q-point: ( 0.13  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.293   (   8.773   15.196    0.000)   17.547
   1.415   (  10.169   17.614    0.000)   20.338
   2.021   (   6.272   10.864    0.000)   12.545
   2.555   (   4.373    7.574    0.000)    8.746
   2.668   (  15.018   26.013    0.000)   30.037
   4.820   (   1.127    1.952    0.000)    2.254
   6.012   (  -4.303   -7.453    0.000)    8.606
   6.253   (   0.233    0.404    0.000)    0.466
   6.322   (   1.243    2.152    0.000)    2.485
   6.419   (   1.302    2.256    0.000)    2.604
   6.705   (   1.507    2.610    0.000)    3.014
   7.191   (  -3.066   -5.311    0.000)    6.133
======================= Grid point 34 (17/135) =======================
q-point: ( 0.20  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.585   (   9.725   12.592    0.000)   15.910
   1.699   (   2.698   18.524    0.000)   18.720
   2.200   (   1.465   11.784    0.000)   11.875
   2.651   (   2.034    1.901    0.000)    2.784
   3.175   (  19.694   17.592    0.000)   26.407
   4.923   (   6.473    6.010    0.000)    8.833
   5.847   (  -5.735   -9.199    0.000)   10.840
   6.266   (   0.798    0.417    0.000)    0.900
   6.364   (   0.402    3.497    0.000)    3.520
   6.452   (  -0.960    3.706    0.000)    3.828
   6.737   (   0.789    0.030    0.000)    0.790
   7.077   (  -2.882   -7.928    0.000)    8.435
======================= Grid point 35 (18/135) =======================
q-point: ( 0.27  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.835   (   1.971   13.581    0.000)   13.723
   1.889   (   3.871   14.714    0.000)   15.214
   2.293   (  -3.696   11.447    0.000)   12.029
   2.669   (  -0.568   -0.048    0.000)    0.570
   3.656   (  20.190   12.243    0.000)   23.612
   5.140   (  11.325    7.063    0.000)   13.347
   5.644   (  -6.670  -10.970    0.000)   12.839
   6.286   (   1.282   -0.187    0.000)    1.295
   6.400   (  -0.319    3.647    0.000)    3.661
   6.461   (  -2.447    4.444    0.000)    5.073
   6.732   (  -0.469   -1.545    0.000)    1.615
   6.942   (  -3.544   -8.485    0.000)    9.195
======================= Grid point 36 (19/135) =======================
q-point: ( 0.33  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.919   (  -6.990   16.155    0.000)   17.603
   2.114   (   7.307    9.081    0.000)   11.656
   2.308   (  -6.937   10.848    0.000)   12.877
   2.639   (  -2.466   -0.493    0.000)    2.515
   4.076   (  17.433    8.283    0.000)   19.301
   5.404   (   3.981   -3.148    0.000)    5.076
   5.431   (   2.780   -2.961    0.000)    4.062
   6.292   (   0.447   -2.649    0.000)    2.687
   6.421   (  -0.728    2.618    0.000)    2.717
   6.465   (  -2.115    4.792    0.000)    5.238
   6.688   (  -2.679   -2.352    0.000)    3.565
   6.810   (  -3.399   -5.257    0.000)    6.260
======================= Grid point 37 (20/135) =======================
q-point: ( 0.40  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.939   (  -8.707   15.838    0.000)   18.073
   2.278   (   1.754    6.699    0.000)    6.925
   2.290   (  -6.527   10.544    0.000)   12.401
   2.589   (  -2.111   -0.692    0.000)    2.222
   4.367   (   9.009    2.968    0.000)    9.485
   5.223   (  -1.844  -10.308    0.000)   10.472
   5.658   (   8.695    1.704    0.000)    8.860
   6.241   (  -0.435   -5.833    0.000)    5.849
   6.430   (  -0.906    2.041    0.000)    2.233
   6.482   (  -1.791    4.736    0.000)    5.063
   6.619   (  -2.556   -0.505    0.000)    2.605
   6.743   (  -1.110   -1.179    0.000)    1.620
======================= Grid point 49 (21/135) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.827   (   7.073   12.250    0.000)   14.145
   2.000   (   6.967   12.067    0.000)   13.933
   2.390   (   4.273    7.401    0.000)    8.546
   2.674   (   0.369    0.640    0.000)    0.738
   3.530   (  10.991   19.036    0.000)   21.981
   5.076   (   5.618    9.731    0.000)   11.236
   5.646   (  -6.588  -11.410    0.000)   13.175
   6.273   (   0.156    0.270    0.000)    0.312
   6.424   (   1.418    2.456    0.000)    2.836
   6.510   (   1.160    2.010    0.000)    2.321
   6.727   (  -0.523   -0.905    0.000)    1.045
   6.923   (  -4.291   -7.432    0.000)    8.582
======================= Grid point 50 (22/135) =======================
q-point: ( 0.27  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.049   (   5.842   10.205    0.000)   11.759
   2.172   (   0.517   12.121    0.000)   12.132
   2.481   (  -1.743    7.445    0.000)    7.647
   2.671   (  -1.319    0.018    0.000)    1.319
   3.902   (  14.153   13.411    0.000)   19.498
   5.291   (   8.845    8.329    0.000)   12.149
   5.413   (  -6.777  -12.346    0.000)   14.084
   6.266   (  -0.413   -2.228    0.000)    2.266
   6.452   (  -0.273    1.489    0.000)    1.514
   6.531   (  -0.902    2.524    0.000)    2.680
   6.702   (  -1.264   -1.446    0.000)    1.921
   6.805   (  -2.294   -4.915    0.000)    5.424
======================= Grid point 51 (23/135) =======================
q-point: ( 0.33  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.180   (  -2.676   10.748    0.000)   11.076
   2.302   (   2.953    9.589    0.000)   10.033
   2.494   (  -4.884    7.723    0.000)    9.138
   2.627   (  -3.180   -0.989    0.000)    3.330
   4.239   (  12.026    9.021    0.000)   15.034
   5.170   (  -5.419  -13.139    0.000)   14.213
   5.537   (   9.562    6.595    0.000)   11.615
   6.199   (  -2.106   -6.885    0.000)    7.200
   6.457   (  -0.378    1.327    0.000)    1.380
   6.539   (  -0.834    2.396    0.000)    2.537
   6.658   (  -2.350   -0.606    0.000)    2.427
   6.753   (  -1.053   -0.662    0.000)    1.244
======================= Grid point 52 (24/135) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.204   (  -6.269   10.858    0.000)   12.538
   2.396   (  -3.156    5.466    0.000)    6.312
   2.487   (  -4.844    8.390    0.000)    9.688
   2.578   (   0.169   -0.293    0.000)    0.339
   4.409   (  -0.929    1.610    0.000)    1.859
   5.037   (   4.907   -8.498    0.000)    9.813
   5.675   (  -0.236    0.409    0.000)    0.472
   6.123   (   3.290   -5.699    0.000)    6.581
   6.462   (  -0.809    1.401    0.000)    1.618
   6.547   (  -0.921    1.595    0.000)    1.842
   6.630   (  -0.782    1.355    0.000)    1.564
   6.743   (  -0.590    1.022    0.000)    1.180
======================= Grid point 66 (25/135) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.237   (   5.137    8.898    0.000)   10.274
   2.343   (   3.461    5.994    0.000)    6.921
   2.575   (   1.447    2.507    0.000)    2.895
   2.656   (  -1.151   -1.993    0.000)    2.301
   4.176   (   8.398   14.545    0.000)   16.796
   5.152   (  -7.683  -13.308    0.000)   15.367
   5.486   (   6.088   10.545    0.000)   12.176
   6.183   (  -3.736   -6.471    0.000)    7.472
   6.467   (   0.373    0.645    0.000)    0.745
   6.557   (   0.194    0.336    0.000)    0.388
   6.676   (  -0.756   -1.310    0.000)    1.513
   6.760   (  -0.149   -0.258    0.000)    0.297
======================= Grid point 67 (26/135) =======================
q-point: ( 0.33  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.360   (  -1.481    8.227    0.000)    8.359
   2.447   (   3.596    4.947    0.000)    6.116
   2.590   (  -3.218   -2.838    0.000)    4.291
   2.597   (  -1.521    2.720    0.000)    3.116
   4.434   (   5.442   11.109    0.000)   12.370
   4.914   (  -3.483  -13.153    0.000)   13.606
   5.682   (   4.492    8.361    0.000)    9.491
   6.040   (  -2.456   -9.168    0.000)    9.491
   6.485   (   0.063    1.695    0.000)    1.696
   6.553   (  -0.089   -0.797    0.000)    0.802
   6.655   (  -0.850    0.119    0.000)    0.859
   6.766   (  -0.287    1.235    0.000)    1.268
======================= Grid point 82 (27/135) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.490   (  -0.000   -0.000    0.000)    0.000
   2.502   (  -2.518   -4.362    0.000)    5.037
   2.502   (   2.518    4.362    0.000)    5.037
   2.623   (  -0.000   -0.000    0.000)    0.000
   4.665   (  -6.269  -10.858    0.000)   12.538
   4.665   (   6.269   10.858    0.000)   12.538
   5.857   (  -4.739   -8.207    0.000)    9.477
   5.857   (   4.739    8.207    0.000)    9.477
   6.522   (  -1.001   -1.734    0.000)    2.002
   6.522   (   1.001    1.734    0.000)    2.002
   6.657   (  -0.000   -0.000    0.000)    0.000
   6.780   (  -0.000   -0.000    0.000)    0.000
======================= Grid point 241 (28/135) =======================
q-point: ( 0.00  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.387   (   0.000    0.000   19.699)   19.699
   0.387   (   0.000   -0.000   19.699)   19.699
   0.924   (   0.000   -0.000   47.732)   47.732
   1.512   (   0.000   -0.000   -0.804)    0.804
   1.512   (  -0.000    0.000   -0.804)    0.804
   4.951   (   0.000    0.000  -15.004)   15.004
   6.244   (   0.000   -0.000    0.098)    0.098
   6.244   (   0.000    0.000    0.098)    0.098
   6.364   (   0.000    0.000   -0.876)    0.876
   6.364   (   0.000    0.000   -0.876)    0.876
   6.520   (   0.000    0.000    8.472)    8.472
   7.350   (   0.000   -0.000   -1.566)    1.566
======================= Grid point 242 (29/135) =======================
q-point: ( 0.07  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.546   (  11.909    6.876   14.004)   19.627
   0.715   (  11.241    6.490   22.972)   26.386
   1.137   (  25.465   14.702   28.561)   40.992
   1.555   (   3.671    2.119   -0.645)    4.287
   1.717   (  16.700    9.642   -1.212)   19.322
   4.893   (  -4.693   -2.709  -14.424)   15.408
   6.239   (  -0.330   -0.190    0.110)    0.396
   6.258   (   0.675    0.390    0.753)    1.084
   6.351   (  -0.915   -0.528    3.311)    3.475
   6.369   (   0.287    0.166   -0.924)    0.982
   6.530   (   1.267    0.732    7.287)    7.432
   7.307   (  -1.879   -1.085   -1.766)    2.798
======================= Grid point 243 (30/135) =======================
q-point: ( 0.13  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.842   (  13.596    7.850    9.126)   18.159
   0.949   (  11.523    6.653   18.349)   22.665
   1.667   (   6.105    3.525   -0.398)    7.061
   1.824   (  31.734   18.322   13.907)   39.194
   2.155   (  20.492   11.831   -1.532)   23.711
   4.774   (  -5.060   -2.922  -13.005)   14.258
   6.236   (   0.210    0.121    0.080)    0.256
   6.240   (  -2.562   -1.479    3.051)    4.250
   6.330   (  -0.893   -0.516    4.901)    5.008
   6.370   (  -0.319   -0.184   -0.961)    1.029
   6.593   (   4.152    2.397    3.842)    6.144
   7.278   (  -1.105   -0.638   -3.375)    3.608
======================= Grid point 244 (31/135) =======================
q-point: ( 0.20  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.127   (  11.515    6.648    6.794)   14.932
   1.225   (  12.836    7.411    9.562)   17.638
   1.812   (   6.447    3.722   -0.329)    7.452
   2.423   (  17.257    9.963    2.165)   20.043
   2.631   (  25.098   14.491    5.600)   29.517
   4.713   (   0.559    0.323  -11.010)   11.029
   6.153   (  -4.910   -2.835    5.109)    7.632
   6.250   (   0.945    0.546   -0.037)    1.092
   6.314   (  -0.431   -0.249    4.362)    4.391
   6.355   (  -0.893   -0.516   -0.870)    1.349
   6.688   (   3.829    2.211    1.325)    4.616
   7.244   (  -2.207   -1.274   -3.497)    4.327
======================= Grid point 245 (32/135) =======================
q-point: ( 0.27  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.344   (   7.805    4.506    5.551)   10.585
   1.511   (  12.514    7.225    3.109)   14.780
   1.939   (   4.601    2.656   -0.372)    5.326
   2.618   (   2.738    1.581   -4.096)    5.174
   3.229   (  26.202   15.128    8.949)   31.552
   4.818   (   8.542    4.932   -8.937)   13.310
   6.028   (  -6.128   -3.538    7.186)   10.085
   6.276   (   1.383    0.798   -0.201)    1.609
   6.315   (   0.550    0.317    3.412)    3.470
   6.336   (  -0.678   -0.392   -0.683)    1.039
   6.747   (   1.209    0.698   -0.236)    1.416
   7.169   (  -4.661   -2.691   -3.175)    6.248
======================= Grid point 246 (33/135) =======================
q-point: ( 0.33  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.476   (   4.059    2.343    4.854)    6.747
   1.777   (  11.201    6.467   -1.731)   13.049
   2.008   (   1.514    0.874   -0.340)    1.781
   2.612   (  -2.196   -1.268   -6.093)    6.599
   3.773   (  22.356   12.907    9.014)   27.343
   5.063   (  12.403    7.161   -8.053)   16.430
   5.884   (  -6.569   -3.793   10.072)   12.609
   6.308   (   1.307    0.755   -0.359)    1.551
   6.330   (   0.251    0.145   -0.513)    0.589
   6.338   (   1.422    0.821    2.712)    3.170
   6.742   (  -1.542   -0.890   -1.546)    2.358
   7.032   (  -7.542   -4.354   -2.815)    9.152
======================= Grid point 247 (34/135) =======================
q-point: ( 0.40  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.536   (   1.579    0.911    4.479)    4.836
   2.003   (   8.762    5.059   -5.120)   11.340
   2.012   (  -0.874   -0.504   -0.142)    1.019
   2.550   (  -2.975   -1.718   -7.595)    8.336
   4.221   (  17.383   10.036    9.194)   22.077
   5.313   (   7.808    4.508   -8.396)   12.320
   5.767   (  -2.118   -1.223   13.287)   13.510
   6.332   (   0.919    0.531   -0.433)    1.147
   6.345   (   0.839    0.485   -0.482)    1.082
   6.369   (   1.010    0.583    2.366)    2.638
   6.681   (  -3.907   -2.256   -2.643)    5.228
   6.850   (  -7.812   -4.510   -2.499)    9.360
======================= Grid point 248 (35/135) =======================
q-point: ( 0.47  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.556   (   0.405    0.234    4.292)    4.318
   1.988   (  -0.819   -0.473    0.095)    0.951
   2.152   (   3.962    2.287   -6.831)    8.221
   2.492   (  -1.816   -1.048   -8.634)    8.885
   4.519   (   8.041    4.642    9.397)   13.210
   5.348   (  -2.144   -1.238   -2.371)    3.428
   5.839   (   5.507    3.179    9.008)   11.026
   6.349   (   0.523    0.302   -0.457)    0.757
   6.357   (  -2.072   -1.196    1.609)    2.883
   6.358   (   0.262    0.151   -0.516)    0.598
   6.589   (  -2.768   -1.598   -2.319)    3.949
   6.731   (  -2.529   -1.460   -2.793)    4.040
======================= Grid point 257 (36/135) =======================
q-point: ( 0.07  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.797   (   9.461   16.387    9.643)   21.237
   0.876   (   6.095   10.557   20.550)   23.894
   1.576   (  12.121   20.995   10.690)   26.495
   1.713   (  10.509   18.201    5.047)   21.615
   1.997   (  10.567   18.302   -1.327)   21.175
   4.805   (  -3.226   -5.587  -13.447)   14.914
   6.245   (   0.219    0.379   -0.002)    0.438
   6.249   (  -0.888   -1.539    2.189)    2.819
   6.334   (  -0.526   -0.911    4.155)    4.286
   6.370   (   0.056    0.096    0.103)    0.151
   6.572   (   2.066    3.579    4.488)    6.101
   7.280   (  -0.935   -1.620   -3.015)    3.548
======================= Grid point 258 (37/135) =======================
q-point: ( 0.13  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.100   (   7.607   14.664    8.661)   18.652
   1.130   (   9.654   13.851   10.146)   19.697
   1.806   (   4.479   10.322   -0.198)   11.253
   2.255   (  17.888   23.311    8.779)   30.667
   2.381   (  19.019   11.881    0.691)   22.436
   4.716   (  -1.630   -2.260  -11.810)   12.134
   6.186   (  -3.222   -4.079    4.472)    6.857
   6.253   (   0.102    1.153    0.223)    1.179
   6.320   (  -0.558   -0.103    3.599)    3.644
   6.378   (  -0.262    1.282   -0.248)    1.332
   6.653   (   3.098    3.333    1.944)    4.948
   7.242   (  -0.980   -3.135   -3.372)    4.707
======================= Grid point 259 (38/135) =======================
q-point: ( 0.20  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.356   (   6.115   13.706    4.235)   15.595
   1.411   (   8.047   13.328    6.286)   16.790
   1.964   (   1.712   11.193   -0.351)   11.329
   2.559   (   6.639    4.076   -2.913)    8.318
   2.939   (  25.683   16.978    9.244)   32.145
   4.742   (   4.819    2.876   -9.786)   11.280
   6.072   (  -4.629   -5.564    6.449)    9.694
   6.272   (   0.717    1.227    0.639)    1.558
   6.317   (  -0.332    1.213    2.247)    2.575
   6.379   (  -1.855    2.783   -0.356)    3.364
   6.723   (   2.260    1.081    0.317)    2.525
   7.180   (  -1.798   -5.308   -3.158)    6.433
======================= Grid point 260 (39/135) =======================
q-point: ( 0.27  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.545   (  -0.142   14.589    3.487)   15.000
   1.669   (   9.671    9.356    0.941)   13.488
   2.070   (  -1.596   10.196   -0.340)   10.326
   2.622   (   0.089   -0.694   -5.258)    5.304
   3.507   (  23.994   13.494    9.232)   29.035
   4.924   (  11.011    5.905   -8.277)   14.987
   5.928   (  -5.210   -6.741    8.804)   12.252
   6.299   (   1.329    0.759    1.066)    1.865
   6.333   (  -0.034    2.487    0.861)    2.632
   6.374   (  -2.277    3.619   -0.108)    4.277
   6.744   (   0.154   -0.975   -0.996)    1.402
   7.073   (  -3.974   -7.014   -2.852)    8.551
======================= Grid point 261 (40/135) =======================
q-point: ( 0.33  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.641   (  -4.242   13.471    3.581)   14.570
   1.909   (   9.037    7.196   -3.562)   12.089
   2.107   (  -4.759    8.710   -0.245)    9.928
   2.584   (  -2.680   -1.396   -6.914)    7.545
   3.997   (  20.566   10.215    9.268)   24.763
   5.184   (  11.997    5.049   -8.306)   15.441
   5.779   (  -3.869   -6.713   12.069)   14.342
   6.328   (   1.390   -0.060    0.803)    1.607
   6.359   (  -0.369    3.422    0.038)    3.442
   6.378   (  -1.140    3.231    0.522)    3.466
   6.711   (  -2.174   -2.158   -2.116)    3.723
   6.915   (  -6.179   -7.173   -2.555)    9.806
======================= Grid point 262 (41/135) =======================
q-point: ( 0.40  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.677   (  -6.118   12.486    3.541)   14.349
   2.086   (  -2.734    7.303   -1.679)    7.977
   2.106   (   2.692    5.981   -4.596)    8.009
   2.521   (  -2.927   -1.297   -8.263)    8.861
   4.379   (  14.560    6.287    9.443)   18.457
   5.331   (   1.910   -2.805   -5.369)    6.352
   5.744   (   5.967   -0.894   11.384)   12.884
   6.331   (   0.090   -2.594    0.374)    2.623
   6.382   (  -1.258    4.009   -0.386)    4.219
   6.394   (  -0.550    2.213    0.959)    2.474
   6.630   (  -4.153   -2.612   -2.579)    5.543
   6.764   (  -4.107   -3.709   -2.525)    6.083
======================= Grid point 263 (42/135) =======================
q-point: ( 0.47  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.685   (  -6.902   11.955    3.507)   14.243
   2.075   (  -5.121    8.869    0.088)   10.242
   2.187   (  -0.753    1.303   -7.075)    7.233
   2.482   (  -0.079    0.136   -8.886)    8.887
   4.549   (   0.926   -1.604    9.560)    9.738
   5.296   (   2.097   -3.633   -0.484)    4.223
   5.837   (   1.976   -3.423    7.475)    8.456
   6.298   (   1.837   -3.183    0.072)    3.676
   6.395   (  -1.222    2.117    0.493)    2.494
   6.395   (  -2.256    3.908   -0.441)    4.534
   6.576   (  -0.193    0.334   -1.891)    1.930
   6.714   (   0.160   -0.277   -2.868)    2.886
======================= Grid point 274 (43/135) =======================
q-point: ( 0.13  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.338   (   7.068   12.243    5.460)   15.154
   1.466   (   9.738   16.867    5.072)   20.125
   2.021   (   6.278   10.875   -0.170)   12.558
   2.520   (   3.935    6.816   -3.505)    8.616
   2.776   (  14.818   25.665   10.757)   31.528
   4.715   (   1.567    2.715  -10.184)   10.656
   6.083   (  -3.721   -6.445    6.332)    9.771
   6.271   (   0.417    0.722    1.340)    1.578
   6.333   (   0.750    1.299    1.358)    2.023
   6.414   (   1.214    2.102   -0.442)    2.467
   6.707   (   1.165    2.018    0.488)    2.381
   7.161   (  -2.941   -5.094   -2.996)    6.601
======================= Grid point 275 (44/135) =======================
q-point: ( 0.20  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.576   (   7.882   10.636    0.055)   13.238
   1.742   (   3.136   17.430    4.147)   18.189
   2.200   (   1.337   11.974   -0.111)   12.048
   2.603   (   1.672    1.031   -4.898)    5.277
   3.277   (  19.261   17.917   10.069)   28.167
   4.830   (   6.734    6.233   -8.697)   12.643
   5.941   (  -4.873   -7.892    8.442)   12.541
   6.285   (   0.720    0.327    1.753)    1.923
   6.364   (   0.297    3.100    0.128)    3.117
   6.446   (  -0.882    3.607   -0.576)    3.758
   6.728   (   0.455   -0.417   -0.683)    0.921
   7.050   (  -2.854   -7.676   -2.652)    8.608
======================= Grid point 276 (45/135) =======================
q-point: ( 0.27  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.795   (   5.589    9.950   -3.593)   11.965
   1.917   (  -0.585   15.828    3.037)   16.127
   2.292   (  -3.801   11.591   -0.154)   12.199
   2.607   (  -1.135   -0.622   -6.275)    6.407
   3.757   (  19.935   12.644    9.805)   25.562
   5.042   (  10.484    6.198   -8.128)   14.643
   5.775   (  -4.793   -8.886   11.127)   15.025
   6.301   (   1.115   -0.433    1.433)    1.867
   6.397   (  -0.247    3.417   -0.204)    3.432
   6.456   (  -2.311    4.311   -0.470)    4.914
   6.715   (  -0.772   -1.756   -1.651)    2.530
   6.918   (  -3.551   -8.233   -2.403)    9.283
======================= Grid point 277 (46/135) =======================
q-point: ( 0.33  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.921   (  -3.888   13.153   -0.921)   13.746
   2.070   (   3.961   10.160   -2.853)   11.272
   2.307   (  -6.919   10.899   -0.158)   12.911
   2.562   (  -3.104   -0.937   -7.673)    8.330
   4.176   (  17.281    8.636    9.806)   21.665
   5.242   (   6.447    0.271   -7.429)    9.840
   5.651   (   1.059   -5.469   13.002)   14.145
   6.301   (   0.175   -2.957    0.799)    3.068
   6.419   (  -0.417    2.327   -0.101)    2.367
   6.460   (  -2.074    4.542   -0.500)    5.018
   6.667   (  -2.802   -2.068   -2.091)    4.062
   6.787   (  -3.513   -5.061   -2.245)    6.557
======================= Grid point 278 (47/135) =======================
q-point: ( 0.40  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.952   (  -7.447   14.023    0.732)   15.894
   2.210   (   1.780    6.620   -5.961)    9.085
   2.290   (  -7.350   11.096   -0.149)   13.310
   2.499   (  -2.430   -0.994   -8.801)    9.184
   4.469   (   8.894    3.227    9.862)   13.667
   5.225   (  -0.216   -7.260    0.147)    7.264
   5.730   (   7.303   -0.784    6.985)   10.136
   6.243   (  -0.641   -6.079    0.056)    6.113
   6.430   (  -0.628    1.615    0.060)    1.734
   6.474   (  -1.712    4.360   -0.767)    4.747
   6.603   (  -2.548    0.059   -1.583)    3.001
   6.719   (  -1.357   -0.909   -2.414)    2.915
======================= Grid point 290 (48/135) =======================
q-point: ( 0.20  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.781   (   6.011   10.412   -3.899)   12.639
   2.027   (   6.645   11.509    2.668)   13.554
   2.393   (   4.323    7.487    0.208)    8.648
   2.612   (   0.006    0.011   -6.211)    6.211
   3.636   (  11.056   19.149   10.347)   24.412
   4.980   (   5.237    9.070   -8.164)   13.279
   5.771   (  -5.382   -9.322   10.895)   15.315
   6.289   (   0.005    0.009    1.464)    1.464
   6.419   (   1.306    2.262   -0.433)    2.648
   6.503   (   1.149    1.989   -0.710)    2.404
   6.712   (  -0.684   -1.184   -1.476)    2.012
   6.901   (  -4.172   -7.226   -2.231)    8.637
======================= Grid point 291 (49/135) =======================
q-point: ( 0.27  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.973   (   5.503    8.860   -6.959)   12.538
   2.190   (   0.323   11.295    1.635)   11.417
   2.484   (  -1.976    7.612    0.114)    7.865
   2.598   (  -1.668   -0.473   -7.058)    7.267
   4.009   (  13.995   13.684   10.391)   22.160
   5.154   (   5.829    4.608   -7.864)   10.819
   5.605   (  -3.017   -7.491   13.200)   15.475
   6.276   (  -0.609   -2.558    0.923)    2.787
   6.444   (  -0.255    1.136   -0.803)    1.414
   6.523   (  -0.913    2.301   -0.844)    2.616
   6.684   (  -1.248   -1.303   -1.726)    2.496
   6.787   (  -2.285   -4.527   -1.825)    5.389
======================= Grid point 292 (50/135) =======================
q-point: ( 0.33  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.118   (   1.112    8.255   -7.298)   11.075
   2.287   (  -0.414   10.326   -0.315)   10.339
   2.492   (  -5.096    7.204   -1.079)    8.890
   2.545   (  -3.744   -0.613   -6.992)    7.955
   4.346   (  11.832    9.214   10.296)   18.191
   5.170   (  -1.985   -7.535   -0.022)    7.792
   5.612   (   6.553    1.837    7.300)    9.980
   6.202   (  -2.314   -7.147    0.253)    7.516
   6.447   (  -0.139    0.884   -0.990)    1.335
   6.526   (  -0.907    1.859   -1.327)    2.457
   6.646   (  -2.159   -0.075   -1.132)    2.438
   6.738   (  -1.264   -0.165   -1.561)    2.015
======================= Grid point 293 (51/135) =======================
q-point: ( 0.40  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.164   (  -4.416    7.650   -5.629)   10.474
   2.360   (  -4.489    7.775   -1.613)    9.122
   2.463   (  -1.537    2.662   -7.702)    8.293
   2.503   (  -2.955    5.118   -1.979)    6.233
   4.514   (  -1.007    1.744   10.201)   10.398
   5.082   (   3.931   -6.809    4.186)    8.907
   5.713   (   0.483   -0.836    3.915)    4.032
   6.123   (   3.287   -5.693   -0.088)    6.575
   6.453   (  -0.524    0.908   -0.997)    1.447
   6.529   (  -0.632    1.094   -1.838)    2.230
   6.624   (  -1.016    1.759   -0.544)    2.103
   6.728   (  -0.888    1.538   -1.489)    2.318
======================= Grid point 307 (52/135) =======================
q-point: ( 0.27  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.139   (   4.658    8.068   -9.230)   13.114
   2.347   (   3.046    5.275   -0.385)    6.103
   2.576   (  -1.239   -2.146   -7.154)    7.570
   2.580   (   1.463    2.534    0.461)    2.962
   4.285   (   8.373   14.503   10.561)   19.798
   5.150   (  -3.588   -6.215   -0.252)    7.181
   5.564   (   2.534    4.389    7.586)    9.123
   6.186   (  -3.894   -6.744    0.323)    7.795
   6.453   (   0.242    0.419   -1.433)    1.512
   6.542   (  -0.059   -0.101   -1.597)    1.601
   6.665   (  -0.457   -0.792   -1.009)    1.361
   6.749   (  -0.002   -0.004   -1.106)    1.106
======================= Grid point 308 (53/135) =======================
q-point: ( 0.33  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.266   (   0.657    7.194   -9.815)   12.186
   2.423   (   2.255    3.025   -4.482)    5.859
   2.510   (  -4.134   -1.232   -4.730)    6.402
   2.602   (  -1.487    3.029    0.525)    3.415
   4.540   (   5.325   10.866   10.238)   15.850
   4.983   (  -3.010  -11.040    6.442)   13.132
   5.703   (   4.192    6.933    2.264)    8.412
   6.041   (  -2.514   -9.136    0.065)    9.476
   6.468   (   0.155    1.360   -1.767)    2.235
   6.531   (  -0.222   -1.084   -2.295)    2.547
   6.653   (  -0.707    0.532   -0.114)    0.892
   6.758   (  -0.351    1.433   -0.860)    1.708
======================= Grid point 323 (54/135) =======================
q-point: ( 0.33  0.33  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.389   (  -2.532   -4.386  -10.976)   12.088
   2.389   (   2.532    4.386  -10.976)   12.088
   2.509   (  -0.000   -0.000    1.855)    1.855
   2.629   (  -0.000   -0.000    0.611)    0.611
   4.760   (  -5.829  -10.096    9.202)   14.852
   4.760   (   5.829   10.096    9.202)   14.852
   5.862   (  -4.523   -7.833    0.520)    9.060
   5.862   (   4.523    7.833    0.520)    9.060
   6.500   (  -0.913   -1.581   -2.290)    2.929
   6.500   (   0.913    1.581   -2.290)    2.929
   6.659   (  -0.000   -0.000    0.218)    0.218
   6.774   (  -0.000   -0.000   -0.706)    0.706
======================= Grid point 482 (55/135) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.744   (   0.000   -0.000   17.304)   17.304
   0.744   (   0.000    0.000   17.304)   17.304
   1.483   (  -0.000    0.000   -2.530)    2.530
   1.483   (  -0.000    0.000   -2.530)    2.530
   1.816   (  -0.000   -0.000   45.205)   45.205
   4.567   (   0.000    0.000  -24.349)   24.349
   6.248   (   0.000   -0.000    0.371)    0.371
   6.248   (   0.000    0.000    0.371)    0.371
   6.341   (   0.000    0.000   -1.480)    1.480
   6.341   (   0.000    0.000   -1.480)    1.480
   6.713   (   0.000    0.000   10.886)   10.886
   7.303   (   0.000   -0.000   -3.391)    3.391
======================= Grid point 483 (56/135) =======================
q-point: ( 0.07  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.838   (   7.809    4.508   15.410)   17.854
   1.046   (  18.998   10.969   14.029)   26.040
   1.529   (   3.984    2.300   -2.274)    5.132
   1.677   (  15.282    8.823   -1.784)   17.737
   1.866   (   7.519    4.341   41.080)   41.988
   4.523   (  -3.551   -2.050  -23.574)   23.927
   6.244   (  -0.308   -0.178    0.417)    0.548
   6.270   (   1.577    0.910    0.613)    1.921
   6.344   (   0.196    0.113   -1.576)    1.592
   6.373   (   2.059    1.189   -0.453)    2.421
   6.706   (  -0.434   -0.251   10.278)   10.291
   7.260   (  -3.012   -1.739   -3.252)    4.761
======================= Grid point 484 (57/135) =======================
q-point: ( 0.13  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.056   (  10.854    6.267   12.273)   17.541
   1.349   (   8.094    4.673   20.040)   22.112
   1.648   (   6.345    3.663   -1.812)    7.547
   2.063   (  17.985   10.384    0.458)   20.773
   2.252   (  24.701   14.261   22.150)   36.113
   4.439   (  -2.986   -1.724  -21.287)   21.564
   6.240   (   0.120    0.069    0.387)    0.411
   6.289   (  -0.501   -0.289    1.962)    2.046
   6.344   (  -0.289   -0.167   -1.640)    1.674
   6.405   (   0.609    0.351    2.119)    2.232
   6.709   (   0.924    0.533    7.940)    8.012
   7.201   (  -2.028   -1.171   -4.370)    4.958
======================= Grid point 485 (58/135) =======================
q-point: ( 0.20  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.293   (   9.913    5.723    9.910)   15.140
   1.492   (   5.878    3.394   17.073)   18.372
   1.796   (   6.482    3.742   -1.546)    7.643
   2.393   (  10.113    5.839   -4.370)   12.468
   2.855   (  27.794   16.047   16.477)   36.076
   4.436   (   3.494    2.017  -17.340)   17.803
   6.246   (  -2.963   -1.711    4.111)    5.348
   6.250   (   0.750    0.433    0.159)    0.881
   6.332   (  -0.694   -0.401   -1.452)    1.658
   6.402   (  -0.673   -0.388    3.620)    3.702
   6.741   (   1.516    0.875    4.500)    4.829
   7.160   (  -2.080   -1.201   -4.899)    5.456
======================= Grid point 486 (59/135) =======================
q-point: ( 0.27  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.483   (   6.851    3.956    8.376)   11.522
   1.628   (   6.243    3.604    9.251)   11.728
   1.923   (   4.601    2.656   -1.472)    5.513
   2.502   (   0.342    0.197   -7.766)    7.776
   3.458   (  25.751   14.867   14.294)   32.992
   4.605   (  11.131    6.427  -12.667)   18.046
   6.169   (  -3.793   -2.190    6.526)    7.860
   6.272   (   1.167    0.674   -0.172)    1.359
   6.318   (  -0.431   -0.249   -1.089)    1.197
   6.390   (  -0.150   -0.087    3.513)    3.517
   6.755   (  -0.597   -0.345    1.623)    1.763
   7.092   (  -4.266   -2.463   -4.543)    6.701
======================= Grid point 487 (60/135) =======================
q-point: ( 0.33  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.598   (   3.512    2.027    7.435)    8.469
   1.768   (   6.250    3.609    1.585)    7.389
   1.993   (   1.667    0.963   -1.335)    2.343
   2.453   (  -3.833   -2.213  -10.086)   11.015
   3.996   (  22.264   12.854   13.466)   29.022
   4.894   (  13.760    7.944   -9.128)   18.324
   6.082   (  -3.871   -2.235    9.418)   10.425
   6.299   (   1.152    0.665   -0.458)    1.407
   6.317   (   0.390    0.225   -0.786)    0.906
   6.400   (   1.017    0.587    3.061)    3.278
   6.711   (  -3.319   -1.916   -1.143)    3.999
   6.964   (  -7.104   -4.102   -3.959)    9.109
======================= Grid point 488 (61/135) =======================
q-point: ( 0.40  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.648   (   1.267    0.732    6.866)    7.020
   1.902   (   5.505    3.178   -4.419)    7.742
   2.003   (  -0.561   -0.324   -0.968)    1.165
   2.359   (  -4.191   -2.420  -11.636)   12.602
   4.445   (  17.509   10.109   13.276)   24.187
   5.175   (  10.272    5.931   -6.206)   13.387
   6.010   (  -1.838   -1.061   11.405)   11.601
   6.322   (   0.887    0.512   -0.569)    1.172
   6.333   (   0.802    0.463   -0.736)    1.183
   6.428   (   1.039    0.600    3.327)    3.537
   6.612   (  -5.314   -3.068   -4.229)    7.453
   6.791   (  -7.577   -4.375   -3.332)    9.362
======================= Grid point 489 (62/135) =======================
q-point: ( 0.47  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.664   (   0.283    0.163    6.545)    6.553
   1.985   (  -0.634   -0.366   -0.562)    0.923
   2.001   (   2.847    1.644   -7.841)    8.502
   2.280   (  -2.338   -1.350  -12.463)   12.752
   4.747   (   8.183    4.725   13.338)   16.346
   5.308   (   1.701    0.982   -1.815)    2.674
   6.031   (   3.084    1.781    9.955)   10.573
   6.338   (   0.491    0.284   -0.629)    0.847
   6.345   (   0.252    0.146   -0.759)    0.813
   6.387   (  -4.745   -2.740    0.881)    5.549
   6.536   (  -0.086   -0.050   -2.638)    2.640
   6.663   (  -3.431   -1.981   -3.966)    5.606
======================= Grid point 498 (63/135) =======================
q-point: ( 0.07  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.020   (   7.414   12.842   12.728)   19.542
   1.279   (   5.724    9.914   18.221)   21.519
   1.628   (   4.306    7.458   -1.326)    8.713
   1.956   (  10.153   17.585   -2.937)   20.517
   2.092   (  11.761   20.370   30.056)   38.165
   4.459   (  -2.184   -3.784  -22.062)   22.490
   6.248   (   0.206    0.357    0.365)    0.551
   6.278   (   0.009    0.015    0.930)    0.930
   6.362   (   0.572    0.990   -0.960)    1.493
   6.390   (   0.232    0.402    1.260)    1.343
   6.705   (   0.241    0.418    8.737)    8.751
   7.211   (  -1.574   -2.725   -3.948)    5.049
======================= Grid point 499 (64/135) =======================
q-point: ( 0.13  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.280   (   5.981   15.434   10.212)   19.449
   1.431   (   5.714    4.235   18.545)   19.862
   1.791   (   4.130   10.389   -1.508)   11.281
   2.279   (  10.682   10.092   -3.785)   15.175
   2.601   (  21.994   21.227   19.322)   36.161
   4.415   (   0.388   -0.048  -18.990)   18.994
   6.255   (  -0.473    0.531    0.927)    1.168
   6.263   (  -0.761   -1.248    2.177)    2.622
   6.361   (  -1.721    1.328   -0.353)    2.203
   6.397   (   0.767   -0.880    2.322)    2.599
   6.726   (   1.158    1.020    5.710)    5.915
   7.162   (  -1.228   -2.958   -4.655)    5.650
======================= Grid point 500 (65/135) =======================
q-point: ( 0.20  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.490   (   3.274   11.474    8.960)   14.922
   1.582   (   5.327    7.744   11.438)   14.805
   1.949   (   1.686   11.345   -1.413)   11.556
   2.467   (   4.549    1.680   -6.556)    8.155
   3.184   (  24.327   17.447   15.495)   33.709
   4.502   (   7.347    4.955  -14.618)   17.095
   6.196   (  -2.973   -3.426    5.627)    7.228
   6.277   (   0.274    2.120    0.139)    2.142
   6.346   (  -1.418    1.453   -0.398)    2.069
   6.396   (  -0.421    0.308    2.267)    2.327
   6.748   (   0.648   -0.440    2.707)    2.818
   7.103   (  -1.733   -4.729   -4.530)    6.774
======================= Grid point 501 (66/135) =======================
q-point: ( 0.27  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.629   (   2.322    8.119    5.421)   10.035
   1.736   (   2.703    9.484    5.983)   11.535
   2.056   (  -1.618   10.364   -1.305)   10.570
   2.482   (  -1.632   -1.857   -9.073)    9.404
   3.737   (  23.326   13.785   14.008)   30.502
   4.739   (  12.822    7.158  -10.574)   18.095
   6.104   (  -3.229   -4.502    8.369)   10.037
   6.305   (   0.392    2.228   -0.170)    2.268
   6.344   (  -1.098    2.497   -0.378)    2.754
   6.395   (  -0.224    1.039    2.168)    2.414
   6.729   (  -1.377   -2.233    0.022)    2.623
   7.003   (  -3.809   -6.495   -4.086)    8.567
======================= Grid point 502 (67/135) =======================
q-point: ( 0.33  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.721   (  -1.240    8.736    3.983)    9.681
   1.861   (   3.037    7.021   -0.042)    7.650
   2.095   (  -4.589    8.810   -1.175)   10.003
   2.408   (  -4.153   -2.135  -10.810)   11.776
   4.223   (  20.352   10.591   13.497)   26.619
   5.031   (  13.177    6.007   -7.580)   16.346
   6.009   (  -2.103   -5.121   10.901)   12.226
   6.328   (   0.357    1.004   -0.545)    1.197
   6.356   (  -1.292    3.848   -0.396)    4.078
   6.412   (   1.043    0.881    2.445)    2.800
   6.660   (  -3.722   -2.945   -2.685)    5.453
   6.853   (  -5.930   -6.825   -3.581)    9.725
======================= Grid point 503 (68/135) =======================
q-point: ( 0.40  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.760   (  -4.228    9.352    4.432)   11.179
   1.974   (   3.105    4.913   -5.498)    8.000
   2.085   (  -5.843    8.048   -0.967)    9.992
   2.318   (  -3.932   -1.740  -11.971)   12.720
   4.607   (  14.518    6.629   13.409)   20.845
   5.245   (   6.744    1.164   -3.754)    7.805
   5.965   (   2.833   -3.322   10.819)   11.667
   6.327   (  -0.320   -2.272   -0.847)    2.446
   6.373   (  -1.346    4.165   -0.414)    4.396
   6.432   (  -0.417    0.219    2.819)    2.858
   6.561   (  -4.032   -1.510   -4.485)    6.218
   6.703   (  -4.591   -4.052   -3.514)    7.060
======================= Grid point 504 (69/135) =======================
q-point: ( 0.47  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.769   (  -5.345    9.257    4.586)   11.632
   2.034   (  -1.275    2.209   -7.457)    7.882
   2.071   (  -4.890    8.469   -0.870)    9.818
   2.268   (  -0.126    0.218  -12.373)   12.376
   4.779   (   0.800   -1.385   13.465)   13.560
   5.289   (   1.452   -2.514   -0.316)    2.920
   6.006   (   2.499   -4.329    9.040)   10.330
   6.285   (   2.032   -3.519   -1.494)    4.330
   6.386   (  -2.199    3.808   -0.415)    4.417
   6.410   (  -1.521    2.635    0.881)    3.168
   6.536   (  -0.020    0.034   -1.883)    1.883
   6.640   (   0.037   -0.065   -4.512)    4.512
======================= Grid point 515 (70/135) =======================
q-point: ( 0.13  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.517   (   3.028    5.245   12.795)   14.156
   1.594   (   8.737   15.132    7.845)   19.153
   2.009   (   6.417   11.115   -1.202)   12.890
   2.424   (   2.787    4.828   -6.488)    8.553
   3.047   (  14.043   24.324   16.802)   32.729
   4.463   (   3.147    5.451  -15.471)   16.702
   6.205   (  -2.221   -3.847    5.398)    6.991
   6.288   (   0.954    1.653    0.513)    1.976
   6.370   (  -0.273   -0.472    1.927)    2.002
   6.402   (   0.949    1.643   -0.886)    2.094
   6.738   (   0.037    0.065    3.142)    3.143
   7.088   (  -2.679   -4.640   -4.323)    6.884
======================= Grid point 516 (71/135) =======================
q-point: ( 0.20  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.624   (   3.823    5.075    5.408)    8.344
   1.844   (   3.115   15.162    6.179)   16.667
   2.191   (   1.079   12.467   -0.934)   12.548
   2.472   (   0.812   -0.660   -8.529)    8.593
   3.526   (  18.242   18.057   15.105)   29.782
   4.628   (   8.380    8.079  -11.723)   16.520
   6.112   (  -3.385   -5.241    8.365)   10.435
   6.315   (   0.488    1.102    1.013)    1.574
   6.371   (   0.005    1.602    0.494)    1.677
   6.432   (  -0.361    2.954   -0.902)    3.110
   6.724   (  -0.732   -1.904    0.656)    2.142
   6.985   (  -2.779   -7.100   -3.865)    8.548
======================= Grid point 517 (72/135) =======================
q-point: ( 0.27  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.740   (   4.259    4.743   -1.081)    6.466
   1.993   (  -1.934   14.079    4.336)   14.858
   2.280   (  -4.146   11.568   -1.368)   12.364
   2.450   (  -2.249   -0.878   -9.240)    9.550
   3.995   (  19.279   13.184   14.207)   27.338
   4.879   (  11.699    7.357   -8.574)   16.264
   5.997   (  -3.333   -6.521   10.947)   13.171
   6.325   (   0.410   -0.462    0.621)    0.876
   6.397   (   0.212    2.638    0.222)    2.656
   6.445   (  -1.699    3.527   -0.627)    3.965
   6.679   (  -1.881   -2.520   -1.624)    3.539
   6.859   (  -3.538   -7.697   -3.529)    9.177
======================= Grid point 518 (73/135) =======================
q-point: ( 0.33  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.849   (   2.833    4.848   -5.030)    7.539
   2.071   (  -3.491   12.561    1.870)   13.171
   2.286   (  -6.563    8.860   -3.256)   11.496
   2.388   (  -5.113    1.404   -8.129)    9.705
   4.412   (  16.915    9.168   13.829)   23.694
   5.124   (   9.650    3.607   -5.261)   11.568
   5.895   (  -0.360   -6.298   11.911)   13.479
   6.310   (  -0.706   -3.211   -0.198)    3.293
   6.422   (   0.452    1.249    0.365)    1.378
   6.449   (  -1.645    3.259   -0.619)    3.703
   6.616   (  -3.375   -1.087   -2.852)    4.551
   6.732   (  -3.762   -4.661   -3.315)    6.846
======================= Grid point 519 (74/135) =======================
q-point: ( 0.40  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.916   (  -1.034    5.333   -5.228)    7.539
   2.134   (  -2.973   11.508   -0.841)   11.916
   2.242   (  -5.376    3.542   -8.287)   10.494
   2.328   (  -5.857    6.503   -4.165)    9.692
   4.704   (   8.708    3.614   13.691)   16.623
   5.228   (   2.532   -2.557    0.016)    3.599
   5.892   (   5.095   -3.994    9.055)   11.132
   6.234   (  -1.375   -6.332   -1.114)    6.574
   6.434   (   0.293    0.033    0.404)    0.501
   6.454   (  -1.355    2.900   -1.245)    3.435
   6.565   (  -2.538    1.808   -2.025)    3.716
   6.659   (  -2.137   -0.100   -3.602)    4.189
======================= Grid point 531 (75/135) =======================
q-point: ( 0.20  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.725   (   3.149    5.454   -0.906)    6.363
   2.087   (   5.322    9.217    2.917)   11.036
   2.394   (   4.497    7.789   -0.261)    8.998
   2.457   (  -0.308   -0.534   -9.107)    9.128
   3.886   (  11.028   19.101   14.799)   26.561
   4.811   (   6.138   10.631   -9.041)   15.246
   5.994   (  -3.988   -6.907   11.182)   13.735
   6.317   (  -0.308   -0.534    1.074)    1.238
   6.410   (   1.031    1.785   -0.403)    2.101
   6.483   (   1.078    1.868   -1.367)    2.553
   6.682   (  -1.259   -2.181   -1.212)    2.795
   6.845   (  -3.911   -6.773   -3.347)    8.507
======================= Grid point 532 (76/135) =======================
q-point: ( 0.27  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.834   (   3.712    5.085   -6.239)    8.864
   2.201   (  -0.050    6.249   -2.549)    6.749
   2.435   (  -3.469    3.440   -5.714)    7.518
   2.496   (  -1.012    6.401   -0.625)    6.510
   4.257   (  13.621   13.999   14.445)   24.293
   5.023   (   7.809    7.392   -6.100)   12.363
   5.859   (  -3.162   -7.305   12.675)   14.967
   6.292   (  -1.098   -3.152    0.405)    3.362
   6.425   (  -0.231   -0.004   -1.011)    1.037
   6.498   (  -0.951    1.283   -1.853)    2.446
   6.646   (  -1.280   -0.708   -1.845)    2.355
   6.741   (  -2.221   -3.559   -2.782)    5.033
======================= Grid point 533 (77/135) =======================
q-point: ( 0.33  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.939   (   2.654    4.704   -9.683)   11.087
   2.228   (  -0.839    1.854   -9.620)    9.833
   2.399   (  -6.260    6.685   -1.156)    9.231
   2.509   (  -4.750    7.445    0.050)    8.832
   4.590   (  11.445    9.409   14.072)   20.433
   5.166   (   3.203    0.391   -0.589)    3.281
   5.787   (   2.314   -4.059   10.092)   11.121
   6.203   (  -2.918   -7.519   -0.396)    8.075
   6.421   (   0.456   -0.458   -1.625)    1.749
   6.487   (  -0.898    0.123   -2.794)    2.937
   6.624   (  -1.862    1.457   -0.973)    2.557
   6.700   (  -1.690    1.068   -2.236)    2.999
======================= Grid point 534 (78/135) =======================
q-point: ( 0.40  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.990   (  -1.718    2.976  -10.624)   11.166
   2.226   (   0.604   -1.047  -14.408)   14.458
   2.375   (  -5.901   10.221    2.508)   12.065
   2.502   (  -4.795    8.305    0.218)    9.592
   4.755   (  -1.064    1.842   13.827)   13.990
   5.167   (   1.824   -3.160    4.014)    5.424
   5.817   (   2.002   -3.467    6.644)    7.757
   6.115   (   3.117   -5.400   -0.807)    6.287
   6.426   (   0.294   -0.509   -1.715)    1.813
   6.478   (   0.087   -0.150   -3.379)    3.384
   6.614   (  -1.555    2.694   -0.378)    3.133
   6.692   (  -1.573    2.725   -2.060)    3.761
======================= Grid point 548 (79/135) =======================
q-point: ( 0.27  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.938   (   3.180    5.508  -10.331)   12.132
   2.242   (  -0.461   -0.798  -11.468)   11.505
   2.493   (   1.285    2.226    0.282)    2.586
   2.589   (   1.557    2.697    0.385)    3.138
   4.534   (   8.241   14.274   14.312)   21.829
   5.139   (   1.908    3.305   -1.019)    3.950
   5.747   (  -1.721   -2.981   10.676)   11.217
   6.190   (  -4.204   -7.281   -0.113)    8.408
   6.415   (  -0.156   -0.269   -2.332)    2.353
   6.493   (  -0.694   -1.202   -3.453)    3.722
   6.647   (   0.250    0.433   -0.613)    0.791
   6.721   (   0.347    0.600   -1.732)    1.866
======================= Grid point 549 (80/135) =======================
q-point: ( 0.33  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.035   (   1.363    5.130  -12.804)   13.861
   2.205   (  -0.829   -2.936  -14.620)   14.935
   2.516   (  -1.992    4.159    2.462)    5.228
   2.613   (  -1.516    3.236    0.557)    3.617
   4.780   (   5.066   10.167   13.661)   17.766
   5.120   (  -1.473   -5.324    6.833)    8.786
   5.769   (   3.032    2.653    4.646)    6.150
   6.040   (  -2.658   -8.876   -0.270)    9.269
   6.421   (   0.347    0.474   -2.988)    3.045
   6.468   (  -0.337   -1.510   -4.030)    4.316
   6.653   (  -0.550    1.267    0.209)    1.397
   6.735   (  -0.454    1.829   -1.387)    2.340
======================= Grid point 564 (81/135) =======================
q-point: ( 0.33  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 192
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.131   (  -2.290   -3.967  -14.623)   15.323
   2.131   (   2.290    3.967  -14.623)   15.323
   2.555   (  -0.000   -0.000    2.745)    2.745
   2.643   (  -0.000   -0.000    0.703)    0.703
   4.973   (  -4.660   -8.071   11.805)   15.040
   4.973   (   4.660    8.071   11.805)   15.040
   5.876   (  -3.794   -6.572    0.900)    7.642
   5.876   (   3.794    6.572    0.900)    7.642
   6.439   (  -0.674   -1.167   -3.838)    4.067
   6.439   (   0.674    1.167   -3.838)    4.067
   6.666   (  -0.000   -0.000    0.536)    0.536
   6.755   (  -0.000   -0.000   -1.228)    1.228
======================= Grid point 723 (82/135) =======================
q-point: ( 0.00  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.042   (   0.000    0.000   13.706)   13.706
   1.042   (  -0.000    0.000   13.706)   13.706
   1.407   (   0.000   -0.000   -5.541)    5.541
   1.407   (   0.000    0.000   -5.541)    5.541
   2.645   (  -0.000   -0.000   41.325)   41.325
   4.036   (   0.000    0.000  -30.906)   30.906
   6.260   (   0.000   -0.000    0.848)    0.848
   6.260   (  -0.000    0.000    0.848)    0.848
   6.310   (  -0.000    0.000   -1.638)    1.638
   6.310   (  -0.000   -0.000   -1.638)    1.638
   6.916   (   0.000   -0.000   10.003)   10.003
   7.218   (   0.000    0.000   -5.567)    5.567
======================= Grid point 724 (83/135) =======================
q-point: ( 0.07  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.112   (   5.942    3.431   12.889)   14.602
   1.289   (  18.137   10.471   11.425)   23.856
   1.459   (   4.420    2.552   -5.229)    7.307
   1.607   (  15.853    9.153   -5.104)   19.003
   2.655   (   1.895    1.094   40.005)   40.065
   4.008   (  -2.128   -1.228  -30.054)   30.154
   6.256   (  -0.236   -0.136    0.944)    0.982
   6.286   (   2.009    1.160    1.132)    2.582
   6.312   (   0.054    0.031   -1.769)    1.770
   6.349   (   2.894    1.671   -1.753)    3.774
   6.899   (  -1.350   -0.780    9.504)    9.631
   7.180   (  -2.904   -1.677   -5.123)    6.123
======================= Grid point 725 (84/135) =======================
q-point: ( 0.13  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.284   (   8.866    5.119   11.139)   15.129
   1.589   (   6.814    3.934   -4.546)    9.087
   1.664   (  12.493    7.213   13.140)   19.513
   1.999   (  16.775    9.685   -5.232)   20.065
   2.804   (  12.960    7.483   32.025)   35.350
   3.975   (   0.262    0.151  -26.984)   26.986
   6.253   (  -0.004   -0.002    0.948)    0.948
   6.310   (  -0.227   -0.131   -1.840)    1.859
   6.324   (   0.791    0.457    1.959)    2.162
   6.408   (   1.861    1.074   -1.432)    2.582
   6.869   (  -1.032   -0.596    8.214)    8.300
   7.114   (  -2.570   -1.484   -4.820)    5.661
======================= Grid point 726 (85/135) =======================
q-point: ( 0.20  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.483   (   8.434    4.870    9.434)   13.559
   1.745   (   6.742    3.892   -3.925)    8.718
   1.822   (   2.875    1.660   16.255)   16.590
   2.277   (   7.177    4.144   -7.712)   11.320
   3.222   (  22.701   13.107   21.257)   33.749
   4.066   (   8.520    4.919  -20.949)   23.144
   6.258   (   0.430    0.248    0.640)    0.810
   6.303   (  -0.371   -0.214   -1.562)    1.620
   6.313   (  -1.490   -0.860    3.223)    3.653
   6.424   (  -0.253   -0.146   -1.167)    1.203
   6.853   (  -0.689   -0.398    6.558)    6.606
   7.067   (  -1.986   -1.146   -4.748)    5.273
======================= Grid point 727 (86/135) =======================
q-point: ( 0.27  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.645   (   5.885    3.398    8.183)   10.636
   1.848   (   0.249    0.144   13.064)   13.067
   1.876   (   4.750    2.742   -3.538)    6.527
   2.326   (  -1.816   -1.048  -10.561)   10.767
   3.752   (  23.791   13.736   15.987)   31.785
   4.345   (  15.477    8.936  -14.393)   22.946
   6.272   (  -1.997   -1.153    4.558)    5.109
   6.272   (   0.823    0.475    0.167)    0.964
   6.297   (  -0.054   -0.031   -1.067)    1.069
   6.413   (  -0.420   -0.242   -1.390)    1.472
   6.820   (  -2.629   -1.518    4.907)    5.770
   7.006   (  -3.825   -2.209   -4.365)    6.210
======================= Grid point 728 (87/135) =======================
q-point: ( 0.33  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.743   (   2.895    1.671    7.349)    8.073
   1.852   (   0.281    0.163    7.064)    7.072
   1.951   (   1.943    1.122   -3.204)    3.911
   2.242   (  -4.876   -2.815  -11.539)   12.839
   4.261   (  21.466   12.393   13.722)   28.331
   4.714   (  16.711    9.648   -9.874)   21.676
   6.229   (  -1.810   -1.045    6.252)    6.592
   6.292   (   0.958    0.553   -0.202)    1.124
   6.303   (   0.569    0.328   -0.686)    0.949
   6.415   (   0.650    0.375   -2.128)    2.257
   6.728   (  -5.558   -3.209    3.314)    7.223
   6.889   (  -6.544   -3.778   -3.836)    8.474
======================= Grid point 729 (88/135) =======================
q-point: ( 0.40  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.782   (   0.846    0.488    6.739)    6.810
   1.865   (   0.943    0.544    1.094)    1.544
   1.969   (  -0.120   -0.070   -2.712)    2.715
   2.133   (  -4.584   -2.647  -11.346)   12.520
   4.698   (  17.159    9.907   12.571)   23.465
   5.059   (  13.324    7.693   -6.503)   16.703
   6.191   (  -1.421   -0.820    7.669)    7.842
   6.313   (   0.838    0.484   -0.311)    1.016
   6.319   (   0.780    0.450   -0.620)    1.094
   6.437   (   0.943    0.544   -4.200)    4.339
   6.580   (  -7.216   -4.166    2.824)    8.798
   6.730   (  -7.127   -4.115   -3.205)    8.832
======================= Grid point 730 (89/135) =======================
q-point: ( 0.47  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.791   (   0.101    0.058    6.328)    6.329
   1.890   (   0.945    0.546   -3.047)    3.236
   1.960   (  -0.386   -0.223   -2.233)    2.277
   2.051   (  -2.273   -1.312  -10.612)   10.931
   4.994   (   8.061    4.654   11.798)   15.027
   5.268   (   4.626    2.671   -2.876)    6.066
   6.190   (   1.790    1.033    6.336)    6.664
   6.328   (   0.426    0.246   -0.397)    0.632
   6.332   (   0.276    0.159   -0.593)    0.673
   6.366   (  -6.559   -3.787   -3.239)    8.238
   6.515   (   1.808    1.044    0.800)    2.236
   6.597   (  -4.429   -2.557   -2.690)    5.779
======================= Grid point 739 (90/135) =======================
q-point: ( 0.07  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.256   (   6.032   10.447   11.417)   16.610
   1.536   (   6.398   11.082    4.289)   13.496
   1.597   (   6.623   11.471    3.041)   13.590
   1.880   (   9.753   16.893   -5.139)   20.172
   2.735   (   5.103    8.839   35.177)   36.628
   3.977   (  -0.648   -1.122  -28.103)   28.133
   6.260   (   0.167    0.289    0.901)    0.961
   6.291   (   0.065    0.112    0.300)    0.327
   6.343   (   1.048    1.816   -0.643)    2.193
   6.386   (   1.034    1.791   -1.322)    2.454
   6.876   (  -0.770   -1.334    8.577)    8.715
   7.129   (  -1.811   -3.137   -4.802)    6.015
======================= Grid point 740 (91/135) =======================
q-point: ( 0.13  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.474   (   5.180   13.434    9.644)   17.329
   1.734   (   5.456    9.996   -0.773)   11.414
   1.774   (   4.122    3.910   12.376)   13.618
   2.178   (   9.520    7.598   -6.731)   13.916
   3.033   (  15.501   15.904   25.054)   33.480
   4.006   (   3.827    3.715  -23.493)   24.091
   6.266   (  -0.194    1.301    0.801)    1.541
   6.294   (   0.379   -0.595    0.581)    0.914
   6.352   (  -1.825    1.502    0.127)    2.367
   6.411   (   1.180   -0.809   -0.900)    1.691
   6.854   (  -0.497   -1.039    7.020)    7.114
   7.073   (  -1.431   -2.805   -4.662)    5.626
======================= Grid point 741 (92/135) =======================
q-point: ( 0.20  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.663   (   1.929   12.134    8.586)   14.990
   1.831   (   1.846    0.540   12.809)   12.953
   1.907   (   1.904   11.695   -2.117)   12.037
   2.313   (   2.553   -0.314   -9.556)    9.896
   3.508   (  21.216   15.906   17.991)   32.044
   4.197   (  11.592    8.494  -17.022)   22.277
   6.274   (  -0.707   -0.310    1.995)    2.139
   6.286   (  -0.595    0.871    1.442)    1.787
   6.331   (  -1.857    1.988   -0.424)    2.753
   6.411   (   0.137   -0.987   -1.155)    1.525
   6.829   (  -1.019   -2.313    5.332)    5.900
   7.017   (  -1.643   -4.135   -4.364)    6.233
======================= Grid point 742 (93/135) =======================
q-point: ( 0.27  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.760   (   0.391    4.474    7.952)    9.132
   1.873   (  -0.414    5.936    7.558)    9.620
   2.018   (  -1.762   10.637   -2.402)   11.046
   2.288   (  -2.997   -2.169  -10.875)   11.487
   4.016   (  21.865   13.361   14.722)   29.552
   4.526   (  16.075    9.523  -11.698)   22.044
   6.231   (  -1.676   -2.590    5.002)    5.876
   6.301   (  -0.139    1.954   -0.190)    1.968
   6.334   (  -1.108    3.280   -0.287)    3.474
   6.408   (   0.357    0.004   -1.471)    1.514
   6.764   (  -3.106   -4.188    3.639)    6.358
   6.926   (  -3.580   -5.965   -3.958)    8.004
======================= Grid point 743 (94/135) =======================
q-point: ( 0.33  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.791   (   0.830    0.918    4.448)    4.617
   1.916   (  -2.801    8.054    4.197)    9.504
   2.055   (  -4.429    8.100   -3.049)    9.723
   2.199   (  -5.403   -0.406  -10.053)   11.420
   4.484   (  19.656   10.513   13.135)   25.873
   4.886   (  15.773    7.940   -7.997)   19.385
   6.177   (  -0.972   -3.742    6.753)    7.781
   6.313   (   0.073    0.351   -1.125)    1.180
   6.352   (  -1.242    4.227    0.193)    4.409
   6.427   (   1.264    1.060   -1.781)    2.427
   6.645   (  -5.590   -4.882    1.980)    7.682
   6.785   (  -5.509   -6.542   -3.518)    9.248
======================= Grid point 744 (95/135) =======================
q-point: ( 0.40  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.812   (   0.540    0.938    1.369)    1.745
   1.941   (  -3.048    7.803    1.466)    8.504
   2.035   (  -4.777    4.060   -5.147)    8.112
   2.117   (  -5.817    3.706   -6.847)    9.719
   4.858   (  14.144    6.616   12.177)   19.801
   5.172   (  10.210    3.889   -4.482)   11.809
   6.136   (   1.277   -4.099    6.932)    8.154
   6.300   (  -0.485   -3.302   -2.070)    3.928
   6.371   (  -1.387    4.358    0.423)    4.593
   6.463   (  -0.303    2.279   -4.079)    4.683
   6.502   (  -4.528   -3.150    3.202)    6.378
   6.639   (  -4.710   -4.613   -3.178)    7.319
======================= Grid point 745 (96/135) =======================
q-point: ( 0.47  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.823   (  -0.595    1.031    0.589)    1.328
   1.965   (  -4.986    8.636    1.169)   10.040
   1.997   (  -0.579    1.002   -8.325)    8.405
   2.086   (  -4.282    7.417   -4.316)    9.590
   5.027   (   0.751   -1.301   11.640)   11.736
   5.274   (   1.025   -1.775   -1.728)    2.681
   6.148   (   2.686   -4.653    5.161)    7.450
   6.246   (   2.223   -3.850   -2.259)    4.987
   6.386   (  -2.248    3.894    0.656)    4.544
   6.416   (  -2.563    4.439   -0.302)    5.135
   6.516   (   0.897   -1.554   -0.268)    1.814
   6.559   (   0.116   -0.200   -3.516)    3.524
======================= Grid point 756 (97/135) =======================
q-point: ( 0.13  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.747   (   7.637   13.228    7.781)   17.142
   1.790   (   0.065    0.112   14.495)   14.495
   1.971   (   6.806   11.789   -2.791)   13.896
   2.273   (   1.358    2.352   -9.366)    9.752
   3.397   (  12.171   21.080   19.365)   31.105
   4.139   (   5.854   10.140  -18.238)   21.673
   6.280   (  -0.677   -1.173    2.710)    3.030
   6.299   (   1.040    1.802    0.717)    2.201
   6.379   (   0.572    0.990   -1.308)    1.737
   6.386   (  -0.592   -1.026   -0.419)    1.256
   6.824   (  -1.285   -2.226    5.305)    5.895
   7.005   (  -2.399   -4.156   -4.193)    6.373
======================= Grid point 757 (98/135) =======================
q-point: ( 0.20  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.778   (  -0.428   -0.262   10.081)   10.093
   1.960   (   2.369   11.255    5.224)   12.632
   2.160   (   0.547   12.980   -2.329)   13.199
   2.291   (   0.859   -0.761   -9.787)    9.854
   3.828   (  16.636   17.122   15.947)   28.709
   4.390   (  11.216   11.357  -13.155)   20.684
   6.239   (  -1.974   -2.562    4.919)    5.887
   6.325   (  -0.314    1.342    0.132)    1.384
   6.372   (  -0.104    1.606   -0.398)    1.658
   6.408   (   0.536    1.698   -1.645)    2.424
   6.764   (  -2.177   -4.188    3.433)    5.836
   6.911   (  -2.690   -6.521   -3.738)    7.984
======================= Grid point 758 (99/135) =======================
q-point: ( 0.27  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.776   (   0.710   -0.432    4.914)    4.984
   2.041   (  -1.998    5.635   -1.596)    6.188
   2.228   (  -4.621    9.491   -2.678)   10.891
   2.310   (  -2.876    7.321   -4.342)    8.984
   4.270   (  18.207   13.043   13.956)   26.389
   4.711   (  14.158    9.700   -9.214)   19.479
   6.169   (  -2.432   -3.944    6.972)    8.371
   6.318   (  -0.594   -0.648   -1.355)    1.615
   6.401   (   0.379    2.774    0.168)    2.805
   6.429   (  -0.369    2.237   -1.145)    2.540
   6.674   (  -3.364   -4.393    1.414)    5.711
   6.791   (  -3.502   -7.242   -3.481)    8.765
======================= Grid point 759 (100/135) =======================
q-point: ( 0.33  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.793   (   1.496   -0.046    0.079)    1.499
   2.029   (  -2.092    0.303   -8.241)    8.507
   2.220   (  -7.670   10.986    0.881)   13.427
   2.313   (  -6.954   10.537   -2.036)   12.788
   4.672   (  16.242    9.196   12.741)   22.599
   5.023   (  12.609    6.377   -5.833)   15.287
   6.089   (  -1.305   -5.444    8.199)    9.928
   6.286   (  -1.482   -3.409   -2.477)    4.467
   6.425   (  -0.237    1.606   -0.028)    1.623
   6.441   (   0.563   -0.110   -0.366)    0.681
   6.580   (  -4.360   -0.456   -0.202)    4.389
   6.668   (  -3.910   -4.236   -3.249)    6.617
======================= Grid point 760 (101/135) =======================
q-point: ( 0.40  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.818   (   1.078    0.062   -3.092)    3.275
   1.996   (  -0.617   -1.482  -11.126)   11.241
   2.199   (  -8.236   12.833    3.048)   15.551
   2.295   (  -7.453   11.489   -1.394)   13.765
   4.956   (   8.362    3.608   11.822)   14.923
   5.214   (   5.072    0.737   -1.920)    5.473
   6.046   (   2.833   -5.072    6.539)    8.747
   6.199   (  -1.875   -6.186   -2.456)    6.915
   6.430   (   0.024    0.183   -0.659)    0.684
   6.442   (   1.018   -2.164   -0.131)    2.396
   6.537   (  -2.975    3.976   -0.579)    5.000
   6.593   (  -3.073    1.166   -3.009)    4.456
======================= Grid point 772 (102/135) =======================
q-point: ( 0.20  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.767   (  -0.168   -0.291    5.423)    5.434
   2.082   (   0.368    0.637   -5.407)    5.457
   2.328   (   3.055    5.292   -2.265)    6.517
   2.380   (   4.789    8.295   -1.302)    9.666
   4.172   (  10.580   18.326   14.487)   25.645
   4.633   (   7.867   13.627   -9.789)   18.531
   6.173   (  -2.530   -4.382    7.291)    8.874
   6.318   (  -0.771   -1.335   -1.060)    1.871
   6.408   (   0.932    1.615    0.230)    1.879
   6.451   (   0.963    1.667   -1.794)    2.631
   6.680   (  -2.326   -4.029    1.184)    4.800
   6.781   (  -3.676   -6.367   -3.276)    8.049
======================= Grid point 773 (103/135) =======================
q-point: ( 0.27  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.773   (   0.872    0.331    0.480)    1.049
   2.050   (  -1.845   -2.291  -11.024)   11.410
   2.412   (  -1.940    8.603    1.942)    9.030
   2.486   (  -1.393    7.960   -0.935)    8.135
   4.529   (  12.991   13.742   13.313)   23.126
   4.906   (  10.078   10.353   -6.692)   15.923
   6.072   (  -3.163   -6.149    9.301)   11.590
   6.279   (  -1.701   -3.449   -2.001)    4.335
   6.412   (  -0.394   -1.656   -0.174)    1.711
   6.456   (  -0.855   -0.575   -2.244)    2.469
   6.626   (  -1.469   -0.185   -0.060)    1.482
   6.687   (  -2.059   -2.385   -2.680)    4.136
======================= Grid point 774 (104/135) =======================
q-point: ( 0.33  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.800   (   1.406    0.905   -3.760)    4.115
   1.994   (  -1.571   -2.706  -12.158)   12.554
   2.433   (  -5.467    9.491    2.757)   11.295
   2.505   (  -4.942    8.240   -0.677)    9.632
   4.848   (  10.835    9.120   12.255)   18.729
   5.136   (   7.011    5.283   -2.917)    9.250
   5.971   (  -0.598   -6.353    8.646)   10.746
   6.179   (  -3.421   -7.295   -2.311)    8.382
   6.392   (   0.832   -2.220   -1.192)    2.654
   6.429   (  -0.450   -2.064   -2.853)    3.550
   6.615   (  -1.813    2.776    0.110)    3.318
   6.658   (  -1.916    2.318   -1.972)    3.596
======================= Grid point 775 (105/135) =======================
q-point: ( 0.40  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.822   (  -0.217    0.376   -5.719)    5.735
   1.962   (   0.953   -1.650  -12.048)   12.198
   2.431   (  -5.821   10.083    2.906)   12.000
   2.500   (  -5.111    8.853   -0.568)   10.238
   5.005   (  -0.949    1.643   11.457)   11.613
   5.215   (   0.163   -0.282    0.488)    0.587
   5.946   (   2.861   -4.956    6.184)    8.426
   6.090   (   2.736   -4.739   -1.903)    5.793
   6.392   (   1.255   -2.173   -1.580)    2.965
   6.411   (   0.879   -1.523   -3.087)    3.553
   6.613   (  -2.035    3.524    0.323)    4.082
   6.655   (  -2.111    3.656   -1.761)    4.574
======================= Grid point 789 (106/135) =======================
q-point: ( 0.27  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.793   (   0.945    1.638   -3.866)    4.304
   1.995   (  -1.797   -3.113  -12.536)   13.041
   2.530   (   2.064    3.575    2.490)    4.821
   2.590   (   1.732    2.999   -0.390)    3.485
   4.797   (   7.859   13.612   12.496)   20.080
   5.100   (   5.194    8.996   -3.432)   10.940
   5.945   (  -3.563   -6.172    9.534)   11.903
   6.172   (  -4.305   -7.456   -2.096)    8.861
   6.374   (  -0.796   -1.379   -1.622)    2.273
   6.421   (  -1.208   -2.092   -3.672)    4.395
   6.644   (   0.743    1.288    0.330)    1.523
   6.687   (   0.696    1.205   -1.710)    2.204
======================= Grid point 790 (107/135) =======================
q-point: ( 0.33  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.830   (   0.851    2.170   -7.346)    7.707
   1.938   (  -0.902   -2.987  -12.024)   12.422
   2.566   (  -1.727    4.251    2.501)    5.226
   2.619   (  -1.568    3.551   -0.151)    3.884
   5.025   (   4.559    8.878   11.148)   14.963
   5.211   (   1.016    1.564    1.983)    2.722
   5.873   (   0.871   -2.691    5.806)    6.459
   6.025   (  -2.589   -8.058   -1.498)    8.595
   6.364   (   0.473   -0.640   -2.536)    2.658
   6.389   (  -0.245   -1.732   -3.781)    4.166
   6.662   (  -0.503    1.819    0.728)    2.023
   6.707   (  -0.499    2.102   -1.486)    2.622
======================= Grid point 805 (108/135) =======================
q-point: ( 0.33  0.33  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.881   (  -1.442   -2.498  -10.284)   10.681
   1.881   (   1.442    2.498  -10.284)   10.681
   2.605   (  -0.000   -0.000    2.378)    2.378
   2.651   (  -0.000   -0.000    0.034)    0.034
   5.172   (  -2.743   -4.751    8.069)    9.757
   5.172   (   2.743    4.751    8.069)    9.757
   5.894   (  -2.327   -4.030    0.879)    4.736
   5.894   (   2.327    4.030    0.879)    4.736
   6.365   (  -0.356   -0.617   -3.429)    3.503
   6.365   (   0.356    0.617   -3.429)    3.503
   6.680   (  -0.000   -0.000    0.945)    0.945
   6.729   (  -0.000   -0.000   -1.435)    1.435
======================= Grid point 964 (109/135) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 81
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.264   (   0.000   -0.000   -9.516)    9.516
   1.264   (   0.000    0.000   -9.516)    9.516
   1.264   (   0.000    0.000    9.516)    9.516
   1.264   (  -0.000   -0.000    9.516)    9.516
   3.391   (   0.000   -0.000  -36.471)   36.471
   3.391   (   0.000   -0.000   36.471)   36.471
   6.281   (   0.000    0.000   -1.361)    1.361
   6.281   (   0.000   -0.000   -1.361)    1.361
   6.281   (   0.000    0.000    1.361)    1.361
   6.281   (   0.000    0.000    1.361)    1.361
   7.089   (   0.000   -0.000   -7.925)    7.925
   7.089   (   0.000   -0.000    7.925)    7.925
======================= Grid point 966 (110/135) =======================
q-point: ( 0.07  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.323   (   5.000    2.887   -9.079)   10.760
   1.323   (   5.000    2.887    9.079)   10.760
   1.479   (  16.650    9.613   -8.334)   20.955
   1.479   (  16.650    9.613    8.334)   20.955
   3.379   (  -0.448   -0.259  -35.525)   35.529
   3.379   (  -0.448   -0.259   35.525)   35.529
   6.280   (  -0.105   -0.060   -1.493)    1.498
   6.280   (  -0.105   -0.060    1.493)    1.498
   6.314   (   2.539    1.466   -1.761)    3.420
   6.314   (   2.539    1.466    1.761)    3.420
   7.061   (  -2.184   -1.261   -7.373)    7.793
   7.061   (  -2.184   -1.261    7.373)    7.793
======================= Grid point 968 (111/135) =======================
q-point: ( 0.13  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.469   (   7.602    4.389   -8.039)   11.903
   1.469   (   7.602    4.389    8.039)   11.903
   1.868   (  15.419    8.902   -8.542)   19.747
   1.868   (  15.419    8.902    8.542)   19.747
   3.417   (   5.089    2.938  -31.118)   31.668
   3.417   (   5.089    2.938   31.118)   31.668
   6.277   (  -0.128   -0.074   -1.541)    1.548
   6.277   (  -0.128   -0.074    1.541)    1.548
   6.367   (   1.622    0.936   -2.481)    3.109
   6.367   (   1.622    0.936    2.481)    3.109
   7.009   (  -2.084   -1.203   -6.325)    6.767
   7.009   (  -2.084   -1.203    6.325)    6.767
======================= Grid point 970 (112/135) =======================
q-point: ( 0.20  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.641   (   7.371    4.255   -6.937)   10.979
   1.641   (   7.371    4.255    6.937)   10.979
   2.092   (   4.349    2.511  -11.825)   12.847
   2.092   (   4.349    2.511   11.825)   12.847
   3.647   (  15.341    8.857  -22.561)   28.685
   3.647   (  15.341    8.857   22.561)   28.685
   6.275   (   0.030    0.017   -1.214)    1.215
   6.275   (   0.030    0.017    1.214)    1.215
   6.376   (  -0.628   -0.363   -3.272)    3.351
   6.376   (  -0.628   -0.363    3.272)    3.351
   6.971   (  -1.612   -0.931   -5.525)    5.830
   6.971   (  -1.612   -0.931    5.525)    5.830
======================= Grid point 972 (113/135) =======================
q-point: ( 0.27  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.784   (   5.163    2.981   -6.121)    8.545
   1.784   (   5.163    2.981    6.121)    8.545
   2.101   (  -2.365   -1.365  -12.830)   13.118
   2.101   (  -2.365   -1.365   12.830)   13.118
   4.057   (  20.160   11.639  -15.649)   28.050
   4.057   (  20.160   11.639   15.649)   28.050
   6.280   (   0.391    0.226   -0.676)    0.813
   6.280   (   0.391    0.226    0.676)    0.813
   6.354   (  -1.038   -0.599   -4.075)    4.248
   6.354   (  -1.038   -0.599    4.075)    4.248
   6.919   (  -3.431   -1.981   -4.939)    6.331
   6.919   (  -3.431   -1.981    4.939)    6.331
======================= Grid point 974 (114/135) =======================
q-point: ( 0.33  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.868   (   2.354    1.359   -5.531)    6.163
   1.868   (   2.354    1.359    5.531)    6.163
   2.026   (  -3.782   -2.184  -10.667)   11.527
   2.026   (  -3.782   -2.184   10.667)   11.527
   4.508   (  19.582   11.305  -11.844)   25.525
   4.508   (  19.582   11.305   11.844)   25.525
   6.293   (   0.760    0.439   -0.274)    0.920
   6.293   (   0.760    0.439    0.274)    0.920
   6.337   (  -0.521   -0.301   -5.196)    5.231
   6.337   (  -0.521   -0.301    5.196)    5.231
   6.811   (  -6.157   -3.555   -4.485)    8.406
   6.811   (  -6.157   -3.555    4.485)    8.406
======================= Grid point 976 (115/135) =======================
q-point: ( 0.40  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.896   (   0.368    0.213   -4.957)    4.975
   1.896   (   0.368    0.213    4.957)    4.975
   1.948   (  -2.931   -1.692   -7.392)    8.130
   1.948   (  -2.931   -1.692    7.392)    8.130
   4.910   (  15.811    9.129   -9.384)   20.528
   4.910   (  15.811    9.129    9.384)   20.528
   6.311   (   0.794    0.459   -0.187)    0.936
   6.311   (   0.794    0.459    0.187)    0.936
   6.323   (  -1.010   -0.583   -6.582)    6.684
   6.323   (  -1.010   -0.583    6.582)    6.684
   6.660   (  -6.662   -3.846   -4.225)    8.776
   6.660   (  -6.662   -3.846    4.225)    8.776
======================= Grid point 978 (116/135) =======================
q-point: ( 0.47  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.896   (  -0.128   -0.074   -4.487)    4.489
   1.896   (  -0.128   -0.074    4.487)    4.489
   1.902   (  -1.088   -0.628   -4.464)    4.637
   1.902   (  -1.088   -0.628    4.464)    4.637
   5.177   (   6.945    4.009   -7.032)   10.665
   5.177   (   6.945    4.009    7.032)   10.665
   6.287   (  -1.698   -0.980   -4.467)    4.878
   6.287   (  -1.698   -0.980    4.467)    4.878
   6.325   (   0.342    0.197   -0.112)    0.410
   6.325   (   0.342    0.197    0.112)    0.410
   6.551   (  -2.485   -1.435   -2.347)    3.708
   6.551   (  -2.485   -1.435    2.347)    3.708
======================= Grid point 996 (117/135) =======================
q-point: ( 0.07  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.444   (   5.176    8.965   -8.034)   13.103
   1.444   (   5.176    8.965    8.034)   13.103
   1.756   (   9.368   16.226   -7.944)   20.351
   1.756   (   9.368   16.226    7.944)   20.351
   3.392   (   1.530    2.650  -32.835)   32.977
   3.392   (   1.530    2.650   32.835)   32.977
   6.281   (   0.089    0.154   -1.130)    1.144
   6.281   (   0.089    0.154    1.130)    1.144
   6.354   (   1.209    2.095   -1.607)    2.904
   6.354   (   1.209    2.095    1.607)    2.904
   7.022   (  -1.429   -2.476   -6.566)    7.161
   7.022   (  -1.429   -2.476    6.566)    7.161
======================= Grid point 998 (118/135) =======================
q-point: ( 0.13  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.635   (   4.619   11.953   -7.003)   14.603
   1.635   (   4.619   11.953    7.003)   14.603
   2.015   (   7.496    5.139  -10.472)   13.866
   2.015   (   7.496    5.139   10.472)   13.866
   3.529   (   8.875    9.141  -26.003)   28.957
   3.529   (   8.875    9.141   26.003)   28.957
   6.285   (  -0.169    0.976   -1.067)    1.455
   6.285   (  -0.169    0.976    1.067)    1.455
   6.377   (   0.150   -0.028   -2.083)    2.088
   6.377   (   0.150   -0.028    2.083)    2.088
   6.976   (  -1.191   -2.202   -5.621)    6.153
   6.976   (  -1.191   -2.202    5.621)    6.153
======================= Grid point 1000 (119/135) =======================
q-point: ( 0.20  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.806   (   1.591   11.575   -6.174)   13.215
   1.806   (   1.591   11.575    6.174)   13.215
   2.100   (   1.324   -1.020  -12.560)   12.670
   2.100   (   1.324   -1.020   12.560)   12.670
   3.858   (  16.645   12.587  -18.259)   27.729
   3.858   (  16.645   12.587   18.259)   27.729
   6.293   (  -0.539    1.607   -0.779)    1.865
   6.293   (  -0.539    1.607    0.779)    1.865
   6.365   (  -0.806   -0.623   -3.049)    3.215
   6.365   (  -0.806   -0.623    3.049)    3.215
   6.930   (  -1.516   -3.464   -4.923)    6.208
   6.930   (  -1.516   -3.464    4.923)    6.208
======================= Grid point 1002 (120/135) =======================
q-point: ( 0.27  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.900   (  -0.875    6.164   -5.554)    8.343
   1.900   (  -0.875    6.164    5.554)    8.343
   2.082   (  -2.874    1.818   -9.774)   10.349
   2.082   (  -2.874    1.818    9.774)   10.349
   4.287   (  19.388   11.871  -13.371)   26.375
   4.287   (  19.388   11.871   13.371)   26.375
   6.292   (  -0.696    0.096   -1.231)    1.418
   6.292   (  -0.696    0.096    1.231)    1.418
   6.356   (  -0.470    1.511   -2.902)    3.305
   6.356   (  -0.470    1.511    2.902)    3.305
   6.847   (  -3.477   -5.394   -4.417)    7.790
   6.847   (  -3.477   -5.394    4.417)    7.790
======================= Grid point 1004 (121/135) =======================
q-point: ( 0.33  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.904   (  -1.189   -0.526   -6.270)    6.403
   1.904   (  -1.189   -0.526    6.270)    6.403
   2.056   (  -5.797    7.631   -4.843)   10.737
   2.056   (  -5.797    7.631    4.843)   10.737
   4.711   (  18.144    9.625  -10.469)   23.053
   4.711   (  18.144    9.625   10.469)   23.053
   6.274   (  -0.125   -2.240   -3.380)    4.057
   6.274   (  -0.125   -2.240    3.380)    4.057
   6.373   (  -0.616    3.510   -2.402)    4.298
   6.373   (  -0.616    3.510    2.402)    4.298
   6.713   (  -5.362   -6.137   -4.066)    9.108
   6.713   (  -5.362   -6.137    4.066)    9.108
======================= Grid point 1006 (122/135) =======================
q-point: ( 0.40  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.873   (  -0.482   -2.008   -5.569)    5.940
   1.873   (  -0.482   -2.008    5.569)    5.940
   2.033   (  -6.601    9.375   -3.584)   12.012
   2.033   (  -6.601    9.375    3.584)   12.012
   5.055   (  12.803    5.770   -8.121)   16.222
   5.055   (  12.803    5.770    8.121)   16.222
   6.240   (  -0.036   -4.613   -4.264)    6.282
   6.240   (  -0.036   -4.613    4.264)    6.282
   6.395   (  -1.490    4.386   -2.287)    5.166
   6.395   (  -1.490    4.386    2.287)    5.166
   6.575   (  -4.007   -4.371   -3.709)    6.995
   6.575   (  -4.007   -4.371    3.709)    6.995
======================= Grid point 1008 (123/135) =======================
q-point: ( 0.47  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.858   (   0.816   -1.413   -4.584)    4.865
   1.858   (   0.816   -1.413    4.584)    4.865
   2.022   (  -5.787   10.023   -3.441)   12.075
   2.022   (  -5.787   10.023    3.441)   12.075
   5.202   (   0.813   -1.409   -6.307)    6.513
   5.202   (   0.813   -1.409    6.307)    6.513
   6.207   (   2.434   -4.215   -1.925)    5.234
   6.207   (   2.434   -4.215    1.925)    5.234
   6.405   (  -2.710    4.693   -0.888)    5.491
   6.405   (  -2.710    4.693    0.888)    5.491
   6.519   (   0.870   -1.507   -0.750)    1.895
   6.519   (   0.870   -1.507    0.750)    1.895
======================= Grid point 1030 (124/135) =======================
q-point: ( 0.13  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.877   (   6.442   11.157   -5.411)   13.974
   1.877   (   6.442   11.157    5.411)   13.974
   2.065   (   0.601    1.040  -12.253)   12.312
   2.065   (   0.601    1.040   12.253)   12.312
   3.774   (   9.166   15.876  -19.656)   26.878
   3.774   (   9.166   15.876   19.656)   26.878
   6.309   (   0.670    1.160   -0.128)    1.346
   6.309   (   0.670    1.160    0.128)    1.346
   6.368   (  -0.346   -0.599   -0.618)    0.928
   6.368   (  -0.346   -0.599    0.618)    0.928
   6.922   (  -2.018   -3.495   -4.708)    6.201
   6.922   (  -2.018   -3.495    4.708)    6.201
======================= Grid point 1032 (125/135) =======================
q-point: ( 0.20  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.965   (  -2.004    0.516   -7.170)    7.462
   1.965   (  -2.004    0.516    7.170)    7.462
   2.147   (   2.556    9.330   -3.629)   10.332
   2.147   (   2.556    9.330    3.629)   10.332
   4.124   (  14.229   14.754  -14.758)   25.258
   4.124   (  14.229   14.754   14.758)   25.258
   6.308   (  -1.162   -0.790   -2.286)    2.684
   6.308   (  -1.162   -0.790    2.286)    2.684
   6.376   (   0.421    2.226   -1.233)    2.579
   6.376   (   0.421    2.226    1.233)    2.579
   6.838   (  -2.634   -5.783   -3.956)    7.485
   6.838   (  -2.634   -5.783    3.956)    7.485
======================= Grid point 1034 (126/135) =======================
q-point: ( 0.27  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.914   (  -1.887   -2.618   -8.703)    9.282
   1.914   (  -1.887   -2.618    8.703)    9.282
   2.246   (  -3.545   12.044   -2.924)   12.891
   2.246   (  -3.545   12.044    2.924)   12.891
   4.515   (  16.512   11.842  -11.494)   23.345
   4.515   (  16.512   11.842   11.494)   23.345
   6.270   (  -1.540   -2.087   -3.785)    4.588
   6.270   (  -1.540   -2.087    3.785)    4.588
   6.408   (   0.154    2.841   -0.742)    2.940
   6.408   (   0.154    2.841    0.742)    2.940
   6.725   (  -3.605   -6.434   -3.376)    8.111
   6.725   (  -3.605   -6.434    3.376)    8.111
======================= Grid point 1036 (127/135) =======================
q-point: ( 0.33  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.864   (  -1.000   -2.616   -7.023)    7.561
   1.864   (  -1.000   -2.616    7.023)    7.561
   2.268   (  -7.296   12.216   -2.800)   14.501
   2.268   (  -7.296   12.216    2.800)   14.501
   4.883   (  14.890    8.308   -9.089)   19.322
   4.883   (  14.890    8.308    9.089)   19.322
   6.213   (  -1.744   -4.201   -5.060)    6.804
   6.213   (  -1.744   -4.201    5.060)    6.804
   6.430   (   0.227    0.775   -0.560)    0.983
   6.430   (   0.227    0.775    0.560)    0.983
   6.609   (  -4.124   -2.839   -2.710)    5.693
   6.609   (  -4.124   -2.839    2.710)    5.693
======================= Grid point 1038 (128/135) =======================
q-point: ( 0.40  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.834   (   0.175   -1.936   -4.961)    5.328
   1.834   (   0.175   -1.936    4.961)    5.328
   2.257   (  -7.822   12.512   -2.647)   14.991
   2.257   (  -7.822   12.512    2.647)   14.991
   5.135   (   7.223    2.748   -6.660)   10.202
   5.135   (   7.223    2.748    6.660)   10.202
   6.141   (  -0.164   -5.450   -3.719)    6.600
   6.141   (  -0.164   -5.450    3.719)    6.600
   6.432   (   1.006   -2.071   -0.696)    2.405
   6.432   (   1.006   -2.071    0.696)    2.405
   6.548   (  -3.337    3.054   -1.601)    4.798
   6.548   (  -3.337    3.054    1.601)    4.798
======================= Grid point 1062 (129/135) =======================
q-point: ( 0.20  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.917   (  -1.954   -3.385   -9.506)   10.278
   1.917   (  -1.954   -3.385    9.506)   10.278
   2.345   (   4.993    8.648   -2.141)   10.213
   2.345   (   4.993    8.648    2.141)   10.213
   4.427   (   9.508   16.468  -12.026)   22.500
   4.427   (   9.508   16.468   12.026)   22.500
   6.275   (  -1.402   -2.428   -3.597)    4.561
   6.275   (  -1.402   -2.428    3.597)    4.561
   6.421   (   0.951    1.647   -1.199)    2.248
   6.421   (   0.951    1.647    1.199)    2.248
   6.722   (  -3.287   -5.694   -2.856)    7.168
   6.722   (  -3.287   -5.694    2.856)    7.168
======================= Grid point 1064 (130/135) =======================
q-point: ( 0.27  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.857   (  -1.438   -2.857   -7.948)    8.568
   1.857   (  -1.438   -2.857    7.948)    8.568
   2.456   (  -1.518    8.704   -2.111)    9.084
   2.456   (  -1.518    8.704    2.111)    9.084
   4.751   (  11.866   12.561   -9.768)   19.849
   4.751   (  11.866   12.561    9.768)   19.849
   6.210   (  -2.601   -4.341   -5.288)    7.319
   6.210   (  -2.601   -4.341    5.288)    7.319
   6.422   (  -0.571   -2.005   -1.182)    2.397
   6.422   (  -0.571   -2.005    1.182)    2.397
   6.644   (  -1.791   -1.037   -1.747)    2.708
   6.644   (  -1.791   -1.037    1.747)    2.708
======================= Grid point 1066 (131/135) =======================
q-point: ( 0.33  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.817   (  -0.488   -2.079   -5.519)    5.918
   1.817   (  -0.488   -2.079    5.519)    5.918
   2.480   (  -5.149    8.968   -1.946)   10.523
   2.480   (  -5.149    8.968    1.946)   10.523
   5.039   (   9.513    7.922   -7.494)   14.471
   5.039   (   9.513    7.922    7.494)   14.471
   6.106   (  -2.641   -6.738   -5.438)    9.053
   6.106   (  -2.641   -6.738    5.438)    9.053
   6.389   (   0.345   -3.200   -1.082)    3.396
   6.389   (   0.345   -3.200    1.082)    3.396
   6.627   (  -1.889    3.003   -1.181)    3.739
   6.627   (  -1.889    3.003    1.181)    3.739
======================= Grid point 1068 (132/135) =======================
q-point: ( 0.40  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.802   (   0.727   -1.260   -3.941)    4.201
   1.802   (   0.727   -1.260    3.941)    4.201
   2.477   (  -5.468    9.471   -1.859)   11.093
   2.477   (  -5.468    9.471    1.859)   11.093
   5.170   (  -0.621    1.075   -5.469)    5.608
   5.170   (  -0.621    1.075    5.469)    5.608
   6.039   (   2.720   -4.711   -3.564)    6.503
   6.039   (   2.720   -4.711    3.564)    6.503
   6.377   (   1.490   -2.582   -0.326)    2.999
   6.377   (   1.490   -2.582    0.326)    2.999
   6.627   (  -2.244    3.886   -1.105)    4.622
   6.627   (  -2.244    3.886    1.105)    4.622
======================= Grid point 1096 (133/135) =======================
q-point: ( 0.27  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.811   (  -1.081   -1.872   -5.779)    6.170
   1.811   (  -1.081   -1.872    5.779)    6.170
   2.571   (   1.944    3.367   -1.631)    4.216
   2.571   (   1.944    3.367    1.631)    4.216
   4.994   (   6.973   12.077   -7.880)   16.018
   4.994   (   6.973   12.077    7.880)   16.018
   6.096   (  -4.091   -7.086   -6.030)   10.165
   6.096   (  -4.091   -7.086    6.030)   10.165
   6.369   (  -1.320   -2.285   -1.376)    2.976
   6.369   (  -1.320   -2.285    1.376)    2.976
   6.659   (   0.866    1.499   -1.168)    2.089
   6.659   (   0.866    1.499    1.168)    2.089
======================= Grid point 1098 (134/135) =======================
q-point: ( 0.33  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.786   (  -0.136   -1.040   -3.076)    3.250
   1.786   (  -0.136   -1.040    3.076)    3.250
   2.604   (  -1.635    3.923   -1.389)    4.471
   2.604   (  -1.635    3.923    1.389)    4.471
   5.183   (   3.309    6.380   -4.987)    8.748
   5.183   (   3.309    6.380    4.987)    8.748
   5.972   (  -1.414   -6.246   -4.126)    7.618
   5.972   (  -1.414   -6.246    4.126)    7.618
   6.343   (   0.170   -1.523   -0.697)    1.684
   6.343   (   0.170   -1.523    0.697)    1.684
   6.681   (  -0.501    2.103   -1.218)    2.481
   6.681   (  -0.501    2.103    1.218)    2.481
======================= Grid point 1128 (135/135) =======================
q-point: ( 0.33  0.33 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 4.63e-04 4.63e-04 
Number of triplets: 119
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.776   (  -0.000   -0.000    0.000)    0.000
   1.776   (  -0.000   -0.000    0.000)    0.000
   2.640   (  -0.000   -0.000   -1.197)    1.197
   2.640   (  -0.000   -0.000    1.197)    1.197
   5.253   (  -0.000   -0.000    0.000)    0.000
   5.253   (  -0.000   -0.000    0.000)    0.000
   5.903   (  -0.000   -0.000    0.000)    0.000
   5.903   (  -0.000   -0.000    0.000)    0.000
   6.328   (  -0.000   -0.000    0.000)    0.000
   6.328   (  -0.000   -0.000    0.000)    0.000
   6.702   (  -0.000   -0.000   -1.305)    1.305
   6.702   (  -0.000   -0.000    1.305)    1.305
=================== End of collection of collisions ===================
----------- Thermal conductivity (W/m-k) with tetrahedron method -----------
#  T(K)        xx         yy         zz         yz         xz         xy        #ipm
    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/21600
   10.0    530.531    530.531    728.228      0.000     -0.000     -0.000 3/21600
   20.0    336.371    336.371    421.910     -0.000     -0.000     -0.000 3/21600
   30.0    266.497    266.497    316.940     -0.000     -0.000     -0.000 3/21600
   40.0    208.070    208.070    238.951     -0.000     -0.000     -0.000 3/21600
   50.0    157.976    157.976    177.101     -0.000     -0.000     -0.000 3/21600
   60.0    121.406    121.406    133.850     -0.000     -0.000     -0.000 3/21600
   70.0     96.405     96.405    105.109     -0.000     -0.000     -0.000 3/21600
   80.0     79.236     79.236     85.760     -0.000     -0.000     -0.000 3/21600
   90.0     67.075     67.075     72.249     -0.000     -0.000     -0.000 3/21600
  100.0     58.143     58.143     62.430     -0.000     -0.000     -0.000 3/21600
  110.0     51.354     51.354     55.025     -0.000     -0.000     -0.000 3/21600
  120.0     46.038     46.038     49.259     -0.000     -0.000     -0.000 3/21600
  130.0     41.768     41.768     44.648     -0.000     -0.000     -0.000 3/21600
  140.0     38.263     38.263     40.876     -0.000     -0.000     -0.000 3/21600
  150.0     35.333     35.333     37.731     -0.000     -0.000     -0.000 3/21600
  160.0     32.846     32.846     35.066     -0.000     -0.000     -0.000 3/21600
  170.0     30.706     30.706     32.776     -0.000     -0.000     -0.000 3/21600
  180.0     28.843     28.843     30.786     -0.000     -0.000     -0.000 3/21600
  190.0     27.207     27.207     29.039     -0.000     -0.000     -0.000 3/21600
  200.0     25.757     25.757     27.492     -0.000     -0.000     -0.000 3/21600
  210.0     24.461     24.461     26.110     -0.000     -0.000     -0.000 3/21600
  220.0     23.297     23.297     24.869     -0.000     -0.000     -0.000 3/21600
  230.0     22.243     22.243     23.746     -0.000     -0.000     -0.000 3/21600
  240.0     21.286     21.286     22.726     -0.000     -0.000     -0.000 3/21600
  250.0     20.411     20.411     21.794     -0.000     -0.000     -0.000 3/21600
  260.0     19.608     19.608     20.939     -0.000     -0.000     -0.000 3/21600
  270.0     18.869     18.869     20.152     -0.000     -0.000     -0.000 3/21600
  280.0     18.185     18.185     19.424     -0.000     -0.000     -0.000 3/21600
  290.0     17.551     17.551     18.749     -0.000     -0.000     -0.000 3/21600
  300.0     16.962     16.962     18.121     -0.000     -0.000     -0.000 3/21600
  310.0     16.412     16.412     17.535     -0.000     -0.000     -0.000 3/21600
  320.0     15.898     15.898     16.988     -0.000     -0.000     -0.000 3/21600
  330.0     15.416     15.416     16.474     -0.000     -0.000     -0.000 3/21600
  340.0     14.963     14.963     15.992     -0.000     -0.000     -0.000 3/21600
  350.0     14.537     14.537     15.538     -0.000     -0.000     -0.000 3/21600
  360.0     14.135     14.135     15.110     -0.000     -0.000     -0.000 3/21600
  370.0     13.755     13.755     14.705     -0.000     -0.000     -0.000 3/21600
  380.0     13.396     13.396     14.322     -0.000     -0.000     -0.000 3/21600
  390.0     13.055     13.055     13.960     -0.000     -0.000     -0.000 3/21600
  400.0     12.732     12.732     13.615     -0.000     -0.000     -0.000 3/21600
  410.0     12.425     12.425     13.287     -0.000     -0.000     -0.000 3/21600
  420.0     12.132     12.132     12.976     -0.000     -0.000     -0.000 3/21600
  430.0     11.853     11.853     12.678     -0.000     -0.000     -0.000 3/21600
  440.0     11.587     11.587     12.395     -0.000     -0.000     -0.000 3/21600
  450.0     11.333     11.333     12.124     -0.000     -0.000     -0.000 3/21600
  460.0     11.090     11.090     11.865     -0.000     -0.000     -0.000 3/21600
  470.0     10.857     10.857     11.617     -0.000     -0.000     -0.000 3/21600
  480.0     10.635     10.635     11.379     -0.000     -0.000     -0.000 3/21600
  490.0     10.421     10.421     11.151     -0.000     -0.000     -0.000 3/21600
  500.0     10.216     10.216     10.932     -0.000     -0.000     -0.000 3/21600
  510.0     10.018     10.018     10.722     -0.000     -0.000     -0.000 3/21600
  520.0      9.829      9.829     10.520     -0.000     -0.000     -0.000 3/21600
  530.0      9.646      9.646     10.325     -0.000     -0.000     -0.000 3/21600
  540.0      9.470      9.470     10.137     -0.000     -0.000     -0.000 3/21600
  550.0      9.301      9.301      9.957     -0.000     -0.000     -0.000 3/21600
  560.0      9.138      9.138      9.782     -0.000     -0.000     -0.000 3/21600
  570.0      8.980      8.980      9.614     -0.000     -0.000     -0.000 3/21600
  580.0      8.828      8.828      9.452     -0.000     -0.000     -0.000 3/21600
  590.0      8.681      8.681      9.295     -0.000     -0.000     -0.000 3/21600
  600.0      8.539      8.539      9.143     -0.000     -0.000     -0.000 3/21600
  610.0      8.401      8.401      8.996     -0.000     -0.000     -0.000 3/21600
  620.0      8.268      8.268      8.854     -0.000     -0.000     -0.000 3/21600
  630.0      8.139      8.139      8.716     -0.000     -0.000     -0.000 3/21600
  640.0      8.014      8.014      8.583     -0.000     -0.000     -0.000 3/21600
  650.0      7.893      7.893      8.453     -0.000     -0.000     -0.000 3/21600
  660.0      7.775      7.775      8.328     -0.000     -0.000     -0.000 3/21600
  670.0      7.661      7.661      8.206     -0.000     -0.000     -0.000 3/21600
  680.0      7.551      7.551      8.088     -0.000     -0.000     -0.000 3/21600
  690.0      7.443      7.443      7.973     -0.000     -0.000     -0.000 3/21600
  700.0      7.339      7.339      7.862     -0.000     -0.000     -0.000 3/21600
  710.0      7.237      7.237      7.753     -0.000     -0.000     -0.000 3/21600
  720.0      7.138      7.138      7.648     -0.000     -0.000     -0.000 3/21600
  730.0      7.042      7.042      7.545     -0.000     -0.000     -0.000 3/21600
  740.0      6.949      6.949      7.445     -0.000     -0.000     -0.000 3/21600
  750.0      6.858      6.858      7.348     -0.000     -0.000     -0.000 3/21600
  760.0      6.769      6.769      7.253     -0.000     -0.000     -0.000 3/21600
  770.0      6.683      6.683      7.161     -0.000     -0.000     -0.000 3/21600
  780.0      6.598      6.598      7.071     -0.000     -0.000     -0.000 3/21600
  790.0      6.516      6.516      6.983     -0.000     -0.000     -0.000 3/21600
  800.0      6.436      6.436      6.897     -0.000     -0.000     -0.000 3/21600
  810.0      6.358      6.358      6.814     -0.000     -0.000     -0.000 3/21600
  820.0      6.282      6.282      6.732     -0.000     -0.000     -0.000 3/21600
  830.0      6.207      6.207      6.653     -0.000     -0.000     -0.000 3/21600
  840.0      6.135      6.135      6.575     -0.000     -0.000     -0.000 3/21600
  850.0      6.064      6.064      6.499     -0.000     -0.000     -0.000 3/21600
  860.0      5.995      5.995      6.425     -0.000     -0.000     -0.000 3/21600
  870.0      5.927      5.927      6.353     -0.000     -0.000     -0.000 3/21600
  880.0      5.861      5.861      6.282     -0.000     -0.000     -0.000 3/21600
  890.0      5.796      5.796      6.213     -0.000     -0.000     -0.000 3/21600
  900.0      5.732      5.732      6.145     -0.000     -0.000     -0.000 3/21600
  910.0      5.670      5.670      6.079     -0.000     -0.000     -0.000 3/21600
  920.0      5.610      5.610      6.014     -0.000     -0.000     -0.000 3/21600
  930.0      5.550      5.550      5.950     -0.000     -0.000     -0.000 3/21600
  940.0      5.492      5.492      5.888     -0.000     -0.000     -0.000 3/21600
  950.0      5.436      5.436      5.828     -0.000     -0.000     -0.000 3/21600
  960.0      5.380      5.380      5.768     -0.000     -0.000     -0.000 3/21600
  970.0      5.325      5.325      5.710     -0.000     -0.000     -0.000 3/21600
  980.0      5.272      5.272      5.652     -0.000     -0.000     -0.000 3/21600
  990.0      5.219      5.219      5.596     -0.000     -0.000     -0.000 3/21600
 1000.0      5.168      5.168      5.541     -0.000     -0.000     -0.000 3/21600

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

