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

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

-------------------------[time 2026-01-08 00:09:48]-------------------------
Compiled with OpenMP support (max 48 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.580347819875001    0.000000000000000    0.000000000000000
  b   -2.290173909937501    3.966697570180421    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.502629120000000
Atomic positions (fractional):
   *1 Cd  0.33333333333333  0.66666666666667  0.00026922740646 112.411
    2 Cd  0.66666666666667  0.33333333333333  0.50026922740646 112.411
   *3 Te  0.33333333333333  0.66666666666667  0.37473077259354 127.600
    4 Te  0.66666666666667  0.33333333333333  0.87473077259354 127.600
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.580347819875001    0.000000000000000    0.000000000000000
  b   -2.290173909937501    3.966697570180421    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.502629120000000
Atomic positions (fractional):
   *1 Cd  0.33333333333333  0.66666666666667  0.00026922740646 112.411 > 1
    2 Cd  0.66666666666667  0.33333333333333  0.50026922740646 112.411 > 2
   *3 Te  0.33333333333333  0.66666666666667  0.37473077259354 127.600 > 3
    4 Te  0.66666666666667  0.33333333333333  0.87473077259354 127.600 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   18.321391279500006    0.000000000000000    0.000000000000000
  b   -9.160695639750005   15.866790280721684    0.000000000000000
  c    0.000000000000000    0.000000000000000   15.005258240000000
Atomic positions (fractional):
   *1 Cd  0.08333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    2 Cd  0.33333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    3 Cd  0.58333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    4 Cd  0.83333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    5 Cd  0.08333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    6 Cd  0.33333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    7 Cd  0.58333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    8 Cd  0.83333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    9 Cd  0.08333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   10 Cd  0.33333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   11 Cd  0.58333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   12 Cd  0.83333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   13 Cd  0.08333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   14 Cd  0.33333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   15 Cd  0.58333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   16 Cd  0.83333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   17 Cd  0.08333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   18 Cd  0.33333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   19 Cd  0.58333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   20 Cd  0.83333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   21 Cd  0.08333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   22 Cd  0.33333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   23 Cd  0.58333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   24 Cd  0.83333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   25 Cd  0.08333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   26 Cd  0.33333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   27 Cd  0.58333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   28 Cd  0.83333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   29 Cd  0.08333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   30 Cd  0.33333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   31 Cd  0.58333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   32 Cd  0.83333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   33 Cd  0.16666666666667  0.08333333333333  0.25013461370323 112.411 > 2
   34 Cd  0.41666666666667  0.08333333333333  0.25013461370323 112.411 > 2
   35 Cd  0.66666666666667  0.08333333333333  0.25013461370323 112.411 > 2
   36 Cd  0.91666666666667  0.08333333333333  0.25013461370323 112.411 > 2
   37 Cd  0.16666666666667  0.33333333333333  0.25013461370323 112.411 > 2
   38 Cd  0.41666666666667  0.33333333333333  0.25013461370323 112.411 > 2
   39 Cd  0.66666666666667  0.33333333333333  0.25013461370323 112.411 > 2
   40 Cd  0.91666666666667  0.33333333333333  0.25013461370323 112.411 > 2
   41 Cd  0.16666666666667  0.58333333333333  0.25013461370323 112.411 > 2
   42 Cd  0.41666666666667  0.58333333333333  0.25013461370323 112.411 > 2
   43 Cd  0.66666666666667  0.58333333333333  0.25013461370323 112.411 > 2
   44 Cd  0.91666666666667  0.58333333333333  0.25013461370323 112.411 > 2
   45 Cd  0.16666666666667  0.83333333333333  0.25013461370323 112.411 > 2
   46 Cd  0.41666666666667  0.83333333333333  0.25013461370323 112.411 > 2
   47 Cd  0.66666666666667  0.83333333333333  0.25013461370323 112.411 > 2
   48 Cd  0.91666666666667  0.83333333333333  0.25013461370323 112.411 > 2
   49 Cd  0.16666666666667  0.08333333333333  0.75013461370323 112.411 > 2
   50 Cd  0.41666666666667  0.08333333333333  0.75013461370323 112.411 > 2
   51 Cd  0.66666666666667  0.08333333333333  0.75013461370323 112.411 > 2
   52 Cd  0.91666666666667  0.08333333333333  0.75013461370323 112.411 > 2
   53 Cd  0.16666666666667  0.33333333333333  0.75013461370323 112.411 > 2
   54 Cd  0.41666666666667  0.33333333333333  0.75013461370323 112.411 > 2
   55 Cd  0.66666666666667  0.33333333333333  0.75013461370323 112.411 > 2
   56 Cd  0.91666666666667  0.33333333333333  0.75013461370323 112.411 > 2
   57 Cd  0.16666666666667  0.58333333333333  0.75013461370323 112.411 > 2
   58 Cd  0.41666666666667  0.58333333333333  0.75013461370323 112.411 > 2
   59 Cd  0.66666666666667  0.58333333333333  0.75013461370323 112.411 > 2
   60 Cd  0.91666666666667  0.58333333333333  0.75013461370323 112.411 > 2
   61 Cd  0.16666666666667  0.83333333333333  0.75013461370323 112.411 > 2
   62 Cd  0.41666666666667  0.83333333333333  0.75013461370323 112.411 > 2
   63 Cd  0.66666666666667  0.83333333333333  0.75013461370323 112.411 > 2
   64 Cd  0.91666666666667  0.83333333333333  0.75013461370323 112.411 > 2
  *65 Te  0.08333333333333  0.16666666666667  0.18736538629677 127.600 > 3
   66 Te  0.33333333333333  0.16666666666667  0.18736538629677 127.600 > 3
   67 Te  0.58333333333333  0.16666666666667  0.18736538629677 127.600 > 3
   68 Te  0.83333333333333  0.16666666666667  0.18736538629677 127.600 > 3
   69 Te  0.08333333333333  0.41666666666667  0.18736538629677 127.600 > 3
   70 Te  0.33333333333333  0.41666666666667  0.18736538629677 127.600 > 3
   71 Te  0.58333333333333  0.41666666666667  0.18736538629677 127.600 > 3
   72 Te  0.83333333333333  0.41666666666667  0.18736538629677 127.600 > 3
   73 Te  0.08333333333333  0.66666666666667  0.18736538629677 127.600 > 3
   74 Te  0.33333333333333  0.66666666666667  0.18736538629677 127.600 > 3
   75 Te  0.58333333333333  0.66666666666667  0.18736538629677 127.600 > 3
   76 Te  0.83333333333333  0.66666666666667  0.18736538629677 127.600 > 3
   77 Te  0.08333333333333  0.91666666666667  0.18736538629677 127.600 > 3
   78 Te  0.33333333333333  0.91666666666667  0.18736538629677 127.600 > 3
   79 Te  0.58333333333333  0.91666666666667  0.18736538629677 127.600 > 3
   80 Te  0.83333333333333  0.91666666666667  0.18736538629677 127.600 > 3
   81 Te  0.08333333333333  0.16666666666667  0.68736538629677 127.600 > 3
   82 Te  0.33333333333333  0.16666666666667  0.68736538629677 127.600 > 3
   83 Te  0.58333333333333  0.16666666666667  0.68736538629677 127.600 > 3
   84 Te  0.83333333333333  0.16666666666667  0.68736538629677 127.600 > 3
   85 Te  0.08333333333333  0.41666666666667  0.68736538629677 127.600 > 3
   86 Te  0.33333333333333  0.41666666666667  0.68736538629677 127.600 > 3
   87 Te  0.58333333333333  0.41666666666667  0.68736538629677 127.600 > 3
   88 Te  0.83333333333333  0.41666666666667  0.68736538629677 127.600 > 3
   89 Te  0.08333333333333  0.66666666666667  0.68736538629677 127.600 > 3
   90 Te  0.33333333333333  0.66666666666667  0.68736538629677 127.600 > 3
   91 Te  0.58333333333333  0.66666666666667  0.68736538629677 127.600 > 3
   92 Te  0.83333333333333  0.66666666666667  0.68736538629677 127.600 > 3
   93 Te  0.08333333333333  0.91666666666667  0.68736538629677 127.600 > 3
   94 Te  0.33333333333333  0.91666666666667  0.68736538629677 127.600 > 3
   95 Te  0.58333333333333  0.91666666666667  0.68736538629677 127.600 > 3
   96 Te  0.83333333333333  0.91666666666667  0.68736538629677 127.600 > 3
   97 Te  0.16666666666667  0.08333333333333  0.43736538629677 127.600 > 4
   98 Te  0.41666666666667  0.08333333333333  0.43736538629677 127.600 > 4
   99 Te  0.66666666666667  0.08333333333333  0.43736538629677 127.600 > 4
  100 Te  0.91666666666667  0.08333333333333  0.43736538629677 127.600 > 4
  101 Te  0.16666666666667  0.33333333333333  0.43736538629677 127.600 > 4
  102 Te  0.41666666666667  0.33333333333333  0.43736538629677 127.600 > 4
  103 Te  0.66666666666667  0.33333333333333  0.43736538629677 127.600 > 4
  104 Te  0.91666666666667  0.33333333333333  0.43736538629677 127.600 > 4
  105 Te  0.16666666666667  0.58333333333333  0.43736538629677 127.600 > 4
  106 Te  0.41666666666667  0.58333333333333  0.43736538629677 127.600 > 4
  107 Te  0.66666666666667  0.58333333333333  0.43736538629677 127.600 > 4
  108 Te  0.91666666666667  0.58333333333333  0.43736538629677 127.600 > 4
  109 Te  0.16666666666667  0.83333333333333  0.43736538629677 127.600 > 4
  110 Te  0.41666666666667  0.83333333333333  0.43736538629677 127.600 > 4
  111 Te  0.66666666666667  0.83333333333333  0.43736538629677 127.600 > 4
  112 Te  0.91666666666667  0.83333333333333  0.43736538629677 127.600 > 4
  113 Te  0.16666666666667  0.08333333333333  0.93736538629677 127.600 > 4
  114 Te  0.41666666666667  0.08333333333333  0.93736538629677 127.600 > 4
  115 Te  0.66666666666667  0.08333333333333  0.93736538629677 127.600 > 4
  116 Te  0.91666666666667  0.08333333333333  0.93736538629677 127.600 > 4
  117 Te  0.16666666666667  0.33333333333333  0.93736538629677 127.600 > 4
  118 Te  0.41666666666667  0.33333333333333  0.93736538629677 127.600 > 4
  119 Te  0.66666666666667  0.33333333333333  0.93736538629677 127.600 > 4
  120 Te  0.91666666666667  0.33333333333333  0.93736538629677 127.600 > 4
  121 Te  0.16666666666667  0.58333333333333  0.93736538629677 127.600 > 4
  122 Te  0.41666666666667  0.58333333333333  0.93736538629677 127.600 > 4
  123 Te  0.66666666666667  0.58333333333333  0.93736538629677 127.600 > 4
  124 Te  0.91666666666667  0.58333333333333  0.93736538629677 127.600 > 4
  125 Te  0.16666666666667  0.83333333333333  0.93736538629677 127.600 > 4
  126 Te  0.41666666666667  0.83333333333333  0.93736538629677 127.600 > 4
  127 Te  0.66666666666667  0.83333333333333  0.93736538629677 127.600 > 4
  128 Te  0.91666666666667  0.83333333333333  0.93736538629677 127.600 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           17.8795723   -0.0000000    0.0000000
            0.0000000   17.8795723    0.0000000
            0.0000000    0.0000000   15.9093710
-------------------------- Born effective charges --------------------------
    1 Cd    2.4657830   -0.0000000    0.0000000
            0.0000000    2.4657830    0.0000000
            0.0000000    0.0000000    2.4850593
    2 Cd    2.4657830   -0.0000000    0.0000000
            0.0000000    2.4657830    0.0000000
            0.0000000    0.0000000    2.4850593
    3 Te   -2.4657830    0.0000000    0.0000000
            0.0000000   -2.4657830    0.0000000
            0.0000000    0.0000000   -2.4850593
    4 Te   -2.4657830    0.0000000    0.0000000
            0.0000000   -2.4657830    0.0000000
            0.0000000    0.0000000   -2.4850593
----------------------------------------------------------------------------
Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
-------------------------------- Symfc start -------------------------------
Symfc version 1.5.4 (https://github.com/symfc/symfc)
Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024)
Computing [2] order force constants.
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 188
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
-----
Solver_atoms: 1 -- 128 / 128
Time (Solver_compr_matrix_reshape): 0.008
Solver_block: 80 / 80
 - Time: 0.121
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.130
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 188
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 00:09:52]-------------------------
Compiled with OpenMP support (max 48 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.580347819875001    0.000000000000000    0.000000000000000
  b   -2.290173909937501    3.966697570180421    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.502629120000000
Atomic positions (fractional):
    1 Cd  0.33333333333333  0.66666666666667  0.00026922740646 112.411
    2 Cd  0.66666666666667  0.33333333333333  0.50026922740646 112.411
    3 Te  0.33333333333333  0.66666666666667  0.37473077259354 127.600
    4 Te  0.66666666666667  0.33333333333333  0.87473077259354 127.600
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   18.321391279500006    0.000000000000000    0.000000000000000
  b   -9.160695639750005   15.866790280721684    0.000000000000000
  c    0.000000000000000    0.000000000000000   15.005258240000000
Atomic positions (fractional):
    1 Cd  0.08333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    2 Cd  0.33333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    3 Cd  0.58333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    4 Cd  0.83333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    5 Cd  0.08333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    6 Cd  0.33333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    7 Cd  0.58333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    8 Cd  0.83333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    9 Cd  0.08333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   10 Cd  0.33333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   11 Cd  0.58333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   12 Cd  0.83333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   13 Cd  0.08333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   14 Cd  0.33333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   15 Cd  0.58333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   16 Cd  0.83333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   17 Cd  0.08333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   18 Cd  0.33333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   19 Cd  0.58333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   20 Cd  0.83333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   21 Cd  0.08333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   22 Cd  0.33333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   23 Cd  0.58333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   24 Cd  0.83333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   25 Cd  0.08333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   26 Cd  0.33333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   27 Cd  0.58333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   28 Cd  0.83333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   29 Cd  0.08333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   30 Cd  0.33333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   31 Cd  0.58333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   32 Cd  0.83333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   33 Cd  0.16666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   34 Cd  0.41666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   35 Cd  0.66666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   36 Cd  0.91666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   37 Cd  0.16666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   38 Cd  0.41666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   39 Cd  0.66666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   40 Cd  0.91666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   41 Cd  0.16666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   42 Cd  0.41666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   43 Cd  0.66666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   44 Cd  0.91666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   45 Cd  0.16666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   46 Cd  0.41666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   47 Cd  0.66666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   48 Cd  0.91666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   49 Cd  0.16666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   50 Cd  0.41666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   51 Cd  0.66666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   52 Cd  0.91666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   53 Cd  0.16666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   54 Cd  0.41666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   55 Cd  0.66666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   56 Cd  0.91666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   57 Cd  0.16666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   58 Cd  0.41666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   59 Cd  0.66666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   60 Cd  0.91666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   61 Cd  0.16666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   62 Cd  0.41666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   63 Cd  0.66666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   64 Cd  0.91666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   65 Te  0.08333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   66 Te  0.33333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   67 Te  0.58333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   68 Te  0.83333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   69 Te  0.08333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   70 Te  0.33333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   71 Te  0.58333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   72 Te  0.83333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   73 Te  0.08333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   74 Te  0.33333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   75 Te  0.58333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   76 Te  0.83333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   77 Te  0.08333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   78 Te  0.33333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   79 Te  0.58333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   80 Te  0.83333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   81 Te  0.08333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   82 Te  0.33333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   83 Te  0.58333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   84 Te  0.83333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   85 Te  0.08333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   86 Te  0.33333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   87 Te  0.58333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   88 Te  0.83333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   89 Te  0.08333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   90 Te  0.33333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   91 Te  0.58333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   92 Te  0.83333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   93 Te  0.08333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   94 Te  0.33333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   95 Te  0.58333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   96 Te  0.83333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   97 Te  0.16666666666667  0.08333333333333  0.43736538629677 127.600 > 97
   98 Te  0.41666666666667  0.08333333333333  0.43736538629677 127.600 > 97
   99 Te  0.66666666666667  0.08333333333333  0.43736538629677 127.600 > 97
  100 Te  0.91666666666667  0.08333333333333  0.43736538629677 127.600 > 97
  101 Te  0.16666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  102 Te  0.41666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  103 Te  0.66666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  104 Te  0.91666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  105 Te  0.16666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  106 Te  0.41666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  107 Te  0.66666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  108 Te  0.91666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  109 Te  0.16666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  110 Te  0.41666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  111 Te  0.66666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  112 Te  0.91666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  113 Te  0.16666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  114 Te  0.41666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  115 Te  0.66666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  116 Te  0.91666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  117 Te  0.16666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  118 Te  0.41666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  119 Te  0.66666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  120 Te  0.91666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  121 Te  0.16666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  122 Te  0.41666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  123 Te  0.66666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  124 Te  0.91666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  125 Te  0.16666666666667  0.83333333333333  0.93736538629677 127.600 > 97
  126 Te  0.41666666666667  0.83333333333333  0.93736538629677 127.600 > 97
  127 Te  0.66666666666667  0.83333333333333  0.93736538629677 127.600 > 97
  128 Te  0.91666666666667  0.83333333333333  0.93736538629677 127.600 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           17.8795723   -0.0000000    0.0000000
            0.0000000   17.8795723    0.0000000
            0.0000000    0.0000000   15.9093710
-------------------------- Born effective charges --------------------------
    1 Cd    2.4657830   -0.0000000    0.0000000
            0.0000000    2.4657830    0.0000000
            0.0000000    0.0000000    2.4850593
    2 Cd    2.4657830   -0.0000000    0.0000000
            0.0000000    2.4657830    0.0000000
            0.0000000    0.0000000    2.4850593
    3 Te   -2.4657830    0.0000000    0.0000000
            0.0000000   -2.4657830    0.0000000
            0.0000000    0.0000000   -2.4850593
    4 Te   -2.4657830    0.0000000    0.0000000
            0.0000000   -2.4657830    0.0000000
            0.0000000    0.0000000   -2.4850593
----------------------------------------------------------------------------
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.00000135 (zzz) -0.00000135 (zzz) -0.00000135 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 00:09:56]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 00:09:57]-------------------------
Compiled with OpenMP support (max 48 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.580347819875001    0.000000000000000    0.000000000000000
  b   -2.290173909937501    3.966697570180421    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.502629120000000
Atomic positions (fractional):
    1 Cd  0.33333333333333  0.66666666666667  0.00026922740646 112.411
    2 Cd  0.66666666666667  0.33333333333333  0.50026922740646 112.411
    3 Te  0.33333333333333  0.66666666666667  0.37473077259354 127.600
    4 Te  0.66666666666667  0.33333333333333  0.87473077259354 127.600
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   18.321391279500006    0.000000000000000    0.000000000000000
  b   -9.160695639750005   15.866790280721684    0.000000000000000
  c    0.000000000000000    0.000000000000000   15.005258240000000
Atomic positions (fractional):
    1 Cd  0.08333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    2 Cd  0.33333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    3 Cd  0.58333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    4 Cd  0.83333333333333  0.16666666666667  0.00013461370323 112.411 > 1
    5 Cd  0.08333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    6 Cd  0.33333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    7 Cd  0.58333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    8 Cd  0.83333333333333  0.41666666666667  0.00013461370323 112.411 > 1
    9 Cd  0.08333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   10 Cd  0.33333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   11 Cd  0.58333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   12 Cd  0.83333333333333  0.66666666666667  0.00013461370323 112.411 > 1
   13 Cd  0.08333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   14 Cd  0.33333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   15 Cd  0.58333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   16 Cd  0.83333333333333  0.91666666666667  0.00013461370323 112.411 > 1
   17 Cd  0.08333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   18 Cd  0.33333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   19 Cd  0.58333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   20 Cd  0.83333333333333  0.16666666666667  0.50013461370323 112.411 > 1
   21 Cd  0.08333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   22 Cd  0.33333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   23 Cd  0.58333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   24 Cd  0.83333333333333  0.41666666666667  0.50013461370323 112.411 > 1
   25 Cd  0.08333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   26 Cd  0.33333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   27 Cd  0.58333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   28 Cd  0.83333333333333  0.66666666666667  0.50013461370323 112.411 > 1
   29 Cd  0.08333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   30 Cd  0.33333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   31 Cd  0.58333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   32 Cd  0.83333333333333  0.91666666666667  0.50013461370323 112.411 > 1
   33 Cd  0.16666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   34 Cd  0.41666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   35 Cd  0.66666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   36 Cd  0.91666666666667  0.08333333333333  0.25013461370323 112.411 > 33
   37 Cd  0.16666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   38 Cd  0.41666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   39 Cd  0.66666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   40 Cd  0.91666666666667  0.33333333333333  0.25013461370323 112.411 > 33
   41 Cd  0.16666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   42 Cd  0.41666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   43 Cd  0.66666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   44 Cd  0.91666666666667  0.58333333333333  0.25013461370323 112.411 > 33
   45 Cd  0.16666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   46 Cd  0.41666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   47 Cd  0.66666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   48 Cd  0.91666666666667  0.83333333333333  0.25013461370323 112.411 > 33
   49 Cd  0.16666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   50 Cd  0.41666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   51 Cd  0.66666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   52 Cd  0.91666666666667  0.08333333333333  0.75013461370323 112.411 > 33
   53 Cd  0.16666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   54 Cd  0.41666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   55 Cd  0.66666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   56 Cd  0.91666666666667  0.33333333333333  0.75013461370323 112.411 > 33
   57 Cd  0.16666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   58 Cd  0.41666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   59 Cd  0.66666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   60 Cd  0.91666666666667  0.58333333333333  0.75013461370323 112.411 > 33
   61 Cd  0.16666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   62 Cd  0.41666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   63 Cd  0.66666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   64 Cd  0.91666666666667  0.83333333333333  0.75013461370323 112.411 > 33
   65 Te  0.08333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   66 Te  0.33333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   67 Te  0.58333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   68 Te  0.83333333333333  0.16666666666667  0.18736538629677 127.600 > 65
   69 Te  0.08333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   70 Te  0.33333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   71 Te  0.58333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   72 Te  0.83333333333333  0.41666666666667  0.18736538629677 127.600 > 65
   73 Te  0.08333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   74 Te  0.33333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   75 Te  0.58333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   76 Te  0.83333333333333  0.66666666666667  0.18736538629677 127.600 > 65
   77 Te  0.08333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   78 Te  0.33333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   79 Te  0.58333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   80 Te  0.83333333333333  0.91666666666667  0.18736538629677 127.600 > 65
   81 Te  0.08333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   82 Te  0.33333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   83 Te  0.58333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   84 Te  0.83333333333333  0.16666666666667  0.68736538629677 127.600 > 65
   85 Te  0.08333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   86 Te  0.33333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   87 Te  0.58333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   88 Te  0.83333333333333  0.41666666666667  0.68736538629677 127.600 > 65
   89 Te  0.08333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   90 Te  0.33333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   91 Te  0.58333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   92 Te  0.83333333333333  0.66666666666667  0.68736538629677 127.600 > 65
   93 Te  0.08333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   94 Te  0.33333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   95 Te  0.58333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   96 Te  0.83333333333333  0.91666666666667  0.68736538629677 127.600 > 65
   97 Te  0.16666666666667  0.08333333333333  0.43736538629677 127.600 > 97
   98 Te  0.41666666666667  0.08333333333333  0.43736538629677 127.600 > 97
   99 Te  0.66666666666667  0.08333333333333  0.43736538629677 127.600 > 97
  100 Te  0.91666666666667  0.08333333333333  0.43736538629677 127.600 > 97
  101 Te  0.16666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  102 Te  0.41666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  103 Te  0.66666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  104 Te  0.91666666666667  0.33333333333333  0.43736538629677 127.600 > 97
  105 Te  0.16666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  106 Te  0.41666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  107 Te  0.66666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  108 Te  0.91666666666667  0.58333333333333  0.43736538629677 127.600 > 97
  109 Te  0.16666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  110 Te  0.41666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  111 Te  0.66666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  112 Te  0.91666666666667  0.83333333333333  0.43736538629677 127.600 > 97
  113 Te  0.16666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  114 Te  0.41666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  115 Te  0.66666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  116 Te  0.91666666666667  0.08333333333333  0.93736538629677 127.600 > 97
  117 Te  0.16666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  118 Te  0.41666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  119 Te  0.66666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  120 Te  0.91666666666667  0.33333333333333  0.93736538629677 127.600 > 97
  121 Te  0.16666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  122 Te  0.41666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  123 Te  0.66666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  124 Te  0.91666666666667  0.58333333333333  0.93736538629677 127.600 > 97
  125 Te  0.16666666666667  0.83333333333333  0.93736538629677 127.600 > 97
  126 Te  0.41666666666667  0.83333333333333  0.93736538629677 127.600 > 97
  127 Te  0.66666666666667  0.83333333333333  0.93736538629677 127.600 > 97
  128 Te  0.91666666666667  0.83333333333333  0.93736538629677 127.600 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           17.8795723   -0.0000000    0.0000000
            0.0000000   17.8795723    0.0000000
            0.0000000    0.0000000   15.9093710
-------------------------- Born effective charges --------------------------
    1 Cd    2.4657830   -0.0000000    0.0000000
            0.0000000    2.4657830    0.0000000
            0.0000000    0.0000000    2.4850593
    2 Cd    2.4657830   -0.0000000    0.0000000
            0.0000000    2.4657830    0.0000000
            0.0000000    0.0000000    2.4850593
    3 Te   -2.4657830    0.0000000    0.0000000
            0.0000000   -2.4657830    0.0000000
            0.0000000    0.0000000   -2.4850593
    4 Te   -2.4657830    0.0000000    0.0000000
            0.0000000   -2.4657830    0.0000000
            0.0000000    0.0000000   -2.4850593
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000135 (zzz) -0.00000135 (zzz) -0.00000135 (zzz)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 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.81, Number of G-points: 313, Lambda: 0.35
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: 2.72e-04 2.72e-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
   0.779   (   0.000   -0.000   -0.000)    0.000
   0.779   (  -0.000    0.000    0.000)    0.000
   3.189   (   0.000    0.000   -0.000)    0.000
   4.282   (   0.000   -0.000   -0.000)    0.000
   4.300   (  -0.000   -0.000    0.000)    0.000
   4.300   (   0.000    0.000   -0.000)    0.000
   4.336   (   0.000    0.000    0.000)    0.000
   4.336   (   0.000   -0.000   -0.000)    0.000
   4.379   (   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: 2.72e-04 2.72e-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.246   (  10.681    6.166    0.000)   12.333
   0.254   (  11.562    6.675    0.000)   13.351
   0.633   (  27.447   15.847    0.000)   31.693
   0.804   (   2.283    1.318    0.000)    2.636
   0.968   (  14.193    8.194    0.000)   16.388
   3.153   (  -2.801   -1.617    0.000)    3.235
   4.272   (  -2.227   -1.286    0.000)    2.572
   4.282   (  -0.110   -0.063    0.000)    0.127
   4.297   (  -0.061   -0.035    0.000)    0.070
   4.346   (   0.628    0.363    0.000)    0.725
   4.424   (   3.615    2.087    0.000)    4.174
   4.763   (   1.708    0.986    0.000)    1.972
======================= Grid point 2 (3/84) =======================
q-point: ( 0.15  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.470   (   9.180    5.300    0.000)   10.601
   0.509   (  11.152    6.439    0.000)   12.878
   0.878   (   4.100    2.367    0.000)    4.735
   1.195   (  19.838   11.453    0.000)   22.907
   1.313   (  17.295    9.985    0.000)   19.971
   3.099   (  -1.088   -0.628    0.000)    1.256
   4.201   (  -3.970   -2.292    0.000)    4.584
   4.283   (   0.238    0.138    0.000)    0.275
   4.303   (   0.674    0.389    0.000)    0.779
   4.355   (   0.156    0.090    0.000)    0.180
   4.518   (   4.279    2.470    0.000)    4.941
   4.805   (   1.646    0.950    0.000)    1.900
======================= Grid point 3 (4/84) =======================
q-point: ( 0.23  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.646   (   6.331    3.655    0.000)    7.311
   0.750   (  10.298    5.945    0.000)   11.891
   0.972   (   3.825    2.208    0.000)    4.417
   1.423   (   3.686    2.128    0.000)    4.256
   1.785   (  21.506   12.416    0.000)   24.833
   3.139   (   5.008    2.892    0.000)    5.783
   4.096   (  -5.444   -3.143    0.000)    6.286
   4.295   (   0.835    0.482    0.000)    0.965
   4.324   (   1.080    0.624    0.000)    1.247
   4.356   (   0.019    0.011    0.000)    0.022
   4.598   (   2.542    1.468    0.000)    2.935
   4.820   (  -0.692   -0.400    0.000)    0.799
======================= Grid point 4 (5/84) =======================
q-point: ( 0.31  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.749   (   2.916    1.683    0.000)    3.367
   0.966   (   8.804    5.083    0.000)   10.166
   1.031   (   1.214    0.701    0.000)    1.402
   1.434   (  -1.545   -0.892    0.000)    1.784
   2.231   (  18.201   10.508    0.000)   21.017
   3.310   (   9.583    5.533    0.000)   11.066
   3.958   (  -6.866   -3.964    0.000)    7.928
   4.319   (   1.286    0.743    0.000)    1.485
   4.344   (   0.606    0.350    0.000)    0.700
   4.366   (   1.021    0.589    0.000)    1.179
   4.625   (  -0.266   -0.153    0.000)    0.307
   4.767   (  -3.853   -2.225    0.000)    4.450
======================= Grid point 5 (6/84) =======================
q-point: ( 0.38  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.786   (   0.712    0.411    0.000)    0.822
   1.024   (  -1.652   -0.954    0.000)    1.908
   1.139   (   6.576    3.797    0.000)    7.594
   1.392   (  -1.794   -1.036    0.000)    2.071
   2.596   (  14.210    8.204    0.000)   16.408
   3.537   (   9.848    5.686    0.000)   11.371
   3.796   (  -7.014   -4.049    0.000)    8.099
   4.342   (   0.231    0.133    0.000)    0.266
   4.359   (   0.759    0.438    0.000)    0.877
   4.391   (   1.029    0.594    0.000)    1.188
   4.575   (  -4.288   -2.476    0.000)    4.952
   4.678   (  -3.078   -1.777    0.000)    3.554
======================= Grid point 6 (7/84) =======================
q-point: ( 0.46  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.795   (   0.231    0.133    0.000)    0.266
   0.980   (  -1.540   -0.889    0.000)    1.778
   1.252   (   3.218    1.858    0.000)    3.716
   1.358   (  -1.269   -0.733    0.000)    1.465
   2.843   (   6.709    3.874    0.000)    7.747
   3.608   (  -3.416   -1.972    0.000)    3.945
   3.779   (   5.171    2.986    0.000)    5.971
   4.303   (  -3.081   -1.779    0.000)    3.558
   4.374   (   0.429    0.248    0.000)    0.495
   4.408   (   0.388    0.224    0.000)    0.448
   4.495   (  -1.495   -0.863    0.000)    1.726
   4.635   (  -0.944   -0.545    0.000)    1.090
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.441   (   6.463   11.194    0.000)   12.926
   0.459   (   7.236   12.532    0.000)   14.471
   0.869   (   3.810    6.598    0.000)    7.619
   1.080   (  13.828   23.951    0.000)   27.656
   1.172   (   7.156   12.395    0.000)   14.312
   3.108   (  -1.279   -2.215    0.000)    2.557
   4.222   (  -2.120   -3.672    0.000)    4.240
   4.280   (   0.023    0.040    0.000)    0.047
   4.318   (   0.615    1.065    0.000)    1.229
   4.346   (   0.496    0.859    0.000)    0.992
   4.495   (   2.456    4.253    0.000)    4.911
   4.787   (   0.704    1.219    0.000)    1.408
======================= Grid point 15 (9/84) =======================
q-point: ( 0.15  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.663   (   8.053    8.907    0.000)   12.008
   0.677   (   3.254   13.557    0.000)   13.943
   0.991   (   2.114    8.237    0.000)    8.504
   1.368   (   5.526    5.248    0.000)    7.621
   1.569   (  18.906   16.350    0.000)   24.996
   3.103   (   1.966    1.668    0.000)    2.578
   4.131   (  -3.718   -4.943    0.000)    6.185
   4.286   (   0.512    0.318    0.000)    0.602
   4.338   (  -0.071    2.525    0.000)    2.526
   4.369   (   0.269    1.676    0.000)    1.697
   4.572   (   2.393    2.571    0.000)    3.512
   4.801   (   0.778   -1.549    0.000)    1.733
======================= Grid point 16 (10/84) =======================
q-point: ( 0.23  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.823   (  -1.047   12.555    0.000)   12.599
   0.877   (   8.050    6.969    0.000)   10.647
   1.086   (  -0.859    8.505    0.000)    8.549
   1.438   (   0.568   -0.173    0.000)    0.593
   2.024   (  18.529   12.157    0.000)   22.161
   3.216   (   7.238    5.036    0.000)    8.817
   4.003   (  -4.937   -6.345    0.000)    8.039
   4.303   (   1.070    0.455    0.000)    1.163
   4.362   (  -0.298    2.826    0.000)    2.841
   4.385   (  -0.793    2.778    0.000)    2.889
   4.616   (   1.031    0.342    0.000)    1.086
   4.773   (  -0.881   -4.205    0.000)    4.296
======================= Grid point 17 (11/84) =======================
q-point: ( 0.31  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.889   (  -4.375   10.999    0.000)   11.837
   1.062   (   6.359    5.286    0.000)    8.269
   1.120   (  -3.778    7.059    0.000)    8.007
   1.416   (  -1.819   -0.600    0.000)    1.915
   2.420   (  16.194    9.090    0.000)   18.571
   3.421   (  10.045    5.764    0.000)   11.581
   3.844   (  -5.735   -7.565    0.000)    9.493
   4.325   (   1.018   -0.296    0.000)    1.060
   4.378   (  -0.749    2.472    0.000)    2.583
   4.404   (  -0.345    2.989    0.000)    3.008
   4.605   (  -1.897   -1.843    0.000)    2.645
   4.698   (  -2.318   -4.091    0.000)    4.702
======================= Grid point 18 (12/84) =======================
q-point: ( 0.38  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.904   (  -5.583   10.067    0.000)   11.512
   1.083   (  -6.242    6.314    0.000)    8.878
   1.211   (   4.651    3.617    0.000)    5.892
   1.380   (  -1.945   -0.207    0.000)    1.956
   2.733   (  11.586    5.974    0.000)   13.035
   3.607   (   2.515   -0.887    0.000)    2.667
   3.709   (   2.013   -2.302    0.000)    3.058
   4.307   (  -1.391   -3.869    0.000)    4.112
   4.390   (  -0.481    1.958    0.000)    2.016
   4.427   (  -0.367    2.885    0.000)    2.908
   4.529   (  -3.449   -1.347    0.000)    3.703
   4.644   (  -1.111   -1.499    0.000)    1.866
======================= Grid point 19 (13/84) =======================
q-point: ( 0.46  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.906   (  -5.632    9.755    0.000)   11.264
   1.055   (  -4.390    7.603    0.000)    8.780
   1.279   (  -0.462    0.799    0.000)    0.923
   1.354   (  -0.301    0.521    0.000)    0.602
   2.875   (   0.289   -0.501    0.000)    0.579
   3.551   (   1.995   -3.455    0.000)    3.990
   3.794   (   0.891   -1.543    0.000)    1.782
   4.249   (   1.698   -2.940    0.000)    3.395
   4.396   (  -0.869    1.505    0.000)    1.737
   4.441   (  -1.613    2.793    0.000)    3.226
   4.496   (  -0.474    0.821    0.000)    0.948
   4.627   (   0.170   -0.294    0.000)    0.340
======================= Grid point 29 (14/84) =======================
q-point: ( 0.15  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.841   (   5.272    9.132    0.000)   10.545
   0.919   (   5.918   10.250    0.000)   11.835
   1.150   (   4.144    7.178    0.000)    8.289
   1.428   (   0.858    1.486    0.000)    1.716
   1.908   (  10.509   18.202    0.000)   21.018
   3.178   (   3.525    6.105    0.000)    7.049
   4.012   (  -4.155   -7.197    0.000)    8.310
   4.296   (   0.409    0.709    0.000)    0.819
   4.385   (   1.250    2.164    0.000)    2.499
   4.412   (   1.275    2.209    0.000)    2.551
   4.609   (   0.659    1.142    0.000)    1.319
   4.749   (  -2.042   -3.537    0.000)    4.085
======================= Grid point 30 (15/84) =======================
q-point: ( 0.23  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.010   (   5.144    6.972    0.000)    8.664
   1.063   (  -0.200   10.402    0.000)   10.404
   1.247   (  -0.932    6.897    0.000)    6.959
   1.437   (  -0.678    0.375    0.000)    0.774
   2.263   (  13.489   12.768    0.000)   18.573
   3.334   (   7.058    7.092    0.000)   10.005
   3.855   (  -5.253   -8.820    0.000)   10.266
   4.309   (   0.504   -0.047    0.000)    0.506
   4.414   (  -0.431    1.984    0.000)    2.030
   4.445   (  -0.219    2.877    0.000)    2.886
   4.613   (  -0.739   -0.527    0.000)    0.908
   4.677   (  -1.022   -4.829    0.000)    4.936
======================= Grid point 31 (16/84) =======================
q-point: ( 0.31  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.098   (  -3.711    8.888    0.000)    9.632
   1.174   (   4.063    6.303    0.000)    7.499
   1.255   (  -4.997    5.869    0.000)    7.708
   1.417   (  -2.110    0.677    0.000)    2.216
   2.594   (  12.688    9.173    0.000)   15.657
   3.529   (   7.445    4.857    0.000)    8.889
   3.679   (  -4.043   -8.512    0.000)    9.423
   4.291   (  -1.313   -3.848    0.000)    4.066
   4.415   (  -0.609    1.168    0.000)    1.317
   4.465   (  -0.517    2.819    0.000)    2.866
   4.572   (  -2.643   -0.579    0.000)    2.706
   4.639   (  -0.118   -1.864    0.000)    1.868
======================= Grid point 32 (17/84) =======================
q-point: ( 0.38  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.105   (  -5.341    9.294    0.000)   10.719
   1.218   (  -5.435    6.225    0.000)    8.264
   1.285   (   2.389    4.083    0.000)    4.730
   1.382   (  -2.488    0.348    0.000)    2.512
   2.836   (   6.281    4.962    0.000)    8.004
   3.489   (  -1.618   -8.291    0.000)    8.447
   3.738   (   5.568    2.845    0.000)    6.253
   4.202   (  -1.404   -6.263    0.000)    6.419
   4.421   (  -0.444    1.372    0.000)    1.442
   4.488   (  -0.088    2.899    0.000)    2.900
   4.530   (  -2.495    0.827    0.000)    2.629
   4.625   (  -0.298   -0.231    0.000)    0.378
======================= Grid point 43 (18/84) =======================
q-point: ( 0.23  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.138   (   3.561    6.167    0.000)    7.121
   1.208   (   2.747    4.757    0.000)    5.493
   1.330   (   0.963    1.669    0.000)    1.927
   1.446   (   0.063    0.109    0.000)    0.126
   2.524   (   8.048   13.939    0.000)   16.096
   3.486   (   4.624    8.009    0.000)    9.248
   3.672   (  -5.384   -9.326    0.000)   10.768
   4.282   (  -1.996   -3.457    0.000)    3.991
   4.427   (  -0.054   -0.093    0.000)    0.108
   4.485   (   0.775    1.342    0.000)    1.549
   4.602   (  -0.621   -1.075    0.000)    1.242
   4.619   (  -0.319   -0.553    0.000)    0.639
======================= Grid point 44 (19/84) =======================
q-point: ( 0.31  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.228   (  -0.478    5.431    0.000)    5.452
   1.281   (   2.184    3.729    0.000)    4.321
   1.327   (  -2.321    1.763    0.000)    2.914
   1.430   (  -2.094   -0.007    0.000)    2.094
   2.785   (   8.026   10.614    0.000)   13.307
   3.451   (  -4.755  -10.255    0.000)   11.304
   3.687   (   5.743    6.951    0.000)    9.016
   4.176   (  -3.409   -7.374    0.000)    8.124
   4.436   (   0.279    1.332    0.000)    1.361
   4.504   (   0.291    1.125    0.000)    1.162
   4.571   (  -1.717   -0.585    0.000)    1.814
   4.630   (   0.308    1.024    0.000)    1.070
======================= Grid point 45 (20/84) =======================
q-point: ( 0.38  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.246   (  -3.162    5.476    0.000)    6.324
   1.299   (  -1.356    2.348    0.000)    2.711
   1.353   (  -1.011    1.751    0.000)    2.022
   1.391   (  -0.398    0.690    0.000)    0.796
   2.933   (  -2.845    4.928    0.000)    5.691
   3.332   (   4.184   -7.248    0.000)    8.369
   3.796   (  -1.773    3.071    0.000)    3.546
   4.086   (   3.041   -5.267    0.000)    6.082
   4.448   (  -0.810    1.404    0.000)    1.621
   4.518   (   0.137   -0.238    0.000)    0.275
   4.544   (  -0.564    0.977    0.000)    1.128
   4.639   (  -0.700    1.212    0.000)    1.400
======================= Grid point 58 (21/84) =======================
q-point: ( 0.31  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.313   (   1.138    1.970    0.000)    2.275
   1.326   (   1.101    1.908    0.000)    2.203
   1.343   (   0.644    1.115    0.000)    1.287
   1.404   (  -1.474   -2.553    0.000)    2.948
   2.996   (   5.962   10.327    0.000)   11.924
   3.247   (  -5.832  -10.101    0.000)   11.663
   3.836   (   4.441    7.691    0.000)    8.881
   4.025   (  -4.402   -7.625    0.000)    8.804
   4.467   (   0.970    1.680    0.000)    1.940
   4.506   (  -0.485   -0.841    0.000)    0.971
   4.556   (  -0.255   -0.441    0.000)    0.510
   4.653   (   0.449    0.777    0.000)    0.898
======================= Grid point 181 (22/84) =======================
q-point: ( 0.00  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.261   (   0.000    0.000   13.065)   13.065
   0.261   (   0.000   -0.000   13.065)   13.065
   0.673   (   0.000    0.000   34.685)   34.685
   0.784   (  -0.000   -0.000    0.384)    0.384
   0.784   (   0.000   -0.000    0.384)    0.384
   3.073   (   0.000    0.000  -11.536)   11.536
   4.300   (  -0.000    0.000    0.027)    0.027
   4.300   (   0.000    0.000    0.027)    0.027
   4.332   (  -0.000    0.000   -0.336)    0.336
   4.332   (   0.000   -0.000   -0.336)    0.336
   4.444   (  -0.000    0.000    6.112)    6.112
   4.756   (   0.000    0.000    0.656)    0.656
======================= Grid point 182 (23/84) =======================
q-point: ( 0.08  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.357   (   7.299    4.214    9.465)   12.673
   0.485   (   6.998    4.040   15.429)   17.417
   0.811   (   2.413    1.393    0.573)    2.844
   0.836   (  18.225   10.522   20.177)   29.154
   0.966   (  14.106    8.144    0.450)   16.295
   3.041   (  -2.434   -1.405  -11.124)   11.474
   4.293   (  -0.912   -0.527    1.520)    1.849
   4.297   (  -0.048   -0.028    0.064)    0.085
   4.339   (  -0.136   -0.078    3.094)    3.098
   4.342   (   0.589    0.340   -0.409)    0.793
   4.464   (   2.011    1.161    4.227)    4.822
   4.743   (   0.455    0.263   -1.025)    1.152
======================= Grid point 183 (24/84) =======================
q-point: ( 0.15  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.537   (   8.033    4.638    6.258)   11.189
   0.624   (   6.705    3.871   11.434)   13.809
   0.886   (   4.094    2.364    0.746)    4.786
   1.234   (  12.502    7.218    0.857)   14.461
   1.372   (  22.236   12.838    8.708)   27.112
   2.998   (  -0.444   -0.256   -9.952)    9.965
   4.245   (  -3.330   -1.923    3.703)    5.339
   4.304   (   0.647    0.374    0.049)    0.749
   4.331   (  -0.340   -0.197    4.185)    4.203
   4.351   (   0.175    0.101   -0.440)    0.484
   4.530   (   3.398    1.962    1.373)    4.157
   4.773   (   1.510    0.872   -2.939)    3.418
======================= Grid point 184 (25/84) =======================
q-point: ( 0.23  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.696   (   5.826    3.364    4.754)    8.238
   0.784   (   7.350    4.243    4.048)    9.402
   0.979   (   3.733    2.155    0.625)    4.356
   1.400   (   2.739    1.581   -2.491)    4.026
   1.867   (  21.203   12.241    8.201)   25.820
   3.054   (   5.617    3.243   -8.298)   10.532
   4.153   (  -4.735   -2.734    5.141)    7.505
   4.324   (   1.045    0.604   -0.036)    1.208
   4.330   (   0.318    0.184    3.320)    3.341
   4.352   (   0.079    0.045   -0.334)    0.346
   4.594   (   1.986    1.147   -0.274)    2.310
   4.787   (  -0.636   -0.367   -3.164)    3.248
======================= Grid point 185 (26/84) =======================
q-point: ( 0.31  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.791   (   2.676    1.545    4.013)    5.065
   0.944   (   6.806    3.930   -1.514)    8.004
   1.037   (   1.212    0.700    0.509)    1.489
   1.394   (  -2.198   -1.269   -4.077)    4.803
   2.308   (  18.090   10.444    7.580)   22.222
   3.232   (   9.505    5.488   -7.175)   13.113
   4.036   (  -5.593   -3.229    7.005)    9.528
   4.343   (   0.633    0.366   -0.051)    0.733
   4.346   (   1.003    0.579    2.517)    2.770
   4.364   (   1.014    0.586   -0.261)    1.200
   4.609   (  -0.803   -0.464   -1.543)    1.800
   4.737   (  -3.634   -2.098   -2.910)    5.107
======================= Grid point 186 (27/84) =======================
q-point: ( 0.38  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.824   (   0.589    0.340    3.645)    3.708
   1.031   (  -1.504   -0.868    0.682)    1.866
   1.083   (   5.441    3.141   -5.174)    8.139
   1.337   (  -2.434   -1.405   -5.485)    6.163
   2.672   (  14.191    8.193    7.412)   17.984
   3.441   (   8.112    4.683   -7.214)   11.823
   3.916   (  -4.397   -2.539    9.463)   10.739
   4.359   (   0.768    0.443    0.005)    0.887
   4.364   (   0.111    0.064    2.235)    2.239
   4.388   (   1.012    0.585   -0.302)    1.207
   4.547   (  -4.608   -2.661   -2.856)    6.039
   4.655   (  -2.853   -1.647   -2.238)    3.982
======================= Grid point 187 (28/84) =======================
q-point: ( 0.46  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.830   (   0.161    0.093    3.418)    3.423
   0.991   (  -1.410   -0.814    1.054)    1.939
   1.180   (   2.890    1.668   -6.919)    7.682
   1.291   (  -1.606   -0.927   -6.609)    6.864
   2.918   (   6.689    3.862    7.338)   10.654
   3.535   (   0.310    0.179   -4.947)    4.960
   3.890   (   1.816    1.049    8.623)    8.874
   4.315   (  -3.727   -2.152    0.971)    4.412
   4.373   (   0.422    0.243   -0.002)    0.487
   4.405   (   0.384    0.222   -0.323)    0.549
   4.476   (  -0.858   -0.496   -1.839)    2.089
   4.613   (  -0.979   -0.565   -2.224)    2.495
======================= Grid point 195 (29/84) =======================
q-point: ( 0.08  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.527   (   6.276   10.870    6.387)   14.084
   0.578   (   3.525    6.106   13.054)   14.837
   0.881   (   3.284    5.687    0.882)    6.626
   1.165   (   6.843   11.853   -0.803)   13.710
   1.202   (  12.963   22.452   12.268)   28.682
   3.004   (  -0.950   -1.645  -10.343)   10.516
   4.261   (  -1.642   -2.844    3.100)    4.516
   4.309   (   0.336    0.583    0.693)    0.966
   4.337   (  -0.065   -0.112    3.181)    3.184
   4.346   (   0.432    0.749    0.079)    0.868
   4.511   (   1.895    3.282    1.901)    4.240
   4.757   (   0.563    0.976   -2.579)    2.815
======================= Grid point 196 (30/84) =======================
q-point: ( 0.15  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.706   (   3.915    8.478    5.014)   10.599
   0.736   (   4.890    9.756    5.622)   12.276
   0.998   (   2.048    8.188    0.598)    8.461
   1.350   (   5.076    4.079   -1.914)    6.788
   1.661   (  17.962   16.479    9.192)   26.052
   3.012   (   2.489    2.211   -8.941)    9.541
   4.183   (  -3.247   -4.372    4.642)    7.155
   4.318   (   0.331    0.288    2.268)    2.310
   4.342   (  -0.195    1.831    0.971)    2.082
   4.367   (   0.191    1.723   -0.154)    1.740
   4.573   (   1.953    2.065    0.250)    2.853
   4.769   (   0.733   -1.449   -3.037)    3.444
======================= Grid point 197 (31/84) =======================
q-point: ( 0.23  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.842   (   1.719    8.749    1.464)    9.036
   0.889   (   3.591    7.855    2.329)    8.946
   1.092   (  -0.911    8.479    0.490)    8.541
   1.404   (  -0.007   -0.663   -3.418)    3.482
   2.106   (  18.114   12.225    8.018)   23.277
   3.136   (   7.511    5.189   -7.566)   11.857
   4.072   (  -4.141   -5.540    6.200)    9.289
   4.329   (   0.828    0.219    2.333)    2.485
   4.362   (  -0.231    2.551    0.017)    2.562
   4.383   (  -0.784    2.786   -0.136)    2.898
   4.606   (   0.586    0.008   -0.976)    1.139
   4.743   (  -0.831   -4.001   -2.950)    5.040
======================= Grid point 198 (32/84) =======================
q-point: ( 0.31  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.915   (  -3.416    9.585    2.241)   10.419
   1.023   (   4.998    4.617   -3.150)    7.498
   1.124   (  -4.245    7.146    0.444)    8.324
   1.368   (  -2.465   -0.996   -4.787)    5.476
   2.498   (  16.032    9.200    7.601)   19.986
   3.336   (   9.319    5.122   -7.121)   12.798
   3.943   (  -4.063   -6.179    8.392)   11.185
   4.345   (   0.744   -0.435    1.824)    2.017
   4.378   (  -0.428    2.162    0.130)    2.207
   4.401   (  -0.409    2.956   -0.270)    2.996
   4.584   (  -2.419   -2.028   -2.165)    3.827
   4.673   (  -2.068   -3.707   -2.405)    4.879
======================= Grid point 199 (33/84) =======================
q-point: ( 0.38  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.930   (  -5.054    9.125    2.490)   10.724
   1.091   (  -5.669    6.274    0.542)    8.473
   1.145   (   3.506    3.385   -5.917)    7.665
   1.318   (  -2.489   -0.623   -6.197)    6.707
   2.810   (  11.488    6.033    7.447)   14.961
   3.497   (   4.408    0.625   -6.544)    7.915
   3.853   (   0.700   -3.321    9.753)   10.327
   4.319   (  -1.743   -4.547    0.933)    4.958
   4.392   (  -0.362    1.764    0.293)    1.824
   4.423   (  -0.312    2.803   -0.337)    2.841
   4.507   (  -3.293   -0.539   -2.115)    3.951
   4.623   (  -1.164   -1.419   -2.113)    2.799
======================= Grid point 200 (34/84) =======================
q-point: ( 0.46  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.932   (  -5.163    8.943    2.478)   10.619
   1.066   (  -4.354    7.542    1.040)    8.771
   1.204   (  -0.517    0.896   -7.018)    7.094
   1.284   (  -0.218    0.378   -6.919)    6.933
   2.951   (   0.259   -0.449    7.369)    7.387
   3.509   (   1.413   -2.447   -3.206)    4.273
   3.877   (   1.405   -2.434    7.189)    7.719
   4.253   (   1.837   -3.182    0.163)    3.678
   4.397   (  -0.840    1.456    0.220)    1.695
   4.437   (  -1.569    2.718   -0.361)    3.159
   4.482   (  -0.628    1.087   -1.283)    1.795
   4.605   (   0.106   -0.184   -2.166)    2.176
======================= Grid point 210 (35/84) =======================
q-point: ( 0.15  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.838   (   3.810    6.599    0.580)    7.642
   0.950   (   5.554    9.620    3.055)   11.520
   1.157   (   4.161    7.207    0.633)    8.346
   1.396   (   0.522    0.904   -3.271)    3.434
   1.997   (  10.341   17.912    8.611)   22.404
   3.096   (   3.716    6.436   -7.823)   10.790
   4.079   (  -3.630   -6.288    6.000)    9.419
   4.323   (   0.239    0.413    2.510)    2.555
   4.382   (   1.121    1.942   -0.272)    2.259
   4.409   (   1.265    2.191   -0.246)    2.541
   4.601   (   0.423    0.733   -0.692)    1.093
   4.720   (  -1.961   -3.397   -2.873)    4.862
======================= Grid point 211 (36/84) =======================
q-point: ( 0.23  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.964   (   4.111    5.405   -3.744)    7.755
   1.085   (  -0.143    9.334    1.990)    9.545
   1.254   (  -1.034    7.081    0.657)    7.186
   1.395   (  -1.056    0.027   -4.153)    4.285
   2.346   (  13.223   12.826    8.085)   20.118
   3.252   (   6.802    6.697   -7.212)   11.964
   3.944   (  -4.212   -7.419    7.856)   11.597
   4.331   (   0.343   -0.269    2.113)    2.157
   4.407   (  -0.423    1.504   -0.776)    1.744
   4.443   (  -0.220    2.873   -0.220)    2.890
   4.597   (  -1.098   -0.682   -1.538)    2.009
   4.653   (  -0.762   -4.386   -2.369)    5.043
======================= Grid point 212 (37/84) =======================
q-point: ( 0.31  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.067   (   1.248    5.196   -4.873)    7.233
   1.151   (  -0.909    7.973   -0.279)    8.029
   1.264   (  -5.017    6.150    0.783)    7.975
   1.362   (  -2.897    0.411   -5.265)    6.024
   2.674   (  12.489    9.226    7.738)   17.348
   3.422   (   5.955    3.286   -7.030)    9.781
   3.815   (  -1.895   -6.146    9.826)   11.744
   4.303   (  -1.807   -4.561    1.074)    5.022
   4.407   (  -0.130    0.819   -0.805)    1.156
   4.462   (  -0.613    2.793   -0.317)    2.877
   4.557   (  -2.542    0.028   -1.477)    2.940
   4.620   (  -0.179   -1.533   -1.811)    2.379
======================= Grid point 213 (38/84) =======================
q-point: ( 0.38  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.101   (  -3.401    6.890   -1.768)    7.885
   1.222   (  -0.637    5.787   -3.129)    6.609
   1.231   (  -3.714    5.907   -0.382)    6.988
   1.312   (  -3.085    0.088   -6.708)    7.384
   2.913   (   6.182    4.925    7.463)   10.870
   3.458   (   0.325   -4.292   -2.395)    4.925
   3.812   (   3.823   -0.632    6.586)    7.641
   4.204   (  -1.543   -6.522    0.021)    6.702
   4.417   (  -0.127    1.015   -0.401)    1.098
   4.482   (  -0.177    2.618   -0.647)    2.702
   4.522   (  -2.387    1.244   -0.760)    2.797
   4.608   (  -0.495    0.217   -1.656)    1.742
======================= Grid point 224 (39/84) =======================
q-point: ( 0.23  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.067   (   2.953    5.115   -6.582)    8.843
   1.213   (   2.311    4.003   -0.177)    4.626
   1.341   (   1.022    1.771    0.971)    2.263
   1.400   (  -0.060   -0.104   -4.131)    4.133
   2.606   (   7.971   13.807    7.966)   17.822
   3.384   (   3.538    6.129   -7.090)   10.018
   3.801   (  -3.732   -6.465    9.715)   12.252
   4.295   (  -2.443   -4.231    1.181)    5.026
   4.413   (  -0.023   -0.040   -1.437)    1.438
   4.483   (   0.738    1.278   -0.265)    1.499
   4.586   (  -0.575   -0.996   -1.537)    1.920
   4.606   (  -0.004   -0.007   -1.305)    1.305
======================= Grid point 225 (40/84) =======================
q-point: ( 0.31  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.156   (   1.807    4.405   -7.691)    9.046
   1.264   (   0.355    3.331   -2.332)    4.082
   1.338   (  -2.558    1.602    1.074)    3.204
   1.376   (  -2.751    0.362   -3.885)    4.774
   2.864   (   7.844   10.366    7.548)   15.032
   3.425   (  -1.263   -4.046   -2.172)    4.763
   3.758   (   2.664    1.527    6.496)    7.185
   4.178   (  -3.576   -7.632    0.072)    8.428
   4.425   (   0.467    1.188   -1.114)    1.694
   4.497   (   0.060    0.717   -0.794)    1.072
   4.561   (  -1.375   -0.183   -0.804)    1.603
   4.622   (   0.192    1.433   -0.851)    1.677
======================= Grid point 226 (41/84) =======================
q-point: ( 0.38  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.196   (  -2.073    3.591   -6.738)    7.912
   1.292   (  -0.489    0.847   -4.944)    5.040
   1.329   (  -1.898    3.287   -0.272)    3.805
   1.330   (  -1.115    1.932   -1.873)    2.913
   3.007   (  -2.688    4.656    7.153)    8.948
   3.357   (   3.132   -5.425    2.181)    6.633
   3.825   (  -0.982    1.701    2.949)    3.543
   4.086   (   3.014   -5.220   -0.070)    6.028
   4.437   (  -0.632    1.095   -1.063)    1.652
   4.502   (   0.279   -0.483   -1.661)    1.752
   4.545   (  -0.720    1.247    0.146)    1.447
   4.631   (  -0.874    1.514   -0.771)    1.911
======================= Grid point 239 (42/84) =======================
q-point: ( 0.31  0.31  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.238   (   2.111    3.656   -8.679)    9.652
   1.294   (  -0.043   -0.074   -6.909)    6.909
   1.342   (  -0.306   -0.529    1.335)    1.469
   1.365   (  -0.575   -0.996   -0.149)    1.159
   3.066   (   5.604    9.707    6.715)   13.066
   3.290   (  -4.798   -8.310    4.020)   10.403
   3.852   (   3.753    6.501    1.616)    7.679
   4.026   (  -4.302   -7.452    0.091)    8.605
   4.453   (   0.882    1.527   -1.464)    2.292
   4.490   (  -0.652   -1.129   -1.668)    2.117
   4.559   (  -0.003   -0.004    0.316)    0.316
   4.648   (   0.476    0.825   -0.547)    1.099
======================= Grid point 362 (43/84) =======================
q-point: ( 0.00  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.487   (   0.000    0.000   10.494)   10.494
   0.487   (  -0.000   -0.000   10.494)   10.494
   0.784   (   0.000   -0.000   -0.606)    0.606
   0.784   (   0.000    0.000   -0.606)    0.606
   1.310   (  -0.000    0.000   31.944)   31.944
   2.780   (   0.000    0.000  -18.670)   18.670
   4.302   (   0.000   -0.000    0.133)    0.133
   4.302   (   0.000    0.000    0.133)    0.133
   4.324   (   0.000   -0.000   -0.523)    0.523
   4.324   (   0.000    0.000   -0.523)    0.523
   4.576   (   0.000   -0.000    6.976)    6.976
   4.766   (  -0.000   -0.000    0.169)    0.169
======================= Grid point 363 (44/84) =======================
q-point: ( 0.08  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.544   (   4.742    2.738    9.363)   10.847
   0.699   (  12.294    7.098    8.494)   16.543
   0.816   (   2.722    1.571   -0.349)    3.162
   0.949   (  12.231    7.062   -0.873)   14.150
   1.352   (   5.998    3.463   29.081)   29.894
   2.758   (  -1.524   -0.880  -18.050)   18.135
   4.300   (   0.007    0.004    0.234)    0.235
   4.314   (   0.763    0.440    0.733)    1.145
   4.331   (   0.474    0.274   -0.667)    0.863
   4.364   (   2.472    1.427   -0.156)    2.858
   4.570   (  -0.262   -0.151    6.137)    6.144
   4.736   (  -1.810   -1.045   -0.067)    2.091
======================= Grid point 364 (45/84) =======================
q-point: ( 0.15  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.675   (   6.317    3.647    7.458)   10.432
   0.867   (   3.259    1.882   12.102)   12.674
   0.896   (   4.165    2.405   -0.005)    4.809
   1.215   (   9.606    5.546   -2.072)   11.284
   1.638   (  18.038   10.414   18.023)   27.543
   2.747   (   1.494    0.863  -15.889)   15.983
   4.306   (   0.569    0.329    0.229)    0.696
   4.309   (  -1.512   -0.873    2.675)    3.194
   4.340   (   0.234    0.135   -0.719)    0.768
   4.398   (   0.341    0.197    1.794)    1.836
   4.577   (   0.898    0.519    3.512)    3.662
   4.718   (   0.180    0.104   -2.296)    2.305
======================= Grid point 365 (46/84) =======================
q-point: ( 0.23  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.804   (   4.878    2.816    6.059)    8.273
   0.910   (   1.348    0.778    8.624)    8.763
   0.987   (   3.598    2.077    0.026)    4.154
   1.332   (   0.868    0.501   -4.503)    4.613
   2.072   (  19.585   11.308   12.678)   25.926
   2.852   (   7.885    4.553  -12.277)   15.285
   4.253   (  -3.128   -1.806    4.769)    5.983
   4.324   (   0.933    0.538    0.053)    1.078
   4.344   (   0.247    0.142   -0.518)    0.591
   4.393   (  -0.345   -0.199    2.233)    2.268
   4.596   (   0.328    0.190    0.794)    0.880
   4.719   (  -0.723   -0.418   -3.393)    3.495
======================= Grid point 366 (47/84) =======================
q-point: ( 0.31  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.884   (   2.176    1.256    5.240)    5.811
   0.942   (   1.655    0.956    1.950)    2.731
   1.044   (   1.245    0.719   -0.009)    1.437
   1.289   (  -3.670   -2.119   -6.592)    7.837
   2.490   (  17.505   10.106   10.795)   22.915
   3.072   (  10.928    6.309   -9.056)   15.532
   4.176   (  -3.676   -2.122    6.924)    8.122
   4.343   (   0.711    0.410    0.010)    0.821
   4.357   (   0.986    0.569   -0.400)    1.207
   4.394   (   0.555    0.321    1.639)    1.760
   4.575   (  -2.474   -1.428   -1.446)    3.202
   4.674   (  -3.110   -1.795   -3.231)    4.831
======================= Grid point 367 (48/84) =======================
q-point: ( 0.38  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.908   (   0.277    0.160    4.737)    4.748
   0.985   (   2.086    1.204   -3.992)    4.663
   1.043   (  -1.123   -0.648    0.308)    1.332
   1.201   (  -3.743   -2.161   -8.182)    9.254
   2.846   (  13.953    8.056   10.000)   18.962
   3.308   (   9.314    5.377   -6.688)   12.664
   4.094   (  -3.337   -1.926    8.630)    9.451
   4.360   (   0.796    0.460    0.113)    0.926
   4.381   (   0.967    0.558   -0.484)    1.217
   4.413   (   0.669    0.386    1.789)    1.949
   4.476   (  -5.900   -3.407   -3.503)    7.661
   4.603   (  -2.710   -1.565   -2.899)    4.265
======================= Grid point 368 (49/84) =======================
q-point: ( 0.46  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.909   (  -0.031   -0.018    4.340)    4.340
   1.012   (  -1.084   -0.626    0.887)    1.534
   1.030   (   1.684    0.972   -7.557)    7.803
   1.132   (  -2.162   -1.248   -9.149)    9.483
   3.088   (   6.572    3.795    9.570)   12.214
   3.451   (   2.968    1.713   -3.973)    5.247
   4.053   (  -0.027   -0.016    7.920)    7.920
   4.319   (  -5.184   -2.993   -1.211)    6.107
   4.374   (   0.404    0.233    0.110)    0.479
   4.396   (   0.378    0.218   -0.515)    0.675
   4.447   (   0.718    0.415   -0.452)    0.944
   4.559   (  -1.140   -0.658   -3.177)    3.438
======================= Grid point 376 (50/84) =======================
q-point: ( 0.08  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.667   (   4.894    8.477    7.562)   12.369
   0.824   (   2.820    4.885   10.260)   11.709
   0.896   (   3.042    5.268    1.077)    6.178
   1.141   (   6.128   10.614   -1.775)   12.384
   1.528   (   8.733   15.127   21.174)   27.448
   2.741   (  -0.009   -0.016  -16.664)   16.664
   4.306   (  -0.401   -0.694    1.426)    1.636
   4.313   (   0.481    0.833    0.136)    0.972
   4.348   (   0.399    0.691    0.101)    0.805
   4.388   (   0.366    0.634    1.144)    1.358
   4.573   (   0.369    0.640    4.296)    4.359
   4.714   (  -0.450   -0.780   -1.573)    1.813
======================= Grid point 377 (51/84) =======================
q-point: ( 0.15  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.823   (   1.799    9.771    6.249)   11.737
   0.890   (   2.546    1.141    9.910)   10.296
   1.006   (   1.993    8.141    0.047)    8.382
   1.297   (   3.899    2.327   -3.658)    5.830
   1.892   (  15.432   15.022   14.458)   25.939
   2.791   (   4.257    4.081  -13.725)   14.938
   4.270   (  -1.893   -2.807    3.855)    5.131
   4.332   (   0.195    1.827    0.154)    1.844
   4.359   (  -0.456    1.746   -0.250)    1.822
   4.391   (   0.284   -0.640    1.877)    2.003
   4.590   (   0.521    0.565    1.697)    1.863
   4.705   (   0.328   -1.595   -3.039)    3.448
======================= Grid point 378 (52/84) =======================
q-point: ( 0.23  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.892   (   1.008    2.159    4.272)    4.892
   0.964   (  -0.071    7.679    4.728)    9.018
   1.099   (  -1.000    8.504    0.135)    8.564
   1.315   (  -1.306   -1.619   -5.656)    6.026
   2.300   (  16.919   11.950   11.646)   23.763
   2.960   (   9.134    6.442  -10.342)   15.228
   4.196   (  -2.718   -3.988    6.079)    7.762
   4.352   (  -0.127    1.964   -0.073)    1.970
   4.370   (  -0.142    1.923    0.092)    1.931
   4.392   (  -0.070    0.664    1.221)    1.392
   4.588   (  -0.795   -1.065   -0.554)    1.440
   4.678   (  -0.733   -3.571   -3.239)    4.877
======================= Grid point 379 (53/84) =======================
q-point: ( 0.31  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.929   (   1.428    1.363   -0.681)    2.088
   1.013   (  -2.495    7.260    2.430)    8.053
   1.131   (  -4.209    7.166    0.028)    8.311
   1.251   (  -4.040   -1.221   -6.872)    8.064
   2.678   (  15.488    9.195   10.423)   20.810
   3.192   (  10.363    5.870   -7.720)   14.193
   4.107   (  -2.603   -4.875    8.073)    9.783
   4.362   (  -0.274    0.450   -0.434)    0.683
   4.390   (   0.090    1.899    0.105)    1.904
   4.403   (   0.273    1.168    1.083)    1.616
   4.534   (  -3.722   -2.272   -2.398)    4.976
   4.620   (  -1.732   -3.110   -2.846)    4.557
======================= Grid point 380 (54/84) =======================
q-point: ( 0.38  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.960   (  -0.319    2.631   -1.784)    3.195
   1.039   (  -1.508    5.856   -1.590)    6.252
   1.109   (  -5.625    6.247    0.044)    8.406
   1.172   (  -4.078   -0.516   -7.781)    8.800
   2.982   (  11.187    6.005    9.789)   16.032
   3.387   (   6.765    2.704   -5.162)    8.929
   4.031   (  -0.312   -4.188    8.528)    9.506
   4.322   (  -2.366   -5.788   -1.198)    6.367
   4.403   (  -0.463    1.274    0.559)    1.467
   4.419   (  -0.127    1.803    0.247)    1.825
   4.463   (  -2.657    1.938   -1.835)    3.766
   4.572   (  -1.396   -1.328   -3.002)    3.567
======================= Grid point 381 (55/84) =======================
q-point: ( 0.46  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.969   (  -2.042    3.538   -0.744)    4.152
   1.071   (  -2.751    4.764   -3.985)    6.793
   1.082   (  -3.550    6.148   -0.877)    7.153
   1.123   (  -0.723    1.252   -7.697)    7.832
   3.121   (   0.240   -0.416    9.489)    9.502
   3.452   (   0.812   -1.406   -2.851)    3.281
   4.021   (   1.824   -3.159    7.181)    8.055
   4.244   (   1.994   -3.453   -1.338)    4.206
   4.405   (  -0.796    1.378    0.548)    1.683
   4.428   (  -1.425    2.468   -0.546)    2.902
   4.459   (  -0.798    1.381   -0.701)    1.742
   4.552   (  -0.069    0.120   -3.142)    3.145
======================= Grid point 391 (56/84) =======================
q-point: ( 0.15  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.891   (   0.400    0.693    5.067)    5.130
   1.021   (   4.566    7.908    3.952)    9.950
   1.168   (   4.260    7.379    0.350)    8.527
   1.311   (  -0.163   -0.283   -5.379)    5.389
   2.204   (   9.739   16.868   12.406)   23.093
   2.911   (   4.760    8.244  -10.984)   14.535
   4.199   (  -2.557   -4.430    5.900)    7.809
   4.367   (   0.770    1.334    0.801)    1.737
   4.379   (   0.044    0.076    0.513)    0.521
   4.403   (   1.232    2.134   -0.408)    2.498
   4.591   (  -0.328   -0.568   -0.100)    0.664
   4.657   (  -1.843   -3.192   -3.155)    4.851
======================= Grid point 392 (57/84) =======================
q-point: ( 0.23  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.911   (   1.262    0.929   -0.846)    1.781
   1.117   (  -0.579    5.203    0.235)    5.240
   1.265   (  -1.486    7.364    0.145)    7.514
   1.304   (  -1.300    1.401   -4.086)    4.511
   2.537   (  12.525   12.616   11.076)   20.946
   3.099   (   7.777    7.778   -8.485)   13.892
   4.102   (  -3.065   -5.677    7.884)   10.187
   4.371   (  -0.398   -0.411    0.154)    0.592
   4.388   (  -0.067    0.019    0.141)    0.157
   4.437   (  -0.129    2.761   -0.339)    2.785
   4.563   (  -2.030   -1.062   -1.738)    2.876
   4.601   (  -0.223   -3.462   -2.602)    4.337
======================= Grid point 393 (58/84) =======================
q-point: ( 0.31  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.950   (   1.916    1.230   -5.539)    5.989
   1.128   (  -1.697    2.300   -4.830)    5.612
   1.265   (  -5.323    5.691   -1.252)    7.893
   1.286   (  -4.894    5.435    0.207)    7.316
   2.853   (  11.979    9.121   10.183)   18.176
   3.299   (   7.548    5.119   -6.025)   10.931
   3.998   (  -1.906   -5.960    8.906)   10.885
   4.310   (  -2.986   -5.957   -1.230)    6.776
   4.395   (   0.893    0.069    0.160)    0.910
   4.454   (  -0.692    2.306   -0.539)    2.467
   4.526   (  -2.180    1.313   -1.346)    2.879
   4.578   (  -0.461   -0.644   -2.396)    2.523
======================= Grid point 394 (59/84) =======================
q-point: ( 0.38  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.990   (   0.982    1.264   -7.769)    7.932
   1.116   (  -0.807    0.098   -9.500)    9.534
   1.229   (  -5.774    7.522    0.853)    9.520
   1.254   (  -5.328    6.638    1.114)    8.585
   3.084   (   5.948    4.693    9.519)   12.166
   3.413   (   2.545    0.041   -2.451)    3.534
   3.950   (   1.975   -3.737    7.228)    8.374
   4.194   (  -1.730   -6.658   -1.280)    6.997
   4.409   (   0.515    0.026   -0.270)    0.582
   4.463   (  -0.297    1.188   -1.389)    1.852
   4.509   (  -2.038    2.392   -0.313)    3.158
   4.570   (  -1.063    1.532   -2.123)    2.826
======================= Grid point 405 (60/84) =======================
q-point: ( 0.23  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.938   (   1.105    1.914   -5.512)    5.939
   1.150   (  -0.229   -0.397   -6.997)    7.012
   1.350   (   1.183    2.050   -0.126)    2.370
   1.361   (   1.218    2.109    1.045)    2.650
   2.790   (   7.701   13.339   10.456)   18.616
   3.256   (   4.657    8.067   -6.386)   11.294
   3.985   (  -3.435   -5.950    9.067)   11.376
   4.306   (  -3.547   -6.143   -1.078)    7.175
   4.387   (   0.265    0.460   -0.333)    0.627
   4.472   (   0.357    0.618   -0.986)    1.217
   4.553   (  -0.224   -0.388   -1.500)    1.566
   4.577   (   0.575    0.996   -1.564)    1.941
======================= Grid point 406 (61/84) =======================
q-point: ( 0.31  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.984   (   1.626    2.333   -8.805)    9.253
   1.130   (  -1.150   -1.060  -10.275)   10.393
   1.358   (  -2.607    3.134    1.342)    4.292
   1.366   (  -2.290    2.096    1.366)    3.392
   3.036   (   7.393    9.704    9.619)   15.536
   3.380   (   2.614    2.623   -2.578)    4.512
   3.899   (  -0.315   -3.399    7.635)    8.363
   4.170   (  -3.742   -7.767   -1.165)    8.699
   4.402   (   0.635    0.605   -1.018)    1.344
   4.468   (  -0.369   -0.637   -2.227)    2.346
   4.552   (  -0.652    1.244   -0.022)    1.405
   4.601   (   0.003    2.109   -1.196)    2.424
======================= Grid point 407 (62/84) =======================
q-point: ( 0.38  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.018   (  -0.878    1.520  -10.041)   10.193
   1.107   (   0.374   -0.647  -11.330)   11.355
   1.350   (  -1.915    3.050    1.647)    3.960
   1.350   (  -2.269    4.198    1.665)    5.055
   3.169   (  -2.292    3.970    8.884)    9.997
   3.395   (   0.952   -1.649    1.351)    2.334
   3.901   (   0.646   -1.118    4.674)    4.849
   4.079   (   2.822   -4.888   -0.824)    5.704
   4.412   (  -0.191    0.331   -1.333)    1.387
   4.460   (   0.496   -0.860   -2.612)    2.794
   4.552   (  -1.008    1.747    0.598)    2.104
   4.612   (  -1.164    2.016   -1.142)    2.593
======================= Grid point 420 (63/84) =======================
q-point: ( 0.31  0.31  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.037   (   1.623    2.811  -10.929)   11.401
   1.098   (  -1.185   -2.052  -11.656)   11.894
   1.377   (  -0.093   -0.160    1.967)    1.976
   1.390   (   0.244    0.423    1.667)    1.737
   3.217   (   4.702    8.144    8.216)   12.487
   3.370   (  -2.110   -3.655    3.630)    5.567
   3.897   (   1.751    3.032    3.000)    4.611
   4.025   (  -3.857   -6.681   -0.388)    7.725
   4.418   (   0.592    1.025   -1.901)    2.239
   4.447   (  -0.692   -1.198   -2.591)    2.937
   4.570   (   0.287    0.496    0.820)    1.000
   4.633   (   0.521    0.902   -0.978)    1.429
======================= Grid point 543 (64/84) =======================
q-point: ( 0.00  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.654   (   0.000    0.000    6.908)    6.908
   0.654   (   0.000   -0.000    6.908)    6.908
   0.750   (   0.000   -0.000   -3.263)    3.263
   0.750   (   0.000    0.000   -3.263)    3.263
   1.882   (   0.000    0.000   28.066)   28.066
   2.375   (   0.000    0.000  -23.655)   23.655
   4.306   (   0.000   -0.000    0.322)    0.322
   4.306   (  -0.000    0.000    0.322)    0.322
   4.314   (  -0.000    0.000   -0.497)    0.497
   4.314   (   0.000    0.000   -0.497)    0.497
   4.690   (  -0.000    0.000    4.730)    4.730
   4.752   (   0.000   -0.000   -1.840)    1.840
======================= Grid point 544 (65/84) =======================
q-point: ( 0.08  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.697   (   3.654    2.109    6.447)    7.705
   0.786   (   3.108    1.794   -2.984)    4.668
   0.833   (  12.362    7.137    5.641)   15.349
   0.914   (  12.009    6.933   -2.936)   14.174
   1.891   (   1.723    0.995   26.927)   27.001
   2.366   (  -0.227   -0.131  -22.880)   22.881
   4.307   (   0.130    0.075    0.491)    0.513
   4.318   (   0.304    0.175   -0.685)    0.769
   4.328   (   1.631    0.942    0.911)    2.092
   4.349   (   2.490    1.438   -1.138)    3.092
   4.672   (  -1.352   -0.781    4.197)    4.478
   4.725   (  -2.022   -1.168   -1.439)    2.743
======================= Grid point 545 (66/84) =======================
q-point: ( 0.15  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.801   (   5.149    2.973    5.436)    8.056
   0.875   (   4.475    2.584   -2.409)    5.701
   1.053   (   5.428    3.134    7.419)    9.712
   1.159   (   7.780    4.492   -4.009)    9.838
   2.018   (  10.434    6.024   20.942)   24.161
   2.408   (   4.912    2.836  -19.456)   20.266
   4.313   (   0.455    0.263    0.515)    0.736
   4.325   (   0.333    0.192   -0.735)    0.829
   4.352   (   0.030    0.017    2.160)    2.160
   4.392   (   0.891    0.515   -1.781)    2.057
   4.645   (  -0.830   -0.479    3.108)    3.253
   4.687   (  -1.036   -0.598   -1.433)    1.867
======================= Grid point 546 (67/84) =======================
q-point: ( 0.23  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.908   (   4.101    2.368    4.527)    6.551
   0.970   (   3.672    2.120   -1.984)    4.682
   1.082   (  -1.669   -0.964    8.752)    8.961
   1.228   (  -1.103   -0.637   -6.516)    6.639
   2.331   (  16.185    9.345   14.178)   23.458
   2.601   (  11.759    6.789  -13.902)   19.433
   4.327   (   0.737    0.426    0.268)    0.892
   4.330   (  -1.623   -0.937    3.540)    4.006
   4.334   (   0.487    0.281   -0.477)    0.737
   4.390   (  -0.722   -0.417   -2.288)    2.435
   4.629   (  -1.102   -0.636    2.274)    2.606
   4.669   (  -1.103   -0.637   -2.026)    2.393
======================= Grid point 547 (68/84) =======================
q-point: ( 0.31  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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   (   1.708    0.986    3.960)    4.424
   1.022   (  -3.078   -1.777    6.035)    7.004
   1.029   (   1.393    0.804   -1.748)    2.376
   1.154   (  -4.609   -2.661   -7.376)    9.095
   2.699   (  16.017    9.247   10.734)   21.384
   2.893   (  13.563    7.831   -9.750)   18.448
   4.290   (  -1.977   -1.141    5.129)    5.614
   4.344   (   0.818    0.472    0.185)    0.962
   4.350   (   0.919    0.531   -0.359)    1.120
   4.375   (  -0.525   -0.303   -3.477)    3.529
   4.580   (  -3.305   -1.908    1.861)    4.246
   4.623   (  -3.024   -1.746   -2.320)    4.192
======================= Grid point 548 (69/84) =======================
q-point: ( 0.38  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.961   (  -2.065   -1.192    1.897)    3.047
   0.989   (  -0.147   -0.085    3.467)    3.471
   1.035   (  -0.635   -0.367   -1.325)    1.514
   1.051   (  -4.177   -2.412   -6.961)    8.469
   3.029   (  13.112    7.570    8.893)   17.558
   3.180   (  11.434    6.602   -7.035)   14.961
   4.237   (  -2.813   -1.624    6.517)    7.282
   4.352   (  -2.308   -1.333   -5.662)    6.258
   4.364   (   0.844    0.487    0.312)    1.024
   4.371   (   0.906    0.523   -0.478)    1.150
   4.495   (  -3.325   -1.920    2.994)    4.869
   4.550   (  -2.975   -1.717   -2.746)    4.398
======================= Grid point 549 (70/84) =======================
q-point: ( 0.46  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.935   (  -0.355   -0.205   -1.638)    1.689
   0.979   (  -1.960   -1.132   -5.963)    6.378
   0.981   (  -0.321   -0.185    2.939)    2.962
   1.015   (  -0.688   -0.397   -0.716)    1.070
   3.256   (   6.077    3.508    7.611)   10.352
   3.372   (   4.854    2.802   -4.764)    7.356
   4.180   (  -1.527   -0.882    5.394)    5.675
   4.266   (  -3.730   -2.154   -3.838)    5.769
   4.378   (   0.385    0.223    0.322)    0.550
   4.386   (   0.376    0.217   -0.503)    0.665
   4.463   (  -0.133   -0.077    1.692)    1.699
   4.503   (  -1.129   -0.652   -2.549)    2.863
======================= Grid point 557 (71/84) =======================
q-point: ( 0.08  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.794   (   4.009    6.945    5.431)    9.685
   0.867   (   3.548    6.145   -2.115)    7.404
   1.001   (   4.200    7.275    6.280)   10.488
   1.092   (   5.312    9.201   -3.489)   11.182
   1.961   (   4.257    7.374   23.148)   24.664
   2.383   (   1.588    2.751  -20.751)   20.992
   4.317   (   0.423    0.733    0.291)    0.895
   4.320   (   0.250    0.433    0.150)    0.522
   4.355   (   0.555    0.961    0.886)    1.420
   4.380   (   0.776    1.344   -1.407)    2.095
   4.650   (  -0.673   -1.166    3.297)    3.561
   4.692   (  -1.041   -1.803   -1.234)    2.421
======================= Grid point 558 (72/84) =======================
q-point: ( 0.15  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.928   (   1.796    8.975    4.538)   10.216
   0.990   (   1.723    8.156   -1.932)    8.557
   1.073   (   1.414   -0.828    8.599)    8.754
   1.209   (   2.417    0.351   -5.716)    6.216
   2.192   (  11.491   11.494   16.508)   23.166
   2.505   (   7.371    7.347  -16.039)   19.119
   4.323   (  -0.294   -0.406    1.605)    1.682
   4.335   (   0.216    1.198    0.053)    1.218
   4.359   (  -0.739    1.121    0.772)    1.549
   4.390   (   0.109   -0.589   -1.632)    1.739
   4.631   (  -0.561   -1.118    2.178)    2.512
   4.664   (  -0.336   -1.829   -1.519)    2.401
======================= Grid point 559 (73/84) =======================
q-point: ( 0.23  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.992   (  -1.399    1.345    5.202)    5.552
   1.041   (  -1.099    3.403    1.615)    3.924
   1.103   (  -1.318    7.187    1.773)    7.518
   1.198   (  -2.192   -1.037   -6.405)    6.849
   2.529   (  14.818   10.747   12.076)   21.929
   2.753   (  11.936    8.643  -11.369)   18.613
   4.292   (  -1.297   -2.459    4.022)    4.889
   4.343   (  -0.298    0.305   -1.145)    1.222
   4.369   (  -0.354    2.642    0.178)    2.672
   4.387   (  -0.246    0.888   -1.562)    1.814
   4.598   (  -1.582   -2.495    1.492)    3.310
   4.631   (  -1.028   -3.197   -1.901)    3.860
======================= Grid point 560 (74/84) =======================
q-point: ( 0.31  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.952   (  -1.116   -2.761    3.449)    4.557
   1.032   (  -2.194   -0.155   -2.788)    3.552
   1.122   (  -4.961    7.139    0.200)    8.696
   1.152   (  -5.172    5.256   -2.621)    7.826
   2.873   (  14.339    8.637    9.699)   19.346
   3.042   (  12.468    7.397   -8.190)   16.651
   4.239   (  -1.513   -3.834    5.784)    7.102
   4.328   (  -0.406   -2.457   -3.412)    4.224
   4.390   (  -0.564    3.017    0.221)    3.077
   4.401   (  -0.208    2.414   -0.956)    2.604
   4.529   (  -3.147   -2.798    1.730)    4.553
   4.571   (  -2.104   -3.062   -2.386)    4.415
======================= Grid point 561 (75/84) =======================
q-point: ( 0.38  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.918   (  -0.063   -1.954    0.115)    1.958
   0.983   (  -1.694   -2.113   -5.871)    6.465
   1.095   (  -6.055    7.609    1.470)    9.835
   1.119   (  -5.882    6.997   -0.720)    9.169
   3.157   (  10.434    5.557    8.156)   14.362
   3.288   (   8.914    4.446   -5.717)   11.485
   4.170   (  -1.185   -4.760    6.067)    7.802
   4.267   (  -2.074   -5.661   -4.223)    7.361
   4.407   (  -0.708    2.548    0.251)    2.657
   4.416   (  -0.822    2.142   -0.707)    2.400
   4.470   (  -1.598    0.098    2.094)    2.636
   4.517   (  -1.621   -1.176   -2.706)    3.366
======================= Grid point 562 (76/84) =======================
q-point: ( 0.46  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.910   (   0.749   -1.298   -1.694)    2.262
   0.954   (   0.497   -0.860   -6.096)    6.177
   1.079   (  -4.830    8.366    1.799)    9.826
   1.101   (  -4.474    7.749   -0.321)    8.954
   3.285   (   0.299   -0.518    7.248)    7.272
   3.390   (   0.465   -0.806   -4.023)    4.130
   4.134   (   1.923   -3.331    4.608)    6.003
   4.203   (   1.930   -3.343   -2.839)    4.792
   4.413   (  -1.008    1.747    0.267)    2.034
   4.419   (  -1.208    2.092   -0.381)    2.445
   4.463   (  -0.491    0.850    1.076)    1.456
   4.497   (  -0.276    0.478   -2.434)    2.495
======================= Grid point 572 (77/84) =======================
q-point: ( 0.15  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.012   (  -1.873   -3.244    6.732)    7.704
   1.079   (   2.028    3.512    1.109)    4.204
   1.167   (   4.399    7.620   -0.259)    8.802
   1.203   (   0.746    1.292   -5.148)    5.360
   2.448   (   8.448   14.633   12.945)   21.286
   2.689   (   6.575   11.388  -12.235)   17.961
   4.292   (  -1.422   -2.464    3.901)    4.828
   4.350   (  -0.319   -0.553   -2.048)    2.146
   4.389   (   0.706    1.222    0.151)    1.419
   4.394   (   1.172    2.030   -0.449)    2.386
   4.599   (  -1.279   -2.215    0.749)    2.666
   4.616   (  -1.749   -3.030   -1.192)    3.696
======================= Grid point 573 (78/84) =======================
q-point: ( 0.23  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.947   (  -1.680   -3.091    4.428)    5.655
   1.056   (  -2.280   -2.716   -5.900)    6.884
   1.262   (  -1.222    7.461    0.616)    7.585
   1.276   (  -0.677    7.101   -0.508)    7.152
   2.746   (  11.340   11.655   10.517)   19.366
   2.932   (   9.616    9.853   -9.135)   16.522
   4.231   (  -2.131   -4.009    5.655)    7.252
   4.321   (  -1.308   -2.395   -3.886)    4.749
   4.414   (   0.173    1.210    1.084)    1.633
   4.429   (   0.064    2.364   -0.592)    2.437
   4.546   (  -1.946   -1.991    0.167)    2.789
   4.563   (  -0.730   -2.969   -1.430)    3.376
======================= Grid point 574 (79/84) =======================
q-point: ( 0.31  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.903   (  -0.546   -2.082    1.211)    2.470
   0.991   (  -1.989   -2.668   -7.273)    7.998
   1.277   (  -5.194    7.555    1.268)    9.255
   1.292   (  -4.971    6.941   -0.158)    8.539
   3.037   (  11.057    8.492    8.766)   16.468
   3.182   (   9.525    7.127   -6.575)   13.593
   4.145   (  -2.172   -5.642    6.541)    8.907
   4.250   (  -2.555   -5.591   -4.529)    7.635
   4.414   (   0.153   -0.869    1.372)    1.632
   4.439   (  -0.450    0.526   -1.100)    1.300
   4.516   (  -1.618    1.763    0.368)    2.421
   4.537   (  -0.996    0.664   -1.728)    2.102
======================= Grid point 575 (80/84) =======================
q-point: ( 0.38  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.888   (   0.435   -1.066   -1.658)    2.019
   0.944   (  -0.786   -1.737   -7.063)    7.316
   1.257   (  -5.529    7.954    1.285)    9.772
   1.268   (  -5.374    7.280    0.166)    9.050
   3.248   (   5.404    4.034    7.290)    9.930
   3.355   (   4.242    2.626   -4.007)    6.399
   4.071   (   0.364   -5.068    5.277)    7.326
   4.151   (  -1.139   -5.957   -3.223)    6.868
   4.412   (   0.574   -0.939    0.608)    1.258
   4.433   (   0.065   -0.622   -1.504)    1.629
   4.513   (  -1.834    3.192    0.705)    3.749
   4.535   (  -1.610    2.858   -1.513)    3.613
======================= Grid point 586 (81/84) =======================
q-point: ( 0.23  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.897   (  -1.100   -1.906    1.699)    2.780
   0.995   (  -1.749   -3.029   -7.794)    8.543
   1.366   (   1.661    2.878    0.987)    3.466
   1.374   (   1.482    2.567    0.159)    2.968
   2.980   (   7.105   12.306    9.100)   16.874
   3.132   (   6.061   10.497   -6.982)   13.988
   4.136   (  -3.207   -5.555    6.778)    9.332
   4.245   (  -3.170   -5.490   -4.701)    7.892
   4.407   (  -0.589   -1.020    1.769)    2.126
   4.443   (  -0.468   -0.811   -1.896)    2.115
   4.539   (   0.511    0.885    0.073)    1.025
   4.551   (   0.798    1.381   -1.052)    1.911
======================= Grid point 587 (82/84) =======================
q-point: ( 0.31  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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|
   0.876   (   0.005   -0.596   -1.369)    1.493
   0.944   (  -1.320   -2.163   -7.833)    8.233
   1.383   (  -2.332    3.177    0.883)    4.039
   1.386   (  -2.278    2.662    0.548)    3.546
   3.204   (   6.579    8.543    7.550)   13.164
   3.317   (   5.090    6.452   -4.385)    9.315
   4.032   (  -1.998   -5.546    6.070)    8.462
   4.125   (  -3.165   -6.844   -3.637)    8.372
   4.396   (   0.209   -0.620    0.586)    0.878
   4.423   (  -0.308   -1.233   -2.143)    2.491
   4.560   (  -0.285    2.147    0.775)    2.300
   4.579   (  -0.127    2.388   -1.142)    2.650
======================= Grid point 588 (83/84) =======================
q-point: ( 0.38  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.878   (   0.030   -0.052   -3.300)    3.300
   0.918   (   0.551   -0.954   -7.119)    7.204
   1.375   (  -2.041    3.535    0.783)    4.157
   1.375   (  -2.302    3.988    0.758)    4.667
   3.317   (  -1.704    2.952    6.148)    7.029
   3.393   (  -0.752    1.303   -1.859)    2.392
   3.988   (   1.800   -3.117    4.105)    5.459
   4.050   (   2.396   -4.150   -2.345)    5.335
   4.396   (   0.305   -0.529   -0.086)    0.617
   4.414   (   0.545   -0.944   -1.871)    2.165
   4.568   (  -1.219    2.112    1.012)    2.640
   4.589   (  -1.287    2.229   -1.201)    2.841
======================= Grid point 601 (84/84) =======================
q-point: ( 0.31  0.31  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.72e-04 2.72e-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.879   (   0.441    0.763   -4.170)    4.263
   0.908   (  -0.864   -1.496   -6.883)    7.097
   1.408   (   0.093    0.162    1.068)    1.084
   1.413   (   0.228    0.395    0.551)    0.715
   3.350   (   3.270    5.664    5.264)    8.395
   3.405   (   0.985    1.706   -0.382)    2.007
   3.958   (  -0.776   -1.345    3.132)    3.496
   4.006   (  -2.744   -4.752   -1.739)    5.756
   4.392   (   0.083    0.143   -0.552)    0.576
   4.403   (  -0.432   -0.748   -1.690)    1.898
   4.589   (   0.447    0.774    1.148)    1.455
   4.612   (   0.521    0.902   -1.202)    1.590
=================== 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    152.531    152.531    188.367     -0.000      0.000      0.000 3/14196
   20.0    109.304    109.304    124.760      0.000      0.000      0.000 3/14196
   30.0     79.872     79.872     87.428      0.000      0.000      0.000 3/14196
   40.0     57.554     57.554     62.359      0.000      0.000      0.000 3/14196
   50.0     43.237     43.237     46.804      0.000      0.000      0.000 3/14196
   60.0     34.180     34.180     37.041      0.000      0.000      0.000 3/14196
   70.0     28.178     28.178     30.577      0.000      0.000      0.000 3/14196
   80.0     23.976     23.976     26.047      0.000      0.000      0.000 3/14196
   90.0     20.890     20.890     22.714      0.000      0.000      0.000 3/14196
  100.0     18.530     18.530     20.162      0.000      0.000      0.000 3/14196
  110.0     16.668     16.668     18.144      0.000      0.000      0.000 3/14196
  120.0     15.159     15.159     16.508      0.000      0.000      0.000 3/14196
  130.0     13.911     13.911     15.152      0.000      0.000      0.000 3/14196
  140.0     12.861     12.861     14.010      0.000      0.000      0.000 3/14196
  150.0     11.963     11.963     13.034      0.000      0.000      0.000 3/14196
  160.0     11.187     11.187     12.189      0.000      0.000      0.000 3/14196
  170.0     10.508     10.508     11.451      0.000      0.000      0.000 3/14196
  180.0      9.910      9.910     10.799      0.000      0.000      0.000 3/14196
  190.0      9.378      9.378     10.219      0.000      0.000      0.000 3/14196
  200.0      8.901      8.901      9.700      0.000      0.000      0.000 3/14196
  210.0      8.472      8.472      9.232      0.000      0.000      0.000 3/14196
  220.0      8.083      8.083      8.808      0.000      0.000      0.000 3/14196
  230.0      7.730      7.730      8.423      0.000      0.000      0.000 3/14196
  240.0      7.406      7.406      8.070      0.000      0.000      0.000 3/14196
  250.0      7.109      7.109      7.746      0.000      0.000      0.000 3/14196
  260.0      6.835      6.835      7.447      0.000      0.000      0.000 3/14196
  270.0      6.582      6.582      7.171      0.000      0.000      0.000 3/14196
  280.0      6.348      6.348      6.915      0.000      0.000      0.000 3/14196
  290.0      6.130      6.130      6.677      0.000      0.000      0.000 3/14196
  300.0      5.926      5.926      6.455      0.000      0.000      0.000 3/14196
  310.0      5.736      5.736      6.248      0.000      0.000      0.000 3/14196
  320.0      5.557      5.557      6.053      0.000      0.000      0.000 3/14196
  330.0      5.390      5.390      5.871      0.000      0.000      0.000 3/14196
  340.0      5.233      5.233      5.699      0.000      0.000      0.000 3/14196
  350.0      5.084      5.084      5.537      0.000      0.000      0.000 3/14196
  360.0      4.944      4.944      5.384      0.000      0.000      0.000 3/14196
  370.0      4.811      4.811      5.240      0.000      0.000      0.000 3/14196
  380.0      4.686      4.686      5.103      0.000      0.000      0.000 3/14196
  390.0      4.567      4.567      4.973      0.000      0.000      0.000 3/14196
  400.0      4.453      4.453      4.849      0.000      0.000      0.000 3/14196
  410.0      4.346      4.346      4.732      0.000      0.000      0.000 3/14196
  420.0      4.243      4.243      4.620      0.000      0.000      0.000 3/14196
  430.0      4.145      4.145      4.513      0.000      0.000      0.000 3/14196
  440.0      4.052      4.052      4.412      0.000      0.000      0.000 3/14196
  450.0      3.963      3.963      4.314      0.000      0.000      0.000 3/14196
  460.0      3.878      3.878      4.221      0.000      0.000      0.000 3/14196
  470.0      3.796      3.796      4.132      0.000      0.000      0.000 3/14196
  480.0      3.718      3.718      4.047      0.000      0.000      0.000 3/14196
  490.0      3.642      3.642      3.965      0.000      0.000      0.000 3/14196
  500.0      3.570      3.570      3.886      0.000      0.000      0.000 3/14196
  510.0      3.501      3.501      3.811      0.000      0.000      0.000 3/14196
  520.0      3.434      3.434      3.738      0.000      0.000      0.000 3/14196
  530.0      3.370      3.370      3.668      0.000      0.000      0.000 3/14196
  540.0      3.308      3.308      3.601      0.000      0.000      0.000 3/14196
  550.0      3.249      3.249      3.536      0.000      0.000      0.000 3/14196
  560.0      3.191      3.191      3.473      0.000      0.000      0.000 3/14196
  570.0      3.136      3.136      3.413      0.000      0.000      0.000 3/14196
  580.0      3.082      3.082      3.354      0.000      0.000      0.000 3/14196
  590.0      3.031      3.031      3.298      0.000      0.000      0.000 3/14196
  600.0      2.981      2.981      3.244      0.000      0.000      0.000 3/14196
  610.0      2.932      2.932      3.191      0.000      0.000      0.000 3/14196
  620.0      2.885      2.885      3.140      0.000      0.000      0.000 3/14196
  630.0      2.840      2.840      3.090      0.000      0.000      0.000 3/14196
  640.0      2.796      2.796      3.042      0.000      0.000      0.000 3/14196
  650.0      2.753      2.753      2.996      0.000      0.000      0.000 3/14196
  660.0      2.712      2.712      2.951      0.000      0.000      0.000 3/14196
  670.0      2.672      2.672      2.907      0.000      0.000      0.000 3/14196
  680.0      2.633      2.633      2.865      0.000      0.000      0.000 3/14196
  690.0      2.595      2.595      2.824      0.000      0.000      0.000 3/14196
  700.0      2.559      2.559      2.784      0.000      0.000      0.000 3/14196
  710.0      2.523      2.523      2.745      0.000      0.000      0.000 3/14196
  720.0      2.488      2.488      2.707      0.000      0.000      0.000 3/14196
  730.0      2.454      2.454      2.670      0.000      0.000      0.000 3/14196
  740.0      2.422      2.422      2.634      0.000      0.000      0.000 3/14196
  750.0      2.390      2.390      2.600      0.000      0.000      0.000 3/14196
  760.0      2.358      2.358      2.566      0.000      0.000      0.000 3/14196
  770.0      2.328      2.328      2.533      0.000      0.000      0.000 3/14196
  780.0      2.298      2.298      2.500      0.000      0.000      0.000 3/14196
  790.0      2.270      2.270      2.469      0.000      0.000      0.000 3/14196
  800.0      2.242      2.242      2.438      0.000      0.000      0.000 3/14196
  810.0      2.214      2.214      2.408      0.000      0.000      0.000 3/14196
  820.0      2.187      2.187      2.379      0.000      0.000      0.000 3/14196
  830.0      2.161      2.161      2.351      0.000      0.000      0.000 3/14196
  840.0      2.136      2.136      2.323      0.000      0.000      0.000 3/14196
  850.0      2.111      2.111      2.296      0.000      0.000      0.000 3/14196
  860.0      2.086      2.086      2.269      0.000      0.000      0.000 3/14196
  870.0      2.063      2.063      2.244      0.000      0.000      0.000 3/14196
  880.0      2.039      2.039      2.218      0.000      0.000      0.000 3/14196
  890.0      2.017      2.017      2.194      0.000      0.000      0.000 3/14196
  900.0      1.995      1.995      2.169      0.000      0.000      0.000 3/14196
  910.0      1.973      1.973      2.146      0.000      0.000      0.000 3/14196
  920.0      1.952      1.952      2.123      0.000      0.000      0.000 3/14196
  930.0      1.931      1.931      2.100      0.000      0.000      0.000 3/14196
  940.0      1.910      1.910      2.078      0.000      0.000      0.000 3/14196
  950.0      1.890      1.890      2.056      0.000      0.000      0.000 3/14196
  960.0      1.871      1.871      2.035      0.000      0.000      0.000 3/14196
  970.0      1.852      1.852      2.014      0.000      0.000      0.000 3/14196
  980.0      1.833      1.833      1.993      0.000      0.000      0.000 3/14196
  990.0      1.815      1.815      1.973      0.000      0.000      0.000 3/14196
 1000.0      1.797      1.797      1.954      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-08 00:10:03]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

