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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [4 4 2]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P6_3mc (186)
Number of symmetry operations in supercell: 384
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.281415999926754    0.000000000000000    0.000000000000000
  b   -2.140707999963377    3.707815020105723    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.034632940000000
Atomic positions (fractional):
   *1 Zn  0.33333333333333  0.66666666666667  0.00064681485892  65.380
    2 Zn  0.66666666666667  0.33333333333333  0.50064681485892  65.380
   *3 Te  0.33333333333333  0.66666666666667  0.37435318514107 127.600
    4 Te  0.66666666666667  0.33333333333333  0.87435318514107 127.600
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.281415999926754    0.000000000000000    0.000000000000000
  b   -2.140707999963377    3.707815020105723    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.034632940000000
Atomic positions (fractional):
   *1 Zn  0.33333333333333  0.66666666666667  0.00064681485892  65.380 > 1
    2 Zn  0.66666666666667  0.33333333333333  0.50064681485892  65.380 > 2
   *3 Te  0.33333333333333  0.66666666666667  0.37435318514107 127.600 > 3
    4 Te  0.66666666666667  0.33333333333333  0.87435318514107 127.600 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   17.125663999707015    0.000000000000000    0.000000000000000
  b   -8.562831999853508   14.831260080422892    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.069265880000000
Atomic positions (fractional):
   *1 Zn  0.08333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    2 Zn  0.33333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    3 Zn  0.58333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    4 Zn  0.83333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    5 Zn  0.08333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    6 Zn  0.33333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    7 Zn  0.58333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    8 Zn  0.83333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    9 Zn  0.08333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   10 Zn  0.33333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   11 Zn  0.58333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   12 Zn  0.83333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   13 Zn  0.08333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   14 Zn  0.33333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   15 Zn  0.58333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   16 Zn  0.83333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   17 Zn  0.08333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   18 Zn  0.33333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   19 Zn  0.58333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   20 Zn  0.83333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   21 Zn  0.08333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   22 Zn  0.33333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   23 Zn  0.58333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   24 Zn  0.83333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   25 Zn  0.08333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   26 Zn  0.33333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   27 Zn  0.58333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   28 Zn  0.83333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   29 Zn  0.08333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   30 Zn  0.33333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   31 Zn  0.58333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   32 Zn  0.83333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   33 Zn  0.16666666666667  0.08333333333333  0.25032340742946  65.380 > 2
   34 Zn  0.41666666666667  0.08333333333333  0.25032340742946  65.380 > 2
   35 Zn  0.66666666666667  0.08333333333333  0.25032340742946  65.380 > 2
   36 Zn  0.91666666666667  0.08333333333333  0.25032340742946  65.380 > 2
   37 Zn  0.16666666666667  0.33333333333333  0.25032340742946  65.380 > 2
   38 Zn  0.41666666666667  0.33333333333333  0.25032340742946  65.380 > 2
   39 Zn  0.66666666666667  0.33333333333333  0.25032340742946  65.380 > 2
   40 Zn  0.91666666666667  0.33333333333333  0.25032340742946  65.380 > 2
   41 Zn  0.16666666666667  0.58333333333333  0.25032340742946  65.380 > 2
   42 Zn  0.41666666666667  0.58333333333333  0.25032340742946  65.380 > 2
   43 Zn  0.66666666666667  0.58333333333333  0.25032340742946  65.380 > 2
   44 Zn  0.91666666666667  0.58333333333333  0.25032340742946  65.380 > 2
   45 Zn  0.16666666666667  0.83333333333333  0.25032340742946  65.380 > 2
   46 Zn  0.41666666666667  0.83333333333333  0.25032340742946  65.380 > 2
   47 Zn  0.66666666666667  0.83333333333333  0.25032340742946  65.380 > 2
   48 Zn  0.91666666666667  0.83333333333333  0.25032340742946  65.380 > 2
   49 Zn  0.16666666666667  0.08333333333333  0.75032340742946  65.380 > 2
   50 Zn  0.41666666666667  0.08333333333333  0.75032340742946  65.380 > 2
   51 Zn  0.66666666666667  0.08333333333333  0.75032340742946  65.380 > 2
   52 Zn  0.91666666666667  0.08333333333333  0.75032340742946  65.380 > 2
   53 Zn  0.16666666666667  0.33333333333333  0.75032340742946  65.380 > 2
   54 Zn  0.41666666666667  0.33333333333333  0.75032340742946  65.380 > 2
   55 Zn  0.66666666666667  0.33333333333333  0.75032340742946  65.380 > 2
   56 Zn  0.91666666666667  0.33333333333333  0.75032340742946  65.380 > 2
   57 Zn  0.16666666666667  0.58333333333333  0.75032340742946  65.380 > 2
   58 Zn  0.41666666666667  0.58333333333333  0.75032340742946  65.380 > 2
   59 Zn  0.66666666666667  0.58333333333333  0.75032340742946  65.380 > 2
   60 Zn  0.91666666666667  0.58333333333333  0.75032340742946  65.380 > 2
   61 Zn  0.16666666666667  0.83333333333333  0.75032340742946  65.380 > 2
   62 Zn  0.41666666666667  0.83333333333333  0.75032340742946  65.380 > 2
   63 Zn  0.66666666666667  0.83333333333333  0.75032340742946  65.380 > 2
   64 Zn  0.91666666666667  0.83333333333333  0.75032340742946  65.380 > 2
  *65 Te  0.08333333333333  0.16666666666667  0.18717659257054 127.600 > 3
   66 Te  0.33333333333333  0.16666666666667  0.18717659257054 127.600 > 3
   67 Te  0.58333333333333  0.16666666666667  0.18717659257054 127.600 > 3
   68 Te  0.83333333333333  0.16666666666667  0.18717659257054 127.600 > 3
   69 Te  0.08333333333333  0.41666666666667  0.18717659257054 127.600 > 3
   70 Te  0.33333333333333  0.41666666666667  0.18717659257054 127.600 > 3
   71 Te  0.58333333333333  0.41666666666667  0.18717659257054 127.600 > 3
   72 Te  0.83333333333333  0.41666666666667  0.18717659257054 127.600 > 3
   73 Te  0.08333333333333  0.66666666666667  0.18717659257054 127.600 > 3
   74 Te  0.33333333333333  0.66666666666667  0.18717659257054 127.600 > 3
   75 Te  0.58333333333333  0.66666666666667  0.18717659257054 127.600 > 3
   76 Te  0.83333333333333  0.66666666666667  0.18717659257054 127.600 > 3
   77 Te  0.08333333333333  0.91666666666667  0.18717659257054 127.600 > 3
   78 Te  0.33333333333333  0.91666666666667  0.18717659257054 127.600 > 3
   79 Te  0.58333333333333  0.91666666666667  0.18717659257054 127.600 > 3
   80 Te  0.83333333333333  0.91666666666667  0.18717659257054 127.600 > 3
   81 Te  0.08333333333333  0.16666666666667  0.68717659257053 127.600 > 3
   82 Te  0.33333333333333  0.16666666666667  0.68717659257053 127.600 > 3
   83 Te  0.58333333333333  0.16666666666667  0.68717659257053 127.600 > 3
   84 Te  0.83333333333333  0.16666666666667  0.68717659257053 127.600 > 3
   85 Te  0.08333333333333  0.41666666666667  0.68717659257053 127.600 > 3
   86 Te  0.33333333333333  0.41666666666667  0.68717659257053 127.600 > 3
   87 Te  0.58333333333333  0.41666666666667  0.68717659257053 127.600 > 3
   88 Te  0.83333333333333  0.41666666666667  0.68717659257053 127.600 > 3
   89 Te  0.08333333333333  0.66666666666667  0.68717659257053 127.600 > 3
   90 Te  0.33333333333333  0.66666666666667  0.68717659257053 127.600 > 3
   91 Te  0.58333333333333  0.66666666666667  0.68717659257053 127.600 > 3
   92 Te  0.83333333333333  0.66666666666667  0.68717659257053 127.600 > 3
   93 Te  0.08333333333333  0.91666666666667  0.68717659257053 127.600 > 3
   94 Te  0.33333333333333  0.91666666666667  0.68717659257053 127.600 > 3
   95 Te  0.58333333333333  0.91666666666667  0.68717659257053 127.600 > 3
   96 Te  0.83333333333333  0.91666666666667  0.68717659257053 127.600 > 3
   97 Te  0.16666666666667  0.08333333333333  0.43717659257053 127.600 > 4
   98 Te  0.41666666666667  0.08333333333333  0.43717659257053 127.600 > 4
   99 Te  0.66666666666667  0.08333333333333  0.43717659257053 127.600 > 4
  100 Te  0.91666666666667  0.08333333333333  0.43717659257053 127.600 > 4
  101 Te  0.16666666666667  0.33333333333333  0.43717659257053 127.600 > 4
  102 Te  0.41666666666667  0.33333333333333  0.43717659257053 127.600 > 4
  103 Te  0.66666666666667  0.33333333333333  0.43717659257053 127.600 > 4
  104 Te  0.91666666666667  0.33333333333333  0.43717659257053 127.600 > 4
  105 Te  0.16666666666667  0.58333333333333  0.43717659257053 127.600 > 4
  106 Te  0.41666666666667  0.58333333333333  0.43717659257053 127.600 > 4
  107 Te  0.66666666666667  0.58333333333333  0.43717659257053 127.600 > 4
  108 Te  0.91666666666667  0.58333333333333  0.43717659257053 127.600 > 4
  109 Te  0.16666666666667  0.83333333333333  0.43717659257053 127.600 > 4
  110 Te  0.41666666666667  0.83333333333333  0.43717659257053 127.600 > 4
  111 Te  0.66666666666667  0.83333333333333  0.43717659257053 127.600 > 4
  112 Te  0.91666666666667  0.83333333333333  0.43717659257053 127.600 > 4
  113 Te  0.16666666666667  0.08333333333333  0.93717659257053 127.600 > 4
  114 Te  0.41666666666667  0.08333333333333  0.93717659257053 127.600 > 4
  115 Te  0.66666666666667  0.08333333333333  0.93717659257053 127.600 > 4
  116 Te  0.91666666666667  0.08333333333333  0.93717659257053 127.600 > 4
  117 Te  0.16666666666667  0.33333333333333  0.93717659257053 127.600 > 4
  118 Te  0.41666666666667  0.33333333333333  0.93717659257053 127.600 > 4
  119 Te  0.66666666666667  0.33333333333333  0.93717659257053 127.600 > 4
  120 Te  0.91666666666667  0.33333333333333  0.93717659257053 127.600 > 4
  121 Te  0.16666666666667  0.58333333333333  0.93717659257053 127.600 > 4
  122 Te  0.41666666666667  0.58333333333333  0.93717659257053 127.600 > 4
  123 Te  0.66666666666667  0.58333333333333  0.93717659257053 127.600 > 4
  124 Te  0.91666666666667  0.58333333333333  0.93717659257053 127.600 > 4
  125 Te  0.16666666666667  0.83333333333333  0.93717659257053 127.600 > 4
  126 Te  0.41666666666667  0.83333333333333  0.93717659257053 127.600 > 4
  127 Te  0.66666666666667  0.83333333333333  0.93717659257053 127.600 > 4
  128 Te  0.91666666666667  0.83333333333333  0.93717659257053 127.600 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           11.9415109    0.0000000    0.0000000
           -0.0000000   11.9415109    0.0000000
            0.0000000    0.0000000   10.7467801
-------------------------- Born effective charges --------------------------
    1 Zn    2.1681442    0.0000000    0.0000000
           -0.0000000    2.1681442    0.0000000
            0.0000000    0.0000000    2.1213866
    2 Zn    2.1681442    0.0000000    0.0000000
           -0.0000000    2.1681442    0.0000000
            0.0000000    0.0000000    2.1213866
    3 Te   -2.1681442   -0.0000000    0.0000000
            0.0000000   -2.1681442    0.0000000
            0.0000000    0.0000000   -2.1213866
    4 Te   -2.1681442   -0.0000000    0.0000000
            0.0000000   -2.1681442    0.0000000
            0.0000000    0.0000000   -2.1213866
----------------------------------------------------------------------------
Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
-------------------------------- Symfc start -------------------------------
Symfc version 1.5.4 (https://github.com/symfc/symfc)
Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024)
Computing [2] order force constants.
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 172
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
-----
Solver_atoms: 1 -- 128 / 128
Time (Solver_compr_matrix_reshape): 0.010
Solver_block: 80 / 80
 - Time: 0.156
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.167
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 172
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-07 23:04:38]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.281415999926754    0.000000000000000    0.000000000000000
  b   -2.140707999963377    3.707815020105723    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.034632940000000
Atomic positions (fractional):
    1 Zn  0.33333333333333  0.66666666666667  0.00064681485892  65.380
    2 Zn  0.66666666666667  0.33333333333333  0.50064681485892  65.380
    3 Te  0.33333333333333  0.66666666666667  0.37435318514107 127.600
    4 Te  0.66666666666667  0.33333333333333  0.87435318514107 127.600
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   17.125663999707015    0.000000000000000    0.000000000000000
  b   -8.562831999853508   14.831260080422892    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.069265880000000
Atomic positions (fractional):
    1 Zn  0.08333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    2 Zn  0.33333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    3 Zn  0.58333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    4 Zn  0.83333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    5 Zn  0.08333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    6 Zn  0.33333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    7 Zn  0.58333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    8 Zn  0.83333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    9 Zn  0.08333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   10 Zn  0.33333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   11 Zn  0.58333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   12 Zn  0.83333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   13 Zn  0.08333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   14 Zn  0.33333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   15 Zn  0.58333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   16 Zn  0.83333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   17 Zn  0.08333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   18 Zn  0.33333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   19 Zn  0.58333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   20 Zn  0.83333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   21 Zn  0.08333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   22 Zn  0.33333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   23 Zn  0.58333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   24 Zn  0.83333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   25 Zn  0.08333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   26 Zn  0.33333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   27 Zn  0.58333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   28 Zn  0.83333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   29 Zn  0.08333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   30 Zn  0.33333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   31 Zn  0.58333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   32 Zn  0.83333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   33 Zn  0.16666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   34 Zn  0.41666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   35 Zn  0.66666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   36 Zn  0.91666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   37 Zn  0.16666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   38 Zn  0.41666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   39 Zn  0.66666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   40 Zn  0.91666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   41 Zn  0.16666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   42 Zn  0.41666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   43 Zn  0.66666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   44 Zn  0.91666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   45 Zn  0.16666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   46 Zn  0.41666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   47 Zn  0.66666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   48 Zn  0.91666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   49 Zn  0.16666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   50 Zn  0.41666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   51 Zn  0.66666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   52 Zn  0.91666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   53 Zn  0.16666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   54 Zn  0.41666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   55 Zn  0.66666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   56 Zn  0.91666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   57 Zn  0.16666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   58 Zn  0.41666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   59 Zn  0.66666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   60 Zn  0.91666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   61 Zn  0.16666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   62 Zn  0.41666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   63 Zn  0.66666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   64 Zn  0.91666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   65 Te  0.08333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   66 Te  0.33333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   67 Te  0.58333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   68 Te  0.83333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   69 Te  0.08333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   70 Te  0.33333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   71 Te  0.58333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   72 Te  0.83333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   73 Te  0.08333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   74 Te  0.33333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   75 Te  0.58333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   76 Te  0.83333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   77 Te  0.08333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   78 Te  0.33333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   79 Te  0.58333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   80 Te  0.83333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   81 Te  0.08333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   82 Te  0.33333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   83 Te  0.58333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   84 Te  0.83333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   85 Te  0.08333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   86 Te  0.33333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   87 Te  0.58333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   88 Te  0.83333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   89 Te  0.08333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   90 Te  0.33333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   91 Te  0.58333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   92 Te  0.83333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   93 Te  0.08333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   94 Te  0.33333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   95 Te  0.58333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   96 Te  0.83333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   97 Te  0.16666666666667  0.08333333333333  0.43717659257053 127.600 > 97
   98 Te  0.41666666666667  0.08333333333333  0.43717659257053 127.600 > 97
   99 Te  0.66666666666667  0.08333333333333  0.43717659257053 127.600 > 97
  100 Te  0.91666666666667  0.08333333333333  0.43717659257053 127.600 > 97
  101 Te  0.16666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  102 Te  0.41666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  103 Te  0.66666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  104 Te  0.91666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  105 Te  0.16666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  106 Te  0.41666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  107 Te  0.66666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  108 Te  0.91666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  109 Te  0.16666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  110 Te  0.41666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  111 Te  0.66666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  112 Te  0.91666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  113 Te  0.16666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  114 Te  0.41666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  115 Te  0.66666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  116 Te  0.91666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  117 Te  0.16666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  118 Te  0.41666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  119 Te  0.66666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  120 Te  0.91666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  121 Te  0.16666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  122 Te  0.41666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  123 Te  0.66666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  124 Te  0.91666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  125 Te  0.16666666666667  0.83333333333333  0.93717659257053 127.600 > 97
  126 Te  0.41666666666667  0.83333333333333  0.93717659257053 127.600 > 97
  127 Te  0.66666666666667  0.83333333333333  0.93717659257053 127.600 > 97
  128 Te  0.91666666666667  0.83333333333333  0.93717659257053 127.600 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           11.9415109    0.0000000    0.0000000
           -0.0000000   11.9415109    0.0000000
            0.0000000    0.0000000   10.7467801
-------------------------- Born effective charges --------------------------
    1 Zn    2.1681442    0.0000000    0.0000000
           -0.0000000    2.1681442    0.0000000
            0.0000000    0.0000000    2.1213866
    2 Zn    2.1681442    0.0000000    0.0000000
           -0.0000000    2.1681442    0.0000000
            0.0000000    0.0000000    2.1213866
    3 Te   -2.1681442   -0.0000000    0.0000000
            0.0000000   -2.1681442    0.0000000
            0.0000000    0.0000000   -2.1213866
    4 Te   -2.1681442   -0.0000000    0.0000000
            0.0000000   -2.1681442    0.0000000
            0.0000000    0.0000000   -2.1213866
----------------------------------------------------------------------------
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.00000095 (xyx) 0.00000095 (xyx) 0.00000095 (xxy)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
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 23:04:43]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 23:04:44]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.281415999926754    0.000000000000000    0.000000000000000
  b   -2.140707999963377    3.707815020105723    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.034632940000000
Atomic positions (fractional):
    1 Zn  0.33333333333333  0.66666666666667  0.00064681485892  65.380
    2 Zn  0.66666666666667  0.33333333333333  0.50064681485892  65.380
    3 Te  0.33333333333333  0.66666666666667  0.37435318514107 127.600
    4 Te  0.66666666666667  0.33333333333333  0.87435318514107 127.600
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   17.125663999707015    0.000000000000000    0.000000000000000
  b   -8.562831999853508   14.831260080422892    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.069265880000000
Atomic positions (fractional):
    1 Zn  0.08333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    2 Zn  0.33333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    3 Zn  0.58333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    4 Zn  0.83333333333333  0.16666666666667  0.00032340742946  65.380 > 1
    5 Zn  0.08333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    6 Zn  0.33333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    7 Zn  0.58333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    8 Zn  0.83333333333333  0.41666666666667  0.00032340742946  65.380 > 1
    9 Zn  0.08333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   10 Zn  0.33333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   11 Zn  0.58333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   12 Zn  0.83333333333333  0.66666666666667  0.00032340742946  65.380 > 1
   13 Zn  0.08333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   14 Zn  0.33333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   15 Zn  0.58333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   16 Zn  0.83333333333333  0.91666666666667  0.00032340742946  65.380 > 1
   17 Zn  0.08333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   18 Zn  0.33333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   19 Zn  0.58333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   20 Zn  0.83333333333333  0.16666666666667  0.50032340742946  65.380 > 1
   21 Zn  0.08333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   22 Zn  0.33333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   23 Zn  0.58333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   24 Zn  0.83333333333333  0.41666666666667  0.50032340742946  65.380 > 1
   25 Zn  0.08333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   26 Zn  0.33333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   27 Zn  0.58333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   28 Zn  0.83333333333333  0.66666666666667  0.50032340742946  65.380 > 1
   29 Zn  0.08333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   30 Zn  0.33333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   31 Zn  0.58333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   32 Zn  0.83333333333333  0.91666666666667  0.50032340742946  65.380 > 1
   33 Zn  0.16666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   34 Zn  0.41666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   35 Zn  0.66666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   36 Zn  0.91666666666667  0.08333333333333  0.25032340742946  65.380 > 33
   37 Zn  0.16666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   38 Zn  0.41666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   39 Zn  0.66666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   40 Zn  0.91666666666667  0.33333333333333  0.25032340742946  65.380 > 33
   41 Zn  0.16666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   42 Zn  0.41666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   43 Zn  0.66666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   44 Zn  0.91666666666667  0.58333333333333  0.25032340742946  65.380 > 33
   45 Zn  0.16666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   46 Zn  0.41666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   47 Zn  0.66666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   48 Zn  0.91666666666667  0.83333333333333  0.25032340742946  65.380 > 33
   49 Zn  0.16666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   50 Zn  0.41666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   51 Zn  0.66666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   52 Zn  0.91666666666667  0.08333333333333  0.75032340742946  65.380 > 33
   53 Zn  0.16666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   54 Zn  0.41666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   55 Zn  0.66666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   56 Zn  0.91666666666667  0.33333333333333  0.75032340742946  65.380 > 33
   57 Zn  0.16666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   58 Zn  0.41666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   59 Zn  0.66666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   60 Zn  0.91666666666667  0.58333333333333  0.75032340742946  65.380 > 33
   61 Zn  0.16666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   62 Zn  0.41666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   63 Zn  0.66666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   64 Zn  0.91666666666667  0.83333333333333  0.75032340742946  65.380 > 33
   65 Te  0.08333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   66 Te  0.33333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   67 Te  0.58333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   68 Te  0.83333333333333  0.16666666666667  0.18717659257054 127.600 > 65
   69 Te  0.08333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   70 Te  0.33333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   71 Te  0.58333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   72 Te  0.83333333333333  0.41666666666667  0.18717659257054 127.600 > 65
   73 Te  0.08333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   74 Te  0.33333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   75 Te  0.58333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   76 Te  0.83333333333333  0.66666666666667  0.18717659257054 127.600 > 65
   77 Te  0.08333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   78 Te  0.33333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   79 Te  0.58333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   80 Te  0.83333333333333  0.91666666666667  0.18717659257054 127.600 > 65
   81 Te  0.08333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   82 Te  0.33333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   83 Te  0.58333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   84 Te  0.83333333333333  0.16666666666667  0.68717659257053 127.600 > 65
   85 Te  0.08333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   86 Te  0.33333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   87 Te  0.58333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   88 Te  0.83333333333333  0.41666666666667  0.68717659257053 127.600 > 65
   89 Te  0.08333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   90 Te  0.33333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   91 Te  0.58333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   92 Te  0.83333333333333  0.66666666666667  0.68717659257053 127.600 > 65
   93 Te  0.08333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   94 Te  0.33333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   95 Te  0.58333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   96 Te  0.83333333333333  0.91666666666667  0.68717659257053 127.600 > 65
   97 Te  0.16666666666667  0.08333333333333  0.43717659257053 127.600 > 97
   98 Te  0.41666666666667  0.08333333333333  0.43717659257053 127.600 > 97
   99 Te  0.66666666666667  0.08333333333333  0.43717659257053 127.600 > 97
  100 Te  0.91666666666667  0.08333333333333  0.43717659257053 127.600 > 97
  101 Te  0.16666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  102 Te  0.41666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  103 Te  0.66666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  104 Te  0.91666666666667  0.33333333333333  0.43717659257053 127.600 > 97
  105 Te  0.16666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  106 Te  0.41666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  107 Te  0.66666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  108 Te  0.91666666666667  0.58333333333333  0.43717659257053 127.600 > 97
  109 Te  0.16666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  110 Te  0.41666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  111 Te  0.66666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  112 Te  0.91666666666667  0.83333333333333  0.43717659257053 127.600 > 97
  113 Te  0.16666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  114 Te  0.41666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  115 Te  0.66666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  116 Te  0.91666666666667  0.08333333333333  0.93717659257053 127.600 > 97
  117 Te  0.16666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  118 Te  0.41666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  119 Te  0.66666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  120 Te  0.91666666666667  0.33333333333333  0.93717659257053 127.600 > 97
  121 Te  0.16666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  122 Te  0.41666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  123 Te  0.66666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  124 Te  0.91666666666667  0.58333333333333  0.93717659257053 127.600 > 97
  125 Te  0.16666666666667  0.83333333333333  0.93717659257053 127.600 > 97
  126 Te  0.41666666666667  0.83333333333333  0.93717659257053 127.600 > 97
  127 Te  0.66666666666667  0.83333333333333  0.93717659257053 127.600 > 97
  128 Te  0.91666666666667  0.83333333333333  0.93717659257053 127.600 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           11.9415109    0.0000000    0.0000000
           -0.0000000   11.9415109    0.0000000
            0.0000000    0.0000000   10.7467801
-------------------------- Born effective charges --------------------------
    1 Zn    2.1681442    0.0000000    0.0000000
           -0.0000000    2.1681442    0.0000000
            0.0000000    0.0000000    2.1213866
    2 Zn    2.1681442    0.0000000    0.0000000
           -0.0000000    2.1681442    0.0000000
            0.0000000    0.0000000    2.1213866
    3 Te   -2.1681442   -0.0000000    0.0000000
            0.0000000   -2.1681442    0.0000000
            0.0000000    0.0000000   -2.1213866
    4 Te   -2.1681442   -0.0000000    0.0000000
            0.0000000   -2.1681442    0.0000000
            0.0000000    0.0000000   -2.1213866
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000095 (xyx) 0.00000095 (xyx) 0.00000095 (xxy)
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 13 13 7 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.86, Number of G-points: 313, Lambda: 0.31
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/84) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 84
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.138   (   0.000    0.000    0.000)    0.000
   1.138   (  -0.000   -0.000   -0.000)    0.000
   4.041   (  -0.000    0.000    0.000)    0.000
   5.301   (   0.000    0.000    0.000)    0.000
   5.301   (   0.000   -0.000   -0.000)    0.000
   5.432   (   0.000    0.000   -0.000)    0.000
   5.482   (   0.000    0.000    0.000)    0.000
   5.510   (   0.000    0.000   -0.000)    0.000
   5.510   (   0.000   -0.000   -0.000)    0.000
======================= Grid point 1 (2/84) =======================
q-point: ( 0.08  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.361   (  14.668    8.469    0.000)   16.937
   0.366   (  15.376    8.878    0.000)   17.755
   0.800   (  32.459   18.740    0.000)   37.480
   1.191   (   4.166    2.405    0.000)    4.810
   1.383   (  17.695   10.216    0.000)   20.432
   3.978   (  -4.789   -2.765    0.000)    5.530
   5.287   (  -0.996   -0.575    0.000)    1.150
   5.344   (   2.621    1.513    0.000)    3.027
   5.394   (  -2.439   -1.408    0.000)    2.816
   5.468   (  -1.000   -0.577    0.000)    1.155
   5.505   (  -0.578   -0.334    0.000)    0.667
   6.118   (  -1.026   -0.592    0.000)    1.185
======================= Grid point 2 (3/84) =======================
q-point: ( 0.15  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.691   (  12.631    7.292    0.000)   14.585
   0.726   (  14.654    8.461    0.000)   16.921
   1.320   (   6.140    3.545    0.000)    7.090
   1.539   (  28.703   16.572    0.000)   33.143
   1.846   (  18.776   10.840    0.000)   21.681
   3.846   (  -5.275   -3.046    0.000)    6.091
   5.262   (  -1.006   -0.581    0.000)    1.161
   5.280   (  -4.976   -2.873    0.000)    5.746
   5.414   (  -0.630   -0.364    0.000)    0.728
   5.471   (   2.897    1.673    0.000)    3.345
   5.479   (  -1.525   -0.881    0.000)    1.761
   6.081   (  -2.050   -1.184    0.000)    2.367
======================= Grid point 3 (4/84) =======================
q-point: ( 0.23  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.952   (   8.899    5.138    0.000)   10.276
   1.065   (  13.614    7.860    0.000)   15.720
   1.461   (   5.225    3.017    0.000)    6.034
   2.093   (  13.291    7.674    0.000)   15.347
   2.273   (  20.163   11.641    0.000)   23.282
   3.780   (   0.638    0.368    0.000)    0.736
   5.159   (  -4.886   -2.821    0.000)    5.642
   5.242   (  -0.623   -0.360    0.000)    0.719
   5.381   (  -1.400   -0.808    0.000)    1.617
   5.438   (  -1.658   -0.957    0.000)    1.915
   5.538   (   1.785    1.031    0.000)    2.061
   6.015   (  -3.618   -2.089    0.000)    4.178
======================= Grid point 4 (5/84) =======================
q-point: ( 0.31  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.111   (   4.423    2.554    0.000)    5.107
   1.374   (  12.036    6.949    0.000)   13.898
   1.553   (   2.208    1.275    0.000)    2.550
   2.216   (   0.137    0.079    0.000)    0.158
   2.776   (  19.805   11.435    0.000)   22.869
   3.876   (   6.402    3.696    0.000)    7.393
   5.057   (  -3.415   -1.971    0.000)    3.943
   5.232   (  -0.218   -0.126    0.000)    0.252
   5.358   (  -0.385   -0.222    0.000)    0.444
   5.409   (  -0.669   -0.386    0.000)    0.772
   5.549   (  -0.674   -0.389    0.000)    0.778
   5.900   (  -6.112   -3.529    0.000)    7.058
======================= Grid point 5 (6/84) =======================
q-point: ( 0.38  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.174   (   1.149    0.663    0.000)    1.327
   1.568   (  -0.682   -0.394    0.000)    0.787
   1.634   (   9.396    5.425    0.000)   10.849
   2.179   (  -2.371   -1.369    0.000)    2.738
   3.196   (  14.993    8.656    0.000)   17.312
   4.020   (   4.317    2.493    0.000)    4.985
   5.013   (   0.188    0.109    0.000)    0.217
   5.232   (   0.179    0.103    0.000)    0.206
   5.363   (   0.722    0.417    0.000)    0.834
   5.406   (   0.244    0.141    0.000)    0.281
   5.518   (  -1.606   -0.927    0.000)    1.854
   5.723   (  -8.423   -4.863    0.000)    9.726
======================= Grid point 6 (7/84) =======================
q-point: ( 0.46  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.185   (   0.067    0.039    0.000)    0.077
   1.543   (  -0.892   -0.515    0.000)    1.030
   1.805   (   4.280    2.471    0.000)    4.942
   2.126   (  -1.654   -0.955    0.000)    1.910
   3.470   (   6.841    3.949    0.000)    7.899
   4.059   (  -0.351   -0.203    0.000)    0.405
   5.063   (   2.838    1.638    0.000)    3.277
   5.237   (   0.164    0.095    0.000)    0.189
   5.391   (   1.691    0.976    0.000)    1.953
   5.413   (   0.249    0.144    0.000)    0.288
   5.486   (  -0.907   -0.524    0.000)    1.048
   5.526   (  -6.759   -3.902    0.000)    7.804
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.631   (   8.543   14.797    0.000)   17.086
   0.661   (   9.625   16.670    0.000)   19.249
   1.256   (   9.342   16.181    0.000)   18.685
   1.399   (  12.121   20.994    0.000)   24.242
   1.680   (  10.125   17.537    0.000)   20.250
   3.881   (  -3.362   -5.823    0.000)    6.724
   5.278   (  -0.470   -0.814    0.000)    0.939
   5.310   (  -2.756   -4.773    0.000)    5.512
   5.401   (   1.613    2.794    0.000)    3.226
   5.446   (  -0.820   -1.421    0.000)    1.641
   5.514   (   0.360    0.624    0.000)    0.721
   6.084   (  -1.481   -2.565    0.000)    2.961
======================= Grid point 15 (9/84) =======================
q-point: ( 0.15  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.945   (  10.564   11.948    0.000)   15.949
   0.979   (   5.168   17.723    0.000)   18.461
   1.477   (   3.570   10.700    0.000)   11.280
   1.928   (  16.069   18.880    0.000)   24.792
   2.048   (  15.870    8.843    0.000)   18.167
   3.780   (  -1.684   -2.501    0.000)    3.015
   5.197   (  -3.833   -4.864    0.000)    6.193
   5.260   (  -0.917    0.192    0.000)    0.937
   5.405   (  -1.905   -0.367    0.000)    1.940
   5.459   (   1.183    0.394    0.000)    1.247
   5.522   (  -0.198    1.303    0.000)    1.318
   6.016   (  -1.653   -4.890    0.000)    5.162
======================= Grid point 16 (10/84) =======================
q-point: ( 0.23  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.199   (  -0.340   16.470    0.000)   16.474
   1.252   (  10.594    9.916    0.000)   14.511
   1.619   (  -0.676   11.126    0.000)   11.146
   2.183   (   4.247    1.546    0.000)    4.519
   2.524   (  20.830   12.169    0.000)   24.124
   3.801   (   4.388    1.611    0.000)    4.675
   5.082   (  -3.321   -4.447    0.000)    5.550
   5.252   (  -0.769    1.130    0.000)    1.367
   5.369   (  -1.522   -0.449    0.000)    1.587
   5.453   (  -2.458    2.091    0.000)    3.227
   5.535   (   0.987   -1.353    0.000)    1.675
   5.920   (  -2.540   -6.790    0.000)    7.249
======================= Grid point 17 (11/84) =======================
q-point: ( 0.31  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.308   (  -5.104   14.689    0.000)   15.551
   1.527   (   9.597    7.757    0.000)   12.340
   1.679   (  -4.667    9.997    0.000)   11.033
   2.201   (  -1.087   -1.066    0.000)    1.523
   2.985   (  17.987    8.767    0.000)   20.010
   3.932   (   6.806    1.321    0.000)    6.933
   4.998   (  -0.846   -3.524    0.000)    3.624
   5.254   (  -0.794    1.947    0.000)    2.103
   5.347   (  -0.088   -0.692    0.000)    0.697
   5.434   (  -1.776    1.936    0.000)    2.627
   5.523   (  -0.359   -1.732    0.000)    1.769
   5.778   (  -4.908   -7.740    0.000)    9.165
======================= Grid point 18 (12/84) =======================
q-point: ( 0.38  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.337   (  -7.371   13.461    0.000)   15.347
   1.670   (  -6.705    9.481    0.000)   11.612
   1.747   (   6.447    5.461    0.000)    8.449
   2.154   (  -2.424   -0.987    0.000)    2.618
   3.339   (  12.594    5.161    0.000)   13.610
   4.018   (   2.884   -2.707    0.000)    3.955
   4.999   (   3.794   -1.380    0.000)    4.037
   5.262   (  -0.931    2.452    0.000)    2.623
   5.350   (   1.438   -1.510    0.000)    2.085
   5.429   (  -0.604    0.992    0.000)    1.162
   5.495   (  -1.003   -0.881    0.000)    1.334
   5.594   (  -6.343   -6.663    0.000)    9.199
======================= Grid point 19 (13/84) =======================
q-point: ( 0.46  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.339   (  -7.633   13.221    0.000)   15.266
   1.655   (  -5.962   10.326    0.000)   11.924
   1.844   (  -0.744    1.288    0.000)    1.487
   2.118   (  -0.163    0.282    0.000)    0.325
   3.495   (   0.885   -1.533    0.000)    1.770
   4.013   (   2.229   -3.861    0.000)    4.458
   5.050   (   1.691   -2.930    0.000)    3.383
   5.268   (  -1.555    2.693    0.000)    3.110
   5.356   (   2.161   -3.743    0.000)    4.323
   5.431   (  -0.067    0.115    0.000)    0.133
   5.468   (   0.261   -0.452    0.000)    0.522
   5.488   (  -0.391    0.677    0.000)    0.782
======================= Grid point 29 (14/84) =======================
q-point: ( 0.15  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.201   (   7.211   12.490    0.000)   14.422
   1.329   (   8.410   14.567    0.000)   16.821
   1.706   (   5.841   10.117    0.000)   11.682
   2.145   (   2.347    4.065    0.000)    4.694
   2.365   (  11.387   19.722    0.000)   22.773
   3.768   (   0.918    1.590    0.000)    1.836
   5.088   (  -3.133   -5.426    0.000)    6.265
   5.266   (   0.216    0.373    0.000)    0.431
   5.390   (  -0.629   -1.090    0.000)    1.259
   5.486   (   0.957    1.658    0.000)    1.914
   5.522   (  -0.695   -1.204    0.000)    1.391
   5.894   (  -3.769   -6.529    0.000)    7.539
======================= Grid point 30 (15/84) =======================
q-point: ( 0.23  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.454   (   7.539   10.023    0.000)   12.542
   1.556   (  -0.360   15.538    0.000)   15.542
   1.864   (  -0.857   11.179    0.000)   11.212
   2.192   (   0.947   -0.038    0.000)    0.948
   2.775   (  15.285   12.298    0.000)   19.618
   3.840   (   4.319    1.987    0.000)    4.754
   4.980   (  -2.370   -5.104    0.000)    5.627
   5.274   (   0.067    0.562    0.000)    0.566
   5.365   (  -1.458    0.167    0.000)    1.468
   5.483   (  -1.586   -1.719    0.000)    2.338
   5.512   (  -0.061    1.560    0.000)    1.562
   5.755   (  -2.999   -8.468    0.000)    8.983
======================= Grid point 31 (16/84) =======================
q-point: ( 0.31  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.625   (  -4.256   13.245    0.000)   13.912
   1.703   (   4.671    9.793    0.000)   10.850
   1.914   (  -6.012   11.268    0.000)   12.771
   2.183   (  -1.469   -0.684    0.000)    1.621
   3.151   (  14.420    7.660    0.000)   16.328
   3.924   (   3.853   -2.092    0.000)    4.384
   4.922   (   1.474   -3.407    0.000)    3.712
   5.278   (   0.510   -0.921    0.000)    1.053
   5.346   (  -1.382    1.168    0.000)    1.809
   5.442   (  -0.894   -1.414    0.000)    1.673
   5.514   (  -1.277    0.904    0.000)    1.564
   5.610   (  -3.679   -7.477    0.000)    8.333
======================= Grid point 32 (17/84) =======================
q-point: ( 0.38  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.646   (  -8.151   14.355    0.000)   16.507
   1.863   (   2.522    5.858    0.000)    6.378
   1.907   (  -7.390   11.762    0.000)   13.891
   2.136   (  -2.036   -0.816    0.000)    2.193
   3.419   (   7.634    2.882    0.000)    8.160
   3.921   (   1.621   -6.132    0.000)    6.343
   4.960   (   4.463   -2.064    0.000)    4.917
   5.266   (   1.562   -4.952    0.000)    5.193
   5.336   (  -1.892    2.471    0.000)    3.113
   5.422   (   0.104   -1.159    0.000)    1.164
   5.472   (  -1.847   -2.185    0.000)    2.861
   5.518   (  -1.626    0.398    0.000)    1.674
======================= Grid point 43 (18/84) =======================
q-point: ( 0.23  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.657   (   5.460    9.457    0.000)   10.920
   1.803   (   4.575    7.923    0.000)    9.149
   2.043   (   3.257    5.642    0.000)    6.515
   2.190   (  -0.197   -0.341    0.000)    0.394
   3.045   (   7.862   13.617    0.000)   15.724
   3.872   (   0.380    0.658    0.000)    0.760
   4.889   (  -1.739   -3.011    0.000)    3.477
   5.264   (  -0.911   -1.577    0.000)    1.822
   5.372   (   0.222    0.385    0.000)    0.445
   5.439   (  -1.035   -1.792    0.000)    2.069
   5.542   (   0.354    0.613    0.000)    0.708
   5.592   (  -3.917   -6.784    0.000)    7.834
======================= Grid point 44 (19/84) =======================
q-point: ( 0.31  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.829   (   2.295    8.168    0.000)    8.484
   1.911   (  -0.060    7.970    0.000)    7.970
   2.104   (  -2.746    6.525    0.000)    7.080
   2.165   (  -1.858   -1.387    0.000)    2.318
   3.320   (   8.767    8.788    0.000)   12.413
   3.848   (  -0.599   -5.140    0.000)    5.175
   4.880   (   2.556   -0.034    0.000)    2.557
   5.206   (  -1.339   -5.581    0.000)    5.739
   5.378   (  -0.431    1.264    0.000)    1.335
   5.412   (  -0.393   -1.139    0.000)    1.205
   5.483   (  -2.133   -2.952    0.000)    3.642
   5.542   (  -0.690    0.116    0.000)    0.699
======================= Grid point 45 (20/84) =======================
q-point: ( 0.38  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.881   (  -4.619    8.001    0.000)    9.239
   1.982   (  -2.759    4.778    0.000)    5.517
   2.106   (  -2.592    4.489    0.000)    5.183
   2.122   (  -0.948    1.643    0.000)    1.897
   3.472   (  -1.449    2.510    0.000)    2.898
   3.789   (   3.370   -5.837    0.000)    6.739
   4.929   (   0.388   -0.673    0.000)    0.777
   5.142   (   3.268   -5.660    0.000)    6.536
   5.383   (  -0.806    1.397    0.000)    1.613
   5.400   (   0.424   -0.735    0.000)    0.849
   5.451   (   0.001   -0.002    0.000)    0.002
   5.534   (  -0.385    0.667    0.000)    0.770
======================= Grid point 58 (21/84) =======================
q-point: ( 0.31  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.982   (   3.401    5.891    0.000)    6.802
   2.013   (   1.583    2.741    0.000)    3.165
   2.110   (  -2.205   -3.820    0.000)    4.411
   2.186   (   0.903    1.564    0.000)    1.806
   3.508   (   4.984    8.632    0.000)    9.967
   3.717   (  -4.074   -7.056    0.000)    8.148
   4.928   (   2.413    4.180    0.000)    4.826
   5.073   (  -3.667   -6.352    0.000)    7.335
   5.390   (   0.114    0.197    0.000)    0.228
   5.397   (  -0.253   -0.438    0.000)    0.506
   5.455   (  -0.142   -0.247    0.000)    0.285
   5.541   (  -0.094   -0.163    0.000)    0.188
======================= Grid point 181 (22/84) =======================
q-point: ( 0.00  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.363   (   0.000    0.000   17.151)   17.151
   0.363   (   0.000   -0.000   17.151)   17.151
   0.833   (   0.000    0.000   40.522)   40.522
   1.141   (   0.000    0.000    0.129)    0.129
   1.141   (  -0.000   -0.000    0.129)    0.129
   3.906   (  -0.000    0.000  -12.808)   12.808
   5.310   (   0.000    0.000    0.881)    0.881
   5.310   (   0.000    0.000    0.881)    0.881
   5.492   (   0.000   -0.000    5.501)    5.501
   5.499   (   0.000    0.000   -1.103)    1.103
   5.499   (   0.000   -0.000   -1.103)    1.103
   6.117   (   0.000   -0.000   -2.408)    2.408
======================= Grid point 182 (23/84) =======================
q-point: ( 0.08  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.513   (  10.360    5.982   12.213)   17.096
   0.652   (  10.098    5.830   18.706)   22.043
   1.055   (  22.057   12.734   24.214)   35.142
   1.195   (   4.236    2.446    0.234)    4.897
   1.378   (  17.363   10.025   -0.294)   20.052
   3.848   (  -4.470   -2.581  -12.412)   13.443
   5.296   (  -1.008   -0.582    0.880)    1.459
   5.336   (   1.282    0.740   -0.057)    1.482
   5.467   (  -1.327   -0.766    6.060)    6.251
   5.493   (  -0.625   -0.361   -1.172)    1.376
   5.497   (  -0.686   -0.396    1.084)    1.343
   6.082   (  -1.781   -1.028   -3.012)    3.647
======================= Grid point 183 (24/84) =======================
q-point: ( 0.15  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.782   (  11.183    6.456    8.065)   15.224
   0.873   (   9.632    5.561   13.930)   17.826
   1.325   (   6.152    3.552    0.353)    7.113
   1.669   (  25.372   14.648   11.043)   31.309
   1.840   (  19.462   11.236    0.249)   22.474
   3.726   (  -4.622   -2.668  -11.251)   12.452
   5.270   (  -1.067   -0.616    0.799)    1.469
   5.303   (  -3.728   -2.152    1.386)    4.522
   5.452   (  -1.579   -0.911    3.564)    4.003
   5.466   (  -1.538   -0.888   -1.247)    2.170
   5.500   (   2.397    1.384    2.795)    3.934
   6.040   (  -1.966   -1.135   -3.731)    4.367
======================= Grid point 184 (25/84) =======================
q-point: ( 0.23  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.019   (   8.270    4.775    6.150)   11.359
   1.114   (  10.229    5.906    5.517)   13.036
   1.466   (   5.199    3.002    0.341)    6.013
   2.092   (   8.693    5.019   -1.494)   10.149
   2.350   (  23.103   13.339    8.209)   27.912
   3.682   (   1.792    1.034   -8.968)    9.203
   5.204   (  -4.039   -2.332    3.383)    5.762
   5.248   (  -0.698   -0.403    0.631)    1.024
   5.409   (  -1.675   -0.967    2.535)    3.189
   5.426   (  -1.637   -0.945   -1.223)    2.252
   5.552   (   1.206    0.696    1.625)    2.140
   5.979   (  -3.339   -1.928   -3.288)    5.067
======================= Grid point 185 (26/84) =======================
q-point: ( 0.31  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.168   (   4.146    2.393    5.242)    7.099
   1.353   (   9.626    5.558   -1.005)   11.161
   1.558   (   2.237    1.292    0.344)    2.606
   2.163   (  -1.187   -0.685   -5.280)    5.455
   2.872   (  19.913   11.497    8.689)   24.581
   3.808   (   7.638    4.410   -5.980)   10.656
   5.122   (  -2.612   -1.508    5.199)    6.011
   5.237   (  -0.261   -0.151    0.493)    0.578
   5.381   (  -0.596   -0.344    2.043)    2.156
   5.397   (  -0.669   -0.386   -1.188)    1.417
   5.550   (  -1.271   -0.734    0.235)    1.487
   5.872   (  -5.705   -3.294   -2.479)    7.039
======================= Grid point 186 (27/84) =======================
q-point: ( 0.38  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.227   (   1.050    0.606    4.861)    5.010
   1.566   (   7.875    4.546   -5.588)   10.673
   1.575   (  -0.568   -0.328    0.516)    0.835
   2.099   (  -3.317   -1.915   -7.587)    8.499
   3.294   (  15.072    8.702    8.744)   19.477
   3.984   (   5.851    3.378   -3.007)    7.395
   5.088   (   0.057    0.033    6.152)    6.153
   5.236   (   0.160    0.092    0.429)    0.467
   5.380   (   0.417    0.241    1.590)    1.661
   5.393   (   0.218    0.126   -1.220)    1.246
   5.505   (  -2.151   -1.242   -1.277)    2.793
   5.708   (  -7.685   -4.437   -1.149)    8.948
======================= Grid point 187 (28/84) =======================
q-point: ( 0.46  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.236   (   0.040    0.023    4.728)    4.728
   1.552   (  -0.812   -0.469    0.769)    1.213
   1.713   (   3.777    2.180   -8.003)    9.114
   2.028   (  -2.074   -1.197   -9.015)    9.328
   3.568   (   6.850    3.955    8.768)   11.809
   4.058   (   0.679    0.392   -0.000)    0.784
   5.127   (   2.435    1.406    5.357)    6.050
   5.241   (   0.157    0.091    0.402)    0.441
   5.387   (  -0.261   -0.151   -0.846)    0.898
   5.400   (   0.234    0.135   -1.273)    1.301
   5.464   (  -0.665   -0.384   -1.540)    1.721
   5.538   (  -5.133   -2.964    1.404)    6.091
======================= Grid point 195 (29/84) =======================
q-point: ( 0.08  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.755   (   8.438   14.615    8.350)   18.829
   0.804   (   5.191    8.991   15.757)   18.869
   1.301   (   5.192    8.992    1.713)   10.524
   1.516   (  14.667   25.404   13.064)   32.111
   1.673   (   9.866   17.088   -0.751)   19.745
   3.758   (  -3.056   -5.293  -11.655)   13.160
   5.284   (  -0.487   -0.843    0.675)    1.185
   5.320   (  -1.606   -2.782    0.539)    3.258
   5.454   (   0.278    0.482    4.518)    4.552
   5.471   (  -0.933   -1.617    1.716)    2.535
   5.505   (   0.479    0.830   -0.274)    0.997
   6.044   (  -1.477   -2.558   -3.578)    4.640
======================= Grid point 196 (30/84) =======================
q-point: ( 0.15  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.007   (   6.541    9.911    6.890)   13.729
   1.054   (   6.313   14.795    6.901)   17.503
   1.484   (   3.255   10.807    0.448)   11.295
   1.979   (   9.873   10.799    0.537)   14.641
   2.094   (  20.633   15.718    8.086)   27.169
   3.673   (  -0.948   -1.736  -10.002)   10.196
   5.236   (  -3.112   -3.897    2.876)    5.756
   5.266   (  -0.891    0.089    0.606)    1.081
   5.428   (  -1.941   -0.894    1.722)    2.744
   5.467   (  -0.814    0.583    0.450)    1.098
   5.526   (   1.086    0.545    1.237)    1.734
   5.979   (  -1.557   -4.608   -3.379)    5.922
======================= Grid point 197 (31/84) =======================
q-point: ( 0.23  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.218   (   4.343   10.402    1.590)   11.384
   1.285   (   3.925   12.570    4.252)   13.838
   1.624   (  -0.750   11.186    0.351)   11.217
   2.144   (   2.951    0.279   -4.015)    4.991
   2.623   (  20.646   12.722    9.056)   25.886
   3.719   (   5.509    2.468   -7.311)    9.481
   5.143   (  -2.740   -3.453    4.855)    6.558
   5.257   (  -0.783    1.106    0.537)    1.458
   5.387   (  -1.508   -0.758    1.442)    2.220
   5.444   (  -2.672    1.959   -0.720)    3.390
   5.539   (   0.722   -1.887    0.686)    2.133
   5.889   (  -2.317   -6.494   -2.789)    7.438
======================= Grid point 198 (32/84) =======================
q-point: ( 0.31  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.345   (  -3.510   12.703    2.843)   13.482
   1.484   (   6.759    7.326   -2.703)   10.327
   1.685   (  -4.629   10.003    0.361)   11.028
   2.132   (  -2.132   -1.802   -6.604)    7.170
   3.084   (  17.887    9.018    8.888)   21.915
   3.881   (   8.081    2.315   -4.363)    9.471
   5.075   (  -0.712   -3.035    6.292)    7.022
   5.259   (  -0.826    1.944    0.483)    2.167
   5.362   (  -0.216   -0.925    1.316)    1.622
   5.424   (  -1.655    1.792   -0.955)    2.620
   5.515   (  -0.826   -2.078   -0.639)    2.326
   5.756   (  -4.459   -7.418   -1.886)    8.858
======================= Grid point 199 (33/84) =======================
q-point: ( 0.38  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.378   (  -6.653   12.394    3.562)   14.510
   1.659   (   3.588    5.645   -5.759)    8.826
   1.681   (  -4.825    8.700   -0.667)    9.971
   2.062   (  -3.113   -1.533   -8.532)    9.211
   3.438   (  12.525    5.281    8.850)   16.220
   4.004   (   4.303   -1.412   -1.073)    4.655
   5.070   (   3.120   -1.718    5.964)    6.946
   5.266   (  -0.923    2.379    0.428)    2.587
   5.355   (   0.633   -2.049    0.272)    2.161
   5.422   (  -0.277    1.184   -0.575)    1.345
   5.477   (  -1.455   -0.829   -1.740)    2.415
   5.588   (  -5.356   -6.391   -0.282)    8.343
======================= Grid point 200 (34/84) =======================
q-point: ( 0.46  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.380   (  -7.100   12.297    3.691)   14.672
   1.663   (  -5.859   10.149    0.627)   11.736
   1.748   (  -0.749    1.297   -8.259)    8.393
   2.017   (  -0.087    0.150   -9.333)    9.335
   3.594   (   0.850   -1.472    8.790)    8.953
   4.022   (   1.992   -3.450    0.899)    4.084
   5.112   (   1.732   -3.000    5.228)    6.271
   5.273   (  -1.544    2.675    0.474)    3.125
   5.336   (   1.858   -3.219   -1.982)    4.212
   5.433   (  -0.785    1.360   -0.472)    1.640
   5.456   (   0.325   -0.563   -0.923)    1.129
   5.489   (   0.555   -0.961    0.888)    1.421
======================= Grid point 210 (35/84) =======================
q-point: ( 0.15  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.201   (   5.461    9.459    1.200)   10.988
   1.377   (   8.049   13.942    4.435)   16.698
   1.715   (   5.863   10.155    0.643)   11.744
   2.105   (   1.713    2.967   -3.957)    5.234
   2.476   (  11.282   19.541   10.006)   24.683
   3.680   (   1.566    2.712   -7.937)    8.533
   5.150   (  -2.484   -4.302    5.005)    7.052
   5.271   (   0.233    0.404    0.536)    0.710
   5.406   (  -0.741   -1.284    1.306)    1.976
   5.486   (   0.514    0.890   -0.515)    1.149
   5.517   (  -0.750   -1.299    0.381)    1.547
   5.863   (  -3.612   -6.256   -2.819)    7.755
======================= Grid point 211 (36/84) =======================
q-point: ( 0.23  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.396   (   5.965    8.077   -4.074)   10.836
   1.598   (   0.065   14.526    3.651)   14.978
   1.872   (  -1.141   11.331    0.532)   11.401
   2.130   (   0.356   -1.114   -5.962)    6.076
   2.882   (  14.939   12.607    9.549)   21.755
   3.779   (   5.349    3.171   -5.303)    8.172
   5.061   (  -2.134   -4.378    6.825)    8.385
   5.280   (   0.057    0.691    0.605)    0.920
   5.377   (  -1.506   -0.176    1.004)    1.819
   5.468   (  -1.661   -1.546   -1.415)    2.674
   5.508   (  -0.221    0.823   -0.257)    0.890
   5.730   (  -2.758   -8.150   -2.194)    8.879
======================= Grid point 212 (37/84) =======================
q-point: ( 0.31  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.568   (   3.177    7.822   -6.528)   10.672
   1.703   (  -3.156   13.013    2.062)   13.548
   1.917   (  -6.060   11.052    0.050)   12.604
   2.099   (  -2.284   -1.242   -7.587)    8.020
   3.256   (  14.151    7.884    9.232)   18.645
   3.900   (   5.306   -0.413   -2.010)    5.688
   5.002   (   0.873   -3.719    6.913)    7.899
   5.284   (   0.273   -0.868    0.491)    1.035
   5.353   (  -1.453    0.682    0.634)    1.725
   5.429   (  -0.638   -1.338   -1.254)    1.941
   5.503   (  -1.403    0.937   -0.950)    1.936
   5.594   (  -3.311   -7.210   -1.319)    8.042
======================= Grid point 213 (38/84) =======================
q-point: ( 0.38  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.645   (  -4.874   10.738   -2.208)   11.997
   1.789   (  -0.379    8.146   -4.117)    9.135
   1.908   (  -7.101   10.935   -0.526)   13.049
   2.037   (  -2.636   -0.681   -8.604)    9.024
   3.520   (   7.469    2.912    8.904)   11.981
   3.936   (   2.264   -4.652    1.377)    5.354
   5.022   (   4.028   -2.689    5.424)    7.271
   5.264   (   0.928   -4.834   -0.405)    4.939
   5.338   (  -1.883    2.299    0.209)    2.979
   5.421   (   0.985   -0.794   -0.117)    1.270
   5.464   (  -2.026   -2.303   -0.553)    3.117
   5.504   (  -1.664    0.342   -1.211)    2.086
======================= Grid point 224 (39/84) =======================
q-point: ( 0.23  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.563   (   4.587    7.945   -7.741)   12.004
   1.828   (   4.254    7.368    2.108)    8.765
   2.053   (   3.243    5.617    0.800)    6.535
   2.106   (  -0.693   -1.200   -7.595)    7.720
   3.155   (   7.835   13.570    9.652)   18.404
   3.843   (   1.531    2.652   -2.444)    3.918
   4.973   (  -2.021   -3.500    7.381)    8.415
   5.272   (  -0.942   -1.632    0.638)    1.990
   5.378   (   0.079    0.137    0.475)    0.500
   5.424   (  -1.021   -1.769   -1.352)    2.450
   5.531   (   0.318    0.550   -0.992)    1.178
   5.575   (  -3.645   -6.313   -1.470)    7.436
======================= Grid point 225 (40/84) =======================
q-point: ( 0.31  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.717   (   3.395    6.929   -9.891)   12.545
   1.923   (  -0.947    7.009    0.337)    7.081
   2.059   (  -3.419   -0.084   -7.565)    8.302
   2.120   (  -2.088    5.284    0.238)    5.687
   3.425   (   8.468    8.584    9.165)   15.145
   3.866   (   0.666   -2.877    1.490)    3.307
   4.943   (   1.783   -1.257    5.711)    6.114
   5.211   (  -1.468   -5.451    0.434)    5.662
   5.379   (  -0.668    0.964   -0.094)    1.177
   5.402   (   0.104   -1.101   -0.856)    1.398
   5.476   (  -2.006   -2.582   -0.580)    3.321
   5.531   (  -0.721    0.282   -1.026)    1.286
======================= Grid point 226 (41/84) =======================
q-point: ( 0.38  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.794   (  -2.797    4.845   -9.686)   11.186
   1.974   (  -4.175    7.232    1.511)    8.486
   1.991   (   0.522   -0.903  -10.948)   10.997
   2.124   (  -3.769    6.529    0.640)    7.566
   3.572   (  -1.305    2.260    8.600)    8.987
   3.835   (   2.629   -4.553    3.851)    6.517
   4.975   (   0.814   -1.410    4.273)    4.573
   5.148   (   2.970   -5.145    0.537)    5.965
   5.377   (  -0.561    0.972   -0.582)    1.264
   5.398   (   0.615   -1.065   -0.178)    1.243
   5.445   (  -0.276    0.478   -0.406)    0.685
   5.524   (  -0.475    0.823   -0.956)    1.348
======================= Grid point 239 (42/84) =======================
q-point: ( 0.31  0.31  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.852   (   3.248    5.625  -11.383)   13.106
   1.967   (  -1.262   -2.187  -10.155)   10.465
   2.057   (   0.363    0.629    1.217)    1.417
   2.194   (   0.895    1.550    0.728)    1.932
   3.603   (   4.519    7.827    8.132)   12.157
   3.778   (  -3.088   -5.348    5.208)    8.078
   4.964   (   1.769    3.064    3.422)    4.923
   5.086   (  -3.327   -5.763    1.256)    6.772
   5.385   (  -0.054   -0.094   -0.512)    0.524
   5.386   (  -0.176   -0.304   -0.999)    1.059
   5.455   (  -0.047   -0.081    0.027)    0.098
   5.533   (  -0.052   -0.090   -0.799)    0.806
======================= Grid point 362 (43/84) =======================
q-point: ( 0.00  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.683   (   0.000    0.000   14.015)   14.015
   0.683   (  -0.000   -0.000   14.015)   14.015
   1.132   (   0.000    0.000   -1.409)    1.409
   1.132   (   0.000    0.000   -1.409)    1.409
   1.635   (   0.000   -0.000   38.193)   38.193
   3.545   (   0.000    0.000  -22.188)   22.188
   5.336   (   0.000   -0.000    1.669)    1.669
   5.336   (   0.000    0.000    1.669)    1.669
   5.467   (  -0.000   -0.000   -1.948)    1.948
   5.467   (   0.000   -0.000   -1.948)    1.948
   5.631   (   0.000   -0.000    7.665)    7.665
   6.045   (   0.000   -0.000   -4.620)    4.620
======================= Grid point 363 (44/84) =======================
q-point: ( 0.08  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.774   (   6.940    4.007   12.443)   14.800
   0.954   (  15.638    9.029   12.092)   21.732
   1.188   (   4.439    2.563   -1.224)    5.270
   1.357   (  16.120    9.307   -1.476)   18.672
   1.691   (   6.910    3.990   34.168)   35.087
   3.499   (  -3.468   -2.002  -21.406)   21.778
   5.322   (  -1.021   -0.589    1.692)    2.062
   5.348   (   0.471    0.272    1.289)    1.399
   5.459   (  -0.745   -0.430   -2.059)    2.232
   5.489   (   1.016    0.586   -1.527)    1.925
   5.615   (  -0.923   -0.533    7.783)    7.856
   6.005   (  -2.638   -1.523   -4.539)    5.466
======================= Grid point 364 (45/84) =======================
q-point: ( 0.15  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.974   (   8.976    5.182    9.926)   14.351
   1.209   (   5.712    3.298   16.788)   18.038
   1.323   (   6.269    3.620   -0.909)    7.296
   1.763   (  15.568    8.988   -2.084)   18.097
   2.045   (  20.979   12.112   19.601)   31.161
   3.416   (  -2.265   -1.308  -18.686)   18.868
   5.294   (  -1.222   -0.706    1.595)    2.130
   5.327   (  -2.322   -1.341    1.362)    3.007
   5.431   (  -1.551   -0.895   -2.158)    2.804
   5.489   (  -1.183   -0.683   -0.556)    1.475
   5.610   (   0.634    0.366    7.180)    7.217
   5.949   (  -1.877   -1.084   -4.703)    5.178
======================= Grid point 365 (46/84) =======================
q-point: ( 0.23  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.172   (   7.073    4.084    8.133)   11.526
   1.304   (   3.212    1.854   12.697)   13.228
   1.465   (   5.218    3.013   -0.724)    6.068
   2.015   (   4.529    2.615   -6.027)    7.979
   2.582   (  22.449   12.961   13.507)   29.230
   3.449   (   5.729    3.308  -13.189)   14.755
   5.260   (  -2.697   -1.557    2.216)    3.822
   5.268   (  -0.899   -0.519    1.322)    1.681
   5.391   (  -1.560   -0.900   -2.067)    2.741
   5.445   (  -1.985   -1.146   -0.313)    2.313
   5.622   (  -0.217   -0.125    5.454)    5.460
   5.902   (  -2.486   -1.435   -3.774)    4.742
======================= Grid point 366 (47/84) =======================
q-point: ( 0.31  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.300   (   3.562    2.056    7.110)    8.214
   1.382   (   3.420    1.974    4.600)    6.063
   1.558   (   2.350    1.357   -0.585)    2.775
   2.007   (  -3.711   -2.143   -9.763)   10.662
   3.088   (  19.505   11.261   11.661)   25.363
   3.666   (  11.040    6.374   -7.467)   14.774
   5.211   (  -1.338   -0.773    3.102)    3.466
   5.253   (  -0.370   -0.214    1.095)    1.175
   5.364   (  -0.655   -0.378   -1.968)    2.109
   5.409   (  -0.899   -0.519   -0.742)    1.276
   5.586   (  -2.789   -1.610    4.161)    5.261
   5.819   (  -4.512   -2.605   -2.347)    5.715
======================= Grid point 367 (48/84) =======================
q-point: ( 0.38  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.349   (   0.804    0.464    6.601)    6.666
   1.467   (   3.583    2.069   -2.760)    4.974
   1.580   (  -0.273   -0.158   -0.270)    0.416
   1.895   (  -4.883   -2.819  -11.699)   12.987
   3.504   (  14.910    8.608   10.924)   20.390
   3.918   (   8.970    5.179   -3.313)   10.874
   5.195   (  -0.060   -0.035    3.372)    3.372
   5.250   (   0.120    0.069    1.004)    1.014
   5.359   (   0.159    0.092   -1.998)    2.006
   5.396   (  -0.736   -0.425   -2.167)    2.328
   5.505   (  -3.207   -1.852    3.583)    5.153
   5.693   (  -5.610   -3.239   -0.294)    6.484
======================= Grid point 368 (49/84) =======================
q-point: ( 0.46  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.355   (  -0.041   -0.024    6.362)    6.363
   1.542   (   2.168    1.252   -7.301)    7.718
   1.564   (  -0.615   -0.355    0.122)    0.721
   1.799   (  -2.623   -1.515  -12.447)   12.810
   3.775   (   6.718    3.879   10.417)   12.988
   4.060   (   2.635    1.521    0.048)    3.043
   5.216   (   1.832    1.058    2.032)    2.933
   5.254   (   0.144    0.083    0.976)    0.990
   5.338   (  -3.522   -2.033   -3.779)    5.551
   5.365   (   0.196    0.113   -2.070)    2.083
   5.474   (   0.363    0.210    3.332)    3.358
   5.580   (  -3.078   -1.777    2.213)    4.187
======================= Grid point 376 (50/84) =======================
q-point: ( 0.08  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.954   (   6.713   11.627   10.166)   16.840
   1.146   (   4.240    7.345   15.270)   17.467
   1.311   (   4.626    8.013   -0.364)    9.260
   1.644   (   9.211   15.955   -2.352)   18.572
   1.906   (   9.997   17.316   24.103)   31.316
   3.434   (  -1.976   -3.423  -19.685)   20.077
   5.307   (  -0.574   -0.994    1.579)    1.952
   5.333   (  -1.027   -1.779    1.060)    2.312
   5.458   (  -0.059   -0.103   -2.053)    2.057
   5.486   (  -0.505   -0.875   -0.649)    1.201
   5.606   (  -0.051   -0.088    7.425)    7.426
   5.956   (  -1.578   -2.733   -4.664)    5.632
======================= Grid point 377 (51/84) =======================
q-point: ( 0.15  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.186   (   3.008   12.834    8.800)   15.849
   1.265   (   4.648    3.254   13.845)   14.962
   1.485   (   3.032   11.003   -0.654)   11.432
   1.928   (   7.989    5.620   -4.626)   10.808
   2.355   (  17.997   16.887   15.665)   29.231
   3.406   (   1.705    1.054  -15.577)   15.705
   5.282   (  -1.737   -2.001    1.734)    3.167
   5.286   (  -1.071   -0.402    1.470)    1.863
   5.431   (  -2.726    0.376   -1.589)    3.178
   5.465   (  -0.754   -1.304   -0.914)    1.762
   5.608   (   0.401   -0.854    6.124)    6.196
   5.898   (  -1.168   -3.872   -4.087)    5.750
======================= Grid point 378 (52/84) =======================
q-point: ( 0.23  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.300   (   2.279    3.719    6.983)    8.233
   1.409   (   0.764   11.583    7.307)   13.717
   1.624   (  -0.825   11.299   -0.581)   11.344
   2.018   (   0.468   -1.879   -8.318)    8.540
   2.852   (  19.451   12.762   12.538)   26.427
   3.539   (   8.897    5.116   -9.779)   14.176
   5.229   (  -1.901   -1.362    3.179)    3.946
   5.274   (  -0.874    1.052    1.167)    1.798
   5.390   (  -1.882    0.340   -1.467)    2.411
   5.430   (  -1.955   -0.218   -1.377)    2.401
   5.587   (  -0.606   -3.148    4.531)    5.551
   5.826   (  -1.586   -5.698   -2.916)    6.595
======================= Grid point 379 (53/84) =======================
q-point: ( 0.31  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.377   (   2.334    3.006    1.185)    3.985
   1.505   (  -2.876   11.048    4.561)   12.294
   1.685   (  -4.450    9.734   -0.670)   10.724
   1.951   (  -4.104   -2.688  -10.596)   11.676
   3.300   (  17.389    9.110   11.316)   22.658
   3.782   (  11.034    4.552   -5.066)   12.966
   5.189   (  -0.822   -1.434    3.998)    4.326
   5.274   (  -0.882    1.883    1.013)    2.313
   5.373   (  -0.987    0.691   -1.117)    1.643
   5.402   (  -0.657   -0.242   -1.704)    1.842
   5.524   (  -2.121   -3.711    2.880)    5.154
   5.719   (  -3.089   -6.611   -1.445)    7.439
======================= Grid point 380 (54/84) =======================
q-point: ( 0.38  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.436   (  -0.761    5.016   -0.035)    5.073
   1.556   (  -2.100    8.252   -0.420)    8.525
   1.680   (  -6.010    8.463   -0.950)   10.423
   1.846   (  -4.550   -1.414  -11.390)   12.346
   3.647   (  12.226    5.180   10.609)   16.996
   3.982   (   7.017    0.926   -1.205)    7.180
   5.178   (   1.519   -1.352    3.392)    3.954
   5.279   (  -0.803    1.807    0.683)    2.092
   5.345   (  -1.783   -2.039   -1.647)    3.170
   5.404   (  -0.550    1.522   -2.077)    2.633
   5.459   (  -0.747   -2.020    2.228)    3.099
   5.594   (  -2.875   -5.769    0.742)    6.489
======================= Grid point 381 (55/84) =======================
q-point: ( 0.46  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.453   (  -3.695    6.399    1.397)    7.520
   1.590   (  -3.036    5.259   -4.161)    7.361
   1.665   (  -4.842    8.386   -1.428)    9.788
   1.788   (  -0.750    1.298  -11.200)   11.300
   3.799   (   0.902   -1.563   10.183)   10.342
   4.044   (   1.601   -2.772    0.967)    3.344
   5.205   (   1.205   -2.087    2.473)    3.454
   5.280   (   1.265   -2.192   -3.189)    4.071
   5.292   (  -1.542    2.670    1.392)    3.383
   5.402   (  -1.480    2.564   -2.147)    3.657
   5.457   (   1.018   -1.763    1.798)    2.716
   5.529   (   1.607   -2.784    2.301)    3.953
======================= Grid point 391 (56/84) =======================
q-point: ( 0.15  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.295   (   1.407    2.437    8.312)    8.775
   1.491   (   7.202   12.474    6.176)   15.672
   1.721   (   6.008   10.406   -0.216)   12.018
   1.984   (   0.147    0.254   -7.953)    7.958
   2.725   (  10.693   18.521   13.484)   25.282
   3.481   (   3.641    6.307  -10.958)   13.158
   5.241   (  -1.039   -1.800    3.383)    3.971
   5.289   (   0.230    0.399    1.270)    1.351
   5.419   (  -0.945   -1.637   -0.599)    1.983
   5.449   (  -0.194   -0.336   -2.613)    2.642
   5.572   (  -1.391   -2.409    4.655)    5.423
   5.798   (  -3.145   -5.448   -3.040)    6.987
======================= Grid point 392 (57/84) =======================
q-point: ( 0.23  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.354   (   2.402    2.563    1.060)    3.669
   1.684   (  -0.084   11.923    4.012)   12.580
   1.872   (  -2.399   10.756   -1.276)   11.094
   1.971   (   0.081   -1.348   -8.382)    8.490
   3.113   (  14.102   12.587   12.041)   22.411
   3.655   (   7.896    6.062   -6.431)   11.851
   5.198   (  -1.794   -1.897    5.524)    6.110
   5.298   (  -0.022    0.712    1.175)    1.374
   5.388   (  -0.840   -0.551   -0.490)    1.117
   5.427   (  -1.776   -0.897   -2.569)    3.249
   5.520   (  -1.376   -2.079    2.071)    3.241
   5.683   (  -2.068   -7.423   -2.016)    7.965
======================= Grid point 393 (58/84) =======================
q-point: ( 0.31  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.432   (   2.995    2.857   -4.918)    6.428
   1.748   (  -3.647    8.084   -1.117)    8.938
   1.868   (  -5.876    5.663   -4.759)    9.447
   1.966   (  -5.018    7.388   -2.320)    9.227
   3.473   (  13.494    7.805   10.984)   19.070
   3.855   (   8.112    2.610   -2.396)    8.852
   5.146   (  -0.617   -2.828    5.990)    6.653
   5.294   (  -0.442   -1.452    0.321)    1.552
   5.367   (  -1.553   -0.486    0.584)    1.729
   5.396   (   0.370   -0.334   -1.586)    1.662
   5.483   (  -1.864    0.318   -0.552)    1.970
   5.570   (  -2.108   -6.891   -0.823)    7.253
======================= Grid point 394 (59/84) =======================
q-point: ( 0.38  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.506   (   1.269    2.629   -8.144)    8.651
   1.746   (  -1.381    1.354  -10.406)   10.584
   1.819   (  -7.909   10.340    1.131)   13.067
   1.945   (  -7.024   10.500   -0.457)   12.641
   3.724   (   7.137    2.513   10.022)   12.557
   3.965   (   3.646   -2.090    1.164)    4.361
   5.138   (   2.686   -2.601    4.732)    6.031
   5.246   (  -0.791   -4.315   -1.296)    4.575
   5.348   (  -1.297    1.657    1.015)    2.336
   5.415   (   1.639   -0.176   -0.210)    1.662
   5.456   (  -1.711    0.148   -0.819)    1.903
   5.483   (  -0.983   -3.123   -0.210)    3.281
======================= Grid point 405 (60/84) =======================
q-point: ( 0.23  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.418   (   2.133    3.694   -4.862)    6.468
   1.808   (  -0.108   -0.187   -8.173)    8.175
   1.976   (   2.118    3.669   -0.686)    4.292
   2.068   (   3.280    5.681    0.389)    6.572
   3.380   (   7.525   13.033   11.328)   18.837
   3.787   (   3.667    6.351   -2.977)    7.915
   5.139   (  -2.094   -3.626    7.754)    8.812
   5.283   (  -1.222   -2.117    0.365)    2.472
   5.383   (   0.096    0.167   -0.228)    0.298
   5.392   (  -1.014   -1.756   -1.548)    2.551
   5.509   (  -0.004   -0.006   -0.979)    0.979
   5.545   (  -3.082   -5.338   -1.184)    6.276
======================= Grid point 406 (61/84) =======================
q-point: ( 0.31  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.503   (   2.502    3.980   -9.339)   10.456
   1.764   (  -1.946   -2.566  -13.455)   13.835
   2.027   (  -2.895    6.881    3.012)    8.049
   2.131   (  -2.434    6.218    0.423)    6.691
   3.633   (   7.758    7.709   10.124)   14.904
   3.894   (   3.054    1.046    0.879)    3.346
   5.079   (   0.141   -2.965    6.733)    7.358
   5.218   (  -1.727   -4.818    0.094)    5.120
   5.365   (  -0.532   -0.241   -1.073)    1.221
   5.391   (   0.774   -0.753   -0.204)    1.099
   5.466   (  -1.675   -1.881   -0.279)    2.535
   5.505   (  -0.814    0.564   -1.452)    1.757
======================= Grid point 407 (62/84) =======================
q-point: ( 0.38  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.559   (  -1.434    2.484  -11.226)   11.587
   1.720   (   1.016   -1.760  -14.210)   14.355
   2.032   (  -4.459    7.723    3.473)    9.570
   2.137   (  -3.940    6.825    0.446)    7.893
   3.763   (  -0.824    1.428    9.047)    9.196
   3.915   (   1.352   -2.342    3.395)    4.340
   5.079   (   1.431   -2.479    5.219)    5.952
   5.159   (   2.034   -3.524    0.487)    4.098
   5.361   (   0.047   -0.082   -0.858)    0.863
   5.394   (   1.201   -2.080   -0.308)    2.422
   5.441   (  -0.763    1.322    0.080)    1.528
   5.499   (  -0.678    1.174   -1.416)    1.960
======================= Grid point 420 (63/84) =======================
q-point: ( 0.31  0.31  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.592   (   2.431    4.211  -12.460)   13.375
   1.696   (  -2.070   -3.586  -14.236)   14.826
   2.114   (   1.044    1.809    3.353)    3.950
   2.209   (   0.906    1.570    0.547)    1.894
   3.781   (   3.413    5.912    8.294)   10.742
   3.886   (  -1.134   -1.965    4.601)    5.130
   5.052   (   0.365    0.633    4.713)    4.769
   5.118   (  -2.381   -4.125    1.711)    5.061
   5.361   (  -0.030   -0.052   -1.388)    1.389
   5.370   (  -0.404   -0.700   -0.930)    1.232
   5.458   (   0.098    0.169    0.299)    0.357
   5.511   (   0.045    0.078   -1.210)    1.213
======================= Grid point 543 (64/84) =======================
q-point: ( 0.00  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.924   (   0.000   -0.000    9.570)    9.570
   0.924   (  -0.000    0.000    9.570)    9.570
   1.070   (   0.000    0.000   -4.937)    4.937
   1.070   (   0.000    0.000   -4.937)    4.937
   2.375   (   0.000    0.000   34.363)   34.363
   3.021   (  -0.000    0.000  -29.081)   29.081
   5.376   (   0.000   -0.000    2.214)    2.214
   5.376   (   0.000    0.000    2.214)    2.214
   5.423   (   0.000   -0.000   -2.339)    2.339
   5.423   (   0.000    0.000   -2.339)    2.339
   5.789   (   0.000    0.000    7.599)    7.599
   5.932   (   0.000   -0.000   -6.432)    6.432
======================= Grid point 544 (65/84) =======================
q-point: ( 0.08  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.994   (   5.470    3.158    8.943)   10.949
   1.132   (   4.795    2.768   -4.635)    7.221
   1.163   (  16.095    9.293    8.480)   20.429
   1.296   (  15.987    9.230   -4.637)   19.033
   2.378   (   1.133    0.654   32.563)   32.589
   2.995   (  -1.707   -0.986  -27.905)   27.975
   5.363   (  -0.987   -0.570    2.286)    2.554
   5.387   (   0.449    0.259    2.483)    2.537
   5.412   (  -0.886   -0.511   -2.450)    2.655
   5.443   (   1.059    0.611   -2.752)    3.012
   5.769   (  -1.265   -0.730    7.113)    7.261
   5.900   (  -2.201   -1.271   -5.802)    6.334
======================= Grid point 545 (66/84) =======================
q-point: ( 0.15  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.157   (   7.489    4.324    7.649)   11.545
   1.276   (   6.658    3.844   -3.977)    8.655
   1.500   (   9.233    5.330   11.529)   15.703
   1.679   (  13.145    7.589   -6.254)   16.416
   2.502   (  10.641    6.144   24.133)   27.081
   2.988   (   2.550    1.472  -23.053)   23.240
   5.333   (  -1.396   -0.806    2.247)    2.765
   5.371   (  -1.934   -1.117    3.056)    3.785
   5.383   (  -1.514   -0.874   -2.499)    3.049
   5.443   (  -1.336   -0.771   -3.510)    3.834
   5.749   (  -0.218   -0.126    6.143)    6.149
   5.857   (  -1.153   -0.665   -4.618)    4.806
======================= Grid point 546 (67/84) =======================
q-point: ( 0.23  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.325   (   6.076    3.508    6.490)    9.558
   1.426   (   5.467    3.156   -3.379)    7.160
   1.588   (  -0.456   -0.263   14.188)   14.198
   1.846   (   0.866    0.500  -10.692)   10.739
   2.871   (  18.442   10.647   14.588)   25.812
   3.166   (  12.104    6.988  -14.403)   20.070
   5.302   (  -1.160   -0.670    1.971)    2.383
   5.313   (  -2.431   -1.403    3.293)    4.327
   5.346   (  -1.399   -0.807   -2.305)    2.814
   5.392   (  -2.327   -1.344   -4.061)    4.869
   5.742   (  -0.919   -0.531    5.437)    5.539
   5.831   (  -1.472   -0.850   -3.528)    3.916
======================= Grid point 547 (68/84) =======================
q-point: ( 0.31  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.435   (   3.002    1.733    5.741)    6.706
   1.524   (   2.594    1.498   -2.964)    4.214
   1.542   (  -2.559   -1.477   10.591)   10.996
   1.781   (  -4.896   -2.827  -12.107)   13.362
   3.317   (  17.983   10.383   10.476)   23.258
   3.507   (  15.033    8.679   -8.359)   19.267
   5.268   (  -1.279   -0.739    3.136)    3.466
   5.281   (  -0.499   -0.288    1.732)    1.826
   5.321   (  -0.602   -0.348   -2.123)    2.234
   5.347   (  -1.419   -0.819   -4.303)    4.605
   5.693   (  -3.076   -1.776    5.167)    6.270
   5.773   (  -3.416   -1.972   -2.690)    4.775
======================= Grid point 548 (69/84) =======================
q-point: ( 0.38  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.473   (   0.477    0.275    5.253)    5.282
   1.488   (  -1.680   -0.970    5.040)    5.400
   1.553   (   0.103    0.059   -2.547)    2.550
   1.657   (  -4.817   -2.781  -10.875)   12.215
   3.706   (  13.969    8.065    8.444)   18.207
   3.840   (  11.935    6.891   -4.866)   14.615
   5.244   (  -0.909   -0.525    2.473)    2.686
   5.277   (   0.092    0.053    1.666)    1.669
   5.315   (  -1.632   -0.942   -4.388)    4.776
   5.316   (   0.106    0.061   -2.110)    2.114
   5.611   (  -3.219   -1.859    4.977)    6.212
   5.679   (  -3.974   -2.295   -1.663)    4.881
======================= Grid point 549 (70/84) =======================
q-point: ( 0.46  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.466   (  -0.263   -0.152    0.562)    0.639
   1.473   (  -0.183   -0.105    4.921)    4.925
   1.546   (  -0.383   -0.221   -2.146)    2.191
   1.569   (  -2.178   -1.257   -9.081)    9.423
   3.956   (   6.040    3.487    6.827)    9.760
   4.045   (   4.540    2.621   -2.097)    5.646
   5.228   (  -0.050   -0.029    0.430)    0.434
   5.262   (  -2.318   -1.338   -3.061)    4.066
   5.281   (   0.141    0.081    1.670)    1.677
   5.321   (   0.159    0.092   -2.159)    2.167
   5.562   (  -0.779   -0.450    3.927)    4.028
   5.605   (  -1.781   -1.028   -0.199)    2.066
======================= Grid point 557 (71/84) =======================
q-point: ( 0.08  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.140   (   5.610    9.717    7.731)   13.625
   1.257   (   5.025    8.703   -3.573)   10.665
   1.415   (   6.565   11.370    9.649)   16.294
   1.567   (   8.236   14.265   -5.441)   17.347
   2.441   (   3.979    6.892   27.281)   28.418
   2.977   (   0.189    0.328  -24.906)   24.909
   5.345   (  -0.691   -1.197    2.092)    2.507
   5.364   (  -1.107   -1.918    1.375)    2.607
   5.417   (  -0.005   -0.009   -1.276)    1.276
   5.447   (  -0.292   -0.506   -2.782)    2.843
   5.748   (  -0.527   -0.913    6.314)    6.402
   5.860   (  -1.196   -2.072   -4.893)    5.447
======================= Grid point 558 (72/84) =======================
q-point: ( 0.15  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.344   (   3.124   12.468    6.589)   14.444
   1.445   (   3.126   11.469   -3.258)   12.325
   1.559   (   3.047    0.167   13.559)   13.899
   1.791   (   5.393    2.354   -9.032)   10.779
   2.698   (  12.846   12.465   17.758)   25.213
   3.063   (   6.649    6.228  -17.793)   19.990
   5.318   (  -1.119   -0.750    1.796)    2.245
   5.324   (  -0.978   -1.684    1.843)    2.681
   5.392   (  -2.705    0.429   -1.426)    3.087
   5.419   (  -1.579   -1.380   -3.125)    3.764
   5.731   (  -0.111   -1.817    5.420)    5.717
   5.822   (  -0.504   -2.805   -3.711)    4.679
======================= Grid point 559 (73/84) =======================
q-point: ( 0.23  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.469   (  -0.964    5.695    7.805)    9.710
   1.558   (  -0.752    5.663    3.219)    6.557
   1.608   (  -0.642    8.564    3.404)    9.238
   1.814   (  -1.316   -2.513  -11.345)   11.695
   3.105   (  17.109   11.446   11.964)   23.809
   3.333   (  13.384    8.643  -10.566)   19.118
   5.287   (  -1.465   -0.676    2.991)    3.398
   5.304   (  -1.042    0.915    1.691)    2.187
   5.350   (  -2.002    0.671   -2.075)    2.960
   5.370   (  -2.241   -0.768   -3.838)    4.510
   5.694   (  -0.899   -4.085    4.962)    6.490
   5.772   (  -0.927   -4.796   -2.846)    5.653
======================= Grid point 560 (74/84) =======================
q-point: ( 0.31  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.457   (  -0.740   -1.753    6.832)    7.092
   1.586   (  -3.096    5.173   -0.620)    6.061
   1.653   (  -4.921    8.910   -0.124)   10.179
   1.746   (  -5.333    2.917   -7.875)    9.948
   3.514   (  16.213    8.419    9.271)   20.486
   3.670   (  14.072    6.819   -6.283)   16.852
   5.250   (  -1.527   -0.971    2.786)    3.323
   5.295   (  -1.241    0.537    0.384)    1.406
   5.331   (  -1.457    0.919   -2.341)    2.907
   5.346   (  -0.486    1.554   -2.407)    2.906
   5.617   (  -1.834   -5.270    4.781)    7.348
   5.688   (  -1.973   -5.896   -2.160)    6.582
======================= Grid point 561 (75/84) =======================
q-point: ( 0.38  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.434   (   0.262   -1.554    2.647)    3.081
   1.555   (  -2.055   -0.882   -7.274)    7.610
   1.634   (  -7.068    9.772    2.458)   12.308
   1.698   (  -6.981    8.856   -3.203)   11.723
   3.836   (  11.350    4.404    7.438)   14.267
   3.942   (   9.549    2.929   -3.201)   10.488
   5.218   (  -0.406   -1.402    1.639)    2.195
   5.264   (  -1.217   -2.028   -2.526)    3.461
   5.323   (  -1.748    2.593    0.752)    3.217
   5.356   (  -1.241    2.716   -2.137)    3.673
   5.542   (  -0.209   -4.474    4.229)    6.160
   5.597   (  -1.099   -5.274   -1.027)    5.484
======================= Grid point 562 (76/84) =======================
q-point: ( 0.46  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.431   (   0.638   -1.105    0.421)    1.344
   1.526   (   0.542   -0.939   -8.479)    8.548
   1.620   (  -6.400   11.086    3.588)   13.294
   1.680   (  -5.950   10.306   -2.273)   12.116
   3.971   (   1.110   -1.923    6.236)    6.620
   4.045   (   1.338   -2.318   -1.309)    2.979
   5.215   (   0.584   -1.011   -0.308)    1.207
   5.226   (   0.703   -1.218   -1.718)    2.220
   5.321   (  -1.693    2.932    1.529)    3.715
   5.358   (  -1.620    2.806   -2.109)    3.867
   5.521   (   1.920   -3.326    3.275)    5.047
   5.557   (   1.973   -3.417   -0.024)    3.946
======================= Grid point 572 (77/84) =======================
q-point: ( 0.15  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.514   (  -1.524   -2.640   12.100)   12.478
   1.605   (   5.808   10.060    4.380)   12.415
   1.699   (   6.298   10.908   -1.873)   12.734
   1.786   (  -0.753   -1.304  -11.014)   11.116
   2.997   (   9.185   15.909   13.008)   22.510
   3.249   (   6.605   11.440  -11.835)   17.736
   5.297   (  -0.347   -0.601    2.503)    2.597
   5.319   (   0.264    0.458    1.248)    1.355
   5.382   (  -1.091   -1.890   -2.192)    3.093
   5.393   (  -0.423   -0.732   -2.557)    2.693
   5.673   (  -2.120   -3.672    4.522)    6.199
   5.744   (  -2.565   -4.442   -2.641)    5.769
======================= Grid point 573 (78/84) =======================
q-point: ( 0.23  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.455   (  -1.459   -2.639    8.439)    8.961
   1.654   (  -2.836   -2.055   -9.145)    9.792
   1.833   (  -0.724   10.966    1.697)   11.120
   1.881   (   0.352   10.249   -2.208)   10.490
   3.341   (  12.759   11.484   10.076)   19.905
   3.517   (  10.671    9.188   -7.431)   15.922
   5.279   (  -1.837   -0.416    3.020)    3.559
   5.323   (  -1.209   -0.285    0.526)    1.349
   5.353   (   0.049    0.427   -1.884)    1.932
   5.375   (  -1.127    0.688   -2.247)    2.606
   5.588   (  -1.850   -5.407    3.725)    6.822
   5.647   (  -1.674   -6.642   -2.024)    7.143
======================= Grid point 574 (79/84) =======================
q-point: ( 0.31  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.417   (  -0.208   -1.632    3.807)    4.148
   1.573   (  -2.482   -3.165  -10.587)   11.325
   1.888   (  -6.141   11.951    3.086)   13.787
   1.938   (  -5.844   11.145   -1.640)   12.690
   3.670   (  12.403    6.857    7.997)   16.273
   3.792   (  10.635    5.098   -4.191)   12.516
   5.226   (  -2.188   -1.428    2.571)    3.665
   5.275   (  -2.178   -2.070   -2.180)    3.712
   5.370   (  -0.202    1.862    0.411)    1.918
   5.384   (   0.843    1.049   -0.543)    1.451
   5.504   (  -1.562   -4.291    2.532)    5.221
   5.551   (  -1.219   -6.509   -1.523)    6.795
======================= Grid point 575 (80/84) =======================
q-point: ( 0.38  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.410   (   0.570   -0.692   -0.098)    0.902
   1.512   (  -0.832   -2.231   -9.528)    9.821
   1.884   (  -7.667   12.375    3.161)   14.897
   1.931   (  -7.404   11.734   -1.370)   13.942
   3.893   (   6.468    1.409    6.116)    9.013
   3.968   (   5.191   -0.126   -1.400)    5.377
   5.191   (   0.274   -1.344    1.201)    1.823
   5.216   (  -1.289   -2.515   -1.453)    3.178
   5.379   (  -0.915    1.912    1.679)    2.704
   5.410   (   0.256    1.673   -1.171)    2.059
   5.451   (  -0.075   -2.879    1.532)    3.262
   5.481   (   0.604   -5.061   -0.708)    5.145
======================= Grid point 586 (81/84) =======================
q-point: ( 0.23  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.410   (  -0.824   -1.428    4.360)    4.662
   1.579   (  -2.237   -3.875  -11.346)   12.197
   2.020   (   3.552    6.153    2.782)    7.630
   2.061   (   3.414    5.914   -1.133)    6.922
   3.583   (   6.744   11.681    8.315)   15.845
   3.713   (   5.470    9.474   -4.666)   11.894
   5.251   (  -1.592   -2.758    3.027)    4.394
   5.282   (  -1.659   -2.874   -0.572)    3.368
   5.372   (  -0.474   -0.821   -0.301)    0.995
   5.374   (   0.467    0.808   -0.375)    1.006
   5.504   (  -1.325   -2.295    0.673)    2.734
   5.525   (  -2.581   -4.471   -1.049)    5.268
======================= Grid point 587 (82/84) =======================
q-point: ( 0.31  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.398   (   0.305   -0.089   -0.026)    0.318
   1.508   (  -1.618   -2.901  -10.182)   10.710
   2.090   (  -2.628    6.900    2.709)    7.865
   2.128   (  -2.537    6.523   -0.955)    7.064
   3.805   (   6.670    6.108    6.354)   11.053
   3.887   (   5.125    3.984   -1.898)    6.763
   5.182   (  -1.627   -2.875    2.986)    4.453
   5.214   (  -1.886   -3.582   -0.603)    4.093
   5.359   (   0.632   -1.632    0.790)    1.920
   5.381   (   1.024   -1.716   -0.887)    2.186
   5.465   (  -1.374   -0.765    0.241)    1.591
   5.478   (  -1.021    0.410   -1.056)    1.525
======================= Grid point 588 (83/84) =======================
q-point: ( 0.38  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.406   (  -0.134    0.232   -2.622)    2.636
   1.473   (   0.847   -1.467   -8.778)    8.939
   2.097   (  -4.281    7.415    2.663)    8.967
   2.134   (  -4.084    7.074   -0.903)    8.219
   3.908   (  -0.164    0.284    4.785)    4.797
   3.954   (   0.477   -0.826    0.096)    0.959
   5.156   (   1.144   -1.982    1.918)    2.986
   5.167   (   1.059   -1.835    0.258)    2.135
   5.357   (   1.058   -1.832    0.755)    2.246
   5.381   (   1.648   -2.855   -1.086)    3.471
   5.450   (  -0.910    1.575    0.771)    1.976
   5.471   (  -0.845    1.464   -1.267)    2.113
======================= Grid point 601 (84/84) =======================
q-point: ( 0.31  0.31  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 5.98e-04 5.98e-04 2.84e-04 2.84e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.412   (   0.705    1.221   -3.982)    4.224
   1.456   (  -1.194   -2.068   -8.027)    8.375
   2.176   (   1.011    1.751    2.471)    3.193
   2.208   (   0.948    1.642   -0.743)    2.036
   3.911   (   2.056    3.561    4.110)    5.813
   3.944   (   0.614    1.063    0.782)    1.456
   5.131   (  -0.722   -1.250    2.512)    2.897
   5.147   (  -1.354   -2.344    0.736)    2.806
   5.340   (  -0.144   -0.249   -0.359)    0.460
   5.349   (  -0.517   -0.895   -1.014)    1.448
   5.468   (   0.154    0.266    0.724)    0.787
   5.487   (   0.130    0.226   -1.118)    1.148
=================== 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/14196
   10.0    540.060    540.060    482.637      0.000      0.000      0.000 3/14196
   20.0    257.106    257.106    271.321     -0.000      0.000      0.000 3/14196
   30.0    199.960    199.960    205.089     -0.000      0.000      0.000 3/14196
   40.0    159.812    159.812    162.154     -0.000      0.000      0.000 3/14196
   50.0    129.156    129.156    130.940     -0.000      0.000      0.000 3/14196
   60.0    106.593    106.593    108.360     -0.000      0.000      0.000 3/14196
   70.0     90.081     90.081     91.903     -0.000      0.000      0.000 3/14196
   80.0     77.790     77.790     79.637     -0.000      0.000      0.000 3/14196
   90.0     68.405     68.405     70.241     -0.000      0.000      0.000 3/14196
  100.0     61.053     61.053     62.850     -0.000      0.000      0.000 3/14196
  110.0     55.156     55.156     56.897     -0.000      0.000      0.000 3/14196
  120.0     50.326     50.326     52.003     -0.000      0.000      0.000 3/14196
  130.0     46.298     46.298     47.908     -0.000      0.000      0.000 3/14196
  140.0     42.889     42.889     44.431     -0.000      0.000      0.000 3/14196
  150.0     39.964     39.964     41.440     -0.000      0.000      0.000 3/14196
  160.0     37.425     37.425     38.838     -0.000      0.000      0.000 3/14196
  170.0     35.201     35.201     36.553     -0.000      0.000      0.000 3/14196
  180.0     33.234     33.234     34.530     -0.000      0.000      0.000 3/14196
  190.0     31.483     31.483     32.725     -0.000      0.000      0.000 3/14196
  200.0     29.912     29.912     31.104     -0.000      0.000      0.000 3/14196
  210.0     28.495     28.495     29.640     -0.000      0.000      0.000 3/14196
  220.0     27.209     27.209     28.311     -0.000      0.000      0.000 3/14196
  230.0     26.038     26.038     27.098     -0.000      0.000      0.000 3/14196
  240.0     24.966     24.966     25.988     -0.000      0.000      0.000 3/14196
  250.0     23.980     23.980     24.966     -0.000      0.000      0.000 3/14196
  260.0     23.071     23.071     24.023     -0.000      0.000      0.000 3/14196
  270.0     22.230     22.230     23.150     -0.000      0.000      0.000 3/14196
  280.0     21.449     21.449     22.340     -0.000      0.000      0.000 3/14196
  290.0     20.723     20.723     21.585     -0.000      0.000      0.000 3/14196
  300.0     20.044     20.044     20.880     -0.000      0.000      0.000 3/14196
  310.0     19.410     19.410     20.220     -0.000      0.000      0.000 3/14196
  320.0     18.815     18.815     19.601     -0.000      0.000      0.000 3/14196
  330.0     18.256     18.256     19.020     -0.000      0.000      0.000 3/14196
  340.0     17.730     17.730     18.472     -0.000      0.000      0.000 3/14196
  350.0     17.233     17.233     17.956     -0.000      0.000      0.000 3/14196
  360.0     16.764     16.764     17.468     -0.000      0.000      0.000 3/14196
  370.0     16.321     16.321     17.006     -0.000      0.000      0.000 3/14196
  380.0     15.900     15.900     16.568     -0.000      0.000      0.000 3/14196
  390.0     15.501     15.501     16.152     -0.000      0.000      0.000 3/14196
  400.0     15.121     15.121     15.757     -0.000      0.000      0.000 3/14196
  410.0     14.760     14.760     15.381     -0.000      0.000      0.000 3/14196
  420.0     14.416     14.416     15.022     -0.000      0.000      0.000 3/14196
  430.0     14.088     14.088     14.680     -0.000      0.000      0.000 3/14196
  440.0     13.774     13.774     14.353     -0.000      0.000      0.000 3/14196
  450.0     13.474     13.474     14.041     -0.000      0.000      0.000 3/14196
  460.0     13.188     13.188     13.742     -0.000      0.000      0.000 3/14196
  470.0     12.913     12.913     13.456     -0.000      0.000      0.000 3/14196
  480.0     12.649     12.649     13.181     -0.000      0.000      0.000 3/14196
  490.0     12.396     12.396     12.918     -0.000      0.000      0.000 3/14196
  500.0     12.153     12.153     12.664     -0.000      0.000      0.000 3/14196
  510.0     11.920     11.920     12.421     -0.000      0.000      0.000 3/14196
  520.0     11.695     11.695     12.187     -0.000      0.000      0.000 3/14196
  530.0     11.479     11.479     11.961     -0.000      0.000      0.000 3/14196
  540.0     11.271     11.271     11.744     -0.000      0.000      0.000 3/14196
  550.0     11.070     11.070     11.535     -0.000      0.000      0.000 3/14196
  560.0     10.876     10.876     11.333     -0.000      0.000      0.000 3/14196
  570.0     10.689     10.689     11.138     -0.000      0.000      0.000 3/14196
  580.0     10.508     10.508     10.949     -0.000      0.000      0.000 3/14196
  590.0     10.334     10.334     10.767     -0.000      0.000      0.000 3/14196
  600.0     10.165     10.165     10.591     -0.000      0.000      0.000 3/14196
  610.0     10.001     10.001     10.421     -0.000      0.000      0.000 3/14196
  620.0      9.843      9.843     10.256     -0.000      0.000      0.000 3/14196
  630.0      9.690      9.690     10.096     -0.000      0.000      0.000 3/14196
  640.0      9.541      9.541      9.941     -0.000      0.000      0.000 3/14196
  650.0      9.397      9.397      9.791     -0.000      0.000      0.000 3/14196
  660.0      9.257      9.257      9.645     -0.000      0.000      0.000 3/14196
  670.0      9.121      9.121      9.503     -0.000      0.000      0.000 3/14196
  680.0      8.990      8.990      9.366     -0.000      0.000      0.000 3/14196
  690.0      8.862      8.862      9.233     -0.000      0.000      0.000 3/14196
  700.0      8.737      8.737      9.103     -0.000      0.000      0.000 3/14196
  710.0      8.617      8.617      8.977     -0.000      0.000      0.000 3/14196
  720.0      8.499      8.499      8.854     -0.000      0.000      0.000 3/14196
  730.0      8.385      8.385      8.735     -0.000      0.000      0.000 3/14196
  740.0      8.273      8.273      8.619     -0.000      0.000      0.000 3/14196
  750.0      8.165      8.165      8.506     -0.000      0.000      0.000 3/14196
  760.0      8.059      8.059      8.396     -0.000      0.000      0.000 3/14196
  770.0      7.956      7.956      8.289     -0.000      0.000      0.000 3/14196
  780.0      7.856      7.856      8.184     -0.000      0.000      0.000 3/14196
  790.0      7.758      7.758      8.082     -0.000      0.000      0.000 3/14196
  800.0      7.663      7.663      7.983     -0.000      0.000      0.000 3/14196
  810.0      7.570      7.570      7.885     -0.000      0.000      0.000 3/14196
  820.0      7.479      7.479      7.791     -0.000      0.000      0.000 3/14196
  830.0      7.390      7.390      7.698     -0.000      0.000      0.000 3/14196
  840.0      7.304      7.304      7.608     -0.000      0.000      0.000 3/14196
  850.0      7.219      7.219      7.520     -0.000      0.000      0.000 3/14196
  860.0      7.136      7.136      7.434     -0.000      0.000      0.000 3/14196
  870.0      7.056      7.056      7.350     -0.000      0.000      0.000 3/14196
  880.0      6.977      6.977      7.267     -0.000      0.000      0.000 3/14196
  890.0      6.899      6.899      7.187     -0.000      0.000      0.000 3/14196
  900.0      6.824      6.824      7.108     -0.000      0.000      0.000 3/14196
  910.0      6.750      6.750      7.031     -0.000      0.000      0.000 3/14196
  920.0      6.678      6.678      6.956     -0.000      0.000      0.000 3/14196
  930.0      6.607      6.607      6.882     -0.000      0.000      0.000 3/14196
  940.0      6.538      6.538      6.810     -0.000      0.000      0.000 3/14196
  950.0      6.470      6.470      6.739     -0.000      0.000      0.000 3/14196
  960.0      6.404      6.404      6.670     -0.000      0.000      0.000 3/14196
  970.0      6.339      6.339      6.602     -0.000      0.000      0.000 3/14196
  980.0      6.275      6.275      6.536     -0.000      0.000      0.000 3/14196
  990.0      6.212      6.212      6.471     -0.000      0.000      0.000 3/14196
 1000.0      6.151      6.151      6.407     -0.000      0.000      0.000 3/14196

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

