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

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

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

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

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

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


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

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

-------------------------[time 2026-01-07 18:42:54]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [4 4 2]
Primitive matrix:
  [ 1. -0.  0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.302581863194725    0.000000000000000    0.000000000000000
  b   -2.151290931597362    3.726145195388813    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.015388610000000
Atomic positions (fractional):
    1 Se  0.66666666666667  0.33333333333333  0.12418947776900  78.960
    2 Se  0.33333333333333  0.66666666666667  0.62418947776900  78.960
    3 Cd  0.66666666666667  0.33333333333333  0.49995052223101 112.411
    4 Cd  0.33333333333333  0.66666666666667  0.99995052223101 112.411
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   17.210327452778898    0.000000000000000    0.000000000000000
  b   -8.605163726389447   14.904580781555252    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.030777219999999
Atomic positions (fractional):
    1 Se  0.16666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    2 Se  0.41666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    3 Se  0.66666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    4 Se  0.91666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    5 Se  0.16666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    6 Se  0.41666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    7 Se  0.66666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    8 Se  0.91666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    9 Se  0.16666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   10 Se  0.41666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   11 Se  0.66666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   12 Se  0.91666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   13 Se  0.16666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   14 Se  0.41666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   15 Se  0.66666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   16 Se  0.91666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   17 Se  0.16666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   18 Se  0.41666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   19 Se  0.66666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   20 Se  0.91666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   21 Se  0.16666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   22 Se  0.41666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   23 Se  0.66666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   24 Se  0.91666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   25 Se  0.16666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   26 Se  0.41666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   27 Se  0.66666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   28 Se  0.91666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   29 Se  0.16666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   30 Se  0.41666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   31 Se  0.66666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   32 Se  0.91666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   33 Se  0.08333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   34 Se  0.33333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   35 Se  0.58333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   36 Se  0.83333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   37 Se  0.08333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   38 Se  0.33333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   39 Se  0.58333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   40 Se  0.83333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   41 Se  0.08333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   42 Se  0.33333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   43 Se  0.58333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   44 Se  0.83333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   45 Se  0.08333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   46 Se  0.33333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   47 Se  0.58333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   48 Se  0.83333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   49 Se  0.08333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   50 Se  0.33333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   51 Se  0.58333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   52 Se  0.83333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   53 Se  0.08333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   54 Se  0.33333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   55 Se  0.58333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   56 Se  0.83333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   57 Se  0.08333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   58 Se  0.33333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   59 Se  0.58333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   60 Se  0.83333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   61 Se  0.08333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   62 Se  0.33333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   63 Se  0.58333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   64 Se  0.83333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   65 Cd  0.16666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   66 Cd  0.41666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   67 Cd  0.66666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   68 Cd  0.91666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   69 Cd  0.16666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   70 Cd  0.41666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   71 Cd  0.66666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   72 Cd  0.91666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   73 Cd  0.16666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   74 Cd  0.41666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   75 Cd  0.66666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   76 Cd  0.91666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   77 Cd  0.16666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   78 Cd  0.41666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   79 Cd  0.66666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   80 Cd  0.91666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   81 Cd  0.16666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   82 Cd  0.41666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   83 Cd  0.66666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   84 Cd  0.91666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   85 Cd  0.16666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   86 Cd  0.41666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   87 Cd  0.66666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   88 Cd  0.91666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   89 Cd  0.16666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   90 Cd  0.41666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   91 Cd  0.66666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   92 Cd  0.91666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   93 Cd  0.16666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   94 Cd  0.41666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   95 Cd  0.66666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   96 Cd  0.91666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   97 Cd  0.08333333333333  0.16666666666667  0.49997526111550 112.411 > 97
   98 Cd  0.33333333333333  0.16666666666667  0.49997526111550 112.411 > 97
   99 Cd  0.58333333333333  0.16666666666667  0.49997526111550 112.411 > 97
  100 Cd  0.83333333333333  0.16666666666667  0.49997526111550 112.411 > 97
  101 Cd  0.08333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  102 Cd  0.33333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  103 Cd  0.58333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  104 Cd  0.83333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  105 Cd  0.08333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  106 Cd  0.33333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  107 Cd  0.58333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  108 Cd  0.83333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  109 Cd  0.08333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  110 Cd  0.33333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  111 Cd  0.58333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  112 Cd  0.83333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  113 Cd  0.08333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  114 Cd  0.33333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  115 Cd  0.58333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  116 Cd  0.83333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  117 Cd  0.08333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  118 Cd  0.33333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  119 Cd  0.58333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  120 Cd  0.83333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  121 Cd  0.08333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  122 Cd  0.33333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  123 Cd  0.58333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  124 Cd  0.83333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  125 Cd  0.08333333333333  0.91666666666667  0.99997526111550 112.411 > 97
  126 Cd  0.33333333333333  0.91666666666667  0.99997526111550 112.411 > 97
  127 Cd  0.58333333333333  0.91666666666667  0.99997526111550 112.411 > 97
  128 Cd  0.83333333333333  0.91666666666667  0.99997526111550 112.411 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           27.4945805    0.0000000    0.0000000
           -0.0000000   27.4945805    0.0000000
            0.0000000    0.0000000   24.1761992
-------------------------- Born effective charges --------------------------
    1 Se   -2.6471517   -0.0000000    0.0000000
           -0.0000000   -2.6471517    0.0000000
            0.0000000    0.0000000   -2.7000492
    2 Se   -2.6471517   -0.0000000    0.0000000
           -0.0000000   -2.6471517    0.0000000
            0.0000000    0.0000000   -2.7000492
    3 Cd    2.6471517    0.0000000    0.0000000
            0.0000000    2.6471517    0.0000000
            0.0000000    0.0000000    2.7000492
    4 Cd    2.6471517    0.0000000    0.0000000
            0.0000000    2.6471517    0.0000000
            0.0000000    0.0000000    2.7000492
----------------------------------------------------------------------------
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.00000064 (yzy) -0.00000064 (yzy) -0.00000064 (yyz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-07 18:43:00]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 18:43:00]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [4 4 2]
Primitive matrix:
  [ 1. -0.  0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.302581863194725    0.000000000000000    0.000000000000000
  b   -2.151290931597362    3.726145195388813    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.015388610000000
Atomic positions (fractional):
    1 Se  0.66666666666667  0.33333333333333  0.12418947776900  78.960
    2 Se  0.33333333333333  0.66666666666667  0.62418947776900  78.960
    3 Cd  0.66666666666667  0.33333333333333  0.49995052223101 112.411
    4 Cd  0.33333333333333  0.66666666666667  0.99995052223101 112.411
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   17.210327452778898    0.000000000000000    0.000000000000000
  b   -8.605163726389447   14.904580781555252    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.030777219999999
Atomic positions (fractional):
    1 Se  0.16666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    2 Se  0.41666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    3 Se  0.66666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    4 Se  0.91666666666667  0.08333333333333  0.06209473888450  78.960 > 1
    5 Se  0.16666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    6 Se  0.41666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    7 Se  0.66666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    8 Se  0.91666666666667  0.33333333333333  0.06209473888450  78.960 > 1
    9 Se  0.16666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   10 Se  0.41666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   11 Se  0.66666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   12 Se  0.91666666666667  0.58333333333333  0.06209473888450  78.960 > 1
   13 Se  0.16666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   14 Se  0.41666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   15 Se  0.66666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   16 Se  0.91666666666667  0.83333333333333  0.06209473888450  78.960 > 1
   17 Se  0.16666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   18 Se  0.41666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   19 Se  0.66666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   20 Se  0.91666666666667  0.08333333333333  0.56209473888450  78.960 > 1
   21 Se  0.16666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   22 Se  0.41666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   23 Se  0.66666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   24 Se  0.91666666666667  0.33333333333333  0.56209473888450  78.960 > 1
   25 Se  0.16666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   26 Se  0.41666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   27 Se  0.66666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   28 Se  0.91666666666667  0.58333333333333  0.56209473888450  78.960 > 1
   29 Se  0.16666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   30 Se  0.41666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   31 Se  0.66666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   32 Se  0.91666666666667  0.83333333333333  0.56209473888450  78.960 > 1
   33 Se  0.08333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   34 Se  0.33333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   35 Se  0.58333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   36 Se  0.83333333333333  0.16666666666667  0.31209473888450  78.960 > 33
   37 Se  0.08333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   38 Se  0.33333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   39 Se  0.58333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   40 Se  0.83333333333333  0.41666666666667  0.31209473888450  78.960 > 33
   41 Se  0.08333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   42 Se  0.33333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   43 Se  0.58333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   44 Se  0.83333333333333  0.66666666666667  0.31209473888450  78.960 > 33
   45 Se  0.08333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   46 Se  0.33333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   47 Se  0.58333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   48 Se  0.83333333333333  0.91666666666667  0.31209473888450  78.960 > 33
   49 Se  0.08333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   50 Se  0.33333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   51 Se  0.58333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   52 Se  0.83333333333333  0.16666666666667  0.81209473888450  78.960 > 33
   53 Se  0.08333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   54 Se  0.33333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   55 Se  0.58333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   56 Se  0.83333333333333  0.41666666666667  0.81209473888450  78.960 > 33
   57 Se  0.08333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   58 Se  0.33333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   59 Se  0.58333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   60 Se  0.83333333333333  0.66666666666667  0.81209473888450  78.960 > 33
   61 Se  0.08333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   62 Se  0.33333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   63 Se  0.58333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   64 Se  0.83333333333333  0.91666666666667  0.81209473888450  78.960 > 33
   65 Cd  0.16666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   66 Cd  0.41666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   67 Cd  0.66666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   68 Cd  0.91666666666667  0.08333333333333  0.24997526111550 112.411 > 65
   69 Cd  0.16666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   70 Cd  0.41666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   71 Cd  0.66666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   72 Cd  0.91666666666667  0.33333333333333  0.24997526111550 112.411 > 65
   73 Cd  0.16666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   74 Cd  0.41666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   75 Cd  0.66666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   76 Cd  0.91666666666667  0.58333333333333  0.24997526111550 112.411 > 65
   77 Cd  0.16666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   78 Cd  0.41666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   79 Cd  0.66666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   80 Cd  0.91666666666667  0.83333333333333  0.24997526111550 112.411 > 65
   81 Cd  0.16666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   82 Cd  0.41666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   83 Cd  0.66666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   84 Cd  0.91666666666667  0.08333333333333  0.74997526111550 112.411 > 65
   85 Cd  0.16666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   86 Cd  0.41666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   87 Cd  0.66666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   88 Cd  0.91666666666667  0.33333333333333  0.74997526111550 112.411 > 65
   89 Cd  0.16666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   90 Cd  0.41666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   91 Cd  0.66666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   92 Cd  0.91666666666667  0.58333333333333  0.74997526111550 112.411 > 65
   93 Cd  0.16666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   94 Cd  0.41666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   95 Cd  0.66666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   96 Cd  0.91666666666667  0.83333333333333  0.74997526111550 112.411 > 65
   97 Cd  0.08333333333333  0.16666666666667  0.49997526111550 112.411 > 97
   98 Cd  0.33333333333333  0.16666666666667  0.49997526111550 112.411 > 97
   99 Cd  0.58333333333333  0.16666666666667  0.49997526111550 112.411 > 97
  100 Cd  0.83333333333333  0.16666666666667  0.49997526111550 112.411 > 97
  101 Cd  0.08333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  102 Cd  0.33333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  103 Cd  0.58333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  104 Cd  0.83333333333333  0.41666666666667  0.49997526111550 112.411 > 97
  105 Cd  0.08333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  106 Cd  0.33333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  107 Cd  0.58333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  108 Cd  0.83333333333333  0.66666666666667  0.49997526111550 112.411 > 97
  109 Cd  0.08333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  110 Cd  0.33333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  111 Cd  0.58333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  112 Cd  0.83333333333333  0.91666666666667  0.49997526111550 112.411 > 97
  113 Cd  0.08333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  114 Cd  0.33333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  115 Cd  0.58333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  116 Cd  0.83333333333333  0.16666666666667  0.99997526111550 112.411 > 97
  117 Cd  0.08333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  118 Cd  0.33333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  119 Cd  0.58333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  120 Cd  0.83333333333333  0.41666666666667  0.99997526111550 112.411 > 97
  121 Cd  0.08333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  122 Cd  0.33333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  123 Cd  0.58333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  124 Cd  0.83333333333333  0.66666666666667  0.99997526111550 112.411 > 97
  125 Cd  0.08333333333333  0.91666666666667  0.99997526111550 112.411 > 97
  126 Cd  0.33333333333333  0.91666666666667  0.99997526111550 112.411 > 97
  127 Cd  0.58333333333333  0.91666666666667  0.99997526111550 112.411 > 97
  128 Cd  0.83333333333333  0.91666666666667  0.99997526111550 112.411 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           27.4945805    0.0000000    0.0000000
           -0.0000000   27.4945805    0.0000000
            0.0000000    0.0000000   24.1761992
-------------------------- Born effective charges --------------------------
    1 Se   -2.6471517   -0.0000000    0.0000000
           -0.0000000   -2.6471517    0.0000000
            0.0000000    0.0000000   -2.7000492
    2 Se   -2.6471517   -0.0000000    0.0000000
           -0.0000000   -2.6471517    0.0000000
            0.0000000    0.0000000   -2.7000492
    3 Cd    2.6471517    0.0000000    0.0000000
            0.0000000    2.6471517    0.0000000
            0.0000000    0.0000000    2.7000492
    4 Cd    2.6471517    0.0000000    0.0000000
            0.0000000    2.6471517    0.0000000
            0.0000000    0.0000000    2.7000492
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000064 (yzy) -0.00000064 (yzy) -0.00000064 (yyz)
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 13 13 7 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.86, Number of G-points: 313, Lambda: 0.46
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: 4.63e-04 4.63e-04 2.72e-04 2.72e-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.959   (   0.000    0.000   -0.000)    0.000
   0.959   (  -0.000   -0.000   -0.000)    0.000
   3.651   (   0.000    0.000    0.000)    0.000
   5.019   (   0.000    0.000   -0.000)    0.000
   5.142   (  -0.000    0.000   -0.000)    0.000
   5.142   (  -0.000    0.000    0.000)    0.000
   5.224   (  -0.000    0.000   -0.000)    0.000
   5.224   (  -0.000    0.000    0.000)    0.000
   5.589   (   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: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.276   (  12.103    6.988    0.000)   13.975
   0.279   (  11.448    6.609    0.000)   13.219
   0.754   (  30.738   17.746    0.000)   35.493
   0.973   (   1.427    0.824    0.000)    1.648
   1.169   (  14.854    8.576    0.000)   17.151
   3.617   (  -2.433   -1.405    0.000)    2.809
   5.029   (   0.811    0.468    0.000)    0.937
   5.151   (   0.739    0.426    0.000)    0.853
   5.246   (   1.653    0.954    0.000)    1.908
   5.247   (   1.478    0.853    0.000)    1.706
   5.599   (   1.302    0.752    0.000)    1.503
   5.635   (   6.585    3.802    0.000)    7.604
======================= Grid point 2 (3/84) =======================
q-point: ( 0.15  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.538   (  10.035    5.794    0.000)   11.587
   0.569   (  12.404    7.162    0.000)   14.323
   1.040   (   4.226    2.440    0.000)    4.879
   1.426   (  22.102   12.761    0.000)   25.521
   1.553   (  18.601   10.739    0.000)   21.478
   3.569   (  -0.680   -0.393    0.000)    0.785
   5.059   (   1.762    1.017    0.000)    2.035
   5.176   (   1.272    0.734    0.000)    1.468
   5.265   (  -0.494   -0.285    0.000)    0.571
   5.294   (   2.258    1.304    0.000)    2.607
   5.654   (   3.289    1.899    0.000)    3.798
   5.822   (   7.994    4.615    0.000)    9.230
======================= Grid point 3 (4/84) =======================
q-point: ( 0.23  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.746   (   7.136    4.120    0.000)    8.240
   0.859   (  11.752    6.785    0.000)   13.570
   1.157   (   4.976    2.873    0.000)    5.745
   1.682   (   3.810    2.200    0.000)    4.399
   2.120   (  24.244   13.997    0.000)   27.994
   3.614   (   4.485    2.589    0.000)    5.179
   5.112   (   2.662    1.537    0.000)    3.074
   5.207   (   1.278    0.738    0.000)    1.476
   5.222   (  -2.865   -1.654    0.000)    3.308
   5.351   (   2.472    1.427    0.000)    2.855
   5.736   (   3.022    1.745    0.000)    3.490
   5.969   (   3.637    2.100    0.000)    4.200
======================= Grid point 4 (5/84) =======================
q-point: ( 0.31  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.872   (   3.469    2.003    0.000)    4.005
   1.120   (   9.974    5.759    0.000)   11.517
   1.247   (   2.080    1.201    0.000)    2.401
   1.694   (  -1.661   -0.959    0.000)    1.918
   2.653   (  20.310   11.726    0.000)   23.452
   3.760   (   6.990    4.036    0.000)    8.071
   5.138   (  -3.437   -1.984    0.000)    3.968
   5.185   (   2.970    1.715    0.000)    3.430
   5.238   (   1.375    0.794    0.000)    1.587
   5.410   (   2.429    1.402    0.000)    2.804
   5.778   (   0.158    0.091    0.000)    0.183
   5.980   (  -2.562   -1.479    0.000)    2.958
======================= Grid point 5 (6/84) =======================
q-point: ( 0.38  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.921   (   1.049    0.606    0.000)    1.212
   1.248   (  -1.695   -0.979    0.000)    1.957
   1.327   (   7.207    4.161    0.000)    8.321
   1.644   (  -2.044   -1.180    0.000)    2.360
   3.081   (  15.295    8.831    0.000)   17.662
   3.908   (   4.662    2.692    0.000)    5.384
   5.070   (  -1.609   -0.929    0.000)    1.858
   5.253   (   2.292    1.323    0.000)    2.646
   5.272   (   1.454    0.840    0.000)    1.679
   5.462   (   1.832    1.058    0.000)    2.115
   5.728   (  -4.618   -2.666    0.000)    5.332
   5.886   (  -3.944   -2.277    0.000)    4.554
======================= Grid point 6 (7/84) =======================
q-point: ( 0.46  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.937   (   0.379    0.219    0.000)    0.438
   1.196   (  -1.801   -1.040    0.000)    2.080
   1.454   (   3.163    1.826    0.000)    3.653
   1.604   (  -1.148   -0.663    0.000)    1.325
   3.355   (   6.701    3.869    0.000)    7.738
   3.966   (   0.546    0.315    0.000)    0.630
   5.076   (   1.541    0.890    0.000)    1.780
   5.276   (  -0.474   -0.274    0.000)    0.547
   5.300   (   0.685    0.395    0.000)    0.791
   5.493   (   0.683    0.394    0.000)    0.788
   5.603   (  -3.821   -2.206    0.000)    4.412
   5.826   (  -1.230   -0.710    0.000)    1.420
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.489   (   7.082   12.266    0.000)   14.163
   0.524   (   7.810   13.527    0.000)   15.619
   1.031   (   3.664    6.347    0.000)    7.328
   1.282   (  15.600   27.020    0.000)   31.200
   1.402   (   7.603   13.168    0.000)   15.205
   3.576   (  -1.010   -1.750    0.000)    2.020
   5.049   (   0.857    1.485    0.000)    1.714
   5.187   (   1.448    2.508    0.000)    2.896
   5.264   (  -0.189   -0.327    0.000)    0.378
   5.272   (   1.374    2.381    0.000)    2.749
   5.639   (   1.667    2.888    0.000)    3.335
   5.758   (   4.304    7.454    0.000)    8.607
======================= Grid point 15 (9/84) =======================
q-point: ( 0.15  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.753   (   9.137    9.988    0.000)   13.537
   0.777   (   3.771   14.614    0.000)   15.093
   1.171   (   2.809    9.321    0.000)    9.735
   1.623   (   5.725    5.677    0.000)    8.063
   1.866   (  21.048   18.843    0.000)   28.250
   3.579   (   1.963    1.754    0.000)    2.633
   5.090   (   1.785    1.990    0.000)    2.673
   5.216   (  -1.438   -1.232    0.000)    1.894
   5.248   (  -0.345    2.612    0.000)    2.635
   5.338   (   2.443    2.993    0.000)    3.864
   5.709   (   2.316    2.872    0.000)    3.689
   5.903   (   4.995    2.661    0.000)    5.660
======================= Grid point 16 (10/84) =======================
q-point: ( 0.23  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.950   (  -0.693   13.464    0.000)   13.481
   1.011   (   9.237    7.777    0.000)   12.075
   1.302   (  -0.052    9.908    0.000)    9.908
   1.702   (   0.360    0.139    0.000)    0.386
   2.411   (  20.554   13.972    0.000)   24.854
   3.684   (   5.796    4.105    0.000)    7.102
   5.139   (  -1.338   -3.489    0.000)    3.736
   5.156   (   1.807    0.777    0.000)    1.967
   5.280   (  -1.362    5.146    0.000)    5.323
   5.411   (   1.546    3.940    0.000)    4.233
   5.764   (   1.356    0.692    0.000)    1.523
   5.963   (   1.450   -2.683    0.000)    3.049
======================= Grid point 17 (11/84) =======================
q-point: ( 0.31  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.033   (  -4.275   11.703    0.000)   12.460
   1.235   (   7.654    5.717    0.000)    9.554
   1.356   (  -4.141    7.907    0.000)    8.926
   1.674   (  -2.261   -0.657    0.000)    2.354
   2.879   (  17.631   10.149    0.000)   20.343
   3.836   (   6.048    2.953    0.000)    6.730
   5.057   (  -1.055   -4.920    0.000)    5.032
   5.224   (   2.637    1.865    0.000)    3.230
   5.305   (  -1.460    4.643    0.000)    4.868
   5.472   (   0.636    3.930    0.000)    3.981
   5.757   (  -1.715   -2.132    0.000)    2.737
   5.910   (  -2.199   -4.549    0.000)    5.053
======================= Grid point 18 (12/84) =======================
q-point: ( 0.38  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.055   (  -5.614   10.557    0.000)   11.957
   1.314   (  -6.450    6.448    0.000)    9.120
   1.411   (   4.641    3.919    0.000)    6.074
   1.627   (  -2.140   -0.534    0.000)    2.206
   3.235   (  11.995    6.129    0.000)   13.470
   3.930   (   2.804   -0.816    0.000)    2.921
   5.025   (   2.264   -2.695    0.000)    3.520
   5.267   (   0.548   -0.340    0.000)    0.645
   5.325   (  -0.924    3.391    0.000)    3.514
   5.514   (  -0.023    3.328    0.000)    3.328
   5.658   (  -4.956   -3.161    0.000)    5.878
   5.839   (  -1.512   -1.933    0.000)    2.454
======================= Grid point 19 (13/84) =======================
q-point: ( 0.46  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.059   (  -5.810   10.062    0.000)   11.619
   1.279   (  -4.553    7.886    0.000)    9.106
   1.484   (  -0.589    1.021    0.000)    1.179
   1.597   (  -0.050    0.087    0.000)    0.100
   3.387   (   0.380   -0.658    0.000)    0.760
   3.940   (   1.520   -2.633    0.000)    3.040
   5.051   (   1.710   -2.962    0.000)    3.421
   5.250   (   0.843   -1.460    0.000)    1.685
   5.336   (  -1.435    2.485    0.000)    2.870
   5.541   (  -2.046    3.544    0.000)    4.093
   5.580   (  -0.063    0.109    0.000)    0.126
   5.815   (   0.131   -0.226    0.000)    0.261
======================= Grid point 29 (14/84) =======================
q-point: ( 0.15  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.964   (   5.930   10.272    0.000)   11.861
   1.054   (   6.413   11.108    0.000)   12.826
   1.370   (   5.055    8.755    0.000)   10.110
   1.696   (   1.008    1.745    0.000)    2.015
   2.278   (  11.905   20.621    0.000)   23.811
   3.651   (   2.984    5.168    0.000)    5.968
   5.135   (  -2.867   -4.966    0.000)    5.735
   5.140   (   1.621    2.807    0.000)    3.241
   5.340   (   2.386    4.133    0.000)    4.773
   5.410   (   2.187    3.788    0.000)    4.374
   5.756   (   0.857    1.485    0.000)    1.714
   5.922   (  -0.476   -0.824    0.000)    0.951
======================= Grid point 30 (15/84) =======================
q-point: ( 0.23  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.166   (   5.571    7.774    0.000)    9.564
   1.224   (   0.494   10.977    0.000)   10.988
   1.498   (  -0.677    7.706    0.000)    7.736
   1.704   (  -1.075    0.317    0.000)    1.120
   2.701   (  14.695   14.348    0.000)   20.538
   3.774   (   4.503    4.249    0.000)    6.191
   5.031   (  -2.274   -5.345    0.000)    5.808
   5.203   (   2.107    2.506    0.000)    3.274
   5.395   (  -1.768    5.462    0.000)    5.741
   5.485   (   2.260    2.405    0.000)    3.300
   5.763   (  -0.852   -0.701    0.000)    1.103
   5.879   (  -0.585   -4.615    0.000)    4.652
======================= Grid point 31 (16/84) =======================
q-point: ( 0.31  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.265   (  -3.665    9.264    0.000)    9.963
   1.365   (   4.712    6.774    0.000)    8.252
   1.507   (  -5.385    5.372    0.000)    7.607
   1.671   (  -2.780    0.339    0.000)    2.800
   3.081   (  13.283    9.726    0.000)   16.463
   3.869   (   2.931   -0.057    0.000)    2.932
   4.965   (   0.532   -3.340    0.000)    3.383
   5.253   (   0.983    0.473    0.000)    1.091
   5.398   (  -2.832    3.805    0.000)    4.743
   5.533   (   1.183    1.219    0.000)    1.698
   5.708   (  -3.421   -1.623    0.000)    3.786
   5.833   (  -0.270   -2.170    0.000)    2.187
======================= Grid point 32 (17/84) =======================
q-point: ( 0.38  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.276   (  -5.371    9.456    0.000)   10.875
   1.452   (  -5.479    5.206    0.000)    7.558
   1.495   (   2.580    4.438    0.000)    5.133
   1.620   (  -2.751   -0.073    0.000)    2.752
   3.343   (   6.196    4.677    0.000)    7.763
   3.875   (   1.203   -4.243    0.000)    4.410
   4.982   (   2.883   -1.241    0.000)    3.139
   5.228   (  -1.027   -3.280    0.000)    3.437
   5.395   (  -1.657    3.330    0.000)    3.720
   5.569   (   1.351    1.850    0.000)    2.291
   5.632   (  -3.788    0.099    0.000)    3.790
   5.819   (  -0.486    0.296    0.000)    0.569
======================= Grid point 43 (18/84) =======================
q-point: ( 0.23  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.318   (   3.901    6.757    0.000)    7.803
   1.388   (   2.988    5.175    0.000)    5.976
   1.587   (   0.480    0.831    0.000)    0.960
   1.707   (  -0.294   -0.509    0.000)    0.587
   3.004   (   8.515   14.748    0.000)   17.029
   3.840   (   0.929    1.610    0.000)    1.859
   4.943   (  -1.504   -2.606    0.000)    3.009
   5.246   (   0.628    1.087    0.000)    1.255
   5.483   (  -0.812   -1.406    0.000)    1.623
   5.492   (   1.700    2.945    0.000)    3.400
   5.737   (  -1.196   -2.072    0.000)    2.393
   5.817   (  -0.349   -0.605    0.000)    0.699
======================= Grid point 44 (19/84) =======================
q-point: ( 0.31  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.415   (  -1.154    5.993    0.000)    6.103
   1.481   (   3.353    3.690    0.000)    4.986
   1.557   (  -3.018   -0.046    0.000)    3.018
   1.674   (  -2.585   -0.695    0.000)    2.677
   3.288   (   7.898   10.363    0.000)   13.029
   3.826   (  -0.885   -4.100    0.000)    4.194
   4.937   (   1.708    0.764    0.000)    1.871
   5.217   (  -2.041   -4.058    0.000)    4.542
   5.455   (  -1.257    2.422    0.000)    2.729
   5.551   (   2.094    1.358    0.000)    2.496
   5.678   (  -2.531   -2.199    0.000)    3.353
   5.833   (   0.373    1.687    0.000)    1.727
======================= Grid point 45 (20/84) =======================
q-point: ( 0.38  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.431   (  -3.440    5.959    0.000)    6.880
   1.507   (  -0.527    0.912    0.000)    1.053
   1.571   (  -0.925    1.602    0.000)    1.850
   1.621   (  -0.130    0.226    0.000)    0.261
   3.437   (  -2.541    4.401    0.000)    5.082
   3.773   (   2.885   -4.997    0.000)    5.770
   4.977   (  -0.374    0.649    0.000)    0.749
   5.154   (   1.938   -3.357    0.000)    3.876
   5.464   (  -1.936    3.353    0.000)    3.872
   5.600   (  -0.348    0.603    0.000)    0.697
   5.618   (   0.654   -1.133    0.000)    1.308
   5.846   (  -1.008    1.746    0.000)    2.016
======================= Grid point 58 (21/84) =======================
q-point: ( 0.31  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.507   (   0.476    0.824    0.000)    0.951
   1.529   (   1.312    2.273    0.000)    2.625
   1.550   (   0.838    1.451    0.000)    1.675
   1.630   (  -1.944   -3.368    0.000)    3.889
   3.496   (   5.228    9.056    0.000)   10.457
   3.707   (  -3.939   -6.822    0.000)    7.878
   4.990   (   2.227    3.857    0.000)    4.453
   5.112   (  -3.001   -5.199    0.000)    6.003
   5.515   (   1.681    2.911    0.000)    3.361
   5.575   (   0.256    0.443    0.000)    0.512
   5.621   (  -1.513   -2.621    0.000)    3.026
   5.867   (   0.623    1.080    0.000)    1.247
======================= Grid point 181 (22/84) =======================
q-point: ( 0.00  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.301   (   0.000   -0.000   14.137)   14.137
   0.301   (   0.000   -0.000   14.137)   14.137
   0.819   (  -0.000   -0.000   39.092)   39.092
   0.958   (   0.000    0.000   -0.220)    0.220
   0.958   (   0.000   -0.000   -0.220)    0.220
   3.521   (   0.000    0.000  -12.020)   12.020
   5.145   (   0.000    0.000    0.296)    0.296
   5.145   (   0.000   -0.000    0.296)    0.296
   5.219   (   0.000   -0.000   -0.475)    0.475
   5.219   (   0.000    0.000   -0.475)    0.475
   5.569   (   0.000   -0.000    5.852)    5.852
   5.668   (   0.000   -0.000    6.832)    6.832
======================= Grid point 182 (23/84) =======================
q-point: ( 0.08  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.410   (   7.795    4.501   10.338)   13.708
   0.584   (   7.956    4.593   16.838)   19.182
   0.976   (   1.705    0.984    0.177)    1.977
   0.994   (  19.722   11.386   23.429)   32.673
   1.160   (  14.672    8.471   -0.209)   16.943
   3.494   (  -1.855   -1.071  -11.416)   11.615
   5.148   (  -0.091   -0.053    6.549)    6.550
   5.154   (   0.752    0.434    0.311)    0.922
   5.240   (   1.592    0.919   -0.556)    1.921
   5.251   (   2.007    1.159    0.162)    2.323
   5.614   (   5.126    2.960    0.937)    5.993
   5.669   (   0.356    0.206    5.712)    5.726
======================= Grid point 183 (24/84) =======================
q-point: ( 0.15  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.617   (   8.779    5.069    6.914)   12.271
   0.735   (   6.724    3.882   14.055)   16.056
   1.047   (   4.227    2.440    0.554)    4.912
   1.459   (  12.990    7.500   -0.267)   15.002
   1.622   (  24.926   14.391    9.944)   30.451
   3.467   (   0.404    0.233   -9.615)    9.626
   5.148   (   0.395    0.228    7.292)    7.306
   5.180   (   1.328    0.767    0.389)    1.582
   5.286   (   0.393    0.227    1.584)    1.648
   5.287   (   2.206    1.274   -0.695)    2.641
   5.684   (   1.789    1.033    2.871)    3.537
   5.790   (   7.535    4.350   -2.473)    9.045
======================= Grid point 184 (25/84) =======================
q-point: ( 0.23  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.803   (   6.524    3.766    5.212)    9.160
   0.913   (   7.941    4.585    5.351)   10.617
   1.161   (   4.787    2.764    0.307)    5.537
   1.640   (   2.795    1.613   -3.926)    5.082
   2.215   (  23.811   13.747    8.979)   28.924
   3.539   (   5.714    3.299   -6.975)    9.601
   5.171   (   1.448    0.836    5.111)    5.378
   5.213   (   1.377    0.795    0.569)    1.689
   5.261   (  -2.253   -1.301    3.295)    4.198
   5.343   (   2.450    1.414   -0.778)    2.934
   5.736   (   2.061    1.190    0.139)    2.384
   5.929   (   3.463    1.999   -3.581)    5.367
======================= Grid point 185 (26/84) =======================
q-point: ( 0.31  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.919   (   3.187    1.840    4.317)    5.673
   1.100   (   7.534    4.350   -1.348)    8.803
   1.248   (   2.023    1.168    0.002)    2.336
   1.634   (  -2.284   -1.319   -5.563)    6.156
   2.741   (  20.108   11.610    8.085)   24.587
   3.713   (   8.061    4.654   -4.432)   10.310
   5.177   (  -2.634   -1.521    3.462)    4.608
   5.232   (   1.728    0.998    4.109)    4.568
   5.246   (   1.455    0.840    0.787)    1.856
   5.402   (   2.407    1.390   -0.829)    2.901
   5.760   (  -0.384   -0.222   -1.556)    1.618
   5.940   (  -2.446   -1.412   -3.750)    4.694
======================= Grid point 186 (27/84) =======================
q-point: ( 0.38  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.965   (   0.946    0.546    3.938)    4.087
   1.251   (  -1.500   -0.866    0.169)    1.740
   1.263   (   5.904    3.408   -5.581)    8.810
   1.570   (  -2.615   -1.510   -6.849)    7.485
   3.165   (  15.174    8.761    7.669)   19.126
   3.887   (   5.790    3.343   -2.066)    6.997
   5.115   (  -1.726   -0.996    4.052)    4.516
   5.282   (   1.466    0.846    0.893)    1.914
   5.284   (   1.929    1.114    2.785)    3.566
   5.453   (   1.808    1.044   -0.884)    2.268
   5.701   (  -4.820   -2.783   -2.485)    6.095
   5.851   (  -3.734   -2.156   -3.345)    5.457
======================= Grid point 187 (28/84) =======================
q-point: ( 0.46  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.978   (   0.325    0.188    3.754)    3.773
   1.204   (  -1.621   -0.936    0.678)    1.991
   1.371   (   2.847    1.644   -7.465)    8.156
   1.519   (  -1.504   -0.868   -7.907)    8.096
   3.436   (   6.597    3.809    7.350)   10.585
   3.969   (   1.259    0.727    0.129)    1.460
   5.113   (   1.197    0.691    3.450)    3.716
   5.301   (  -0.563   -0.325    2.420)    2.506
   5.309   (   0.673    0.388    0.880)    1.174
   5.483   (   0.672    0.388   -0.925)    1.207
   5.576   (  -3.690   -2.130   -2.483)    4.931
   5.792   (  -1.223   -0.706   -3.214)    3.511
======================= Grid point 195 (29/84) =======================
q-point: ( 0.08  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.604   (   6.761   11.711    7.055)   15.252
   0.686   (   3.468    6.006   15.801)   17.256
   1.041   (   3.218    5.573    0.630)    6.466
   1.383   (   7.171   12.420   -1.752)   14.449
   1.424   (  14.518   25.146   13.734)   32.120
   3.467   (  -0.462   -0.800  -10.220)   10.261
   5.145   (   0.076    0.131    7.502)    7.504
   5.188   (   1.336    2.314    0.155)    2.677
   5.267   (   1.307    2.264   -0.447)    2.652
   5.281   (   0.443    0.768    1.128)    1.435
   5.680   (   0.733    1.270    3.813)    4.086
   5.728   (   4.075    7.058   -1.957)    8.381
======================= Grid point 196 (30/84) =======================
q-point: ( 0.15  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.818   (   3.740   10.427    5.734)   12.474
   0.850   (   5.810    8.939    7.190)   12.859
   1.176   (   2.708    9.178    0.344)    9.575
   1.589   (   5.139    4.436   -3.272)    7.537
   1.972   (  20.025   18.807   10.098)   29.269
   3.493   (   2.903    2.763   -8.010)    8.957
   5.159   (   0.886    0.935    5.900)    6.039
   5.225   (  -0.375    1.455    0.542)    1.597
   5.275   (  -0.812    0.377    2.344)    2.509
   5.331   (   2.408    2.909   -0.618)    3.827
   5.718   (   1.433    1.859    1.004)    2.553
   5.867   (   4.747    2.600   -3.128)    6.251
======================= Grid point 197 (31/84) =======================
q-point: ( 0.23  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.976   (   1.619    9.829    2.033)   10.166
   1.030   (   4.689    7.897    2.570)    9.537
   1.303   (  -0.169    9.763    0.018)    9.765
   1.649   (  -0.209   -0.375   -4.930)    4.949
   2.504   (  20.060   13.980    8.639)   25.932
   3.626   (   6.785    4.909   -5.498)   10.018
   5.175   (  -1.401   -3.291    3.177)    4.784
   5.202   (   1.578    0.401    3.910)    4.236
   5.291   (  -1.250    4.761    1.153)    5.056
   5.403   (   1.517    3.882   -0.764)    4.237
   5.752   (   0.789    0.184   -0.966)    1.260
   5.925   (   1.329   -2.464   -3.536)    4.510
======================= Grid point 198 (32/84) =======================
q-point: ( 0.31  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.062   (  -3.264   10.094    2.330)   10.862
   1.192   (   5.727    4.986   -3.319)    8.287
   1.355   (  -4.250    7.895   -0.148)    8.967
   1.607   (  -2.834   -0.980   -6.189)    6.877
   2.966   (  17.376   10.210    7.989)   21.679
   3.803   (   7.064    3.743   -3.131)    8.586
   5.100   (  -0.978   -4.837    3.902)    6.291
   5.254   (   2.230    1.363    2.648)    3.721
   5.316   (  -1.312    4.470    1.155)    4.800
   5.462   (   0.616    3.831   -0.889)    3.981
   5.734   (  -2.039   -2.283   -2.094)    3.709
   5.874   (  -2.134   -4.160   -3.318)    5.733
======================= Grid point 199 (33/84) =======================
q-point: ( 0.38  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.085   (  -4.933    9.399    2.543)   10.915
   1.310   (  -2.600    5.652   -1.937)    6.516
   1.342   (   0.191    4.542   -4.330)    6.278
   1.546   (  -2.665   -0.862   -7.431)    7.942
   3.318   (  11.788    6.145    7.594)   15.310
   3.925   (   3.782    0.093   -0.635)    3.836
   5.064   (   1.947   -3.075    3.645)    5.150
   5.288   (   0.381   -0.603    1.911)    2.040
   5.335   (  -0.789    3.040    1.108)    3.330
   5.503   (   0.001    3.147   -1.047)    3.317
   5.634   (  -4.955   -2.683   -2.194)    6.047
   5.806   (  -1.575   -1.698   -3.038)    3.820
======================= Grid point 200 (34/84) =======================
q-point: ( 0.46  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.089   (  -5.207    9.018    2.554)   10.722
   1.287   (  -4.523    7.835    0.664)    9.071
   1.400   (  -0.716    1.240   -7.392)    7.529
   1.509   (   0.000   -0.000   -8.168)    8.168
   3.468   (   0.335   -0.581    7.354)    7.384
   3.951   (   1.346   -2.330    0.810)    2.810
   5.083   (   1.816   -3.146    3.030)    4.730
   5.269   (   1.001   -1.734    1.659)    2.600
   5.346   (  -1.269    2.199    0.971)    2.718
   5.531   (  -1.961    3.397   -1.005)    4.049
   5.559   (  -0.329    0.570   -1.876)    1.988
   5.783   (  -0.001    0.002   -3.008)    3.008
======================= Grid point 210 (35/84) =======================
q-point: ( 0.15  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.973   (   4.090    7.084    1.343)    8.289
   1.089   (   5.971   10.341    3.237)   12.372
   1.372   (   4.993    8.649    0.149)    9.988
   1.644   (   0.647    1.121   -4.782)    4.954
   2.379   (  11.693   20.253    9.273)   25.158
   3.588   (   3.609    6.251   -5.976)    9.371
   5.172   (  -2.555   -4.425    3.116)    5.984
   5.188   (   1.093    1.894    4.251)    4.781
   5.348   (   2.303    3.988    0.720)    4.661
   5.401   (   2.123    3.678   -0.859)    4.333
   5.747   (   0.448    0.775   -0.723)    1.151
   5.886   (  -0.414   -0.717   -3.225)    3.329
======================= Grid point 211 (36/84) =======================
q-point: ( 0.23  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.118   (   4.458    5.859   -3.812)    8.291
   1.247   (   0.364    9.845    1.957)   10.044
   1.499   (  -0.777    7.846    0.002)    7.885
   1.643   (  -1.384    0.047   -5.580)    5.749
   2.795   (  14.374   14.336    8.554)   22.030
   3.734   (   5.321    5.245   -3.832)    8.397
   5.075   (  -2.217   -5.176    3.924)    6.863
   5.234   (   1.719    1.900    2.835)    3.821
   5.401   (  -1.744    5.260    0.524)    5.567
   5.475   (   2.313    2.249   -0.936)    3.359
   5.742   (  -1.079   -1.016   -1.835)    2.359
   5.847   (  -0.533   -4.126   -2.895)    5.068
======================= Grid point 212 (37/84) =======================
q-point: ( 0.31  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.233   (   0.840    5.686   -4.909)    7.559
   1.334   (   0.151    8.117   -0.841)    8.161
   1.510   (  -5.322    5.810    0.239)    7.882
   1.599   (  -3.479    0.207   -6.410)    7.296
   3.169   (  13.006    9.711    8.001)   18.097
   3.857   (   3.955    1.208   -1.266)    4.325
   5.006   (   0.237   -3.708    3.796)    5.312
   5.274   (   0.774    0.107    1.985)    2.133
   5.400   (  -2.653    2.912    0.036)    3.939
   5.522   (   1.061    1.478   -0.919)    2.039
   5.687   (  -3.250   -1.456   -1.949)    4.060
   5.807   (  -0.396   -1.606   -2.451)    2.957
======================= Grid point 213 (38/84) =======================
q-point: ( 0.38  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.267   (  -3.416    6.933   -2.212)    8.039
   1.428   (   1.376    5.891   -4.719)    7.672
   1.461   (  -5.374    5.293    0.725)    7.578
   1.534   (  -3.518   -0.146   -7.684)    8.452
   3.426   (   6.033    4.591    7.514)   10.674
   3.890   (   1.616   -3.159    1.256)    3.764
   5.011   (   2.666   -1.742    2.843)    4.269
   5.242   (  -1.096   -3.648    1.311)    4.028
   5.394   (  -1.231    2.665   -0.184)    2.941
   5.558   (   1.118    2.321   -0.869)    2.719
   5.616   (  -3.568   -0.023   -1.485)    3.864
   5.795   (  -0.742    0.897   -2.202)    2.491
======================= Grid point 224 (39/84) =======================
q-point: ( 0.23  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.237   (   3.186    5.519   -7.096)    9.537
   1.392   (   2.562    4.437   -0.250)    5.129
   1.592   (   0.647    1.120    0.457)    1.372
   1.643   (  -0.335   -0.580   -5.382)    5.423
   3.095   (   8.387   14.527    8.261)   18.698
   3.825   (   1.798    3.114   -1.530)    3.908
   4.985   (  -1.742   -3.017    3.871)    5.208
   5.267   (   0.405    0.701    1.985)    2.144
   5.462   (  -1.026   -1.778   -1.941)    2.825
   5.500   (   1.716    2.972    0.707)    3.503
   5.713   (  -1.226   -2.124   -2.266)    3.339
   5.798   (  -0.065   -0.113   -1.836)    1.841
======================= Grid point 225 (40/84) =======================
q-point: ( 0.31  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.336   (   1.512    4.754   -8.194)    9.593
   1.458   (   1.012    3.492   -2.595)    4.466
   1.566   (  -3.184    0.106    0.889)    3.307
   1.607   (  -3.006   -0.192   -4.688)    5.572
   3.373   (   7.657   10.007    7.624)   14.727
   3.845   (  -0.043   -2.499    1.479)    2.904
   4.966   (   1.323    0.048    2.790)    3.088
   5.230   (  -2.249   -4.484    1.186)    5.155
   5.438   (  -0.609    1.911   -1.750)    2.662
   5.553   (   1.493    1.616    0.082)    2.202
   5.658   (  -2.252   -1.955   -1.783)    3.474
   5.818   (   0.281    2.056   -1.361)    2.482
======================= Grid point 226 (41/84) =======================
q-point: ( 0.38  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.376   (  -2.386    4.132   -7.050)    8.513
   1.497   (  -0.121    0.210   -5.216)    5.222
   1.544   (  -1.784    3.091   -0.451)    3.597
   1.549   (  -0.562    0.973   -2.447)    2.692
   3.516   (  -2.348    4.068    7.076)    8.493
   3.810   (   2.354   -4.077    3.203)    5.694
   4.997   (  -0.153    0.264    1.989)    2.012
   5.163   (   1.993   -3.451    0.904)    4.086
   5.450   (  -1.597    2.767   -1.425)    3.498
   5.592   (   0.771   -1.336   -2.513)    2.948
   5.606   (  -0.674    1.168    0.755)    1.546
   5.833   (  -1.174    2.033   -1.276)    2.672
======================= Grid point 239 (42/84) =======================
q-point: ( 0.31  0.31  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.429   (   2.272    3.936   -9.334)   10.381
   1.494   (   0.119    0.207   -7.205)    7.209
   1.543   (  -0.864   -1.496    1.629)    2.375
   1.581   (  -0.913   -1.581   -0.493)    1.891
   3.570   (   4.828    8.363    6.570)   11.680
   3.756   (  -3.248   -5.626    4.260)    7.769
   5.006   (   1.955    3.385    1.574)    4.214
   5.121   (  -2.973   -5.150    0.878)    6.011
   5.491   (   1.595    2.763   -2.246)    3.902
   5.561   (  -0.772   -1.338   -1.949)    2.487
   5.624   (  -0.387   -0.671    0.811)    1.121
   5.857   (   0.655    1.134   -1.050)    1.678
======================= Grid point 362 (43/84) =======================
q-point: ( 0.00  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.564   (   0.000   -0.000   11.507)   11.507
   0.564   (   0.000    0.000   11.507)   11.507
   0.942   (   0.000   -0.000   -1.532)    1.532
   0.942   (   0.000    0.000   -1.532)    1.532
   1.573   (   0.000    0.000   34.699)   34.699
   3.198   (  -0.000    0.000  -18.988)   18.988
   5.154   (   0.000    0.000    0.593)    0.593
   5.154   (   0.000    0.000    0.593)    0.593
   5.206   (  -0.000    0.000   -0.819)    0.819
   5.206   (   0.000    0.000   -0.819)    0.819
   5.717   (   0.000    0.000    7.848)    7.848
   5.815   (   0.000   -0.000    6.417)    6.417
======================= Grid point 363 (44/84) =======================
q-point: ( 0.08  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.629   (   5.074    2.929   10.306)   11.855
   0.820   (  13.948    8.053    8.606)   18.261
   0.970   (   2.354    1.359   -1.004)    2.897
   1.124   (  12.628    7.291   -1.831)   14.696
   1.619   (   6.302    3.638   32.111)   32.925
   3.185   (  -0.643   -0.371  -18.273)   18.288
   5.164   (   0.816    0.471    0.663)    1.152
   5.216   (   3.715    2.145    0.980)    4.401
   5.224   (   1.423    0.822   -0.970)    1.908
   5.250   (   3.221    1.860   -0.369)    3.737
   5.700   (  -0.176   -0.101    6.532)    6.535
   5.797   (  -1.330   -0.768    5.747)    5.948
======================= Grid point 364 (45/84) =======================
q-point: ( 0.15  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.779   (   6.902    3.985    8.240)   11.464
   1.021   (   3.423    1.976   12.185)   12.810
   1.052   (   4.339    2.505   -0.357)    5.023
   1.413   (   9.599    5.542   -3.349)   11.578
   1.949   (  19.828   11.448   20.583)   30.788
   3.203   (   3.000    1.732  -15.722)   16.100
   5.192   (   1.477    0.853    0.823)    1.893
   5.267   (   2.068    1.194   -1.206)    2.676
   5.278   (   0.991    0.572    3.716)    3.888
   5.314   (   1.782    1.029    1.388)    2.482
   5.743   (   2.164    1.249    1.670)    3.006
   5.775   (   1.019    0.588    2.431)    2.701
======================= Grid point 365 (46/84) =======================
q-point: ( 0.23  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.931   (   5.406    3.121    6.673)    9.138
   1.064   (   1.163    0.672    8.770)    8.872
   1.162   (   4.420    2.552   -0.469)    5.126
   1.535   (   0.745    0.430   -6.089)    6.149
   2.460   (  21.706   12.532   14.362)   28.887
   3.347   (   8.812    5.088  -11.515)   15.367
   5.230   (   1.631    0.942    1.108)    2.185
   5.280   (  -0.348   -0.201    4.492)    4.509
   5.320   (   2.368    1.367   -1.368)    3.058
   5.329   (  -0.425   -0.245    2.923)    2.964
   5.748   (  -0.236   -0.136    1.048)    1.083
   5.848   (   2.279    1.316   -3.667)    4.513
======================= Grid point 366 (47/84) =======================
q-point: ( 0.31  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.027   (   2.622    1.514    5.755)    6.503
   1.097   (   1.697    0.980    1.596)    2.528
   1.241   (   1.921    1.109   -0.782)    2.353
   1.487   (  -3.732   -2.155   -8.269)    9.324
   2.951   (  19.152   11.057   11.822)   25.076
   3.590   (  10.721    6.190   -7.471)   14.459
   5.259   (  -1.611   -0.930    4.139)    4.538
   5.269   (   1.668    0.963    1.456)    2.414
   5.323   (   0.553    0.319    4.003)    4.053
   5.377   (   2.330    1.345   -1.506)    3.083
   5.726   (  -1.958   -1.130   -1.387)    2.652
   5.846   (  -2.445   -1.412   -4.896)    5.652
======================= Grid point 367 (48/84) =======================
q-point: ( 0.38  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.063   (   0.659    0.380    5.239)    5.294
   1.145   (   2.233    1.289   -4.815)    5.462
   1.250   (  -0.999   -0.577   -0.470)    1.245
   1.392   (  -3.830   -2.212   -9.827)   10.776
   3.358   (  14.588    8.423   10.451)   19.824
   3.824   (   8.322    4.805   -4.045)   10.427
   5.216   (  -1.535   -0.886    5.409)    5.692
   5.308   (   1.504    0.868    1.638)    2.388
   5.351   (   1.383    0.798    3.239)    3.611
   5.427   (   1.745    1.007   -1.623)    2.587
   5.640   (  -5.436   -3.138   -3.049)    6.978
   5.762   (  -3.557   -2.054   -4.969)    6.447
======================= Grid point 368 (49/84) =======================
q-point: ( 0.46  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.071   (   0.151    0.087    4.850)    4.853
   1.195   (   1.671    0.965   -8.522)    8.738
   1.216   (  -1.170   -0.676    0.316)    1.388
   1.317   (  -2.119   -1.224  -10.791)   11.065
   3.617   (   6.246    3.606    9.555)   11.971
   3.960   (   2.806    1.620   -1.236)    3.467
   5.205   (   0.521    0.301    5.167)    5.202
   5.335   (   0.647    0.374    1.633)    1.796
   5.366   (  -0.457   -0.264    3.669)    3.707
   5.456   (   0.649    0.375   -1.684)    1.843
   5.511   (  -3.537   -2.042   -3.553)    5.413
   5.702   (  -1.328   -0.767   -5.326)    5.543
======================= Grid point 376 (50/84) =======================
q-point: ( 0.08  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.770   (   5.255    9.101    8.356)   13.426
   0.973   (   3.163    5.479   10.239)   12.036
   1.052   (   3.019    5.229    0.695)    6.078
   1.334   (   6.227   10.785   -2.973)   12.803
   1.822   (   9.563   16.563   23.993)   30.683
   3.188   (   0.808    1.399  -16.636)   16.715
   5.195   (   1.238    2.144    0.582)    2.544
   5.244   (   1.278    2.214   -0.041)    2.556
   5.281   (   0.902    1.562    2.785)    3.318
   5.293   (   1.316    2.280    0.021)    2.633
   5.724   (   1.673    2.898    2.195)    4.002
   5.773   (  -0.744   -1.288    4.436)    4.679
======================= Grid point 377 (51/84) =======================
q-point: ( 0.15  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.950   (   2.275   10.692    6.662)   12.801
   1.042   (   2.205    0.492   10.258)   10.504
   1.178   (   2.563    9.011   -0.317)    9.374
   1.500   (   3.743    2.424   -5.187)    6.841
   2.250   (  17.054   16.856   16.438)   29.072
   3.273   (   5.254    5.275  -13.200)   15.155
   5.237   (   0.364    2.602    0.749)    2.732
   5.279   (  -0.052   -0.004    3.952)    3.952
   5.309   (   1.109    2.236    0.442)    2.535
   5.328   (   1.295    0.841    1.110)    1.902
   5.750   (  -0.525   -0.300    1.926)    2.019
   5.804   (   3.268    1.389   -2.277)    4.218
======================= Grid point 378 (52/84) =======================
q-point: ( 0.23  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.036   (   0.919    2.420    4.244)    4.971
   1.114   (   0.374    7.505    4.892)    8.967
   1.298   (  -0.483    9.529   -0.532)    9.556
   1.519   (  -1.501   -1.273   -7.240)    7.503
   2.732   (  18.515   13.439   12.955)   26.292
   3.473   (   9.282    6.920   -9.165)   14.766
   5.241   (  -0.593   -2.240    2.707)    3.563
   5.280   (   0.303    1.149    2.938)    3.169
   5.333   (  -0.177    2.444    2.769)    3.698
   5.383   (   1.254    3.799   -0.993)    4.122
   5.735   (  -0.722   -1.258   -0.503)    1.536
   5.839   (   0.697   -2.235   -4.294)    4.891
======================= Grid point 379 (53/84) =======================
q-point: ( 0.31  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.076   (   1.573    1.181   -1.094)    2.251
   1.175   (  -1.988    7.468    2.275)    8.056
   1.345   (  -4.121    7.552   -1.115)    8.675
   1.451   (  -4.394   -0.587   -8.016)    9.160
   3.170   (  16.492   10.007   11.155)   22.283
   3.714   (   9.477    5.562   -5.559)   12.315
   5.193   (  -0.364   -4.216    4.685)    6.313
   5.306   (   1.053    1.199    1.879)    2.465
   5.357   (  -0.391    3.006    2.810)    4.134
   5.438   (   0.414    3.755   -1.342)    4.010
   5.684   (  -2.951   -2.733   -2.376)    4.672
   5.788   (  -2.219   -3.442   -4.659)    6.204
======================= Grid point 380 (54/84) =======================
q-point: ( 0.38  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.110   (   0.150    2.006   -2.551)    3.249
   1.213   (  -1.263    6.336   -1.911)    6.738
   1.321   (  -5.609    5.540   -2.103)    8.159
   1.366   (  -4.814    0.814   -7.428)    8.889
   3.507   (  11.160    5.960    9.995)   16.124
   3.897   (   5.937    2.052   -2.222)    6.663
   5.159   (   1.488   -3.656    5.184)    6.516
   5.328   (   0.122   -0.297    1.457)    1.491
   5.373   (  -0.200    0.734    2.457)    2.572
   5.474   (  -0.139    2.989   -1.608)    3.397
   5.581   (  -4.995   -1.321   -2.557)    5.765
   5.723   (  -1.935   -1.133   -4.813)    5.310
======================= Grid point 381 (55/84) =======================
q-point: ( 0.46  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.121   (  -1.396    2.418   -1.878)    3.364
   1.256   (  -3.513    6.084   -3.398)    7.804
   1.280   (  -1.799    3.116   -6.419)    7.359
   1.320   (  -2.492    4.317   -4.341)    6.610
   3.649   (   0.296   -0.512    9.405)    9.424
   3.958   (   0.952   -1.650   -0.367)    1.940
   5.164   (   2.026   -3.510    4.580)    6.115
   5.307   (   1.281   -2.219    1.725)    3.089
   5.374   (  -0.328    0.568    1.729)    1.849
   5.501   (  -1.644    2.847   -1.858)    3.776
   5.517   (  -1.399    2.423   -1.706)    3.277
   5.700   (  -0.445    0.771   -4.913)    4.993
======================= Grid point 391 (56/84) =======================
q-point: ( 0.15  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.036   (   0.205    0.356    4.911)    4.928
   1.169   (   4.737    8.205    4.079)   10.315
   1.372   (   4.920    8.522   -0.241)    9.843
   1.520   (   0.001    0.002   -6.896)    6.896
   2.622   (  10.876   18.838   13.837)   25.780
   3.422   (   5.178    8.969   -9.954)   14.365
   5.235   (  -1.438   -2.491    2.662)    3.919
   5.282   (   0.159    0.275    3.910)    3.923
   5.368   (   2.044    3.539    1.141)    4.243
   5.380   (   2.015    3.490   -0.823)    4.114
   5.736   (  -0.677   -1.173   -0.225)    1.373
   5.810   (  -0.523   -0.905   -3.701)    3.846
======================= Grid point 392 (57/84) =======================
q-point: ( 0.23  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.054   (   1.280    0.809   -1.484)    2.120
   1.277   (  -0.186    5.269   -0.289)    5.281
   1.490   (  -1.866    7.149   -1.747)    7.593
   1.520   (  -0.866    2.676   -4.095)    4.968
   3.013   (  13.436   13.857   12.013)   22.734
   3.625   (   7.292    7.581   -6.663)   12.452
   5.168   (  -1.706   -4.198    4.613)    6.467
   5.297   (   0.860    0.780    2.441)    2.703
   5.414   (  -1.300    3.727    0.652)    4.001
   5.456   (   2.092    2.808   -0.489)    3.536
   5.699   (  -1.684   -2.005   -2.070)    3.338
   5.774   (  -0.610   -3.175   -3.927)    5.087
======================= Grid point 393 (58/84) =======================
q-point: ( 0.31  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.096   (   1.962    1.218   -6.446)    6.847
   1.297   (  -1.321    2.306   -5.849)    6.424
   1.475   (  -5.500    5.477   -2.220)    8.073
   1.519   (  -5.087    5.874   -0.059)    7.771
   3.368   (  12.214    9.369   10.618)   18.700
   3.813   (   6.268    3.909   -3.159)    8.034
   5.104   (  -0.150   -4.000    5.390)    6.714
   5.319   (   0.486   -0.157    1.606)    1.685
   5.395   (  -2.017   -0.618   -0.216)    2.120
   5.505   (   0.326    2.665   -0.670)    2.768
   5.639   (  -2.776   -1.110   -2.231)    3.730
   5.742   (  -0.802   -0.151   -3.657)    3.746
======================= Grid point 394 (59/84) =======================
q-point: ( 0.38  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.138   (   0.839    1.263   -8.572)    8.705
   1.291   (  -0.419   -0.017  -11.035)   11.043
   1.435   (  -5.891    7.735    0.779)    9.754
   1.480   (  -5.425    6.280    0.779)    8.335
   3.609   (   5.610    4.175    9.494)   11.792
   3.909   (   2.620   -0.753    0.143)    2.730
   5.090   (   2.096   -2.672    4.585)    5.706
   5.273   (  -1.248   -4.309    1.317)    4.676
   5.383   (   0.059    0.790   -0.740)    1.084
   5.531   (   0.274    1.557   -2.291)    2.783
   5.588   (  -2.887    1.473   -0.364)    3.261
   5.737   (  -1.423    2.462   -3.314)    4.367
======================= Grid point 405 (60/84) =======================
q-point: ( 0.23  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.083   (   1.160    2.008   -6.524)    6.924
   1.318   (  -0.082   -0.142   -7.937)    7.939
   1.566   (   1.161    2.011   -0.724)    2.433
   1.602   (   1.060    1.837    0.331)    2.146
   3.301   (   7.930   13.735   10.983)   19.292
   3.774   (   3.684    6.380   -3.604)    8.201
   5.085   (  -1.938   -3.357    5.532)    6.755
   5.313   (   0.198    0.342    1.731)    1.776
   5.414   (  -1.748   -3.027   -2.036)    4.046
   5.514   (   1.399    2.424    0.350)    2.821
   5.658   (  -1.062   -1.840   -2.477)    3.263
   5.748   (   0.510    0.884   -2.843)    3.021
======================= Grid point 406 (61/84) =======================
q-point: ( 0.31  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.134   (   1.613    2.473   -9.931)   10.361
   1.299   (  -1.003   -1.077  -11.560)   11.653
   1.576   (  -2.618    3.273    1.266)    4.379
   1.595   (  -2.551    1.180    1.231)    3.068
   3.559   (   6.982    8.983    9.639)   14.912
   3.867   (   1.938    1.083    0.248)    2.234
   5.045   (   0.472   -1.363    4.787)    4.999
   5.259   (  -2.724   -5.347    1.298)    6.140
   5.390   (   0.437    1.261   -2.500)    2.834
   5.533   (  -0.217   -0.649   -2.780)    2.863
   5.636   (  -0.863    0.938    0.495)    1.367
   5.780   (   0.130    2.765   -2.329)    3.617
======================= Grid point 407 (62/84) =======================
q-point: ( 0.38  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.170   (  -1.006    1.743  -11.106)   11.287
   1.277   (   0.509   -0.882  -12.766)   12.807
   1.567   (  -2.489    4.310    1.853)    5.311
   1.569   (  -1.348    2.335    1.797)    3.240
   3.686   (  -1.840    3.187    8.581)    9.337
   3.877   (   1.134   -1.965    2.698)    3.525
   5.056   (   0.337   -0.583    3.694)    3.755
   5.188   (   2.000   -3.464    1.345)    4.220
   5.407   (  -0.864    1.497   -2.480)    3.023
   5.520   (   0.964   -1.669   -4.287)    4.700
   5.635   (  -1.260    2.183    2.073)    3.263
   5.797   (  -1.486    2.574   -2.237)    3.720
======================= Grid point 420 (63/84) =======================
q-point: ( 0.31  0.31  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.194   (   1.730    2.997  -12.163)   12.646
   1.264   (  -1.258   -2.180  -13.001)   13.243
   1.589   (  -0.412   -0.713    2.471)    2.605
   1.611   (   0.277    0.479    1.933)    2.011
   3.726   (   3.796    6.574    7.783)   10.872
   3.849   (  -1.548   -2.682    4.125)    5.158
   5.055   (   1.211    2.098    3.162)    3.983
   5.147   (  -2.743   -4.750    1.567)    5.705
   5.430   (   1.300    2.251   -3.390)    4.272
   5.495   (  -1.336   -2.314   -4.150)    4.936
   5.659   (   0.449    0.777    2.449)    2.608
   5.825   (   0.716    1.239   -2.026)    2.481
======================= Grid point 543 (64/84) =======================
q-point: ( 0.00  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.763   (   0.000    0.000    7.893)    7.893
   0.763   (   0.000   -0.000    7.893)    7.893
   0.886   (   0.000    0.000   -4.256)    4.256
   0.886   (   0.000    0.000   -4.256)    4.256
   2.224   (   0.000    0.000   29.113)   29.113
   2.762   (   0.000    0.000  -23.840)   23.840
   5.169   (   0.000   -0.000    0.841)    0.841
   5.169   (   0.000    0.000    0.841)    0.841
   5.188   (   0.000   -0.000   -0.942)    0.942
   5.188   (  -0.000    0.000   -0.942)    0.942
   5.854   (   0.000    0.000    4.918)    4.918
   5.896   (   0.000   -0.000    1.044)    1.044
======================= Grid point 544 (65/84) =======================
q-point: ( 0.08  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.811   (   3.851    2.223    7.328)    8.571
   0.924   (   3.078    1.777   -3.777)    5.186
   0.969   (  13.242    7.645    6.090)   16.459
   1.069   (  12.447    7.186   -3.739)   14.851
   2.239   (   2.360    1.362   28.230)   28.361
   2.762   (   0.676    0.390  -23.184)   23.197
   5.181   (   0.965    0.557    0.986)    1.488
   5.203   (   1.190    0.687   -1.122)    1.774
   5.225   (   3.825    2.208    0.354)    4.431
   5.237   (   3.811    2.200   -0.728)    4.460
   5.826   (  -1.927   -1.113    4.763)    5.258
   5.869   (  -2.044   -1.180    0.807)    2.494
======================= Grid point 545 (66/84) =======================
q-point: ( 0.15  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.929   (   5.576    3.219    6.103)    8.871
   1.021   (   4.752    2.744   -2.914)    6.213
   1.215   (   5.528    3.192    7.231)    9.646
   1.333   (   7.608    4.392   -4.633)    9.932
   2.402   (  11.946    6.897   22.627)   26.500
   2.839   (   6.610    3.816  -19.902)   21.315
   5.213   (   1.674    0.966    1.218)    2.284
   5.240   (   1.880    1.085   -1.389)    2.577
   5.310   (   2.697    1.557    0.885)    3.237
   5.328   (   2.835    1.637   -0.548)    3.319
   5.791   (  -0.632   -0.365    2.679)    2.777
   5.817   (  -1.824   -1.053    0.544)    2.175
======================= Grid point 546 (67/84) =======================
q-point: ( 0.23  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.055   (   4.577    2.643    5.124)    7.361
   1.133   (   4.263    2.461   -2.542)    5.539
   1.243   (  -1.920   -1.109    8.409)    8.696
   1.400   (  -1.323   -0.764   -7.083)    7.246
   2.772   (  17.728   10.235   15.725)   25.813
   3.082   (  13.032    7.524  -14.343)   20.789
   5.257   (   1.939    1.119    1.510)    2.701
   5.289   (   2.200    1.270   -1.627)    3.017
   5.340   (   0.077    0.044    1.877)    1.879
   5.362   (   0.393    0.227    0.024)    0.454
   5.773   (  -1.389   -0.802    1.253)    2.035
   5.798   (  -0.301   -0.174   -1.426)    1.468
======================= Grid point 547 (68/84) =======================
q-point: ( 0.31  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.136   (   2.151    1.242    4.627)    5.252
   1.173   (  -3.304   -1.907    5.722)    6.877
   1.209   (   1.924    1.111   -2.532)    3.368
   1.317   (  -4.759   -2.748   -8.070)    9.763
   3.195   (  17.065    9.853   11.649)   22.891
   3.414   (  14.022    8.096   -9.733)   18.891
   5.303   (   1.935    1.117    1.862)    2.909
   5.332   (  -0.492   -0.284    2.866)    2.921
   5.342   (   2.172    1.254   -1.883)    3.137
   5.366   (   0.294    0.170   -0.025)    0.340
   5.719   (  -3.389   -1.957    0.863)    4.008
   5.761   (  -3.097   -1.788   -3.194)    4.795
======================= Grid point 548 (69/84) =======================
q-point: ( 0.38  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.105   (  -2.140   -1.236    1.484)    2.882
   1.161   (   0.208    0.120    4.125)    4.132
   1.204   (  -4.280   -2.471   -7.762)    9.201
   1.225   (  -0.375   -0.216   -2.105)    2.149
   3.562   (  13.219    7.632    9.104)   17.773
   3.720   (  11.026    6.366   -6.354)   14.229
   5.320   (  -0.574   -0.331    4.480)    4.528
   5.346   (   1.571    0.907    2.066)    2.750
   5.382   (   0.899    0.519   -0.960)    1.414
   5.389   (   1.658    0.957   -2.058)    2.811
   5.608   (  -5.642   -3.258    0.736)    6.557
   5.665   (  -4.257   -2.458   -4.204)    6.468
======================= Grid point 549 (70/84) =======================
q-point: ( 0.46  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.076   (  -0.369   -0.213   -2.140)    2.181
   1.126   (  -1.980   -1.143   -6.798)    7.172
   1.159   (  -0.170   -0.098    3.508)    3.514
   1.208   (  -0.635   -0.366   -1.316)    1.506
   3.795   (   5.519    3.186    7.401)    9.767
   3.910   (   4.340    2.506   -3.916)    6.360
   5.310   (  -0.082   -0.047    4.958)    4.959
   5.374   (   0.627    0.362    2.090)    2.212
   5.397   (   0.062    0.036   -2.837)    2.838
   5.417   (   0.629    0.363   -2.110)    2.232
   5.487   (  -3.251   -1.877    3.377)    5.049
   5.588   (  -1.817   -1.049   -5.641)    6.018
======================= Grid point 557 (71/84) =======================
q-point: ( 0.08  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.921   (   4.261    7.381    6.133)   10.500
   1.011   (   3.717    6.439   -2.648)    7.892
   1.158   (   4.343    7.522    6.250)   10.700
   1.262   (   5.319    9.212   -4.184)   11.431
   2.331   (   5.030    8.712   24.752)   26.718
   2.802   (   2.493    4.318  -21.170)   21.749
   5.211   (   1.192    2.065    0.904)    2.550
   5.230   (   1.210    2.097   -0.883)    2.577
   5.292   (   1.786    3.094    0.127)    3.574
   5.300   (   1.953    3.383   -0.396)    3.927
   5.795   (  -0.747   -1.293    3.510)    3.815
   5.829   (  -1.329   -2.303    0.618)    2.730
======================= Grid point 558 (72/84) =======================
q-point: ( 0.15  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.073   (   2.055    9.325    5.073)   10.813
   1.149   (   1.823    8.111   -2.308)    8.627
   1.237   (   1.519   -0.389    8.179)    8.328
   1.384   (   2.246    0.373   -6.315)    6.712
   2.610   (  12.725   12.909   18.197)   25.685
   2.968   (   8.560    8.733  -16.534)   20.565
   5.255   (   0.843    2.599    1.075)    2.936
   5.279   (   1.247    2.391   -1.174)    2.941
   5.339   (   0.874    1.333    0.998)    1.881
   5.350   (   1.349    0.615    0.240)    1.502
   5.779   (  -0.685   -1.010    0.933)    1.536
   5.789   (  -0.181   -1.544    0.110)    1.558
======================= Grid point 559 (73/84) =======================
q-point: ( 0.23  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.140   (  -1.362    0.120    5.257)    5.432
   1.195   (  -0.699    2.981    1.572)    3.442
   1.288   (  -1.301    8.339    0.767)    8.475
   1.374   (  -2.321    0.194   -6.514)    6.918
   3.004   (  15.888   11.848   13.321)   23.880
   3.260   (  12.559    9.458  -11.640)   19.562
   5.283   (   0.540   -0.128    1.672)    1.762
   5.307   (   0.524    0.337   -0.319)    0.700
   5.371   (   0.406    2.994    0.655)    3.092
   5.376   (   0.414    2.413    0.347)    2.472
   5.738   (  -1.684   -2.599    0.841)    3.209
   5.770   (  -0.740   -2.644   -2.367)    3.625
======================= Grid point 560 (74/84) =======================
q-point: ( 0.31  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.090   (  -0.989   -3.377    3.125)    4.706
   1.182   (  -1.832   -1.484   -3.944)    4.595
   1.310   (  -5.284    7.617    0.161)    9.272
   1.347   (  -5.037    7.273   -2.473)    9.186
   3.393   (  14.782    9.122   10.262)   20.175
   3.578   (  12.294    7.544   -7.845)   16.420
   5.276   (   0.766   -2.811    3.300)    4.402
   5.321   (   1.187   -0.662   -0.758)    1.556
   5.403   (  -0.249    3.078    1.248)    3.331
   5.418   (  -0.140    4.088   -0.522)    4.124
   5.660   (  -3.636   -3.310    0.395)    4.933
   5.702   (  -2.922   -3.303   -3.478)    5.616
======================= Grid point 561 (75/84) =======================
q-point: ( 0.38  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.051   (   0.151   -2.420   -0.280)    2.441
   1.125   (  -1.464   -2.629   -6.821)    7.455
   1.286   (  -6.184    8.356    1.952)   10.577
   1.323   (  -6.021    7.565   -1.300)    9.755
   3.696   (   9.976    5.263    8.090)   13.880
   3.827   (   8.170    3.934   -4.814)   10.266
   5.260   (   1.384   -3.723    4.481)    5.988
   5.329   (   1.083   -2.404   -1.814)    3.200
   5.417   (  -0.429    0.595    1.725)    1.875
   5.448   (  -0.630    2.984   -1.130)    3.253
   5.559   (  -4.546   -0.109    1.179)    4.698
   5.624   (  -2.895   -0.486   -4.542)    5.408
======================= Grid point 562 (76/84) =======================
q-point: ( 0.46  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.042   (   1.001   -1.734   -2.047)    2.863
   1.095   (   0.677   -1.173   -7.044)    7.173
   1.272   (  -5.117    8.863    2.187)   10.465
   1.305   (  -4.692    8.126   -0.824)    9.420
   3.821   (   0.364   -0.630    6.990)    7.027
   3.924   (   0.583   -1.010   -3.222)    3.426
   5.255   (   2.161   -3.742    4.022)    5.903
   5.315   (   1.919   -3.324   -1.450)    4.102
   5.412   (   0.281   -0.487    2.026)    2.102
   5.458   (  -0.907    1.571   -2.309)    2.937
   5.519   (  -2.454    4.251    2.283)    5.413
   5.597   (  -1.351    2.340   -4.860)    5.560
======================= Grid point 572 (77/84) =======================
q-point: ( 0.15  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.156   (  -2.097   -3.633    6.090)    7.395
   1.227   (   1.503    2.604    0.508)    3.049
   1.362   (   4.809    8.329   -0.591)    9.635
   1.390   (   1.452    2.515   -4.724)    5.545
   2.915   (   9.290   16.090   14.326)   23.461
   3.191   (   7.222   12.508  -12.619)   19.179
   5.279   (  -0.354   -0.614    1.765)    1.902
   5.309   (   0.041    0.070   -0.919)    0.923
   5.384   (   1.750    3.031    0.737)    3.577
   5.386   (   1.900    3.292    0.513)    3.835
   5.739   (  -1.561   -2.703    0.406)    3.148
   5.755   (  -1.211   -2.098   -1.525)    2.862
======================= Grid point 573 (78/84) =======================
q-point: ( 0.23  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.080   (  -1.763   -3.467    3.984)    5.568
   1.197   (  -2.173   -2.943   -6.475)    7.437
   1.466   (  -1.152    7.783    0.762)    7.905
   1.489   (  -0.500    7.828   -1.113)    7.922
   3.256   (  11.836   12.468   11.368)   20.610
   3.465   (   9.660   10.231   -9.080)   16.746
   5.250   (  -0.720   -2.390    3.398)    4.216
   5.304   (   0.070   -0.584   -1.646)    1.748
   5.436   (   0.076    2.486    1.430)    2.868
   5.454   (   0.958    3.667   -0.232)    3.797
   5.674   (  -2.015   -3.044    0.013)    3.651
   5.703   (  -1.251   -2.932   -2.682)    4.166
======================= Grid point 574 (79/84) =======================
q-point: ( 0.31  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.030   (  -0.465   -2.308    0.769)    2.477
   1.128   (  -1.859   -2.840   -8.099)    8.782
   1.487   (  -5.158    7.929    1.602)    9.593
   1.513   (  -4.979    7.172   -0.792)    8.766
   3.572   (  10.813    8.338    8.945)   16.324
   3.723   (   8.753    6.491   -5.814)   12.351
   5.211   (  -0.124   -3.315    4.892)    5.911
   5.297   (   0.109   -1.657   -3.306)    3.700
   5.425   (  -1.249   -1.704    2.804)    3.510
   5.481   (  -0.721    0.848   -2.159)    2.429
   5.622   (  -2.003    0.648    0.966)    2.316
   5.669   (  -1.358    1.204   -3.255)    3.727
======================= Grid point 575 (80/84) =======================
q-point: ( 0.38  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.013   (   0.511   -1.234   -2.055)    2.451
   1.079   (  -0.612   -1.877   -8.020)    8.260
   1.468   (  -5.585    8.231    1.575)   10.071
   1.487   (  -5.434    7.298   -0.221)    9.102
   3.782   (   4.882    3.216    6.949)    9.081
   3.883   (   3.815    1.574   -2.959)    5.078
   5.186   (   1.233   -3.172    4.556)    5.687
   5.263   (  -0.165   -3.580   -2.621)    4.441
   5.395   (   0.214   -1.432    2.427)    2.826
   5.464   (  -0.186   -1.205   -3.815)    4.005
   5.611   (  -2.357    4.012    2.248)    5.167
   5.669   (  -2.076    3.906   -3.213)    5.467
======================= Grid point 586 (81/84) =======================
q-point: ( 0.23  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.022   (  -1.155   -2.001    1.224)    2.615
   1.130   (  -1.769   -3.064   -8.621)    9.319
   1.582   (   1.800    3.118    1.293)    3.826
   1.601   (   1.528    2.646   -0.469)    3.091
   3.513   (   7.010   12.141    9.412)   16.886
   3.674   (   5.565    9.638   -6.334)   12.805
   5.196   (  -1.540   -2.668    5.123)    5.978
   5.288   (  -0.637   -1.103   -3.721)    3.933
   5.428   (  -1.598   -2.768    2.885)    4.305
   5.491   (  -0.408   -0.706   -2.698)    2.818
   5.643   (   0.316    0.547    1.211)    1.365
   5.689   (   0.835    1.446   -2.840)    3.294
======================= Grid point 587 (82/84) =======================
q-point: ( 0.31  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.999   (   0.048   -0.594   -1.897)    1.989
   1.075   (  -1.253   -2.223   -8.736)    9.101
   1.605   (  -2.262    3.369    1.091)    4.202
   1.613   (  -2.330    2.445    0.339)    3.395
   3.735   (   5.815    7.239    7.167)   11.730
   3.840   (   4.082    4.615   -3.116)    6.905
   5.150   (  -0.368   -2.325    5.199)    5.707
   5.243   (  -1.595   -3.588   -3.405)    5.198
   5.381   (  -0.319   -1.345    2.267)    2.655
   5.454   (  -0.946   -2.308   -4.273)    4.948
   5.671   (  -0.104    2.925    2.518)    3.861
   5.727   (   0.035    3.186   -2.785)    4.232
======================= Grid point 588 (83/84) =======================
q-point: ( 0.38  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.001   (   0.022   -0.038   -3.811)    3.811
   1.048   (   0.626   -1.084   -8.056)    8.152
   1.598   (  -2.361    4.089    0.892)    4.805
   1.599   (  -1.930    3.343    0.828)    3.948
   3.834   (  -1.099    1.903    5.505)    5.928
   3.899   (  -0.117    0.203   -0.876)    0.906
   5.138   (   0.771   -1.335    4.021)    4.306
   5.198   (   1.441   -2.496   -1.139)    3.100
   5.377   (   0.187   -0.323    0.414)    0.557
   5.430   (   0.944   -1.634   -4.140)    4.550
   5.686   (  -1.603    2.776    2.815)    4.265
   5.744   (  -1.667    2.886   -2.799)    4.353
======================= Grid point 601 (84/84) =======================
q-point: ( 0.31  0.31  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 2.72e-04 2.72e-04 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.003   (   0.483    0.837   -4.826)    4.922
   1.036   (  -0.901   -1.560   -7.721)    7.928
   1.631   (   0.007    0.012    1.362)    1.362
   1.639   (   0.271    0.469    0.598)    0.807
   3.856   (   2.312    4.004    4.525)    6.469
   3.899   (   0.471    0.815    0.354)    1.006
   5.127   (   0.117    0.203    3.489)    3.497
   5.170   (  -1.683   -2.915   -0.128)    3.368
   5.377   (   0.362    0.627   -0.970)    1.210
   5.410   (  -1.007   -1.745   -3.647)    4.166
   5.716   (   0.674    1.167    2.989)    3.279
   5.776   (   0.740    1.282   -2.762)    3.134
=================== 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    183.159    183.159    233.725      0.000     -0.000     -0.000 3/14196
   20.0    129.896    129.896    159.751      0.000     -0.000     -0.000 3/14196
   30.0    104.208    104.208    121.460      0.000     -0.000     -0.000 3/14196
   40.0     84.509     84.509     95.636      0.000     -0.000     -0.000 3/14196
   50.0     69.954     69.954     77.768      0.000     -0.000     -0.000 3/14196
   60.0     59.363     59.363     65.239      0.000     -0.000     -0.000 3/14196
   70.0     51.521     51.521     56.191      0.000     -0.000     -0.000 3/14196
   80.0     45.544     45.544     49.424      0.000     -0.000     -0.000 3/14196
   90.0     40.858     40.858     44.190      0.000     -0.000     -0.000 3/14196
  100.0     37.087     37.087     40.023      0.000     -0.000     -0.000 3/14196
  110.0     33.987     33.987     36.622      0.000     -0.000     -0.000 3/14196
  120.0     31.390     31.390     33.790      0.000     -0.000     -0.000 3/14196
  130.0     29.182     29.182     31.391      0.000     -0.000     -0.000 3/14196
  140.0     27.278     27.278     29.331      0.000     -0.000     -0.000 3/14196
  150.0     25.619     25.619     27.540      0.000     -0.000     -0.000 3/14196
  160.0     24.159     24.159     25.966      0.000     -0.000     -0.000 3/14196
  170.0     22.863     22.863     24.571      0.000     -0.000     -0.000 3/14196
  180.0     21.704     21.704     23.326      0.000     -0.000     -0.000 3/14196
  190.0     20.662     20.662     22.206      0.000     -0.000     -0.000 3/14196
  200.0     19.718     19.718     21.193      0.000     -0.000     -0.000 3/14196
  210.0     18.860     18.860     20.272      0.000     -0.000     -0.000 3/14196
  220.0     18.075     18.075     19.431      0.000     -0.000     -0.000 3/14196
  230.0     17.355     17.355     18.659      0.000     -0.000     -0.000 3/14196
  240.0     16.692     16.692     17.948      0.000     -0.000     -0.000 3/14196
  250.0     16.079     16.079     17.291      0.000     -0.000     -0.000 3/14196
  260.0     15.511     15.511     16.682      0.000     -0.000     -0.000 3/14196
  270.0     14.982     14.982     16.115      0.000     -0.000     -0.000 3/14196
  280.0     14.489     14.489     15.587      0.000     -0.000     -0.000 3/14196
  290.0     14.028     14.028     15.093      0.000     -0.000     -0.000 3/14196
  300.0     13.596     13.596     14.630      0.000     -0.000     -0.000 3/14196
  310.0     13.190     13.190     14.195      0.000     -0.000     -0.000 3/14196
  320.0     12.808     12.808     13.786      0.000     -0.000     -0.000 3/14196
  330.0     12.449     12.449     13.401      0.000     -0.000     -0.000 3/14196
  340.0     12.109     12.109     13.036      0.000     -0.000     -0.000 3/14196
  350.0     11.787     11.787     12.692      0.000     -0.000     -0.000 3/14196
  360.0     11.483     11.483     12.365      0.000     -0.000     -0.000 3/14196
  370.0     11.194     11.194     12.055      0.000     -0.000     -0.000 3/14196
  380.0     10.919     10.919     11.761      0.000     -0.000     -0.000 3/14196
  390.0     10.658     10.658     11.481      0.000     -0.000     -0.000 3/14196
  400.0     10.409     10.409     11.214      0.000     -0.000     -0.000 3/14196
  410.0     10.172     10.172     10.959      0.000     -0.000     -0.000 3/14196
  420.0      9.945      9.945     10.716      0.000     -0.000     -0.000 3/14196
  430.0      9.729      9.729     10.484      0.000     -0.000     -0.000 3/14196
  440.0      9.522      9.522     10.262      0.000     -0.000     -0.000 3/14196
  450.0      9.323      9.323     10.049      0.000     -0.000     -0.000 3/14196
  460.0      9.133      9.133      9.844      0.000     -0.000     -0.000 3/14196
  470.0      8.950      8.950      9.648      0.000     -0.000     -0.000 3/14196
  480.0      8.775      8.775      9.460      0.000     -0.000     -0.000 3/14196
  490.0      8.606      8.606      9.279      0.000     -0.000     -0.000 3/14196
  500.0      8.444      8.444      9.105      0.000     -0.000     -0.000 3/14196
  510.0      8.288      8.288      8.938      0.000     -0.000     -0.000 3/14196
  520.0      8.138      8.138      8.776      0.000     -0.000     -0.000 3/14196
  530.0      7.993      7.993      8.620      0.000     -0.000     -0.000 3/14196
  540.0      7.853      7.853      8.470      0.000     -0.000     -0.000 3/14196
  550.0      7.718      7.718      8.325      0.000     -0.000     -0.000 3/14196
  560.0      7.588      7.588      8.185      0.000     -0.000     -0.000 3/14196
  570.0      7.462      7.462      8.050      0.000     -0.000     -0.000 3/14196
  580.0      7.340      7.340      7.919      0.000     -0.000     -0.000 3/14196
  590.0      7.222      7.222      7.793      0.000     -0.000     -0.000 3/14196
  600.0      7.108      7.108      7.670      0.000     -0.000     -0.000 3/14196
  610.0      6.998      6.998      7.551      0.000     -0.000     -0.000 3/14196
  620.0      6.891      6.891      7.436      0.000     -0.000     -0.000 3/14196
  630.0      6.787      6.787      7.324      0.000     -0.000     -0.000 3/14196
  640.0      6.686      6.686      7.216      0.000     -0.000     -0.000 3/14196
  650.0      6.588      6.588      7.111      0.000     -0.000     -0.000 3/14196
  660.0      6.493      6.493      7.009      0.000     -0.000     -0.000 3/14196
  670.0      6.401      6.401      6.910      0.000     -0.000     -0.000 3/14196
  680.0      6.311      6.311      6.813      0.000     -0.000     -0.000 3/14196
  690.0      6.224      6.224      6.720      0.000     -0.000     -0.000 3/14196
  700.0      6.140      6.140      6.629      0.000     -0.000     -0.000 3/14196
  710.0      6.057      6.057      6.540      0.000     -0.000     -0.000 3/14196
  720.0      5.977      5.977      6.454      0.000     -0.000     -0.000 3/14196
  730.0      5.899      5.899      6.370      0.000     -0.000     -0.000 3/14196
  740.0      5.823      5.823      6.288      0.000     -0.000     -0.000 3/14196
  750.0      5.749      5.749      6.208      0.000     -0.000     -0.000 3/14196
  760.0      5.676      5.676      6.130      0.000     -0.000     -0.000 3/14196
  770.0      5.606      5.606      6.054      0.000     -0.000     -0.000 3/14196
  780.0      5.537      5.537      5.980      0.000     -0.000     -0.000 3/14196
  790.0      5.470      5.470      5.908      0.000     -0.000     -0.000 3/14196
  800.0      5.405      5.405      5.838      0.000     -0.000     -0.000 3/14196
  810.0      5.341      5.341      5.769      0.000     -0.000     -0.000 3/14196
  820.0      5.279      5.279      5.702      0.000     -0.000     -0.000 3/14196
  830.0      5.218      5.218      5.636      0.000     -0.000     -0.000 3/14196
  840.0      5.158      5.158      5.572      0.000     -0.000     -0.000 3/14196
  850.0      5.100      5.100      5.510      0.000     -0.000     -0.000 3/14196
  860.0      5.043      5.043      5.449      0.000     -0.000     -0.000 3/14196
  870.0      4.988      4.988      5.389      0.000     -0.000     -0.000 3/14196
  880.0      4.933      4.933      5.330      0.000     -0.000     -0.000 3/14196
  890.0      4.880      4.880      5.273      0.000     -0.000     -0.000 3/14196
  900.0      4.828      4.828      5.217      0.000     -0.000     -0.000 3/14196
  910.0      4.777      4.777      5.162      0.000     -0.000     -0.000 3/14196
  920.0      4.727      4.727      5.108      0.000     -0.000     -0.000 3/14196
  930.0      4.678      4.678      5.056      0.000     -0.000     -0.000 3/14196
  940.0      4.631      4.631      5.004      0.000     -0.000     -0.000 3/14196
  950.0      4.584      4.584      4.954      0.000     -0.000     -0.000 3/14196
  960.0      4.538      4.538      4.904      0.000     -0.000     -0.000 3/14196
  970.0      4.493      4.493      4.856      0.000     -0.000     -0.000 3/14196
  980.0      4.449      4.449      4.808      0.000     -0.000     -0.000 3/14196
  990.0      4.405      4.405      4.762      0.000     -0.000     -0.000 3/14196
 1000.0      4.363      4.363      4.716      0.000     -0.000     -0.000 3/14196

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

