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

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

-------------------------[time 2026-01-08 19:08:04]-------------------------
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.495121562572470    0.000000000000000    0.000000000000000
  b   -2.247560781286235    3.892889466286960    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.257937420000000
Atomic positions (fractional):
   *1 Li  0.33333333333333  0.66666666666667  0.99831627125308   6.941
    2 Li  0.66666666666667  0.33333333333333  0.49831627125308   6.941
   *3 I   0.33333333333333  0.66666666666667  0.37668372874692 126.904
    4 I   0.66666666666667  0.33333333333333  0.87668372874692 126.904
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.495121562572470    0.000000000000000    0.000000000000000
  b   -2.247560781286235    3.892889466286960    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.257937420000000
Atomic positions (fractional):
   *1 Li  0.33333333333333  0.66666666666667  0.99831627125308   6.941 > 1
    2 Li  0.66666666666667  0.33333333333333  0.49831627125308   6.941 > 2
   *3 I   0.33333333333333  0.66666666666667  0.37668372874692 126.904 > 3
    4 I   0.66666666666667  0.33333333333333  0.87668372874692 126.904 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   17.980486250289879    0.000000000000000    0.000000000000000
  b   -8.990243125144939   15.571557865147840    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.515874840000000
Atomic positions (fractional):
   *1 Li  0.08333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    3 Li  0.58333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    4 Li  0.83333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    5 Li  0.08333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    6 Li  0.33333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    7 Li  0.58333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    8 Li  0.83333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    9 Li  0.08333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   10 Li  0.33333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   11 Li  0.58333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   12 Li  0.83333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   13 Li  0.08333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   14 Li  0.33333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   15 Li  0.58333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   16 Li  0.83333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   17 Li  0.08333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   18 Li  0.33333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   19 Li  0.58333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   20 Li  0.83333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   21 Li  0.08333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   22 Li  0.33333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   23 Li  0.58333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   24 Li  0.83333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   25 Li  0.08333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   26 Li  0.33333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   27 Li  0.58333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   28 Li  0.83333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   29 Li  0.08333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   30 Li  0.33333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   31 Li  0.58333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   32 Li  0.83333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   33 Li  0.16666666666667  0.08333333333333  0.24915813562654   6.941 > 2
   34 Li  0.41666666666667  0.08333333333333  0.24915813562654   6.941 > 2
   35 Li  0.66666666666667  0.08333333333333  0.24915813562654   6.941 > 2
   36 Li  0.91666666666667  0.08333333333333  0.24915813562654   6.941 > 2
   37 Li  0.16666666666667  0.33333333333333  0.24915813562654   6.941 > 2
   38 Li  0.41666666666667  0.33333333333333  0.24915813562654   6.941 > 2
   39 Li  0.66666666666667  0.33333333333333  0.24915813562654   6.941 > 2
   40 Li  0.91666666666667  0.33333333333333  0.24915813562654   6.941 > 2
   41 Li  0.16666666666667  0.58333333333333  0.24915813562654   6.941 > 2
   42 Li  0.41666666666667  0.58333333333333  0.24915813562654   6.941 > 2
   43 Li  0.66666666666667  0.58333333333333  0.24915813562654   6.941 > 2
   44 Li  0.91666666666667  0.58333333333333  0.24915813562654   6.941 > 2
   45 Li  0.16666666666667  0.83333333333333  0.24915813562654   6.941 > 2
   46 Li  0.41666666666667  0.83333333333333  0.24915813562654   6.941 > 2
   47 Li  0.66666666666667  0.83333333333333  0.24915813562654   6.941 > 2
   48 Li  0.91666666666667  0.83333333333333  0.24915813562654   6.941 > 2
   49 Li  0.16666666666667  0.08333333333333  0.74915813562654   6.941 > 2
   50 Li  0.41666666666667  0.08333333333333  0.74915813562654   6.941 > 2
   51 Li  0.66666666666667  0.08333333333333  0.74915813562654   6.941 > 2
   52 Li  0.91666666666667  0.08333333333333  0.74915813562654   6.941 > 2
   53 Li  0.16666666666667  0.33333333333333  0.74915813562654   6.941 > 2
   54 Li  0.41666666666667  0.33333333333333  0.74915813562654   6.941 > 2
   55 Li  0.66666666666667  0.33333333333333  0.74915813562654   6.941 > 2
   56 Li  0.91666666666667  0.33333333333333  0.74915813562654   6.941 > 2
   57 Li  0.16666666666667  0.58333333333333  0.74915813562654   6.941 > 2
   58 Li  0.41666666666667  0.58333333333333  0.74915813562654   6.941 > 2
   59 Li  0.66666666666667  0.58333333333333  0.74915813562654   6.941 > 2
   60 Li  0.91666666666667  0.58333333333333  0.74915813562654   6.941 > 2
   61 Li  0.16666666666667  0.83333333333333  0.74915813562654   6.941 > 2
   62 Li  0.41666666666667  0.83333333333333  0.74915813562654   6.941 > 2
   63 Li  0.66666666666667  0.83333333333333  0.74915813562654   6.941 > 2
   64 Li  0.91666666666667  0.83333333333333  0.74915813562654   6.941 > 2
  *65 I   0.08333333333333  0.16666666666667  0.18834186437346 126.904 > 3
   66 I   0.33333333333333  0.16666666666667  0.18834186437346 126.904 > 3
   67 I   0.58333333333333  0.16666666666667  0.18834186437346 126.904 > 3
   68 I   0.83333333333333  0.16666666666667  0.18834186437346 126.904 > 3
   69 I   0.08333333333333  0.41666666666667  0.18834186437346 126.904 > 3
   70 I   0.33333333333333  0.41666666666667  0.18834186437346 126.904 > 3
   71 I   0.58333333333333  0.41666666666667  0.18834186437346 126.904 > 3
   72 I   0.83333333333333  0.41666666666667  0.18834186437346 126.904 > 3
   73 I   0.08333333333333  0.66666666666667  0.18834186437346 126.904 > 3
   74 I   0.33333333333333  0.66666666666667  0.18834186437346 126.904 > 3
   75 I   0.58333333333333  0.66666666666667  0.18834186437346 126.904 > 3
   76 I   0.83333333333333  0.66666666666667  0.18834186437346 126.904 > 3
   77 I   0.08333333333333  0.91666666666667  0.18834186437346 126.904 > 3
   78 I   0.33333333333333  0.91666666666667  0.18834186437346 126.904 > 3
   79 I   0.58333333333333  0.91666666666667  0.18834186437346 126.904 > 3
   80 I   0.83333333333333  0.91666666666667  0.18834186437346 126.904 > 3
   81 I   0.08333333333333  0.16666666666667  0.68834186437346 126.904 > 3
   82 I   0.33333333333333  0.16666666666667  0.68834186437346 126.904 > 3
   83 I   0.58333333333333  0.16666666666667  0.68834186437346 126.904 > 3
   84 I   0.83333333333333  0.16666666666667  0.68834186437346 126.904 > 3
   85 I   0.08333333333333  0.41666666666667  0.68834186437346 126.904 > 3
   86 I   0.33333333333333  0.41666666666667  0.68834186437346 126.904 > 3
   87 I   0.58333333333333  0.41666666666667  0.68834186437346 126.904 > 3
   88 I   0.83333333333333  0.41666666666667  0.68834186437346 126.904 > 3
   89 I   0.08333333333333  0.66666666666667  0.68834186437346 126.904 > 3
   90 I   0.33333333333333  0.66666666666667  0.68834186437346 126.904 > 3
   91 I   0.58333333333333  0.66666666666667  0.68834186437346 126.904 > 3
   92 I   0.83333333333333  0.66666666666667  0.68834186437346 126.904 > 3
   93 I   0.08333333333333  0.91666666666667  0.68834186437346 126.904 > 3
   94 I   0.33333333333333  0.91666666666667  0.68834186437346 126.904 > 3
   95 I   0.58333333333333  0.91666666666667  0.68834186437346 126.904 > 3
   96 I   0.83333333333333  0.91666666666667  0.68834186437346 126.904 > 3
   97 I   0.16666666666667  0.08333333333333  0.43834186437346 126.904 > 4
   98 I   0.41666666666667  0.08333333333333  0.43834186437346 126.904 > 4
   99 I   0.66666666666667  0.08333333333333  0.43834186437346 126.904 > 4
  100 I   0.91666666666667  0.08333333333333  0.43834186437346 126.904 > 4
  101 I   0.16666666666667  0.33333333333333  0.43834186437346 126.904 > 4
  102 I   0.41666666666667  0.33333333333333  0.43834186437346 126.904 > 4
  103 I   0.66666666666667  0.33333333333333  0.43834186437346 126.904 > 4
  104 I   0.91666666666667  0.33333333333333  0.43834186437346 126.904 > 4
  105 I   0.16666666666667  0.58333333333333  0.43834186437346 126.904 > 4
  106 I   0.41666666666667  0.58333333333333  0.43834186437346 126.904 > 4
  107 I   0.66666666666667  0.58333333333333  0.43834186437346 126.904 > 4
  108 I   0.91666666666667  0.58333333333333  0.43834186437346 126.904 > 4
  109 I   0.16666666666667  0.83333333333333  0.43834186437346 126.904 > 4
  110 I   0.41666666666667  0.83333333333333  0.43834186437346 126.904 > 4
  111 I   0.66666666666667  0.83333333333333  0.43834186437346 126.904 > 4
  112 I   0.91666666666667  0.83333333333333  0.43834186437346 126.904 > 4
  113 I   0.16666666666667  0.08333333333333  0.93834186437346 126.904 > 4
  114 I   0.41666666666667  0.08333333333333  0.93834186437346 126.904 > 4
  115 I   0.66666666666667  0.08333333333333  0.93834186437346 126.904 > 4
  116 I   0.91666666666667  0.08333333333333  0.93834186437346 126.904 > 4
  117 I   0.16666666666667  0.33333333333333  0.93834186437346 126.904 > 4
  118 I   0.41666666666667  0.33333333333333  0.93834186437346 126.904 > 4
  119 I   0.66666666666667  0.33333333333333  0.93834186437346 126.904 > 4
  120 I   0.91666666666667  0.33333333333333  0.93834186437346 126.904 > 4
  121 I   0.16666666666667  0.58333333333333  0.93834186437346 126.904 > 4
  122 I   0.41666666666667  0.58333333333333  0.93834186437346 126.904 > 4
  123 I   0.66666666666667  0.58333333333333  0.93834186437346 126.904 > 4
  124 I   0.91666666666667  0.58333333333333  0.93834186437346 126.904 > 4
  125 I   0.16666666666667  0.83333333333333  0.93834186437346 126.904 > 4
  126 I   0.41666666666667  0.83333333333333  0.93834186437346 126.904 > 4
  127 I   0.66666666666667  0.83333333333333  0.93834186437346 126.904 > 4
  128 I   0.91666666666667  0.83333333333333  0.93834186437346 126.904 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.2940445    0.0000000    0.0000000
            0.0000000    3.2940445    0.0000000
            0.0000000    0.0000000    3.3276319
-------------------------- Born effective charges --------------------------
    1 Li    1.0201279    0.0000000    0.0000000
            0.0000000    1.0201279    0.0000000
            0.0000000    0.0000000    1.0443918
    2 Li    1.0201279    0.0000000    0.0000000
            0.0000000    1.0201279    0.0000000
            0.0000000    0.0000000    1.0443918
    3 I    -1.0201279   -0.0000000    0.0000000
            0.0000000   -1.0201279    0.0000000
            0.0000000    0.0000000   -1.0443918
    4 I    -1.0201279   -0.0000000    0.0000000
            0.0000000   -1.0201279    0.0000000
            0.0000000    0.0000000   -1.0443918
----------------------------------------------------------------------------
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: 156
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.009
Solver_block: 80 / 80
 - Time: 0.160
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.171
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 156
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 19:08:08]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [4 4 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3mc (186)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.495121562572470    0.000000000000000    0.000000000000000
  b   -2.247560781286235    3.892889466286960    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.257937420000000
Atomic positions (fractional):
    1 Li  0.33333333333333  0.66666666666667  0.99831627125308   6.941
    2 Li  0.66666666666667  0.33333333333333  0.49831627125308   6.941
    3 I   0.33333333333333  0.66666666666667  0.37668372874692 126.904
    4 I   0.66666666666667  0.33333333333333  0.87668372874692 126.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   17.980486250289879    0.000000000000000    0.000000000000000
  b   -8.990243125144939   15.571557865147840    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.515874840000000
Atomic positions (fractional):
    1 Li  0.08333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    3 Li  0.58333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    4 Li  0.83333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    5 Li  0.08333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    6 Li  0.33333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    7 Li  0.58333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    8 Li  0.83333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    9 Li  0.08333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   10 Li  0.33333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   11 Li  0.58333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   12 Li  0.83333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   13 Li  0.08333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   14 Li  0.33333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   15 Li  0.58333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   16 Li  0.83333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   17 Li  0.08333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   18 Li  0.33333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   19 Li  0.58333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   20 Li  0.83333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   21 Li  0.08333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   22 Li  0.33333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   23 Li  0.58333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   24 Li  0.83333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   25 Li  0.08333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   26 Li  0.33333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   27 Li  0.58333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   28 Li  0.83333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   29 Li  0.08333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   30 Li  0.33333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   31 Li  0.58333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   32 Li  0.83333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   33 Li  0.16666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   34 Li  0.41666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   35 Li  0.66666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   36 Li  0.91666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   37 Li  0.16666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   38 Li  0.41666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   39 Li  0.66666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   40 Li  0.91666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   41 Li  0.16666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   42 Li  0.41666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   43 Li  0.66666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   44 Li  0.91666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   45 Li  0.16666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   46 Li  0.41666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   47 Li  0.66666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   48 Li  0.91666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   49 Li  0.16666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   50 Li  0.41666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   51 Li  0.66666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   52 Li  0.91666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   53 Li  0.16666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   54 Li  0.41666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   55 Li  0.66666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   56 Li  0.91666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   57 Li  0.16666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   58 Li  0.41666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   59 Li  0.66666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   60 Li  0.91666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   61 Li  0.16666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   62 Li  0.41666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   63 Li  0.66666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   64 Li  0.91666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   65 I   0.08333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   66 I   0.33333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   67 I   0.58333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   68 I   0.83333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   69 I   0.08333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   70 I   0.33333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   71 I   0.58333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   72 I   0.83333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   73 I   0.08333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   74 I   0.33333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   75 I   0.58333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   76 I   0.83333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   77 I   0.08333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   78 I   0.33333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   79 I   0.58333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   80 I   0.83333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   81 I   0.08333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   82 I   0.33333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   83 I   0.58333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   84 I   0.83333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   85 I   0.08333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   86 I   0.33333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   87 I   0.58333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   88 I   0.83333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   89 I   0.08333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   90 I   0.33333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   91 I   0.58333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   92 I   0.83333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   93 I   0.08333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   94 I   0.33333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   95 I   0.58333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   96 I   0.83333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   97 I   0.16666666666667  0.08333333333333  0.43834186437346 126.904 > 97
   98 I   0.41666666666667  0.08333333333333  0.43834186437346 126.904 > 97
   99 I   0.66666666666667  0.08333333333333  0.43834186437346 126.904 > 97
  100 I   0.91666666666667  0.08333333333333  0.43834186437346 126.904 > 97
  101 I   0.16666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  102 I   0.41666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  103 I   0.66666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  104 I   0.91666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  105 I   0.16666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  106 I   0.41666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  107 I   0.66666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  108 I   0.91666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  109 I   0.16666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  110 I   0.41666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  111 I   0.66666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  112 I   0.91666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  113 I   0.16666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  114 I   0.41666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  115 I   0.66666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  116 I   0.91666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  117 I   0.16666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  118 I   0.41666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  119 I   0.66666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  120 I   0.91666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  121 I   0.16666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  122 I   0.41666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  123 I   0.66666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  124 I   0.91666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  125 I   0.16666666666667  0.83333333333333  0.93834186437346 126.904 > 97
  126 I   0.41666666666667  0.83333333333333  0.93834186437346 126.904 > 97
  127 I   0.66666666666667  0.83333333333333  0.93834186437346 126.904 > 97
  128 I   0.91666666666667  0.83333333333333  0.93834186437346 126.904 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.2940445    0.0000000    0.0000000
            0.0000000    3.2940445    0.0000000
            0.0000000    0.0000000    3.3276319
-------------------------- Born effective charges --------------------------
    1 Li    1.0201279    0.0000000    0.0000000
            0.0000000    1.0201279    0.0000000
            0.0000000    0.0000000    1.0443918
    2 Li    1.0201279    0.0000000    0.0000000
            0.0000000    1.0201279    0.0000000
            0.0000000    0.0000000    1.0443918
    3 I    -1.0201279   -0.0000000    0.0000000
            0.0000000   -1.0201279    0.0000000
            0.0000000    0.0000000   -1.0443918
    4 I    -1.0201279   -0.0000000    0.0000000
            0.0000000   -1.0201279    0.0000000
            0.0000000    0.0000000   -1.0443918
----------------------------------------------------------------------------
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.00000041 (zzz) -0.00000041 (zzz) -0.00000041 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 19:08:13]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 19:08:14]-------------------------
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.495121562572470    0.000000000000000    0.000000000000000
  b   -2.247560781286235    3.892889466286960    0.000000000000000
  c    0.000000000000000    0.000000000000000    7.257937420000000
Atomic positions (fractional):
    1 Li  0.33333333333333  0.66666666666667  0.99831627125308   6.941
    2 Li  0.66666666666667  0.33333333333333  0.49831627125308   6.941
    3 I   0.33333333333333  0.66666666666667  0.37668372874692 126.904
    4 I   0.66666666666667  0.33333333333333  0.87668372874692 126.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   17.980486250289879    0.000000000000000    0.000000000000000
  b   -8.990243125144939   15.571557865147840    0.000000000000000
  c    0.000000000000000    0.000000000000000   14.515874840000000
Atomic positions (fractional):
    1 Li  0.08333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    3 Li  0.58333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    4 Li  0.83333333333333  0.16666666666667  0.49915813562654   6.941 > 1
    5 Li  0.08333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    6 Li  0.33333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    7 Li  0.58333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    8 Li  0.83333333333333  0.41666666666667  0.49915813562654   6.941 > 1
    9 Li  0.08333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   10 Li  0.33333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   11 Li  0.58333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   12 Li  0.83333333333333  0.66666666666667  0.49915813562654   6.941 > 1
   13 Li  0.08333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   14 Li  0.33333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   15 Li  0.58333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   16 Li  0.83333333333333  0.91666666666667  0.49915813562654   6.941 > 1
   17 Li  0.08333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   18 Li  0.33333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   19 Li  0.58333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   20 Li  0.83333333333333  0.16666666666667  0.99915813562654   6.941 > 1
   21 Li  0.08333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   22 Li  0.33333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   23 Li  0.58333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   24 Li  0.83333333333333  0.41666666666667  0.99915813562654   6.941 > 1
   25 Li  0.08333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   26 Li  0.33333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   27 Li  0.58333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   28 Li  0.83333333333333  0.66666666666667  0.99915813562654   6.941 > 1
   29 Li  0.08333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   30 Li  0.33333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   31 Li  0.58333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   32 Li  0.83333333333333  0.91666666666667  0.99915813562654   6.941 > 1
   33 Li  0.16666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   34 Li  0.41666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   35 Li  0.66666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   36 Li  0.91666666666667  0.08333333333333  0.24915813562654   6.941 > 33
   37 Li  0.16666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   38 Li  0.41666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   39 Li  0.66666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   40 Li  0.91666666666667  0.33333333333333  0.24915813562654   6.941 > 33
   41 Li  0.16666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   42 Li  0.41666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   43 Li  0.66666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   44 Li  0.91666666666667  0.58333333333333  0.24915813562654   6.941 > 33
   45 Li  0.16666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   46 Li  0.41666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   47 Li  0.66666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   48 Li  0.91666666666667  0.83333333333333  0.24915813562654   6.941 > 33
   49 Li  0.16666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   50 Li  0.41666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   51 Li  0.66666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   52 Li  0.91666666666667  0.08333333333333  0.74915813562654   6.941 > 33
   53 Li  0.16666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   54 Li  0.41666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   55 Li  0.66666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   56 Li  0.91666666666667  0.33333333333333  0.74915813562654   6.941 > 33
   57 Li  0.16666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   58 Li  0.41666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   59 Li  0.66666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   60 Li  0.91666666666667  0.58333333333333  0.74915813562654   6.941 > 33
   61 Li  0.16666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   62 Li  0.41666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   63 Li  0.66666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   64 Li  0.91666666666667  0.83333333333333  0.74915813562654   6.941 > 33
   65 I   0.08333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   66 I   0.33333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   67 I   0.58333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   68 I   0.83333333333333  0.16666666666667  0.18834186437346 126.904 > 65
   69 I   0.08333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   70 I   0.33333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   71 I   0.58333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   72 I   0.83333333333333  0.41666666666667  0.18834186437346 126.904 > 65
   73 I   0.08333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   74 I   0.33333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   75 I   0.58333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   76 I   0.83333333333333  0.66666666666667  0.18834186437346 126.904 > 65
   77 I   0.08333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   78 I   0.33333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   79 I   0.58333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   80 I   0.83333333333333  0.91666666666667  0.18834186437346 126.904 > 65
   81 I   0.08333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   82 I   0.33333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   83 I   0.58333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   84 I   0.83333333333333  0.16666666666667  0.68834186437346 126.904 > 65
   85 I   0.08333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   86 I   0.33333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   87 I   0.58333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   88 I   0.83333333333333  0.41666666666667  0.68834186437346 126.904 > 65
   89 I   0.08333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   90 I   0.33333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   91 I   0.58333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   92 I   0.83333333333333  0.66666666666667  0.68834186437346 126.904 > 65
   93 I   0.08333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   94 I   0.33333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   95 I   0.58333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   96 I   0.83333333333333  0.91666666666667  0.68834186437346 126.904 > 65
   97 I   0.16666666666667  0.08333333333333  0.43834186437346 126.904 > 97
   98 I   0.41666666666667  0.08333333333333  0.43834186437346 126.904 > 97
   99 I   0.66666666666667  0.08333333333333  0.43834186437346 126.904 > 97
  100 I   0.91666666666667  0.08333333333333  0.43834186437346 126.904 > 97
  101 I   0.16666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  102 I   0.41666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  103 I   0.66666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  104 I   0.91666666666667  0.33333333333333  0.43834186437346 126.904 > 97
  105 I   0.16666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  106 I   0.41666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  107 I   0.66666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  108 I   0.91666666666667  0.58333333333333  0.43834186437346 126.904 > 97
  109 I   0.16666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  110 I   0.41666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  111 I   0.66666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  112 I   0.91666666666667  0.83333333333333  0.43834186437346 126.904 > 97
  113 I   0.16666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  114 I   0.41666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  115 I   0.66666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  116 I   0.91666666666667  0.08333333333333  0.93834186437346 126.904 > 97
  117 I   0.16666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  118 I   0.41666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  119 I   0.66666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  120 I   0.91666666666667  0.33333333333333  0.93834186437346 126.904 > 97
  121 I   0.16666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  122 I   0.41666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  123 I   0.66666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  124 I   0.91666666666667  0.58333333333333  0.93834186437346 126.904 > 97
  125 I   0.16666666666667  0.83333333333333  0.93834186437346 126.904 > 97
  126 I   0.41666666666667  0.83333333333333  0.93834186437346 126.904 > 97
  127 I   0.66666666666667  0.83333333333333  0.93834186437346 126.904 > 97
  128 I   0.91666666666667  0.83333333333333  0.93834186437346 126.904 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.2940445    0.0000000    0.0000000
            0.0000000    3.2940445    0.0000000
            0.0000000    0.0000000    3.3276319
-------------------------- Born effective charges --------------------------
    1 Li    1.0201279    0.0000000    0.0000000
            0.0000000    1.0201279    0.0000000
            0.0000000    0.0000000    1.0443918
    2 Li    1.0201279    0.0000000    0.0000000
            0.0000000    1.0201279    0.0000000
            0.0000000    0.0000000    1.0443918
    3 I    -1.0201279   -0.0000000    0.0000000
            0.0000000   -1.0201279    0.0000000
            0.0000000    0.0000000   -1.0443918
    4 I    -1.0201279   -0.0000000    0.0000000
            0.0000000   -1.0201279    0.0000000
            0.0000000    0.0000000   -1.0443918
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000041 (zzz) -0.00000041 (zzz) -0.00000041 (zzz)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 13 13 7 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.83, Number of G-points: 311, Lambda: 0.16
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: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
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.872   (   0.000   -0.000    0.000)    0.000
   0.872   (  -0.000    0.000   -0.000)    0.000
   2.384   (   0.000    0.000    0.000)    0.000
   7.926   (   0.000    0.000    0.000)    0.000
   8.080   (  -0.000   -0.000   -0.000)    0.000
   8.080   (   0.000   -0.000    0.000)    0.000
   8.203   (   0.000    0.000   -0.000)    0.000
   8.203   (  -0.000    0.000    0.000)    0.000
   9.262   (   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: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.207   (   8.917    5.148    0.000)   10.297
   0.220   (   9.764    5.637    0.000)   11.274
   0.500   (  21.366   12.335    0.000)   24.671
   0.895   (   2.012    1.161    0.000)    2.323
   1.008   (  10.600    6.120    0.000)   12.240
   2.369   (  -1.277   -0.738    0.000)    1.475
   7.920   (  -0.478   -0.276    0.000)    0.552
   8.038   (  -2.453   -1.417    0.000)    2.833
   8.106   (   1.780    1.028    0.000)    2.056
   8.228   (   0.891    0.514    0.000)    1.029
   9.231   (  -2.379   -1.374    0.000)    2.747
  10.014   (  -1.846   -1.066    0.000)    2.131
======================= Grid point 2 (3/84) =======================
q-point: ( 0.15  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.402   (   8.028    4.635    0.000)    9.270
   0.442   (   9.696    5.598    0.000)   11.196
   0.964   (   3.888    2.245    0.000)    4.489
   0.967   (  19.370   11.183    0.000)   22.367
   1.286   (  12.532    7.236    0.000)   14.471
   2.329   (  -2.036   -1.176    0.000)    2.351
   7.906   (  -0.779   -0.450    0.000)    0.900
   8.000   (  -0.764   -0.441    0.000)    0.882
   8.129   (  -0.398   -0.230    0.000)    0.460
   8.217   (  -1.809   -1.044    0.000)    2.088
   9.181   (  -1.281   -0.740    0.000)    1.480
   9.962   (  -2.481   -1.432    0.000)    2.865
======================= Grid point 3 (4/84) =======================
q-point: ( 0.23  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.568   (   6.370    3.678    0.000)    7.355
   0.662   (   9.550    5.514    0.000)   11.028
   1.062   (   4.402    2.541    0.000)    5.083
   1.375   (  16.246    9.380    0.000)   18.759
   1.540   (   9.145    5.280    0.000)   10.560
   2.286   (  -1.475   -0.852    0.000)    1.703
   7.886   (  -0.890   -0.514    0.000)    1.028
   7.993   (  -0.267   -0.154    0.000)    0.309
   8.080   (  -3.468   -2.002    0.000)    4.004
   8.161   (  -2.513   -1.451    0.000)    2.901
   9.187   (   1.587    0.916    0.000)    1.833
   9.906   (  -2.449   -1.414    0.000)    2.827
======================= Grid point 4 (5/84) =======================
q-point: ( 0.31  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.689   (   4.295    2.480    0.000)    4.959
   0.875   (   9.046    5.223    0.000)   10.446
   1.151   (   3.145    1.816    0.000)    3.631
   1.657   (   4.841    2.795    0.000)    5.590
   1.731   (  10.985    6.342    0.000)   12.684
   2.276   (   0.786    0.454    0.000)    0.907
   7.867   (  -0.767   -0.443    0.000)    0.886
   7.973   (  -1.401   -0.809    0.000)    1.618
   8.004   (  -2.652   -1.531    0.000)    3.062
   8.123   (  -0.787   -0.454    0.000)    0.909
   9.240   (   2.947    1.701    0.000)    3.403
   9.842   (  -3.378   -1.950    0.000)    3.900
======================= Grid point 5 (6/84) =======================
q-point: ( 0.38  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.765   (   2.396    1.383    0.000)    2.766
   1.068   (   7.556    4.362    0.000)    8.725
   1.200   (   1.231    0.710    0.000)    1.421
   1.677   (  -1.432   -0.827    0.000)    1.653
   1.967   (   8.567    4.946    0.000)    9.892
   2.312   (   1.933    1.116    0.000)    2.232
   7.857   (  -0.050   -0.029    0.000)    0.057
   7.942   (  -1.127   -0.651    0.000)    1.302
   7.965   (  -0.921   -0.532    0.000)    1.064
   8.115   (  -0.081   -0.047    0.000)    0.094
   9.322   (   4.316    2.492    0.000)    4.983
   9.743   (  -5.342   -3.084    0.000)    6.168
======================= Grid point 6 (7/84) =======================
q-point: ( 0.46  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.801   (   0.791    0.457    0.000)    0.914
   1.201   (   3.505    2.024    0.000)    4.048
   1.214   (   0.149    0.086    0.000)    0.172
   1.632   (  -1.654   -0.955    0.000)    1.909
   2.105   (   3.234    1.867    0.000)    3.735
   2.345   (   0.759    0.438    0.000)    0.876
   7.860   (   0.207    0.120    0.000)    0.239
   7.924   (  -0.450   -0.260    0.000)    0.519
   7.954   (  -0.192   -0.111    0.000)    0.222
   8.116   (   0.096    0.055    0.000)    0.111
   9.435   (   5.156    2.977    0.000)    5.954
   9.609   (  -5.739   -3.314    0.000)    6.627
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.378   (   5.855   10.140    0.000)   11.709
   0.382   (   5.563    9.635    0.000)   11.126
   0.829   (  10.472   18.137    0.000)   20.943
   0.990   (   3.691    6.392    0.000)    7.381
   1.173   (   6.577   11.392    0.000)   13.155
   2.341   (  -1.117   -1.934    0.000)    2.234
   7.910   (  -0.406   -0.703    0.000)    0.812
   8.016   (  -1.651   -2.859    0.000)    3.302
   8.058   (  -0.457   -0.792    0.000)    0.915
   8.301   (   1.872    3.242    0.000)    3.744
   9.192   (  -1.199   -2.078    0.000)    2.399
   9.963   (  -2.108   -3.651    0.000)    4.216
======================= Grid point 15 (9/84) =======================
q-point: ( 0.15  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.570   (   3.991   10.895    0.000)   11.603
   0.582   (   7.293    8.190    0.000)   10.967
   1.065   (   2.807    7.336    0.000)    7.855
   1.234   (  12.081   14.390    0.000)   18.789
   1.419   (  10.259    6.453    0.000)   12.120
   2.298   (  -1.459   -1.794    0.000)    2.312
   7.894   (  -0.687   -0.619    0.000)    0.924
   7.968   (   0.451   -3.674    0.000)    3.702
   8.046   (  -0.428   -1.555    0.000)    1.613
   8.314   (  -4.934    4.894    0.000)    6.949
   9.169   (   0.384   -0.401    0.000)    0.555
   9.886   (  -1.089   -5.820    0.000)    5.921
======================= Grid point 16 (10/84) =======================
q-point: ( 0.23  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.722   (   1.169   10.843    0.000)   10.906
   0.789   (   7.885    7.106    0.000)   10.615
   1.170   (   0.950    7.728    0.000)    7.787
   1.561   (  10.799    8.171    0.000)   13.542
   1.623   (   8.370    4.062    0.000)    9.304
   2.268   (  -0.164   -0.821    0.000)    0.837
   7.875   (  -0.817   -0.512    0.000)    0.964
   7.951   (   1.640   -4.039    0.000)    4.360
   8.015   (  -2.249    0.262    0.000)    2.264
   8.251   (  -7.383    5.197    0.000)    9.028
   9.195   (   2.638   -0.231    0.000)    2.648
   9.815   (   0.173   -7.379    0.000)    7.381
======================= Grid point 17 (11/84) =======================
q-point: ( 0.31  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.824   (  -1.365   10.140    0.000)   10.232
   0.987   (   7.539    6.013    0.000)    9.644
   1.244   (  -1.067    7.195    0.000)    7.274
   1.672   (   1.315   -0.729    0.000)    1.503
   1.860   (  10.039    6.198    0.000)   11.798
   2.282   (   1.993    0.061    0.000)    1.994
   7.861   (  -0.414    0.098    0.000)    0.425
   7.941   (   1.080   -2.490    0.000)    2.714
   7.979   (  -2.562    1.007    0.000)    2.753
   8.187   (  -5.217    4.196    0.000)    6.695
   9.251   (   4.520   -0.657    0.000)    4.567
   9.735   (  -0.710   -8.337    0.000)    8.367
======================= Grid point 18 (12/84) =======================
q-point: ( 0.38  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.879   (  -3.268    9.277    0.000)    9.836
   1.154   (   5.419    4.399    0.000)    6.980
   1.276   (  -2.949    6.513    0.000)    7.149
   1.648   (  -1.593   -1.945    0.000)    2.514
   2.048   (   6.157    3.091    0.000)    6.889
   2.319   (   2.078   -0.507    0.000)    2.139
   7.863   (  -0.028    0.646    0.000)    0.647
   7.938   (   0.610   -1.264    0.000)    1.404
   7.950   (  -1.858    0.964    0.000)    2.093
   8.156   (  -2.712    3.016    0.000)    4.056
   9.340   (   6.616   -0.723    0.000)    6.656
   9.623   (  -2.527   -8.464    0.000)    8.833
======================= Grid point 19 (13/84) =======================
q-point: ( 0.46  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.897   (  -4.842    8.387    0.000)    9.685
   1.231   (  -0.576    0.998    0.000)    1.153
   1.280   (  -3.556    6.160    0.000)    7.112
   1.616   (   0.335   -0.580    0.000)    0.670
   2.118   (   0.304   -0.526    0.000)    0.607
   2.335   (   0.859   -1.488    0.000)    1.719
   7.867   (  -0.302    0.523    0.000)    0.604
   7.937   (  -0.872    1.511    0.000)    1.744
   7.938   (   0.649   -1.124    0.000)    1.298
   8.149   (  -1.558    2.699    0.000)    3.116
   9.431   (   2.415   -4.183    0.000)    4.830
   9.520   (   2.545   -4.408    0.000)    5.089
======================= Grid point 29 (14/84) =======================
q-point: ( 0.15  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.754   (   5.262    9.115    0.000)   10.525
   0.784   (   5.757    9.972    0.000)   11.514
   1.214   (   4.175    7.232    0.000)    8.351
   1.514   (   7.811   13.529    0.000)   15.622
   1.543   (   3.399    5.887    0.000)    6.798
   2.266   (  -0.743   -1.287    0.000)    1.486
   7.884   (  -0.256   -0.444    0.000)    0.513
   7.907   (  -1.029   -1.782    0.000)    2.058
   8.041   (  -0.040   -0.069    0.000)    0.079
   8.362   (  -0.105   -0.182    0.000)    0.210
   9.167   (   0.116    0.201    0.000)    0.232
   9.759   (  -3.696   -6.401    0.000)    7.391
======================= Grid point 30 (15/84) =======================
q-point: ( 0.23  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.935   (   5.434    8.160    0.000)    9.803
   0.950   (   1.862   10.678    0.000)   10.839
   1.338   (   0.870    8.647    0.000)    8.691
   1.628   (   2.671    0.421    0.000)    2.704
   1.767   (   8.507    8.678    0.000)   12.152
   2.256   (   0.654   -0.447    0.000)    0.792
   7.873   (  -0.677    0.337    0.000)    0.757
   7.900   (   1.269   -1.210    0.000)    1.754
   8.031   (  -1.620    0.728    0.000)    1.776
   8.320   (  -4.805    1.431    0.000)    5.014
   9.180   (   1.994   -1.074    0.000)    2.264
   9.646   (  -0.818   -8.825    0.000)    8.863
======================= Grid point 31 (16/84) =======================
q-point: ( 0.31  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.044   (  -2.489   11.178    0.000)   11.451
   1.111   (   5.593    6.485    0.000)    8.563
   1.410   (  -2.638    8.838    0.000)    9.224
   1.640   (   0.203   -2.301    0.000)    2.310
   1.969   (   6.981    4.870    0.000)    8.511
   2.274   (   1.951   -0.977    0.000)    2.182
   7.872   (  -0.148    0.885    0.000)    0.898
   7.912   (   1.009   -0.581    0.000)    1.164
   8.003   (  -2.447    1.000    0.000)    2.644
   8.250   (  -4.639    1.658    0.000)    4.927
   9.221   (   4.780   -2.007    0.000)    5.184
   9.547   (  -0.299   -9.829    0.000)    9.834
======================= Grid point 32 (17/84) =======================
q-point: ( 0.38  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.087   (  -5.016   10.762    0.000)   11.874
   1.241   (   2.404    4.281    0.000)    4.910
   1.426   (  -4.482    7.448    0.000)    8.693
   1.611   (  -0.731   -1.272    0.000)    1.467
   2.090   (   3.124    1.321    0.000)    3.392
   2.291   (   1.797   -2.326    0.000)    2.940
   7.880   (  -0.246    1.089    0.000)    1.117
   7.919   (   0.434   -0.519    0.000)    0.676
   7.976   (  -1.705    1.319    0.000)    2.156
   8.206   (  -2.258    1.643    0.000)    2.793
   9.300   (   6.698   -3.007    0.000)    7.342
   9.442   (  -0.541   -8.926    0.000)    8.943
======================= Grid point 43 (18/84) =======================
q-point: ( 0.23  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.092   (   4.504    7.801    0.000)    9.008
   1.134   (   4.293    7.436    0.000)    8.586
   1.493   (   3.825    6.625    0.000)    7.650
   1.612   (  -0.985   -1.706    0.000)    1.970
   1.929   (   4.303    7.454    0.000)    8.607
   2.249   (  -0.142   -0.246    0.000)    0.284
   7.881   (   0.180    0.312    0.000)    0.360
   7.898   (   0.474    0.821    0.000)    0.948
   8.032   (  -0.359   -0.622    0.000)    0.718
   8.307   (  -1.404   -2.432    0.000)    2.809
   9.165   (  -0.004   -0.006    0.000)    0.007
   9.490   (  -3.912   -6.776    0.000)    7.824
======================= Grid point 44 (19/84) =======================
q-point: ( 0.31  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 340
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.232   (   1.660    7.216    0.000)    7.404
   1.245   (   1.535    7.078    0.000)    7.243
   1.542   (  -3.441    0.662    0.000)    3.504
   1.613   (   1.535    3.003    0.000)    3.373
   2.054   (   3.525    3.747    0.000)    5.144
   2.242   (   0.657   -2.212    0.000)    2.308
   7.891   (   0.296    1.010    0.000)    1.052
   7.914   (   0.291    0.541    0.000)    0.614
   8.010   (  -1.374   -0.315    0.000)    1.410
   8.249   (  -2.324   -1.433    0.000)    2.730
   9.188   (   3.213   -0.907    0.000)    3.339
   9.374   (  -1.691   -7.294    0.000)    7.488
======================= Grid point 45 (20/84) =======================
q-point: ( 0.38  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.271   (  -4.097    7.096    0.000)    8.194
   1.319   (  -1.967    3.407    0.000)    3.934
   1.497   (   0.330   -0.571    0.000)    0.659
   1.645   (  -2.565    4.444    0.000)    5.131
   2.108   (  -0.407    0.704    0.000)    0.813
   2.234   (   1.836   -3.180    0.000)    3.672
   7.904   (  -0.662    1.147    0.000)    1.324
   7.916   (  -0.072    0.125    0.000)    0.144
   7.994   (  -0.215    0.372    0.000)    0.430
   8.219   (   0.031   -0.054    0.000)    0.062
   9.233   (   1.931   -3.345    0.000)    3.863
   9.303   (   2.623   -4.543    0.000)    5.245
======================= Grid point 58 (21/84) =======================
q-point: ( 0.31  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.347   (   1.596    2.765    0.000)    3.193
   1.356   (   3.114    5.394    0.000)    6.228
   1.479   (  -2.510   -4.348    0.000)    5.021
   1.700   (   1.653    2.864    0.000)    3.307
   2.120   (   1.657    2.870    0.000)    3.314
   2.191   (  -1.534   -2.656    0.000)    3.067
   7.913   (   0.616    1.067    0.000)    1.232
   7.923   (   0.086    0.149    0.000)    0.172
   7.995   (  -0.493   -0.854    0.000)    0.986
   8.216   (  -0.697   -1.208    0.000)    1.395
   9.195   (   0.879    1.523    0.000)    1.759
   9.261   (  -2.404   -4.163    0.000)    4.807
======================= Grid point 181 (22/84) =======================
q-point: ( 0.00  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.215   (   0.000   -0.000   10.656)   10.656
   0.215   (   0.000    0.000   10.656)   10.656
   0.519   (  -0.000    0.000   25.988)   25.988
   0.856   (   0.000    0.000   -1.713)    1.713
   0.856   (  -0.000    0.000   -1.713)    1.713
   2.321   (   0.000   -0.000   -6.349)    6.349
   8.087   (  -0.000   -0.000    0.694)    0.694
   8.087   (   0.000   -0.000    0.694)    0.694
   8.198   (   0.000    0.000   -0.524)    0.524
   8.198   (   0.000   -0.000   -0.524)    0.524
   9.295   (  -0.000   -0.000    3.282)    3.282
   9.847   (   0.000   -0.000   -2.806)    2.806
======================= Grid point 182 (23/84) =======================
q-point: ( 0.08  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.300   (   6.247    3.607    7.742)   10.582
   0.331   (   1.706    0.985   15.192)   15.319
   0.671   (  16.543    9.551   14.143)   23.768
   0.880   (   2.121    1.225   -1.562)    2.905
   0.983   (  10.259    5.923   -2.067)   12.025
   2.307   (  -1.186   -0.685   -6.247)    6.395
   8.020   (  -5.309   -3.065    4.318)    7.498
   8.048   (  -2.294   -1.324    1.017)    2.837
   8.139   (   0.212    0.123    2.762)    2.773
   8.219   (   0.720    0.416   -0.857)    1.195
   9.273   (  -1.565   -0.904    4.158)    4.533
   9.915   (   2.006    1.158   -6.718)    7.106
======================= Grid point 183 (24/84) =======================
q-point: ( 0.15  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.447   (   7.938    4.583    5.688)   10.787
   0.459   (   7.170    4.140    5.274)    9.816
   0.950   (   3.933    2.271   -1.365)    4.743
   1.076   (  17.581   10.150    8.922)   22.175
   1.256   (  12.347    7.128   -2.693)   14.509
   2.272   (  -1.671   -0.965   -5.723)    6.040
   7.936   (  -2.448   -1.413    2.466)    3.751
   8.011   (  -0.829   -0.479    1.132)    1.483
   8.139   (  -1.167   -0.674    0.833)    1.584
   8.207   (  -1.771   -1.022   -1.015)    2.283
   9.248   (  -0.164   -0.094    6.297)    6.300
   9.919   (  -0.919   -0.531   -4.012)    4.150
======================= Grid point 184 (25/84) =======================
q-point: ( 0.23  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.611   (   6.004    3.466    4.184)    8.097
   0.645   (   8.908    5.143    0.597)   10.304
   1.049   (   4.390    2.535   -1.335)    5.242
   1.445   (  13.481    7.783    4.065)   16.089
   1.510   (  10.263    5.926   -1.015)   11.895
   2.243   (  -0.533   -0.308   -4.198)    4.243
   7.891   (  -1.615   -0.932    0.380)    1.903
   8.002   (  -0.395   -0.228    0.854)    0.968
   8.079   (  -3.610   -2.084    0.117)    4.170
   8.153   (  -2.419   -1.396   -0.806)    2.907
   9.274   (   2.252    1.300    8.144)    8.549
   9.890   (  -1.501   -0.866   -1.538)    2.317
======================= Grid point 185 (26/84) =======================
q-point: ( 0.31  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.727   (   4.114    2.375    3.653)    5.992
   0.842   (   8.240    4.758   -2.075)    9.739
   1.138   (   3.157    1.823   -1.341)    3.884
   1.603   (   2.166    1.250   -5.950)    6.454
   1.788   (  12.093    6.982    5.626)   15.054
   2.260   (   2.061    1.190   -1.557)    2.844
   7.859   (  -1.217   -0.703   -0.834)    1.634
   7.980   (  -1.391   -0.803    0.727)    1.763
   8.004   (  -2.518   -1.454    0.232)    2.916
   8.116   (  -0.836   -0.483   -0.756)    1.226
   9.336   (   3.066    1.770    8.907)    9.585
   9.847   (  -2.527   -1.459    0.420)    2.948
======================= Grid point 186 (27/84) =======================
q-point: ( 0.38  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.800   (   2.308    1.333    3.384)    4.307
   1.016   (   6.819    3.937   -4.274)    8.959
   1.188   (   1.287    0.743   -1.260)    1.948
   1.595   (  -2.118   -1.223   -7.901)    8.271
   2.023   (   8.293    4.788    5.039)   10.820
   2.321   (   2.708    1.564    0.591)    3.183
   7.841   (  -0.278   -0.161   -1.337)    1.375
   7.950   (  -1.095   -0.632    0.795)    1.494
   7.967   (  -0.887   -0.512    0.202)    1.043
   8.106   (  -0.152   -0.088   -0.910)    0.927
   9.415   (   3.937    2.273    8.600)    9.728
   9.768   (  -4.391   -2.535    2.404)    5.611
======================= Grid point 187 (28/84) =======================
q-point: ( 0.46  0.00  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.835   (   0.760    0.439    3.255)    3.371
   1.137   (   3.232    1.866   -5.628)    6.753
   1.203   (   0.188    0.108   -1.136)    1.156
   1.540   (  -1.876   -1.083   -8.687)    8.953
   2.156   (   3.089    1.784    4.535)    5.769
   2.365   (   1.023    0.590    1.607)    1.994
   7.841   (   0.131    0.076   -1.503)    1.510
   7.932   (  -0.432   -0.249    0.858)    0.992
   7.956   (  -0.194   -0.112    0.054)    0.230
   8.106   (   0.060    0.035   -1.039)    1.042
   9.517   (   4.738    2.735    7.491)    9.276
   9.654   (  -5.125   -2.959    4.368)    7.355
======================= Grid point 195 (29/84) =======================
q-point: ( 0.08  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.397   (   3.799    6.580    7.967)   11.009
   0.438   (   5.155    8.930    5.515)   11.693
   0.928   (   5.676    9.831    3.622)   11.916
   0.992   (   7.432   12.872    4.363)   15.490
   1.153   (   6.422   11.123   -2.084)   13.012
   2.282   (  -0.970   -1.680   -5.939)    6.248
   7.958   (  -1.436   -2.487    3.767)    4.737
   8.015   (  -1.907   -3.304    0.055)    3.815
   8.083   (  -0.846   -1.466    2.186)    2.764
   8.284   (   1.628    2.821   -1.678)    3.664
   9.251   (  -0.538   -0.932    5.676)    5.777
   9.913   (  -0.933   -1.616   -4.609)    4.973
======================= Grid point 196 (30/84) =======================
q-point: ( 0.15  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.563   (   6.508    7.415    1.484)    9.977
   0.617   (   4.218   10.150    4.245)   11.782
   1.054   (   2.176    7.362   -1.342)    7.794
   1.317   (  11.542   13.182    6.894)   18.829
   1.390   (  10.051    6.448   -2.451)   12.190
   2.248   (  -0.928   -1.232   -4.937)    5.172
   7.912   (  -2.067   -1.029    1.312)    2.656
   7.966   (   1.055   -4.486    0.201)    4.612
   8.060   (  -0.936   -1.379    1.231)    2.072
   8.295   (  -4.463    4.438   -1.845)    6.559
   9.251   (   1.014    0.352    7.608)    7.684
   9.866   (  -0.406   -4.586   -1.993)    5.017
======================= Grid point 197 (31/84) =======================
q-point: ( 0.23  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.739   (   4.307    7.805    0.193)    8.917
   0.779   (   4.289    8.866    1.978)   10.046
   1.158   (   0.891    7.830   -1.322)    7.991
   1.548   (   5.731    2.365   -4.342)    7.569
   1.655   (  11.920    8.768    5.611)   15.825
   2.240   (   0.926    0.071   -2.762)    2.914
   7.872   (  -1.349   -0.989   -0.591)    1.774
   7.949   (   1.690   -4.150    0.273)    4.489
   8.025   (  -2.426    0.436    0.981)    2.653
   8.238   (  -6.860    4.914   -1.302)    8.538
   9.288   (   2.912   -0.032    8.674)    9.150
   9.815   (   0.657   -6.420   -0.007)    6.454
======================= Grid point 198 (32/84) =======================
q-point: ( 0.31  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.851   (  -0.616    9.138    2.240)    9.429
   0.946   (   6.278    5.812   -2.691)    8.969
   1.232   (  -1.089    7.210   -1.324)    7.411
   1.600   (   0.302   -1.313   -7.083)    7.210
   1.919   (   9.813    6.074    5.349)   12.720
   2.281   (   2.969    0.804   -0.217)    3.083
   7.849   (  -0.669   -0.074   -1.228)    1.400
   7.940   (   1.172   -2.642    0.035)    2.890
   7.988   (  -2.630    1.066    0.913)    2.981
   8.178   (  -5.030    4.041   -0.966)    6.524
   9.346   (   4.410   -0.841    8.784)    9.865
   9.753   (  -0.097   -7.505    1.736)    7.704
======================= Grid point 199 (33/84) =======================
q-point: ( 0.38  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.907   (  -2.914    8.687    2.602)    9.525
   1.094   (   4.865    4.050   -5.010)    8.073
   1.263   (  -2.893    6.314   -1.390)    7.083
   1.560   (  -2.057   -1.989   -8.218)    8.701
   2.100   (   5.916    2.979    4.750)    8.151
   2.336   (   2.520   -0.097    1.322)    2.847
   7.845   (  -0.206    0.581   -1.469)    1.593
   7.937   (   0.859   -1.582   -0.118)    1.804
   7.959   (  -1.953    1.088    0.876)    2.401
   8.147   (  -2.716    2.987   -0.945)    4.147
   9.427   (   6.145   -1.112    7.999)   10.148
   9.660   (  -1.779   -7.778    3.617)    8.761
======================= Grid point 200 (34/84) =======================
q-point: ( 0.46  0.08  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.925   (  -4.552    7.884    2.648)    9.480
   1.165   (  -0.535    0.926   -5.806)    5.903
   1.266   (  -3.340    5.786   -1.470)    6.841
   1.525   (   0.146   -0.253   -8.460)    8.465
   2.168   (   0.287   -0.498    4.460)    4.497
   2.357   (   0.793   -1.374    1.834)    2.425
   7.847   (  -0.297    0.514   -1.566)    1.674
   7.937   (   0.784   -1.358   -0.221)    1.584
   7.947   (  -0.892    1.544    0.930)    2.011
   8.139   (  -1.580    2.737   -0.975)    3.307
   9.512   (   2.468   -4.275    7.359)    8.861
   9.568   (   2.466   -4.271    4.634)    6.767
======================= Grid point 210 (35/84) =======================
q-point: ( 0.15  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.721   (   4.768    8.258   -1.475)    9.649
   0.820   (   5.590    9.682    3.453)   11.701
   1.205   (   4.225    7.317   -1.105)    8.521
   1.499   (   2.827    4.897   -4.509)    7.232
   1.585   (   7.504   12.997    6.467)   16.342
   2.232   (  -0.112   -0.193   -3.268)    3.275
   7.888   (  -1.207   -2.091   -1.540)    2.864
   7.900   (  -0.495   -0.857    1.416)    1.728
   8.054   (  -0.086   -0.149    1.218)    1.230
   8.340   (  -0.026   -0.044   -2.093)    2.093
   9.259   (   0.302    0.524    8.540)    8.561
   9.762   (  -3.125   -5.414    0.168)    6.253
======================= Grid point 211 (36/84) =======================
q-point: ( 0.23  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.882   (   5.085    7.034   -3.848)    9.495
   0.983   (   1.742   10.459    3.084)   11.043
   1.329   (   0.584    8.853   -1.065)    8.936
   1.563   (   2.060   -0.654   -6.426)    6.779
   1.829   (   8.212    8.526    5.608)   13.099
   2.247   (   1.552    0.611   -0.914)    1.902
   7.861   (  -0.731    0.145   -1.438)    1.620
   7.897   (   1.080   -1.390    0.310)    1.787
   8.042   (  -1.710    0.710    1.087)    2.147
   8.303   (  -4.437    1.348   -1.672)    4.929
   9.274   (   1.966   -1.163    8.722)    9.016
   9.665   (  -0.456   -7.903    1.834)    8.126
======================= Grid point 212 (37/84) =======================
q-point: ( 0.31  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.024   (   2.735    7.032   -4.596)    8.835
   1.093   (   0.117    9.380    1.619)    9.519
   1.394   (  -2.975    8.179   -1.908)    8.910
   1.562   (  -0.127   -2.070   -7.222)    7.514
   2.025   (   6.754    4.720    4.992)    9.634
   2.287   (   2.541   -0.283    1.018)    2.752
   7.859   (  -0.419    0.999   -1.360)    1.739
   7.905   (   1.139   -0.972   -0.340)    1.536
   8.012   (  -2.486    0.981    0.957)    2.839
   8.237   (  -4.370    1.544   -1.263)    4.804
   9.311   (   4.481   -2.355    8.256)    9.684
   9.581   (   0.165   -9.021    3.302)    9.608
======================= Grid point 213 (38/84) =======================
q-point: ( 0.38  0.15  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.095   (  -3.360    8.856   -0.909)    9.515
   1.182   (   0.951    5.279   -3.357)    6.328
   1.393   (  -3.855    4.705   -4.063)    7.315
   1.539   (  -1.556    1.231   -5.876)    6.202
   2.141   (   2.981    1.268    4.535)    5.573
   2.316   (   1.915   -1.922    2.010)    3.377
   7.862   (  -0.487    1.078   -1.546)    1.947
   7.913   (   0.681   -0.750   -0.554)    1.155
   7.985   (  -1.693    1.299    0.960)    2.340
   8.195   (  -2.157    1.586   -1.059)    2.879
   9.379   (   6.383   -3.320    7.298)   10.248
   9.492   (  -0.192   -8.435    4.734)    9.674
======================= Grid point 224 (39/84) =======================
q-point: ( 0.23  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.020   (   3.952    6.844   -6.070)    9.965
   1.162   (   4.210    7.291    2.699)    8.841
   1.487   (   3.812    6.603   -0.772)    7.663
   1.529   (  -1.409   -2.441   -7.971)    8.454
   1.986   (   4.166    7.215    5.093)    9.765
   2.261   (   0.394    0.683    0.852)    1.161
   7.875   (   0.491    0.850   -1.886)    2.126
   7.884   (   0.000    0.000    0.203)    0.203
   8.042   (  -0.405   -0.702    0.982)    1.274
   8.290   (  -1.326   -2.297   -1.618)    3.107
   9.256   (  -0.193   -0.335    8.322)    8.331
   9.525   (  -3.495   -6.053    3.351)    7.751
======================= Grid point 225 (40/84) =======================
q-point: ( 0.31  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.147   (   3.167    5.877   -7.673)   10.171
   1.269   (  -0.092    7.459    2.422)    7.843
   1.458   (  -2.959   -2.652   -8.332)    9.231
   1.594   (   0.663    5.917   -1.367)    6.109
   2.107   (   3.392    3.576    4.614)    6.751
   2.269   (   0.954   -1.557    2.164)    2.831
   7.882   (  -0.376    1.169   -1.156)    1.687
   7.898   (   0.834    0.100   -1.042)    1.339
   8.019   (  -1.370   -0.350    0.903)    1.678
   8.236   (  -2.155   -1.366   -1.296)    2.862
   9.270   (   2.867   -1.349    7.531)    8.171
   9.423   (  -1.316   -6.670    4.590)    8.203
======================= Grid point 226 (41/84) =======================
q-point: ( 0.38  0.23  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.214   (  -2.211    3.829   -7.858)    9.017
   1.309   (  -3.499    6.060    2.147)    7.320
   1.403   (   1.031   -1.786   -8.905)    9.141
   1.629   (  -3.157    5.468   -1.392)    6.466
   2.158   (  -0.361    0.626    4.344)    4.404
   2.266   (   1.597   -2.766    2.714)    4.191
   7.885   (  -0.685    1.187   -1.577)    2.089
   7.906   (   0.070   -0.122   -0.923)    0.933
   8.003   (  -0.217    0.376    0.939)    1.035
   8.208   (   0.042   -0.072   -1.115)    1.118
   9.308   (   2.055   -3.560    6.861)    7.998
   9.360   (   2.453   -4.249    5.365)    7.270
======================= Grid point 239 (42/84) =======================
q-point: ( 0.31  0.31  0.14)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.258   (   2.920    5.058   -8.903)   10.647
   1.364   (   0.219    0.379   -2.751)    2.785
   1.386   (  -1.212   -2.099   -4.044)    4.715
   1.690   (   1.594    2.761   -1.064)    3.361
   2.168   (   1.504    2.604    4.136)    5.114
   2.230   (  -1.242   -2.151    3.197)    4.049
   7.901   (   0.368    0.637   -1.071)    1.299
   7.903   (   0.186    0.322   -1.660)    1.701
   8.004   (  -0.431   -0.747    0.979)    1.304
   8.205   (  -0.653   -1.131   -1.119)    1.719
   9.267   (   0.599    1.038    6.679)    6.786
   9.321   (  -2.108   -3.652    5.598)    7.009
======================= Grid point 362 (43/84) =======================
q-point: ( 0.00  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.414   (   0.000    0.000    9.435)    9.435
   0.414   (   0.000   -0.000    9.435)    9.435
   0.804   (   0.000   -0.000   -3.576)    3.576
   0.804   (  -0.000   -0.000   -3.576)    3.576
   1.015   (   0.000    0.000   24.203)   24.203
   2.137   (   0.000   -0.000  -12.212)   12.212
   8.106   (   0.000   -0.000    1.203)    1.203
   8.106   (   0.000    0.000    1.203)    1.203
   8.183   (  -0.000   -0.000   -0.991)    0.991
   8.183   (  -0.000    0.000   -0.991)    0.991
   9.386   (  -0.000   -0.000    5.760)    5.760
   9.767   (   0.000    0.000   -5.169)    5.169
======================= Grid point 363 (44/84) =======================
q-point: ( 0.08  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.465   (   4.176    2.411    8.513)    9.784
   0.567   (   7.735    4.466    8.903)   12.611
   0.832   (   2.400    1.386   -3.347)    4.345
   0.922   (   9.137    5.275   -2.988)   10.966
   1.053   (   5.857    3.381   21.365)   22.410
   2.127   (  -0.787   -0.454  -11.898)   11.933
   8.074   (  -3.090   -1.784    2.081)    4.130
   8.076   (  -1.843   -1.064    1.789)    2.780
   8.172   (  -0.439   -0.253    0.330)    0.605
   8.195   (   0.245    0.141   -1.590)    1.615
   9.383   (   0.043    0.025    6.642)    6.642
   9.787   (   1.433    0.827   -6.381)    6.592
======================= Grid point 364 (45/84) =======================
q-point: ( 0.15  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.585   (   5.871    3.390    7.022)    9.761
   0.668   (   2.869    1.657   13.316)   13.722
   0.908   (   4.092    2.363   -2.987)    5.590
   1.168   (  11.063    6.387   -4.032)   13.396
   1.305   (  14.101    8.141   12.106)   20.290
   2.111   (  -0.287   -0.166  -10.426)   10.432
   7.987   (  -3.935   -2.272    2.787)    5.331
   8.043   (  -1.005   -0.580    2.010)    2.321
   8.149   (  -2.113   -1.220   -0.298)    2.458
   8.178   (  -1.657   -0.956   -1.864)    2.671
   9.403   (   1.903    1.099    8.757)    9.028
   9.821   (   1.400    0.808   -5.596)    5.825
======================= Grid point 365 (46/84) =======================
q-point: ( 0.23  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.716   (   5.330    3.077    6.017)    8.607
   0.748   (   4.243    2.450    9.518)   10.705
   1.008   (   4.395    2.538   -2.842)    5.817
   1.372   (   6.071    3.505   -7.375)   10.175
   1.631   (  13.815    7.976    8.760)   18.199
   2.131   (   2.351    1.358   -6.900)    7.415
   7.906   (  -3.094   -1.786    1.452)    3.856
   8.026   (  -0.745   -0.430    1.527)    1.753
   8.077   (  -3.689   -2.130   -0.820)    4.338
   8.130   (  -2.146   -1.239   -1.467)    2.879
   9.467   (   3.391    1.958   10.438)   11.148
   9.848   (   0.852    0.492   -2.765)    2.935
======================= Grid point 366 (47/84) =======================
q-point: ( 0.31  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.821   (   3.758    2.170    5.467)    6.980
   0.851   (   4.562    2.634    3.839)    6.519
   1.097   (   3.198    1.846   -2.803)    4.636
   1.436   (  -0.075   -0.043  -10.431)   10.431
   1.919   (  11.265    6.504    7.049)   14.795
   2.217   (   4.818    2.782   -2.769)    6.215
   7.849   (  -1.814   -1.047    0.445)    2.141
   8.001   (  -1.349   -0.779    1.320)    2.042
   8.002   (  -2.604   -1.503   -0.959)    3.156
   8.094   (  -0.964   -0.556   -1.357)    1.755
   9.544   (   3.182    1.837   10.966)   11.565
   9.854   (  -0.417   -0.241    0.067)    0.487
======================= Grid point 367 (48/84) =======================
q-point: ( 0.38  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.887   (   2.107    1.216    5.144)    5.690
   0.950   (   4.047    2.337   -1.210)    4.828
   1.150   (   1.438    0.830   -2.656)    3.132
   1.395   (  -2.981   -1.721  -11.664)   12.161
   2.136   (   7.565    4.367    5.852)   10.514
   2.325   (   4.142    2.391   -0.333)    4.794
   7.824   (  -0.422   -0.244    0.145)    0.509
   7.961   (  -1.154   -0.666   -1.218)    1.805
   7.973   (  -0.995   -0.574    1.468)    1.864
   8.081   (  -0.341   -0.197   -1.612)    1.659
   9.613   (   3.015    1.741   10.308)   10.880
   9.825   (  -2.152   -1.242    2.807)    3.749
======================= Grid point 368 (49/84) =======================
q-point: ( 0.46  0.00  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.919   (   0.679    0.392    4.946)    5.008
   1.026   (   2.136    1.233   -4.497)    5.129
   1.168   (   0.291    0.168   -2.447)    2.470
   1.329   (  -2.067   -1.193  -11.768)   12.008
   2.256   (   2.746    1.585    5.078)    5.987
   2.391   (   1.503    0.868    0.736)    1.885
   7.823   (   0.087    0.050    0.019)    0.102
   7.945   (  -0.317   -0.183   -1.398)    1.445
   7.957   (  -0.376   -0.217    1.601)    1.659
   8.077   (  -0.035   -0.020   -1.827)    1.827
   9.685   (   3.298    1.904    8.479)    9.295
   9.761   (  -3.231   -1.866    5.680)    6.796
======================= Grid point 376 (50/84) =======================
q-point: ( 0.08  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.569   (   4.163    7.211    7.228)   11.026
   0.641   (   1.491    2.583   12.826)   13.168
   0.900   (   3.020    5.231   -3.148)    6.811
   1.092   (   5.957   10.318   -4.098)   12.599
   1.213   (   7.350   12.730   14.895)   20.927
   2.113   (  -0.383   -0.663  -11.007)   11.034
   8.024   (  -2.345   -4.061    1.004)    4.796
   8.031   (  -1.504   -2.604    3.479)    4.598
   8.120   (  -1.262   -2.187    1.104)    2.756
   8.240   (   1.122    1.944   -2.692)    3.504
   9.393   (   0.679    1.176    8.095)    8.208
   9.805   (   0.365    0.633   -5.966)    6.011
======================= Grid point 377 (51/84) =======================
q-point: ( 0.15  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.693   (   2.394    4.145    9.691)   10.808
   0.729   (   4.249    7.420    7.471)   11.354
   1.012   (   2.324    7.448   -2.926)    8.333
   1.291   (   7.138    5.336   -6.015)   10.752
   1.503   (  11.224   11.357    9.953)   18.816
   2.114   (   0.862    0.688   -8.471)    8.543
   7.931   (  -3.093   -3.403    0.949)    4.696
   7.991   (   0.443   -4.427    1.787)    4.795
   8.086   (  -1.631   -0.567    1.248)    2.130
   8.247   (  -3.463    3.490   -2.942)    5.729
   9.432   (   2.174    1.572    9.835)   10.194
   9.812   (   1.224   -1.828   -3.442)    4.085
======================= Grid point 378 (52/84) =======================
q-point: ( 0.23  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.782   (   3.714    3.108    5.509)    7.335
   0.865   (   1.717    9.118    5.688)   10.883
   1.117   (   0.860    7.937   -2.804)    8.462
   1.410   (   3.180    0.364   -9.120)    9.666
   1.796   (  11.143    8.402    7.718)   15.948
   2.167   (   3.654    2.308   -4.423)    6.183
   7.862   (  -1.624   -2.155    0.208)    2.706
   7.958   (   0.802   -4.253    0.065)    4.329
   8.051   (  -2.748    1.082    1.619)    3.368
   8.202   (  -5.741    4.165   -2.365)    7.477
   9.490   (   3.412    0.156   10.708)   11.239
   9.811   (   1.988   -4.093   -0.611)    4.591
======================= Grid point 379 (53/84) =======================
q-point: ( 0.31  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.874   (   3.429    2.771    1.215)    4.573
   0.965   (  -0.256    8.607    4.272)    9.612
   1.190   (  -1.070    6.917   -2.967)    7.602
   1.415   (  -1.138   -1.683  -10.900)   11.088
   2.040   (   9.062    5.514    6.330)   12.353
   2.268   (   4.838    2.223   -1.211)    5.460
   7.831   (  -0.662   -0.453   -0.058)    0.805
   7.936   (   0.768   -2.850   -0.778)    3.053
   8.014   (  -2.839    1.344    1.660)    3.553
   8.149   (  -4.586    3.562   -1.979)    6.135
   9.548   (   4.143   -1.375   10.576)   11.441
   9.793   (   1.396   -5.599    1.863)    6.064
======================= Grid point 380 (54/84) =======================
q-point: ( 0.38  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.946   (   1.226    3.294   -0.249)    3.524
   1.033   (  -0.574    6.817    1.255)    6.955
   1.216   (  -2.466    4.574   -3.779)    6.425
   1.364   (  -2.877   -0.395  -10.466)   10.861
   2.205   (   5.377    2.574    5.328)    7.995
   2.356   (   3.332    0.649    0.434)    3.422
   7.825   (  -0.184    0.334   -0.169)    0.417
   7.926   (   0.712   -1.909   -1.190)    2.359
   7.985   (  -2.027    1.349    1.711)    2.976
   8.119   (  -2.649    2.794   -1.924)    4.304
   9.609   (   4.995   -2.073    9.371)   10.819
   9.746   (   0.077   -6.243    4.418)    7.648
======================= Grid point 381 (55/84) =======================
q-point: ( 0.46  0.08  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.973   (  -1.891    3.275    0.437)    3.807
   1.068   (  -2.451    4.245   -1.261)    5.062
   1.210   (  -1.610    2.788   -5.005)    5.951
   1.334   (  -1.561    2.704   -9.133)    9.652
   2.266   (   0.326   -0.565    4.892)    4.935
   2.388   (   0.672   -1.163    0.976)    1.661
   7.826   (  -0.174    0.302   -0.240)    0.423
   7.922   (   0.886   -1.534   -1.359)    2.233
   7.974   (  -0.988    1.711    1.794)    2.669
   8.110   (  -1.558    2.698   -1.946)    3.674
   9.675   (   2.552   -4.420    7.964)    9.459
   9.682   (   2.383   -4.128    6.287)    7.890
======================= Grid point 391 (56/84) =======================
q-point: ( 0.15  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.765   (   2.237    3.875    6.346)    7.765
   0.908   (   5.224    9.048    5.185)   11.663
   1.169   (   4.431    7.674   -2.591)    9.232
   1.369   (   1.434    2.484   -8.541)    9.010
   1.735   (   6.782   11.746    8.110)   15.803
   2.146   (   1.577    2.731   -5.227)    6.104
   7.866   (  -1.494   -2.588   -0.321)    3.006
   7.926   (  -1.154   -1.998    0.768)    2.432
   8.086   (  -0.040   -0.070    2.128)    2.129
   8.285   (   0.094    0.163   -3.314)    3.319
   9.458   (   0.567    0.982   10.432)   10.494
   9.760   (  -1.692   -2.930   -0.468)    3.416
======================= Grid point 392 (57/84) =======================
q-point: ( 0.23  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.845   (   2.930    3.343    1.216)    4.609
   1.062   (   1.592    9.849    4.626)   10.997
   1.290   (  -0.447    8.633   -3.185)    9.213
   1.394   (   1.624   -1.618   -9.855)   10.118
   1.956   (   7.530    7.743    6.672)   12.695
   2.218   (   3.371    2.747   -2.007)    4.789
   7.837   (  -0.406   -0.324   -0.536)    0.747
   7.902   (   0.247   -1.625   -0.207)    1.657
   8.073   (  -1.929    0.734    2.119)    2.958
   8.257   (  -3.699    1.163   -2.882)    4.831
   9.474   (   1.893   -1.523   10.372)   10.653
   9.706   (   0.469   -5.793    1.898)    6.114
======================= Grid point 393 (58/84) =======================
q-point: ( 0.31  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.928   (   3.075    2.834   -3.240)    5.290
   1.157   (  -1.723    9.732    3.832)   10.600
   1.300   (  -3.145    0.998   -8.418)    9.042
   1.420   (  -0.811    4.989   -5.197)    7.249
   2.135   (   6.192    4.122    5.552)    9.282
   2.301   (   3.619    0.977    0.103)    3.750
   7.835   (  -0.290    0.792   -0.565)    1.015
   7.897   (   0.658   -1.264   -0.753)    1.612
   8.040   (  -2.583    1.002    1.929)    3.376
   8.201   (  -3.836    1.334   -2.386)    4.710
   9.499   (   3.858   -3.200    9.668)   10.890
   9.658   (   1.290   -7.309    3.838)    8.356
======================= Grid point 394 (59/84) =======================
q-point: ( 0.38  0.15  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.997   (   1.448    2.144   -5.929)    6.469
   1.203   (  -2.932    8.281    1.782)    8.964
   1.244   (  -2.743   -0.623   -8.521)    8.973
   1.445   (  -3.764    7.698   -3.619)    9.302
   2.239   (   2.744    0.886    4.782)    5.584
   2.350   (   2.143   -1.175    1.168)    2.709
   7.837   (  -0.481    0.872   -0.634)    1.181
   7.897   (   0.571   -1.016   -1.163)    1.646
   8.014   (  -1.645    1.356    1.945)    2.885
   8.164   (  -1.961    1.490   -2.120)    3.250
   9.544   (   5.424   -4.093    8.256)   10.693
   9.604   (   0.966   -7.423    5.844)    9.497
======================= Grid point 405 (60/84) =======================
q-point: ( 0.23  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.914   (   2.116    3.665   -3.261)    5.343
   1.226   (   3.492    6.048    2.660)    7.473
   1.333   (  -1.449   -2.509  -10.189)   10.593
   1.459   (   3.821    6.617   -2.107)    7.926
   2.098   (   3.721    6.445    5.637)    9.336
   2.271   (   1.372    2.376   -0.076)    2.744
   7.844   (   0.438    0.759   -0.884)    1.245
   7.883   (  -0.260   -0.451   -0.377)    0.643
   8.071   (  -0.499   -0.865    2.023)    2.256
   8.246   (  -1.143   -1.980   -2.747)    3.574
   9.445   (  -0.621   -1.076    9.672)    9.751
   9.602   (  -2.596   -4.496    3.811)    6.440
======================= Grid point 406 (61/84) =======================
q-point: ( 0.31  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.989   (   2.324    3.366   -6.971)    8.083
   1.236   (  -2.051   -3.241  -11.144)   11.785
   1.346   (  -0.694    6.558    3.299)    7.374
   1.560   (  -0.070    6.845   -2.280)    7.215
   2.204   (   3.029    2.902    4.697)    6.298
   2.306   (   1.516   -0.354    1.323)    2.043
   7.855   (  -0.304    1.076   -0.985)    1.490
   7.881   (   0.421   -0.396   -0.815)    0.999
   8.046   (  -1.320   -0.385    1.899)    2.345
   8.200   (  -1.840   -1.166   -2.310)    3.176
   9.441   (   2.131   -2.284    8.655)    9.201
   9.529   (  -0.427   -5.427    5.332)    7.620
======================= Grid point 407 (62/84) =======================
q-point: ( 0.38  0.23  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.037   (  -1.090    1.888   -8.644)    8.915
   1.189   (   1.453   -2.517  -11.698)   12.054
   1.374   (  -3.711    6.428    4.037)    8.449
   1.592   (  -3.532    6.118   -2.379)    7.454
   2.247   (  -0.114    0.198    4.201)    4.207
   2.316   (   1.152   -1.996    1.948)    3.018
   7.857   (  -0.614    1.063   -1.003)    1.585
   7.884   (   0.205   -0.356   -1.154)    1.225
   8.031   (  -0.208    0.360    1.936)    1.980
   8.176   (   0.024   -0.041   -2.086)    2.086
   9.462   (   2.250   -3.898    7.696)    8.915
   9.485   (   2.232   -3.866    6.430)    7.827
======================= Grid point 420 (63/84) =======================
q-point: ( 0.31  0.31  0.29)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.061   (   2.129    3.687   -9.773)   10.660
   1.158   (  -2.306   -3.994  -11.611)   12.493
   1.441   (   1.464    2.536    4.021)    4.974
   1.657   (   1.528    2.646   -2.265)    3.803
   2.251   (   1.071    1.856    3.771)    4.337
   2.290   (  -0.667   -1.155    2.491)    2.825
   7.872   (   0.268    0.464   -1.256)    1.366
   7.878   (   0.023    0.039   -1.076)    1.077
   8.033   (  -0.350   -0.607    1.915)    2.039
   8.174   (  -0.552   -0.956   -2.028)    2.309
   9.418   (   0.039    0.068    7.644)    7.644
   9.449   (  -1.523   -2.638    6.456)    7.138
======================= Grid point 543 (64/84) =======================
q-point: ( 0.00  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.583   (   0.000   -0.000    7.653)    7.653
   0.583   (  -0.000    0.000    7.653)    7.653
   0.714   (   0.000   -0.000   -5.608)    5.608
   0.714   (   0.000    0.000   -5.608)    5.608
   1.465   (   0.000    0.000   21.272)   21.272
   1.845   (  -0.000   -0.000  -17.246)   17.246
   8.133   (   0.000    0.000    1.415)    1.415
   8.133   (  -0.000    0.000    1.415)    1.415
   8.160   (   0.000    0.000   -1.321)    1.321
   8.160   (   0.000    0.000   -1.321)    1.321
   9.513   (   0.000   -0.000    6.925)    6.925
   9.649   (  -0.000   -0.000   -6.663)    6.663
======================= Grid point 544 (65/84) =======================
q-point: ( 0.08  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.622   (   3.267    1.886    7.252)    8.175
   0.718   (   9.276    5.356    6.797)   12.686
   0.747   (   2.770    1.599   -5.343)    6.227
   0.837   (   9.280    5.358   -5.271)   11.942
   1.476   (   1.700    0.981   20.483)   20.576
   1.844   (   0.108    0.062  -16.700)   16.700
   8.116   (  -1.602   -0.925    2.219)    2.889
   8.116   (  -1.302   -0.752    2.171)    2.641
   8.157   (  -0.615   -0.355   -1.688)    1.832
   8.158   (  -0.436   -0.252   -2.064)    2.124
   9.522   (   1.010    0.583    7.155)    7.250
   9.659   (   1.087    0.627   -6.737)    6.853
======================= Grid point 545 (66/84) =======================
q-point: ( 0.15  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.720   (   4.972    2.871    6.423)    8.615
   0.831   (   4.416    2.550   -4.815)    7.013
   0.898   (   5.411    3.124    9.779)   11.604
   1.059   (   8.735    5.043   -6.772)   12.148
   1.582   (   7.888    4.554   15.053)   17.594
   1.870   (   2.716    1.568  -13.774)   14.127
   8.050   (  -3.866   -2.232    3.549)    5.703
   8.087   (  -1.240   -0.716    2.450)    2.838
   8.117   (  -2.960   -1.709   -2.861)    4.457
   8.136   (  -1.472   -0.850   -2.385)    2.929
   9.569   (   3.156    1.822    7.781)    8.592
   9.705   (   2.906    1.678   -6.211)    7.060
======================= Grid point 546 (67/84) =======================
q-point: ( 0.23  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.835   (   4.802    2.773    5.768)    8.001
   0.936   (   4.505    2.601   -4.456)    6.850
   0.971   (   1.657    0.956   12.037)   12.188
   1.195   (   3.063    1.768  -10.378)   10.964
   1.808   (  10.966    6.331    9.017)   15.545
   1.980   (   6.675    3.854   -8.313)   11.336
   7.957   (  -3.858   -2.227    3.644)    5.755
   8.033   (  -3.897   -2.250   -3.528)    5.718
   8.060   (  -1.226   -0.708    1.882)    2.355
   8.097   (  -1.726   -0.997   -1.855)    2.723
   9.654   (   3.876    2.238    8.056)    9.216
   9.777   (   3.047    1.759   -4.642)    5.825
======================= Grid point 547 (68/84) =======================
q-point: ( 0.31  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.930   (   3.465    2.000    5.420)    6.736
   0.992   (   0.366    0.211    9.879)    9.888
   1.027   (   3.289    1.899   -4.297)    5.735
   1.212   (  -1.105   -0.638  -11.835)   11.904
   2.049   (   9.743    5.625    5.943)   12.724
   2.149   (   7.559    4.364   -4.213)    9.692
   7.886   (  -2.256   -1.303    3.160)    4.096
   7.957   (  -2.610   -1.507   -3.455)    4.585
   8.031   (  -1.260   -0.727    1.663)    2.210
   8.064   (  -1.122   -0.648   -1.679)    2.121
   9.732   (   2.844    1.642    7.595)    8.274
   9.833   (   1.668    0.963   -2.738)    3.348
======================= Grid point 548 (69/84) =======================
q-point: ( 0.38  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.991   (   1.875    1.083    5.154)    5.590
   0.995   (   0.032    0.018    5.837)    5.837
   1.083   (   1.645    0.950   -4.095)    4.514
   1.166   (  -2.563   -1.480  -10.855)   11.251
   2.238   (   6.583    3.801    4.295)    8.731
   2.302   (   5.446    3.144   -2.153)    6.647
   7.854   (  -0.695   -0.401    2.679)    2.797
   7.916   (  -1.104   -0.637   -3.137)    3.386
   8.006   (  -0.820   -0.473    1.892)    2.115
   8.045   (  -0.586   -0.338   -1.956)    2.070
   9.784   (   1.763    1.018    6.379)    6.696
   9.852   (   0.098    0.057   -0.651)    0.661
======================= Grid point 549 (70/84) =======================
q-point: ( 0.46  0.00  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.998   (   0.187    0.108    2.211)    2.221
   1.018   (   0.564    0.325    4.922)    4.965
   1.106   (   0.427    0.246   -3.854)    3.885
   1.114   (  -1.494   -0.863   -9.141)    9.302
   2.342   (   2.341    1.352    3.424)    4.362
   2.387   (   1.937    1.118   -1.227)    2.551
   7.847   (  -0.052   -0.030    2.295)    2.295
   7.901   (  -0.266   -0.154   -2.825)    2.842
   7.994   (  -0.282   -0.163    2.092)    2.117
   8.037   (  -0.160   -0.093   -2.192)    2.200
   9.817   (   1.225    0.707    4.350)    4.574
   9.842   (  -0.821   -0.474    1.800)    2.034
======================= Grid point 557 (71/84) =======================
q-point: ( 0.08  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.707   (   3.514    6.086    6.521)    9.587
   0.816   (   3.043    5.270   -3.414)    6.977
   0.854   (   3.727    6.455    7.063)   10.269
   0.992   (   5.137    8.898   -6.027)   11.911
   1.537   (   3.415    5.916   17.055)   18.372
   1.856   (   0.938    1.625  -14.927)   15.045
   8.056   (  -2.504   -4.337    2.042)    5.408
   8.082   (  -2.115   -3.664    0.611)    4.274
   8.131   (  -0.240   -0.416    1.237)    1.327
   8.182   (   0.583    1.010   -3.074)    3.288
   9.549   (   1.336    2.314    7.477)    7.940
   9.684   (   1.120    1.941   -6.330)    6.715
======================= Grid point 558 (72/84) =======================
q-point: ( 0.15  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.835   (   2.564    8.132    5.896)   10.367
   0.927   (   3.156    4.415    2.933)    6.169
   0.952   (   2.324    4.471    4.054)    6.467
   1.143   (   4.802    2.967   -8.922)   10.557
   1.712   (   8.033    8.243   10.982)   15.909
   1.925   (   4.070    4.104  -10.456)   11.947
   7.963   (  -2.402   -5.017    2.239)    5.996
   8.002   (  -1.394   -5.441   -1.090)    5.721
   8.120   (  -2.063    1.597    2.611)    3.691
   8.183   (  -2.685    2.611   -3.457)    5.097
   9.609   (   2.912    1.866    7.712)    8.452
   9.731   (   2.628    0.674   -4.893)    5.595
======================= Grid point 559 (73/84) =======================
q-point: ( 0.23  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.929   (   0.821    3.926    7.561)    8.559
   0.997   (   1.519    4.231    7.813)    9.014
   1.048   (   0.767    7.724   -3.688)    8.594
   1.207   (   1.417   -0.324  -11.080)   11.175
   1.943   (   9.404    6.989    7.019)   13.658
   2.068   (   6.793    4.859   -5.693)   10.108
   7.889   (  -1.236   -3.328    2.322)    4.242
   7.938   (  -0.403   -4.254   -2.111)    4.766
   8.091   (  -3.424    2.261    2.531)    4.821
   8.147   (  -4.575    3.302   -2.999)    6.390
   9.675   (   3.684   -0.313    7.512)    8.373
   9.778   (   3.271   -1.760   -3.099)    4.838
======================= Grid point 560 (74/84) =======================
q-point: ( 0.31  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.954   (   0.862   -0.515    6.850)    6.923
   1.061   (  -1.008    6.715    3.949)    7.856
   1.110   (  -1.043    4.640   -4.529)    6.567
   1.206   (  -1.623    2.008   -9.137)    9.495
   2.152   (   7.972    4.629    4.902)   10.441
   2.229   (   6.563    3.525   -2.893)    7.992
   7.854   (  -0.372   -1.299    2.204)    2.585
   7.903   (   0.128   -2.403   -2.425)    3.416
   8.053   (  -3.268    2.022    2.312)    4.485
   8.103   (  -3.936    2.862   -2.584)    5.510
   9.725   (   3.702   -2.310    6.758)    8.044
   9.804   (   2.849   -3.713   -1.338)    4.867
======================= Grid point 561 (75/84) =======================
q-point: ( 0.38  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.960   (   1.069   -1.281    3.495)    3.873
   1.080   (  -0.900   -0.067   -5.448)    5.522
   1.117   (  -3.289    6.356    1.912)    7.408
   1.207   (  -3.501    6.677   -5.416)    9.283
   2.296   (   4.746    1.975    3.668)    6.315
   2.347   (   4.076    1.333   -1.485)    4.538
   7.843   (   0.002   -0.408    1.926)    1.969
   7.889   (   0.308   -1.426   -2.374)    2.786
   8.025   (  -2.156    1.807    2.346)    3.663
   8.074   (  -2.411    2.368   -2.506)    4.207
   9.759   (   3.609   -3.307    5.262)    7.186
   9.804   (   2.051   -4.676    0.713)    5.156
======================= Grid point 562 (76/84) =======================
q-point: ( 0.46  0.08  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.965   (   0.853   -1.477    1.561)    2.312
   1.063   (   0.993   -1.720   -7.852)    8.099
   1.123   (  -4.470    7.743    4.378)    9.955
   1.210   (  -4.241    7.345   -4.573)    9.636
   2.347   (   0.433   -0.750    3.201)    3.316
   2.389   (   0.553   -0.958   -1.007)    1.496
   7.842   (   0.171   -0.296    1.755)    1.788
   7.885   (   0.630   -1.091   -2.299)    2.621
   8.016   (  -1.178    2.040    2.414)    3.373
   8.065   (  -1.406    2.436   -2.528)    3.782
   9.783   (   2.445   -4.235    3.210)    5.849
   9.788   (   2.546   -4.410    2.951)    5.885
======================= Grid point 572 (77/84) =======================
q-point: ( 0.15  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.945   (   0.138    0.239   10.900)   10.903
   1.011   (   4.782    8.282    4.890)   10.741
   1.105   (   4.636    8.030   -3.763)   10.007
   1.173   (   0.248    0.430  -10.890)   10.901
   1.891   (   5.528    9.574    7.616)   13.425
   2.031   (   3.678    6.370   -6.495)    9.813
   7.881   (  -1.701   -2.947    1.704)    3.806
   7.918   (  -1.693   -2.933   -1.551)    3.725
   8.141   (   0.124    0.215    3.377)    3.386
   8.214   (   0.151    0.261   -3.842)    3.854
   9.635   (   0.472    0.817    7.123)    7.185
   9.732   (  -0.275   -0.476   -2.831)    2.884
======================= Grid point 573 (78/84) =======================
q-point: ( 0.23  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.942   (  -0.024   -0.632    8.181)    8.206
   1.123   (  -0.659    1.918   -5.570)    5.927
   1.180   (   0.524    6.481   -0.429)    6.516
   1.253   (   2.332    7.525   -4.395)    9.021
   2.077   (   6.505    6.464    5.379)   10.631
   2.165   (   5.108    4.781   -3.553)    7.846
   7.847   (  -0.266   -0.909    1.465)    1.744
   7.882   (  -0.161   -1.528   -1.722)    2.308
   8.126   (  -2.369    0.893    3.187)    4.070
   8.193   (  -3.005    1.036   -3.507)    4.733
   9.646   (   1.755   -2.333    6.438)    7.069
   9.719   (   1.345   -3.744   -1.159)    4.143
======================= Grid point 574 (79/84) =======================
q-point: ( 0.31  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.936   (   0.426   -1.040    4.358)    4.501
   1.085   (  -1.686   -2.610   -9.780)   10.262
   1.257   (  -2.068    9.488    4.180)   10.572
   1.339   (  -1.766    8.939   -3.940)    9.927
   2.230   (   5.461    3.187    3.891)    7.424
   2.285   (   4.609    2.133   -1.761)    5.375
   7.843   (  -0.061    0.072    1.188)    1.192
   7.873   (   0.199   -0.847   -1.603)    1.823
   8.088   (  -2.857    1.116    2.798)    4.152
   8.146   (  -3.313    1.222   -3.020)    4.646
   9.654   (   3.200   -4.327    5.504)    7.698
   9.705   (   2.416   -5.678    0.290)    6.178
======================= Grid point 575 (80/84) =======================
q-point: ( 0.38  0.15  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.938   (   0.855   -0.950    0.876)    1.550
   1.036   (  -0.374   -2.666   -8.531)    8.945
   1.294   (  -4.505    9.509    4.335)   11.380
   1.374   (  -4.358    9.156   -3.787)   10.824
   2.317   (   2.543    0.241    3.004)    3.943
   2.354   (   2.355   -0.456   -0.815)    2.534
   7.841   (  -0.212    0.212    0.985)    1.029
   7.869   (   0.186   -0.625   -1.568)    1.699
   8.060   (  -1.673    1.439    2.662)    3.457
   8.115   (  -1.788    1.466   -2.764)    3.604
   9.670   (   4.043   -5.165    3.980)    7.672
   9.688   (   2.587   -6.318    2.056)    7.130
======================= Grid point 586 (81/84) =======================
q-point: ( 0.23  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.927   (  -0.407   -0.705    4.781)    4.850
   1.086   (  -1.820   -3.152  -10.541)   11.152
   1.329   (   3.845    6.660    4.093)    8.711
   1.405   (   3.850    6.669   -3.382)    8.411
   2.194   (   3.052    5.286    3.986)    7.290
   2.252   (   2.266    3.925   -1.903)    4.915
   7.844   (   0.184    0.318    0.806)    0.886
   7.867   (  -0.170   -0.295   -1.217)    1.263
   8.122   (  -0.663   -1.149    3.031)    3.308
   8.186   (  -0.905   -1.567   -3.300)    3.763
   9.599   (  -1.137   -1.969    5.426)    5.883
   9.649   (  -1.761   -3.050    0.314)    3.536
======================= Grid point 587 (82/84) =======================
q-point: ( 0.31  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 592
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.920   (   0.352   -0.430    0.794)    0.969
   1.023   (  -1.501   -3.285   -9.034)    9.730
   1.427   (  -0.431    7.116    4.211)    8.280
   1.503   (  -0.328    7.033   -3.429)    7.831
   2.279   (   2.562    1.896    2.841)    4.269
   2.314   (   2.061    0.802   -0.659)    2.307
   7.850   (  -0.144    0.502    0.350)    0.629
   7.864   (   0.016   -0.188   -0.920)    0.939
   8.092   (  -1.367   -0.537    2.685)    3.061
   8.148   (  -1.562   -0.848   -2.842)    3.352
   9.575   (   1.342   -3.283    4.384)    5.639
   9.602   (   0.524   -4.288    1.546)    4.588
======================= Grid point 588 (83/84) =======================
q-point: ( 0.38  0.23  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.925   (   0.177   -0.307   -1.633)    1.671
   0.988   (   1.186   -2.055   -7.588)    7.951
   1.458   (  -3.694    6.399    4.229)    8.513
   1.535   (  -3.647    6.316   -3.474)    8.079
   2.312   (   0.280   -0.485    2.277)    2.345
   2.335   (   0.709   -1.229   -0.076)    1.421
   7.850   (  -0.329    0.570    0.246)    0.703
   7.863   (   0.034   -0.059   -0.941)    0.944
   8.076   (  -0.153    0.265    2.577)    2.595
   8.129   (  -0.057    0.098   -2.638)    2.640
   9.576   (   2.298   -3.981    3.290)    5.653
   9.581   (   2.223   -3.850    2.695)    5.199
======================= Grid point 601 (84/84) =======================
q-point: ( 0.31  0.31  0.43)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 0.00e+00 0.00e+00 
Number of triplets: 322
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.924   (   0.485    0.841   -2.969)    3.124
   0.965   (  -1.362   -2.358   -6.817)    7.341
   1.524   (   1.475    2.555    4.206)    5.138
   1.601   (   1.491    2.582   -3.434)    4.548
   2.307   (   0.514    0.890    1.769)    2.046
   2.320   (  -0.069   -0.120    0.445)    0.466
   7.858   (   0.153    0.265   -0.197)    0.364
   7.862   (  -0.020   -0.034   -0.553)    0.555
   8.077   (  -0.355   -0.614    2.497)    2.596
   8.128   (  -0.438   -0.759   -2.537)    2.684
   9.534   (  -0.495   -0.857    3.462)    3.601
   9.543   (  -0.982   -1.700    2.487)    3.169
=================== 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    147.521    147.521    159.449      0.000     -0.000      0.000 3/14196
   20.0     72.054     72.054     77.962      0.000     -0.000      0.000 3/14196
   30.0     46.278     46.278     50.818      0.000     -0.000      0.000 3/14196
   40.0     34.235     34.235     37.884      0.000     -0.000      0.000 3/14196
   50.0     27.065     27.065     30.069      0.000     -0.000      0.000 3/14196
   60.0     22.156     22.156     24.662      0.000     -0.000      0.000 3/14196
   70.0     18.563     18.563     20.668      0.000     -0.000      0.000 3/14196
   80.0     15.855     15.855     17.635      0.000     -0.000      0.000 3/14196
   90.0     13.777     13.777     15.295      0.000     -0.000      0.000 3/14196
  100.0     12.151     12.151     13.462      0.000     -0.000      0.000 3/14196
  110.0     10.855     10.855     12.002      0.000     -0.000      0.000 3/14196
  120.0      9.804      9.804     10.819      0.000     -0.000      0.000 3/14196
  130.0      8.937      8.937      9.845      0.000     -0.000      0.000 3/14196
  140.0      8.211      8.211      9.032      0.000     -0.000      0.000 3/14196
  150.0      7.595      7.595      8.343      0.000     -0.000      0.000 3/14196
  160.0      7.066      7.066      7.754      0.000     -0.000      0.000 3/14196
  170.0      6.606      6.606      7.243      0.000     -0.000      0.000 3/14196
  180.0      6.204      6.204      6.796      0.000     -0.000      0.000 3/14196
  190.0      5.849      5.849      6.402      0.000     -0.000      0.000 3/14196
  200.0      5.532      5.532      6.053      0.000     -0.000      0.000 3/14196
  210.0      5.249      5.249      5.740      0.000     -0.000      0.000 3/14196
  220.0      4.994      4.994      5.459      0.000     -0.000      0.000 3/14196
  230.0      4.763      4.763      5.204      0.000     -0.000      0.000 3/14196
  240.0      4.553      4.553      4.973      0.000     -0.000      0.000 3/14196
  250.0      4.361      4.361      4.762      0.000     -0.000      0.000 3/14196
  260.0      4.185      4.185      4.568      0.000     -0.000      0.000 3/14196
  270.0      4.022      4.022      4.390      0.000     -0.000      0.000 3/14196
  280.0      3.872      3.872      4.225      0.000     -0.000      0.000 3/14196
  290.0      3.733      3.733      4.072      0.000     -0.000      0.000 3/14196
  300.0      3.604      3.604      3.931      0.000     -0.000      0.000 3/14196
  310.0      3.483      3.483      3.799      0.000     -0.000      0.000 3/14196
  320.0      3.371      3.371      3.675      0.000     -0.000      0.000 3/14196
  330.0      3.265      3.265      3.560      0.000     -0.000      0.000 3/14196
  340.0      3.166      3.166      3.452      0.000     -0.000      0.000 3/14196
  350.0      3.073      3.073      3.350      0.000     -0.000      0.000 3/14196
  360.0      2.985      2.985      3.254      0.000     -0.000      0.000 3/14196
  370.0      2.903      2.903      3.164      0.000     -0.000      0.000 3/14196
  380.0      2.824      2.824      3.078      0.000     -0.000      0.000 3/14196
  390.0      2.750      2.750      2.997      0.000     -0.000      0.000 3/14196
  400.0      2.680      2.680      2.920      0.000     -0.000      0.000 3/14196
  410.0      2.613      2.613      2.847      0.000     -0.000      0.000 3/14196
  420.0      2.550      2.550      2.778      0.000     -0.000      0.000 3/14196
  430.0      2.489      2.489      2.712      0.000     -0.000      0.000 3/14196
  440.0      2.432      2.432      2.649      0.000     -0.000      0.000 3/14196
  450.0      2.377      2.377      2.589      0.000     -0.000      0.000 3/14196
  460.0      2.324      2.324      2.532      0.000     -0.000      0.000 3/14196
  470.0      2.274      2.274      2.477      0.000     -0.000      0.000 3/14196
  480.0      2.226      2.226      2.425      0.000     -0.000      0.000 3/14196
  490.0      2.180      2.180      2.374      0.000     -0.000      0.000 3/14196
  500.0      2.136      2.136      2.326      0.000     -0.000      0.000 3/14196
  510.0      2.093      2.093      2.280      0.000     -0.000      0.000 3/14196
  520.0      2.052      2.052      2.235      0.000     -0.000      0.000 3/14196
  530.0      2.013      2.013      2.193      0.000     -0.000      0.000 3/14196
  540.0      1.976      1.976      2.151      0.000     -0.000      0.000 3/14196
  550.0      1.939      1.939      2.112      0.000     -0.000      0.000 3/14196
  560.0      1.904      1.904      2.074      0.000     -0.000      0.000 3/14196
  570.0      1.870      1.870      2.037      0.000     -0.000      0.000 3/14196
  580.0      1.838      1.838      2.001      0.000     -0.000      0.000 3/14196
  590.0      1.806      1.806      1.967      0.000     -0.000      0.000 3/14196
  600.0      1.776      1.776      1.934      0.000     -0.000      0.000 3/14196
  610.0      1.747      1.747      1.902      0.000     -0.000      0.000 3/14196
  620.0      1.718      1.718      1.871      0.000     -0.000      0.000 3/14196
  630.0      1.691      1.691      1.841      0.000     -0.000      0.000 3/14196
  640.0      1.664      1.664      1.812      0.000     -0.000      0.000 3/14196
  650.0      1.638      1.638      1.784      0.000     -0.000      0.000 3/14196
  660.0      1.613      1.613      1.757      0.000     -0.000      0.000 3/14196
  670.0      1.589      1.589      1.730      0.000     -0.000      0.000 3/14196
  680.0      1.565      1.565      1.704      0.000     -0.000      0.000 3/14196
  690.0      1.543      1.543      1.680      0.000     -0.000      0.000 3/14196
  700.0      1.520      1.520      1.655      0.000     -0.000      0.000 3/14196
  710.0      1.499      1.499      1.632      0.000     -0.000      0.000 3/14196
  720.0      1.478      1.478      1.609      0.000     -0.000      0.000 3/14196
  730.0      1.458      1.458      1.587      0.000     -0.000      0.000 3/14196
  740.0      1.438      1.438      1.565      0.000     -0.000      0.000 3/14196
  750.0      1.418      1.418      1.544      0.000     -0.000      0.000 3/14196
  760.0      1.400      1.400      1.524      0.000     -0.000      0.000 3/14196
  770.0      1.381      1.381      1.504      0.000     -0.000      0.000 3/14196
  780.0      1.364      1.364      1.485      0.000     -0.000      0.000 3/14196
  790.0      1.346      1.346      1.466      0.000     -0.000      0.000 3/14196
  800.0      1.329      1.329      1.447      0.000     -0.000      0.000 3/14196
  810.0      1.313      1.313      1.429      0.000     -0.000      0.000 3/14196
  820.0      1.297      1.297      1.412      0.000     -0.000      0.000 3/14196
  830.0      1.281      1.281      1.395      0.000     -0.000      0.000 3/14196
  840.0      1.266      1.266      1.378      0.000     -0.000      0.000 3/14196
  850.0      1.251      1.251      1.362      0.000     -0.000      0.000 3/14196
  860.0      1.236      1.236      1.346      0.000     -0.000      0.000 3/14196
  870.0      1.222      1.222      1.330      0.000     -0.000      0.000 3/14196
  880.0      1.208      1.208      1.315      0.000     -0.000      0.000 3/14196
  890.0      1.194      1.194      1.300      0.000     -0.000      0.000 3/14196
  900.0      1.181      1.181      1.286      0.000     -0.000      0.000 3/14196
  910.0      1.168      1.168      1.272      0.000     -0.000      0.000 3/14196
  920.0      1.155      1.155      1.258      0.000     -0.000      0.000 3/14196
  930.0      1.143      1.143      1.244      0.000     -0.000      0.000 3/14196
  940.0      1.131      1.131      1.231      0.000     -0.000      0.000 3/14196
  950.0      1.119      1.119      1.218      0.000     -0.000      0.000 3/14196
  960.0      1.107      1.107      1.205      0.000     -0.000      0.000 3/14196
  970.0      1.095      1.095      1.193      0.000     -0.000      0.000 3/14196
  980.0      1.084      1.084      1.181      0.000     -0.000      0.000 3/14196
  990.0      1.073      1.073      1.169      0.000     -0.000      0.000 3/14196
 1000.0      1.063      1.063      1.157      0.000     -0.000      0.000 3/14196

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

