
-----------------------------
------- 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    3.862004957877586    0.000000000000000    0.000000000000000
  b   -1.931002478938793    3.344594403063441    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.338390410000000
Atomic positions (fractional):
   *1 Al  0.33333333333333  0.66666666666667  0.00005672915003  26.982
    2 Al  0.66666666666667  0.33333333333333  0.50005672915003  26.982
   *3 P   0.33333333333333  0.66666666666667  0.37494327084997  30.974
    4 P   0.66666666666667  0.33333333333333  0.87494327084997  30.974
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.862004957877586    0.000000000000000    0.000000000000000
  b   -1.931002478938793    3.344594403063441    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.338390410000000
Atomic positions (fractional):
   *1 Al  0.33333333333333  0.66666666666667  0.00005672915003  26.982 > 1
    2 Al  0.66666666666667  0.33333333333333  0.50005672915003  26.982 > 2
   *3 P   0.33333333333333  0.66666666666667  0.37494327084997  30.974 > 3
    4 P   0.66666666666667  0.33333333333333  0.87494327084997  30.974 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   15.448019831510344    0.000000000000000    0.000000000000000
  b   -7.724009915755172   13.378377612253765    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.676780819999999
Atomic positions (fractional):
   *1 Al  0.08333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    2 Al  0.33333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    3 Al  0.58333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    4 Al  0.83333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    5 Al  0.08333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    6 Al  0.33333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    7 Al  0.58333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    8 Al  0.83333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    9 Al  0.08333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   10 Al  0.33333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   11 Al  0.58333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   12 Al  0.83333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   13 Al  0.08333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   14 Al  0.33333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   15 Al  0.58333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   16 Al  0.83333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   17 Al  0.08333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   18 Al  0.33333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   19 Al  0.58333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   20 Al  0.83333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   21 Al  0.08333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   22 Al  0.33333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   23 Al  0.58333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   24 Al  0.83333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   25 Al  0.08333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   26 Al  0.33333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   27 Al  0.58333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   28 Al  0.83333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   29 Al  0.08333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   30 Al  0.33333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   31 Al  0.58333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   32 Al  0.83333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   33 Al  0.16666666666667  0.08333333333333  0.25002836457501  26.982 > 2
   34 Al  0.41666666666667  0.08333333333333  0.25002836457501  26.982 > 2
   35 Al  0.66666666666667  0.08333333333333  0.25002836457501  26.982 > 2
   36 Al  0.91666666666667  0.08333333333333  0.25002836457501  26.982 > 2
   37 Al  0.16666666666667  0.33333333333333  0.25002836457501  26.982 > 2
   38 Al  0.41666666666667  0.33333333333333  0.25002836457501  26.982 > 2
   39 Al  0.66666666666667  0.33333333333333  0.25002836457501  26.982 > 2
   40 Al  0.91666666666667  0.33333333333333  0.25002836457501  26.982 > 2
   41 Al  0.16666666666667  0.58333333333333  0.25002836457501  26.982 > 2
   42 Al  0.41666666666667  0.58333333333333  0.25002836457501  26.982 > 2
   43 Al  0.66666666666667  0.58333333333333  0.25002836457501  26.982 > 2
   44 Al  0.91666666666667  0.58333333333333  0.25002836457501  26.982 > 2
   45 Al  0.16666666666667  0.83333333333333  0.25002836457501  26.982 > 2
   46 Al  0.41666666666667  0.83333333333333  0.25002836457501  26.982 > 2
   47 Al  0.66666666666667  0.83333333333333  0.25002836457501  26.982 > 2
   48 Al  0.91666666666667  0.83333333333333  0.25002836457501  26.982 > 2
   49 Al  0.16666666666667  0.08333333333333  0.75002836457501  26.982 > 2
   50 Al  0.41666666666667  0.08333333333333  0.75002836457501  26.982 > 2
   51 Al  0.66666666666667  0.08333333333333  0.75002836457501  26.982 > 2
   52 Al  0.91666666666667  0.08333333333333  0.75002836457501  26.982 > 2
   53 Al  0.16666666666667  0.33333333333333  0.75002836457501  26.982 > 2
   54 Al  0.41666666666667  0.33333333333333  0.75002836457501  26.982 > 2
   55 Al  0.66666666666667  0.33333333333333  0.75002836457501  26.982 > 2
   56 Al  0.91666666666667  0.33333333333333  0.75002836457501  26.982 > 2
   57 Al  0.16666666666667  0.58333333333333  0.75002836457501  26.982 > 2
   58 Al  0.41666666666667  0.58333333333333  0.75002836457501  26.982 > 2
   59 Al  0.66666666666667  0.58333333333333  0.75002836457501  26.982 > 2
   60 Al  0.91666666666667  0.58333333333333  0.75002836457501  26.982 > 2
   61 Al  0.16666666666667  0.83333333333333  0.75002836457501  26.982 > 2
   62 Al  0.41666666666667  0.83333333333333  0.75002836457501  26.982 > 2
   63 Al  0.66666666666667  0.83333333333333  0.75002836457501  26.982 > 2
   64 Al  0.91666666666667  0.83333333333333  0.75002836457501  26.982 > 2
  *65 P   0.08333333333333  0.16666666666667  0.18747163542499  30.974 > 3
   66 P   0.33333333333333  0.16666666666667  0.18747163542499  30.974 > 3
   67 P   0.58333333333333  0.16666666666667  0.18747163542499  30.974 > 3
   68 P   0.83333333333333  0.16666666666667  0.18747163542499  30.974 > 3
   69 P   0.08333333333333  0.41666666666667  0.18747163542499  30.974 > 3
   70 P   0.33333333333333  0.41666666666667  0.18747163542499  30.974 > 3
   71 P   0.58333333333333  0.41666666666667  0.18747163542499  30.974 > 3
   72 P   0.83333333333333  0.41666666666667  0.18747163542499  30.974 > 3
   73 P   0.08333333333333  0.66666666666667  0.18747163542499  30.974 > 3
   74 P   0.33333333333333  0.66666666666667  0.18747163542499  30.974 > 3
   75 P   0.58333333333333  0.66666666666667  0.18747163542499  30.974 > 3
   76 P   0.83333333333333  0.66666666666667  0.18747163542499  30.974 > 3
   77 P   0.08333333333333  0.91666666666667  0.18747163542499  30.974 > 3
   78 P   0.33333333333333  0.91666666666667  0.18747163542499  30.974 > 3
   79 P   0.58333333333333  0.91666666666667  0.18747163542499  30.974 > 3
   80 P   0.83333333333333  0.91666666666667  0.18747163542499  30.974 > 3
   81 P   0.08333333333333  0.16666666666667  0.68747163542499  30.974 > 3
   82 P   0.33333333333333  0.16666666666667  0.68747163542499  30.974 > 3
   83 P   0.58333333333333  0.16666666666667  0.68747163542499  30.974 > 3
   84 P   0.83333333333333  0.16666666666667  0.68747163542499  30.974 > 3
   85 P   0.08333333333333  0.41666666666667  0.68747163542499  30.974 > 3
   86 P   0.33333333333333  0.41666666666667  0.68747163542499  30.974 > 3
   87 P   0.58333333333333  0.41666666666667  0.68747163542499  30.974 > 3
   88 P   0.83333333333333  0.41666666666667  0.68747163542499  30.974 > 3
   89 P   0.08333333333333  0.66666666666667  0.68747163542499  30.974 > 3
   90 P   0.33333333333333  0.66666666666667  0.68747163542499  30.974 > 3
   91 P   0.58333333333333  0.66666666666667  0.68747163542499  30.974 > 3
   92 P   0.83333333333333  0.66666666666667  0.68747163542499  30.974 > 3
   93 P   0.08333333333333  0.91666666666667  0.68747163542499  30.974 > 3
   94 P   0.33333333333333  0.91666666666667  0.68747163542499  30.974 > 3
   95 P   0.58333333333333  0.91666666666667  0.68747163542499  30.974 > 3
   96 P   0.83333333333333  0.91666666666667  0.68747163542499  30.974 > 3
   97 P   0.16666666666667  0.08333333333333  0.43747163542499  30.974 > 4
   98 P   0.41666666666667  0.08333333333333  0.43747163542499  30.974 > 4
   99 P   0.66666666666667  0.08333333333333  0.43747163542499  30.974 > 4
  100 P   0.91666666666667  0.08333333333333  0.43747163542499  30.974 > 4
  101 P   0.16666666666667  0.33333333333333  0.43747163542499  30.974 > 4
  102 P   0.41666666666667  0.33333333333333  0.43747163542499  30.974 > 4
  103 P   0.66666666666667  0.33333333333333  0.43747163542499  30.974 > 4
  104 P   0.91666666666667  0.33333333333333  0.43747163542499  30.974 > 4
  105 P   0.16666666666667  0.58333333333333  0.43747163542499  30.974 > 4
  106 P   0.41666666666667  0.58333333333333  0.43747163542499  30.974 > 4
  107 P   0.66666666666667  0.58333333333333  0.43747163542499  30.974 > 4
  108 P   0.91666666666667  0.58333333333333  0.43747163542499  30.974 > 4
  109 P   0.16666666666667  0.83333333333333  0.43747163542499  30.974 > 4
  110 P   0.41666666666667  0.83333333333333  0.43747163542499  30.974 > 4
  111 P   0.66666666666667  0.83333333333333  0.43747163542499  30.974 > 4
  112 P   0.91666666666667  0.83333333333333  0.43747163542499  30.974 > 4
  113 P   0.16666666666667  0.08333333333333  0.93747163542499  30.974 > 4
  114 P   0.41666666666667  0.08333333333333  0.93747163542499  30.974 > 4
  115 P   0.66666666666667  0.08333333333333  0.93747163542499  30.974 > 4
  116 P   0.91666666666667  0.08333333333333  0.93747163542499  30.974 > 4
  117 P   0.16666666666667  0.33333333333333  0.93747163542499  30.974 > 4
  118 P   0.41666666666667  0.33333333333333  0.93747163542499  30.974 > 4
  119 P   0.66666666666667  0.33333333333333  0.93747163542499  30.974 > 4
  120 P   0.91666666666667  0.33333333333333  0.93747163542499  30.974 > 4
  121 P   0.16666666666667  0.58333333333333  0.93747163542499  30.974 > 4
  122 P   0.41666666666667  0.58333333333333  0.93747163542499  30.974 > 4
  123 P   0.66666666666667  0.58333333333333  0.93747163542499  30.974 > 4
  124 P   0.91666666666667  0.58333333333333  0.93747163542499  30.974 > 4
  125 P   0.16666666666667  0.83333333333333  0.93747163542499  30.974 > 4
  126 P   0.41666666666667  0.83333333333333  0.93747163542499  30.974 > 4
  127 P   0.66666666666667  0.83333333333333  0.93747163542499  30.974 > 4
  128 P   0.91666666666667  0.83333333333333  0.93747163542499  30.974 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            9.0303717   -0.0000000    0.0000000
            0.0000000    9.0303717    0.0000000
            0.0000000    0.0000000    9.2322622
-------------------------- Born effective charges --------------------------
    1 Al    2.2508269   -0.0000000    0.0000000
            0.0000000    2.2508269    0.0000000
            0.0000000    0.0000000    2.4009123
    2 Al    2.2508269   -0.0000000    0.0000000
            0.0000000    2.2508269    0.0000000
            0.0000000    0.0000000    2.4009123
    3 P    -2.2508269    0.0000000    0.0000000
            0.0000000   -2.2508269    0.0000000
            0.0000000    0.0000000   -2.4009123
    4 P    -2.2508269    0.0000000    0.0000000
            0.0000000   -2.2508269    0.0000000
            0.0000000    0.0000000   -2.4009123
----------------------------------------------------------------------------
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.010
Solver_block: 80 / 80
 - Time: 0.158
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.169
--------------------------------- Symfc end --------------------------------
Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 156
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
Max drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) 
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    3.862004957877586    0.000000000000000    0.000000000000000
  b   -1.931002478938793    3.344594403063441    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.338390410000000
Atomic positions (fractional):
    1 Al  0.33333333333333  0.66666666666667  0.00005672915003  26.982
    2 Al  0.66666666666667  0.33333333333333  0.50005672915003  26.982
    3 P   0.33333333333333  0.66666666666667  0.37494327084997  30.974
    4 P   0.66666666666667  0.33333333333333  0.87494327084997  30.974
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.448019831510344    0.000000000000000    0.000000000000000
  b   -7.724009915755172   13.378377612253765    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.676780819999999
Atomic positions (fractional):
    1 Al  0.08333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    2 Al  0.33333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    3 Al  0.58333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    4 Al  0.83333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    5 Al  0.08333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    6 Al  0.33333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    7 Al  0.58333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    8 Al  0.83333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    9 Al  0.08333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   10 Al  0.33333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   11 Al  0.58333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   12 Al  0.83333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   13 Al  0.08333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   14 Al  0.33333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   15 Al  0.58333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   16 Al  0.83333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   17 Al  0.08333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   18 Al  0.33333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   19 Al  0.58333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   20 Al  0.83333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   21 Al  0.08333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   22 Al  0.33333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   23 Al  0.58333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   24 Al  0.83333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   25 Al  0.08333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   26 Al  0.33333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   27 Al  0.58333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   28 Al  0.83333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   29 Al  0.08333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   30 Al  0.33333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   31 Al  0.58333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   32 Al  0.83333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   33 Al  0.16666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   34 Al  0.41666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   35 Al  0.66666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   36 Al  0.91666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   37 Al  0.16666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   38 Al  0.41666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   39 Al  0.66666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   40 Al  0.91666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   41 Al  0.16666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   42 Al  0.41666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   43 Al  0.66666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   44 Al  0.91666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   45 Al  0.16666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   46 Al  0.41666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   47 Al  0.66666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   48 Al  0.91666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   49 Al  0.16666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   50 Al  0.41666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   51 Al  0.66666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   52 Al  0.91666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   53 Al  0.16666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   54 Al  0.41666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   55 Al  0.66666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   56 Al  0.91666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   57 Al  0.16666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   58 Al  0.41666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   59 Al  0.66666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   60 Al  0.91666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   61 Al  0.16666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   62 Al  0.41666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   63 Al  0.66666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   64 Al  0.91666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   65 P   0.08333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   66 P   0.33333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   67 P   0.58333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   68 P   0.83333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   69 P   0.08333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   70 P   0.33333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   71 P   0.58333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   72 P   0.83333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   73 P   0.08333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   74 P   0.33333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   75 P   0.58333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   76 P   0.83333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   77 P   0.08333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   78 P   0.33333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   79 P   0.58333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   80 P   0.83333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   81 P   0.08333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   82 P   0.33333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   83 P   0.58333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   84 P   0.83333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   85 P   0.08333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   86 P   0.33333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   87 P   0.58333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   88 P   0.83333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   89 P   0.08333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   90 P   0.33333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   91 P   0.58333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   92 P   0.83333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   93 P   0.08333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   94 P   0.33333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   95 P   0.58333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   96 P   0.83333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   97 P   0.16666666666667  0.08333333333333  0.43747163542499  30.974 > 97
   98 P   0.41666666666667  0.08333333333333  0.43747163542499  30.974 > 97
   99 P   0.66666666666667  0.08333333333333  0.43747163542499  30.974 > 97
  100 P   0.91666666666667  0.08333333333333  0.43747163542499  30.974 > 97
  101 P   0.16666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  102 P   0.41666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  103 P   0.66666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  104 P   0.91666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  105 P   0.16666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  106 P   0.41666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  107 P   0.66666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  108 P   0.91666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  109 P   0.16666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  110 P   0.41666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  111 P   0.66666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  112 P   0.91666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  113 P   0.16666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  114 P   0.41666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  115 P   0.66666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  116 P   0.91666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  117 P   0.16666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  118 P   0.41666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  119 P   0.66666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  120 P   0.91666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  121 P   0.16666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  122 P   0.41666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  123 P   0.66666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  124 P   0.91666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  125 P   0.16666666666667  0.83333333333333  0.93747163542499  30.974 > 97
  126 P   0.41666666666667  0.83333333333333  0.93747163542499  30.974 > 97
  127 P   0.66666666666667  0.83333333333333  0.93747163542499  30.974 > 97
  128 P   0.91666666666667  0.83333333333333  0.93747163542499  30.974 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            9.0303717   -0.0000000    0.0000000
            0.0000000    9.0303717    0.0000000
            0.0000000    0.0000000    9.2322622
-------------------------- Born effective charges --------------------------
    1 Al    2.2508269   -0.0000000    0.0000000
            0.0000000    2.2508269    0.0000000
            0.0000000    0.0000000    2.4009123
    2 Al    2.2508269   -0.0000000    0.0000000
            0.0000000    2.2508269    0.0000000
            0.0000000    0.0000000    2.4009123
    3 P    -2.2508269    0.0000000    0.0000000
            0.0000000   -2.2508269    0.0000000
            0.0000000    0.0000000   -2.4009123
    4 P    -2.2508269    0.0000000    0.0000000
            0.0000000   -2.2508269    0.0000000
            0.0000000    0.0000000   -2.4009123
----------------------------------------------------------------------------
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.00000048 (yyy) 0.00000048 (yyy) 0.00000048 (yyy)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-07 23:04: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    3.862004957877586    0.000000000000000    0.000000000000000
  b   -1.931002478938793    3.344594403063441    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.338390410000000
Atomic positions (fractional):
    1 Al  0.33333333333333  0.66666666666667  0.00005672915003  26.982
    2 Al  0.66666666666667  0.33333333333333  0.50005672915003  26.982
    3 P   0.33333333333333  0.66666666666667  0.37494327084997  30.974
    4 P   0.66666666666667  0.33333333333333  0.87494327084997  30.974
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.448019831510344    0.000000000000000    0.000000000000000
  b   -7.724009915755172   13.378377612253765    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.676780819999999
Atomic positions (fractional):
    1 Al  0.08333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    2 Al  0.33333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    3 Al  0.58333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    4 Al  0.83333333333333  0.16666666666667  0.00002836457501  26.982 > 1
    5 Al  0.08333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    6 Al  0.33333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    7 Al  0.58333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    8 Al  0.83333333333333  0.41666666666667  0.00002836457501  26.982 > 1
    9 Al  0.08333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   10 Al  0.33333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   11 Al  0.58333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   12 Al  0.83333333333333  0.66666666666667  0.00002836457501  26.982 > 1
   13 Al  0.08333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   14 Al  0.33333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   15 Al  0.58333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   16 Al  0.83333333333333  0.91666666666667  0.00002836457501  26.982 > 1
   17 Al  0.08333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   18 Al  0.33333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   19 Al  0.58333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   20 Al  0.83333333333333  0.16666666666667  0.50002836457501  26.982 > 1
   21 Al  0.08333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   22 Al  0.33333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   23 Al  0.58333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   24 Al  0.83333333333333  0.41666666666667  0.50002836457501  26.982 > 1
   25 Al  0.08333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   26 Al  0.33333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   27 Al  0.58333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   28 Al  0.83333333333333  0.66666666666667  0.50002836457501  26.982 > 1
   29 Al  0.08333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   30 Al  0.33333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   31 Al  0.58333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   32 Al  0.83333333333333  0.91666666666667  0.50002836457501  26.982 > 1
   33 Al  0.16666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   34 Al  0.41666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   35 Al  0.66666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   36 Al  0.91666666666667  0.08333333333333  0.25002836457501  26.982 > 33
   37 Al  0.16666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   38 Al  0.41666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   39 Al  0.66666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   40 Al  0.91666666666667  0.33333333333333  0.25002836457501  26.982 > 33
   41 Al  0.16666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   42 Al  0.41666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   43 Al  0.66666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   44 Al  0.91666666666667  0.58333333333333  0.25002836457501  26.982 > 33
   45 Al  0.16666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   46 Al  0.41666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   47 Al  0.66666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   48 Al  0.91666666666667  0.83333333333333  0.25002836457501  26.982 > 33
   49 Al  0.16666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   50 Al  0.41666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   51 Al  0.66666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   52 Al  0.91666666666667  0.08333333333333  0.75002836457501  26.982 > 33
   53 Al  0.16666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   54 Al  0.41666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   55 Al  0.66666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   56 Al  0.91666666666667  0.33333333333333  0.75002836457501  26.982 > 33
   57 Al  0.16666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   58 Al  0.41666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   59 Al  0.66666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   60 Al  0.91666666666667  0.58333333333333  0.75002836457501  26.982 > 33
   61 Al  0.16666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   62 Al  0.41666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   63 Al  0.66666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   64 Al  0.91666666666667  0.83333333333333  0.75002836457501  26.982 > 33
   65 P   0.08333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   66 P   0.33333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   67 P   0.58333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   68 P   0.83333333333333  0.16666666666667  0.18747163542499  30.974 > 65
   69 P   0.08333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   70 P   0.33333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   71 P   0.58333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   72 P   0.83333333333333  0.41666666666667  0.18747163542499  30.974 > 65
   73 P   0.08333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   74 P   0.33333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   75 P   0.58333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   76 P   0.83333333333333  0.66666666666667  0.18747163542499  30.974 > 65
   77 P   0.08333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   78 P   0.33333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   79 P   0.58333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   80 P   0.83333333333333  0.91666666666667  0.18747163542499  30.974 > 65
   81 P   0.08333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   82 P   0.33333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   83 P   0.58333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   84 P   0.83333333333333  0.16666666666667  0.68747163542499  30.974 > 65
   85 P   0.08333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   86 P   0.33333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   87 P   0.58333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   88 P   0.83333333333333  0.41666666666667  0.68747163542499  30.974 > 65
   89 P   0.08333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   90 P   0.33333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   91 P   0.58333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   92 P   0.83333333333333  0.66666666666667  0.68747163542499  30.974 > 65
   93 P   0.08333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   94 P   0.33333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   95 P   0.58333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   96 P   0.83333333333333  0.91666666666667  0.68747163542499  30.974 > 65
   97 P   0.16666666666667  0.08333333333333  0.43747163542499  30.974 > 97
   98 P   0.41666666666667  0.08333333333333  0.43747163542499  30.974 > 97
   99 P   0.66666666666667  0.08333333333333  0.43747163542499  30.974 > 97
  100 P   0.91666666666667  0.08333333333333  0.43747163542499  30.974 > 97
  101 P   0.16666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  102 P   0.41666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  103 P   0.66666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  104 P   0.91666666666667  0.33333333333333  0.43747163542499  30.974 > 97
  105 P   0.16666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  106 P   0.41666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  107 P   0.66666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  108 P   0.91666666666667  0.58333333333333  0.43747163542499  30.974 > 97
  109 P   0.16666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  110 P   0.41666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  111 P   0.66666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  112 P   0.91666666666667  0.83333333333333  0.43747163542499  30.974 > 97
  113 P   0.16666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  114 P   0.41666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  115 P   0.66666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  116 P   0.91666666666667  0.08333333333333  0.93747163542499  30.974 > 97
  117 P   0.16666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  118 P   0.41666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  119 P   0.66666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  120 P   0.91666666666667  0.33333333333333  0.93747163542499  30.974 > 97
  121 P   0.16666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  122 P   0.41666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  123 P   0.66666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  124 P   0.91666666666667  0.58333333333333  0.93747163542499  30.974 > 97
  125 P   0.16666666666667  0.83333333333333  0.93747163542499  30.974 > 97
  126 P   0.41666666666667  0.83333333333333  0.93747163542499  30.974 > 97
  127 P   0.66666666666667  0.83333333333333  0.93747163542499  30.974 > 97
  128 P   0.91666666666667  0.83333333333333  0.93747163542499  30.974 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            9.0303717   -0.0000000    0.0000000
            0.0000000    9.0303717    0.0000000
            0.0000000    0.0000000    9.2322622
-------------------------- Born effective charges --------------------------
    1 Al    2.2508269   -0.0000000    0.0000000
            0.0000000    2.2508269    0.0000000
            0.0000000    0.0000000    2.4009123
    2 Al    2.2508269   -0.0000000    0.0000000
            0.0000000    2.2508269    0.0000000
            0.0000000    0.0000000    2.4009123
    3 P    -2.2508269    0.0000000    0.0000000
            0.0000000   -2.2508269    0.0000000
            0.0000000    0.0000000   -2.4009123
    4 P    -2.2508269    0.0000000    0.0000000
            0.0000000   -2.2508269    0.0000000
            0.0000000    0.0000000   -2.4009123
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000048 (yyy) 0.00000048 (yyy) 0.00000048 (yyy)
Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) 
--------------------------- 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 ... [ 15 15 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.96, Number of G-points: 313, Lambda: 0.30
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/135) =======================
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: 135
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
   2.917   (   0.000   -0.000    0.000)    0.000
   2.917   (   0.000    0.000   -0.000)    0.000
   9.918   (   0.000    0.000    0.000)    0.000
  11.785   (   0.000   -0.000    0.000)    0.000
  12.577   (   0.000    0.000   -0.000)    0.000
  12.577   (   0.000    0.000    0.000)    0.000
  12.772   (   0.000   -0.000    0.000)    0.000
  13.026   (   0.000    0.000    0.000)    0.000
  13.026   (   0.000   -0.000   -0.000)    0.000
======================= Grid point 1 (2/135) =======================
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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.749   (  33.427   19.299    0.000)   38.599
   0.758   (  32.387   18.699    0.000)   37.397
   1.606   (  68.249   39.404    0.000)   78.808
   3.014   (   8.143    4.702    0.000)    9.403
   3.339   (  33.350   19.254    0.000)   38.509
   9.797   (  -9.998   -5.772    0.000)   11.545
  11.721   (  -5.368   -3.099    0.000)    6.198
  12.537   (  -3.150   -1.819    0.000)    3.637
  12.635   (   4.226    2.440    0.000)    4.880
  12.745   (  -1.930   -1.114    0.000)    2.229
  13.010   (  -1.593   -0.920    0.000)    1.839
  14.517   (  -2.901   -1.675    0.000)    3.350
======================= Grid point 2 (3/135) =======================
q-point: ( 0.13  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.475   (  29.457   17.007    0.000)   34.014
   1.507   (  32.519   18.775    0.000)   37.550
   3.131   (  63.738   36.799    0.000)   73.598
   3.267   (  13.162    7.599    0.000)   15.198
   4.259   (  42.877   24.755    0.000)   49.510
   9.502   ( -14.338   -8.278    0.000)   16.556
  11.553   (  -8.767   -5.062    0.000)   10.124
  12.448   (  -4.327   -2.498    0.000)    4.996
  12.634   (  -4.282   -2.472    0.000)    4.945
  12.792   (   5.472    3.159    0.000)    6.319
  12.945   (  -4.078   -2.355    0.000)    4.709
  14.416   (  -5.879   -3.394    0.000)    6.789
======================= Grid point 3 (4/135) =======================
q-point: ( 0.20  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.092   (  23.694   13.680    0.000)   27.359
   2.244   (  31.460   18.164    0.000)   36.327
   3.584   (  13.545    7.820    0.000)   15.641
   4.511   (  55.064   31.791    0.000)   63.582
   5.181   (  35.420   20.450    0.000)   40.900
   9.219   (  -8.118   -4.687    0.000)    9.374
  11.332   ( -10.202   -5.890    0.000)   11.780
  12.349   (  -4.141   -2.391    0.000)    4.781
  12.518   (  -5.379   -3.105    0.000)    6.211
  12.829   (  -5.661   -3.268    0.000)    6.537
  12.896   (   2.902    1.675    0.000)    3.351
  14.242   (  -9.437   -5.448    0.000)   10.896
======================= Grid point 4 (5/135) =======================
q-point: ( 0.27  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.548   (  15.703    9.066    0.000)   18.132
   2.949   (  29.659   17.124    0.000)   34.247
   3.856   (   9.512    5.492    0.000)   10.984
   5.494   (  23.240   13.418    0.000)   26.836
   5.969   (  39.451   22.777    0.000)   45.555
   9.214   (   8.562    4.943    0.000)    9.886
  11.091   ( -10.636   -6.141    0.000)   12.281
  12.261   (  -3.442   -1.987    0.000)    3.974
  12.402   (  -4.463   -2.576    0.000)    5.153
  12.704   (  -4.807   -2.775    0.000)    5.550
  12.913   (  -1.428   -0.825    0.000)    1.649
  13.973   ( -14.155   -8.172    0.000)   16.344
======================= Grid point 5 (6/135) =======================
q-point: ( 0.33  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.815   (   7.815    4.512    0.000)    9.024
   3.599   (  26.467   15.281    0.000)   30.562
   4.003   (   3.157    1.823    0.000)    3.645
   5.687   (  -0.987   -0.570    0.000)    1.140
   6.930   (  40.418   23.335    0.000)   46.670
   9.566   (  19.166   11.066    0.000)   22.131
  10.861   (  -8.271   -4.775    0.000)    9.550
  12.194   (  -2.304   -1.330    0.000)    2.661
  12.321   (  -2.397   -1.384    0.000)    2.768
  12.621   (  -2.335   -1.348    0.000)    2.696
  12.842   (  -4.405   -2.543    0.000)    5.087
  13.584   ( -19.682  -11.363    0.000)   22.726
======================= Grid point 6 (7/135) =======================
q-point: ( 0.40  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.929   (   2.681    1.548    0.000)    3.096
   4.013   (  -1.652   -0.954    0.000)    1.908
   4.149   (  20.771   11.992    0.000)   23.985
   5.582   (  -6.580   -3.799    0.000)    7.598
   7.771   (  31.959   18.451    0.000)   36.903
   9.879   (   2.865    1.654    0.000)    3.308
  10.844   (  10.600    6.120    0.000)   12.240
  12.156   (  -0.971   -0.561    0.000)    1.121
  12.291   (  -0.396   -0.229    0.000)    0.457
  12.592   (  -0.387   -0.223    0.000)    0.446
  12.732   (  -4.650   -2.684    0.000)    5.369
  13.078   ( -23.620  -13.637    0.000)   27.274
======================= Grid point 7 (8/135) =======================
q-point: ( 0.47  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.963   (   0.568    0.328    0.000)    0.656
   3.966   (  -1.580   -0.912    0.000)    1.824
   4.513   (   9.360    5.404    0.000)   10.808
   5.439   (  -4.459   -2.575    0.000)    5.149
   8.344   (  15.348    8.861    0.000)   17.722
   9.760   (  -7.170   -4.139    0.000)    8.279
  11.238   (  16.051    9.267    0.000)   18.534
  12.144   (  -0.185   -0.107    0.000)    0.213
  12.289   (  -0.200   -0.115    0.000)    0.231
  12.571   ( -15.980   -9.226    0.000)   18.452
  12.591   (   0.110    0.063    0.000)    0.127
  12.652   (  -2.393   -1.381    0.000)    2.763
======================= Grid point 16 (9/135) =======================
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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.305   (  18.852   32.652    0.000)   37.704
   1.369   (  20.717   35.883    0.000)   41.435
   2.677   (  34.477   59.716    0.000)   68.955
   3.282   (  11.285   19.545    0.000)   22.569
   3.929   (  22.454   38.892    0.000)   44.909
   9.589   (  -8.116  -14.057    0.000)   16.231
  11.605   (  -4.705   -8.149    0.000)    9.409
  12.486   (  -2.506   -4.340    0.000)    5.011
  12.692   (   2.215    3.837    0.000)    4.431
  12.697   (  -1.781   -3.085    0.000)    3.562
  12.997   (  -0.564   -0.977    0.000)    1.128
  14.440   (  -3.477   -6.022    0.000)    6.954
======================= Grid point 17 (10/135) =======================
q-point: ( 0.13  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.982   (  24.035   27.492    0.000)   36.517
   2.060   (  15.649   37.292    0.000)   40.442
   3.548   (  13.301   22.014    0.000)   25.721
   4.029   (  39.884   47.129    0.000)   61.740
   4.768   (  33.786   25.209    0.000)   42.154
   9.298   (  -9.014  -10.760    0.000)   14.036
  11.411   (  -7.462   -8.941    0.000)   11.646
  12.383   (  -3.708   -4.261    0.000)    5.648
  12.603   (  -5.800   -0.807    0.000)    5.856
  12.784   (   4.959   -1.510    0.000)    5.184
  12.954   (  -4.259    0.829    0.000)    4.338
  14.281   (  -5.557  -10.128    0.000)   11.553
======================= Grid point 18 (11/135) =======================
q-point: ( 0.20  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.616   (   7.304   35.205    0.000)   35.955
   2.671   (  24.611   24.522    0.000)   34.743
   3.912   (   3.818   23.772    0.000)   24.076
   5.145   (  37.197   29.304    0.000)   47.353
   5.527   (  33.527   15.159    0.000)   36.794
   9.149   (   1.107   -1.415    0.000)    1.797
  11.179   (  -8.670   -9.417    0.000)   12.801
  12.287   (  -2.811   -3.473    0.000)    4.468
  12.489   (  -6.640   -0.296    0.000)    6.647
  12.788   (  -3.906   -0.150    0.000)    3.908
  12.898   (  -0.724   -2.393    0.000)    2.501
  14.044   (  -8.013  -14.265    0.000)   16.361
======================= Grid point 19 (12/135) =======================
q-point: ( 0.27  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.988   (  -3.060   32.507    0.000)   32.651
   3.327   (  24.863   20.333    0.000)   32.119
   4.146   (  -3.702   22.413    0.000)   22.717
   5.640   (   7.569    2.378    0.000)    7.933
   6.422   (  43.606   21.960    0.000)   48.824
   9.338   (  16.270    6.987    0.000)   17.706
  10.932   (  -8.409   -9.846    0.000)   12.948
  12.221   (  -1.776   -1.645    0.000)    2.421
  12.377   (  -5.682   -0.377    0.000)    5.694
  12.699   (  -6.058    1.774    0.000)    6.312
  12.849   (  -1.032   -5.311    0.000)    5.410
  13.712   ( -12.220  -17.565    0.000)   21.398
======================= Grid point 20 (13/135) =======================
q-point: ( 0.33  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.173   ( -10.364   29.163    0.000)   30.950
   3.919   (  21.868   16.682    0.000)   27.504
   4.231   ( -10.585   20.077    0.000)   22.697
   5.643   (  -3.938   -3.099    0.000)    5.011
   7.336   (  38.477   17.379    0.000)   42.220
   9.727   (  16.367    3.820    0.000)   16.807
  10.749   (   0.385   -5.874    0.000)    5.887
  12.191   (  -1.390    0.682    0.000)    1.549
  12.295   (  -2.635   -1.248    0.000)    2.916
  12.631   (  -3.498    0.984    0.000)    3.634
  12.762   (  -2.925   -4.337    0.000)    5.232
  13.272   ( -17.320  -19.088    0.000)   25.775
======================= Grid point 21 (14/135) =======================
q-point: ( 0.40  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.238   ( -13.982   27.265    0.000)   30.641
   4.209   ( -13.823   19.585    0.000)   23.972
   4.387   (  14.658   11.859    0.000)   18.855
   5.507   (  -6.561   -3.484    0.000)    7.429
   8.063   (  27.867   10.906    0.000)   29.925
   9.786   (  -3.733  -10.229    0.000)   10.889
  10.961   (  21.626    4.439    0.000)   22.077
  12.181   (  -1.542    2.180    0.000)    2.670
  12.254   (  -0.050   -3.939    0.000)    3.939
  12.585   (  -1.738   -2.469    0.000)    3.019
  12.678   (  -3.368   -2.942    0.000)    4.472
  12.779   ( -17.239  -12.337    0.000)   21.199
======================= Grid point 22 (15/135) =======================
q-point: ( 0.47  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.251   ( -15.286   26.477    0.000)   30.573
   4.181   ( -12.218   21.163    0.000)   24.436
   4.593   (  -1.520    2.633    0.000)    3.040
   5.414   (  -0.119    0.206    0.000)    0.238
   8.398   (   2.079   -3.600    0.000)    4.157
   9.644   (   3.909   -6.770    0.000)    7.818
  11.264   (   4.232   -7.331    0.000)    8.465
  12.184   (  -2.147    3.718    0.000)    4.293
  12.191   (   4.958   -8.587    0.000)    9.916
  12.472   (  -0.132    0.229    0.000)    0.265
  12.584   (   1.460   -2.529    0.000)    2.920
  12.668   (  -2.210    3.828    0.000)    4.420
======================= Grid point 33 (16/135) =======================
q-point: ( 0.13  0.13  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.561   (  17.312   29.985    0.000)   34.624
   2.806   (  20.134   34.872    0.000)   40.267
   4.029   (  13.941   24.146    0.000)   27.881
   4.975   (  27.009   46.781    0.000)   54.018
   5.235   (  12.262   21.239    0.000)   24.525
   9.135   (  -2.301   -3.985    0.000)    4.601
  11.206   (  -6.609  -11.446    0.000)   13.217
  12.293   (  -2.457   -4.256    0.000)    4.914
  12.568   (  -1.844   -3.194    0.000)    3.689
  12.781   (   0.521    0.902    0.000)    1.041
  12.926   (  -2.124   -3.678    0.000)    4.248
  14.027   (  -8.654  -14.989    0.000)   17.308
======================= Grid point 34 (17/135) =======================
q-point: ( 0.20  0.13  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.158   (  19.485   25.410    0.000)   32.021
   3.390   (   4.990   37.770    0.000)   38.099
   4.445   (   3.256   27.394    0.000)   27.586
   5.531   (   7.913    8.467    0.000)   11.589
   5.935   (  37.292   27.863    0.000)   46.552
   9.184   (   7.739    5.137    0.000)    9.289
  10.960   (  -7.930  -12.439    0.000)   14.751
  12.224   (  -1.461   -2.683    0.000)    3.055
  12.483   (  -5.398   -0.571    0.000)    5.428
  12.789   (   0.692   -4.432    0.000)    4.485
  12.824   (  -5.392   -0.496    0.000)    5.415
  13.695   (  -9.406  -19.999    0.000)   22.101
======================= Grid point 35 (18/135) =======================
q-point: ( 0.27  0.13  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.684   (   2.314   29.749    0.000)   29.839
   3.777   (   8.715   29.261    0.000)   30.532
   4.673   (  -8.692   28.302    0.000)   29.607
   5.641   (   1.420   -1.313    0.000)    1.934
   6.835   (  38.530   19.847    0.000)   43.341
   9.471   (  16.344    5.765    0.000)   17.331
  10.706   (  -5.596  -12.261    0.000)   13.477
  12.188   (  -0.170   -2.258    0.000)    2.264
  12.382   (  -6.123    0.882    0.000)    6.186
  12.679   (  -3.089   -6.939    0.000)    7.596
  12.786   (  -3.103    2.057    0.000)    3.723
  13.305   ( -11.307  -22.184    0.000)   24.900
======================= Grid point 36 (19/135) =======================
q-point: ( 0.33  0.13  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.853   ( -15.964   35.696    0.000)   39.103
   4.264   (  15.844   18.250    0.000)   24.168
   4.733   ( -15.982   28.101    0.000)   32.328
   5.579   (  -5.119   -3.191    0.000)    6.032
   7.634   (  34.107   13.088    0.000)   36.532
   9.684   (   4.189   -8.364    0.000)    9.354
  10.656   (  12.880   -2.879    0.000)   13.198
  12.155   (   0.203   -6.145    0.000)    6.149
  12.303   (  -4.583    2.276    0.000)    5.117
  12.574   (  -1.507   -6.107    0.000)    6.290
  12.730   (  -5.554   -1.761    0.000)    5.826
  12.896   ( -11.466  -14.292    0.000)   18.324
======================= Grid point 37 (20/135) =======================
q-point: ( 0.40  0.13  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.894   ( -19.574   35.427    0.000)   40.474
   4.623   (   5.987   11.993    0.000)   13.405
   4.718   ( -17.781   29.038    0.000)   34.049
   5.446   (  -4.798   -2.577    0.000)    5.446
   8.209   (  19.004    3.944    0.000)   19.409
   9.514   (  -2.299  -15.832    0.000)   15.998
  10.989   (  19.539   -1.972    0.000)   19.639
  12.026   (  -0.908  -16.130    0.000)   16.155
  12.269   (  -2.810    3.436    0.000)    4.438
  12.506   (  -1.425   -3.143    0.000)    3.451
  12.534   (  -3.702   -5.368    0.000)    6.520
  12.761   (  -3.527    3.321    0.000)    4.845
======================= Grid point 49 (21/135) =======================
q-point: ( 0.20  0.20  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.665   (  14.539   25.182    0.000)   29.078
   4.037   (  14.923   25.847    0.000)   29.846
   4.928   (  11.311   19.592    0.000)   22.623
   5.611   (   0.769    1.332    0.000)    1.538
   6.555   (  19.234   33.314    0.000)   38.468
   9.348   (   6.119   10.598    0.000)   12.238
  10.694   (  -7.862  -13.617    0.000)   15.724
  12.176   (  -1.428   -2.474    0.000)    2.857
  12.459   (  -1.173   -2.031    0.000)    2.345
  12.647   (  -4.324   -7.490    0.000)    8.649
  12.846   (   0.416    0.721    0.000)    0.833
  13.276   ( -12.300  -21.304    0.000)   24.600
======================= Grid point 50 (22/135) =======================
q-point: ( 0.27  0.20  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.146   (  13.346   21.387    0.000)   25.209
   4.418   (  -0.747   28.308    0.000)   28.317
   5.196   (  -4.359   22.603    0.000)   23.019
   5.604   (  -1.065   -2.461    0.000)    2.682
   7.240   (  27.487   21.522    0.000)   34.911
   9.532   (   5.474   -1.106    0.000)    5.584
  10.490   (   0.423   -7.216    0.000)    7.229
  12.104   (  -2.418   -7.147    0.000)    7.545
  12.402   (  -4.328    0.635    0.000)    4.374
  12.523   (  -1.102   -7.546    0.000)    7.626
  12.821   (  -7.117   -2.969    0.000)    7.711
  12.897   (  -5.175  -13.538    0.000)   14.494
======================= Grid point 51 (23/135) =======================
q-point: ( 0.33  0.20  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.468   (  -4.705   23.984    0.000)   24.441
   4.677   (   3.523   23.287    0.000)   23.552
   5.263   ( -12.644   22.501    0.000)   25.810
   5.516   (  -5.716   -2.649    0.000)    6.300
   7.886   (  25.444   13.058    0.000)   28.599
   9.412   (  -5.607  -17.930    0.000)   18.787
  10.652   (  17.414    2.830    0.000)   17.642
  11.907   (  -5.553  -18.557    0.000)   19.370
  12.340   (  -3.100    0.546    0.000)    3.148
  12.451   (   0.250   -5.585    0.000)    5.590
  12.597   (  -5.992   -4.137    0.000)    7.282
  12.837   (  -2.250    1.522    0.000)    2.717
======================= Grid point 52 (24/135) =======================
q-point: ( 0.40  0.20  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.531   ( -14.795   25.626    0.000)   29.590
   4.874   (  -7.631   13.218    0.000)   15.262
   5.250   ( -11.003   19.057    0.000)   22.006
   5.427   (  -2.311    4.003    0.000)    4.622
   8.230   (   0.819   -1.418    0.000)    1.637
   9.205   (   7.991  -13.841    0.000)   15.983
  10.896   (   3.696   -6.401    0.000)    7.391
  11.694   (   8.585  -14.870    0.000)   17.170
  12.315   (  -0.690    1.196    0.000)    1.381
  12.425   (   2.300   -3.984    0.000)    4.600
  12.534   (  -1.669    2.890    0.000)    3.337
  12.825   (  -1.560    2.702    0.000)    3.120
======================= Grid point 66 (25/135) =======================
q-point: ( 0.27  0.27  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: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.557   (  11.158   19.326    0.000)   22.316
   4.847   (   8.763   15.179    0.000)   17.527
   5.528   (  -3.087   -5.348    0.000)    6.175
   5.537   (   6.352   11.002    0.000)   12.704
   7.707   (  14.476   25.073    0.000)   28.952
   9.367   (  -8.853  -15.334    0.000)   17.707
  10.507   (   5.387    9.330    0.000)   10.773
  11.860   ( -10.006  -17.332    0.000)   20.013
  12.389   (  -1.362   -2.359    0.000)    2.723
  12.400   (  -2.349   -4.069    0.000)    4.698
  12.631   (  -3.185   -5.517    0.000)    6.371
  12.865   (   0.070    0.122    0.000)    0.141
======================= Grid point 67 (26/135) =======================
q-point: ( 0.33  0.27  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: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.862   (   0.303   17.088    0.000)   17.090
   5.081   (   3.229   14.639    0.000)   14.991
   5.369   (  -6.972   -6.894    0.000)    9.804
   5.668   (  -3.769   11.907    0.000)   12.490
   8.184   (  12.347   17.208    0.000)   21.179
   9.008   (  -6.544  -21.836    0.000)   22.796
  10.784   (  10.324   10.697    0.000)   14.866
  11.485   (  -6.386  -21.982    0.000)   22.891
  12.333   (  -1.448   -0.596    0.000)    1.566
  12.364   (   0.778   -3.039    0.000)    3.137
  12.586   (  -1.222    1.121    0.000)    1.658
  12.862   (  -0.673    0.755    0.000)    1.012
======================= Grid point 82 (27/135) =======================
q-point: ( 0.33  0.33  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: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   5.157   (  -5.777  -10.006    0.000)   11.554
   5.157   (   5.777   10.006    0.000)   11.554
   5.246   (  -0.000   -0.000    0.000)    0.000
   5.789   (  -0.000   -0.000    0.000)    0.000
   8.571   ( -10.612  -18.380    0.000)   21.223
   8.571   (  10.612   18.380    0.000)   21.223
  11.076   (  -9.109  -15.778    0.000)   18.218
  11.076   (   9.109   15.778    0.000)   18.218
  12.330   (  -0.253   -0.438    0.000)    0.505
  12.330   (   0.253    0.438    0.000)    0.505
  12.601   (  -0.000   -0.000    0.000)    0.000
  12.869   (  -0.000   -0.000    0.000)    0.000
======================= Grid point 241 (28/135) =======================
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: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.738   (   0.000   -0.000   36.406)   36.406
   0.738   (   0.000    0.000   36.406)   36.406
   1.679   (   0.000    0.000   84.319)   84.319
   2.903   (  -0.000   -0.000   -1.551)    1.551
   2.903   (   0.000   -0.000   -1.551)    1.551
   9.518   (   0.000   -0.000  -35.914)   35.914
  12.077   (  -0.000    0.000   24.940)   24.940
  12.587   (   0.000   -0.000    1.045)    1.045
  12.587   (   0.000   -0.000    1.045)    1.045
  13.003   (   0.000   -0.000   -2.282)    2.282
  13.003   (  -0.000   -0.000   -2.282)    2.282
  14.435   (   0.000    0.000   -5.405)    5.405
======================= Grid point 242 (29/135) =======================
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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.058   (  23.260   13.429   25.461)   37.009
   1.348   (  22.937   13.243   39.765)   47.778
   2.114   (  46.219   26.685   50.980)   73.806
   3.002   (   8.306    4.795   -1.336)    9.683
   3.314   (  32.856   18.969   -2.108)   37.997
   9.412   (  -8.798   -5.080  -34.792)   36.245
  11.990   (  -6.861   -3.961   22.711)   24.053
  12.548   (  -3.130   -1.807    1.086)    3.774
  12.650   (   4.617    2.666    2.476)    5.879
  12.878   (  -6.256   -3.612    5.631)    9.159
  12.986   (  -1.673   -0.966   -2.393)    3.075
  14.422   (  -1.925   -1.112   -7.746)    8.059
======================= Grid point 243 (30/135) =======================
q-point: ( 0.13  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.651   (  26.396   15.240   16.487)   34.653
   1.835   (  22.618   13.059   31.915)   41.240
   3.259   (  13.273    7.663   -0.987)   15.358
   3.403   (  58.936   34.027   25.520)   72.681
   4.226   (  42.557   24.570   -2.990)   49.231
   9.153   ( -12.390   -7.153  -31.743)   34.818
  11.792   (  -9.848   -5.686   20.150)   23.137
  12.458   (  -4.365   -2.520    1.079)    5.155
  12.683   (  -4.132   -2.385    4.441)    6.518
  12.835   (   3.763    2.173    3.704)    5.710
  12.919   (  -4.128   -2.383   -2.559)    5.410
  14.339   (  -5.230   -3.020   -7.405)    9.555
======================= Grid point 244 (31/135) =======================
q-point: ( 0.20  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.221   (  22.331   12.893   12.354)   28.593
   2.391   (  25.201   14.550   16.050)   33.232
   3.577   (  13.544    7.820   -0.851)   15.662
   4.694   (  50.378   29.086   15.734)   60.262
   5.151   (  37.062   21.398   -1.585)   42.825
   8.923   (  -5.457   -3.150  -26.852)   27.582
  11.556   ( -10.406   -6.008   19.079)   22.547
  12.357   (  -4.238   -2.447    0.916)    4.979
  12.556   (  -5.878   -3.393    3.623)    7.693
  12.803   (  -5.649   -3.262   -2.593)    7.020
  12.921   (   2.343    1.353    2.517)    3.695
  14.179   (  -8.863   -5.117   -6.211)   11.971
======================= Grid point 245 (32/135) =======================
q-point: ( 0.27  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.654   (  14.993    8.656   10.250)   20.119
   2.975   (  25.165   14.529    4.076)   29.342
   3.849   (   9.501    5.485   -0.867)   11.005
   5.454   (  14.197    8.197   -6.555)   17.655
   6.093   (  45.629   26.344   13.987)   54.513
   8.986   (  11.530    6.657  -20.041)   24.060
  11.322   (  -9.796   -5.656   19.744)   22.755
  12.267   (  -3.533   -2.040    0.689)    4.137
  12.431   (  -4.751   -2.743    2.826)    6.171
  12.678   (  -4.798   -2.770   -2.550)    6.099
  12.923   (  -2.138   -1.235    1.141)    2.720
  13.923   ( -13.576   -7.838   -4.896)   16.424
======================= Grid point 246 (33/135) =======================
q-point: ( 0.33  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.909   (   7.431    4.290    9.102)   12.508
   3.534   (  23.096   13.334   -4.960)   27.126
   3.997   (   3.242    1.872   -0.785)    3.825
   5.537   (  -3.701   -2.137  -14.884)   15.485
   7.108   (  41.268   23.826   17.050)   50.610
   9.397   (  21.438   12.377  -13.765)   28.324
  11.117   (  -7.278   -4.202   21.572)   23.151
  12.198   (  -2.356   -1.360    0.529)    2.771
  12.346   (  -2.613   -1.509    2.357)    3.829
  12.595   (  -2.381   -1.375   -2.588)    3.776
  12.833   (  -5.288   -3.053   -0.775)    6.155
  13.548   ( -18.911  -10.918   -3.349)   22.092
======================= Grid point 247 (34/135) =======================
q-point: ( 0.40  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.016   (   2.465    1.423    8.461)    8.927
   4.010   (  -1.469   -0.848   -0.461)    1.758
   4.020   (  18.563   10.717  -11.461)   24.306
   5.385   (  -8.081   -4.666  -18.858)   21.040
   7.964   (  32.532   18.782   18.269)   41.771
   9.804   (   9.871    5.699   -5.722)   12.754
  11.067   (   5.590    3.228   19.142)   20.201
  12.160   (  -1.001   -0.578    0.444)    1.238
  12.309   (  -0.845   -0.488    1.805)    2.052
  12.564   (  -0.457   -0.264   -2.731)    2.781
  12.698   (  -6.301   -3.638   -3.531)    8.088
  13.069   ( -21.673  -12.513   -0.552)   25.032
======================= Grid point 248 (35/135) =======================
q-point: ( 0.47  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.046   (   0.484    0.279    8.119)    8.138
   3.967   (  -1.470   -0.849   -0.086)    1.700
   4.349   (   8.649    4.994  -14.925)   17.959
   5.218   (  -5.066   -2.925  -20.879)   21.683
   8.548   (  15.661    9.042   19.208)   26.381
   9.820   (  -4.195   -2.422    5.357)    7.222
  11.360   (  14.240    8.221   11.204)   19.897
  12.148   (  -0.197   -0.114    0.397)    0.458
  12.274   (  -3.551   -2.050   -1.885)    4.513
  12.506   (  -9.300   -5.369   -4.806)   11.765
  12.563   (   0.077    0.044   -2.856)    2.858
  12.701   (  -6.315   -3.646    3.778)    8.213
======================= Grid point 257 (36/135) =======================
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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.557   (  18.273   31.650   17.485)   40.514
   1.687   (  12.127   21.005   35.642)   43.111
   2.969   (  27.687   47.955   25.090)   60.793
   3.287   (  14.710   25.479    3.204)   29.594
   3.906   (  22.150   38.365   -2.376)   44.364
   9.229   (  -7.096  -12.290  -32.725)   35.670
  11.850   (  -5.444   -9.429   20.662)   23.355
  12.490   (  -2.438   -4.222    0.761)    4.934
  12.724   (   2.369    4.103    3.607)    5.955
  12.772   (  -3.187   -5.520    5.107)    8.167
  12.972   (  -0.595   -1.031   -2.535)    2.801
  14.360   (  -3.019   -5.230   -7.486)    9.618
======================= Grid point 258 (37/135) =======================
q-point: ( 0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.150   (  15.891   25.629   16.604)   34.424
   2.212   (  17.982   29.882   16.231)   38.467
   3.554   (  10.535   21.589   -0.233)   24.023
   4.243   (  39.821   46.116   20.422)   64.261
   4.730   (  33.610   24.596   -3.433)   41.790
   8.981   (  -7.243   -8.812  -28.916)   31.084
  11.638   (  -7.810   -9.298   19.344)   22.839
  12.391   (  -3.576   -4.311    0.865)    5.667
  12.643   (  -6.318   -1.825    3.444)    7.424
  12.808   (   3.337   -1.733    2.164)    4.338
  12.934   (  -3.079    0.673   -1.602)    3.535
  14.216   (  -5.229   -9.420   -6.410)   12.537
======================= Grid point 259 (38/135) =======================
q-point: ( 0.20  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.670   (  13.302   25.545    6.767)   29.585
   2.782   (  14.681   28.000   10.960)   33.461
   3.906   (   3.689   23.822   -0.857)   24.121
   5.221   (  25.246   19.422    1.548)   31.890
   5.571   (  42.523   23.573    9.447)   49.529
   8.890   (   3.697    0.631  -23.182)   23.483
  11.407   (  -8.374   -8.766   19.575)   23.025
  12.295   (  -2.817   -3.456    0.871)    4.543
  12.515   (  -6.730   -1.138    2.347)    7.217
  12.772   (  -5.605    0.296   -1.869)    5.916
  12.905   (   0.423   -3.371    1.152)    3.587
  13.990   (  -7.694  -13.653   -5.228)   16.521
======================= Grid point 260 (39/135) =======================
q-point: ( 0.27  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.056   (  -0.446   28.765    6.138)   29.416
   3.311   (  19.351   19.533    0.497)   27.500
   4.139   (  -3.690   22.399   -0.890)   22.719
   5.525   (   4.011   -0.246  -11.989)   12.644
   6.594   (  44.483   23.610   16.862)   53.108
   9.142   (  18.688    8.551  -16.548)   26.386
  11.180   (  -7.536   -8.605   21.230)   24.116
  12.229   (  -1.919   -1.484    0.824)    2.562
  12.395   (  -5.500   -1.115    1.785)    5.889
  12.672   (  -6.111    1.485   -2.646)    6.823
  12.849   (  -1.671   -5.726    0.154)    5.967
  13.671   ( -11.767  -17.032   -3.982)   21.081
======================= Grid point 261 (40/135) =======================
q-point: ( 0.33  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.244   (  -9.223   27.212    6.733)   29.511
   3.819   (  18.749   15.202   -8.294)   25.523
   4.225   ( -10.403   20.016   -0.768)   22.571
   5.466   (  -5.887   -4.299  -17.193)   18.674
   7.527   (  38.794   18.171   18.064)   46.492
   9.593   (  20.108    6.441  -10.352)   23.515
  11.010   (  -1.083   -6.413   21.963)   22.905
  12.197   (  -1.681    0.861    0.620)    1.988
  12.311   (  -2.395   -1.896    1.660)    3.476
  12.600   (  -3.495    0.356   -3.086)    4.676
  12.746   (  -3.853   -4.347   -1.650)    6.038
  13.247   ( -16.392  -18.547   -2.267)   24.856
======================= Grid point 262 (41/135) =======================
q-point: ( 0.40  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.309   ( -13.217   25.846    6.782)   29.811
   4.203   ( -11.383   18.628   -1.746)   21.900
   4.238   (  11.008   11.697  -12.395)   20.289
   5.295   (  -7.617   -4.154  -20.106)   21.899
   8.266   (  28.113   11.479   19.047)   35.845
   9.799   (   2.233   -5.893    1.212)    6.417
  11.121   (  16.799    1.450   14.545)   22.269
  12.183   (  -1.640    1.865    0.239)    2.495
  12.260   (  -1.096   -4.889    0.623)    5.049
  12.537   (  -2.880   -4.334   -4.627)    6.963
  12.651   (  -4.657    0.624   -2.913)    5.528
  12.791   ( -13.766  -13.210    1.479)   19.137
======================= Grid point 263 (42/135) =======================
q-point: ( 0.47  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.321   ( -14.566   25.228    6.747)   29.902
   4.179   ( -12.030   20.836   -0.354)   24.062
   4.424   (  -1.461    2.531  -15.393)   15.668
   5.189   (  -0.089    0.154  -21.193)   21.194
   8.607   (   1.923   -3.331   19.598)   19.972
   9.728   (   3.508   -6.076    7.455)   10.238
  11.375   (   4.408   -7.635   10.340)   13.588
  12.152   (   4.608   -7.981   -4.151)   10.108
  12.187   (  -2.090    3.620    0.385)    4.198
  12.450   (  -0.091    0.158   -1.843)    1.852
  12.597   (   0.019   -0.034    0.143)    0.148
  12.650   (  -0.446    0.772   -0.658)    1.108
======================= Grid point 274 (43/135) =======================
q-point: ( 0.13  0.13  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.624   (  14.095   24.413    8.219)   29.363
   2.903   (  19.419   33.634    9.458)   39.972
   4.030   (  13.816   23.930   -0.397)   27.635
   5.165   (  11.124   19.267   -7.024)   23.330
   5.184   (  26.914   46.616   19.909)   57.391
   8.864   (  -0.815   -1.411  -24.468)   24.522
  11.435   (  -6.261  -10.845   19.750)   23.385
  12.301   (  -2.439   -4.224    0.835)    4.948
  12.588   (  -2.242   -3.883    1.702)    4.796
  12.800   (   0.234    0.405    1.756)    1.818
  12.901   (  -2.280   -3.949   -2.081)    5.013
  13.975   (  -8.305  -14.385   -5.131)   17.385
======================= Grid point 275 (44/135) =======================
q-point: ( 0.20  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.118   (  15.943   21.720   -1.861)   27.007
   3.476   (   6.077   35.900    8.095)   37.300
   4.441   (   2.898   27.619   -0.588)   27.777
   5.429   (   6.751    3.974  -10.835)   13.370
   6.129   (  36.350   30.830   18.778)   51.229
   8.967   (   9.746    7.263  -18.785)   22.374
  11.209   (  -7.221  -11.047   21.524)   25.248
  12.233   (  -1.380   -2.622    0.938)    3.108
  12.494   (  -5.480   -1.215    0.835)    5.675
  12.787   (  -4.331   -5.348   -2.041)    7.179
  12.802   (  -0.911   -0.465   -0.138)    1.032
  13.654   (  -9.157  -19.416   -4.022)   21.841
======================= Grid point 276 (45/135) =======================
q-point: ( 0.27  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.581   (  11.527   20.991   -8.759)   25.499
   3.842   (  -2.049   33.436    6.434)   34.111
   4.665   (  -8.906   28.239   -0.991)   29.626
   5.484   (  -0.292   -3.078  -15.398)   15.705
   7.036   (  38.160   21.204   19.005)   47.613
   9.309   (  18.683    7.825  -12.982)   24.059
  10.978   (  -5.496  -11.470   23.248)   26.499
  12.198   (  -0.134   -2.197    1.120)    2.470
  12.386   (  -6.023    0.225    0.367)    6.038
  12.649   (  -3.819   -6.527   -3.037)    8.149
  12.775   (  -2.874    1.004   -1.037)    3.217
  13.274   ( -10.943  -21.687   -2.932)   24.468
======================= Grid point 277 (46/135) =======================
q-point: ( 0.33  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.854   (  -9.082   29.332   -2.357)   30.796
   4.163   (   8.297   21.051   -6.067)   23.427
   4.721   ( -15.775   27.469   -1.478)   31.711
   5.383   (  -6.592   -3.743  -18.494)   19.987
   7.843   (  33.972   14.099   19.558)   41.658
   9.633   (  10.975   -2.603   -3.937)   11.947
  10.874   (   7.480   -6.659   19.476)   21.900
  12.163   (  -0.198   -6.581    0.797)    6.632
  12.303   (  -4.475    1.594   -0.026)    4.751
  12.535   (  -2.224   -5.726   -3.848)    7.248
  12.722   (  -4.516   -0.516   -0.925)    4.639
  12.873   ( -11.319  -15.271   -1.915)   19.105
======================= Grid point 278 (47/135) =======================
q-point: ( 0.40  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.921   ( -17.052   31.904    1.549)   36.209
   4.467   (   4.101   13.211  -13.030)   19.004
   4.703   ( -17.205   27.611   -2.062)   32.598
   5.228   (  -5.781   -2.105  -19.993)   20.918
   8.425   (  18.869    4.641   20.100)   27.957
   9.603   (  -0.180  -12.607    7.748)   14.798
  11.099   (  17.895   -3.915   10.604)   21.166
  12.010   (  -2.315  -16.395   -1.854)   16.661
  12.266   (  -2.514    2.785   -0.328)    3.767
  12.471   (  -0.925   -0.592   -2.573)    2.798
  12.555   (  -3.231   -7.592    1.489)    8.384
  12.730   (  -3.580    3.508   -2.964)    5.823
======================= Grid point 290 (48/135) =======================
q-point: ( 0.20  0.20  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.551   (  12.496   21.644   -9.333)   26.678
   4.097   (  14.450   25.027    5.786)   29.473
   4.928   (  11.337   19.637   -0.187)   22.676
   5.452   (  -0.429   -0.744  -15.515)   15.539
   6.775   (  19.560   33.879   20.437)   44.137
   9.175   (   7.567   13.106  -14.097)   20.682
  10.969   (  -7.327  -12.691   23.730)   27.890
  12.186   (  -1.417   -2.455    1.021)    3.013
  12.459   (  -1.411   -2.444   -0.194)    2.829
  12.627   (  -4.632   -8.022   -1.963)    9.469
  12.825   (   0.329    0.569   -1.917)    2.027
  13.246   ( -12.006  -20.794   -2.950)   24.192
======================= Grid point 291 (49/135) =======================
q-point: ( 0.27  0.20  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.969   (  11.929   18.834  -15.416)   27.104
   4.467   (  -0.339   26.621    4.360)   26.978
   5.184   (  -5.997   22.004   -2.094)   22.903
   5.420   (  -0.685   -3.032  -16.565)   16.854
   7.466   (  27.114   22.571   20.949)   41.030
   9.440   (   9.976    4.577   -7.140)   13.094
  10.747   (  -2.491  -10.461   22.467)   24.908
  12.112   (  -2.566   -7.503    0.744)    7.965
  12.395   (  -4.464    0.144   -0.764)    4.531
  12.496   (  -1.281   -7.318   -2.624)    7.880
  12.803   (  -6.034   -2.217   -1.701)    6.650
  12.872   (  -5.726  -13.958   -2.322)   15.264
======================= Grid point 292 (50/135) =======================
q-point: ( 0.33  0.20  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.301   (   3.950   17.617  -17.282)   24.993
   4.673   (  -4.349   25.915    1.545)   26.323
   5.179   ( -11.693    9.954  -13.264)   20.292
   5.375   (  -8.069    9.594   -7.513)   14.615
   8.115   (  25.002   13.922   21.107)   35.559
   9.488   (  -0.994  -11.709    6.329)   13.347
  10.780   (  13.533   -1.787   12.542)   18.537
  11.905   (  -6.099  -18.709   -0.332)   19.681
  12.325   (  -3.047   -0.201   -1.620)    3.456
  12.434   (   0.781   -4.694   -1.530)    4.999
  12.591   (  -5.802   -4.429   -0.500)    7.316
  12.811   (  -2.352    1.736   -2.591)    3.906
======================= Grid point 293 (51/135) =======================
q-point: ( 0.40  0.20  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.421   (  -9.335   16.169  -14.256)   23.490
   4.812   ( -11.966   20.725   -1.088)   23.956
   5.077   (  -2.055    3.560  -20.869)   21.270
   5.350   ( -10.870   18.827   -2.690)   21.905
   8.460   (   0.461   -0.799   21.123)   21.143
   9.349   (   6.675  -11.562   12.669)   18.405
  10.973   (   4.565   -7.907    7.630)   11.898
  11.684   (   8.195  -14.194   -1.057)   16.424
  12.297   (  -0.381    0.660   -1.746)    1.905
  12.420   (   2.693   -4.664   -0.548)    5.413
  12.525   (  -1.921    3.326   -0.733)    3.910
  12.798   (  -1.670    2.893   -2.571)    4.216
======================= Grid point 307 (52/135) =======================
q-point: ( 0.27  0.27  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.333   (  10.032   17.376  -20.060)   28.372
   4.861   (   7.559   13.093   -0.437)   15.125
   5.321   (  -2.828   -4.899  -17.525)   18.416
   5.539   (   6.382   11.054    0.047)   12.764
   7.947   (  14.637   25.352   21.961)   36.596
   9.431   (  -4.064   -7.040    5.169)    9.633
  10.648   (   1.583    2.742   13.887)   14.244
  11.862   ( -10.093  -17.482    0.046)   20.186
  12.374   (  -2.314   -4.008   -2.591)    5.304
  12.377   (  -1.501   -2.599   -1.175)    3.223
  12.623   (  -2.936   -5.086   -0.692)    5.913
  12.840   (   0.093    0.162   -2.382)    2.390
======================= Grid point 308 (53/135) =======================
q-point: ( 0.33  0.27  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.623   (   2.645   14.976  -22.264)   26.962
   5.020   (   2.316    3.260  -15.227)   15.743
   5.205   (  -8.462    4.081   -4.930)   10.609
   5.671   (  -3.758   12.101    0.132)   12.672
   8.423   (  12.016   17.233   21.775)   30.258
   9.184   (  -5.473  -18.136   15.655)   24.576
  10.838   (   9.556    8.028    5.497)   13.638
  11.488   (  -6.518  -21.481    0.170)   22.449
  12.310   (  -1.446   -0.696   -2.286)    2.793
  12.347   (   0.985   -3.320   -1.591)    3.811
  12.582   (  -1.244    1.266   -0.308)    1.801
  12.838   (  -0.700    0.867   -2.299)    2.555
======================= Grid point 323 (54/135) =======================
q-point: ( 0.33  0.33  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: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.890   (  -5.638   -9.765  -24.543)   27.009
   4.890   (   5.638    9.765  -24.543)   27.009
   5.292   (  -0.000   -0.000    4.530)    4.530
   5.792   (  -0.000   -0.000    0.210)    0.210
   8.797   (  -9.901  -17.149   20.426)   28.449
   8.797   (   9.901   17.149   20.426)   28.449
  11.094   (  -8.539  -14.791    1.775)   17.171
  11.094   (   8.539   14.791    1.775)   17.171
  12.308   (  -0.367   -0.635   -2.122)    2.245
  12.308   (   0.367    0.635   -2.122)    2.245
  12.598   (  -0.000   -0.000   -0.200)    0.200
  12.846   (  -0.000   -0.000   -2.183)    2.183
======================= Grid point 482 (55/135) =======================
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: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.417   (   0.000   -0.000   32.086)   32.086
   1.417   (   0.000    0.000   32.086)   32.086
   2.844   (   0.000    0.000   -4.844)    4.844
   2.844   (   0.000   -0.000   -4.844)    4.844
   3.310   (   0.000    0.000   80.761)   80.761
   8.644   (  -0.000    0.000  -50.988)   50.988
  12.619   (   0.000   -0.000    2.259)    2.259
  12.619   (   0.000    0.000    2.259)    2.259
  12.623   (  -0.000   -0.000   28.571)   28.571
  12.940   (   0.000    0.000   -4.014)    4.014
  12.940   (   0.000   -0.000   -4.014)    4.014
  14.275   (   0.000   -0.000  -10.850)   10.850
======================= Grid point 483 (56/135) =======================
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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.610   (  15.475    8.934   28.427)   33.576
   1.972   (  34.728   20.050   26.861)   48.266
   2.949   (   8.735    5.043   -4.494)   11.043
   3.234   (  29.364   16.953   -0.753)   33.915
   3.439   (  16.461    9.504   69.560)   72.111
   8.561   (  -6.744   -3.894  -49.692)   50.298
  12.428   (  -9.416   -5.436   17.002)   20.181
  12.581   (  -3.049   -1.760    2.363)    4.240
  12.780   (   6.213    3.587   13.406)   15.205
  12.909   (  -2.637   -1.522   -1.561)    3.421
  12.919   (  -1.900   -1.097   -4.217)    4.754
  14.233   (  -3.120   -1.801  -11.508)   12.059
======================= Grid point 484 (57/135) =======================
q-point: ( 0.13  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.053   (  21.407   12.359   22.574)   33.475
   2.570   (  16.524    9.540   37.190)   41.799
   3.216   (  13.639    7.875   -3.831)   16.209
   4.006   (  39.499   22.805   17.339)   48.795
   4.235   (  45.067   26.020   21.586)   56.338
   8.374   (  -8.124   -4.690  -45.451)   46.409
  12.191   ( -10.875   -6.279   17.684)   21.689
  12.492   (  -4.453   -2.571    2.408)    5.678
  12.785   (  -5.731   -3.309    3.647)    7.556
  12.848   (  -4.257   -2.458   -4.506)    6.668
  12.936   (   4.823    2.784    8.733)   10.357
  14.145   (  -4.554   -2.629  -11.799)   12.917
======================= Grid point 485 (58/135) =======================
q-point: ( 0.20  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.536   (  19.548   11.286   18.366)   29.100
   2.877   (  12.349    7.130   31.515)   34.591
   3.539   (  13.655    7.884   -3.379)   16.126
   4.838   (  27.915   16.117   -6.060)   32.799
   5.351   (  50.529   29.173   25.525)   63.685
   8.277   (   1.975    1.140  -36.792)   36.863
  11.943   ( -10.350   -5.975   17.926)   21.544
  12.387   (  -4.503   -2.600    2.147)    5.626
  12.633   (  -6.811   -3.932    2.736)    8.327
  12.731   (  -5.594   -3.230   -4.529)    7.889
  13.012   (   1.164    0.672    7.669)    7.786
  14.012   (  -7.344   -4.240  -10.189)   13.256
======================= Grid point 486 (59/135) =======================
q-point: ( 0.27  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.920   (  13.408    7.741   15.762)   22.094
   3.173   (  13.403    7.738   16.748)   22.804
   3.812   (   9.553    5.516   -3.179)   11.480
   5.203   (   4.358    2.516  -18.290)   18.970
   6.497   (  48.102   27.772   25.094)   60.949
   8.534   (  20.114   11.613  -24.219)   33.556
  11.725   (  -8.517   -4.917   18.787)   21.205
  12.290   (  -3.778   -2.181    1.755)    4.702
  12.490   (  -5.299   -3.059    1.880)    6.401
  12.608   (  -4.755   -2.745   -4.402)    7.037
  12.977   (  -4.148   -2.395    5.374)    7.198
  13.794   ( -11.830   -6.830   -7.747)   15.704
======================= Grid point 487 (60/135) =======================
q-point: ( 0.33  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.147   (   6.523    3.766   14.209)   16.082
   3.487   (  13.690    7.904    2.099)   15.947
   3.963   (   3.505    2.024   -2.911)    4.985
   5.131   (  -8.232   -4.753  -25.114)   26.853
   7.539   (  42.070   24.289   24.921)   54.598
   9.116   (  27.858   16.084  -14.081)   35.115
  11.553   (  -6.355   -3.669   20.070)   21.370
  12.217   (  -2.481   -1.432    1.501)    3.234
  12.396   (  -2.866   -1.655    1.424)    3.603
  12.524   (  -2.495   -1.441   -4.431)    5.285
  12.831   (  -8.188   -4.727    1.969)    9.657
  13.466   ( -16.500   -9.526   -4.751)   19.636
======================= Grid point 488 (61/135) =======================
q-point: ( 0.40  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.238   (   1.920    1.109   13.257)   13.441
   3.789   (  12.170    7.027   -9.543)   16.987
   3.986   (  -0.973   -0.562   -2.322)    2.580
   4.905   ( -10.271   -5.930  -27.947)   30.360
   8.415   (  33.440   19.307   25.525)   46.287
   9.695   (  19.724   11.388   -5.387)   23.404
  11.460   (  -0.294   -0.170   18.147)   18.151
  12.177   (  -1.063   -0.614    1.396)    1.859
  12.348   (  -2.146   -1.239    0.193)    2.485
  12.490   (  -0.638   -0.369   -4.651)    4.709
  12.618   (  -9.987   -5.766   -1.935)   11.693
  13.065   ( -16.850   -9.728   -0.023)   19.457
======================= Grid point 489 (62/135) =======================
q-point: ( 0.47  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.259   (   0.263    0.152   12.701)   12.705
   3.952   (  -1.180   -0.681   -1.689)    2.170
   4.017   (   6.344    3.663  -16.381)   17.945
   4.705   (  -5.769   -3.331  -28.707)   29.470
   9.018   (  16.327    9.427   26.338)   32.390
   9.934   (   1.931    1.115    5.286)    5.737
  11.595   (  10.394    6.001   10.047)   15.652
  12.164   (  -0.222   -0.128    1.350)    1.374
  12.188   ( -12.067   -6.967   -7.591)   15.868
  12.457   (  -2.015   -1.163    2.109)    3.140
  12.485   (  -0.008   -0.005   -4.851)    4.851
  12.783   (  -6.343   -3.662    4.023)    8.356
======================= Grid point 498 (63/135) =======================
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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.978   (  14.541   25.185   23.448)   37.356
   2.427   (  11.395   19.737   34.096)   41.012
   3.161   (   9.101   15.763   -2.433)   18.364
   3.829   (  21.393   37.054   -5.688)   43.162
   3.877   (  21.711   37.604   52.083)   67.809
   8.426   (  -5.032   -8.716  -46.950)   48.017
  12.256   (  -6.174  -10.693   17.565)   21.470
  12.521   (  -2.471   -4.279    2.367)    5.479
  12.831   (  -3.081   -5.337    0.410)    6.176
  12.855   (   0.743    1.287    6.886)    7.044
  12.915   (   0.900    1.559    0.713)    1.936
  14.166   (  -2.766   -4.791  -11.812)   13.044
======================= Grid point 499 (64/135) =======================
q-point: ( 0.13  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.504   (  11.827   30.173   19.058)   37.597
   2.749   (  11.915    9.249   34.155)   37.337
   3.522   (   9.254   21.744   -3.207)   23.848
   4.555   (  24.623   26.350   -4.523)   36.347
   4.855   (  41.629   37.163   30.498)   63.594
   8.276   (  -2.608   -3.746  -40.657)   40.913
  12.030   (  -7.998   -9.162   18.006)   21.729
  12.421   (  -3.560   -4.454    2.280)    6.141
  12.702   (  -5.776   -4.407    1.264)    7.374
  12.815   (  -4.136   -1.159   -3.117)    5.307
  12.971   (   2.485    0.162    8.137)    8.509
  14.042   (  -4.403   -7.885  -10.703)   14.004
======================= Grid point 500 (65/135) =======================
q-point: ( 0.20  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.917   (   7.413   19.653   17.410)   27.282
   3.091   (  10.354   19.637   20.417)   30.161
   3.869   (   3.611   23.978   -3.195)   24.458
   5.075   (  14.190    6.621  -14.232)   21.160
   5.964   (  45.625   31.669   26.709)   61.627
   8.350   (  11.247    6.725  -29.792)   32.546
  11.810   (  -7.929   -7.056   18.953)   21.722
  12.323   (  -2.961   -3.415    2.070)    4.971
  12.556   (  -6.309   -3.292    0.818)    7.163
  12.711   (  -6.934    0.264   -4.267)    8.146
  12.972   (  -0.093   -4.693    6.580)    8.082
  13.850   (  -6.597  -12.045   -8.583)   16.194
======================= Grid point 501 (66/135) =======================
q-point: ( 0.27  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.202   (   5.792   15.072    9.487)   18.728
   3.431   (   4.624   21.774   11.632)   25.115
   4.102   (  -3.611   22.361   -3.099)   22.862
   5.178   (  -1.644   -4.133  -22.412)   22.849
   7.028   (  44.259   25.183   25.410)   56.910
   8.788   (  25.609   13.230  -18.327)   34.157
  11.617   (  -6.873   -6.067   20.559)   22.511
  12.256   (  -2.334   -1.089    1.888)    3.194
  12.430   (  -4.506   -2.698    0.769)    5.308
  12.598   (  -6.160    0.572   -4.803)    7.832
  12.875   (  -3.944   -7.289    3.515)    9.002
  13.568   ( -10.237  -15.573   -6.134)   19.620
======================= Grid point 502 (67/135) =======================
q-point: ( 0.33  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.398   (  -2.686   17.884    7.288)   19.498
   3.704   (   6.171   15.645   -0.320)   16.821
   4.191   (  -9.878   19.606   -2.997)   22.157
   5.019   (  -9.148   -5.688  -26.462)   28.570
   7.974   (  39.090   19.322   25.458)   50.492
   9.393   (  27.071   11.175   -9.735)   30.862
  11.459   (  -3.588   -5.998   21.027)   22.158
  12.217   (  -2.335    1.072    1.448)    2.949
  12.353   (  -1.529   -2.641    1.634)    3.462
  12.512   (  -3.571   -1.520   -5.562)    6.782
  12.711   (  -6.783   -5.428   -0.609)    8.709
  13.195   ( -13.362  -16.935   -2.778)   21.750
======================= Grid point 503 (68/135) =======================
q-point: ( 0.40  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.475   (  -9.573   20.120    8.767)   23.944
   3.954   (   6.808   10.511  -11.627)   17.088
   4.178   ( -12.578   18.322   -2.884)   22.410
   4.799   (  -9.380   -4.101  -27.727)   29.557
   8.731   (  28.608   12.341   26.059)   40.618
   9.822   (  12.442    1.197    0.565)   12.512
  11.438   (   8.125   -2.329   15.200)   17.392
  12.186   (  -2.488   -1.270   -0.535)    2.844
  12.270   (  -5.660   -6.145   -0.669)    8.381
  12.438   (  -0.617   -3.479   -3.052)    4.669
  12.576   (  -5.983    2.252   -3.822)    7.448
  12.832   (  -8.838  -12.792    2.075)   15.686
======================= Grid point 504 (69/135) =======================
q-point: ( 0.47  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.489   ( -11.783   20.408    9.202)   25.299
   4.089   (  -2.507    4.343  -15.683)   16.465
   4.152   ( -10.963   18.989   -3.024)   22.134
   4.678   (  -0.718    1.244  -27.711)   27.748
   9.084   (   1.763   -3.054   26.587)   26.820
   9.888   (   2.856   -4.947    7.686)    9.576
  11.600   (   4.119   -7.135    9.957)   12.923
  12.027   (   3.664   -6.346   -8.117)   10.935
  12.206   (  -2.103    3.642    1.753)    4.556
  12.431   (   0.660   -1.143    1.176)    1.768
  12.538   (  -2.582    4.471   -4.878)    7.103
  12.692   (   2.962   -5.130    3.369)    6.814
======================= Grid point 515 (70/135) =======================
q-point: ( 0.13  0.13  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.938   (   6.528   11.307   23.422)   26.815
   3.152   (  17.777   30.791   14.902)   38.551
   4.001   (  13.914   24.100   -2.772)   27.966
   4.956   (   7.968   13.800  -14.249)   21.377
   5.690   (  25.916   44.888   29.747)   59.761
   8.287   (   3.590    6.218  -32.174)   32.966
  11.845   (  -5.288   -9.160   19.408)   22.103
  12.330   (  -2.457   -4.256    2.122)    5.353
  12.610   (  -2.929   -5.072   -0.187)    5.860
  12.783   (  -1.446   -2.504   -4.605)    5.438
  12.937   (  -1.987   -3.441    7.475)    8.465
  13.835   (  -7.330  -12.695   -8.686)   17.040
======================= Grid point 516 (71/135) =======================
q-point: ( 0.20  0.13  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.178   (   8.363   11.098    9.401)   16.778
   3.685   (   6.291   32.418   12.435)   35.287
   4.412   (   2.174   28.216   -2.650)   28.424
   5.112   (   3.877   -1.884  -20.625)   21.071
   6.603   (  34.441   32.532   27.344)   54.701
   8.553   (  15.531   13.594  -21.837)   30.048
  11.663   (  -6.180   -7.905   21.962)   24.146
  12.264   (  -1.214   -2.494    2.193)    3.536
  12.498   (  -5.160   -2.809   -0.955)    5.952
  12.701   (  -5.153   -2.668   -5.700)    8.134
  12.850   (  -2.135   -5.942    4.801)    7.932
  13.546   (  -8.264  -17.846   -6.601)   20.744
======================= Grid point 517 (72/135) =======================
q-point: ( 0.27  0.13  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.441   (   9.179   10.708   -3.087)   14.437
   4.010   (  -4.930   31.460    9.770)   33.309
   4.619   (  -9.639   26.678   -4.437)   28.711
   5.085   (  -2.992   -3.634  -23.004)   23.481
   7.504   (  37.243   22.953   26.428)   51.110
   9.047   (  24.407   12.774  -12.970)   30.449
  11.470   (  -6.198   -9.101   24.024)   26.427
  12.233   (  -0.054   -2.026    2.487)    3.208
  12.388   (  -4.866   -1.470   -0.669)    5.127
  12.563   (  -5.348   -5.790   -5.514)    9.619
  12.757   (  -3.670   -2.167   -0.099)    4.264
  13.200   (  -9.547  -20.358   -4.397)   22.911
======================= Grid point 518 (73/135) =======================
q-point: ( 0.33  0.13  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.689   (   6.069   11.293  -10.759)   16.737
   4.178   (  -9.036   28.690    5.294)   30.542
   4.634   ( -14.221   19.212   -9.778)   25.826
   4.964   ( -10.621    4.459  -19.050)   22.262
   8.317   (  33.694   15.493   26.374)   45.507
   9.555   (  20.507    5.145   -4.214)   21.558
  11.305   (   0.056   -9.472   21.723)   23.698
  12.184   (  -2.260   -8.518    0.777)    8.846
  12.302   (  -3.421   -1.199   -0.061)    3.625
  12.437   (  -2.497   -3.131   -5.167)    6.538
  12.685   (  -4.282    1.440   -2.760)    5.294
  12.838   (  -9.255  -16.971   -1.528)   19.391
======================= Grid point 519 (74/135) =======================
q-point: ( 0.40  0.13  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.844   (  -2.708   12.838  -10.905)   17.061
   4.304   (  -7.961   25.722   -0.568)   26.931
   4.540   ( -11.475    8.844  -18.429)   23.442
   4.848   ( -13.476   16.254  -11.443)   24.015
   8.908   (  18.849    5.494   26.590)   33.053
   9.760   (   5.008   -6.438    7.004)   10.751
  11.353   (  13.139   -6.550   13.252)   19.777
  11.947   (  -6.135  -16.588   -4.652)   18.288
  12.264   (  -1.067    0.581    0.717)    1.411
  12.427   (   0.455    2.314   -1.392)    2.739
  12.584   (  -2.266   -9.782    1.006)   10.091
  12.654   (  -3.524    3.366   -4.294)    6.495
======================= Grid point 531 (75/135) =======================
q-point: ( 0.20  0.20  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.408   (   7.029   12.174   -2.747)   14.323
   4.241   (  12.503   21.656    7.704)   26.166
   4.912   (  11.517   19.949   -1.638)   23.093
   5.045   (  -1.972   -3.416  -23.860)   24.184
   7.269   (  19.863   34.404   27.551)   48.345
   8.886   (  11.260   19.502  -14.492)   26.780
  11.481   (  -6.171  -10.689   25.640)   28.456
  12.217   (  -1.509   -2.614    2.198)    3.734
  12.437   (  -1.886   -3.266   -2.377)    4.458
  12.572   (  -4.956   -8.584   -3.670)   10.569
  12.788   (  -0.749   -1.298   -1.080)    1.847
  13.168   ( -11.164  -19.337   -4.753)   22.828
======================= Grid point 532 (76/135) =======================
q-point: ( 0.27  0.20  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.660   (   8.168   11.468  -13.105)   19.234
   4.523   (  -1.000   15.771   -3.599)   16.207
   4.960   (  -7.679    6.912  -15.622)   18.729
   5.217   (  -0.269   18.468   -2.933)   18.701
   7.965   (  26.475   23.823   27.403)   44.938
   9.300   (  16.432   12.464   -7.067)   21.802
  11.247   (  -5.980  -12.853   25.993)   29.607
  12.130   (  -3.082   -9.150    1.054)    9.713
  12.364   (  -4.794   -1.598   -2.681)    5.720
  12.428   (  -1.343   -5.785   -3.799)    7.050
  12.757   (  -3.558   -0.375   -2.795)    4.540
  12.817   (  -6.895  -15.095   -3.029)   16.869
======================= Grid point 533 (77/135) =======================
q-point: ( 0.33  0.20  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.900   (   5.991   10.488  -20.086)   23.437
   4.563   (  -3.395    2.771  -22.328)   22.754
   4.907   ( -13.570   19.403   -0.339)   23.680
   5.311   ( -10.975   22.073   -1.871)   24.722
   8.614   (  24.324   14.819   27.047)   39.279
   9.601   (   8.387   -0.644    3.937)    9.288
  11.104   (   5.066   -9.460   18.692)   21.553
  11.887   (  -7.812  -18.995   -1.635)   20.603
  12.277   (  -2.083   -1.641   -2.856)    3.897
  12.400   (   2.010   -3.610   -1.783)    4.501
  12.581   (  -5.226   -4.733   -0.408)    7.062
  12.742   (  -2.625    2.160   -4.096)    5.323
======================= Grid point 534 (78/135) =======================
q-point: ( 0.40  0.20  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.020   (  -3.679    6.372  -22.261)   23.445
   4.530   (   1.255   -2.174  -30.980)   31.082
   4.885   ( -15.149   26.238    6.838)   31.059
   5.319   ( -13.691   23.713   -1.453)   27.420
   8.956   (   0.200   -0.347   26.733)   26.736
   9.612   (   4.136   -7.163   12.121)   14.674
  11.176   (   6.096  -10.558   11.879)   17.022
  11.648   (   6.935  -12.011   -2.693)   14.128
  12.254   (   0.519   -0.898   -2.256)    2.483
  12.403   (   3.793   -6.570   -1.280)    7.693
  12.514   (  -2.343    4.058   -0.167)    4.688
  12.730   (  -1.931    3.344   -4.136)    5.658
======================= Grid point 548 (79/135) =======================
q-point: ( 0.27  0.27  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: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.897   (   6.951   12.040  -21.143)   25.304
   4.599   (  -2.225   -3.854  -26.929)   27.294
   5.151   (   5.547    9.608    3.008)   11.495
   5.532   (   6.497   11.254   -0.914)   13.027
   8.460   (  14.658   25.389   27.554)   40.233
   9.513   (   4.096    7.094    2.163)    8.472
  11.012   (  -4.813   -8.336   21.313)   23.386
  11.856   ( -10.318  -17.871   -0.748)   20.649
  12.306   (  -2.239   -3.878   -3.994)    6.001
  12.340   (  -1.737   -3.009   -2.656)    4.373
  12.610   (  -2.312   -4.005   -0.491)    4.651
  12.777   (   0.120    0.208   -3.831)    3.839
======================= Grid point 549 (80/135) =======================
q-point: ( 0.33  0.27  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: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.113   (   3.143   10.753  -26.308)   28.594
   4.486   (  -2.327   -6.986  -31.393)   32.245
   5.274   (  -4.590   13.671    6.872)   15.975
   5.667   (  -3.818   12.440   -0.723)   13.033
   8.924   (  11.419   16.467   26.528)   33.246
   9.511   (  -2.084   -8.866   14.954)   17.510
  10.998   (   6.594    0.470   10.428)   12.347
  11.488   (  -6.811  -19.847   -0.327)   20.985
  12.249   (  -1.156   -0.884   -3.625)    3.907
  12.303   (   1.363   -4.181   -2.914)    5.276
  12.579   (  -1.245    1.523    0.089)    1.969
  12.777   (  -0.765    1.146   -3.746)    3.992
======================= Grid point 564 (81/135) =======================
q-point: ( 0.33  0.33  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: 192
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.319   (  -4.799   -8.312  -30.306)   31.790
   4.319   (   4.799    8.312  -30.306)   31.790
   5.412   (  -0.000   -0.000    7.214)    7.214
   5.791   (  -0.000   -0.000   -0.536)    0.536
   9.255   (  -7.902  -13.687   23.277)   28.135
   9.255   (   7.902   13.687   23.277)   28.135
  11.143   (  -6.714  -11.629    3.110)   13.784
  11.143   (   6.714   11.629    3.110)   13.784
  12.250   (  -0.597   -1.034   -3.584)    3.778
  12.250   (   0.597    1.034   -3.584)    3.778
  12.597   (  -0.000   -0.000    0.168)    0.168
  12.788   (  -0.000   -0.000   -3.568)    3.568
======================= Grid point 723 (82/135) =======================
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: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.989   (  -0.000    0.000   25.556)   25.556
   1.989   (  -0.000    0.000   25.556)   25.556
   2.696   (   0.000   -0.000  -10.513)   10.513
   2.696   (   0.000    0.000  -10.513)   10.513
   4.852   (   0.000   -0.000   75.270)   75.270
   7.541   (   0.000   -0.000  -60.417)   60.417
  12.676   (   0.000   -0.000    3.559)    3.559
  12.676   (   0.000    0.000    3.559)    3.559
  12.851   (   0.000    0.000   -4.805)    4.805
  12.851   (   0.000    0.000   -4.805)    4.805
  13.165   (  -0.000    0.000   25.924)   25.924
  14.007   (  -0.000    0.000  -16.274)   16.274
======================= Grid point 724 (83/135) =======================
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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.134   (  11.949    6.899   24.037)   27.715
   2.457   (  33.891   19.567   22.286)   45.034
   2.809   (   9.357    5.402  -10.056)   14.760
   3.099   (  31.398   18.128   -9.979)   37.604
   4.882   (   4.163    2.404   72.179)   72.339
   7.487   (  -4.195   -2.422  -58.828)   59.027
  12.626   (  -4.697   -2.712    6.295)    8.310
  12.642   (  -2.867   -1.655    3.745)    4.998
  12.826   (  -2.226   -1.285   -5.060)    5.675
  12.842   (  -1.303   -0.752   -4.611)    4.850
  13.177   (   1.225    0.708   22.944)   22.987
  13.964   (  -3.204   -1.850  -15.846)   16.272
======================= Grid point 725 (84/135) =======================
q-point: ( 0.13  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.490   (  17.818   10.287   20.878)   29.313
   3.092   (  14.364    8.293   -9.017)   18.879
   3.203   (  25.114   14.499   26.833)   39.509
   3.927   (  36.122   20.855  -10.841)   43.096
   5.172   (  23.963   13.835   56.356)   62.783
   7.399   (  -1.673   -0.966  -52.870)   52.905
  12.461   (  -9.137   -5.275   10.504)   14.888
  12.554   (  -4.521   -2.610    3.901)    6.517
  12.749   (  -4.406   -2.544   -5.389)    7.411
  12.767   (  -5.389   -3.111   -4.472)    7.663
  13.220   (   2.050    1.184   17.859)   18.015
  13.885   (  -3.482   -2.010  -14.515)   15.061
======================= Grid point 726 (85/135) =======================
q-point: ( 0.20  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.902   (  16.998    9.814   17.917)   26.575
   3.429   (  14.129    8.158   -8.001)   18.172
   3.542   (   6.618    3.821   33.742)   34.596
   4.590   (  18.686   10.788  -17.674)   27.891
   5.980   (  43.131   24.902   35.944)   61.420
   7.516   (  13.918    8.036  -39.759)   42.884
  12.237   (  -9.683   -5.591   12.390)   16.689
  12.444   (  -4.858   -2.805    3.644)    6.689
  12.611   (  -7.499   -4.329   -5.155)   10.077
  12.632   (  -5.454   -3.149   -5.335)    8.254
  13.241   (  -0.876   -0.506   14.546)   14.581
  13.791   (  -5.189   -2.996  -12.098)   13.500
======================= Grid point 727 (86/135) =======================
q-point: ( 0.27  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.239   (  11.798    6.811   15.804)   20.865
   3.619   (   1.679    0.970   27.355)   27.424
   3.711   (   9.859    5.692   -7.303)   13.525
   4.759   (  -2.258   -1.304  -26.155)   26.285
   7.024   (  45.768   26.424   27.389)   59.523
   8.053   (  30.945   17.866  -24.629)   43.398
  12.034   (  -7.740   -4.469   13.138)   15.890
  12.339   (  -4.102   -2.368    3.199)    5.716
  12.449   (  -6.143   -3.547   -6.338)    9.513
  12.512   (  -4.635   -2.676   -5.081)    7.380
  13.161   (  -6.152   -3.552   12.509)   14.386
  13.626   (  -9.430   -5.444   -9.387)   14.376
======================= Grid point 728 (87/135) =======================
q-point: ( 0.33  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.437   (   5.510    3.181   14.429)   15.770
   3.656   (   1.881    1.086   15.101)   15.256
   3.871   (   3.963    2.288   -6.718)    8.129
   4.594   ( -10.135   -5.852  -28.243)   30.572
   8.037   (  41.594   24.014   24.585)   53.955
   8.831   (  34.575   19.962  -15.497)   42.825
  11.878   (  -6.089   -3.515   13.586)   15.297
  12.261   (  -2.605   -1.504    2.950)    4.213
  12.332   (  -4.257   -2.458   -8.061)    9.441
  12.428   (  -2.610   -1.507   -5.034)    5.867
  12.966   ( -10.293   -5.943   11.067)   16.241
  13.359   ( -13.400   -7.737   -6.493)   16.781
======================= Grid point 729 (88/135) =======================
q-point: ( 0.40  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.508   (   1.220    0.704   13.443)   13.517
   3.713   (   3.014    1.740    2.828)    4.485
   3.907   (  -0.288   -0.166   -5.942)    5.951
   4.353   (  -9.949   -5.744  -26.491)   28.875
   8.910   (  33.548   19.369   23.582)   45.351
   9.563   (  27.240   15.727   -8.989)   32.713
  11.749   (  -4.736   -2.734   11.993)   13.181
  12.219   (  -1.091   -0.630    2.916)    3.176
  12.222   (  -6.595   -3.807  -11.304)   13.629
  12.390   (  -0.857   -0.495   -5.231)    5.324
  12.727   (  -9.145   -5.280   10.837)   15.131
  13.045   ( -12.653   -7.305   -2.791)   14.874
======================= Grid point 730 (89/135) =======================
q-point: ( 0.47  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.517   (  -0.042   -0.024   12.780)   12.780
   3.783   (   2.420    1.397   -5.762)    6.404
   3.887   (  -0.796   -0.460   -5.207)    5.287
   4.170   (  -4.984   -2.877  -23.734)   24.421
   9.517   (  16.500    9.526   22.792)   29.706
   9.991   (   8.543    4.932   -0.587)    9.882
  11.717   (   3.229    1.864    3.620)    5.197
  11.995   ( -11.492   -6.635  -10.688)   17.038
  12.205   (  -0.229   -0.132    2.934)    2.946
  12.381   (  -0.111   -0.064   -5.427)    5.428
  12.597   (  -2.214   -1.278   10.091)   10.410
  12.832   (  -5.010   -2.893    0.259)    5.791
======================= Grid point 739 (90/135) =======================
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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.429   (  12.026   20.830   21.372)   32.176
   2.962   (  13.883   24.046   13.609)   30.922
   3.079   (  12.256   21.228    1.550)   24.561
   3.672   (  20.419   35.367  -10.542)   42.177
   5.040   (   9.404   16.288   62.555)   65.322
   7.415   (  -2.085   -3.612  -55.091)   55.248
  12.509   (  -4.878   -8.449    9.213)   13.419
  12.583   (  -2.474   -4.286    3.839)    6.263
  12.791   (  -2.476   -4.288   -4.152)    6.462
  12.793   (  -1.426   -2.469   -5.280)    6.000
  13.202   (   0.925    1.601   19.141)   19.230
  13.902   (  -2.330   -4.036  -14.963)   15.672
======================= Grid point 740 (91/135) =======================
q-point: ( 0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.880   (  10.614   26.927   18.285)   34.235
   3.387   (  12.371   16.485    8.558)   22.317
   3.451   (   8.042   12.923   15.560)   21.767
   4.346   (  22.199   18.511  -14.740)   32.446
   5.609   (  29.206   29.705   43.028)   59.890
   7.420   (   4.989    4.515  -45.540)   46.035
  12.316   (  -7.096   -8.234   11.452)   15.790
  12.480   (  -3.729   -4.838    3.427)    7.004
  12.671   (  -5.680   -4.607   -3.924)    8.300
  12.722   (  -4.533   -1.914   -5.593)    7.449
  13.222   (   0.636   -1.463   15.557)   15.638
  13.807   (  -3.046   -6.005  -12.959)   14.604
======================= Grid point 741 (92/135) =======================
q-point: ( 0.20  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.273   (   3.910   23.926   17.102)   29.669
   3.588   (   4.940    3.025   28.501)   29.083
   3.770   (   3.762   24.638   -6.640)   25.793
   4.704   (   7.695    0.535  -23.010)   24.269
   6.538   (  40.778   29.943   30.484)   59.065
   7.751   (  22.028   15.491  -30.649)   40.799
  12.120   (  -7.372   -5.520   12.851)   15.810
  12.375   (  -3.134   -3.682    2.773)    5.574
  12.519   (  -6.582   -3.818   -4.235)    8.709
  12.606   (  -6.777   -0.891   -6.072)    9.143
  13.175   (  -2.055   -6.306   13.036)   14.626
  13.663   (  -4.839  -10.037  -10.363)   15.217
======================= Grid point 742 (93/135) =======================
q-point: ( 0.27  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.466   (   1.232    8.091   16.596)   18.505
   3.719   (  -1.040   15.649   16.462)   22.737
   4.008   (  -3.537   22.223   -6.263)   23.358
   4.680   (  -5.249   -5.132  -27.049)   28.027
   7.542   (  42.561   25.029   25.720)   55.672
   8.424   (  33.359   18.644  -19.032)   42.692
  11.953   (  -6.856   -3.774   13.907)   15.958
  12.304   (  -2.638   -0.787    2.839)    3.955
  12.382   (  -4.986   -3.056   -5.961)    8.350
  12.488   (  -5.545   -0.638   -5.848)    8.084
  13.022   (  -5.887  -10.026   10.943)   15.966
  13.433   (  -8.051  -13.764   -7.787)   17.746
======================= Grid point 743 (94/135) =======================
q-point: ( 0.33  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.539   (   1.800    2.918    9.371)    9.979
   3.825   (  -6.507   18.455    9.554)   21.777
   4.092   (  -9.056   17.132   -7.188)   20.668
   4.495   ( -10.877   -1.517  -24.704)   27.035
   8.473   (  38.700   19.518   23.996)   49.543
   9.186   (  33.163   15.349  -12.018)   38.468
  11.799   (  -6.028   -4.786   13.955)   15.937
  12.240   (  -4.218   -2.350   -2.105)    5.267
  12.299   (  -3.001    0.387   -4.445)    5.377
  12.416   (  -2.041   -0.278   -4.048)    4.542
  12.794   (  -7.771  -10.289    9.020)   15.736
  13.127   (  -9.782  -15.047   -4.750)   18.565
======================= Grid point 744 (95/135) =======================
q-point: ( 0.40  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.593   (   0.237    3.995    3.740)    5.477
   3.880   (  -6.704   16.861    3.280)   18.439
   4.058   ( -10.014   10.306  -10.073)   17.549
   4.313   ( -12.118    6.648  -18.178)   22.836
   9.228   (  28.608   12.353   23.073)   38.773
   9.790   (  21.048    6.816   -4.754)   22.630
  11.675   (  -0.891   -4.733    9.483)   10.636
  12.076   (  -8.360  -10.007  -10.301)   16.617
  12.266   (  -3.408    3.828    1.981)    5.495
  12.412   (   0.041    2.217   -3.416)    4.072
  12.597   (  -4.040   -6.673    7.677)   10.945
  12.846   (  -5.437  -12.048   -1.281)   13.280
======================= Grid point 745 (96/135) =======================
q-point: ( 0.47  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.616   (  -2.768    4.795    2.672)    6.148
   3.930   (  -9.804   16.981    1.536)   19.668
   3.994   (  -2.948    5.106  -14.986)   16.104
   4.238   (  -8.419   14.582  -13.375)   21.504
   9.582   (   1.891   -3.275   22.328)   22.646
   9.992   (   2.498   -4.327    1.779)    5.303
  11.706   (   2.640   -4.573    1.529)    5.497
  11.857   (   2.616   -4.530   -8.209)    9.734
  12.256   (  -2.594    4.493    3.168)    6.079
  12.433   (  -2.660    4.607   -5.476)    7.635
  12.530   (   2.887   -5.001    8.086)    9.936
  12.735   (   3.742   -6.481    0.230)    7.487
======================= Grid point 756 (97/135) =======================
q-point: ( 0.13  0.13  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.456   (  16.038   27.778   15.167)   35.480
   3.491   (   0.876    1.518   30.597)   30.647
   3.911   (  14.389   24.923   -6.420)   29.486
   4.600   (   3.826    6.626  -21.924)   23.221
   6.314   (  22.929   39.714   32.870)   56.421
   7.634   (  10.249   17.751  -33.578)   39.340
  12.159   (  -4.187   -7.253   12.523)   15.066
  12.380   (  -2.598   -4.499    2.457)    5.747
  12.571   (  -3.292   -5.703   -3.285)    7.359
  12.669   (  -2.103   -3.642   -6.252)    7.535
  13.151   (  -3.271   -5.666   13.147)   14.685
  13.643   (  -5.992  -10.378  -10.633)   16.021
======================= Grid point 757 (98/135) =======================
q-point: ( 0.20  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.493   (  -0.088    0.470   21.562)   21.568
   3.935   (   5.592   26.640   11.847)   29.686
   4.325   (   0.391   28.400   -6.479)   29.132
   4.651   (   3.058   -3.675  -24.676)   25.135
   7.152   (  32.044   31.610   27.442)   52.717
   8.121   (  22.536   21.396  -22.140)   38.155
  12.024   (  -5.253   -4.417   14.462)   16.008
  12.316   (  -0.967   -2.094    2.394)    3.324
  12.445   (  -5.286   -4.135   -3.936)    7.781
  12.575   (  -4.979   -2.703   -6.745)    8.809
  13.001   (  -4.581  -10.263   10.185)   15.167
  13.400   (  -6.866  -15.828   -8.285)   19.139
======================= Grid point 758 (99/135) =======================
q-point: ( 0.27  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.512   (   2.019    0.401   10.636)   10.834
   4.147   (  -5.401   15.400   -1.959)   16.437
   4.460   (  -9.530   19.619   -7.718)   23.137
   4.710   (  -4.691   15.862  -12.295)   20.611
   8.019   (  35.950   23.154   24.800)   49.432
   8.785   (  29.999   17.710  -14.337)   37.671
  11.869   (  -7.481   -5.274   16.188)   18.597
  12.293   (  -2.021   -2.371    2.673)    4.105
  12.323   (  -2.829   -2.088   -6.013)    6.966
  12.456   (  -4.526   -3.431   -4.612)    7.316
  12.808   (  -5.737   -9.354    5.778)   12.402
  13.100   (  -7.524  -18.572   -5.929)   20.897
======================= Grid point 759 (100/135) =======================
q-point: ( 0.33  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.570   (   3.423    1.337    0.508)    3.710
   4.112   (  -5.718    0.837  -18.535)   19.415
   4.468   ( -15.775   26.232    3.414)   30.800
   4.744   ( -14.907   25.001   -6.812)   29.894
   8.820   (  33.107   15.615   23.404)   43.447
   9.442   (  27.398   10.477   -8.063)   30.421
  11.670   (  -6.292   -8.591   14.969)   18.370
  12.120   (  -8.643  -10.299   -8.055)   15.673
  12.325   (   1.073   -1.357    2.678)    3.188
  12.387   (   0.167   -1.063    0.039)    1.076
  12.645   (  -5.753   -2.486    0.080)    6.268
  12.798   (  -6.176  -16.332   -3.076)   17.730
======================= Grid point 760 (101/135) =======================
q-point: ( 0.40  0.13  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.639   (   2.035    1.463   -6.089)    6.585
   4.025   (  -2.026   -2.630  -23.291)   23.527
   4.450   ( -18.211   30.015    7.541)   35.908
   4.722   ( -17.431   28.125   -5.083)   33.477
   9.402   (  18.673    5.116   22.059)   29.351
   9.840   (  11.276   -1.216    0.190)   11.343
  11.561   (   4.749   -6.815    7.562)   11.233
  11.828   (  -8.046  -13.965   -7.053)   17.593
  12.318   (   0.644   -1.831    4.625)    5.015
  12.428   (   2.965    1.542    1.253)    3.569
  12.549   (  -5.175   -0.161   -3.499)    6.249
  12.606   (  -0.144   -8.623   -1.421)    8.740
======================= Grid point 772 (102/135) =======================
q-point: ( 0.20  0.20  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.495   (   0.397    0.687   11.761)   11.788
   4.263   (   0.884    1.532  -12.830)   12.951
   4.676   (   7.074   12.253   -5.540)   15.194
   4.855   (  11.886   20.587   -4.244)   24.148
   7.801   (  19.411   33.621   25.482)   46.438
   8.596   (  15.051   26.069  -15.471)   33.845
  11.918   (  -4.154   -7.195   17.938)   19.769
  12.274   (  -1.944   -3.367    3.457)    5.202
  12.356   (  -2.137   -3.701   -5.706)    7.129
  12.481   (  -3.940   -6.824   -5.370)    9.535
  12.813   (  -4.357   -7.547    3.929)    9.559
  13.062   (  -9.917  -17.177   -6.060)   20.740
======================= Grid point 773 (103/135) =======================
q-point: ( 0.27  0.20  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.533   (   2.424    2.013    1.400)    3.448
   4.169   (  -4.778   -5.834  -25.314)   26.413
   4.917   (  -2.934   22.878    6.155)   23.872
   5.157   (  -1.717   21.921   -3.938)   22.338
   8.482   (  25.629   23.678   23.939)   42.315
   9.140   (  21.176   18.104   -9.882)   29.560
  11.708   (  -8.829  -11.724   19.646)   24.523
  12.136   (  -5.835  -12.531   -1.897)   13.953
  12.296   (  -3.472   -2.193   -2.541)    4.829
  12.365   (   0.606   -4.820   -2.750)    5.582
  12.707   (  -3.307   -1.679   -1.752)    4.102
  12.757   (  -6.420  -13.409   -3.020)   15.171
======================= Grid point 774 (104/135) =======================
q-point: ( 0.33  0.20  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.611   (   3.228    2.980   -7.340)    8.555
   4.033   (  -3.741   -5.730  -25.709)   26.604
   5.020   ( -12.226   25.613    7.689)   29.405
   5.260   ( -11.502   23.866   -3.614)   26.738
   9.113   (  23.425   14.220   22.317)   35.341
   9.615   (  16.320    7.091   -3.027)   18.049
  11.451   (  -3.554  -12.879   15.276)   20.294
  11.828   (  -9.968  -17.753   -4.772)   20.911
  12.243   (   0.960   -4.034    0.340)    4.161
  12.360   (   3.473   -5.590   -2.579)    7.068
  12.578   (  -4.595   -2.895    0.207)    5.435
  12.661   (  -3.053    2.107   -3.885)    5.371
======================= Grid point 775 (105/135) =======================
q-point: ( 0.40  0.20  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.666   (  -0.799    1.384  -11.391)   11.503
   3.960   (   1.769   -3.063  -24.662)   24.914
   5.038   ( -15.230   26.379    7.926)   31.474
   5.274   ( -14.329   24.819   -3.432)   28.863
   9.438   (   0.616   -1.067   20.708)   20.744
   9.776   (   2.379   -4.120    3.695)    6.024
  11.390   (   5.910  -10.237    8.703)   14.678
  11.583   (   5.261   -9.111   -3.583)   11.114
  12.234   (   2.387   -4.134    0.852)    4.849
  12.362   (   4.849   -8.398   -2.992)   10.149
  12.526   (  -2.464    4.268    1.427)    5.131
  12.645   (  -2.167    3.754   -4.220)    6.050
======================= Grid point 789 (106/135) =======================
q-point: ( 0.27  0.27  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: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.597   (   2.502    4.333   -7.449)    8.974
   4.037   (  -3.922   -6.793  -26.509)   27.645
   5.272   (   7.090   12.280    7.330)   15.962
   5.494   (   6.704   11.612   -3.104)   13.763
   8.965   (  14.027   24.295   22.570)   36.005
   9.495   (   9.462   16.390   -4.425)   19.435
  11.420   (  -8.688  -15.049   18.586)   25.444
  11.822   ( -10.620  -18.395   -3.105)   21.467
  12.233   (  -2.100   -3.637   -2.802)    5.048
  12.275   (  -1.806   -3.129   -3.712)    5.180
  12.610   (  -1.434   -2.484    0.599)    2.931
  12.699   (  -0.002   -0.004   -3.757)    3.757
======================= Grid point 790 (107/135) =======================
q-point: ( 0.33  0.27  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: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.690   (   1.905    4.852  -14.422)   15.335
   3.915   (  -1.967   -5.937  -24.281)   25.073
   5.421   (  -4.170   13.813    7.432)   16.230
   5.634   (  -3.917   12.881   -2.819)   13.756
   9.395   (  10.324   13.951   19.906)   26.410
   9.711   (   3.099    1.532    4.538)    5.705
  11.219   (   1.261   -8.029   10.936)   13.626
  11.467   (  -6.567  -17.102   -1.997)   18.428
  12.185   (  -0.279   -1.707   -2.235)    2.826
  12.237   (   1.685   -4.791   -3.568)    6.207
  12.591   (  -1.149    1.669    1.292)    2.403
  12.699   (  -0.853    1.450   -3.870)    4.220
======================= Grid point 805 (108/135) =======================
q-point: ( 0.33  0.33  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: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.798   (  -2.878   -4.984  -20.286)   21.087
   3.798   (   2.878    4.984  -20.286)   21.087
   5.559   (  -0.000   -0.000    7.273)    7.273
   5.762   (  -0.000   -0.000   -2.568)    2.568
   9.637   (  -4.485   -7.768   14.126)   16.733
   9.637   (   4.485    7.768   14.126)   16.733
  11.207   (  -3.640   -6.304    2.963)    7.859
  11.207   (   3.640    6.304    2.963)    7.859
  12.179   (  -0.579   -1.003   -3.209)    3.411
  12.179   (   0.579    1.003   -3.209)    3.411
  12.610   (  -0.000   -0.000    1.277)    1.277
  12.714   (  -0.000   -0.000   -3.702)    3.702
======================= Grid point 964 (109/135) =======================
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: 81
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.418   (  -0.000    0.000  -17.886)   17.886
   2.418   (  -0.000    0.000  -17.886)   17.886
   2.418   (  -0.000    0.000   17.886)   17.886
   2.418   (   0.000    0.000   17.886)   17.886
   6.270   (   0.000    0.000  -68.363)   68.363
   6.270   (   0.000    0.000   68.363)   68.363
  12.757   (   0.000    0.000   -4.564)    4.564
  12.757   (   0.000    0.000   -4.564)    4.564
  12.757   (   0.000   -0.000    4.564)    4.564
  12.757   (   0.000   -0.000    4.564)    4.564
  13.634   (   0.000    0.000  -21.452)   21.452
  13.634   (   0.000   -0.000   21.452)   21.452
======================= Grid point 966 (110/135) =======================
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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.542   (  10.274    5.932  -17.151)   20.855
   2.542   (  10.274    5.932   17.151)   20.855
   2.840   (  32.202   18.592  -16.344)   40.617
   2.840   (  32.202   18.592   16.344)   40.617
   6.251   (  -0.840   -0.485  -66.311)   66.319
   6.251   (  -0.840   -0.485   66.311)   66.319
  12.727   (  -2.576   -1.487   -4.812)    5.657
  12.727   (  -2.576   -1.487    4.812)    5.657
  12.739   (  -2.064   -1.191   -5.553)    6.043
  12.739   (  -2.064   -1.191    5.553)    6.043
  13.608   (  -1.800   -1.040  -20.193)   20.299
  13.608   (  -1.800   -1.040   20.193)   20.299
======================= Grid point 968 (111/135) =======================
q-point: ( 0.13  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.852   (  15.640    9.030  -15.417)   23.745
   2.852   (  15.640    9.030   15.417)   23.745
   3.643   (  32.260   18.625  -18.131)   41.428
   3.643   (  32.260   18.625   18.131)   41.428
   6.307   (   8.164    4.714  -57.351)   58.120
   6.307   (   8.164    4.714   57.351)   58.120
  12.640   (  -6.669   -3.850   -7.949)   11.067
  12.640   (  -6.669   -3.850    7.949)   11.067
  12.644   (  -4.509   -2.603   -5.086)    7.279
  12.644   (  -4.509   -2.603    5.086)    7.279
  13.573   (  -1.206   -0.696  -17.132)   17.188
  13.573   (  -1.206   -0.696   17.132)   17.188
======================= Grid point 970 (112/135) =======================
q-point: ( 0.20  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.217   (  15.194    8.772  -13.576)   22.184
   3.217   (  15.194    8.772   13.576)   22.184
   4.149   (  10.557    6.095  -26.973)   29.600
   4.149   (  10.557    6.095   26.973)   29.600
   6.729   (  28.962   16.721  -39.471)   51.734
   6.729   (  28.962   16.721   39.471)   51.734
  12.455   (  -8.623   -4.978   -9.885)   14.031
  12.455   (  -8.623   -4.978    9.885)   14.031
  12.529   (  -5.201   -3.003   -4.918)    7.762
  12.529   (  -5.201   -3.003    4.918)    7.762
  13.533   (  -2.946   -1.701  -14.216)   14.618
  13.533   (  -2.946   -1.701   14.216)   14.618
======================= Grid point 972 (113/135) =======================
q-point: ( 0.27  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.519   (  10.584    6.111  -12.192)   17.263
   3.519   (  10.584    6.111   12.192)   17.263
   4.197   (  -3.834   -2.213  -29.960)   30.285
   4.197   (  -3.834   -2.213   29.960)   30.285
   7.554   (  40.052   23.124  -26.116)   53.113
   7.554   (  40.052   23.124   26.116)   53.113
  12.270   (  -7.162   -4.135  -10.961)   13.730
  12.270   (  -7.162   -4.135   10.961)   13.730
  12.416   (  -4.411   -2.547   -4.535)    6.820
  12.416   (  -4.411   -2.547    4.535)    6.820
  13.415   (  -7.506   -4.334  -12.144)   14.920
  13.415   (  -7.506   -4.334   12.144)   14.920
======================= Grid point 974 (114/135) =======================
q-point: ( 0.33  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.693   (   4.633    2.675  -11.207)   12.419
   3.693   (   4.633    2.675   11.207)   12.419
   4.060   (  -6.937   -4.005  -24.837)   26.097
   4.060   (  -6.937   -4.005   24.837)   26.097
   8.482   (  39.258   22.666  -20.180)   49.620
   8.482   (  39.258   22.666   20.180)   49.620
  12.125   (  -5.764   -3.328  -11.901)   13.636
  12.125   (  -5.764   -3.328   11.901)   13.636
  12.334   (  -2.658   -1.535   -4.374)    5.344
  12.334   (  -2.658   -1.535    4.374)    5.344
  13.193   ( -11.291   -6.519  -10.486)   16.732
  13.193   ( -11.291   -6.519   10.486)   16.732
======================= Grid point 976 (115/135) =======================
q-point: ( 0.40  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.746   (   0.467    0.270  -10.302)   10.316
   3.746   (   0.467    0.270   10.302)   10.316
   3.912   (  -5.401   -3.118  -17.047)   18.152
   3.912   (  -5.401   -3.118   17.047)   18.152
   9.313   (  31.807   18.364  -16.692)   40.342
   9.313   (  31.807   18.364   16.692)   40.342
  11.981   (  -7.223   -4.170  -12.172)   14.756
  11.981   (  -7.223   -4.170   12.172)   14.756
  12.292   (  -1.025   -0.592   -4.466)    4.620
  12.292   (  -1.025   -0.592    4.466)    4.620
  12.934   ( -10.061   -5.809   -8.602)   14.456
  12.934   ( -10.061   -5.809    8.602)   14.456
======================= Grid point 978 (116/135) =======================
q-point: ( 0.47  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.741   (  -0.402   -0.232   -9.577)    9.588
   3.741   (  -0.402   -0.232    9.577)    9.588
   3.825   (  -1.993   -1.150  -10.248)   10.504
   3.825   (  -1.993   -1.150   10.248)   10.504
   9.870   (  14.233    8.217  -12.229)   20.485
   9.870   (  14.233    8.217   12.229)   20.485
  11.811   (  -5.664   -3.270   -7.095)    9.649
  11.811   (  -5.664   -3.270    7.095)    9.649
  12.280   (  -0.194   -0.112   -4.585)    4.590
  12.280   (  -0.194   -0.112    4.585)    4.590
  12.773   (  -3.549   -2.049   -6.441)    7.634
  12.773   (  -3.549   -2.049    6.441)    7.634
======================= Grid point 996 (117/135) =======================
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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.796   (  10.521   18.223  -15.364)   26.054
   2.796   (  10.521   18.223   15.364)   26.054
   3.405   (  19.117   33.111  -16.577)   41.672
   3.405   (  19.117   33.111   16.577)   41.672
   6.268   (   2.272    3.935  -60.833)   61.002
   6.268   (   2.272    3.935   60.833)   61.002
  12.643   (  -3.648   -6.319   -1.853)    7.528
  12.643   (  -3.648   -6.319    1.853)    7.528
  12.708   (  -1.551   -2.686   -0.148)    3.105
  12.708   (  -1.551   -2.686    0.148)    3.105
  13.576   (  -1.088   -1.885  -18.017)   18.148
  13.576   (  -1.088   -1.885   18.017)   18.148
======================= Grid point 998 (118/135) =======================
q-point: ( 0.13  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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.199   (   9.775   24.391  -13.634)   29.603
   3.199   (   9.775   24.391   13.634)   29.603
   3.974   (  16.584   12.070  -23.219)   30.981
   3.974   (  16.584   12.070   23.219)   30.981
   6.504   (  15.923   16.227  -46.710)   51.949
   6.504   (  15.923   16.227   46.710)   51.949
  12.494   (  -4.951   -7.843   -5.449)   10.757
  12.494   (  -4.951   -7.843    5.449)   10.757
  12.623   (  -5.485   -1.867   -2.664)    6.377
  12.623   (  -5.485   -1.867    2.664)    6.377
  13.531   (  -1.272   -3.813  -15.073)   15.599
  13.531   (  -1.272   -3.813   15.073)   15.599
======================= Grid point 1000 (119/135) =======================
q-point: ( 0.20  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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.568   (   3.489   24.481  -12.530)   27.722
   3.568   (   3.489   24.481   12.530)   27.722
   4.184   (   3.667   -1.699  -29.147)   29.426
   4.184   (   3.667   -1.699   29.147)   29.426
   7.145   (  32.740   24.099  -30.844)   51.029
   7.145   (  32.740   24.099   30.844)   51.029
  12.335   (  -5.429   -5.528   -8.565)   11.549
  12.335   (  -5.429   -5.528    8.565)   11.549
  12.497   (  -6.663   -1.440   -4.107)    7.958
  12.497   (  -6.663   -1.440    4.107)    7.958
  13.437   (  -3.268   -8.152  -12.674)   15.419
  13.437   (  -3.268   -8.152   12.674)   15.419
======================= Grid point 1002 (120/135) =======================
q-point: ( 0.27  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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.777   (  -2.324   14.537  -12.668)   19.422
   3.777   (  -2.324   14.537   12.668)   19.422
   4.161   (  -4.747    3.088  -23.780)   24.445
   4.161   (  -4.747    3.088   23.780)   24.445
   8.019   (  39.141   22.828  -22.334)   50.517
   8.019   (  39.141   22.828   22.334)   50.517
  12.195   (  -6.019   -3.144  -11.051)   12.971
  12.195   (  -6.019   -3.144   11.051)   12.971
  12.386   (  -4.402   -0.767   -4.483)    6.330
  12.386   (  -4.402   -0.767    4.483)    6.330
  13.250   (  -6.556  -12.100  -10.879)   17.543
  13.250   (  -6.556  -12.100   10.879)   17.543
======================= Grid point 1004 (121/135) =======================
q-point: ( 0.33  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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.793   (  -3.096    0.844  -14.332)   14.687
   3.793   (  -3.096    0.844   14.332)   14.687
   4.118   ( -11.195   15.476  -13.097)   23.159
   4.118   ( -11.195   15.476   13.097)   23.159
   8.892   (  36.990   18.219  -18.149)   45.051
   8.892   (  36.990   18.219   18.149)   45.051
  12.049   (  -6.777   -4.721  -12.022)   14.586
  12.049   (  -6.777   -4.721   12.022)   14.586
  12.335   (  -2.215    1.681   -3.852)    4.750
  12.335   (  -2.215    1.681    3.852)    4.750
  12.990   (  -7.745  -13.302   -9.142)   17.903
  12.990   (  -7.745  -13.302    9.142)   17.903
======================= Grid point 1006 (122/135) =======================
q-point: ( 0.40  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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.738   (  -1.352   -2.640  -12.078)   12.437
   3.738   (  -1.352   -2.640   12.078)   12.437
   4.076   ( -13.861   20.040   -8.950)   25.958
   4.076   ( -13.861   20.040    8.950)   25.958
   9.604   (  26.595   10.647  -14.426)   32.074
   9.604   (  26.595   10.647   14.426)   32.074
  11.862   (  -7.009   -7.620  -10.297)   14.602
  11.862   (  -7.009   -7.620   10.297)   14.602
  12.329   (  -2.105    3.956   -4.096)    6.071
  12.329   (  -2.105    3.956    4.096)    6.071
  12.764   (  -3.340  -10.759   -7.100)   13.317
  12.764   (  -3.340  -10.759    7.100)   13.317
======================= Grid point 1008 (123/135) =======================
q-point: ( 0.47  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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.711   (   1.041   -1.803   -9.832)   10.050
   3.711   (   1.041   -1.803    9.832)   10.050
   4.054   ( -12.595   21.814   -8.118)   26.465
   4.054   ( -12.595   21.814    8.118)   26.465
   9.914   (   2.229   -3.861  -10.437)   11.349
   9.914   (   2.229   -3.861   10.437)   11.349
  11.738   (   2.031   -3.518   -3.363)    5.274
  11.738   (   2.031   -3.518    3.363)    5.274
  12.332   (  -2.717    4.706   -4.547)    7.085
  12.332   (  -2.717    4.706    4.547)    7.085
  12.681   (   3.915   -6.781   -5.810)    9.750
  12.681   (   3.915   -6.781    5.810)    9.750
======================= Grid point 1030 (124/135) =======================
q-point: ( 0.13  0.13 -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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.723   (  14.367   24.885  -11.217)   30.846
   3.723   (  14.367   24.885   11.217)   30.846
   4.098   (   1.103    1.911  -28.371)   28.457
   4.098   (   1.103    1.911   28.371)   28.457
   6.971   (  17.395   30.130  -33.578)   48.351
   6.971   (  17.395   30.130   33.578)   48.351
  12.352   (  -3.490   -6.044   -6.674)    9.657
  12.352   (  -3.490   -6.044    6.674)    9.657
  12.560   (  -2.583   -4.474   -3.735)    6.375
  12.560   (  -2.583   -4.474    3.735)    6.375
  13.414   (  -4.642   -8.039  -12.716)   15.744
  13.414   (  -4.642   -8.039   12.716)   15.744
======================= Grid point 1032 (125/135) =======================
q-point: ( 0.20  0.13 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.944   (  -4.757    1.839  -19.252)   19.916
   3.944   (  -4.757    1.839   19.252)   19.916
   4.272   (   6.829   20.548   -9.750)   23.747
   4.272   (   6.829   20.548    9.750)   23.747
   7.665   (  28.272   27.842  -24.464)   46.615
   7.665   (  28.272   27.842   24.464)   46.615
  12.251   (  -3.521   -2.988   -8.844)    9.977
  12.251   (  -3.521   -2.988    8.844)    9.977
  12.454   (  -4.718   -3.094   -4.947)    7.504
  12.454   (  -4.718   -3.094    4.947)    7.504
  13.214   (  -5.701  -13.603  -10.522)   18.118
  13.214   (  -5.701  -13.603   10.522)   18.118
======================= Grid point 1034 (126/135) =======================
q-point: ( 0.27  0.13 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.835   (  -4.144   -4.808  -20.490)   21.451
   3.835   (  -4.144   -4.808   20.490)   21.451
   4.524   (  -6.640   27.591   -8.271)   29.559
   4.524   (  -6.640   27.591    8.271)   29.559
   8.457   (  33.799   21.351  -19.330)   44.406
   8.457   (  33.799   21.351   19.330)   44.406
  12.132   (  -7.364   -3.766  -11.305)   14.008
  12.132   (  -7.364   -3.766   11.305)   14.008
  12.373   (  -1.425   -1.065   -4.048)    4.422
  12.373   (  -1.425   -1.065    4.048)    4.422
  12.959   (  -6.355  -15.743   -8.228)   18.866
  12.959   (  -6.355  -15.743    8.228)   18.866
======================= Grid point 1036 (127/135) =======================
q-point: ( 0.33  0.13 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.732   (  -2.137   -4.487  -15.453)   16.233
   3.732   (  -2.137   -4.487   15.453)   16.233
   4.606   ( -15.553   29.000   -7.629)   33.780
   4.606   ( -15.553   29.000    7.629)   33.780
   9.210   (  31.401   13.938  -15.796)   37.813
   9.210   (  31.401   13.938   15.796)   37.813
  11.921   (  -9.412   -8.597  -11.301)   17.036
  11.921   (  -9.412   -8.597   11.301)   17.036
  12.374   (   1.120    0.686   -1.910)    2.318
  12.374   (   1.120    0.686    1.910)    2.318
  12.710   (  -5.231  -12.498   -5.282)   14.542
  12.710   (  -5.231  -12.498    5.282)   14.542
======================= Grid point 1038 (128/135) =======================
q-point: ( 0.40  0.13 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.677   (   0.108   -3.008  -10.363)   10.792
   3.677   (   0.108   -3.008   10.363)   10.792
   4.602   ( -17.988   30.065   -7.149)   35.757
   4.602   ( -17.988   30.065    7.149)   35.757
   9.736   (  16.530    2.804  -11.219)   20.173
   9.736   (  16.530    2.804   11.219)   20.173
  11.690   (  -4.081   -9.328   -6.533)   12.097
  11.690   (  -4.081   -9.328    6.533)   12.097
  12.409   (   1.313    1.233   -3.800)    4.205
  12.409   (   1.313    1.233    3.800)    4.205
  12.552   (  -1.189   -7.420   -3.471)    8.277
  12.552   (  -1.189   -7.420    3.471)    8.277
======================= Grid point 1062 (129/135) =======================
q-point: ( 0.20  0.20 -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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.846   (  -3.917   -6.785  -22.414)   23.744
   3.846   (  -3.917   -6.785   22.414)   23.744
   4.746   (  12.202   21.134   -6.480)   25.249
   4.746   (  12.202   21.134    6.480)   25.249
   8.251   (  17.850   30.917  -19.988)   40.914
   8.251   (  17.850   30.917   19.988)   40.914
  12.190   (  -2.594   -4.493  -10.280)   11.515
  12.190   (  -2.594   -4.493   10.280)   11.515
  12.372   (  -2.650   -4.589   -5.258)    7.465
  12.372   (  -2.650   -4.589    5.258)    7.465
  12.930   (  -7.983  -13.826   -7.030)   17.444
  12.930   (  -7.983  -13.826    7.030)   17.444
======================= Grid point 1064 (130/135) =======================
q-point: ( 0.27  0.20 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.725   (  -2.772   -5.219  -17.497)   18.468
   3.725   (  -2.772   -5.219   17.497)   18.468
   5.053   (  -2.138   23.202   -6.487)   24.186
   5.053   (  -2.138   23.202    6.487)   24.186
   8.884   (  24.084   21.776  -16.547)   36.442
   8.884   (  24.084   21.776   16.547)   36.442
  12.006   ( -10.288  -10.355  -10.859)   18.193
  12.006   ( -10.288  -10.355   10.859)   18.193
  12.310   (   1.369   -5.390   -2.703)    6.183
  12.310   (   1.369   -5.390    2.703)    6.183
  12.706   (  -4.962   -7.629   -1.674)    9.253
  12.706   (  -4.962   -7.629    1.674)    9.253
======================= Grid point 1066 (131/135) =======================
q-point: ( 0.33  0.20 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.651   (  -0.884   -3.355  -11.473)   11.986
   3.651   (  -0.884   -3.355   11.473)   11.986
   5.162   ( -11.818   24.970   -6.272)   28.328
   5.162   ( -11.818   24.970    6.272)   28.328
   9.462   (  21.197   11.707  -12.769)   27.376
   9.462   (  21.197   11.707   12.769)   27.376
  11.690   (  -9.366  -14.444   -9.275)   19.555
  11.690   (  -9.366  -14.444    9.275)   19.555
  12.293   (   3.656   -7.269   -3.833)    8.994
  12.293   (   3.656   -7.269    3.833)    8.994
  12.599   (  -3.791    0.517   -2.148)    4.387
  12.599   (  -3.791    0.517    2.148)    4.387
======================= Grid point 1068 (132/135) =======================
q-point: ( 0.40  0.20 -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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.627   (   1.090   -1.887   -7.886)    8.181
   3.627   (   1.090   -1.887    7.886)    8.181
   5.179   ( -14.818   25.666   -6.139)   30.265
   5.179   ( -14.818   25.666    6.139)   30.265
   9.732   (   1.432   -2.480   -8.590)    9.055
   9.732   (   1.432   -2.480    8.590)    9.055
  11.510   (   4.677   -8.101   -4.163)   10.238
  11.510   (   4.677   -8.101    4.163)   10.238
  12.288   (   4.521   -7.831   -4.028)    9.899
  12.288   (   4.521   -7.831    4.028)    9.899
  12.571   (  -2.331    4.038   -3.135)    5.618
  12.571   (  -2.331    4.038    3.135)    5.618
======================= Grid point 1096 (133/135) =======================
q-point: ( 0.27  0.27 -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: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.642   (  -1.772   -3.070  -12.011)   12.523
   3.642   (  -1.772   -3.070   12.011)   12.523
   5.406   (   6.967   12.067   -5.818)   15.100
   5.406   (   6.967   12.067    5.818)   15.100
   9.322   (  12.436   21.540  -13.379)   28.242
   9.322   (  12.436   21.540   13.379)   28.242
  11.704   ( -10.445  -18.091   -9.813)   23.080
  11.704   ( -10.445  -18.091    9.813)   23.080
  12.217   (  -1.821   -3.155   -1.390)    3.899
  12.217   (  -1.821   -3.155    1.390)    3.899
  12.638   (  -0.519   -0.899   -2.325)    2.547
  12.638   (  -0.519   -0.899    2.325)    2.547
======================= Grid point 1098 (134/135) =======================
q-point: ( 0.33  0.27 -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: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.602   (  -0.203   -1.531   -6.118)    6.310
   3.602   (  -0.203   -1.531    6.118)    6.310
   5.551   (  -4.027   13.384   -5.523)   15.029
   5.551   (  -4.027   13.384    5.523)   15.029
   9.676   (   7.710    9.159   -8.149)   14.482
   9.676   (   7.710    9.159    8.149)   14.482
  11.392   (  -4.059  -13.473   -6.142)   15.353
  11.392   (  -4.059  -13.473    6.142)   15.353
  12.179   (   1.144   -3.679   -1.679)    4.203
  12.179   (   1.144   -3.679    1.679)    4.203
  12.632   (  -0.987    1.646   -2.822)    3.413
  12.632   (  -0.987    1.646    2.822)    3.413
======================= Grid point 1128 (135/135) =======================
q-point: ( 0.33  0.33 -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: 119
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.588   (  -0.000   -0.000    0.000)    0.000
   3.588   (  -0.000   -0.000    0.000)    0.000
   5.685   (  -0.000   -0.000   -5.266)    5.266
   5.685   (  -0.000   -0.000    5.266)    5.266
   9.780   (  -0.000   -0.000    0.000)    0.000
   9.780   (  -0.000   -0.000    0.000)    0.000
  11.241   (  -0.000   -0.000    0.000)    0.000
  11.241   (  -0.000   -0.000    0.000)    0.000
  12.143   (  -0.000   -0.000    0.000)    0.000
  12.143   (  -0.000   -0.000    0.000)    0.000
  12.649   (  -0.000   -0.000   -2.717)    2.717
  12.649   (  -0.000   -0.000    2.717)    2.717
=================== 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/21600
   10.0  16738.726  16738.726  16331.418     -0.000      0.000     -0.000 3/21600
   20.0   8204.962   8204.962   8204.316      0.000      0.000     -0.000 3/21600
   30.0   5106.057   5106.057   4711.299      0.000     -0.000     -0.000 3/21600
   40.0   3143.925   3143.925   2849.488     -0.000     -0.000     -0.000 3/21600
   50.0   2034.889   2034.889   1840.882     -0.000     -0.000     -0.000 3/21600
   60.0   1385.676   1385.676   1255.380     -0.000     -0.000     -0.000 3/21600
   70.0    985.383    985.383    893.250     -0.000     -0.000     -0.000 3/21600
   80.0    729.054    729.054    660.567     -0.000     -0.000     -0.000 3/21600
   90.0    559.607    559.607    506.687     -0.000     -0.000     -0.000 3/21600
  100.0    444.037    444.037    401.920     -0.000     -0.000     -0.000 3/21600
  110.0    362.716    362.716    328.392     -0.000     -0.000     -0.000 3/21600
  120.0    303.747    303.747    275.208     -0.000     -0.000     -0.000 3/21600
  130.0    259.771    259.771    235.627     -0.000     -0.000     -0.000 3/21600
  140.0    226.128    226.128    205.391     -0.000     -0.000     -0.000 3/21600
  150.0    199.794    199.794    181.745     -0.000     -0.000     -0.000 3/21600
  160.0    178.756    178.756    162.862     -0.000     -0.000     -0.000 3/21600
  170.0    161.640    161.640    147.500     -0.000     -0.000     -0.000 3/21600
  180.0    147.491    147.491    134.795     -0.000     -0.000     -0.000 3/21600
  190.0    135.625    135.625    124.134     -0.000     -0.000     -0.000 3/21600
  200.0    125.547    125.547    115.072     -0.000     -0.000     -0.000 3/21600
  210.0    116.892    116.892    107.281     -0.000     -0.000     -0.000 3/21600
  220.0    109.382    109.382    100.515     -0.000     -0.000     -0.000 3/21600
  230.0    102.808    102.808     94.584     -0.000     -0.000     -0.000 3/21600
  240.0     97.006     97.006     89.343     -0.000     -0.000     -0.000 3/21600
  250.0     91.849     91.849     84.679     -0.000     -0.000     -0.000 3/21600
  260.0     87.234     87.234     80.500     -0.000     -0.000     -0.000 3/21600
  270.0     83.079     83.079     76.734     -0.000     -0.000     -0.000 3/21600
  280.0     79.320     79.320     73.321     -0.000     -0.000     -0.000 3/21600
  290.0     75.901     75.901     70.213     -0.000     -0.000     -0.000 3/21600
  300.0     72.778     72.778     67.372     -0.000     -0.000     -0.000 3/21600
  310.0     69.913     69.913     64.762     -0.000     -0.000     -0.000 3/21600
  320.0     67.275     67.275     62.356     -0.000     -0.000     -0.000 3/21600
  330.0     64.838     64.838     60.131     -0.000     -0.000     -0.000 3/21600
  340.0     62.580     62.580     58.067     -0.000     -0.000     -0.000 3/21600
  350.0     60.480     60.480     56.146     -0.000     -0.000     -0.000 3/21600
  360.0     58.522     58.522     54.354     -0.000     -0.000     -0.000 3/21600
  370.0     56.693     56.693     52.677     -0.000     -0.000     -0.000 3/21600
  380.0     54.979     54.979     51.105     -0.000     -0.000     -0.000 3/21600
  390.0     53.370     53.370     49.629     -0.000     -0.000     -0.000 3/21600
  400.0     51.857     51.857     48.238     -0.000     -0.000     -0.000 3/21600
  410.0     50.430     50.430     46.926     -0.000     -0.000     -0.000 3/21600
  420.0     49.083     49.083     45.687     -0.000     -0.000     -0.000 3/21600
  430.0     47.808     47.808     44.513     -0.000     -0.000     -0.000 3/21600
  440.0     46.600     46.600     43.401     -0.000     -0.000     -0.000 3/21600
  450.0     45.454     45.454     42.344     -0.000     -0.000     -0.000 3/21600
  460.0     44.365     44.365     41.340     -0.000     -0.000     -0.000 3/21600
  470.0     43.329     43.329     40.383     -0.000     -0.000     -0.000 3/21600
  480.0     42.341     42.341     39.472     -0.000     -0.000     -0.000 3/21600
  490.0     41.399     41.399     38.601     -0.000     -0.000     -0.000 3/21600
  500.0     40.500     40.500     37.770     -0.000     -0.000     -0.000 3/21600
  510.0     39.640     39.640     36.974     -0.000     -0.000     -0.000 3/21600
  520.0     38.816     38.816     36.213     -0.000     -0.000     -0.000 3/21600
  530.0     38.027     38.027     35.483     -0.000     -0.000     -0.000 3/21600
  540.0     37.271     37.271     34.782     -0.000     -0.000     -0.000 3/21600
  550.0     36.545     36.545     34.110     -0.000     -0.000     -0.000 3/21600
  560.0     35.847     35.847     33.463     -0.000     -0.000     -0.000 3/21600
  570.0     35.176     35.176     32.841     -0.000     -0.000     -0.000 3/21600
  580.0     34.530     34.530     32.243     -0.000     -0.000     -0.000 3/21600
  590.0     33.909     33.909     31.666     -0.000     -0.000     -0.000 3/21600
  600.0     33.309     33.309     31.110     -0.000     -0.000     -0.000 3/21600
  610.0     32.731     32.731     30.574     -0.000     -0.000     -0.000 3/21600
  620.0     32.174     32.174     30.056     -0.000     -0.000     -0.000 3/21600
  630.0     31.635     31.635     29.556     -0.000     -0.000     -0.000 3/21600
  640.0     31.115     31.115     29.072     -0.000     -0.000     -0.000 3/21600
  650.0     30.611     30.611     28.605     -0.000     -0.000     -0.000 3/21600
  660.0     30.124     30.124     28.152     -0.000     -0.000     -0.000 3/21600
  670.0     29.653     29.653     27.714     -0.000     -0.000     -0.000 3/21600
  680.0     29.196     29.196     27.290     -0.000     -0.000     -0.000 3/21600
  690.0     28.754     28.754     26.878     -0.000     -0.000     -0.000 3/21600
  700.0     28.325     28.325     26.479     -0.000     -0.000     -0.000 3/21600
  710.0     27.909     27.909     26.092     -0.000     -0.000     -0.000 3/21600
  720.0     27.505     27.505     25.716     -0.000     -0.000     -0.000 3/21600
  730.0     27.113     27.113     25.352     -0.000     -0.000     -0.000 3/21600
  740.0     26.732     26.732     24.997     -0.000     -0.000     -0.000 3/21600
  750.0     26.362     26.362     24.652     -0.000     -0.000     -0.000 3/21600
  760.0     26.002     26.002     24.317     -0.000     -0.000     -0.000 3/21600
  770.0     25.652     25.652     23.991     -0.000     -0.000     -0.000 3/21600
  780.0     25.311     25.311     23.674     -0.000     -0.000     -0.000 3/21600
  790.0     24.980     24.980     23.365     -0.000     -0.000     -0.000 3/21600
  800.0     24.657     24.657     23.064     -0.000     -0.000     -0.000 3/21600
  810.0     24.342     24.342     22.771     -0.000     -0.000     -0.000 3/21600
  820.0     24.036     24.036     22.486     -0.000     -0.000     -0.000 3/21600
  830.0     23.737     23.737     22.207     -0.000     -0.000     -0.000 3/21600
  840.0     23.446     23.446     21.936     -0.000     -0.000     -0.000 3/21600
  850.0     23.162     23.162     21.671     -0.000     -0.000     -0.000 3/21600
  860.0     22.884     22.884     21.412     -0.000     -0.000     -0.000 3/21600
  870.0     22.614     22.614     21.160     -0.000     -0.000     -0.000 3/21600
  880.0     22.350     22.350     20.914     -0.000     -0.000     -0.000 3/21600
  890.0     22.092     22.092     20.673     -0.000     -0.000     -0.000 3/21600
  900.0     21.840     21.840     20.438     -0.000     -0.000     -0.000 3/21600
  910.0     21.593     21.593     20.208     -0.000     -0.000     -0.000 3/21600
  920.0     21.353     21.353     19.983     -0.000     -0.000     -0.000 3/21600
  930.0     21.117     21.117     19.764     -0.000     -0.000     -0.000 3/21600
  940.0     20.887     20.887     19.549     -0.000     -0.000     -0.000 3/21600
  950.0     20.662     20.662     19.339     -0.000     -0.000     -0.000 3/21600
  960.0     20.442     20.442     19.133     -0.000     -0.000     -0.000 3/21600
  970.0     20.226     20.226     18.932     -0.000     -0.000     -0.000 3/21600
  980.0     20.015     20.015     18.735     -0.000     -0.000     -0.000 3/21600
  990.0     19.808     19.808     18.542     -0.000     -0.000     -0.000 3/21600
 1000.0     19.606     19.606     18.353     -0.000     -0.000     -0.000 3/21600

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

