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

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

-------------------------[time 2026-01-07 23:04: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.001504419534202    0.000000000000000    0.000000000000000
  b   -2.000752209767101    3.465404480672323    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.574457220000000
Atomic positions (fractional):
   *1 Al  0.33333333333333  0.66666666666667  0.00072316479453  26.982
    2 Al  0.66666666666667  0.33333333333333  0.50072316479453  26.982
   *3 As  0.33333333333333  0.66666666666667  0.37527683520547  74.922
    4 As  0.66666666666667  0.33333333333333  0.87527683520547  74.922
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.001504419534202    0.000000000000000    0.000000000000000
  b   -2.000752209767101    3.465404480672323    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.574457220000000
Atomic positions (fractional):
   *1 Al  0.33333333333333  0.66666666666667  0.00072316479453  26.982 > 1
    2 Al  0.66666666666667  0.33333333333333  0.50072316479453  26.982 > 2
   *3 As  0.33333333333333  0.66666666666667  0.37527683520547  74.922 > 3
    4 As  0.66666666666667  0.33333333333333  0.87527683520547  74.922 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   16.006017678136807    0.000000000000000    0.000000000000000
  b   -8.003008839068404   13.861617922689293    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.148914440000000
Atomic positions (fractional):
   *1 Al  0.08333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    2 Al  0.33333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    3 Al  0.58333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    4 Al  0.83333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    5 Al  0.08333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    6 Al  0.33333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    7 Al  0.58333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    8 Al  0.83333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    9 Al  0.08333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   10 Al  0.33333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   11 Al  0.58333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   12 Al  0.83333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   13 Al  0.08333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   14 Al  0.33333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   15 Al  0.58333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   16 Al  0.83333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   17 Al  0.08333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   18 Al  0.33333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   19 Al  0.58333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   20 Al  0.83333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   21 Al  0.08333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   22 Al  0.33333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   23 Al  0.58333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   24 Al  0.83333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   25 Al  0.08333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   26 Al  0.33333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   27 Al  0.58333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   28 Al  0.83333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   29 Al  0.08333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   30 Al  0.33333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   31 Al  0.58333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   32 Al  0.83333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   33 Al  0.16666666666667  0.08333333333333  0.25036158239727  26.982 > 2
   34 Al  0.41666666666667  0.08333333333333  0.25036158239727  26.982 > 2
   35 Al  0.66666666666667  0.08333333333333  0.25036158239727  26.982 > 2
   36 Al  0.91666666666667  0.08333333333333  0.25036158239727  26.982 > 2
   37 Al  0.16666666666667  0.33333333333333  0.25036158239727  26.982 > 2
   38 Al  0.41666666666667  0.33333333333333  0.25036158239727  26.982 > 2
   39 Al  0.66666666666667  0.33333333333333  0.25036158239727  26.982 > 2
   40 Al  0.91666666666667  0.33333333333333  0.25036158239727  26.982 > 2
   41 Al  0.16666666666667  0.58333333333333  0.25036158239727  26.982 > 2
   42 Al  0.41666666666667  0.58333333333333  0.25036158239727  26.982 > 2
   43 Al  0.66666666666667  0.58333333333333  0.25036158239727  26.982 > 2
   44 Al  0.91666666666667  0.58333333333333  0.25036158239727  26.982 > 2
   45 Al  0.16666666666667  0.83333333333333  0.25036158239727  26.982 > 2
   46 Al  0.41666666666667  0.83333333333333  0.25036158239727  26.982 > 2
   47 Al  0.66666666666667  0.83333333333333  0.25036158239727  26.982 > 2
   48 Al  0.91666666666667  0.83333333333333  0.25036158239727  26.982 > 2
   49 Al  0.16666666666667  0.08333333333333  0.75036158239727  26.982 > 2
   50 Al  0.41666666666667  0.08333333333333  0.75036158239727  26.982 > 2
   51 Al  0.66666666666667  0.08333333333333  0.75036158239727  26.982 > 2
   52 Al  0.91666666666667  0.08333333333333  0.75036158239727  26.982 > 2
   53 Al  0.16666666666667  0.33333333333333  0.75036158239727  26.982 > 2
   54 Al  0.41666666666667  0.33333333333333  0.75036158239727  26.982 > 2
   55 Al  0.66666666666667  0.33333333333333  0.75036158239727  26.982 > 2
   56 Al  0.91666666666667  0.33333333333333  0.75036158239727  26.982 > 2
   57 Al  0.16666666666667  0.58333333333333  0.75036158239727  26.982 > 2
   58 Al  0.41666666666667  0.58333333333333  0.75036158239727  26.982 > 2
   59 Al  0.66666666666667  0.58333333333333  0.75036158239727  26.982 > 2
   60 Al  0.91666666666667  0.58333333333333  0.75036158239727  26.982 > 2
   61 Al  0.16666666666667  0.83333333333333  0.75036158239727  26.982 > 2
   62 Al  0.41666666666667  0.83333333333333  0.75036158239727  26.982 > 2
   63 Al  0.66666666666667  0.83333333333333  0.75036158239727  26.982 > 2
   64 Al  0.91666666666667  0.83333333333333  0.75036158239727  26.982 > 2
  *65 As  0.08333333333333  0.16666666666667  0.18763841760273  74.922 > 3
   66 As  0.33333333333333  0.16666666666667  0.18763841760273  74.922 > 3
   67 As  0.58333333333333  0.16666666666667  0.18763841760273  74.922 > 3
   68 As  0.83333333333333  0.16666666666667  0.18763841760273  74.922 > 3
   69 As  0.08333333333333  0.41666666666667  0.18763841760273  74.922 > 3
   70 As  0.33333333333333  0.41666666666667  0.18763841760273  74.922 > 3
   71 As  0.58333333333333  0.41666666666667  0.18763841760273  74.922 > 3
   72 As  0.83333333333333  0.41666666666667  0.18763841760273  74.922 > 3
   73 As  0.08333333333333  0.66666666666667  0.18763841760273  74.922 > 3
   74 As  0.33333333333333  0.66666666666667  0.18763841760273  74.922 > 3
   75 As  0.58333333333333  0.66666666666667  0.18763841760273  74.922 > 3
   76 As  0.83333333333333  0.66666666666667  0.18763841760273  74.922 > 3
   77 As  0.08333333333333  0.91666666666667  0.18763841760273  74.922 > 3
   78 As  0.33333333333333  0.91666666666667  0.18763841760273  74.922 > 3
   79 As  0.58333333333333  0.91666666666667  0.18763841760273  74.922 > 3
   80 As  0.83333333333333  0.91666666666667  0.18763841760273  74.922 > 3
   81 As  0.08333333333333  0.16666666666667  0.68763841760273  74.922 > 3
   82 As  0.33333333333333  0.16666666666667  0.68763841760273  74.922 > 3
   83 As  0.58333333333333  0.16666666666667  0.68763841760273  74.922 > 3
   84 As  0.83333333333333  0.16666666666667  0.68763841760273  74.922 > 3
   85 As  0.08333333333333  0.41666666666667  0.68763841760273  74.922 > 3
   86 As  0.33333333333333  0.41666666666667  0.68763841760273  74.922 > 3
   87 As  0.58333333333333  0.41666666666667  0.68763841760273  74.922 > 3
   88 As  0.83333333333333  0.41666666666667  0.68763841760273  74.922 > 3
   89 As  0.08333333333333  0.66666666666667  0.68763841760273  74.922 > 3
   90 As  0.33333333333333  0.66666666666667  0.68763841760273  74.922 > 3
   91 As  0.58333333333333  0.66666666666667  0.68763841760273  74.922 > 3
   92 As  0.83333333333333  0.66666666666667  0.68763841760273  74.922 > 3
   93 As  0.08333333333333  0.91666666666667  0.68763841760273  74.922 > 3
   94 As  0.33333333333333  0.91666666666667  0.68763841760273  74.922 > 3
   95 As  0.58333333333333  0.91666666666667  0.68763841760273  74.922 > 3
   96 As  0.83333333333333  0.91666666666667  0.68763841760273  74.922 > 3
   97 As  0.16666666666667  0.08333333333333  0.43763841760273  74.922 > 4
   98 As  0.41666666666667  0.08333333333333  0.43763841760273  74.922 > 4
   99 As  0.66666666666667  0.08333333333333  0.43763841760273  74.922 > 4
  100 As  0.91666666666667  0.08333333333333  0.43763841760273  74.922 > 4
  101 As  0.16666666666667  0.33333333333333  0.43763841760273  74.922 > 4
  102 As  0.41666666666667  0.33333333333333  0.43763841760273  74.922 > 4
  103 As  0.66666666666667  0.33333333333333  0.43763841760273  74.922 > 4
  104 As  0.91666666666667  0.33333333333333  0.43763841760273  74.922 > 4
  105 As  0.16666666666667  0.58333333333333  0.43763841760273  74.922 > 4
  106 As  0.41666666666667  0.58333333333333  0.43763841760273  74.922 > 4
  107 As  0.66666666666667  0.58333333333333  0.43763841760273  74.922 > 4
  108 As  0.91666666666667  0.58333333333333  0.43763841760273  74.922 > 4
  109 As  0.16666666666667  0.83333333333333  0.43763841760273  74.922 > 4
  110 As  0.41666666666667  0.83333333333333  0.43763841760273  74.922 > 4
  111 As  0.66666666666667  0.83333333333333  0.43763841760273  74.922 > 4
  112 As  0.91666666666667  0.83333333333333  0.43763841760273  74.922 > 4
  113 As  0.16666666666667  0.08333333333333  0.93763841760273  74.922 > 4
  114 As  0.41666666666667  0.08333333333333  0.93763841760273  74.922 > 4
  115 As  0.66666666666667  0.08333333333333  0.93763841760273  74.922 > 4
  116 As  0.91666666666667  0.08333333333333  0.93763841760273  74.922 > 4
  117 As  0.16666666666667  0.33333333333333  0.93763841760273  74.922 > 4
  118 As  0.41666666666667  0.33333333333333  0.93763841760273  74.922 > 4
  119 As  0.66666666666667  0.33333333333333  0.93763841760273  74.922 > 4
  120 As  0.91666666666667  0.33333333333333  0.93763841760273  74.922 > 4
  121 As  0.16666666666667  0.58333333333333  0.93763841760273  74.922 > 4
  122 As  0.41666666666667  0.58333333333333  0.93763841760273  74.922 > 4
  123 As  0.66666666666667  0.58333333333333  0.93763841760273  74.922 > 4
  124 As  0.91666666666667  0.58333333333333  0.93763841760273  74.922 > 4
  125 As  0.16666666666667  0.83333333333333  0.93763841760273  74.922 > 4
  126 As  0.41666666666667  0.83333333333333  0.93763841760273  74.922 > 4
  127 As  0.66666666666667  0.83333333333333  0.93763841760273  74.922 > 4
  128 As  0.91666666666667  0.83333333333333  0.93763841760273  74.922 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           10.7659988   -0.0000000    0.0000000
            0.0000000   10.7659988    0.0000000
            0.0000000    0.0000000   10.7905166
-------------------------- Born effective charges --------------------------
    1 Al    2.1721112   -0.0000000    0.0000000
            0.0000000    2.1721112    0.0000000
            0.0000000    0.0000000    2.3242496
    2 Al    2.1721112   -0.0000000    0.0000000
            0.0000000    2.1721112    0.0000000
            0.0000000    0.0000000    2.3242496
    3 As   -2.1721112    0.0000000    0.0000000
            0.0000000   -2.1721112    0.0000000
            0.0000000    0.0000000   -2.3242496
    4 As   -2.1721112    0.0000000    0.0000000
            0.0000000   -2.1721112    0.0000000
            0.0000000    0.0000000   -2.3242496
----------------------------------------------------------------------------
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.167
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.178
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 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 (zz) -0.00000000 (zz) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-07 23:04:09]-------------------------
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.001504419534202    0.000000000000000    0.000000000000000
  b   -2.000752209767101    3.465404480672323    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.574457220000000
Atomic positions (fractional):
    1 Al  0.33333333333333  0.66666666666667  0.00072316479453  26.982
    2 Al  0.66666666666667  0.33333333333333  0.50072316479453  26.982
    3 As  0.33333333333333  0.66666666666667  0.37527683520547  74.922
    4 As  0.66666666666667  0.33333333333333  0.87527683520547  74.922
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.006017678136807    0.000000000000000    0.000000000000000
  b   -8.003008839068404   13.861617922689293    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.148914440000000
Atomic positions (fractional):
    1 Al  0.08333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    2 Al  0.33333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    3 Al  0.58333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    4 Al  0.83333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    5 Al  0.08333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    6 Al  0.33333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    7 Al  0.58333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    8 Al  0.83333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    9 Al  0.08333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   10 Al  0.33333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   11 Al  0.58333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   12 Al  0.83333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   13 Al  0.08333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   14 Al  0.33333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   15 Al  0.58333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   16 Al  0.83333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   17 Al  0.08333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   18 Al  0.33333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   19 Al  0.58333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   20 Al  0.83333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   21 Al  0.08333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   22 Al  0.33333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   23 Al  0.58333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   24 Al  0.83333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   25 Al  0.08333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   26 Al  0.33333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   27 Al  0.58333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   28 Al  0.83333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   29 Al  0.08333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   30 Al  0.33333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   31 Al  0.58333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   32 Al  0.83333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   33 Al  0.16666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   34 Al  0.41666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   35 Al  0.66666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   36 Al  0.91666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   37 Al  0.16666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   38 Al  0.41666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   39 Al  0.66666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   40 Al  0.91666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   41 Al  0.16666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   42 Al  0.41666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   43 Al  0.66666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   44 Al  0.91666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   45 Al  0.16666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   46 Al  0.41666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   47 Al  0.66666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   48 Al  0.91666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   49 Al  0.16666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   50 Al  0.41666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   51 Al  0.66666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   52 Al  0.91666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   53 Al  0.16666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   54 Al  0.41666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   55 Al  0.66666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   56 Al  0.91666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   57 Al  0.16666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   58 Al  0.41666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   59 Al  0.66666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   60 Al  0.91666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   61 Al  0.16666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   62 Al  0.41666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   63 Al  0.66666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   64 Al  0.91666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   65 As  0.08333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   66 As  0.33333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   67 As  0.58333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   68 As  0.83333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   69 As  0.08333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   70 As  0.33333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   71 As  0.58333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   72 As  0.83333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   73 As  0.08333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   74 As  0.33333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   75 As  0.58333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   76 As  0.83333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   77 As  0.08333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   78 As  0.33333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   79 As  0.58333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   80 As  0.83333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   81 As  0.08333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   82 As  0.33333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   83 As  0.58333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   84 As  0.83333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   85 As  0.08333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   86 As  0.33333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   87 As  0.58333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   88 As  0.83333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   89 As  0.08333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   90 As  0.33333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   91 As  0.58333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   92 As  0.83333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   93 As  0.08333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   94 As  0.33333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   95 As  0.58333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   96 As  0.83333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   97 As  0.16666666666667  0.08333333333333  0.43763841760273  74.922 > 97
   98 As  0.41666666666667  0.08333333333333  0.43763841760273  74.922 > 97
   99 As  0.66666666666667  0.08333333333333  0.43763841760273  74.922 > 97
  100 As  0.91666666666667  0.08333333333333  0.43763841760273  74.922 > 97
  101 As  0.16666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  102 As  0.41666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  103 As  0.66666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  104 As  0.91666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  105 As  0.16666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  106 As  0.41666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  107 As  0.66666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  108 As  0.91666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  109 As  0.16666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  110 As  0.41666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  111 As  0.66666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  112 As  0.91666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  113 As  0.16666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  114 As  0.41666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  115 As  0.66666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  116 As  0.91666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  117 As  0.16666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  118 As  0.41666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  119 As  0.66666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  120 As  0.91666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  121 As  0.16666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  122 As  0.41666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  123 As  0.66666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  124 As  0.91666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  125 As  0.16666666666667  0.83333333333333  0.93763841760273  74.922 > 97
  126 As  0.41666666666667  0.83333333333333  0.93763841760273  74.922 > 97
  127 As  0.66666666666667  0.83333333333333  0.93763841760273  74.922 > 97
  128 As  0.91666666666667  0.83333333333333  0.93763841760273  74.922 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           10.7659988   -0.0000000    0.0000000
            0.0000000   10.7659988    0.0000000
            0.0000000    0.0000000   10.7905166
-------------------------- Born effective charges --------------------------
    1 Al    2.1721112   -0.0000000    0.0000000
            0.0000000    2.1721112    0.0000000
            0.0000000    0.0000000    2.3242496
    2 Al    2.1721112   -0.0000000    0.0000000
            0.0000000    2.1721112    0.0000000
            0.0000000    0.0000000    2.3242496
    3 As   -2.1721112    0.0000000    0.0000000
            0.0000000   -2.1721112    0.0000000
            0.0000000    0.0000000   -2.3242496
    4 As   -2.1721112    0.0000000    0.0000000
            0.0000000   -2.1721112    0.0000000
            0.0000000    0.0000000   -2.3242496
----------------------------------------------------------------------------
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.00000119 (zzz) -0.00000119 (zzz) -0.00000119 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-07 23:04:14]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 23:04: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.001504419534202    0.000000000000000    0.000000000000000
  b   -2.000752209767101    3.465404480672323    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.574457220000000
Atomic positions (fractional):
    1 Al  0.33333333333333  0.66666666666667  0.00072316479453  26.982
    2 Al  0.66666666666667  0.33333333333333  0.50072316479453  26.982
    3 As  0.33333333333333  0.66666666666667  0.37527683520547  74.922
    4 As  0.66666666666667  0.33333333333333  0.87527683520547  74.922
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.006017678136807    0.000000000000000    0.000000000000000
  b   -8.003008839068404   13.861617922689293    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.148914440000000
Atomic positions (fractional):
    1 Al  0.08333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    2 Al  0.33333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    3 Al  0.58333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    4 Al  0.83333333333333  0.16666666666667  0.00036158239727  26.982 > 1
    5 Al  0.08333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    6 Al  0.33333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    7 Al  0.58333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    8 Al  0.83333333333333  0.41666666666667  0.00036158239727  26.982 > 1
    9 Al  0.08333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   10 Al  0.33333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   11 Al  0.58333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   12 Al  0.83333333333333  0.66666666666667  0.00036158239727  26.982 > 1
   13 Al  0.08333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   14 Al  0.33333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   15 Al  0.58333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   16 Al  0.83333333333333  0.91666666666667  0.00036158239727  26.982 > 1
   17 Al  0.08333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   18 Al  0.33333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   19 Al  0.58333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   20 Al  0.83333333333333  0.16666666666667  0.50036158239727  26.982 > 1
   21 Al  0.08333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   22 Al  0.33333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   23 Al  0.58333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   24 Al  0.83333333333333  0.41666666666667  0.50036158239727  26.982 > 1
   25 Al  0.08333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   26 Al  0.33333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   27 Al  0.58333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   28 Al  0.83333333333333  0.66666666666667  0.50036158239727  26.982 > 1
   29 Al  0.08333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   30 Al  0.33333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   31 Al  0.58333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   32 Al  0.83333333333333  0.91666666666667  0.50036158239727  26.982 > 1
   33 Al  0.16666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   34 Al  0.41666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   35 Al  0.66666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   36 Al  0.91666666666667  0.08333333333333  0.25036158239727  26.982 > 33
   37 Al  0.16666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   38 Al  0.41666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   39 Al  0.66666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   40 Al  0.91666666666667  0.33333333333333  0.25036158239727  26.982 > 33
   41 Al  0.16666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   42 Al  0.41666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   43 Al  0.66666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   44 Al  0.91666666666667  0.58333333333333  0.25036158239727  26.982 > 33
   45 Al  0.16666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   46 Al  0.41666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   47 Al  0.66666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   48 Al  0.91666666666667  0.83333333333333  0.25036158239727  26.982 > 33
   49 Al  0.16666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   50 Al  0.41666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   51 Al  0.66666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   52 Al  0.91666666666667  0.08333333333333  0.75036158239727  26.982 > 33
   53 Al  0.16666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   54 Al  0.41666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   55 Al  0.66666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   56 Al  0.91666666666667  0.33333333333333  0.75036158239727  26.982 > 33
   57 Al  0.16666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   58 Al  0.41666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   59 Al  0.66666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   60 Al  0.91666666666667  0.58333333333333  0.75036158239727  26.982 > 33
   61 Al  0.16666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   62 Al  0.41666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   63 Al  0.66666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   64 Al  0.91666666666667  0.83333333333333  0.75036158239727  26.982 > 33
   65 As  0.08333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   66 As  0.33333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   67 As  0.58333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   68 As  0.83333333333333  0.16666666666667  0.18763841760273  74.922 > 65
   69 As  0.08333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   70 As  0.33333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   71 As  0.58333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   72 As  0.83333333333333  0.41666666666667  0.18763841760273  74.922 > 65
   73 As  0.08333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   74 As  0.33333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   75 As  0.58333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   76 As  0.83333333333333  0.66666666666667  0.18763841760273  74.922 > 65
   77 As  0.08333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   78 As  0.33333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   79 As  0.58333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   80 As  0.83333333333333  0.91666666666667  0.18763841760273  74.922 > 65
   81 As  0.08333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   82 As  0.33333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   83 As  0.58333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   84 As  0.83333333333333  0.16666666666667  0.68763841760273  74.922 > 65
   85 As  0.08333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   86 As  0.33333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   87 As  0.58333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   88 As  0.83333333333333  0.41666666666667  0.68763841760273  74.922 > 65
   89 As  0.08333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   90 As  0.33333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   91 As  0.58333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   92 As  0.83333333333333  0.66666666666667  0.68763841760273  74.922 > 65
   93 As  0.08333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   94 As  0.33333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   95 As  0.58333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   96 As  0.83333333333333  0.91666666666667  0.68763841760273  74.922 > 65
   97 As  0.16666666666667  0.08333333333333  0.43763841760273  74.922 > 97
   98 As  0.41666666666667  0.08333333333333  0.43763841760273  74.922 > 97
   99 As  0.66666666666667  0.08333333333333  0.43763841760273  74.922 > 97
  100 As  0.91666666666667  0.08333333333333  0.43763841760273  74.922 > 97
  101 As  0.16666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  102 As  0.41666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  103 As  0.66666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  104 As  0.91666666666667  0.33333333333333  0.43763841760273  74.922 > 97
  105 As  0.16666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  106 As  0.41666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  107 As  0.66666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  108 As  0.91666666666667  0.58333333333333  0.43763841760273  74.922 > 97
  109 As  0.16666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  110 As  0.41666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  111 As  0.66666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  112 As  0.91666666666667  0.83333333333333  0.43763841760273  74.922 > 97
  113 As  0.16666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  114 As  0.41666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  115 As  0.66666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  116 As  0.91666666666667  0.08333333333333  0.93763841760273  74.922 > 97
  117 As  0.16666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  118 As  0.41666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  119 As  0.66666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  120 As  0.91666666666667  0.33333333333333  0.93763841760273  74.922 > 97
  121 As  0.16666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  122 As  0.41666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  123 As  0.66666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  124 As  0.91666666666667  0.58333333333333  0.93763841760273  74.922 > 97
  125 As  0.16666666666667  0.83333333333333  0.93763841760273  74.922 > 97
  126 As  0.41666666666667  0.83333333333333  0.93763841760273  74.922 > 97
  127 As  0.66666666666667  0.83333333333333  0.93763841760273  74.922 > 97
  128 As  0.91666666666667  0.83333333333333  0.93763841760273  74.922 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           10.7659988   -0.0000000    0.0000000
            0.0000000   10.7659988    0.0000000
            0.0000000    0.0000000   10.7905166
-------------------------- Born effective charges --------------------------
    1 Al    2.1721112   -0.0000000    0.0000000
            0.0000000    2.1721112    0.0000000
            0.0000000    0.0000000    2.3242496
    2 Al    2.1721112   -0.0000000    0.0000000
            0.0000000    2.1721112    0.0000000
            0.0000000    0.0000000    2.3242496
    3 As   -2.1721112    0.0000000    0.0000000
            0.0000000   -2.1721112    0.0000000
            0.0000000    0.0000000   -2.3242496
    4 As   -2.1721112    0.0000000    0.0000000
            0.0000000   -2.1721112    0.0000000
            0.0000000    0.0000000   -2.3242496
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000119 (zzz) -0.00000119 (zzz) -0.00000119 (zzz)
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 14 14 8 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.92, Number of G-points: 313, Lambda: 0.32
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/120) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 120
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   1.959   (  -0.000    0.000    0.000)    0.000
   1.959   (  -0.000   -0.000    0.000)    0.000
   6.195   (  -0.000   -0.000    0.000)    0.000
  10.268   (  -0.000    0.000   -0.000)    0.000
  10.268   (   0.000    0.000    0.000)    0.000
  10.437   (   0.000    0.000    0.000)    0.000
  10.606   (  -0.000   -0.000    0.000)    0.000
  10.606   (   0.000    0.000    0.000)    0.000
  10.721   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/120) =======================
q-point: ( 0.07  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.563   (  23.293   13.448    0.000)   26.896
   0.571   (  24.559   14.179    0.000)   28.358
   1.213   (  49.538   28.601    0.000)   57.202
   2.040   (   6.497    3.751    0.000)    7.502
   2.307   (  26.159   15.103    0.000)   30.206
   6.133   (  -5.077   -2.931    0.000)    5.862
  10.212   (  -4.179   -2.413    0.000)    4.826
  10.283   (   0.289    0.167    0.000)    0.334
  10.406   (  -2.162   -1.248    0.000)    2.496
  10.595   (  -1.265   -0.730    0.000)    1.461
  10.621   (  -7.669   -4.428    0.000)    8.855
  11.568   (  -0.281   -0.162    0.000)    0.324
======================= Grid point 2 (3/120) =======================
q-point: ( 0.14  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.095   (  20.953   12.097    0.000)   24.194
   1.146   (  23.785   13.732    0.000)   27.464
   2.245   (  10.090    5.825    0.000)   11.651
   2.343   (  44.916   25.932    0.000)   51.864
   3.032   (  31.719   18.313    0.000)   36.626
   5.969   (  -8.080   -4.665    0.000)    9.330
  10.096   (  -5.114   -2.952    0.000)    5.905
  10.193   (  -8.664   -5.002    0.000)   10.004
  10.323   (  -6.738   -3.890    0.000)    7.780
  10.483   (  -0.352   -0.203    0.000)    0.406
  10.533   (  -3.961   -2.287    0.000)    4.574
  11.557   (  -0.784   -0.452    0.000)    0.905
======================= Grid point 3 (4/120) =======================
q-point: ( 0.21  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.540   (  16.006    9.241    0.000)   18.482
   1.701   (  22.774   13.148    0.000)   26.297
   2.487   (   9.549    5.513    0.000)   11.026
   3.323   (  36.555   21.105    0.000)   42.210
   3.710   (  23.942   13.823    0.000)   27.646
   5.792   (  -5.569   -3.215    0.000)    6.430
   9.936   ( -10.662   -6.155    0.000)   12.311
   9.981   (  -4.460   -2.575    0.000)    5.150
  10.163   (  -6.143   -3.547    0.000)    7.093
  10.419   (  -5.091   -2.939    0.000)    5.879
  10.512   (   1.122    0.648    0.000)    1.295
  11.522   (  -2.375   -1.371    0.000)    2.742
======================= Grid point 4 (5/120) =======================
q-point: ( 0.29  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.842   (   9.269    5.351    0.000)   10.703
   2.224   (  20.981   12.113    0.000)   24.227
   2.670   (   5.385    3.109    0.000)    6.218
   3.927   (   9.967    5.754    0.000)   11.509
   4.256   (  25.904   14.955    0.000)   29.911
   5.757   (   2.877    1.661    0.000)    3.322
   9.750   (  -4.394   -2.537    0.000)    5.074
   9.886   (  -3.476   -2.007    0.000)    4.013
  10.036   (  -4.328   -2.499    0.000)    4.997
  10.314   (  -3.368   -1.944    0.000)    3.889
  10.514   (  -0.888   -0.513    0.000)    1.025
  11.433   (  -5.360   -3.095    0.000)    6.190
======================= Grid point 5 (6/120) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.987   (   3.327    1.921    0.000)    3.842
   2.689   (  17.758   10.253    0.000)   20.505
   2.734   (   0.165    0.095    0.000)    0.191
   3.961   (  -3.521   -2.033    0.000)    4.066
   4.840   (  21.436   12.376    0.000)   24.752
   5.877   (   5.579    3.221    0.000)    6.442
   9.726   (   1.907    1.101    0.000)    2.202
   9.819   (  -2.136   -1.233    0.000)    2.467
   9.967   (  -1.237   -0.714    0.000)    1.428
  10.266   (  -0.774   -0.447    0.000)    0.894
  10.478   (  -1.924   -1.111    0.000)    2.221
  11.266   (  -8.455   -4.881    0.000)    9.763
======================= Grid point 6 (7/120) =======================
q-point: ( 0.43  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.026   (   0.467    0.270    0.000)    0.540
   2.701   (  -2.171   -1.253    0.000)    2.507
   3.050   (  11.933    6.890    0.000)   13.780
   3.847   (  -5.023   -2.900    0.000)    5.800
   5.249   (  12.561    7.252    0.000)   14.504
   5.971   (   1.939    1.119    0.000)    2.239
   9.784   (  -0.835   -0.482    0.000)    0.964
   9.790   (   2.208    1.275    0.000)    2.550
   9.993   (   3.387    1.955    0.000)    3.911
  10.265   (   0.405    0.234    0.000)    0.467
  10.430   (  -2.006   -1.158    0.000)    2.316
  11.057   (  -8.037   -4.640    0.000)    9.280
======================= Grid point 7 (8/120) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 160
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.029   (   0.000    0.000    0.000)    0.000
   2.669   (   0.000    0.000    0.000)    0.000
   3.207   (   0.000    0.000    0.000)    0.000
   3.773   (   0.000    0.000    0.000)    0.000
   5.407   (   0.000    0.000    0.000)    0.000
   5.986   (   0.000    0.000    0.000)    0.000
   9.775   (   0.000    0.000    0.000)    0.000
   9.814   (   0.000    0.000    0.000)    0.000
  10.059   (   0.000    0.000    0.000)    0.000
  10.272   (   0.000    0.000    0.000)    0.000
  10.398   (   0.000    0.000    0.000)    0.000
  10.947   (   0.000    0.000    0.000)    0.000
======================= Grid point 16 (9/120) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.993   (  13.790   23.885    0.000)   27.579
   1.031   (  15.272   26.452    0.000)   30.544
   1.977   (  21.708   37.599    0.000)   43.415
   2.283   (  11.463   19.854    0.000)   22.925
   2.774   (  16.754   29.019    0.000)   33.508
   6.019   (  -4.426   -7.665    0.000)    8.851
  10.150   (  -3.027   -5.243    0.000)    6.054
  10.202   (  -3.704   -6.415    0.000)    7.408
  10.354   (  -1.829   -3.168    0.000)    3.658
  10.483   (  -4.006   -6.939    0.000)    8.012
  10.612   (   0.314    0.545    0.000)    0.629
  11.540   (  -1.364   -2.363    0.000)    2.729
======================= Grid point 17 (10/120) =======================
q-point: ( 0.14  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.503   (  17.385   19.922    0.000)   26.441
   1.548   (  10.081   27.993    0.000)   29.753
   2.475   (   8.538   16.568    0.000)   18.639
   2.980   (  28.138   32.309    0.000)   42.845
   3.404   (  24.348   16.911    0.000)   29.645
   5.845   (  -5.643   -6.757    0.000)    8.803
  10.007   (  -6.790   -9.509    0.000)   11.684
  10.029   (  -4.531   -4.912    0.000)    6.682
  10.255   (  -6.542   -2.002    0.000)    6.842
  10.403   (   1.415   -5.255    0.000)    5.442
  10.589   (  -4.929    3.255    0.000)    5.906
  11.480   (   0.001   -6.810    0.000)    6.810
======================= Grid point 18 (11/120) =======================
q-point: ( 0.21  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.937   (   2.216   26.506    0.000)   26.599
   2.016   (  18.103   16.951    0.000)   24.801
   2.743   (   1.296   18.170    0.000)   18.217
   3.726   (  21.728   15.365    0.000)   26.612
   3.945   (  23.567   10.637    0.000)   25.856
   5.731   (  -0.510   -2.357    0.000)    2.411
   9.802   (  -5.648   -5.632    0.000)    7.976
   9.926   (  -3.056   -2.551    0.000)    3.981
  10.112   (  -6.885   -1.625    0.000)    7.074
  10.364   (  -1.332   -3.211    0.000)    3.476
  10.523   (  -3.939    0.203    0.000)    3.944
  11.395   (   0.290  -10.799    0.000)   10.803
======================= Grid point 19 (12/120) =======================
q-point: ( 0.29  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.164   (  -5.547   23.724    0.000)   24.364
   2.494   (  17.353   13.763    0.000)   22.148
   2.886   (  -5.498   16.880    0.000)   17.753
   3.963   (   2.145   -1.094    0.000)    2.408
   4.543   (  25.526   12.473    0.000)   28.410
   5.784   (   6.208    0.525    0.000)    6.230
   9.715   (  -0.417   -0.659    0.000)    0.780
   9.862   (  -2.689   -0.263    0.000)    2.702
   9.991   (  -3.963   -2.124    0.000)    4.496
  10.306   (  -2.710   -0.030    0.000)    2.710
  10.458   (  -0.981   -4.516    0.000)    4.621
  11.258   (  -1.761  -13.485    0.000)   13.600
======================= Grid point 20 (13/120) =======================
q-point: ( 0.36  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.248   ( -10.391   21.413    0.000)   23.801
   2.904   (  12.646   10.847    0.000)   16.660
   2.905   (  -9.108   15.286    0.000)   17.793
   3.905   (  -4.350   -2.933    0.000)    5.246
   5.046   (  18.599    7.750    0.000)   20.149
   5.898   (   5.655   -1.348    0.000)    5.813
   9.755   (   3.147    1.586    0.000)    3.524
   9.818   (  -2.590    0.847    0.000)    2.725
   9.948   (   1.606   -1.200    0.000)    2.004
  10.276   (  -1.033    0.656    0.000)    1.224
  10.394   (   0.374   -6.722    0.000)    6.732
  11.062   (  -3.473  -14.297    0.000)   14.713
======================= Grid point 21 (14/120) =======================
q-point: ( 0.43  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.261   ( -11.839   20.641    0.000)   23.795
   2.868   ( -11.144   16.386    0.000)   19.816
   3.180   (   6.041    5.690    0.000)    8.299
   3.797   (  -3.524   -1.789    0.000)    3.953
   5.344   (   8.122    1.934    0.000)    8.349
   5.942   (   2.669   -3.804    0.000)    4.647
   9.790   (  -1.133    0.943    0.000)    1.474
   9.811   (   0.822    0.741    0.000)    1.107
  10.017   (   4.692    0.192    0.000)    4.696
  10.274   (  -0.022    0.188    0.000)    0.189
  10.333   (   1.765   -8.160    0.000)    8.349
  10.885   (   0.254  -11.374    0.000)   11.377
======================= Grid point 31 (15/120) =======================
q-point: ( 0.14  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.932   (  12.295   21.296    0.000)   24.591
   2.116   (  14.362   24.876    0.000)   28.725
   2.852   (  10.465   18.125    0.000)   20.929
   3.630   (  17.294   29.954    0.000)   34.588
   3.708   (   7.302   12.647    0.000)   14.604
   5.725   (  -2.347   -4.065    0.000)    4.694
   9.824   (  -4.169   -7.221    0.000)    8.338
   9.939   (  -2.181   -3.777    0.000)    4.362
  10.213   (  -2.060   -3.567    0.000)    4.119
  10.295   (  -2.349   -4.068    0.000)    4.697
  10.610   (  -0.886   -1.535    0.000)    1.772
  11.302   (  -5.597   -9.695    0.000)   11.194
======================= Grid point 32 (16/120) =======================
q-point: ( 0.21  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.366   (  13.543   17.520    0.000)   22.145
   2.524   (   1.137   27.050    0.000)   27.074
   3.161   (   0.810   20.567    0.000)   20.583
   3.883   (   4.973    1.960    0.000)    5.346
   4.239   (  22.426   16.915    0.000)   28.090
   5.700   (   2.347   -0.717    0.000)    2.454
   9.722   (  -1.697   -2.231    0.000)    2.803
   9.883   (  -1.033   -1.597    0.000)    1.902
  10.088   (  -7.926   -0.849    0.000)    7.971
  10.251   (   2.986   -7.662    0.000)    8.223
  10.529   (  -6.773    0.167    0.000)    6.775
  11.104   (  -1.330  -16.516    0.000)   16.570
======================= Grid point 33 (17/120) =======================
q-point: ( 0.29  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 500
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.698   (  -5.937   24.146    0.000)   24.865
   2.795   (   9.825   16.871    0.000)   19.523
   3.298   (  -9.008   21.284    0.000)   23.112
   3.921   (  -0.139   -2.384    0.000)    2.388
   4.770   (  20.986    9.937    0.000)   23.219
   5.768   (   5.758   -2.124    0.000)    6.138
   9.718   (   2.059    0.677    0.000)    2.167
   9.865   (  -0.413    0.215    0.000)    0.466
   9.949   (  -5.263   -1.831    0.000)    5.572
  10.225   (   3.101   -8.178    0.000)    8.746
  10.403   (  -6.631   -0.263    0.000)    6.636
  10.909   (   0.027  -19.331    0.000)   19.331
======================= Grid point 34 (18/120) =======================
q-point: ( 0.36  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.762   ( -13.850   25.979    0.000)   29.441
   3.125   (   9.081   11.123    0.000)   14.359
   3.311   ( -13.480   21.830    0.000)   25.656
   3.848   (  -4.218   -2.578    0.000)    4.943
   5.169   (  14.445    4.469    0.000)   15.121
   5.824   (   4.477   -5.694    0.000)    7.243
   9.789   (   3.228    1.641    0.000)    3.621
   9.837   (  -2.993    0.737    0.000)    3.083
   9.917   (   4.278   -1.746    0.000)    4.620
  10.186   (   2.472  -10.004    0.000)   10.304
  10.300   (  -3.346   -1.864    0.000)    3.830
  10.720   (   0.970  -17.237    0.000)   17.264
======================= Grid point 35 (19/120) =======================
q-point: ( 0.43  0.14  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.769   ( -15.055   26.077    0.000)   30.111
   3.284   (  -2.566    4.444    0.000)    5.131
   3.298   ( -13.315   23.063    0.000)   26.631
   3.777   (   0.028   -0.048    0.000)    0.056
   5.333   (   1.587   -2.749    0.000)    3.174
   5.825   (   4.072   -7.053    0.000)    8.144
   9.813   (  -0.624    1.081    0.000)    1.248
   9.826   (  -0.482    0.835    0.000)    0.964
   9.995   (   1.543   -2.672    0.000)    3.085
  10.119   (   6.309  -10.928    0.000)   12.618
  10.260   (   1.194   -2.068    0.000)    2.388
  10.632   (   6.802  -11.781    0.000)   13.603
======================= Grid point 47 (20/120) =======================
q-point: ( 0.21  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.724   (   9.821   17.011    0.000)   19.642
   2.983   (   9.566   16.568    0.000)   19.131
   3.525   (   7.830   13.562    0.000)   15.660
   3.890   (  -0.236   -0.409    0.000)    0.472
   4.604   (  10.514   18.211    0.000)   21.028
   5.697   (   0.163    0.282    0.000)    0.326
   9.702   (   0.102    0.176    0.000)    0.203
   9.861   (  -0.261   -0.452    0.000)    0.522
  10.069   (  -1.936   -3.354    0.000)    3.872
  10.076   (  -3.664   -6.346    0.000)    7.327
  10.510   (  -1.515   -2.624    0.000)    3.030
  10.758   (  -9.159  -15.863    0.000)   18.317
======================= Grid point 48 (21/120) =======================
q-point: ( 0.29  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.052   (   7.857   14.335    0.000)   16.347
   3.217   (  -2.077   18.585    0.000)   18.700
   3.697   (  -5.246   15.983    0.000)   16.822
   3.871   (  -0.618   -2.489    0.000)    2.564
   4.972   (  13.530   10.078    0.000)   16.871
   5.701   (   2.128   -4.185    0.000)    4.694
   9.730   (   2.977    0.181    0.000)    2.983
   9.863   (   0.723   -0.288    0.000)    0.778
   9.918   (  -5.916   -1.262    0.000)    6.049
  10.046   (   2.573   -7.686    0.000)    8.105
  10.373   (  -9.288   -5.743    0.000)   10.920
  10.526   (   1.714  -14.267    0.000)   14.369
======================= Grid point 49 (22/120) =======================
q-point: ( 0.36  0.21  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.224   (  -6.797   17.068    0.000)   18.371
   3.378   (   0.810   13.584    0.000)   13.609
   3.708   (  -7.836   10.756    0.000)   13.308
   3.815   (  -4.980    3.469    0.000)    6.069
   5.248   (   7.555    3.687    0.000)    8.407
   5.677   (   2.775   -8.107    0.000)    8.569
   9.813   (   5.165   -0.261    0.000)    5.172
   9.839   (  -1.896   -0.068    0.000)    1.897
   9.886   (   2.944   -0.856    0.000)    3.066
   9.988   (  -1.078   -7.808    0.000)    7.883
  10.176   (  -0.175  -11.287    0.000)   11.288
  10.457   (   1.739   -6.299    0.000)    6.535
======================= Grid point 62 (23/120) =======================
q-point: ( 0.29  0.29  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 285
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.330   (   7.020   12.158    0.000)   14.039
   3.492   (   4.902    8.491    0.000)    9.804
   3.787   (  -2.839   -4.917    0.000)    5.678
   3.946   (   4.006    6.939    0.000)    8.013
   5.188   (   5.997   10.387    0.000)   11.994
   5.602   (  -3.112   -5.390    0.000)    6.224
   9.736   (   0.736    1.275    0.000)    1.472
   9.865   (   0.230    0.398    0.000)    0.460
   9.885   (  -1.008   -1.745    0.000)    2.015
   9.941   (  -1.647   -2.852    0.000)    3.293
  10.106   (  -7.963  -13.793    0.000)   15.927
  10.415   (  -0.973   -1.685    0.000)    1.945
======================= Grid point 63 (24/120) =======================
q-point: ( 0.36  0.29  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.484   (  -5.304    9.187    0.000)   10.609
   3.608   (  -3.304    5.722    0.000)    6.608
   3.679   (   2.542   -4.402    0.000)    5.083
   4.013   (  -3.750    6.495    0.000)    7.500
   5.336   (  -2.712    4.698    0.000)    5.425
   5.512   (   4.031   -6.982    0.000)    8.062
   9.786   (   0.233   -0.404    0.000)    0.466
   9.859   (  -0.862    1.492    0.000)    1.723
   9.880   (   1.476   -3.250    0.000)    3.570
   9.880   (   0.484   -0.145    0.000)    0.505
   9.949   (   4.421   -7.657    0.000)    8.841
  10.398   (   0.566   -0.981    0.000)    1.133
======================= Grid point 211 (25/120) =======================
q-point: ( 0.00  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 96
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.522   (   0.000   -0.000   26.668)   26.668
   0.522   (  -0.000   -0.000   26.668)   26.668
   1.198   (   0.000   -0.000   62.145)   62.145
   1.957   (  -0.000    0.000   -0.377)    0.377
   1.957   (  -0.000    0.000   -0.377)    0.377
   6.067   (   0.000    0.000  -13.331)   13.331
  10.277   (   0.000   -0.000    0.934)    0.934
  10.277   (  -0.000   -0.000    0.934)    0.934
  10.589   (  -0.000    0.000   -1.678)    1.678
  10.589   (   0.000   -0.000   -1.678)    1.678
  10.768   (   0.000   -0.000    4.785)    4.785
  11.530   (   0.000   -0.000   -2.180)    2.180
======================= Grid point 212 (26/120) =======================
q-point: ( 0.07  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.769   (  17.102    9.874   18.168)   26.834
   0.995   (  17.108    9.877   29.054)   35.134
   1.549   (  34.709   20.039   35.558)   53.579
   2.039   (   6.584    3.801   -0.246)    7.607
   2.297   (  25.846   14.922   -0.825)   29.856
   6.005   (  -5.051   -2.916  -13.281)   14.506
  10.222   (  -4.118   -2.377    1.050)    4.869
  10.274   (  -0.817   -0.472    0.059)    0.946
  10.496   (  -4.651   -2.685    4.156)    6.791
  10.576   (  -1.414   -0.816   -1.901)    2.506
  10.678   (  -6.774   -3.911    5.754)    9.711
  11.525   (  -0.260   -0.150   -3.772)    3.784
======================= Grid point 213 (27/120) =======================
q-point: ( 0.14  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.214   (  18.911   10.918   11.563)   24.709
   1.372   (  16.682    9.631   22.318)   29.481
   2.245   (  10.118    5.842   -0.078)   11.683
   2.517   (  41.737   24.097   17.061)   51.124
   3.017   (  31.561   18.222   -1.340)   36.468
   5.845   (  -7.781   -4.492  -12.848)   15.677
  10.107   (  -5.147   -2.971    1.109)    6.045
  10.174   (  -8.100   -4.677   -1.226)    9.433
  10.372   (  -7.292   -4.210    3.678)    9.188
  10.511   (  -4.074   -2.352   -2.245)    5.212
  10.554   (  -1.042   -0.601    7.186)    7.286
  11.518   (  -0.452   -0.261   -3.931)    3.966
======================= Grid point 214 (28/120) =======================
q-point: ( 0.21  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.626   (  15.098    8.717    8.579)   19.431
   1.791   (  18.202   10.509   10.123)   23.328
   2.487   (   9.508    5.489   -0.089)   10.979
   3.432   (  32.732   18.898    9.166)   38.892
   3.701   (  25.596   14.778    0.372)   29.559
   5.684   (  -4.316   -2.492  -10.891)   11.977
   9.940   (  -9.823   -5.672    0.131)   11.344
   9.990   (  -4.568   -2.638    0.922)    5.355
  10.193   (  -6.940   -4.007    2.626)    8.433
  10.395   (  -5.147   -2.972   -2.434)    6.423
  10.573   (   1.079    0.623    6.474)    6.593
  11.494   (  -1.750   -1.010   -2.809)    3.461
======================= Grid point 215 (29/120) =======================
q-point: ( 0.29  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.912   (   8.810    5.087    7.085)   12.397
   2.223   (  17.756   10.252    0.916)   20.523
   2.669   (   5.372    3.102   -0.158)    6.206
   3.875   (   5.217    3.012   -6.479)    8.846
   4.353   (  28.349   16.367   10.344)   34.330
   5.693   (   5.124    2.958   -6.155)    8.537
   9.766   (  -4.267   -2.463    1.041)    5.035
   9.892   (  -3.554   -2.052    0.683)    4.160
  10.053   (  -4.687   -2.706    1.350)    5.578
  10.289   (  -3.449   -1.991   -2.582)    4.746
  10.580   (  -0.500   -0.289    7.067)    7.090
  11.424   (  -4.371   -2.523   -0.856)    5.119
======================= Grid point 216 (30/120) =======================
q-point: ( 0.36  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.050   (   3.095    1.787    6.310)    7.252
   2.622   (  15.465    8.929   -5.814)   18.780
   2.735   (   0.252    0.146   -0.066)    0.299
   3.844   (  -5.231   -3.020  -11.907)   13.351
   4.956   (  21.629   12.488   11.257)   27.395
   5.864   (   7.424    4.286   -1.245)    8.662
   9.735   (   1.167    0.674    0.412)    1.409
   9.824   (  -2.160   -1.247    0.570)    2.558
   9.978   (  -1.486   -0.858    0.569)    1.808
  10.238   (  -0.894   -0.516   -2.832)    3.014
  10.557   (  -1.135   -0.655    8.517)    8.618
  11.285   (  -7.095   -4.096    1.952)    8.422
======================= Grid point 217 (31/120) =======================
q-point: ( 0.43  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.084   (   0.351    0.203    5.909)    5.922
   2.704   (  -2.038   -1.177    0.231)    2.365
   2.940   (  10.674    6.163  -10.218)   16.010
   3.701   (  -5.861   -3.384  -14.505)   16.006
   5.367   (  12.555    7.248   11.258)   18.355
   5.993   (   3.115    1.798    2.163)    4.197
   9.779   (   1.607    0.928   -1.093)    2.153
   9.789   (  -0.837   -0.483    0.552)    1.113
   9.986   (   1.926    1.112   -1.284)    2.568
  10.235   (   0.314    0.181   -3.103)    3.124
  10.536   (  -0.556   -0.321   10.957)   10.976
  11.110   (  -6.744   -3.894    5.170)    9.347
======================= Grid point 218 (32/120) =======================
q-point: (-0.50  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 228
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.086   (   0.000    0.000    5.761)    5.761
   2.674   (   0.000    0.000    0.423)    0.423
   3.082   (   0.000    0.000  -11.778)   11.778
   3.616   (  -0.000   -0.000  -15.382)   15.382
   5.524   (   0.000    0.000   11.136)   11.136
   6.024   (  -0.000   -0.000    3.554)    3.554
   9.780   (   0.000    0.000    0.553)    0.553
   9.797   (   0.000    0.000   -1.472)    1.472
  10.022   (   0.000    0.000   -3.640)    3.640
  10.241   (   0.000    0.000   -3.220)    3.220
  10.530   (   0.000    0.000   12.882)   12.882
  11.018   (   0.000    0.000    6.831)    6.831
======================= Grid point 227 (33/120) =======================
q-point: ( 0.07  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.156   (  13.628   23.605   12.159)   29.846
   1.258   (   9.015   15.614   25.061)   30.873
   2.127   (  13.668   23.674   10.375)   29.238
   2.335   (  17.310   29.981    9.717)   35.957
   2.764   (  16.532   28.635   -1.121)   33.084
   5.893   (  -4.326   -7.492  -13.047)   15.655
  10.151   (  -2.931   -5.076    0.543)    5.886
  10.198   (  -3.872   -6.707   -0.368)    7.754
  10.400   (  -2.863   -4.959    3.269)    6.594
  10.553   (  -3.278   -5.678    5.939)    8.847
  10.593   (  -0.057   -0.098   -0.887)    0.895
  11.502   (  -1.146   -1.984   -3.779)    4.419
======================= Grid point 228 (34/120) =======================
q-point: ( 0.14  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.604   (  10.636   19.736   10.433)   24.728
   1.655   (  12.749   21.489   11.511)   27.510
   2.483   (   7.338   16.707    0.460)   18.254
   3.117   (  26.968   31.230   13.286)   43.349
   3.386   (  24.585   16.501   -1.553)   29.650
   5.728   (  -5.049   -6.138  -11.947)   14.349
  10.009   (  -6.563   -8.957   -0.044)   11.104
  10.033   (  -4.155   -4.832    0.614)    6.403
  10.281   (  -6.284   -3.169    2.047)    7.330
  10.457   (  -1.744   -4.046    4.628)    6.390
  10.580   (  -2.663    1.998    0.253)    3.339
  11.451   (   0.197   -6.109   -2.951)    6.787
======================= Grid point 229 (35/120) =======================
q-point: ( 0.21  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.974   (   6.339   20.357    4.014)   21.696
   2.072   (  11.202   18.671    6.507)   22.725
   2.743   (   1.092   18.256   -0.056)   18.288
   3.740   (  14.418    8.362   -1.915)   16.777
   4.005   (  28.644   16.354    8.702)   34.113
   5.644   (   1.231   -1.026   -8.559)    8.707
   9.819   (  -5.376   -4.948    1.118)    7.392
   9.931   (  -2.952   -2.782    0.572)    4.097
  10.129   (  -7.072   -2.141    1.384)    7.517
  10.374   (  -4.888   -0.683    0.051)    4.936
  10.551   (  -0.346   -2.473    4.190)    4.878
  11.380   (   0.786   -9.948   -1.494)   10.090
======================= Grid point 230 (36/120) =======================
q-point: ( 0.29  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.210   (  -4.037   21.504    4.324)   22.303
   2.460   (  13.854   12.822   -2.085)   18.992
   2.886   (  -5.501   16.854   -0.188)   17.730
   3.869   (  -0.116   -2.615   -9.854)   10.196
   4.658   (  25.724   13.037   11.351)   30.993
   5.748   (   8.249    1.923   -3.376)    9.118
   9.735   (  -0.969   -0.309    1.373)    1.709
   9.867   (  -2.663   -0.314    0.545)    2.736
  10.002   (  -4.030   -2.354    0.786)    4.733
  10.290   (  -4.189    1.099   -1.919)    4.737
  10.514   (   1.095   -5.838    6.619)    8.894
  11.265   (  -0.812  -12.543    0.740)   12.591
======================= Grid point 231 (37/120) =======================
q-point: ( 0.36  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.295   (  -9.662   20.184    4.688)   22.863
   2.811   (  11.904    9.502   -8.227)   17.311
   2.906   ( -10.112   15.479   -0.131)   18.489
   3.770   (  -5.560   -3.696  -13.470)   15.034
   5.165   (  18.535    7.943   11.430)   23.179
   5.907   (   7.104   -0.196    0.862)    7.158
   9.758   (   1.825    1.320   -0.080)    2.254
   9.822   (  -2.608    0.828    0.524)    2.786
   9.955   (   1.067   -1.381    0.023)    1.745
  10.252   (  -1.565    1.207   -2.581)    3.250
  10.476   (   2.222   -7.063    8.913)   11.587
  11.099   (  -2.185  -13.402    3.708)   14.076
======================= Grid point 232 (38/120) =======================
q-point: ( 0.43  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.308   ( -11.315   19.710    4.733)   23.214
   2.869   ( -10.872   16.108    0.021)   19.434
   3.058   (   5.391    5.252  -11.399)   13.660
   3.642   (  -3.945   -2.077  -15.170)   15.812
   5.462   (   8.012    2.012   11.244)   13.952
   5.979   (   3.166   -3.113    3.456)    5.627
   9.787   (  -0.081    0.316   -0.730)    0.799
   9.803   (  -0.481    1.267   -0.036)    1.356
   9.990   (   2.786   -0.936   -2.922)    4.144
  10.249   (  -0.353    0.879   -2.700)    2.862
  10.452   (   3.680   -7.655   11.956)   14.666
  10.949   (   1.044  -11.050    6.328)   12.776
======================= Grid point 242 (39/120) =======================
q-point: ( 0.14  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.958   (   9.951   17.236    4.000)   20.300
   2.179   (  13.901   24.077    6.365)   28.521
   2.861   (  10.404   18.020    0.581)   20.816
   3.661   (   6.385   11.059   -4.962)   13.700
   3.763   (  17.192   29.778   13.123)   36.804
   5.629   (  -1.456   -2.522   -9.525)    9.960
   9.842   (  -3.685   -6.383    1.309)    7.486
   9.941   (  -2.241   -3.882    0.192)    4.487
  10.218   (  -2.337   -4.048    0.317)    4.685
  10.363   (  -2.453   -4.248    6.812)    8.395
  10.588   (  -0.986   -1.709   -1.831)    2.691
  11.290   (  -5.097   -8.829   -1.289)   10.276
======================= Grid point 243 (40/120) =======================
q-point: ( 0.21  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.314   (  10.916   14.985   -3.684)   18.902
   2.582   (   2.136   25.680    5.534)   26.356
   3.165   (   0.309   20.782    0.206)   20.785
   3.802   (   3.915   -0.523   -8.564)    9.431
   4.365   (  21.882   17.908   12.352)   30.855
   5.646   (   4.035    1.090   -5.197)    6.669
   9.755   (  -1.700   -1.352    2.481)    3.297
   9.884   (  -0.842   -1.743    0.140)    1.941
  10.096   (  -7.501   -1.318    0.483)    7.631
  10.300   (   1.438   -6.981    5.155)    8.796
  10.511   (  -5.493   -1.006   -1.346)    5.744
  11.110   (  -0.908  -15.271    0.650)   15.312
======================= Grid point 244 (41/120) =======================
q-point: ( 0.29  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.629   (   5.187   15.472   -7.672)   18.032
   2.811   (  -2.264   22.485    3.431)   22.857
   3.295   (  -9.227   21.010   -0.562)   22.954
   3.801   (  -1.395   -3.404  -12.010)   12.561
   4.896   (  20.657   10.395   12.038)   26.071
   5.763   (   7.425   -0.500   -0.461)    7.456
   9.747   (   0.679    1.128    2.003)    2.397
   9.869   (  -0.275    0.079    0.376)    0.472
   9.958   (  -5.184   -1.760    0.603)    5.507
  10.260   (   1.961   -6.637    3.365)    7.696
  10.401   (  -4.289   -2.387    0.772)    4.969
  10.937   (   0.720  -18.238    2.875)   18.477
======================= Grid point 245 (42/120) =======================
q-point: ( 0.36  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.768   ( -10.123   22.302   -1.070)   24.515
   3.029   (   5.121   12.641   -7.116)   15.383
   3.302   ( -13.117   20.758   -1.332)   24.591
   3.704   (  -5.198   -2.420  -13.832)   14.973
   5.292   (  14.150    4.674   11.646)   18.913
   5.858   (   5.417   -4.340    3.198)    7.642
   9.782   (   1.192    0.772   -0.621)    1.550
   9.844   (  -2.413    1.043    0.597)    2.696
   9.921   (   2.534   -1.820   -0.705)    3.198
  10.238   (   3.185   -6.250    4.103)    8.127
  10.315   (  -1.431   -4.475    3.298)    5.740
  10.770   (   1.913  -17.206    5.141)   18.060
======================= Grid point 246 (43/120) =======================
q-point: ( 0.43  0.14  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.786   ( -13.531   23.437    0.790)   27.074
   3.164   (  -3.334    5.774  -10.211)   12.195
   3.285   ( -12.370   21.425   -1.998)   24.820
   3.625   (  -0.598    1.035  -14.178)   14.228
   5.454   (   1.482   -2.567   11.355)   11.735
   5.876   (   3.677   -6.369    4.720)    8.739
   9.790   (   0.058   -0.101   -1.440)    1.444
   9.835   (  -1.002    1.735    0.635)    2.102
   9.946   (   1.867   -3.234   -4.131)    5.569
  10.247   (   6.301  -10.913   12.280)   17.595
  10.262   (  -0.098    0.169   -0.995)    1.014
  10.691   (   7.506  -13.001    6.103)   16.205
======================= Grid point 258 (44/120) =======================
q-point: ( 0.21  0.21  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.619   (   8.434   14.608   -9.182)   19.205
   3.022   (   9.288   16.087    3.918)   18.984
   3.533   (   7.809   13.526    0.542)   15.628
   3.768   (  -1.200   -2.079  -12.302)   12.533
   4.739   (  10.570   18.308   12.814)   24.720
   5.684   (   1.416    2.453   -1.209)    3.080
   9.750   (   0.278    0.481    3.644)    3.686
   9.860   (  -0.304   -0.527   -0.109)    0.619
  10.055   (  -2.258   -3.911   -1.288)    4.696
  10.137   (  -3.697   -6.404    6.937)   10.139
  10.480   (  -1.593   -2.759   -2.921)    4.322
  10.794   (  -8.246  -14.282    3.541)   16.868
======================= Grid point 259 (45/120) =======================
q-point: ( 0.29  0.21  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.906   (   7.382   12.616  -13.308)   19.768
   3.250   (  -2.040   17.523    2.894)   17.877
   3.657   (  -6.966    8.355   -8.324)   13.698
   3.765   (   0.088    4.333   -5.912)    7.330
   5.105   (  13.129   10.237   12.414)   20.767
   5.733   (   3.356   -2.254    2.893)    4.971
   9.763   (   0.924   -0.168    2.243)    2.432
   9.859   (   0.548   -0.828   -0.424)    1.080
   9.936   (  -5.462   -0.568    0.524)    5.517
  10.075   (   3.378   -8.679    5.589)   10.862
  10.382   (  -6.893   -1.258   -0.768)    7.049
  10.561   (   0.713  -16.875    4.585)   17.501
======================= Grid point 260 (46/120) =======================
q-point: ( 0.36  0.21  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.106   (  -0.373   11.567  -13.632)   17.882
   3.363   (  -4.409   17.079    0.860)   17.660
   3.569   (  -5.888    0.107  -15.023)   16.136
   3.785   (  -7.960   13.699   -0.883)   15.868
   5.374   (   7.233    3.681   11.643)   14.193
   5.742   (   3.000   -6.587    5.853)    9.308
   9.776   (   1.147   -1.510   -1.312)    2.306
   9.856   (   1.080   -0.512   -1.037)    1.583
   9.884   (  -0.557   -1.485   -1.320)    2.063
  10.055   (   3.578   -7.274   10.047)   12.910
  10.255   (  -1.804   -5.163    4.443)    7.047
  10.462   (   3.338  -10.326    1.958)   11.027
======================= Grid point 273 (47/120) =======================
q-point: ( 0.29  0.29  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.155   (   6.406   11.095  -16.362)   20.781
   3.488   (   2.862    4.957   -4.735)    7.428
   3.643   (  -1.293   -2.240   -9.330)    9.682
   3.951   (   3.965    6.867    0.384)    7.939
   5.318   (   5.742    9.945   11.896)   16.534
   5.673   (  -2.102   -3.640    6.397)    7.654
   9.741   (  -0.569   -0.986    0.594)    1.284
   9.840   (  -1.098   -1.902   -2.976)    3.698
   9.899   (  -1.526   -2.644   -0.959)    3.200
   9.981   (  -0.032   -0.056    9.447)    9.447
  10.235   (  -6.218  -10.771   10.169)   16.065
  10.382   (  -0.995   -1.723   -3.391)    3.931
======================= Grid point 274 (48/120) =======================
q-point: ( 0.36  0.29  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.305   (  -4.440    7.691  -17.293)   19.440
   3.486   (   2.631   -4.557  -18.287)   19.029
   3.635   (  -3.775    6.539    3.324)    8.249
   4.018   (  -3.667    6.351    0.344)    7.342
   5.456   (  -2.408    4.171   10.925)   11.940
   5.607   (   3.327   -5.763    8.482)   10.781
   9.743   (   0.442   -0.765   -1.626)    1.851
   9.793   (   2.049   -3.548   -3.954)    5.694
   9.878   (  -0.045    0.078    0.484)    0.492
  10.000   (  -0.572    0.991   10.765)   10.825
  10.106   (   3.642   -6.308   10.844)   13.064
  10.365   (   0.753   -1.304   -3.274)    3.604
======================= Grid point 422 (49/120) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 120
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.000   (   0.000    0.000   23.212)   23.212
   1.000   (   0.000    0.000   23.212)   23.212
   1.934   (   0.000    0.000   -2.334)    2.334
   1.934   (   0.000    0.000   -2.334)    2.334
   2.348   (   0.000    0.000   58.433)   58.433
   5.696   (   0.000    0.000  -25.450)   25.450
  10.304   (  -0.000   -0.000    1.921)    1.921
  10.304   (   0.000    0.000    1.921)    1.921
  10.544   (   0.000    0.000   -2.977)    2.977
  10.544   (   0.000    0.000   -2.977)    2.977
  10.892   (   0.000    0.000    7.932)    7.932
  11.468   (   0.000    0.000   -4.396)    4.396
======================= Grid point 423 (50/120) =======================
q-point: ( 0.07  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.149   (  11.541    6.663   20.320)   24.300
   1.434   (  25.633   14.799   19.366)   35.371
   2.020   (   6.821    3.938   -2.116)    8.155
   2.258   (  23.512   13.575   -0.520)   27.155
   2.451   (  12.419    7.170   50.229)   52.236
   5.638   (  -4.674   -2.699  -25.043)   25.618
  10.252   (  -3.916   -2.261    2.168)    5.015
  10.287   (  -2.114   -1.220    1.322)    2.776
  10.519   (  -2.416   -1.395   -1.186)    3.032
  10.525   (  -1.826   -1.054   -3.375)    3.980
  10.825   (  -4.803   -2.773    9.172)   10.719
  11.438   (  -1.590   -0.918   -5.267)    5.578
======================= Grid point 424 (51/120) =======================
q-point: ( 0.14  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.486   (  15.476    8.935   15.853)   23.889
   1.873   (  11.289    6.518   26.927)   29.917
   2.231   (  10.257    5.922   -1.749)   11.972
   2.882   (  29.040   16.766    8.931)   34.701
   3.057   (  32.969   19.035   17.217)   41.782
   5.497   (  -6.245   -3.605  -23.183)   24.278
  10.139   (  -5.205   -3.005    2.351)    6.454
  10.163   (  -8.304   -4.794    0.265)    9.593
  10.412   (  -6.886   -3.976   -0.096)    7.952
  10.451   (  -4.369   -2.523   -3.963)    6.415
  10.738   (  -1.115   -0.644   11.475)   11.547
  11.420   (   0.025    0.015   -5.987)    5.987
======================= Grid point 425 (52/120) =======================
q-point: ( 0.21  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.837   (  13.220    7.632   12.739)   19.882
   2.092   (   8.641    4.989   20.775)   23.047
   2.474   (   9.485    5.476   -1.575)   11.065
   3.483   (  17.862   10.313   -5.325)   21.302
   3.864   (  33.701   19.458   17.012)   42.471
   5.406   (   0.275    0.159  -17.551)   17.554
   9.941   (  -9.050   -5.225    0.307)   10.454
  10.018   (  -4.840   -2.794    2.093)    5.968
  10.222   (  -8.105   -4.679   -0.588)    9.377
  10.331   (  -5.276   -3.046   -4.241)    7.423
  10.754   (   1.410    0.814   12.014)   12.124
  11.424   (  -0.078   -0.045   -4.365)    4.366
======================= Grid point 426 (53/120) =======================
q-point: ( 0.29  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.090   (   7.782    4.493   10.880)   14.111
   2.307   (   9.363    5.406    8.751)   13.909
   2.656   (   5.389    3.111   -1.527)    6.407
   3.670   (  -1.071   -0.619  -14.744)   14.796
   4.617   (  28.873   16.670   16.033)   36.995
   5.551   (  10.981    6.340   -8.242)   15.123
   9.777   (  -4.280   -2.471    0.375)    4.956
   9.914   (  -3.741   -2.160    1.738)    4.656
  10.056   (  -5.586   -3.225   -1.837)    6.706
  10.220   (  -3.649   -2.107   -4.439)    6.120
  10.780   (   0.652    0.377   13.257)   13.279
  11.403   (  -1.913   -1.105   -1.421)    2.627
======================= Grid point 427 (54/120) =======================
q-point: ( 0.36  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.209   (   2.541    1.467    9.826)   10.255
   2.530   (   9.191    5.306   -2.290)   10.857
   2.724   (   0.498    0.288   -1.299)    1.421
   3.538   (  -8.012   -4.626  -19.322)   21.422
   5.218   (  21.366   12.336   15.028)   28.888
   5.837   (  11.276    6.510   -1.612)   13.120
   9.734   (   0.192    0.111   -0.328)    0.396
   9.843   (  -2.192   -1.265    1.616)    3.002
   9.961   (  -2.434   -1.405   -2.748)    3.931
  10.163   (  -1.208   -0.697   -4.833)    5.031
  10.791   (   0.550    0.317   15.039)   15.052
  11.331   (  -3.976   -2.296    2.431)    5.195
======================= Grid point 428 (55/120) =======================
q-point: ( 0.43  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.233   (   0.055    0.032    9.220)    9.220
   2.702   (  -1.689   -0.975   -0.795)    2.106
   2.730   (   7.118    4.110  -10.108)   13.028
   3.349   (  -6.978   -4.029  -20.863)   22.365
   5.621   (  12.150    7.015   14.076)   19.874
   6.040   (   5.416    3.127    2.299)    6.663
   9.755   (   1.048    0.605   -1.130)    1.655
   9.809   (  -0.820   -0.473    1.652)    1.904
   9.933   (  -0.144   -0.083   -4.078)    4.082
  10.153   (   0.073    0.042   -5.282)    5.283
  10.813   (   1.167    0.674   16.614)   16.669
  11.229   (  -4.031   -2.327    6.495)    7.991
======================= Grid point 429 (56/120) =======================
q-point: (-0.50  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 235
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.231   (   0.000    0.000    8.982)    8.982
   2.676   (   0.000    0.000   -0.482)    0.482
   2.829   (   0.000    0.000  -13.106)   13.106
   3.252   (  -0.000   -0.000  -21.095)   21.095
   5.772   (   0.000    0.000   13.514)   13.514
   6.100   (  -0.000   -0.000    3.858)    3.858
   9.769   (   0.000    0.000   -1.294)    1.294
   9.800   (   0.000    0.000    1.684)    1.684
   9.935   (   0.000    0.000   -5.101)    5.101
  10.156   (   0.000    0.000   -5.479)    5.479
  10.832   (   0.000    0.000   17.200)   17.200
  11.173   (   0.000    0.000    8.475)    8.475
======================= Grid point 438 (57/120) =======================
q-point: ( 0.07  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.440   (  11.005   19.061   16.394)   27.444
   1.766   (   7.795   13.502   24.739)   29.242
   2.195   (   7.108   12.312   -0.482)   14.224
   2.726   (  15.951   27.627   -3.196)   32.061
   2.797   (  16.100   27.887   35.047)   47.594
   5.538   (  -3.717   -6.438  -23.934)   25.062
  10.177   (  -3.008   -5.210    2.138)    6.385
  10.196   (  -4.288   -7.427    0.424)    8.587
  10.435   (  -3.164   -5.480    0.026)    6.328
  10.523   (  -0.372   -0.644   -4.049)    4.116
  10.744   (  -2.219   -3.843   11.075)   11.931
  11.408   (  -0.996   -1.725   -5.830)    6.160
======================= Grid point 439 (58/120) =======================
q-point: ( 0.14  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.837   (   7.225   22.611   13.269)   27.194
   1.999   (   8.839    6.630   23.053)   25.564
   2.479   (   6.293   17.039   -1.191)   18.203
   3.271   (  17.734   16.974   -3.674)   24.821
   3.513   (  28.640   25.461   20.370)   43.399
   5.414   (  -2.618   -3.573  -20.346)   20.823
  10.007   (  -6.405   -8.436    0.145)   10.593
  10.055   (  -3.878   -5.380    1.707)    6.849
  10.301   (  -6.185   -4.394   -0.362)    7.595
  10.462   (  -6.492    0.760   -3.695)    7.508
  10.705   (   1.501   -2.604   11.656)   12.037
  11.375   (   0.888   -4.426   -4.790)    6.582
======================= Grid point 440 (59/120) =======================
q-point: ( 0.21  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.116   (   4.355   13.227   10.972)   17.728
   2.259   (   5.949   15.561   12.790)   21.003
   2.732   (   0.800   18.385   -1.468)   18.461
   3.608   (   7.733    1.357  -11.303)   13.763
   4.270   (  29.138   19.387   17.147)   38.973
   5.437   (   6.526    3.026  -12.423)   14.356
   9.832   (  -5.112   -4.182    0.459)    6.621
   9.948   (  -2.714   -3.604    1.145)    4.654
  10.137   (  -7.491   -2.750   -0.946)    8.036
  10.336   (  -8.175    1.704   -3.845)    9.194
  10.714   (   3.133   -4.739   12.160)   13.422
  11.342   (   2.169   -7.849   -2.435)    8.500
======================= Grid point 441 (60/120) =======================
q-point: ( 0.29  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.300   (   1.675   11.712    5.127)   12.894
   2.488   (   2.387   14.740    5.851)   16.038
   2.871   (  -5.428   16.628   -1.603)   17.565
   3.602   (  -3.970   -4.921  -17.456)   18.566
   4.927   (  25.160   13.290   15.650)   32.474
   5.673   (  12.863    5.034   -4.304)   14.467
   9.750   (  -1.833    0.006    0.310)    1.859
   9.884   (  -2.282   -0.677    1.171)    2.653
   9.996   (  -4.840   -2.454   -1.756)    5.704
  10.228   (  -5.500    1.772   -4.458)    7.298
  10.715   (   3.685   -6.513   13.690)   15.602
  11.282   (   1.469  -10.334    0.782)   10.467
======================= Grid point 442 (61/120) =======================
q-point: ( 0.36  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.400   (  -5.602   14.531    5.208)   16.421
   2.671   (   4.389    9.443   -4.227)   11.239
   2.892   (  -9.666   14.782   -1.640)   17.737
   3.439   (  -7.576   -4.250  -19.864)   21.681
   5.425   (  18.003    7.835   14.555)   24.440
   5.925   (  10.014    2.008    0.773)   10.243
   9.749   (   0.420    1.144   -0.585)    1.352
   9.841   (  -2.551    0.824    1.613)    3.128
   9.926   (  -0.561   -1.886   -3.010)    3.597
  10.179   (  -2.199    1.649   -4.938)    5.652
  10.720   (   4.851   -7.066   15.480)   17.694
  11.185   (   0.644  -11.402    4.688)   12.345
======================= Grid point 443 (62/120) =======================
q-point: ( 0.43  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.421   (  -9.041   16.042    6.249)   19.446
   2.820   (   2.134    6.220  -10.432)   12.331
   2.862   (  -9.592   14.079   -2.224)   17.181
   3.286   (  -4.878   -1.525  -20.300)   20.933
   5.711   (   7.681    1.808   13.585)   15.711
   6.052   (   4.192   -1.856    3.672)    5.874
   9.766   (   0.305    0.360   -1.057)    1.158
   9.819   (  -1.189    1.332    1.607)    2.402
   9.916   (   0.916   -1.655   -4.424)    4.811
  10.171   (  -0.849    1.433   -5.268)    5.525
  10.740   (   5.361   -7.515   16.726)   19.105
  11.094   (   2.605  -10.241    7.911)   13.200
======================= Grid point 453 (63/120) =======================
q-point: ( 0.14  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.127   (   4.448    7.704   14.057)   16.636
   2.341   (  12.827   22.217   10.038)   27.548
   2.860   (  10.424   18.055   -0.866)   20.866
   3.514   (   3.807    6.593  -10.704)   13.135
   4.080   (  16.405   28.415   18.900)   37.865
   5.394   (   1.443    2.499  -14.383)   14.670
   9.858   (  -3.133   -5.426    0.464)    6.283
   9.948   (  -2.538   -4.396    0.533)    5.104
  10.215   (  -2.739   -4.745   -0.687)    5.522
  10.443   (  -1.628   -2.819   -2.379)    4.032
  10.642   (  -1.891   -3.276   10.578)   11.234
  11.255   (  -3.811   -6.601   -2.339)    7.973
======================= Grid point 454 (64/120) =======================
q-point: ( 0.21  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.296   (   5.781    7.483    3.143)    9.965
   2.718   (   2.319   23.411    8.312)   24.951
   3.158   (  -0.866   21.020   -1.341)   21.081
   3.565   (   1.745   -4.381  -15.752)   16.443
   4.656   (  20.622   18.180   16.792)   32.214
   5.527   (   8.173    5.462   -6.860)   11.987
   9.786   (  -1.876   -0.365    0.832)    2.084
   9.887   (  -0.504   -2.181    0.059)    2.239
  10.093   (  -6.694   -2.102   -0.682)    7.050
  10.356   (  -5.502   -0.575   -2.538)    6.086
  10.578   (   1.136   -7.512   10.877)   13.268
  11.124   (   0.256  -12.426    0.489)   12.438
======================= Grid point 455 (65/120) =======================
q-point: ( 0.29  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 788
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.481   (   6.000    7.366   -5.793)   11.127
   2.911   (  -6.259   22.355    6.127)   24.010
   3.254   (  -9.495   16.779   -5.186)   19.965
   3.519   (  -3.623   -0.458  -14.952)   15.391
   5.169   (  19.688   10.488   15.271)   27.034
   5.751   (  10.828    2.652   -0.884)   11.184
   9.769   (  -1.053    1.440    0.414)    1.832
   9.876   (   0.336   -0.623    0.037)    0.709
   9.956   (  -5.860   -1.321   -0.491)    6.027
  10.261   (  -4.657    0.777   -3.844)    6.088
  10.534   (   4.254  -10.036   12.921)   16.905
  11.004   (   2.428  -15.649    3.580)   16.236
======================= Grid point 456 (66/120) =======================
q-point: ( 0.36  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.643   (   2.468    8.147  -10.052)   13.173
   2.994   (  -7.795   20.018    2.152)   21.589
   3.204   (  -9.989    9.011  -11.210)   17.511
   3.439   (  -9.988    9.647   -9.342)   16.736
   5.548   (  13.386    4.386   13.879)   19.776
   5.924   (   7.333   -1.931    3.144)    8.209
   9.770   (  -0.247    0.606   -0.532)    0.843
   9.859   (  -1.475   -0.098    0.498)    1.560
   9.887   (  -0.669   -1.270   -1.699)    2.224
  10.213   (  -2.108    1.405   -4.529)    5.189
  10.524   (   6.613  -11.009   15.026)   19.766
  10.893   (   3.829  -15.962    6.766)   17.754
======================= Grid point 457 (67/120) =======================
q-point: ( 0.43  0.14  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.711   (  -4.502    7.798   -9.641)   13.191
   3.053   ( -10.723   18.574    1.123)   21.476
   3.129   (  -3.051    5.284  -16.804)   17.877
   3.394   (  -9.378   16.244   -6.053)   19.709
   5.700   (   1.543   -2.672   13.083)   13.442
   5.974   (   3.037   -5.260    4.845)    7.770
   9.770   (   0.027   -0.047   -0.721)    0.723
   9.855   (  -1.106    1.916    1.624)    2.745
   9.868   (   1.613   -2.794   -3.752)    4.948
  10.204   (  -0.947    1.641   -4.857)    5.213
  10.530   (   6.770  -11.725   16.035)   20.986
  10.838   (   7.885  -13.657    8.297)   17.819
======================= Grid point 469 (68/120) =======================
q-point: ( 0.21  0.21  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.459   (   4.821    8.351   -5.831)   11.269
   3.110   (   7.526   13.035    3.728)   15.507
   3.467   (  -1.810   -3.136  -17.063)   17.443
   3.536   (   7.820   13.544   -0.483)   15.647
   5.026   (  10.257   17.765   15.911)   25.961
   5.655   (   3.933    6.812   -1.842)    8.079
   9.796   (   0.325    0.564    1.097)    1.276
   9.857   (  -0.453   -0.784   -0.328)    0.963
  10.031   (  -2.440   -4.226   -0.815)    4.948
  10.295   (  -2.527   -4.377    3.699)    6.263
  10.450   (  -2.683   -4.647    4.878)    7.252
  10.871   (  -6.447  -11.167    3.883)   13.466
======================= Grid point 470 (69/120) =======================
q-point: ( 0.29  0.21  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.637   (   5.399    7.923  -12.897)   16.071
   3.220   (  -2.860    2.435  -13.765)   14.269
   3.446   (  -4.812   10.669   -2.087)   11.889
   3.728   (  -2.590   14.461   -0.888)   14.718
   5.374   (  12.151    9.683   14.423)   21.200
   5.789   (   5.702    1.180    2.524)    6.346
   9.788   (  -1.267   -0.170    0.544)    1.390
   9.844   (  -0.040   -2.271   -1.147)    2.545
   9.936   (  -3.757   -1.146    0.551)    3.966
  10.235   (   0.615   -3.315    4.040)    5.262
  10.367   (  -1.227   -4.788    4.910)    6.967
  10.681   (   0.051  -14.678    6.680)   16.126
======================= Grid point 471 (70/120) =======================
q-point: ( 0.36  0.21  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.790   (   1.990    6.736  -16.935)   18.334
   3.156   (  -1.906   -2.786  -22.125)   22.381
   3.464   (  -9.672   17.369    5.006)   20.501
   3.775   (  -8.955   16.536   -0.721)   18.819
   5.620   (   6.656    2.900   12.776)   14.695
   5.859   (   3.624   -3.986    5.657)    7.812
   9.765   (  -0.180   -1.212   -0.086)    1.229
   9.811   (   0.607   -3.280   -2.211)    4.002
   9.892   (  -1.322    0.491    1.503)    2.061
  10.232   (   0.053    0.391   -2.231)    2.266
  10.319   (   3.090   -7.515   12.397)   14.823
  10.573   (   4.729  -14.094    8.038)   16.900
======================= Grid point 484 (71/120) =======================
q-point: ( 0.29  0.29  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 429
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.805   (   4.671    8.091  -18.129)   20.395
   3.159   (  -2.977   -5.156  -22.811)   23.575
   3.655   (   4.047    7.009    5.234)    9.638
   3.953   (   3.961    6.861   -0.347)    7.930
   5.566   (   4.892    8.473   12.739)   16.063
   5.799   (  -0.292   -0.506    5.945)    5.974
   9.756   (  -1.222   -2.116    0.849)    2.587
   9.790   (  -1.718   -2.975   -1.728)    3.846
   9.907   (  -0.652   -1.129    1.428)    1.934
  10.194   (   0.195    0.337    9.515)    9.523
  10.306   (  -1.345   -2.330   -2.061)    3.389
  10.433   (  -5.093   -8.821    9.690)   14.058
======================= Grid point 485 (72/120) =======================
q-point: ( 0.36  0.29  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 427
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.922   (  -3.232    5.599  -20.646)   21.634
   3.073   (   2.535   -4.392  -23.002)   23.554
   3.725   (  -3.913    6.777    5.736)    9.703
   4.020   (  -3.670    6.356   -0.312)    7.346
   5.679   (  -1.597    2.766   11.108)   11.558
   5.776   (   2.077   -3.597    8.168)    9.163
   9.734   (   0.644   -1.116    0.374)    1.341
   9.749   (   1.241   -2.150   -0.917)    2.647
   9.896   (   0.000   -0.001    1.557)    1.557
  10.230   (   0.088   -0.152   12.521)   12.523
  10.252   (  -0.433    0.749   -3.019)    3.141
  10.342   (   4.529   -7.844    7.913)   12.027
======================= Grid point 633 (73/120) =======================
q-point: ( 0.00  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 96
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.394   (   0.000    0.000   18.037)   18.037
   1.394   (   0.000    0.000   18.037)   18.037
   1.854   (   0.000    0.000   -6.417)    6.417
   1.854   (  -0.000    0.000   -6.417)    6.417
   3.406   (   0.000    0.000   52.605)   52.605
   5.109   (   0.000    0.000  -35.984)   35.984
  10.349   (   0.000    0.000    2.873)    2.873
  10.349   (   0.000    0.000    2.873)    2.873
  10.480   (   0.000    0.000   -3.622)    3.622
  10.480   (   0.000    0.000   -3.622)    3.622
  11.056   (   0.000    0.000    8.944)    8.944
  11.363   (   0.000   -0.000   -6.536)    6.536
======================= Grid point 634 (74/120) =======================
q-point: ( 0.07  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.507   (   8.979    5.184   16.851)   19.785
   1.768   (  25.635   14.800   15.651)   33.484
   1.944   (   7.187    4.150   -6.111)   10.307
   2.184   (  24.486   14.137   -6.118)   28.928
   3.427   (   2.909    1.680   49.826)   49.939
   5.064   (  -3.470   -2.004  -35.004)   35.233
  10.304   (  -3.548   -2.048    3.254)    5.232
  10.326   (  -2.642   -1.525    2.852)    4.176
  10.453   (  -2.406   -1.389   -4.110)    4.960
  10.469   (  -1.633   -0.943   -3.675)    4.131
  11.008   (  -3.183   -1.837    9.679)   10.353
  11.325   (  -2.342   -1.352   -6.728)    7.251
======================= Grid point 635 (75/120) =======================
q-point: ( 0.14  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.780   (  12.957    7.481   14.439)   20.793
   2.166   (  10.642    6.144   -5.435)   13.437
   2.324   (  16.711    9.648   20.117)   27.875
   2.832   (  26.395   15.239   -7.208)   31.319
   3.648   (  17.790   10.271   36.242)   41.659
   4.980   (  -2.263   -1.307  -30.784)   30.895
  10.192   (  -8.502   -4.909    2.929)   10.245
  10.196   (  -5.195   -3.000    3.632)    7.013
  10.366   (  -4.735   -2.734   -4.775)    7.259
  10.367   (  -7.151   -4.128   -4.304)    9.311
  10.961   (  -0.059   -0.034   11.259)   11.259
  11.298   (   0.457    0.264   -6.855)    6.875
======================= Grid point 636 (76/120) =======================
q-point: ( 0.21  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.080   (  11.495    6.637   12.260)   18.069
   2.416   (   9.690    5.595   -4.825)   12.185
   2.532   (   3.012    1.739   24.127)   24.377
   3.296   (  10.419    6.015  -13.779)   18.292
   4.240   (  29.057   16.776   21.181)   39.678
   5.042   (   8.979    5.184  -20.254)   22.753
   9.968   (  -9.118   -5.264    2.884)   10.916
  10.070   (  -5.144   -2.970    3.464)    6.876
  10.163   (  -8.920   -5.150   -5.194)   11.535
  10.241   (  -5.375   -3.103   -4.990)    7.964
  10.994   (   2.146    1.239   12.228)   12.476
  11.330   (   1.681    0.970   -5.682)    6.004
======================= Grid point 637 (77/120) =======================
q-point: ( 0.29  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.301   (   6.745    3.894   10.781)   13.301
   2.559   (   0.328    0.189   17.302)   17.306
   2.601   (   5.556    3.208   -4.429)    7.795
   3.333   (  -5.195   -2.999  -20.034)   20.913
   4.923   (  27.000   15.588   15.427)   34.785
   5.387   (  17.794   10.274   -9.252)   22.534
   9.798   (  -4.769   -2.754    2.307)    5.971
   9.961   (  -3.918   -2.262    3.172)    5.526
   9.980   (  -5.999   -3.464   -5.562)    8.884
  10.128   (  -3.863   -2.230   -5.092)    6.769
  11.041   (   1.708    0.986   13.109)   13.257
  11.359   (   0.542    0.313   -3.588)    3.642
======================= Grid point 638 (78/120) =======================
q-point: ( 0.36  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.401   (   1.911    1.104    9.850)   10.094
   2.573   (   1.025    0.592    7.168)    7.265
   2.676   (   0.870    0.502   -4.004)    4.128
   3.150   (  -8.680   -5.011  -20.383)   22.713
   5.489   (  20.105   11.608   12.760)   26.491
   5.792   (  14.962    8.638   -3.507)   17.629
   9.739   (  -0.598   -0.345    1.247)    1.425
   9.881   (  -2.548   -1.471   -5.109)    5.895
   9.888   (  -2.144   -1.238    3.179)    4.030
  10.063   (  -1.604   -0.926   -5.468)    5.773
  11.077   (   1.398    0.807   13.708)   13.803
  11.353   (  -0.954   -0.551   -0.756)    1.336
======================= Grid point 639 (79/120) =======================
q-point: ( 0.43  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.413   (  -0.333   -0.192    9.201)    9.209
   2.610   (   1.891    1.092   -1.571)    2.690
   2.664   (  -1.234   -0.712   -3.401)    3.687
   2.963   (  -6.433   -3.714  -18.443)   19.883
   5.863   (  11.076    6.395   10.769)   16.720
   6.064   (   7.569    4.370   -0.205)    8.742
   9.745   (   0.672    0.388    0.333)    0.845
   9.847   (  -0.571   -0.330   -4.553)    4.601
   9.856   (  -0.732   -0.423    3.363)    3.468
  10.044   (  -0.242   -0.140   -5.945)    5.952
  11.111   (   1.370    0.791   13.518)   13.610
  11.322   (  -1.441   -0.832    2.406)    2.925
======================= Grid point 640 (80/120) =======================
q-point: (-0.50  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 228
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.406   (   0.000    0.000    8.915)    8.915
   2.641   (   0.000    0.000   -5.483)    5.483
   2.645   (   0.000    0.000   -3.056)    3.056
   2.877   (  -0.000   -0.000  -16.962)   16.962
   5.999   (   0.000    0.000    9.752)    9.752
   6.152   (   0.000    0.000    1.184)    1.184
   9.755   (   0.000    0.000    0.050)    0.050
   9.842   (   0.000    0.000   -4.414)    4.414
   9.848   (   0.000    0.000    3.460)    3.460
  10.042   (   0.000    0.000   -6.160)    6.160
  11.130   (   0.000    0.000   13.097)   13.097
  11.301   (   0.000    0.000    4.040)    4.040
======================= Grid point 649 (81/120) =======================
q-point: ( 0.07  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.742   (   9.178   15.897   14.714)   23.525
   2.101   (   8.340   14.446    1.692)   16.766
   2.207   (  10.584   18.332   10.753)   23.743
   2.633   (  15.102   26.158   -6.809)   30.962
   3.546   (   6.918   11.982   41.466)   43.714
   4.998   (  -2.045   -3.541  -32.448)   32.705
  10.221   (  -4.681   -8.107    2.212)    9.619
  10.229   (  -3.234   -5.602    3.081)    7.164
  10.403   (  -2.791   -4.834   -2.947)    6.312
  10.438   (  -0.770   -1.333   -4.612)    4.862
  10.957   (  -1.270   -2.199   10.817)   11.111
  11.288   (  -0.913   -1.582   -6.834)    7.074
======================= Grid point 650 (82/120) =======================
q-point: ( 0.14  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.087   (   6.665   20.624   12.455)   24.998
   2.406   (   7.945   14.078    2.201)   16.314
   2.476   (   5.127    6.639   16.623)   18.620
   3.127   (  15.119   11.235  -10.845)   21.736
   3.970   (  20.655   20.361   26.111)   39.026
   4.978   (   2.755    2.106  -25.038)   25.277
  10.026   (  -6.171   -8.867    2.024)   10.991
  10.090   (  -4.038   -7.316    1.358)    8.466
  10.273   (  -6.937   -2.963   -1.699)    7.732
  10.371   (  -6.243    0.384   -5.214)    8.144
  10.939   (   1.774   -2.147   11.900)   12.221
  11.274   (   1.760   -2.778   -6.015)    6.855
======================= Grid point 651 (83/120) =======================
q-point: ( 0.21  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.351   (   0.745   15.475   12.277)   19.768
   2.557   (   2.560    3.882   18.118)   18.706
   2.679   (   0.670   18.559   -3.821)   18.960
   3.328   (   2.679   -2.648  -17.801)   18.195
   4.606   (  26.075   17.910   17.484)   36.144
   5.190   (  14.068    8.799  -13.488)   21.384
   9.852   (  -4.920   -4.567    2.000)    7.005
   9.960   (  -2.484   -5.751   -0.703)    6.304
  10.110   (  -8.648   -0.699   -0.564)    8.694
  10.244   (  -8.276    1.416   -5.403)    9.984
  10.961   (   4.023   -4.505   12.589)   13.963
  11.281   (   3.657   -5.764   -4.299)    8.067
======================= Grid point 652 (84/120) =======================
q-point: ( 0.29  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.442   (   0.095    3.845   10.181)   10.883
   2.648   (  -3.356   13.236    9.620)   16.703
   2.816   (  -5.231   15.703   -4.084)   17.047
   3.241   (  -6.087   -4.333  -19.534)   20.914
   5.215   (  23.602   12.452   13.862)   30.071
   5.581   (  17.520    8.141   -5.754)   20.157
   9.764   (  -2.341   -0.332    1.549)    2.827
   9.887   (  -1.124   -2.548   -2.393)    3.672
   9.967   (  -6.577    0.078    0.932)    6.643
  10.130   (  -5.689    1.437   -5.556)    8.081
  10.984   (   4.921   -6.464   13.270)   15.559
  11.276   (   3.724   -8.204   -1.979)    9.224
======================= Grid point 653 (85/120) =======================
q-point: ( 0.36  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.475   (   0.228    2.581    4.312)    5.030
   2.694   (  -6.304   13.126    4.545)   15.254
   2.828   (  -8.377   11.198   -5.262)   14.942
   3.079   (  -9.113    0.904  -16.177)   18.589
   5.680   (  16.810    6.964   11.576)   21.565
   5.922   (  12.762    3.942   -1.542)   13.445
   9.747   (  -0.319    0.951    0.690)    1.217
   9.852   (  -0.994   -1.203   -4.086)    4.374
   9.889   (  -2.769    0.890    3.188)    4.316
  10.074   (  -2.496    1.468   -5.840)    6.519
  11.007   (   5.747   -7.339   13.461)   16.374
  11.247   (   3.372   -9.493    0.975)   10.121
======================= Grid point 654 (86/120) =======================
q-point: ( 0.43  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.500   (  -1.710    4.392    1.771)    5.035
   2.718   (  -4.814   11.684   -0.412)   12.644
   2.775   (  -7.259    6.846   -7.255)   12.337
   2.962   (  -8.367    7.934  -11.297)   16.143
   5.940   (   7.108    1.097    9.842)   12.190
   6.100   (   5.258   -0.798    0.950)    5.402
   9.756   (   0.105    0.496    0.130)    0.523
   9.834   (   0.215   -0.929   -3.718)    3.838
   9.867   (  -1.320    1.279    3.260)    3.742
  10.060   (  -1.068    1.436   -6.127)    6.383
  11.033   (   5.694   -8.017   13.007)   16.306
  11.212   (   4.099   -9.414    3.576)   10.873
======================= Grid point 664 (87/120) =======================
q-point: ( 0.14  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.469   (  -0.138   -0.238   20.687)   20.689
   2.538   (  11.670   20.214   10.073)   25.422
   2.820   (  10.612   18.381   -3.541)   21.517
   3.251   (   0.626    1.085  -16.828)   16.875
   4.446   (  14.359   24.870   19.029)   34.449
   5.104   (   6.034   10.452  -15.828)   19.905
   9.873   (  -3.167   -5.486    1.405)    6.489
   9.950   (  -3.158   -5.470   -0.805)    6.367
  10.205   (  -2.467   -4.272    0.413)    4.950
  10.351   (  -1.581   -2.739   -5.767)    6.578
  10.880   (  -1.831   -3.172   12.415)   12.944
  11.195   (  -2.495   -4.321   -4.252)    6.556
======================= Grid point 665 (88/120) =======================
q-point: ( 0.21  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.452   (  -0.246   -0.603   12.875)   12.892
   2.872   (   1.387   18.519    6.473)   19.667
   3.093   (  -3.969   17.217   -6.546)   18.843
   3.250   (   3.282   -0.203  -14.702)   15.065
   4.964   (  18.821   16.844   14.980)   29.366
   5.388   (  12.696   10.208   -7.920)   18.114
   9.801   (  -1.841   -0.585    1.054)    2.201
   9.875   (  -0.750   -2.566   -1.573)    3.102
  10.093   (  -6.405   -1.422    1.391)    6.707
  10.262   (  -6.098   -0.131   -5.974)    8.537
  10.824   (   2.184   -8.048   12.745)   15.230
  11.112   (   1.519   -9.815   -2.233)   10.180
======================= Grid point 666 (89/120) =======================
q-point: ( 0.29  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.460   (   1.530    0.143    4.188)    4.461
   2.924   (  -5.580    2.257  -11.671)   13.131
   3.135   (  -8.341   17.989    1.310)   19.871
   3.335   (  -6.800   16.560   -5.641)   18.769
   5.438   (  18.252    9.443   12.321)   23.961
   5.719   (  13.987    5.403   -2.949)   15.281
   9.777   (  -1.522    1.219    0.782)    2.101
   9.852   (  -0.320   -1.268   -2.404)    2.736
   9.973   (  -5.538    0.309    2.681)    6.161
  10.160   (  -5.171    1.130   -6.003)    8.003
  10.796   (   5.660  -10.864   12.900)   17.789
  11.046   (   4.234  -13.241    0.068)   13.902
======================= Grid point 667 (90/120) =======================
q-point: ( 0.36  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.504   (   2.230    1.114   -2.871)    3.802
   2.828   (  -3.685   -2.419  -16.982)   17.544
   3.147   ( -13.094   21.772    5.462)   25.986
   3.340   ( -12.947   20.239   -3.823)   24.328
   5.782   (  12.294    3.204   10.126)   16.246
   5.960   (   9.240    0.010    0.308)    9.245
   9.767   (  -0.750    0.659    0.298)    1.042
   9.824   (  -0.720   -1.661   -2.527)    3.108
   9.916   (  -2.386    1.416    3.307)    4.317
  10.108   (  -2.491    1.572   -6.068)    6.745
  10.799   (   7.603  -11.933   12.593)   18.942
  10.991   (   5.710  -14.386    2.663)   15.705
======================= Grid point 668 (91/120) =======================
q-point: ( 0.43  0.14  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.537   (  -0.206    0.356   -5.869)    5.884
   2.775   (   0.663   -1.148  -17.231)   17.282
   3.138   ( -13.333   23.094    5.932)   27.318
   3.324   ( -12.801   22.172   -3.276)   25.810
   5.915   (   1.896   -3.283    8.938)    9.709
   6.041   (   2.613   -4.526    1.804)    5.529
   9.763   (  -0.135    0.234    0.110)    0.292
   9.809   (   0.799   -1.384   -2.421)    2.901
   9.902   (  -1.101    1.906    3.362)    4.019
  10.095   (  -1.077    1.865   -6.147)    6.513
  10.809   (   7.338  -12.709   12.209)   19.090
  10.964   (   7.865  -13.622    4.032)   16.238
======================= Grid point 680 (92/120) =======================
q-point: ( 0.21  0.21  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.445   (   0.293    0.508    4.954)    4.988
   2.976   (  -3.056   -5.293  -17.785)   18.806
   3.336   (   7.353   12.735    3.925)   15.220
   3.508   (   7.929   13.733   -2.611)   16.071
   5.304   (   9.324   16.150   12.758)   22.595
   5.606   (   6.194   10.728   -3.720)   12.934
   9.806   (   0.213    0.369    0.349)    0.551
   9.842   (  -0.540   -0.935   -1.188)    1.606
  10.042   (  -1.990   -3.447    2.292)    4.593
  10.230   (  -1.629   -2.821   -6.105)    6.920
  10.668   (  -3.609   -6.251   12.536)   14.466
  10.915   (  -4.932   -8.542    0.046)    9.863
======================= Grid point 681 (93/120) =======================
q-point: ( 0.29  0.21  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.478   (   1.936    1.840   -2.639)    3.754
   2.853   (  -3.642   -5.153  -19.592)   20.583
   3.523   (  -3.463   16.359    5.955)   17.750
   3.698   (  -3.172   15.512   -2.523)   16.033
   5.616   (  10.872    8.105   10.481)   17.139
   5.813   (   7.786    3.931   -0.371)    8.730
   9.794   (  -1.591   -0.051    0.262)    1.613
   9.820   (  -0.771   -2.023   -1.190)    2.470
   9.972   (  -2.779   -0.599    3.217)    4.293
  10.168   (  -2.769   -0.350   -6.080)    6.690
  10.582   (   2.590   -8.821   11.880)   15.022
  10.774   (   1.013  -12.118    2.308)   12.377
======================= Grid point 682 (94/120) =======================
q-point: ( 0.36  0.21  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 784
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.537   (   1.455    2.263   -8.155)    8.587
   2.755   (  -1.492   -3.770  -18.218)   18.664
   3.576   (  -9.457   18.007    6.067)   21.225
   3.747   (  -9.232   17.238   -2.398)   19.701
   5.825   (   6.039    1.335    8.386)   10.420
   5.936   (   4.461   -1.971    2.157)    5.332
   9.767   (  -0.564   -0.696    0.257)    0.932
   9.783   (  -0.183   -2.131   -0.851)    2.302
   9.938   (  -1.242    0.625    3.364)    3.640
  10.133   (  -1.470    0.811   -6.161)    6.385
  10.564   (   5.519   -9.585   11.096)   15.667
  10.699   (   5.097  -12.829    4.192)   14.427
======================= Grid point 695 (95/120) =======================
q-point: ( 0.29  0.29  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.536   (   2.030    3.516   -8.764)    9.659
   2.750   (  -2.691   -4.661  -18.508)   19.274
   3.767   (   4.233    7.332    5.950)   10.348
   3.931   (   4.031    6.982   -2.160)    8.346
   5.770   (   3.690    6.391    8.252)   11.071
   5.879   (   1.227    2.125    2.178)    3.281
   9.770   (  -1.088   -1.885    0.485)    2.230
   9.773   (  -1.268   -2.197   -0.199)    2.545
   9.953   (  -0.532   -0.922    3.451)    3.611
  10.153   (  -0.494   -0.856   -6.267)    6.344
  10.460   (  -1.676   -2.902   10.087)   10.629
  10.575   (  -3.992   -6.914    4.620)    9.224
======================= Grid point 696 (96/120) =======================
q-point: ( 0.36  0.29  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 420
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.593   (  -1.603    2.777  -12.418)   12.825
   2.682   (   1.777   -3.078  -16.527)   16.905
   3.839   (  -3.894    6.745    5.909)    9.776
   3.999   (  -3.743    6.484   -2.066)    7.767
   5.849   (  -0.637    1.103    6.508)    6.632
   5.894   (   1.040   -1.801    3.989)    4.499
   9.746   (   0.530   -0.919    0.640)    1.239
   9.746   (   0.640   -1.109    0.351)    1.328
   9.942   (   0.008   -0.014    3.428)    3.428
  10.143   (  -0.125    0.216   -6.445)    6.449
  10.441   (   1.221   -2.115    8.849)    9.180
  10.497   (   3.507   -6.074    6.000)    9.230
======================= Grid point 844 (97/120) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.680   (  -0.000    0.000  -12.042)   12.042
   1.680   (  -0.000    0.000  -12.042)   12.042
   1.680   (   0.000    0.000   12.042)   12.042
   1.680   (  -0.000    0.000   12.042)   12.042
   4.337   (  -0.000   -0.000  -45.037)   45.037
   4.337   (   0.000    0.000   45.037)   45.037
  10.411   (   0.000    0.000   -3.537)    3.537
  10.411   (  -0.000   -0.000   -3.537)    3.537
  10.411   (   0.000   -0.000    3.537)    3.537
  10.411   (   0.000    0.000    3.537)    3.537
  11.222   (  -0.000    0.000   -8.248)    8.248
  11.222   (  -0.000   -0.000    8.248)    8.248
======================= Grid point 846 (98/120) =======================
q-point: ( 0.07  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.778   (   7.791    4.498  -11.499)   14.600
   1.778   (   7.791    4.498   11.499)   14.600
   2.022   (  24.813   14.326  -10.987)   30.685
   2.022   (  24.813   14.326   10.987)   30.685
   4.317   (  -1.134   -0.655  -43.370)   43.389
   4.317   (  -1.134   -0.655   43.370)   43.389
  10.374   (  -3.021   -1.744   -4.013)    5.317
  10.374   (  -3.021   -1.744    4.013)    5.317
  10.393   (  -2.222   -1.283   -4.023)    4.772
  10.393   (  -2.222   -1.283    4.023)    4.772
  11.181   (  -2.528   -1.460   -8.391)    8.885
  11.181   (  -2.528   -1.460    8.391)    8.885
======================= Grid point 848 (99/120) =======================
q-point: ( 0.14  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.017   (  11.455    6.614  -10.235)   16.724
   2.017   (  11.455    6.614   10.235)   16.724
   2.639   (  22.752   13.136  -13.191)   29.397
   2.639   (  22.752   13.136   13.191)   29.397
   4.343   (   5.314    3.068  -35.776)   36.298
   4.343   (   5.314    3.068   35.776)   36.298
  10.271   (  -8.022   -4.632   -5.119)   10.584
  10.271   (  -8.022   -4.632    5.119)   10.584
  10.275   (  -5.039   -2.909   -4.586)    7.409
  10.275   (  -5.039   -2.909    4.586)    7.409
  11.152   (   0.594    0.343   -8.768)    8.795
  11.152   (   0.594    0.343    8.768)    8.795
======================= Grid point 850 (100/120) =======================
q-point: ( 0.21  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.285   (  10.319    5.958   -8.951)   14.903
   2.285   (  10.319    5.958    8.951)   14.903
   2.966   (   4.393    2.536  -20.664)   21.277
   2.966   (   4.393    2.536   20.664)   21.277
   4.646   (  19.829   11.448  -21.245)   31.235
   4.646   (  19.829   11.448   21.245)   31.235
  10.052   (  -9.245   -5.338   -5.683)   12.094
  10.052   (  -9.245   -5.338    5.683)   12.094
  10.148   (  -5.344   -3.085   -4.616)    7.706
  10.148   (  -5.344   -3.085    4.616)    7.706
  11.196   (   2.460    1.420   -8.783)    9.231
  11.196   (   2.460    1.420    8.783)    9.231
======================= Grid point 852 (101/120) =======================
q-point: ( 0.29  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.483   (   5.987    3.457   -8.042)   10.605
   2.483   (   5.987    3.457    8.042)   10.605
   2.934   (  -4.975   -2.872  -21.193)   21.957
   2.934   (  -4.975   -2.872   21.193)   21.957
   5.186   (  23.317   13.462  -12.123)   29.527
   5.186   (  23.317   13.462   12.123)   29.527
   9.871   (  -5.539   -3.198   -5.243)    8.270
   9.871   (  -5.539   -3.198    5.243)    8.270
  10.035   (  -3.972   -2.293   -4.518)    6.438
  10.035   (  -3.972   -2.293    4.518)    6.438
  11.249   (   1.806    1.043   -8.352)    8.609
  11.249   (   1.806    1.043    8.352)    8.609
======================= Grid point 854 (102/120) =======================
q-point: ( 0.36  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.567   (   1.343    0.776   -7.368)    7.529
   2.567   (   1.343    0.776    7.368)    7.529
   2.798   (  -5.659   -3.267  -15.908)   17.198
   2.798   (  -5.659   -3.267   15.908)   17.198
   5.687   (  17.918   10.345   -7.900)   22.147
   5.687   (  17.918   10.345    7.900)   22.147
   9.789   (  -1.637   -0.945   -4.014)    4.437
   9.789   (  -1.637   -0.945    4.014)    4.437
   9.964   (  -1.947   -1.124   -4.731)    5.238
   9.964   (  -1.947   -1.124    4.731)    5.238
  11.280   (   0.910    0.526   -7.286)    7.361
  11.280   (   0.910    0.526    7.286)    7.361
======================= Grid point 856 (103/120) =======================
q-point: ( 0.43  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.568   (  -0.767   -0.443   -6.731)    6.789
   2.568   (  -0.767   -0.443    6.731)    6.789
   2.688   (  -3.297   -1.903   -9.764)   10.479
   2.688   (  -3.297   -1.903    9.764)   10.479
   6.016   (   9.465    5.465   -5.159)   12.086
   6.016   (   9.465    5.465    5.159)   12.086
   9.775   (   0.028    0.016   -2.780)    2.780
   9.775   (   0.028    0.016    2.780)    2.780
   9.937   (  -0.534   -0.309   -5.095)    5.132
   9.937   (  -0.534   -0.309    5.095)    5.132
  11.296   (   0.422    0.244   -5.547)    5.568
  11.296   (   0.422    0.244    5.547)    5.568
======================= Grid point 858 (104/120) =======================
q-point: (-0.50  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.555   (   0.000    0.000   -6.408)    6.408
   2.555   (   0.000    0.000    6.408)    6.408
   2.648   (   0.000    0.000   -6.449)    6.449
   2.648   (   0.000    0.000    6.449)    6.449
   6.129   (   0.000    0.000   -3.872)    3.872
   6.129   (   0.000    0.000    3.872)    3.872
   9.777   (   0.000    0.000   -2.299)    2.299
   9.777   (   0.000    0.000    2.299)    2.299
   9.932   (   0.000    0.000   -5.267)    5.267
   9.932   (   0.000    0.000    5.267)    5.267
  11.301   (   0.000    0.000   -4.452)    4.452
  11.301   (   0.000    0.000    4.452)    4.452
======================= Grid point 876 (105/120) =======================
q-point: ( 0.07  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.981   (   8.085   14.004  -10.068)   19.049
   1.981   (   8.085   14.004   10.068)   19.049
   2.459   (  13.824   23.944  -11.617)   29.990
   2.459   (  13.824   23.944   11.617)   29.990
   4.317   (   1.208    2.092  -38.779)   38.854
   4.317   (   1.208    2.092   38.779)   38.854
  10.263   (  -4.611   -7.986   -0.909)    9.266
  10.263   (  -4.611   -7.986    0.909)    9.266
  10.369   (  -1.285   -2.225   -1.330)    2.893
  10.369   (  -1.285   -2.225    1.330)    2.893
  11.143   (  -0.866   -1.500   -8.610)    8.782
  11.143   (  -0.866   -1.500    8.610)    8.782
======================= Grid point 878 (106/120) =======================
q-point: ( 0.14  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.293   (   6.336   18.854   -8.867)   21.777
   2.293   (   6.336   18.854    8.867)   21.777
   2.857   (  10.368    6.131  -17.650)   21.368
   2.857   (  10.368    6.131   17.650)   21.368
   4.479   (  11.157   10.893  -26.934)   31.122
   4.479   (  11.157   10.893   26.934)   31.122
  10.077   (  -5.242   -9.513   -2.640)   11.178
  10.077   (  -5.242   -9.513    2.640)   11.178
  10.284   (  -6.650   -0.199   -3.255)    7.407
  10.284   (  -6.650   -0.199    3.255)    7.407
  11.137   (   2.103   -1.970   -8.678)    9.144
  11.137   (   2.103   -1.970    8.678)    9.144
======================= Grid point 880 (107/120) =======================
q-point: ( 0.21  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.551   (   0.357   17.845   -8.469)   19.756
   2.551   (   0.357   17.845    8.469)   19.756
   2.951   (   0.540   -2.722  -21.036)   21.218
   2.951   (   0.540   -2.722   21.036)   21.218
   4.919   (  21.068   14.200  -15.262)   29.639
   4.919   (  21.068   14.200   15.262)   29.639
   9.910   (  -4.076   -5.875   -3.702)    8.052
   9.910   (  -4.076   -5.875    3.702)    8.052
  10.148   (  -8.258    0.894   -4.139)    9.280
  10.148   (  -8.258    0.894    4.139)    9.280
  11.164   (   4.312   -4.636   -8.397)   10.516
  11.164   (   4.312   -4.636    8.397)   10.516
======================= Grid point 882 (108/120) =======================
q-point: ( 0.29  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.651   (  -3.600    7.334   -9.785)   12.747
   2.651   (  -3.600    7.334    9.785)   12.747
   2.905   (  -5.561    4.972  -14.510)   16.316
   2.905   (  -5.561    4.972   14.510)   16.316
   5.438   (  21.106   10.650   -9.530)   25.490
   5.438   (  21.106   10.650    9.530)   25.490
   9.816   (  -2.375   -1.541   -3.806)    4.744
   9.816   (  -2.375   -1.541    3.806)    4.744
  10.030   (  -5.670    0.902   -4.691)    7.414
  10.030   (  -5.670    0.902    4.691)    7.414
  11.187   (   4.952   -6.913   -7.706)   11.476
  11.187   (   4.952   -6.913    7.706)   11.476
======================= Grid point 884 (109/120) =======================
q-point: ( 0.36  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.615   (  -2.725   -0.690   -9.966)   10.355
   2.615   (  -2.725   -0.690    9.966)   10.355
   2.861   ( -10.181   13.196   -7.865)   18.430
   2.861   ( -10.181   13.196    7.865)   18.430
   5.851   (  15.029    5.566   -6.264)   17.207
   5.851   (  15.029    5.566    6.264)   17.207
   9.782   (  -0.892    0.100   -2.933)    3.067
   9.782   (  -0.892    0.100    2.933)    3.067
   9.968   (  -2.596    1.132   -5.033)    5.775
   9.968   (  -2.596    1.132    5.033)    5.775
  11.199   (   5.169   -8.144   -6.309)   11.526
  11.199   (   5.169   -8.144    6.309)   11.526
======================= Grid point 886 (110/120) =======================
q-point: ( 0.43  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.572   (  -0.705   -1.145   -7.429)    7.550
   2.572   (  -0.705   -1.145    7.429)    7.550
   2.820   ( -10.572   15.744   -5.856)   19.847
   2.820   ( -10.572   15.744    5.856)   19.847
   6.073   (   6.267    0.149   -4.062)    7.470
   6.073   (   6.267    0.149    4.062)    7.470
   9.776   (  -0.034    0.048   -2.081)    2.082
   9.776   (  -0.034    0.048    2.081)    2.082
   9.949   (  -1.181    1.335   -5.225)    5.521
   9.949   (  -1.181    1.335    5.225)    5.521
  11.205   (   5.197   -8.690   -4.685)   11.157
  11.205   (   5.197   -8.690    4.685)   11.157
======================= Grid point 906 (111/120) =======================
q-point: ( 0.14  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.703   (  10.178   17.628   -6.400)   21.338
   2.703   (  10.178   17.628    6.400)   21.338
   2.889   (  -0.516   -0.894  -20.187)   20.214
   2.889   (  -0.516   -0.894   20.187)   20.214
   4.791   (  10.787   18.684  -17.204)   27.594
   4.791   (  10.787   18.684   17.204)   27.594
   9.914   (  -3.425   -5.933   -2.557)    7.312
   9.914   (  -3.425   -5.933    2.557)    7.312
  10.251   (  -1.886   -3.267   -4.259)    5.689
  10.251   (  -1.886   -3.267    4.259)    5.689
  11.079   (  -1.828   -3.165   -8.227)    9.003
  11.079   (  -1.828   -3.165    8.227)    9.003
======================= Grid point 908 (112/120) =======================
q-point: ( 0.21  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.751   (  -4.549   -3.328  -16.885)   17.801
   2.751   (  -4.549   -3.328   16.885)   17.801
   3.075   (   3.044   18.310   -6.158)   19.556
   3.075   (   3.044   18.310    6.158)   19.556
   5.212   (  16.280   14.101  -10.989)   24.179
   5.212   (  16.280   14.101   10.989)   24.179
   9.833   (  -1.483   -1.726   -2.309)    3.242
   9.833   (  -1.483   -1.726    2.309)    3.242
  10.155   (  -6.036   -0.620   -4.843)    7.763
  10.155   (  -6.036   -0.620    4.843)    7.763
  11.021   (   2.190   -8.388   -7.614)   11.538
  11.021   (   2.190   -8.388    7.614)   11.538
======================= Grid point 910 (113/120) =======================
q-point: ( 0.29  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 492
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.638   (  -3.037   -4.285  -13.843)   14.806
   2.638   (  -3.037   -4.285   13.843)   14.806
   3.232   (  -7.815   20.875   -5.805)   23.034
   3.232   (  -7.815   20.875    5.805)   23.034
   5.625   (  16.406    7.621   -7.260)   19.492
   5.625   (  16.406    7.621    7.260)   19.492
   9.805   (  -1.287    0.118   -2.162)    2.519
   9.805   (  -1.287    0.118    2.162)    2.519
  10.051   (  -5.055    0.818   -5.203)    7.301
  10.051   (  -5.055    0.818    5.203)    7.301
  10.986   (   5.412  -11.696   -6.621)   14.489
  10.986   (   5.412  -11.696    6.621)   14.489
======================= Grid point 912 (114/120) =======================
q-point: ( 0.36  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.564   (  -1.052   -2.909   -9.217)    9.722
   2.564   (  -1.052   -2.909    9.217)    9.722
   3.255   ( -13.217   21.961   -5.280)   26.170
   3.255   ( -13.217   21.961    5.280)   26.170
   5.922   (  10.905    1.630   -4.561)   11.932
   5.922   (  10.905    1.630    4.561)   11.932
   9.784   (  -0.928   -0.095   -1.539)    1.799
   9.784   (  -0.928   -0.095    1.539)    1.799
   9.998   (  -2.456    1.461   -5.236)    5.965
   9.998   (  -2.456    1.461    5.236)    5.965
  10.971   (   7.083  -13.069   -5.070)   15.705
  10.971   (   7.083  -13.069    5.070)   15.705
======================= Grid point 914 (115/120) =======================
q-point: ( 0.43  0.14 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.540   (   0.893   -1.547   -6.562)    6.800
   2.540   (   0.893   -1.547    6.562)    6.800
   3.246   ( -13.322   23.075   -5.032)   27.116
   3.246   ( -13.322   23.075    5.032)   27.116
   6.030   (   2.289   -3.965   -3.181)    5.575
   6.030   (   2.289   -3.965    3.181)    5.575
   9.775   (   0.084   -0.145   -1.151)    1.164
   9.775   (   0.084   -0.145    1.151)    1.164
   9.984   (  -1.083    1.876   -5.240)    5.671
   9.984   (  -1.083    1.876    5.240)    5.671
  10.967   (   7.708  -13.350   -4.092)   15.949
  10.967   (   7.708  -13.350    4.092)   15.949
======================= Grid point 938 (116/120) =======================
q-point: ( 0.21  0.21 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.642   (  -3.064   -5.307  -15.110)   16.305
   2.642   (  -3.064   -5.307   15.110)   16.305
   3.434   (   8.103   14.035   -4.977)   16.953
   3.434   (   8.103   14.035    4.977)   16.953
   5.500   (   7.945   13.761   -7.751)   17.679
   5.500   (   7.945   13.761    7.751)   17.679
   9.819   (  -0.166   -0.287   -1.050)    1.101
   9.819   (  -0.166   -0.287    1.050)    1.101
  10.117   (  -1.678   -2.906   -5.285)    6.260
  10.117   (  -1.678   -2.906    5.285)    6.260
  10.855   (  -4.022   -6.966   -6.608)   10.410
  10.855   (  -4.022   -6.966    6.608)   10.410
======================= Grid point 940 (117/120) =======================
q-point: ( 0.29  0.21 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.555   (  -1.619   -3.301  -10.552)   11.174
   2.555   (  -1.619   -3.301   10.552)   11.174
   3.628   (  -3.339   16.089   -4.791)   17.116
   3.628   (  -3.339   16.089    4.791)   17.116
   5.764   (   9.441    6.102   -5.040)   12.320
   5.764   (   9.441    6.102    5.040)   12.320
   9.802   (  -1.233   -0.732   -0.683)    1.588
   9.802   (  -1.233   -0.732    0.683)    1.588
  10.055   (  -2.675   -0.430   -5.418)    6.058
  10.055   (  -2.675   -0.430    5.418)    6.058
  10.749   (   2.045  -10.270   -5.169)   11.678
  10.749   (   2.045  -10.270    5.169)   11.678
======================= Grid point 942 (118/120) =======================
q-point: ( 0.36  0.21 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 490
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.515   (  -0.065   -1.703   -6.107)    6.340
   2.515   (  -0.065   -1.703    6.107)    6.340
   3.680   (  -9.372   17.693   -4.614)   20.546
   3.680   (  -9.372   17.693    4.614)   20.546
   5.932   (   5.309   -0.348   -2.808)    6.016
   5.932   (   5.309   -0.348    2.808)    6.016
   9.773   (  -0.545   -0.964   -0.380)    1.170
   9.773   (  -0.545   -0.964    0.380)    1.170
  10.021   (  -1.352    0.696   -5.306)    5.519
  10.021   (  -1.352    0.696    5.306)    5.519
  10.708   (   5.509  -11.297   -3.582)   13.070
  10.708   (   5.509  -11.297    3.582)   13.070
======================= Grid point 968 (119/120) =======================
q-point: ( 0.29  0.29 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 267
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.507   (  -0.807   -1.397   -6.037)    6.249
   2.507   (  -0.807   -1.397    6.037)    6.249
   3.868   (   4.141    7.173   -4.437)    9.396
   3.868   (   4.141    7.173    4.437)    9.396
   5.874   (   2.467    4.272   -2.777)    5.661
   5.874   (   2.467    4.272    2.777)    5.661
   9.774   (  -0.992   -1.719   -0.028)    1.985
   9.774   (  -0.992   -1.719    0.028)    1.985
  10.038   (  -0.534   -0.925   -5.468)    5.572
  10.038   (  -0.534   -0.925    5.468)    5.572
  10.591   (  -2.814   -4.875   -3.123)    6.437
  10.591   (  -2.814   -4.875    3.123)    6.437
======================= Grid point 970 (120/120) =======================
q-point: ( 0.36  0.29 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 
Number of triplets: 266
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.493   (   0.270   -0.468   -2.510)    2.567
   2.493   (   0.270   -0.468    2.510)    2.567
   3.938   (  -3.828    6.631   -4.325)    8.794
   3.938   (  -3.828    6.631    4.325)    8.794
   5.922   (   0.216   -0.375   -1.119)    1.200
   5.922   (   0.216   -0.375    1.119)    1.200
   9.752   (   0.435   -0.754   -0.044)    0.872
   9.752   (   0.435   -0.754    0.044)    0.872
  10.027   (  -0.032    0.055   -5.430)    5.431
  10.027   (  -0.032    0.055    5.430)    5.431
  10.543   (   2.417   -4.186   -1.524)    5.068
  10.543   (   2.417   -4.186    1.524)    5.068
=================== 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/18816
   10.0  10015.749  10015.749  12305.263     -0.000      0.000      0.000 3/18816
   20.0   3496.346   3496.346   3657.321      0.000      0.000      0.000 3/18816
   30.0   1679.700   1679.700   1710.971     -0.000     -0.000      0.000 3/18816
   40.0    978.737    978.737    994.063     -0.000     -0.000      0.000 3/18816
   50.0    666.394    666.394    676.122     -0.000     -0.000      0.000 3/18816
   60.0    499.017    499.017    505.043     -0.000     -0.000      0.000 3/18816
   70.0    395.909    395.909    398.829     -0.000     -0.000      0.000 3/18816
   80.0    326.323    326.323    326.758     -0.000     -0.000      0.000 3/18816
   90.0    276.517    276.517    275.183     -0.000      0.000      0.000 3/18816
  100.0    239.391    239.391    236.925     -0.000      0.000      0.000 3/18816
  110.0    210.847    210.847    207.726     -0.000      0.000      0.000 3/18816
  120.0    188.337    188.337    184.881     -0.000      0.000      0.000 3/18816
  130.0    170.195    170.195    166.608     -0.000      0.000      0.000 3/18816
  140.0    155.295    155.295    151.700     -0.000      0.000      0.000 3/18816
  150.0    142.854    142.854    139.324     -0.000      0.000      0.000 3/18816
  160.0    132.316    132.316    128.890     -0.000      0.000      0.000 3/18816
  170.0    123.276    123.276    119.975     -0.000      0.000      0.000 3/18816
  180.0    115.435    115.435    112.266     -0.000      0.000      0.000 3/18816
  190.0    108.567    108.567    105.531     -0.000      0.000      0.000 3/18816
  200.0    102.498    102.498     99.593     -0.000      0.000      0.000 3/18816
  210.0     97.096     97.096     94.316     -0.000      0.000      0.000 3/18816
  220.0     92.254     92.254     89.593     -0.000      0.000      0.000 3/18816
  230.0     87.887     87.887     85.338     -0.000      0.000      0.000 3/18816
  240.0     83.929     83.929     81.484     -0.000      0.000      0.000 3/18816
  250.0     80.321     80.321     77.975     -0.000      0.000      0.000 3/18816
  260.0     77.020     77.020     74.766     -0.000      0.000      0.000 3/18816
  270.0     73.987     73.987     71.818     -0.000      0.000      0.000 3/18816
  280.0     71.190     71.190     69.101     -0.000      0.000      0.000 3/18816
  290.0     68.601     68.601     66.587     -0.000      0.000      0.000 3/18816
  300.0     66.199     66.199     64.255     -0.000      0.000      0.000 3/18816
  310.0     63.963     63.963     62.084     -0.000      0.000      0.000 3/18816
  320.0     61.875     61.875     60.058     -0.000      0.000      0.000 3/18816
  330.0     59.923     59.923     58.164     -0.000      0.000      0.000 3/18816
  340.0     58.092     58.092     56.387     -0.000      0.000      0.000 3/18816
  350.0     56.372     56.372     54.718     -0.000      0.000      0.000 3/18816
  360.0     54.752     54.752     53.147     -0.000      0.000      0.000 3/18816
  370.0     53.224     53.224     51.664     -0.000      0.000      0.000 3/18816
  380.0     51.780     51.780     50.264     -0.000      0.000      0.000 3/18816
  390.0     50.414     50.414     48.939     -0.000      0.000      0.000 3/18816
  400.0     49.119     49.119     47.682     -0.000      0.000      0.000 3/18816
  410.0     47.890     47.890     46.490     -0.000      0.000      0.000 3/18816
  420.0     46.721     46.721     45.356     -0.000      0.000      0.000 3/18816
  430.0     45.609     45.609     44.277     -0.000      0.000      0.000 3/18816
  440.0     44.549     44.549     43.249     -0.000      0.000      0.000 3/18816
  450.0     43.538     43.538     42.268     -0.000      0.000      0.000 3/18816
  460.0     42.572     42.572     41.331     -0.000      0.000      0.000 3/18816
  470.0     41.648     41.648     40.435     -0.000      0.000      0.000 3/18816
  480.0     40.764     40.764     39.577     -0.000      0.000      0.000 3/18816
  490.0     39.917     39.917     38.755     -0.000      0.000      0.000 3/18816
  500.0     39.105     39.105     37.968     -0.000      0.000      0.000 3/18816
  510.0     38.326     38.326     37.211     -0.000      0.000      0.000 3/18816
  520.0     37.577     37.577     36.485     -0.000      0.000      0.000 3/18816
  530.0     36.857     36.857     35.786     -0.000      0.000      0.000 3/18816
  540.0     36.164     36.164     35.114     -0.000      0.000      0.000 3/18816
  550.0     35.498     35.498     34.467     -0.000      0.000      0.000 3/18816
  560.0     34.855     34.855     33.844     -0.000      0.000      0.000 3/18816
  570.0     34.236     34.236     33.243     -0.000      0.000      0.000 3/18816
  580.0     33.638     33.638     32.663     -0.000      0.000      0.000 3/18816
  590.0     33.061     33.061     32.102     -0.000      0.000      0.000 3/18816
  600.0     32.503     32.503     31.561     -0.000      0.000      0.000 3/18816
  610.0     31.964     31.964     31.038     -0.000      0.000      0.000 3/18816
  620.0     31.443     31.443     30.533     -0.000      0.000      0.000 3/18816
  630.0     30.939     30.939     30.043     -0.000      0.000      0.000 3/18816
  640.0     30.450     30.450     29.569     -0.000      0.000      0.000 3/18816
  650.0     29.977     29.977     29.110     -0.000      0.000      0.000 3/18816
  660.0     29.518     29.518     28.665     -0.000      0.000      0.000 3/18816
  670.0     29.074     29.074     28.233     -0.000      0.000      0.000 3/18816
  680.0     28.642     28.642     27.814     -0.000      0.000      0.000 3/18816
  690.0     28.223     28.223     27.408     -0.000      0.000      0.000 3/18816
  700.0     27.817     27.817     27.013     -0.000      0.000      0.000 3/18816
  710.0     27.422     27.422     26.630     -0.000      0.000      0.000 3/18816
  720.0     27.038     27.038     26.257     -0.000      0.000      0.000 3/18816
  730.0     26.664     26.664     25.895     -0.000      0.000      0.000 3/18816
  740.0     26.301     26.301     25.542     -0.000      0.000      0.000 3/18816
  750.0     25.948     25.948     25.199     -0.000      0.000      0.000 3/18816
  760.0     25.604     25.604     24.865     -0.000      0.000      0.000 3/18816
  770.0     25.269     25.269     24.540     -0.000      0.000      0.000 3/18816
  780.0     24.943     24.943     24.224     -0.000      0.000      0.000 3/18816
  790.0     24.625     24.625     23.915     -0.000      0.000      0.000 3/18816
  800.0     24.315     24.315     23.614     -0.000      0.000      0.000 3/18816
  810.0     24.013     24.013     23.321     -0.000      0.000      0.000 3/18816
  820.0     23.718     23.718     23.035     -0.000      0.000      0.000 3/18816
  830.0     23.431     23.431     22.756     -0.000      0.000      0.000 3/18816
  840.0     23.150     23.150     22.483     -0.000      0.000      0.000 3/18816
  850.0     22.876     22.876     22.217     -0.000      0.000      0.000 3/18816
  860.0     22.609     22.609     21.958     -0.000      0.000      0.000 3/18816
  870.0     22.348     22.348     21.704     -0.000      0.000      0.000 3/18816
  880.0     22.092     22.092     21.456     -0.000      0.000      0.000 3/18816
  890.0     21.843     21.843     21.214     -0.000      0.000      0.000 3/18816
  900.0     21.599     21.599     20.977     -0.000      0.000      0.000 3/18816
  910.0     21.360     21.360     20.745     -0.000      0.000      0.000 3/18816
  920.0     21.127     21.127     20.519     -0.000      0.000      0.000 3/18816
  930.0     20.899     20.899     20.297     -0.000      0.000      0.000 3/18816
  940.0     20.675     20.675     20.080     -0.000      0.000      0.000 3/18816
  950.0     20.457     20.457     19.868     -0.000      0.000      0.000 3/18816
  960.0     20.243     20.243     19.660     -0.000      0.000      0.000 3/18816
  970.0     20.033     20.033     19.457     -0.000      0.000      0.000 3/18816
  980.0     19.828     19.828     19.257     -0.000      0.000      0.000 3/18816
  990.0     19.627     19.627     19.062     -0.000      0.000      0.000 3/18816
 1000.0     19.429     19.429     18.871     -0.000      0.000      0.000 3/18816

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

