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

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

-------------------------[time 2026-01-09 04:43:34]-------------------------
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.898379686525482    0.000000000000000    0.000000000000000
  b   -1.949189843262741    3.376095842128284    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.433165160000000
Atomic positions (fractional):
   *1 Si  0.33333333333333  0.66666666666667  0.37461647148081  28.085
    2 Si  0.66666666666667  0.33333333333333  0.87461647148081  28.085
   *3 Ge  0.33333333333333  0.66666666666667  0.00038252851920  72.640
    4 Ge  0.66666666666667  0.33333333333333  0.50038252851920  72.640
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.898379686525482    0.000000000000000    0.000000000000000
  b   -1.949189843262741    3.376095842128284    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.433165160000000
Atomic positions (fractional):
   *1 Si  0.33333333333333  0.66666666666667  0.37461647148081  28.085 > 1
    2 Si  0.66666666666667  0.33333333333333  0.87461647148081  28.085 > 2
   *3 Ge  0.33333333333333  0.66666666666667  0.00038252851920  72.640 > 3
    4 Ge  0.66666666666667  0.33333333333333  0.50038252851920  72.640 > 4
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   15.593518746101928    0.000000000000000    0.000000000000000
  b   -7.796759373050964   13.504383368513135    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.866330319999999
Atomic positions (fractional):
   *1 Si  0.08333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    2 Si  0.33333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    3 Si  0.58333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    4 Si  0.83333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    5 Si  0.08333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    6 Si  0.33333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    7 Si  0.58333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    8 Si  0.83333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    9 Si  0.08333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   10 Si  0.33333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   11 Si  0.58333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   12 Si  0.83333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   13 Si  0.08333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   14 Si  0.33333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   15 Si  0.58333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   16 Si  0.83333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   17 Si  0.08333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   18 Si  0.33333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   19 Si  0.58333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   20 Si  0.83333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   21 Si  0.08333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   22 Si  0.33333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   23 Si  0.58333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   24 Si  0.83333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   25 Si  0.08333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   26 Si  0.33333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   27 Si  0.58333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   28 Si  0.83333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   29 Si  0.08333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   30 Si  0.33333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   31 Si  0.58333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   32 Si  0.83333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   33 Si  0.16666666666667  0.08333333333333  0.43730823574040  28.085 > 2
   34 Si  0.41666666666667  0.08333333333333  0.43730823574040  28.085 > 2
   35 Si  0.66666666666667  0.08333333333333  0.43730823574040  28.085 > 2
   36 Si  0.91666666666667  0.08333333333333  0.43730823574040  28.085 > 2
   37 Si  0.16666666666667  0.33333333333333  0.43730823574040  28.085 > 2
   38 Si  0.41666666666667  0.33333333333333  0.43730823574040  28.085 > 2
   39 Si  0.66666666666667  0.33333333333333  0.43730823574040  28.085 > 2
   40 Si  0.91666666666667  0.33333333333333  0.43730823574040  28.085 > 2
   41 Si  0.16666666666667  0.58333333333333  0.43730823574040  28.085 > 2
   42 Si  0.41666666666667  0.58333333333333  0.43730823574040  28.085 > 2
   43 Si  0.66666666666667  0.58333333333333  0.43730823574040  28.085 > 2
   44 Si  0.91666666666667  0.58333333333333  0.43730823574040  28.085 > 2
   45 Si  0.16666666666667  0.83333333333333  0.43730823574040  28.085 > 2
   46 Si  0.41666666666667  0.83333333333333  0.43730823574040  28.085 > 2
   47 Si  0.66666666666667  0.83333333333333  0.43730823574040  28.085 > 2
   48 Si  0.91666666666667  0.83333333333333  0.43730823574040  28.085 > 2
   49 Si  0.16666666666667  0.08333333333333  0.93730823574040  28.085 > 2
   50 Si  0.41666666666667  0.08333333333333  0.93730823574040  28.085 > 2
   51 Si  0.66666666666667  0.08333333333333  0.93730823574040  28.085 > 2
   52 Si  0.91666666666667  0.08333333333333  0.93730823574040  28.085 > 2
   53 Si  0.16666666666667  0.33333333333333  0.93730823574040  28.085 > 2
   54 Si  0.41666666666667  0.33333333333333  0.93730823574040  28.085 > 2
   55 Si  0.66666666666667  0.33333333333333  0.93730823574040  28.085 > 2
   56 Si  0.91666666666667  0.33333333333333  0.93730823574040  28.085 > 2
   57 Si  0.16666666666667  0.58333333333333  0.93730823574040  28.085 > 2
   58 Si  0.41666666666667  0.58333333333333  0.93730823574040  28.085 > 2
   59 Si  0.66666666666667  0.58333333333333  0.93730823574040  28.085 > 2
   60 Si  0.91666666666667  0.58333333333333  0.93730823574040  28.085 > 2
   61 Si  0.16666666666667  0.83333333333333  0.93730823574040  28.085 > 2
   62 Si  0.41666666666667  0.83333333333333  0.93730823574040  28.085 > 2
   63 Si  0.66666666666667  0.83333333333333  0.93730823574040  28.085 > 2
   64 Si  0.91666666666667  0.83333333333333  0.93730823574040  28.085 > 2
  *65 Ge  0.08333333333333  0.16666666666667  0.00019126425960  72.640 > 3
   66 Ge  0.33333333333333  0.16666666666667  0.00019126425960  72.640 > 3
   67 Ge  0.58333333333333  0.16666666666667  0.00019126425960  72.640 > 3
   68 Ge  0.83333333333333  0.16666666666667  0.00019126425960  72.640 > 3
   69 Ge  0.08333333333333  0.41666666666667  0.00019126425960  72.640 > 3
   70 Ge  0.33333333333333  0.41666666666667  0.00019126425960  72.640 > 3
   71 Ge  0.58333333333333  0.41666666666667  0.00019126425960  72.640 > 3
   72 Ge  0.83333333333333  0.41666666666667  0.00019126425960  72.640 > 3
   73 Ge  0.08333333333333  0.66666666666667  0.00019126425960  72.640 > 3
   74 Ge  0.33333333333333  0.66666666666667  0.00019126425960  72.640 > 3
   75 Ge  0.58333333333333  0.66666666666667  0.00019126425960  72.640 > 3
   76 Ge  0.83333333333333  0.66666666666667  0.00019126425960  72.640 > 3
   77 Ge  0.08333333333333  0.91666666666667  0.00019126425960  72.640 > 3
   78 Ge  0.33333333333333  0.91666666666667  0.00019126425960  72.640 > 3
   79 Ge  0.58333333333333  0.91666666666667  0.00019126425960  72.640 > 3
   80 Ge  0.83333333333333  0.91666666666667  0.00019126425960  72.640 > 3
   81 Ge  0.08333333333333  0.16666666666667  0.50019126425960  72.640 > 3
   82 Ge  0.33333333333333  0.16666666666667  0.50019126425960  72.640 > 3
   83 Ge  0.58333333333333  0.16666666666667  0.50019126425960  72.640 > 3
   84 Ge  0.83333333333333  0.16666666666667  0.50019126425960  72.640 > 3
   85 Ge  0.08333333333333  0.41666666666667  0.50019126425960  72.640 > 3
   86 Ge  0.33333333333333  0.41666666666667  0.50019126425960  72.640 > 3
   87 Ge  0.58333333333333  0.41666666666667  0.50019126425960  72.640 > 3
   88 Ge  0.83333333333333  0.41666666666667  0.50019126425960  72.640 > 3
   89 Ge  0.08333333333333  0.66666666666667  0.50019126425960  72.640 > 3
   90 Ge  0.33333333333333  0.66666666666667  0.50019126425960  72.640 > 3
   91 Ge  0.58333333333333  0.66666666666667  0.50019126425960  72.640 > 3
   92 Ge  0.83333333333333  0.66666666666667  0.50019126425960  72.640 > 3
   93 Ge  0.08333333333333  0.91666666666667  0.50019126425960  72.640 > 3
   94 Ge  0.33333333333333  0.91666666666667  0.50019126425960  72.640 > 3
   95 Ge  0.58333333333333  0.91666666666667  0.50019126425960  72.640 > 3
   96 Ge  0.83333333333333  0.91666666666667  0.50019126425960  72.640 > 3
   97 Ge  0.16666666666667  0.08333333333333  0.25019126425960  72.640 > 4
   98 Ge  0.41666666666667  0.08333333333333  0.25019126425960  72.640 > 4
   99 Ge  0.66666666666667  0.08333333333333  0.25019126425960  72.640 > 4
  100 Ge  0.91666666666667  0.08333333333333  0.25019126425960  72.640 > 4
  101 Ge  0.16666666666667  0.33333333333333  0.25019126425960  72.640 > 4
  102 Ge  0.41666666666667  0.33333333333333  0.25019126425960  72.640 > 4
  103 Ge  0.66666666666667  0.33333333333333  0.25019126425960  72.640 > 4
  104 Ge  0.91666666666667  0.33333333333333  0.25019126425960  72.640 > 4
  105 Ge  0.16666666666667  0.58333333333333  0.25019126425960  72.640 > 4
  106 Ge  0.41666666666667  0.58333333333333  0.25019126425960  72.640 > 4
  107 Ge  0.66666666666667  0.58333333333333  0.25019126425960  72.640 > 4
  108 Ge  0.91666666666667  0.58333333333333  0.25019126425960  72.640 > 4
  109 Ge  0.16666666666667  0.83333333333333  0.25019126425960  72.640 > 4
  110 Ge  0.41666666666667  0.83333333333333  0.25019126425960  72.640 > 4
  111 Ge  0.66666666666667  0.83333333333333  0.25019126425960  72.640 > 4
  112 Ge  0.91666666666667  0.83333333333333  0.25019126425960  72.640 > 4
  113 Ge  0.16666666666667  0.08333333333333  0.75019126425960  72.640 > 4
  114 Ge  0.41666666666667  0.08333333333333  0.75019126425960  72.640 > 4
  115 Ge  0.66666666666667  0.08333333333333  0.75019126425960  72.640 > 4
  116 Ge  0.91666666666667  0.08333333333333  0.75019126425960  72.640 > 4
  117 Ge  0.16666666666667  0.33333333333333  0.75019126425960  72.640 > 4
  118 Ge  0.41666666666667  0.33333333333333  0.75019126425960  72.640 > 4
  119 Ge  0.66666666666667  0.33333333333333  0.75019126425960  72.640 > 4
  120 Ge  0.91666666666667  0.33333333333333  0.75019126425960  72.640 > 4
  121 Ge  0.16666666666667  0.58333333333333  0.75019126425960  72.640 > 4
  122 Ge  0.41666666666667  0.58333333333333  0.75019126425960  72.640 > 4
  123 Ge  0.66666666666667  0.58333333333333  0.75019126425960  72.640 > 4
  124 Ge  0.91666666666667  0.58333333333333  0.75019126425960  72.640 > 4
  125 Ge  0.16666666666667  0.83333333333333  0.75019126425960  72.640 > 4
  126 Ge  0.41666666666667  0.83333333333333  0.75019126425960  72.640 > 4
  127 Ge  0.66666666666667  0.83333333333333  0.75019126425960  72.640 > 4
  128 Ge  0.91666666666667  0.83333333333333  0.75019126425960  72.640 > 4
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           18.9308102    0.0000000    0.0000000
            0.0000000   18.9308102    0.0000000
            0.0000000    0.0000000   17.3356274
-------------------------- Born effective charges --------------------------
    1 Si    0.1134643    0.0000000    0.0000000
            0.0000000    0.1134643    0.0000000
            0.0000000    0.0000000    0.1443076
    2 Si    0.1134643    0.0000000    0.0000000
            0.0000000    0.1134643    0.0000000
            0.0000000    0.0000000    0.1443076
    3 Ge   -0.1134643   -0.0000000    0.0000000
            0.0000000   -0.1134643    0.0000000
            0.0000000    0.0000000   -0.1443076
    4 Ge   -0.1134643   -0.0000000    0.0000000
            0.0000000   -0.1134643    0.0000000
            0.0000000    0.0000000   -0.1443076
----------------------------------------------------------------------------
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.159
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.170
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) 
Permutation basis: 384/384
Permutation basis: 2856/2856
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 224
Number of blocks in projector: 156
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 98
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 80
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 46
Use standard eigh solver.
Tree of FC basis block matrices:
- (224, 220), data: False
|-- (46, 44), data: True
|-- (80, 80), data: True
|-- (98, 96), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-09 04:43:39]-------------------------
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.898379686525482    0.000000000000000    0.000000000000000
  b   -1.949189843262741    3.376095842128284    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.433165160000000
Atomic positions (fractional):
    1 Si  0.33333333333333  0.66666666666667  0.37461647148081  28.085
    2 Si  0.66666666666667  0.33333333333333  0.87461647148081  28.085
    3 Ge  0.33333333333333  0.66666666666667  0.00038252851920  72.640
    4 Ge  0.66666666666667  0.33333333333333  0.50038252851920  72.640
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.593518746101928    0.000000000000000    0.000000000000000
  b   -7.796759373050964   13.504383368513135    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.866330319999999
Atomic positions (fractional):
    1 Si  0.08333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    2 Si  0.33333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    3 Si  0.58333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    4 Si  0.83333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    5 Si  0.08333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    6 Si  0.33333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    7 Si  0.58333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    8 Si  0.83333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    9 Si  0.08333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   10 Si  0.33333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   11 Si  0.58333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   12 Si  0.83333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   13 Si  0.08333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   14 Si  0.33333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   15 Si  0.58333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   16 Si  0.83333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   17 Si  0.08333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   18 Si  0.33333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   19 Si  0.58333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   20 Si  0.83333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   21 Si  0.08333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   22 Si  0.33333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   23 Si  0.58333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   24 Si  0.83333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   25 Si  0.08333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   26 Si  0.33333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   27 Si  0.58333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   28 Si  0.83333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   29 Si  0.08333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   30 Si  0.33333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   31 Si  0.58333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   32 Si  0.83333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   33 Si  0.16666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   34 Si  0.41666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   35 Si  0.66666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   36 Si  0.91666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   37 Si  0.16666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   38 Si  0.41666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   39 Si  0.66666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   40 Si  0.91666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   41 Si  0.16666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   42 Si  0.41666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   43 Si  0.66666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   44 Si  0.91666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   45 Si  0.16666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   46 Si  0.41666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   47 Si  0.66666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   48 Si  0.91666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   49 Si  0.16666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   50 Si  0.41666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   51 Si  0.66666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   52 Si  0.91666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   53 Si  0.16666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   54 Si  0.41666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   55 Si  0.66666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   56 Si  0.91666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   57 Si  0.16666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   58 Si  0.41666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   59 Si  0.66666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   60 Si  0.91666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   61 Si  0.16666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   62 Si  0.41666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   63 Si  0.66666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   64 Si  0.91666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   65 Ge  0.08333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   66 Ge  0.33333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   67 Ge  0.58333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   68 Ge  0.83333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   69 Ge  0.08333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   70 Ge  0.33333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   71 Ge  0.58333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   72 Ge  0.83333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   73 Ge  0.08333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   74 Ge  0.33333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   75 Ge  0.58333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   76 Ge  0.83333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   77 Ge  0.08333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   78 Ge  0.33333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   79 Ge  0.58333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   80 Ge  0.83333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   81 Ge  0.08333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   82 Ge  0.33333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   83 Ge  0.58333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   84 Ge  0.83333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   85 Ge  0.08333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   86 Ge  0.33333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   87 Ge  0.58333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   88 Ge  0.83333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   89 Ge  0.08333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   90 Ge  0.33333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   91 Ge  0.58333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   92 Ge  0.83333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   93 Ge  0.08333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   94 Ge  0.33333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   95 Ge  0.58333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   96 Ge  0.83333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   97 Ge  0.16666666666667  0.08333333333333  0.25019126425960  72.640 > 97
   98 Ge  0.41666666666667  0.08333333333333  0.25019126425960  72.640 > 97
   99 Ge  0.66666666666667  0.08333333333333  0.25019126425960  72.640 > 97
  100 Ge  0.91666666666667  0.08333333333333  0.25019126425960  72.640 > 97
  101 Ge  0.16666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  102 Ge  0.41666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  103 Ge  0.66666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  104 Ge  0.91666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  105 Ge  0.16666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  106 Ge  0.41666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  107 Ge  0.66666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  108 Ge  0.91666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  109 Ge  0.16666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  110 Ge  0.41666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  111 Ge  0.66666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  112 Ge  0.91666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  113 Ge  0.16666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  114 Ge  0.41666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  115 Ge  0.66666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  116 Ge  0.91666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  117 Ge  0.16666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  118 Ge  0.41666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  119 Ge  0.66666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  120 Ge  0.91666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  121 Ge  0.16666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  122 Ge  0.41666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  123 Ge  0.66666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  124 Ge  0.91666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  125 Ge  0.16666666666667  0.83333333333333  0.75019126425960  72.640 > 97
  126 Ge  0.41666666666667  0.83333333333333  0.75019126425960  72.640 > 97
  127 Ge  0.66666666666667  0.83333333333333  0.75019126425960  72.640 > 97
  128 Ge  0.91666666666667  0.83333333333333  0.75019126425960  72.640 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           18.9308102    0.0000000    0.0000000
            0.0000000   18.9308102    0.0000000
            0.0000000    0.0000000   17.3356274
-------------------------- Born effective charges --------------------------
    1 Si    0.1134643    0.0000000    0.0000000
            0.0000000    0.1134643    0.0000000
            0.0000000    0.0000000    0.1443076
    2 Si    0.1134643    0.0000000    0.0000000
            0.0000000    0.1134643    0.0000000
            0.0000000    0.0000000    0.1443076
    3 Ge   -0.1134643   -0.0000000    0.0000000
            0.0000000   -0.1134643    0.0000000
            0.0000000    0.0000000   -0.1443076
    4 Ge   -0.1134643   -0.0000000    0.0000000
            0.0000000   -0.1134643    0.0000000
            0.0000000    0.0000000   -0.1443076
----------------------------------------------------------------------------
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.00000102 (zzz) -0.00000102 (zzz) -0.00000102 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-09 04:43:44]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-09 04:43:44]-------------------------
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.898379686525482    0.000000000000000    0.000000000000000
  b   -1.949189843262741    3.376095842128284    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.433165160000000
Atomic positions (fractional):
    1 Si  0.33333333333333  0.66666666666667  0.37461647148081  28.085
    2 Si  0.66666666666667  0.33333333333333  0.87461647148081  28.085
    3 Ge  0.33333333333333  0.66666666666667  0.00038252851920  72.640
    4 Ge  0.66666666666667  0.33333333333333  0.50038252851920  72.640
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.593518746101928    0.000000000000000    0.000000000000000
  b   -7.796759373050964   13.504383368513135    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.866330319999999
Atomic positions (fractional):
    1 Si  0.08333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    2 Si  0.33333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    3 Si  0.58333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    4 Si  0.83333333333333  0.16666666666667  0.18730823574040  28.085 > 1
    5 Si  0.08333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    6 Si  0.33333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    7 Si  0.58333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    8 Si  0.83333333333333  0.41666666666667  0.18730823574040  28.085 > 1
    9 Si  0.08333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   10 Si  0.33333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   11 Si  0.58333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   12 Si  0.83333333333333  0.66666666666667  0.18730823574040  28.085 > 1
   13 Si  0.08333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   14 Si  0.33333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   15 Si  0.58333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   16 Si  0.83333333333333  0.91666666666667  0.18730823574040  28.085 > 1
   17 Si  0.08333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   18 Si  0.33333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   19 Si  0.58333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   20 Si  0.83333333333333  0.16666666666667  0.68730823574040  28.085 > 1
   21 Si  0.08333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   22 Si  0.33333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   23 Si  0.58333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   24 Si  0.83333333333333  0.41666666666667  0.68730823574040  28.085 > 1
   25 Si  0.08333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   26 Si  0.33333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   27 Si  0.58333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   28 Si  0.83333333333333  0.66666666666667  0.68730823574040  28.085 > 1
   29 Si  0.08333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   30 Si  0.33333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   31 Si  0.58333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   32 Si  0.83333333333333  0.91666666666667  0.68730823574040  28.085 > 1
   33 Si  0.16666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   34 Si  0.41666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   35 Si  0.66666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   36 Si  0.91666666666667  0.08333333333333  0.43730823574040  28.085 > 33
   37 Si  0.16666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   38 Si  0.41666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   39 Si  0.66666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   40 Si  0.91666666666667  0.33333333333333  0.43730823574040  28.085 > 33
   41 Si  0.16666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   42 Si  0.41666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   43 Si  0.66666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   44 Si  0.91666666666667  0.58333333333333  0.43730823574040  28.085 > 33
   45 Si  0.16666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   46 Si  0.41666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   47 Si  0.66666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   48 Si  0.91666666666667  0.83333333333333  0.43730823574040  28.085 > 33
   49 Si  0.16666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   50 Si  0.41666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   51 Si  0.66666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   52 Si  0.91666666666667  0.08333333333333  0.93730823574040  28.085 > 33
   53 Si  0.16666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   54 Si  0.41666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   55 Si  0.66666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   56 Si  0.91666666666667  0.33333333333333  0.93730823574040  28.085 > 33
   57 Si  0.16666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   58 Si  0.41666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   59 Si  0.66666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   60 Si  0.91666666666667  0.58333333333333  0.93730823574040  28.085 > 33
   61 Si  0.16666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   62 Si  0.41666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   63 Si  0.66666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   64 Si  0.91666666666667  0.83333333333333  0.93730823574040  28.085 > 33
   65 Ge  0.08333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   66 Ge  0.33333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   67 Ge  0.58333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   68 Ge  0.83333333333333  0.16666666666667  0.00019126425960  72.640 > 65
   69 Ge  0.08333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   70 Ge  0.33333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   71 Ge  0.58333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   72 Ge  0.83333333333333  0.41666666666667  0.00019126425960  72.640 > 65
   73 Ge  0.08333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   74 Ge  0.33333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   75 Ge  0.58333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   76 Ge  0.83333333333333  0.66666666666667  0.00019126425960  72.640 > 65
   77 Ge  0.08333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   78 Ge  0.33333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   79 Ge  0.58333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   80 Ge  0.83333333333333  0.91666666666667  0.00019126425960  72.640 > 65
   81 Ge  0.08333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   82 Ge  0.33333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   83 Ge  0.58333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   84 Ge  0.83333333333333  0.16666666666667  0.50019126425960  72.640 > 65
   85 Ge  0.08333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   86 Ge  0.33333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   87 Ge  0.58333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   88 Ge  0.83333333333333  0.41666666666667  0.50019126425960  72.640 > 65
   89 Ge  0.08333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   90 Ge  0.33333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   91 Ge  0.58333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   92 Ge  0.83333333333333  0.66666666666667  0.50019126425960  72.640 > 65
   93 Ge  0.08333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   94 Ge  0.33333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   95 Ge  0.58333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   96 Ge  0.83333333333333  0.91666666666667  0.50019126425960  72.640 > 65
   97 Ge  0.16666666666667  0.08333333333333  0.25019126425960  72.640 > 97
   98 Ge  0.41666666666667  0.08333333333333  0.25019126425960  72.640 > 97
   99 Ge  0.66666666666667  0.08333333333333  0.25019126425960  72.640 > 97
  100 Ge  0.91666666666667  0.08333333333333  0.25019126425960  72.640 > 97
  101 Ge  0.16666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  102 Ge  0.41666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  103 Ge  0.66666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  104 Ge  0.91666666666667  0.33333333333333  0.25019126425960  72.640 > 97
  105 Ge  0.16666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  106 Ge  0.41666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  107 Ge  0.66666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  108 Ge  0.91666666666667  0.58333333333333  0.25019126425960  72.640 > 97
  109 Ge  0.16666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  110 Ge  0.41666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  111 Ge  0.66666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  112 Ge  0.91666666666667  0.83333333333333  0.25019126425960  72.640 > 97
  113 Ge  0.16666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  114 Ge  0.41666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  115 Ge  0.66666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  116 Ge  0.91666666666667  0.08333333333333  0.75019126425960  72.640 > 97
  117 Ge  0.16666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  118 Ge  0.41666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  119 Ge  0.66666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  120 Ge  0.91666666666667  0.33333333333333  0.75019126425960  72.640 > 97
  121 Ge  0.16666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  122 Ge  0.41666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  123 Ge  0.66666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  124 Ge  0.91666666666667  0.58333333333333  0.75019126425960  72.640 > 97
  125 Ge  0.16666666666667  0.83333333333333  0.75019126425960  72.640 > 97
  126 Ge  0.41666666666667  0.83333333333333  0.75019126425960  72.640 > 97
  127 Ge  0.66666666666667  0.83333333333333  0.75019126425960  72.640 > 97
  128 Ge  0.91666666666667  0.83333333333333  0.75019126425960  72.640 > 97
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           18.9308102    0.0000000    0.0000000
            0.0000000   18.9308102    0.0000000
            0.0000000    0.0000000   17.3356274
-------------------------- Born effective charges --------------------------
    1 Si    0.1134643    0.0000000    0.0000000
            0.0000000    0.1134643    0.0000000
            0.0000000    0.0000000    0.1443076
    2 Si    0.1134643    0.0000000    0.0000000
            0.0000000    0.1134643    0.0000000
            0.0000000    0.0000000    0.1443076
    3 Ge   -0.1134643   -0.0000000    0.0000000
            0.0000000   -0.1134643    0.0000000
            0.0000000    0.0000000   -0.1443076
    4 Ge   -0.1134643   -0.0000000    0.0000000
            0.0000000   -0.1134643    0.0000000
            0.0000000    0.0000000   -0.1443076
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000102 (zzz) -0.00000102 (zzz) -0.00000102 (zzz)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 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.95, Number of G-points: 313, Lambda: 0.42
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: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
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
   1.617   (   0.000    0.000   -0.000)    0.000
   1.617   (   0.000   -0.000    0.000)    0.000
   6.886   (   0.000    0.000    0.000)    0.000
  11.476   (   0.000   -0.000    0.000)    0.000
  12.017   (  -0.000    0.000   -0.000)    0.000
  12.017   (  -0.000    0.000    0.000)    0.000
  12.098   (   0.000    0.000   -0.000)    0.000
  12.337   (   0.000   -0.000   -0.000)    0.000
  12.337   (   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: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.605   (  26.335   15.205    0.000)   30.409
   0.638   (  27.618   15.945    0.000)   31.891
   1.246   (  53.963   31.156    0.000)   62.312
   1.739   (  10.121    5.843    0.000)   11.686
   2.071   (  34.524   19.933    0.000)   39.865
   6.812   (  -6.395   -3.692    0.000)    7.385
  11.345   ( -10.633   -6.139    0.000)   12.278
  11.969   (  -4.175   -2.410    0.000)    4.821
  12.022   (  -2.731   -1.577    0.000)    3.154
  12.070   (   0.040    0.023    0.000)    0.046
  12.292   (  -3.848   -2.222    0.000)    4.443
  12.317   (  -1.773   -1.024    0.000)    2.047
======================= Grid point 2 (3/135) =======================
q-point: ( 0.13  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.190   (  24.587   14.195    0.000)   28.391
   1.257   (  26.617   15.367    0.000)   30.735
   2.036   (  14.725    8.502    0.000)   17.003
   2.444   (  50.738   29.293    0.000)   58.587
   2.964   (  40.413   23.333    0.000)   46.665
   6.607   ( -11.038   -6.373    0.000)   12.746
  11.039   ( -15.212   -8.783    0.000)   17.565
  11.835   (  -7.275   -4.200    0.000)    8.401
  11.883   (  -8.038   -4.641    0.000)    9.281
  12.083   (  -0.752   -0.434    0.000)    0.868
  12.168   (  -6.879   -3.972    0.000)    7.943
  12.268   (  -2.120   -1.224    0.000)    2.448
======================= Grid point 3 (4/135) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.704   (  19.873   11.474    0.000)   22.948
   1.848   (  25.100   14.492    0.000)   28.983
   2.356   (  12.154    7.017    0.000)   14.035
   3.541   (  45.053   26.011    0.000)   52.023
   3.803   (  31.208   18.018    0.000)   36.036
   6.349   ( -10.268   -5.928    0.000)   11.856
  10.695   ( -14.064   -8.120    0.000)   16.240
  11.652   (  -8.370   -4.833    0.000)    9.665
  11.683   (  -9.111   -5.260    0.000)   10.520
  11.994   (  -7.956   -4.593    0.000)    9.187
  12.005   (  -6.007   -3.468    0.000)    6.937
  12.232   (  -1.222   -0.705    0.000)    1.411
======================= Grid point 4 (5/135) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.071   (  11.880    6.859    0.000)   13.718
   2.394   (  22.652   13.078    0.000)   26.156
   2.541   (   3.281    1.894    0.000)    3.789
   4.306   (  12.010    6.934    0.000)   13.868
   4.482   (  36.972   21.346    0.000)   42.692
   6.213   (  -0.445   -0.257    0.000)    0.514
  10.432   (  -8.488   -4.900    0.000)    9.801
  11.472   (  -7.078   -4.086    0.000)    8.173
  11.487   (  -7.771   -4.487    0.000)    8.974
  11.824   (  -9.601   -5.543    0.000)   11.086
  11.826   (  -6.335   -3.657    0.000)    7.315
  12.199   (  -1.852   -1.070    0.000)    2.139
======================= Grid point 5 (6/135) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.238   (   2.967    1.713    0.000)    3.426
   2.485   (  -8.087   -4.669    0.000)    9.338
   2.869   (  18.736   10.817    0.000)   21.634
   4.385   (  -2.754   -1.590    0.000)    3.181
   5.217   (  27.355   15.793    0.000)   31.587
   6.309   (   6.919    3.995    0.000)    7.990
  10.323   (  -0.955   -0.551    0.000)    1.103
  11.333   (  -5.900   -3.406    0.000)    6.813
  11.342   (  -4.229   -2.441    0.000)    4.883
  11.584   ( -10.823   -6.249    0.000)   12.497
  11.719   (  -2.969   -1.714    0.000)    3.428
  12.149   (  -2.395   -1.383    0.000)    2.765
======================= Grid point 6 (7/135) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.208   ( -14.064   -8.120    0.000)   16.240
   2.234   (  -3.087   -1.782    0.000)    3.564
   3.234   (  12.910    7.454    0.000)   14.907
   4.281   (  -4.911   -2.835    0.000)    5.671
   5.725   (  17.164    9.910    0.000)   19.819
   6.448   (   4.400    2.540    0.000)    5.080
  10.381   (   5.659    3.267    0.000)    6.534
  11.173   (  -9.448   -5.455    0.000)   10.910
  11.277   (  -1.602   -0.925    0.000)    1.850
  11.385   (  -5.198   -3.001    0.000)    6.002
  11.685   (  -0.294   -0.170    0.000)    0.339
  12.097   (  -2.014   -1.163    0.000)    2.326
======================= Grid point 7 (8/135) =======================
q-point: ( 0.47  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.921   (  -8.039   -4.641    0.000)    9.283
   2.170   (  -1.952   -1.127    0.000)    2.255
   3.441   (   4.814    2.780    0.000)    5.559
   4.196   (  -2.087   -1.205    0.000)    2.410
   5.995   (   6.208    3.584    0.000)    7.168
   6.503   (   0.763    0.440    0.000)    0.881
  10.546   (   7.246    4.183    0.000)    8.367
  10.937   (  -9.016   -5.206    0.000)   10.411
  11.258   (  -0.295   -0.170    0.000)    0.341
  11.329   (  -0.757   -0.437    0.000)    0.874
  11.688   (   0.266    0.154    0.000)    0.307
  12.064   (  -0.787   -0.455    0.000)    0.909
======================= Grid point 16 (9/135) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.087   (  15.297   26.496    0.000)   30.594
   1.099   (  17.149   29.703    0.000)   34.298
   1.983   (  10.707   18.545    0.000)   21.414
   2.105   (  28.164   48.782    0.000)   56.329
   2.690   (  22.399   38.796    0.000)   44.798
   6.672   (  -5.851  -10.134    0.000)   11.702
  11.131   (  -8.377  -14.509    0.000)   16.754
  11.848   (  -4.916   -8.514    0.000)    9.831
  11.999   (  -2.225   -3.855    0.000)    4.451
  12.039   (  -1.103   -1.911    0.000)    2.207
  12.236   (  -1.921   -3.327    0.000)    3.842
  12.258   (  -2.665   -4.616    0.000)    5.330
======================= Grid point 17 (10/135) =======================
q-point: ( 0.13  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.631   (  19.722   21.641    0.000)   29.279
   1.675   (  13.382   31.432    0.000)   34.162
   2.363   (   8.589   22.923    0.000)   24.479
   3.101   (  37.535   35.313    0.000)   51.535
   3.475   (  28.775   26.353    0.000)   39.019
   6.436   (  -8.760  -10.261    0.000)   13.492
  10.815   ( -11.460  -13.057    0.000)   17.373
  11.661   (  -5.578  -10.107    0.000)   11.545
  11.855   ( -10.205   -1.354    0.000)   10.294
  11.985   (  -0.434   -8.034    0.000)    8.045
  12.104   (  -7.957   -4.601    0.000)    9.192
  12.207   (   0.821   -4.061    0.000)    4.143
======================= Grid point 18 (11/135) =======================
q-point: ( 0.20  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.141   (   4.913   30.175    0.000)   30.572
   2.172   (  20.153   17.970    0.000)   27.001
   2.676   (  -0.364   23.989    0.000)   23.992
   4.046   (  37.685   24.191    0.000)   44.782
   4.105   (  20.817   11.577    0.000)   23.819
   6.228   (  -4.436   -5.685    0.000)    7.211
  10.518   (  -9.433   -9.344    0.000)   13.278
  11.474   (  -4.766   -9.223    0.000)   10.382
  11.651   ( -11.127   -1.886    0.000)   11.286
  11.872   (  -3.552   -7.046    0.000)    7.891
  11.898   (  -8.391   -8.007    0.000)   11.598
  12.180   (   0.444   -4.204    0.000)    4.227
======================= Grid point 19 (12/135) =======================
q-point: ( 0.27  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.418   (  -5.191   26.296    0.000)   26.804
   2.663   (  17.556   14.437    0.000)   22.730
   2.794   ( -12.210   22.044    0.000)   25.200
   4.370   (   4.395    0.052    0.000)    4.395
   4.851   (  32.652   15.658    0.000)   36.212
   6.215   (   5.614    0.151    0.000)    5.616
  10.330   (  -3.548   -5.423    0.000)    6.481
  11.321   (  -3.262   -8.180    0.000)    8.806
  11.458   (  -8.750   -1.997    0.000)    8.975
  11.665   (  -8.014   -9.540    0.000)   12.460
  11.753   (  -3.557   -4.006    0.000)    5.357
  12.137   (  -0.188   -4.819    0.000)    4.823
======================= Grid point 20 (13/135) =======================
q-point: ( 0.33  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.494   ( -13.346   22.508    0.000)   26.167
   2.663   ( -25.228   21.743    0.000)   33.305
   3.086   (  14.705   10.995    0.000)   18.361
   4.337   (  -4.114   -2.648    0.000)    4.893
   5.461   (  24.319    8.897    0.000)   25.896
   6.351   (   7.969   -0.051    0.000)    7.969
  10.296   (   4.099   -2.570    0.000)    4.838
  11.172   (  -3.773  -12.270    0.000)   12.837
  11.323   (  -4.807   -1.719    0.000)    5.105
  11.466   (  -6.198   -3.820    0.000)    7.281
  11.684   (  -0.884   -2.203    0.000)    2.374
  12.084   (  -0.248   -4.702    0.000)    4.709
======================= Grid point 21 (14/135) =======================
q-point: ( 0.40  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.414   ( -29.102   27.257    0.000)   39.873
   2.449   ( -16.563   21.793    0.000)   27.373
   3.377   (   8.289    7.084    0.000)   10.903
   4.237   (  -4.088   -1.525    0.000)    4.363
   5.852   (  14.291    3.148    0.000)   14.633
   6.444   (   4.473   -3.001    0.000)    5.386
  10.408   (  10.272   -0.882    0.000)   10.310
  10.942   (  -5.442  -16.405    0.000)   17.284
  11.256   (  -1.473   -1.393    0.000)    2.028
  11.388   (  -2.956    2.054    0.000)    3.600
  11.669   (   1.027   -1.839    0.000)    2.106
  12.042   (   0.438   -3.860    0.000)    3.885
======================= Grid point 22 (15/135) =======================
q-point: ( 0.47  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.282   ( -21.068   36.490    0.000)   42.135
   2.410   ( -13.587   23.533    0.000)   27.174
   3.487   (  -1.129    1.955    0.000)    2.258
   4.191   (  -0.350    0.607    0.000)    0.701
   5.995   (   2.026   -3.510    0.000)    4.053
   6.462   (   2.598   -4.499    0.000)    5.195
  10.555   (   2.384   -4.129    0.000)    4.768
  10.749   (   6.807  -11.790    0.000)   13.614
  11.237   (   0.544   -0.942    0.000)    1.088
  11.376   (  -1.955    3.385    0.000)    3.909
  11.671   (   1.249   -2.164    0.000)    2.499
  12.026   (   1.663   -2.881    0.000)    3.326
======================= Grid point 33 (16/135) =======================
q-point: ( 0.13  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.082   (  13.514   23.408    0.000)   27.029
   2.302   (  17.100   29.619    0.000)   34.201
   2.845   (  13.573   23.508    0.000)   27.145
   3.804   (  20.456   35.431    0.000)   40.912
   3.947   (  11.718   20.296    0.000)   23.436
   6.240   (  -4.899   -8.485    0.000)    9.798
  10.555   (  -7.264  -12.582    0.000)   14.529
  11.468   (  -5.319   -9.213    0.000)   10.639
  11.792   (  -4.925   -8.530    0.000)    9.849
  11.816   (  -3.014   -5.221    0.000)    6.029
  11.939   (  -6.578  -11.393    0.000)   13.156
  12.142   (  -1.578   -2.732    0.000)    3.155
======================= Grid point 34 (17/135) =======================
q-point: ( 0.20  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.537   (  14.967   18.738    0.000)   23.982
   2.780   (   2.456   31.640    0.000)   31.735
   3.215   (  -1.551   28.218    0.000)   28.260
   4.247   (   9.379    3.245    0.000)    9.924
   4.488   (  26.669   20.640    0.000)   33.723
   6.144   (   0.801   -2.572    0.000)    2.693
  10.338   (  -5.392   -8.781    0.000)   10.304
  11.299   (  -4.318   -8.985    0.000)    9.969
  11.590   (  -9.378   -4.450    0.000)   10.380
  11.689   (  -6.994  -11.081    0.000)   13.104
  11.737   (  -1.318   -7.401    0.000)    7.518
  12.088   (  -0.576   -4.527    0.000)    4.564
======================= Grid point 35 (18/135) =======================
q-point: ( 0.27  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.953   (  14.003   14.800    0.000)   20.375
   3.007   ( -11.365   30.817    0.000)   32.846
   3.335   ( -17.666   30.422    0.000)   35.179
   4.339   (   1.721   -2.513    0.000)    3.046
   5.101   (  26.657   10.105    0.000)   28.508
   6.203   (   7.124   -1.738    0.000)    7.333
  10.210   (  -0.073   -6.899    0.000)    6.900
  11.113   (  -5.104  -13.904    0.000)   14.811
  11.402   (  -7.295   -3.287    0.000)    8.001
  11.505   (  -3.968   -5.538    0.000)    6.813
  11.658   (  -0.360   -5.751    0.000)    5.762
  12.035   (  -0.196   -4.726    0.000)    4.730
======================= Grid point 36 (19/135) =======================
q-point: ( 0.33  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.034   ( -19.099   29.924    0.000)   35.499
   3.235   ( -20.583   30.304    0.000)   36.633
   3.316   (   3.733   15.322    0.000)   15.770
   4.292   (  -3.530   -1.933    0.000)    4.025
   5.578   (  20.785    3.430    0.000)   21.066
   6.301   (   7.086   -5.035    0.000)    8.692
  10.209   (   7.658   -6.407    0.000)    9.985
  10.853   (  -5.602  -18.816    0.000)   19.632
  11.285   (  -3.758   -1.834    0.000)    4.182
  11.437   (  -1.439   -0.072    0.000)    1.441
  11.613   (   1.406   -5.255    0.000)    5.439
  11.993   (   0.296   -3.715    0.000)    3.726
======================= Grid point 37 (20/135) =======================
q-point: ( 0.40  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.992   ( -19.822   31.042    0.000)   36.831
   3.129   ( -27.881   40.721    0.000)   49.351
   3.520   (   2.706    7.493    0.000)    7.967
   4.217   (  -2.584   -0.711    0.000)    2.681
   5.858   (  10.406   -2.269    0.000)   10.650
   6.333   (   4.795   -7.980    0.000)    9.310
  10.341   (  14.522   -5.903    0.000)   15.676
  10.575   (  -4.181  -19.492    0.000)   19.935
  11.231   (  -0.972   -0.926    0.000)    1.342
  11.430   (  -0.992    1.904    0.000)    2.147
  11.596   (   2.804   -5.627    0.000)    6.287
  11.973   (   0.907   -2.531    0.000)    2.689
======================= Grid point 49 (21/135) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.902   (  10.494   18.177    0.000)   20.988
   3.305   (  11.575   20.048    0.000)   23.149
   3.699   (  10.722   18.571    0.000)   21.445
   4.270   (   0.204    0.353    0.000)    0.407
   4.882   (  11.489   19.899    0.000)   22.978
   6.123   (   0.152    0.263    0.000)    0.303
  10.175   (  -4.303   -7.454    0.000)    8.607
  11.086   (  -7.659  -13.266    0.000)   15.318
  11.481   (  -3.704   -6.415    0.000)    7.408
  11.495   (  -4.057   -7.027    0.000)    8.114
  11.613   (  -2.813   -4.873    0.000)    5.627
  12.013   (  -1.775   -3.074    0.000)    3.549
======================= Grid point 50 (22/135) =======================
q-point: ( 0.27  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.248   (  10.377   15.065    0.000)   18.293
   3.556   (  -6.291   23.229    0.000)   24.066
   3.905   ( -10.929   25.401    0.000)   27.652
   4.286   (   2.238   -2.689    0.000)    3.498
   5.277   (  17.635    8.627    0.000)   19.632
   6.145   (   3.992   -3.880    0.000)    5.567
  10.060   (   0.397   -7.852    0.000)    7.862
  10.785   (  -8.320  -18.218    0.000)   20.028
  11.347   (  -5.555   -2.323    0.000)    6.021
  11.430   (   0.625   -3.081    0.000)    3.144
  11.539   (  -0.193   -5.298    0.000)    5.302
  11.967   (  -0.552   -2.078    0.000)    2.150
======================= Grid point 51 (23/135) =======================
q-point: ( 0.33  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.536   (   6.503   12.638    0.000)   14.213
   3.600   ( -15.236   25.747    0.000)   29.917
   3.905   ( -18.539   29.080    0.000)   34.487
   4.258   (  -3.072   -1.176    0.000)    3.290
   5.621   (  15.332    1.667    0.000)   15.422
   6.161   (   4.675   -8.718    0.000)    9.892
  10.054   (   8.585   -8.440    0.000)   12.039
  10.472   (  -6.151  -18.850    0.000)   19.828
  11.255   (  -3.044   -1.094    0.000)    3.235
  11.434   (   2.087   -1.592    0.000)    2.624
  11.495   (   0.707   -4.979    0.000)    5.029
  11.950   (   0.047   -0.889    0.000)    0.890
======================= Grid point 52 (24/135) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.590   ( -15.918   27.570    0.000)   31.835
   3.669   (  -4.267    7.391    0.000)    8.535
   3.885   ( -18.654   32.309    0.000)   37.307
   4.210   (  -0.449    0.778    0.000)    0.899
   5.777   (   3.173   -5.495    0.000)    6.345
   6.147   (   5.893  -10.207    0.000)   11.786
  10.149   (   8.164  -14.140    0.000)   16.327
  10.262   (   5.572   -9.650    0.000)   11.143
  11.224   (  -0.141    0.244    0.000)    0.282
  11.458   (  -0.498    0.862    0.000)    0.996
  11.464   (   4.063   -7.037    0.000)    8.126
  11.947   (   0.237   -0.411    0.000)    0.475
======================= Grid point 66 (25/135) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.542   (   8.333   14.434    0.000)   16.667
   3.889   (   6.027   10.440    0.000)   12.055
   4.222   (  -1.546   -2.678    0.000)    3.093
   4.291   (   6.804   11.785    0.000)   13.608
   5.469   (   6.303   10.917    0.000)   12.606
   6.062   (  -2.506   -4.340    0.000)    5.011
   9.932   (  -2.418   -4.188    0.000)    4.836
  10.429   (  -9.918  -17.179    0.000)   19.837
  11.297   (  -1.737   -3.008    0.000)    3.474
  11.368   (  -1.425   -2.468    0.000)    2.850
  11.478   (  -0.685   -1.187    0.000)    1.371
  11.946   (  -0.192   -0.333    0.000)    0.385
======================= Grid point 67 (26/135) =======================
q-point: ( 0.33  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.787   (   2.902   12.361    0.000)   12.697
   4.002   (  -5.596   14.127    0.000)   15.195
   4.143   (  -2.487   -3.778    0.000)    4.523
   4.426   (  -6.003   15.475    0.000)   16.599
   5.674   (   6.690    4.341    0.000)    7.975
   5.975   (   0.939   -9.608    0.000)    9.654
   9.920   (   5.450   -3.895    0.000)    6.699
  10.138   (  -4.411  -14.258    0.000)   14.925
  11.245   (  -1.966    0.136    0.000)    1.971
  11.342   (   2.338   -5.127    0.000)    5.635
  11.469   (  -0.194    0.149    0.000)    0.244
  11.945   (   0.001    0.156    0.000)    0.156
======================= Grid point 82 (27/135) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.006   (  -4.675   -8.097    0.000)    9.350
   4.006   (   4.675    8.097    0.000)    9.350
   4.143   (   0.000    0.000    0.000)    0.000
   4.587   (   0.000    0.000    0.000)    0.000
   5.798   (  -3.970   -6.877    0.000)    7.940
   5.798   (   3.970    6.877    0.000)    7.940
   9.935   (  -2.865   -4.963    0.000)    5.731
   9.935   (   2.865    4.963    0.000)    5.731
  11.267   (  -1.170   -2.026    0.000)    2.339
  11.267   (   1.170    2.026    0.000)    2.339
  11.471   (  -0.000   -0.000    0.000)    0.000
  11.948   (  -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: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.583   (   0.000    0.000   28.892)   28.892
   0.583   (   0.000    0.000   28.892)   28.892
   1.354   (   0.000    0.000   68.696)   68.696
   1.668   (   0.000    0.000    4.737)    4.737
   1.668   (  -0.000    0.000    4.737)    4.737
   6.754   (   0.000    0.000  -13.574)   13.574
  11.508   (   0.000   -0.000    3.141)    3.141
  12.020   (   0.000    0.000    0.341)    0.341
  12.020   (   0.000    0.000    0.341)    0.341
  12.083   (   0.000    0.000   -1.827)    1.827
  12.317   (   0.000    0.000   -2.017)    2.017
  12.317   (   0.000    0.000   -2.017)    2.017
======================= Grid point 242 (29/135) =======================
q-point: ( 0.07  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.836   (  18.729   10.813   19.749)   29.287
   1.122   (  25.558   14.756   25.788)   39.191
   1.673   (  31.691   18.297   46.093)   58.853
   1.784   (   9.636    5.564    4.154)   11.877
   2.108   (  33.678   19.444    3.507)   39.045
   6.681   (  -6.284   -3.628  -13.440)   15.273
  11.379   ( -10.396   -6.002    3.446)   12.489
  11.971   (  -4.196   -2.422    0.318)    4.855
  12.012   (  -4.204   -2.427   -1.327)    5.032
  12.062   (   1.052    0.607   -0.625)    1.366
  12.272   (  -3.894   -2.248   -2.071)    4.950
  12.298   (  -1.611   -0.930   -1.895)    2.655
======================= Grid point 243 (30/135) =======================
q-point: ( 0.13  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.316   (  21.877   12.631   11.902)   27.925
   1.577   (  18.120   10.462   27.507)   34.560
   2.068   (  14.239    8.221    3.036)   16.720
   2.639   (  46.714   26.970   19.828)   57.470
   2.983   (  39.671   22.904    1.781)   45.843
   6.481   ( -10.665   -6.158  -12.878)   17.819
  11.080   ( -14.874   -8.588    4.112)   17.661
  11.837   (  -7.317   -4.225    0.249)    8.453
  11.858   (  -8.442   -4.874   -2.628)   10.096
  12.079   (  -0.679   -0.392   -0.290)    0.836
  12.146   (  -6.949   -4.012   -2.204)    8.321
  12.253   (  -1.968   -1.136   -1.526)    2.737
======================= Grid point 244 (31/135) =======================
q-point: ( 0.20  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.788   (  18.647   10.766    8.123)   23.012
   1.999   (  18.839   10.877   15.045)   26.450
   2.380   (  11.988    6.921    2.302)   14.032
   3.681   (  43.434   25.077   13.310)   51.889
   3.803   (  30.292   17.489    0.105)   34.978
   6.239   (  -9.035   -5.216  -11.080)   15.219
  10.745   ( -13.636   -7.873    4.962)   16.509
  11.650   (  -9.267   -5.350   -3.271)   11.189
  11.653   (  -8.424   -4.864    0.142)    9.728
  11.970   (  -8.035   -4.639   -2.372)    9.576
  12.007   (  -5.614   -3.241    0.217)    6.487
  12.219   (  -1.230   -0.710   -1.383)    1.982
======================= Grid point 245 (32/135) =======================
q-point: ( 0.27  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.135   (  11.300    6.524    6.256)   14.471
   2.426   (  18.288   10.559    3.957)   21.485
   2.567   (   3.586    2.070    2.454)    4.813
   4.257   (   9.441    5.451   -5.281)   12.113
   4.609   (  36.583   21.121   12.195)   43.968
   6.150   (   2.339    1.351   -6.188)    6.752
  10.492   (  -8.030   -4.636    5.957)   11.021
  11.455   (  -7.557   -4.363   -3.208)    9.297
  11.471   (  -7.129   -4.116    0.025)    8.231
  11.800   (  -6.429   -3.712   -2.565)    7.854
  11.837   (  -9.032   -5.214    1.331)   10.513
  12.184   (  -1.893   -1.093   -1.473)    2.636
======================= Grid point 246 (33/135) =======================
q-point: ( 0.33  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.293   (   2.726    1.574    5.377)    6.231
   2.525   (  -7.228   -4.173    3.751)    9.151
   2.817   (  15.719    9.075   -4.219)   18.635
   4.279   (  -4.798   -2.770  -10.665)   12.019
   5.341   (  27.435   15.840   11.808)   33.808
   6.309   (   9.376    5.413   -0.043)   10.827
  10.392   (  -0.709   -0.409    6.772)    6.821
  11.312   (  -5.079   -2.932   -2.188)    6.260
  11.340   (  -4.265   -2.463   -0.069)    4.926
  11.609   ( -10.485   -6.054    2.558)   12.375
  11.691   (  -3.082   -1.779   -2.802)    4.530
  12.134   (  -2.424   -1.400   -1.580)    3.215
======================= Grid point 247 (34/135) =======================
q-point: ( 0.40  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.270   ( -11.445   -6.608    5.895)   14.471
   2.288   (  -4.489   -2.591    5.318)    7.426
   3.126   (  10.959    6.327   -9.771)   15.988
   4.141   (  -5.925   -3.421  -13.702)   15.315
   5.853   (  17.393   10.042   12.045)   23.418
   6.494   (   5.989    3.458    4.214)    8.098
  10.450   (   5.423    3.131    6.802)    9.246
  11.191   (  -6.919   -3.995    1.342)    8.102
  11.275   (  -1.623   -0.937   -0.127)    1.878
  11.399   (  -6.446   -3.721    1.842)    7.667
  11.654   (  -0.395   -0.228   -3.057)    3.090
  12.081   (  -2.030   -1.172   -1.659)    2.871
======================= Grid point 248 (35/135) =======================
q-point: ( 0.47  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.011   (  -7.257   -4.190    8.620)   12.022
   2.222   (  -1.952   -1.127    5.007)    5.491
   3.301   (   4.065    2.347  -12.810)   13.643
   4.042   (  -2.351   -1.357  -14.891)   15.136
   6.127   (   6.323    3.651   12.325)   14.325
   6.574   (   1.354    0.782    6.567)    6.750
  10.603   (   6.331    3.655    5.516)    9.158
  11.000   (  -7.561   -4.365    6.113)   10.658
  11.256   (  -0.301   -0.174   -0.153)    0.380
  11.328   (  -0.920   -0.531    0.001)    1.062
  11.656   (   0.225    0.130   -3.228)    3.238
  12.048   (  -0.792   -0.457   -1.701)    1.931
======================= Grid point 257 (36/135) =======================
q-point: ( 0.07  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.237   (  15.042   26.054   12.888)   32.728
   1.453   (  10.163   17.603   29.209)   35.585
   2.019   (   9.979   17.285    3.326)   20.234
   2.337   (  25.262   43.755   24.315)   56.070
   2.714   (  21.957   38.031    2.109)   43.965
   6.544   (  -5.698   -9.869  -13.101)   17.364
  11.170   (  -8.189  -14.183    3.922)   16.841
  11.841   (  -5.160   -8.937   -0.793)   10.350
  11.982   (  -2.204   -3.817   -1.567)    4.678
  12.037   (  -1.265   -2.191   -0.229)    2.540
  12.216   (  -1.874   -3.246   -2.026)    4.261
  12.242   (  -2.533   -4.387   -1.596)    5.311
======================= Grid point 258 (37/135) =======================
q-point: ( 0.13  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.755   (  11.873   28.380    8.924)   32.032
   1.836   (  15.390   16.262   18.474)   29.028
   2.387   (   8.407   22.362    2.248)   23.995
   3.264   (  35.590   34.349   16.029)   51.995
   3.482   (  28.134   25.516    0.516)   37.985
   6.318   (  -8.163   -9.632  -12.020)   17.433
  10.862   ( -11.158  -12.652    4.687)   17.508
  11.648   (  -5.695  -10.372   -1.456)   11.922
  11.840   ( -10.286   -1.290   -1.432)   10.465
  11.974   (  -0.687   -8.263   -1.071)    8.360
  12.091   (  -7.709   -4.057   -1.191)    8.793
  12.191   (   0.972   -4.272   -1.559)    4.650
======================= Grid point 259 (38/135) =======================
q-point: ( 0.20  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.174   (   7.973   22.849    4.943)   24.700
   2.271   (  12.942   19.826    8.931)   25.304
   2.695   (  -0.157   23.590    1.805)   23.659
   4.077   (  19.195   10.469   -2.684)   22.029
   4.186   (  36.541   23.852   13.035)   45.542
   6.139   (  -2.549   -4.212   -8.812)   10.094
  10.575   (  -9.072   -8.803    5.702)   13.867
  11.459   (  -4.746   -9.164   -1.669)   10.454
  11.635   ( -11.197   -1.808   -1.483)   11.438
  11.853   (  -4.355   -6.755   -2.094)    8.306
  11.900   (  -7.270   -8.008    0.420)   10.824
  12.164   (   0.528   -4.414   -1.600)    4.725
======================= Grid point 260 (39/135) =======================
q-point: ( 0.27  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.452   (  -2.856   22.911    3.010)   23.284
   2.661   (  13.200   14.142    0.914)   19.367
   2.815   ( -12.299   21.780    2.163)   25.106
   4.288   (   2.010   -1.621   -8.413)    8.801
   4.979   (  32.266   15.986   12.226)   38.028
   6.186   (   8.328    1.934   -2.800)    8.996
  10.399   (  -3.320   -4.833    6.772)    8.958
  11.311   (  -3.017   -7.492   -1.132)    8.156
  11.442   (  -8.630   -2.074   -1.450)    8.993
  11.677   (  -7.067   -9.295    0.731)   11.699
  11.733   (  -4.140   -4.208   -1.564)    6.107
  12.120   (  -0.115   -5.010   -1.694)    5.290
======================= Grid point 261 (40/135) =======================
q-point: ( 0.33  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.532   ( -12.242   20.974    3.571)   24.546
   2.702   ( -23.864   21.215    3.738)   32.149
   3.006   (  11.759   10.139   -6.708)   16.913
   4.209   (  -5.510   -3.719  -12.636)   14.278
   5.590   (  24.285    9.291   12.140)   28.696
   6.378   (   9.919    1.273    2.443)   10.295
  10.370   (   3.913   -2.152    7.326)    8.580
  11.183   (  -2.835  -10.486    1.025)   10.911
  11.310   (  -4.414   -1.941   -1.321)    5.000
  11.483   (  -6.201   -4.704    1.766)    7.981
  11.657   (  -1.200   -2.120   -2.667)    3.612
  12.067   (  -0.171   -4.877   -1.773)    5.193
======================= Grid point 262 (41/135) =======================
q-point: ( 0.40  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.472   ( -26.668   25.852    5.517)   37.549
   2.491   ( -16.561   21.058    4.042)   27.093
   3.252   (   6.676    6.574  -11.273)   14.659
   4.086   (  -4.564   -2.113  -14.678)   15.516
   5.985   (  14.332    3.445   12.367)   19.241
   6.508   (   5.507   -2.115    5.825)    8.290
  10.477   (   9.476   -0.795    6.764)   11.670
  10.988   (  -4.117  -15.467    4.549)   16.640
  11.247   (  -1.204   -1.459   -0.876)    2.085
  11.394   (  -3.240    1.877    0.777)    3.824
  11.639   (   0.852   -1.754   -3.052)    3.621
  12.025   (   0.524   -4.037   -1.829)    4.462
======================= Grid point 263 (42/135) =======================
q-point: ( 0.47  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.352   ( -19.916   34.495    6.652)   40.383
   2.453   ( -13.060   22.621    4.042)   26.431
   3.345   (  -1.273    2.205  -12.931)   13.179
   4.033   (  -0.211    0.365  -15.206)   15.212
   6.130   (   1.938   -3.357   12.480)   13.068
   6.539   (   2.450   -4.244    7.091)    8.619
  10.605   (   2.288   -3.962    4.991)    6.771
  10.823   (   6.741  -11.675    7.261)   15.312
  11.231   (   0.583   -1.009   -0.657)    1.338
  11.380   (  -2.038    3.530    0.464)    4.103
  11.639   (   1.176   -2.037   -3.171)    3.948
  12.009   (   1.765   -3.057   -1.852)    3.987
======================= Grid point 274 (43/135) =======================
q-point: ( 0.13  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.158   (  10.329   17.891    8.355)   22.284
   2.363   (  16.544   28.656    5.941)   33.618
   2.859   (  13.408   23.224    1.417)   26.854
   3.930   (  10.745   18.611   -2.013)   21.585
   3.954   (  20.188   34.967   14.437)   42.880
   6.142   (  -3.962   -6.863   -9.747)   12.562
  10.612   (  -6.899  -11.949    5.687)   14.924
  11.452   (  -5.277   -9.140   -1.593)   10.674
  11.781   (  -4.981   -8.628   -1.077)   10.020
  11.798   (  -3.051   -5.285   -1.762)    6.351
  11.936   (  -6.342  -10.985   -0.272)   12.687
  12.123   (  -1.634   -2.830   -1.924)    3.792
======================= Grid point 275 (44/135) =======================
q-point: ( 0.20  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.512   (  11.501   15.344   -1.084)   19.206
   2.834   (   3.614   30.254    5.153)   30.902
   3.228   (  -1.400   27.970    1.210)   28.031
   4.185   (   7.478    1.104   -6.539)    9.995
   4.631   (  26.054   21.118   13.540)   36.168
   6.092   (   2.869   -0.392   -5.054)    5.824
  10.408   (  -5.038   -7.985    7.017)   11.763
  11.290   (  -3.957   -8.321   -0.932)    9.260
  11.574   (  -9.611   -4.441   -1.515)   10.695
  11.688   (  -5.416  -10.827   -0.863)   12.137
  11.725   (  -2.608   -7.620   -0.404)    8.064
  12.068   (  -0.496   -4.762   -2.074)    5.218
======================= Grid point 276 (45/135) =======================
q-point: ( 0.27  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.836   (   8.127   14.418   -8.179)   18.462
   3.084   (  -6.753   27.794    5.243)   29.080
   3.349   ( -17.193   29.823    1.306)   34.448
   4.226   (   0.009   -3.968  -11.237)   11.917
   5.240   (  26.128   10.841   13.057)   31.155
   6.209   (   9.298    0.075    0.533)    9.314
  10.293   (  -0.067   -6.168    8.117)   10.195
  11.125   (  -4.305  -12.555    1.180)   13.325
  11.385   (  -7.240   -3.290   -1.786)    8.151
  11.510   (  -3.627   -6.268    0.379)    7.252
  11.639   (  -0.806   -5.442   -1.689)    5.755
  12.014   (  -0.073   -5.015   -2.175)    5.467
======================= Grid point 277 (46/135) =======================
q-point: ( 0.33  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.008   ( -10.648   22.996   -4.828)   25.798
   3.228   (  -5.420   20.391   -1.138)   21.130
   3.284   ( -19.459   28.704    0.350)   34.680
   4.145   (  -4.522   -2.639  -14.080)   15.022
   5.717   (  20.481    4.031   12.912)   24.545
   6.354   (   8.534   -3.620    4.858)   10.466
  10.292   (   7.111   -5.963    8.232)   12.405
  10.890   (  -4.766  -17.948    3.704)   18.936
  11.272   (  -3.445   -1.620   -1.310)    4.026
  11.440   (  -1.235   -0.530    0.347)    1.388
  11.590   (   0.951   -4.866   -2.266)    5.451
  11.971   (   0.443   -4.052   -2.267)    4.664
======================= Grid point 278 (47/135) =======================
q-point: ( 0.40  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.000   ( -17.461   28.109   -0.230)   33.092
   3.159   ( -26.529   38.958    2.645)   47.207
   3.399   (   0.248    9.116   -9.736)   13.340
   4.055   (  -2.801   -0.982  -15.297)   15.582
   5.998   (  10.202   -1.894   12.847)   16.514
   6.415   (   5.149   -7.062    7.392)   11.447
  10.413   (  13.284   -5.584    7.088)   16.059
  10.634   (  -3.055  -19.111    5.932)   20.242
  11.225   (  -0.869   -0.786   -0.716)    1.373
  11.433   (  -0.898    1.803    0.432)    2.060
  11.570   (   2.495   -5.264   -2.656)    6.403
  11.950   (   1.087   -2.889   -2.342)    3.875
======================= Grid point 290 (48/135) =======================
q-point: ( 0.20  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.814   (   8.758   15.169   -7.468)   19.041
   3.342   (  11.382   19.715    3.693)   23.062
   3.709   (  10.688   18.512    0.990)   21.399
   4.163   (  -1.122   -1.944  -10.801)   11.032
   5.035   (  11.697   20.259   14.177)   27.354
   6.118   (   1.603    2.776   -0.529)    3.249
  10.261   (  -3.918   -6.787    8.519)   11.575
  11.098   (  -6.925  -11.995    1.163)   13.899
  11.467   (  -3.715   -6.435   -1.446)    7.570
  11.486   (  -4.313   -7.470   -1.277)    8.720
  11.605   (  -2.707   -4.689   -0.543)    5.442
  11.988   (  -1.862   -3.225   -2.453)    4.459
======================= Grid point 291 (49/135) =======================
q-point: ( 0.27  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.104   (   8.284   13.462  -12.491)   20.146
   3.598   (  -4.726   22.001    3.637)   22.795
   3.900   ( -11.745   23.534   -1.364)   26.337
   4.153   (   2.146   -2.591  -11.917)   12.383
   5.432   (  17.098    9.372   14.229)   24.138
   6.189   (   5.544   -1.854    3.953)    7.057
  10.154   (   0.266   -7.470    9.282)   11.917
  10.819   (  -7.722  -17.353    3.365)   19.290
  11.332   (  -5.420   -2.030   -1.512)    5.982
  11.419   (   1.127   -3.853   -1.254)    4.206
  11.533   (  -0.656   -4.489   -0.485)    4.563
  11.940   (  -0.537   -2.321   -2.706)    3.606
======================= Grid point 292 (50/135) =======================
q-point: ( 0.33  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.340   (   2.853   13.188  -15.627)   20.646
   3.668   ( -11.363   23.038    4.055)   26.005
   3.874   ( -15.438   21.965   -5.316)   27.369
   4.128   (  -6.158    4.516   -9.774)   12.403
   5.772   (  14.723    2.222   13.765)   20.277
   6.245   (   5.517   -7.098    7.474)   11.691
  10.144   (   7.972   -8.336    8.814)   14.516
  10.523   (  -5.564  -18.240    5.070)   19.732
  11.248   (  -2.818   -0.876   -0.836)    3.068
  11.424   (   2.784   -2.743   -1.183)    4.083
  11.488   (  -0.123   -3.507   -0.465)    3.540
  11.921   (   0.129   -1.164   -2.883)    3.112
======================= Grid point 293 (51/135) =======================
q-point: ( 0.40  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.431   (  -7.221   12.507  -15.462)   21.157
   3.691   ( -12.298   21.300    3.318)   24.818
   3.843   ( -12.793   22.157   -7.225)   26.586
   4.083   (  -5.396    9.347   -8.642)   13.826
   5.925   (   2.907   -5.034   13.394)   14.602
   6.247   (   5.343   -9.254    8.963)   13.946
  10.232   (   8.060  -13.960    8.232)   18.100
  10.322   (   5.402   -9.357    5.932)   12.326
  11.219   (  -0.174    0.302   -0.503)    0.612
  11.445   (   3.812   -6.602   -1.902)    7.857
  11.460   (  -0.459    0.796    0.303)    0.967
  11.917   (   0.383   -0.664   -2.952)    3.050
======================= Grid point 307 (52/135) =======================
q-point: ( 0.27  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.369   (   7.495   12.983  -15.585)   21.624
   3.907   (   5.112    8.855   -0.340)   10.230
   4.056   (  -1.350   -2.339  -13.443)   13.712
   4.298   (   6.755   11.701    0.650)   13.527
   5.632   (   6.303   10.917   14.694)   19.360
   6.146   (  -1.362   -2.360    7.511)    7.990
  10.027   (  -2.554   -4.424    9.413)   10.710
  10.478   (  -9.488  -16.433    4.875)   19.592
  11.289   (  -1.604   -2.778   -0.930)    3.340
  11.353   (  -1.282   -2.220   -1.643)    3.045
  11.480   (  -0.658   -1.139    0.268)    1.342
  11.915   (  -0.275   -0.476   -3.084)    3.132
======================= Grid point 308 (53/135) =======================
q-point: ( 0.33  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.587   (   1.970   11.955  -17.698)   21.448
   3.936   (  -1.426   -3.416  -17.598)   17.983
   4.053   (  -5.509   11.919    3.064)   13.483
   4.432   (  -5.860   15.506    0.538)   16.586
   5.831   (   6.221    4.256   13.982)   15.884
   6.092   (   1.222   -7.944   10.302)   13.066
  10.007   (   5.086   -4.258    8.546)   10.818
  10.202   (  -4.095  -13.612    6.258)   15.531
  11.240   (  -1.851    0.210   -0.529)    1.937
  11.331   (   2.266   -4.747   -1.165)    5.388
  11.470   (  -0.203    0.101    0.210)    0.309
  11.913   (   0.020    0.029   -3.270)    3.270
======================= Grid point 323 (54/135) =======================
q-point: ( 0.33  0.33  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.800   (  -4.610   -7.985  -18.698)   20.847
   3.800   (   4.610    7.985  -18.698)   20.847
   4.179   (   0.000    0.000    3.610)    3.610
   4.590   (   0.000    0.000    0.233)    0.233
   5.943   (  -3.424   -5.930   12.766)   14.486
   5.943   (   3.424    5.930   12.766)   14.486
  10.012   (  -2.545   -4.408    7.526)    9.086
  10.012   (   2.545    4.408    7.526)    9.086
  11.262   (  -1.092   -1.891   -0.604)    2.265
  11.262   (   1.092    1.891   -0.604)    2.265
  11.472   (  -0.000   -0.000    0.156)    0.156
  11.914   (  -0.000   -0.000   -3.391)    3.391
======================= Grid point 482 (55/135) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.102   (   0.000    0.000   24.108)   24.108
   1.102   (   0.000    0.000   24.108)   24.108
   1.770   (  -0.000    0.000    4.861)    4.861
   1.770   (   0.000    0.000    4.861)    4.861
   2.651   (   0.000   -0.000   64.456)   64.456
   6.363   (  -0.000   -0.000  -26.499)   26.499
  11.593   (   0.000   -0.000    5.406)    5.406
  12.030   (   0.000    0.000   -3.593)    3.593
  12.033   (   0.000   -0.000    1.055)    1.055
  12.033   (   0.000    0.000    1.055)    1.055
  12.262   (  -0.000    0.000   -3.444)    3.444
  12.262   (  -0.000    0.000   -3.444)    3.444
======================= Grid point 483 (56/135) =======================
q-point: ( 0.07  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.249   (  12.021    6.940   20.945)   25.127
   1.546   (  30.159   17.412   19.842)   40.080
   1.873   (   8.631    4.983    4.207)   10.818
   2.179   (  31.530   18.204    3.348)   36.562
   2.748   (   9.851    5.688   58.374)   59.472
   6.295   (  -5.813   -3.356  -26.084)   26.934
  11.473   (  -9.737   -5.622    5.967)   12.728
  11.962   (  -5.827   -3.364   -3.846)    7.750
  11.984   (  -4.242   -2.449    1.029)    5.005
  12.058   (   1.570    0.906    0.364)    1.849
  12.216   (  -4.009   -2.315   -3.527)    5.820
  12.247   (  -1.255   -0.724   -3.156)    3.472
======================= Grid point 484 (57/135) =======================
q-point: ( 0.13  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.597   (  17.321   10.000   15.658)   25.401
   2.127   (  17.275    9.974   26.263)   32.980
   2.133   (  13.211    7.627    2.923)   15.532
   3.003   (  37.367   21.574    1.501)   43.173
   3.209   (  31.401   18.130   35.846)   50.987
   6.115   (  -9.203   -5.313  -24.451)   26.661
  11.191   ( -14.006   -8.086    7.093)   17.660
  11.782   (  -9.397   -5.425   -5.099)   11.989
  11.848   (  -7.411   -4.279    0.946)    8.610
  12.075   (  -0.590   -0.340   -0.015)    0.681
  12.086   (  -7.124   -4.113   -3.733)    9.034
  12.212   (  -1.605   -0.927   -2.566)    3.166
======================= Grid point 485 (58/135) =======================
q-point: ( 0.20  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.989   (  16.079    9.283   11.736)   21.965
   2.395   (   8.696    5.020   24.381)   26.368
   2.429   (  11.687    6.748    2.128)   13.662
   3.758   (  26.015   15.020   -3.576)   30.252
   4.049   (  39.111   22.581   21.792)   50.144
   5.935   (  -4.651   -2.685  -19.702)   20.420
  10.877   ( -12.627   -7.290    8.364)   16.809
  11.559   (  -9.637   -5.564   -6.063)   12.672
  11.661   (  -8.545   -4.934    0.811)    9.900
  11.906   (  -8.235   -4.754   -3.990)   10.312
  12.015   (  -4.631   -2.673    0.702)    5.393
  12.181   (  -1.314   -0.759   -2.439)    2.872
======================= Grid point 486 (59/135) =======================
q-point: ( 0.27  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.293   (   9.925    5.730    9.322)   14.773
   2.582   (   8.130    4.694   12.609)   15.720
   2.621   (   4.405    2.543    2.554)    5.691
   4.088   (   2.747    1.586  -12.391)   12.790
   4.908   (  35.183   20.313   17.397)   44.194
   5.991   (   9.957    5.749   -9.610)   14.984
  10.648   (  -7.017   -4.051    9.670)   12.616
  11.365   (  -6.981   -4.030   -5.927)   10.005
  11.477   (  -7.236   -4.178    0.673)    8.382
  11.732   (  -6.675   -3.854   -4.294)    8.823
  11.874   (  -7.557   -4.363    2.467)    9.068
  12.142   (  -2.107   -1.216   -2.788)    3.700
======================= Grid point 487 (60/135) =======================
q-point: ( 0.33  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.428   (   2.049    1.183    7.927)    8.272
   2.612   (  -5.130   -2.962    4.566)    7.479
   2.767   (   7.969    4.601    0.563)    9.219
   3.987   (  -8.572   -4.949  -18.690)   21.149
   5.623   (  27.116   15.655   15.833)   35.086
   6.301   (  14.649    8.458   -0.984)   16.944
  10.566   (  -0.221   -0.128   10.607)   10.610
  11.248   (  -3.395   -1.960   -4.487)    5.958
  11.344   (  -4.329   -2.499    0.587)    5.033
  11.616   (  -3.386   -1.955   -4.687)    6.103
  11.681   (  -9.130   -5.271    4.678)   11.534
  12.087   (  -2.637   -1.522   -3.247)    4.451
======================= Grid point 488 (61/135) =======================
q-point: ( 0.40  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.406   (  -4.596   -2.653    7.405)    9.111
   2.423   (  -8.685   -5.014    7.743)   12.670
   2.931   (   6.066    3.502   -8.266)   10.834
   3.792   (  -7.469   -4.313  -20.998)   22.701
   6.133   (  17.404   10.048   15.229)   25.214
   6.583   (   9.282    5.359    4.130)   11.486
  10.624   (   4.734    2.733   10.392)   11.742
  11.195   (  -1.848   -1.067   -1.429)    2.568
  11.278   (  -1.647   -0.951    0.567)    1.984
  11.476   (  -8.140   -4.700    6.217)   11.270
  11.573   (  -0.670   -0.387   -5.120)    5.178
  12.030   (  -2.196   -1.268   -3.690)    4.477
======================= Grid point 489 (62/135) =======================
q-point: ( 0.47  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.215   (  -5.591   -3.228   11.169)   12.900
   2.344   (  -2.015   -1.163    6.846)    7.231
   3.030   (   2.299    1.327  -12.903)   13.174
   3.673   (  -2.694   -1.555  -21.605)   21.828
   6.407   (   6.314    3.645   14.869)   16.560
   6.716   (   2.625    1.515    7.019)    7.646
  10.744   (   4.346    2.509    8.468)    9.843
  11.134   (  -3.108   -1.795    5.556)    6.615
  11.258   (  -0.306   -0.177    0.576)    0.676
  11.351   (  -2.402   -1.387    3.770)    4.680
  11.570   (   0.112    0.065   -5.415)    5.417
  11.994   (  -0.860   -0.496   -3.986)    4.107
======================= Grid point 498 (63/135) =======================
q-point: ( 0.07  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.537   (  11.724   20.307   16.519)   28.683
   1.979   (  11.636   20.154   22.662)   32.483
   2.093   (   9.542   16.527    4.130)   19.526
   2.754   (  20.837   36.090    1.418)   41.697
   3.024   (  13.967   24.192   43.789)   51.940
   6.170   (  -5.082   -8.802  -25.082)   27.063
  11.276   (  -7.697  -13.332    6.778)   16.821
  11.808   (  -5.716   -9.900   -2.999)   11.818
  11.951   (  -2.225   -3.853   -1.107)    4.585
  12.032   (  -1.830   -3.169   -0.203)    3.665
  12.163   (  -1.621   -2.807   -3.206)    4.559
  12.199   (  -2.179   -3.774   -2.630)    5.090
======================= Grid point 499 (64/135) =======================
q-point: ( 0.13  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.969   (  10.232   25.721   12.298)   30.291
   2.287   (   9.738    7.517   25.183)   28.027
   2.438   (   7.786   20.791    3.039)   22.408
   3.472   (  25.554   23.393   -1.510)   34.677
   3.690   (  29.749   29.143   26.194)   49.198
   5.982   (  -5.946   -7.260  -22.092)   24.003
  10.987   ( -10.418  -11.668    7.962)   17.552
  11.598   (  -6.371  -10.514   -3.900)   12.898
  11.811   (  -9.966   -1.798   -1.146)   10.191
  11.943   (  -2.283   -8.396   -2.184)    8.971
  12.065   (  -6.088   -3.167   -1.297)    6.984
  12.151   (   1.139   -4.658   -2.471)    5.394
======================= Grid point 500 (65/135) =======================
q-point: ( 0.20  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.332   (   3.896   20.989   10.323)   23.712
   2.515   (   8.006   10.160   16.315)   20.820
   2.734   (   0.334   22.861    1.776)   22.933
   3.975   (  13.553    5.805   -8.310)   16.925
   4.513   (  34.108   23.366   19.258)   45.609
   5.906   (   3.401    0.486  -14.484)   14.886
  10.726   (  -8.253   -7.523    9.382)   14.586
  11.403   (  -5.170   -8.220   -4.303)   10.621
  11.606   ( -10.760   -2.463   -1.043)   11.087
  11.794   (  -5.759   -6.600   -3.817)    9.554
  11.916   (  -5.062   -7.441    1.271)    9.089
  12.120   (   0.575   -4.929   -2.778)    5.687
======================= Grid point 501 (66/135) =======================
q-point: ( 0.27  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.530   (   0.421   13.199    5.627)   14.355
   2.740   (   2.548   15.073    7.187)   16.892
   2.866   ( -10.660   21.128    2.652)   23.813
   4.045   (  -3.043   -5.032  -16.332)   17.359
   5.273   (  31.082   16.173   16.690)   38.810
   6.114   (  14.625    6.098   -4.511)   16.475
  10.575   (  -2.879   -3.406   10.768)   11.655
  11.269   (  -2.857   -4.947   -3.561)    6.732
  11.415   (  -7.884   -3.231   -0.846)    8.562
  11.642   (  -4.353   -5.722   -4.033)    8.243
  11.747   (  -5.754   -7.596    2.884)    9.956
  12.072   (  -0.099   -5.505   -3.209)    6.373
======================= Grid point 502 (67/135) =======================
q-point: ( 0.33  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.607   (  -7.413   14.165    3.435)   16.352
   2.787   ( -20.138   19.589    4.279)   28.418
   2.909   (   1.446   11.374   -1.305)   11.540
   3.881   (  -7.911   -5.213  -20.073)   22.197
   5.874   (  23.840    9.671   15.668)   30.122
   6.426   (  13.960    3.945    1.975)   14.641
  10.559   (   3.316   -1.099   11.364)   11.889
  11.196   (  -1.079   -4.494   -0.575)    4.657
  11.280   (  -3.569   -4.226   -1.052)    5.630
  11.525   (  -3.299   -6.208    0.710)    7.066
  11.601   (  -3.922   -2.518   -1.008)    4.768
  12.015   (  -0.151   -5.317   -3.642)    6.446
======================= Grid point 503 (68/135) =======================
q-point: ( 0.40  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.576   ( -12.848   16.897    4.378)   21.673
   2.607   ( -23.384   23.205    7.049)   33.690
   3.029   (   1.326    7.144   -9.258)   11.769
   3.721   (  -5.383   -2.640  -21.288)   22.116
   6.267   (  14.104    3.670   15.118)   20.999
   6.632   (   7.740   -0.395    6.001)    9.802
  10.649   (   7.357   -0.600   10.282)   12.657
  11.099   (  -0.970  -11.652    6.194)   13.232
  11.219   (  -0.260   -2.016   -2.167)    2.971
  11.430   (  -3.480    0.338    3.231)    4.761
  11.558   (   0.117   -1.613   -4.906)    5.166
  11.969   (   0.590   -4.411   -3.990)    5.977
======================= Grid point 504 (69/135) =======================
q-point: ( 0.47  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.509   ( -17.206   29.802    8.516)   35.450
   2.545   ( -11.131   19.279    4.679)   22.748
   3.078   (  -2.464    4.268  -12.049)   13.018
   3.660   (  -0.189    0.328  -21.498)   21.501
   6.411   (   1.850   -3.205   14.831)   15.286
   6.692   (   2.225   -3.854    7.589)    8.798
  10.736   (   2.018   -3.495    8.057)    9.011
  11.003   (   6.289  -10.892   10.263)   16.233
  11.209   (   0.802   -1.388   -1.769)    2.387
  11.406   (  -2.062    3.572    2.743)    4.953
  11.553   (   1.093   -1.893   -5.496)    5.915
  11.952   (   1.939   -3.359   -4.131)    5.667
======================= Grid point 515 (70/135) =======================
q-point: ( 0.13  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.419   (   3.996    6.921   18.217)   19.893
   2.514   (  15.267   26.443    9.031)   31.841
   2.890   (  13.110   22.708    1.356)   26.256
   3.850   (   7.813   13.532   -6.706)   17.004
   4.310   (  18.866   32.677   20.928)   43.147
   5.881   (  -0.827   -1.432  -16.468)   16.551
  10.763   (  -6.027  -10.440    9.384)   15.277
  11.403   (  -5.060   -8.765   -3.569)   10.732
  11.753   (  -5.123   -8.873   -1.715)   10.388
  11.753   (  -3.448   -5.972   -2.571)    7.360
  11.932   (  -5.534   -9.585   -0.063)   11.068
  12.073   (  -1.835   -3.179   -3.051)    4.773
======================= Grid point 516 (71/135) =======================
q-point: ( 0.20  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.556   (   4.601    6.688    6.659)   10.499
   2.962   (   4.378   27.715    7.498)   29.044
   3.254   (  -1.186   27.448    1.131)   27.497
   3.986   (   3.166   -3.866  -13.957)   14.825
   4.953   (  24.527   21.250   18.243)   37.229
   5.964   (   8.172    5.245   -7.708)   12.397
  10.592   (  -4.220   -6.028   11.265)   13.455
  11.260   (  -3.100   -6.320   -2.313)    7.410
  11.536   (  -9.904   -4.903   -2.220)   11.272
  11.638   (  -3.151   -9.569   -3.744)   10.747
  11.744   (  -4.313   -8.744    1.927)    9.938
  12.012   (  -0.419   -5.390   -3.528)    6.455
======================= Grid point 517 (72/135) =======================
q-point: ( 0.27  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.707   (   4.760    7.213   -3.454)    9.307
   3.199   (  -7.625   26.683    6.153)   28.425
   3.372   ( -15.814   27.356    0.593)   31.604
   3.932   (  -2.995   -5.074  -17.730)   18.684
   5.545   (  24.912   11.611   16.752)   32.187
   6.214   (  13.806    3.820   -0.317)   14.329
  10.503   (  -0.158   -4.290   12.811)   13.511
  11.153   (  -2.193   -8.406    1.369)    8.794
  11.335   (  -7.318   -4.216   -3.215)    9.037
  11.498   (  -1.750   -7.632   -2.460)    8.208
  11.619   (  -2.651   -5.231    0.710)    5.907
  11.953   (   0.085   -5.700   -3.928)    6.923
======================= Grid point 518 (73/135) =======================
q-point: ( 0.33  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.843   (   0.934    9.307   -9.048)   13.013
   3.256   ( -14.773   24.708    2.079)   28.863
   3.314   ( -19.998   27.318    1.981)   33.913
   3.810   (  -7.396   -0.471  -18.093)   19.552
   6.011   (  19.670    4.596   15.793)   25.641
   6.455   (  11.519   -0.971    4.639)   12.456
  10.504   (   5.568   -4.813   12.794)   14.760
  10.987   (  -2.704  -15.559    5.866)   16.847
  11.231   (  -2.525   -0.990   -3.107)    4.125
  11.441   (   0.700   -2.717   -1.255)    3.074
  11.543   (  -1.377   -3.421   -1.289)    3.906
  11.907   (   0.656   -4.764   -4.263)    6.426
======================= Grid point 519 (74/135) =======================
q-point: ( 0.40  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.909   (  -5.792   12.826   -9.004)   16.708
   3.201   ( -21.276   31.607    0.171)   38.101
   3.295   ( -12.125   21.746    0.684)   24.908
   3.703   (  -5.153    2.266  -18.324)   19.170
   6.285   (   9.748   -1.663   15.039)   17.999
   6.572   (   6.083   -5.453    7.652)   11.193
  10.592   (   9.714   -5.094   10.563)   15.228
  10.791   (   0.149  -17.729    9.544)   20.135
  11.201   (  -0.531   -0.352   -1.821)    1.929
  11.449   (   0.671    0.838    0.657)    1.259
  11.499   (   0.562   -4.103   -3.817)    5.632
  11.883   (   1.400   -3.584   -4.474)    5.901
======================= Grid point 531 (75/135) =======================
q-point: ( 0.20  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.695   (   4.479    7.759   -2.848)    9.401
   3.432   (  10.458   18.114    4.804)   21.461
   3.730   (  10.593   18.347    0.952)   21.207
   3.874   (  -3.131   -5.423  -17.821)   18.889
   5.361   (  11.691   20.249   17.755)   29.359
   6.100   (   4.677    8.101   -1.570)    9.485
  10.483   (  -2.951   -5.112   13.560)   14.789
  11.127   (  -4.825   -8.357    1.649)    9.790
  11.426   (  -3.865   -6.695   -2.713)    8.192
  11.431   (  -4.774   -8.269   -4.528)   10.568
  11.610   (  -2.841   -4.921    1.391)    5.850
  11.922   (  -2.105   -3.645   -4.187)    5.937
======================= Grid point 532 (76/135) =======================
q-point: ( 0.27  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.861   (   5.109    8.611  -10.388)   14.428
   3.624   (  -4.607    8.988   -9.214)   13.671
   3.779   (  -7.590   13.070   -4.212)   15.690
   3.999   (  -2.202   16.626   -1.844)   16.873
   5.751   (  16.035    9.909   16.930)   25.337
   6.268   (   8.530    1.932    3.474)    9.411
  10.395   (  -0.147   -6.548   14.634)   16.033
  10.908   (  -6.254  -15.212    5.452)   17.328
  11.284   (  -4.947   -1.169   -3.592)    6.224
  11.371   (   2.114   -5.278   -3.886)    6.887
  11.535   (  -1.527   -3.178    1.065)    3.683
  11.867   (  -0.549   -2.908   -4.649)    5.511
======================= Grid point 533 (77/135) =======================
q-point: ( 0.33  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.020   (   2.455    9.271  -15.095)   17.884
   3.568   (  -4.923    1.717  -20.029)   20.697
   3.789   ( -13.025   22.235    4.857)   26.223
   4.038   ( -14.320   22.892   -0.687)   27.011
   6.075   (  13.574    2.413   15.714)   20.905
   6.399   (   7.344   -4.264    7.313)   11.207
  10.370   (   6.268   -8.155   13.591)   17.045
  10.655   (  -3.990  -16.778    7.909)   18.973
  11.220   (  -2.116   -0.232   -2.073)    2.971
  11.376   (   3.406   -4.231   -3.867)    6.667
  11.491   (  -1.048   -1.276    1.115)    1.992
  11.843   (   0.285   -1.799   -4.972)    5.295
======================= Grid point 534 (78/135) =======================
q-point: ( 0.40  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.095   (  -4.284    7.420  -16.690)   18.760
   3.530   (  -1.080    1.871  -21.449)   21.558
   3.783   ( -13.976   24.207    5.687)   28.524
   4.019   ( -15.525   26.891    0.027)   31.051
   6.216   (   2.769   -4.797   14.888)   15.885
   6.434   (   4.499   -7.793    9.082)   12.785
  10.442   (   7.795  -13.502   12.391)   19.915
  10.476   (   5.044   -8.736    9.368)   13.766
  11.201   (  -0.255    0.441   -1.384)    1.475
  11.388   (   3.138   -5.435   -4.095)    7.494
  11.475   (  -0.313    0.542    1.420)    1.552
  11.837   (   0.712   -1.234   -5.094)    5.289
======================= Grid point 548 (79/135) =======================
q-point: ( 0.27  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.044   (   5.513    9.549  -15.487)   19.011
   3.605   (  -2.513   -4.352  -22.914)   23.459
   4.037   (   5.548    9.609    5.243)   12.272
   4.309   (   6.560   11.362    0.284)   13.123
   5.951   (   5.906   10.230   16.329)   20.153
   6.299   (   0.679    1.175    7.127)    7.255
  10.269   (  -2.974   -5.150   14.575)   15.742
  10.603   (  -8.457  -14.647    7.494)   18.499
  11.258   (  -1.239   -2.147   -2.256)    3.352
  11.301   (  -0.761   -1.318   -3.716)    4.015
  11.494   (  -0.659   -1.141    1.302)    1.852
  11.833   (  -0.487   -0.844   -5.198)    5.288
======================= Grid point 549 (80/135) =======================
q-point: ( 0.33  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.210   (   1.374    9.647  -18.799)   21.175
   3.516   (  -1.441   -4.999  -23.006)   23.587
   4.147   (  -4.770   12.782    5.865)   14.851
   4.440   (  -5.536   14.971    0.028)   15.962
   6.127   (   5.406    3.288   14.814)   16.109
   6.303   (   1.960   -5.067   10.112)   11.479
  10.225   (   4.082   -5.226   12.952)   14.551
  10.362   (  -3.211  -12.076    9.533)   15.717
  11.221   (  -1.494    0.406   -1.489)    2.148
  11.293   (   2.039   -3.644   -2.738)    4.994
  11.482   (  -0.213   -0.047    1.164)    1.184
  11.825   (   0.060   -0.287   -5.492)    5.500
======================= Grid point 564 (81/135) =======================
q-point: ( 0.33  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 192
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.387   (  -3.978   -6.891  -21.503)   22.928
   3.387   (   3.978    6.891  -21.503)   22.928
   4.274   (   0.000    0.000    5.882)    5.882
   4.590   (   0.000    0.000   -0.423)    0.423
   6.207   (  -2.294   -3.974   12.799)   13.597
   6.207   (   2.294    3.974   12.799)   13.597
  10.203   (  -1.761   -3.050   11.346)   11.880
  10.203   (   1.761    3.050   11.346)   11.880
  11.241   (  -0.861   -1.491   -1.657)    2.390
  11.241   (   0.861    1.491   -1.657)    2.390
  11.482   (  -0.000   -0.000    1.053)    1.053
  11.823   (  -0.000   -0.000   -5.671)    5.671
======================= Grid point 723 (82/135) =======================
q-point: ( 0.00  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.503   (   0.000    0.000   16.788)   16.788
   1.503   (  -0.000    0.000   16.788)   16.788
   1.823   (   0.000    0.000   -0.118)    0.118
   1.823   (   0.000   -0.000   -0.118)    0.118
   3.842   (   0.000    0.000   57.716)   57.716
   5.731   (   0.000    0.000  -38.362)   38.362
  11.709   (   0.000    0.000    6.398)    6.398
  11.945   (   0.000    0.000   -5.144)    5.144
  12.064   (   0.000   -0.000    2.184)    2.184
  12.064   (   0.000    0.000    2.184)    2.184
  12.190   (  -0.000    0.000   -3.887)    3.887
  12.190   (   0.000    0.000   -3.887)    3.887
======================= Grid point 724 (83/135) =======================
q-point: ( 0.07  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.610   (   8.977    5.183   15.548)   18.687
   1.891   (  28.912   16.693   15.150)   36.662
   1.916   (   7.862    4.539   -0.423)    9.088
   2.206   (  29.720   17.159   -1.110)   34.336
   3.856   (   1.858    1.072   54.706)   54.748
   5.675   (  -4.644   -2.681  -37.442)   37.824
  11.602   (  -8.787   -5.073    7.129)   12.400
  11.867   (  -6.684   -3.859   -5.794)    9.651
  12.014   (  -4.270   -2.465    2.184)    5.393
  12.077   (   0.822    0.474    1.592)    1.854
  12.141   (  -4.143   -2.392   -3.964)    6.213
  12.182   (  -0.706   -0.408   -3.453)    3.548
======================= Grid point 725 (84/135) =======================
q-point: ( 0.13  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.883   (  14.121    8.153   12.938)   20.815
   2.157   (  12.454    7.191   -0.981)   14.414
   2.561   (  24.762   14.296   18.403)   34.003
   2.981   (  34.523   19.932   -3.761)   40.041
   3.996   (  12.948    7.475   42.803)   45.339
   5.545   (  -5.759   -3.325  -33.907)   34.553
  11.345   ( -12.888   -7.441    8.472)   17.124
  11.663   ( -10.539   -6.085   -7.112)   14.095
  11.877   (  -7.483   -4.320    2.163)    8.907
  12.008   (  -7.325   -4.229   -4.150)    9.421
  12.083   (  -0.642   -0.371    0.892)    1.160
  12.159   (  -1.193   -0.689   -2.740)    3.067
======================= Grid point 726 (85/135) =======================
q-point: ( 0.20  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.211   (  13.748    7.938   10.418)   18.988
   2.443   (  11.663    6.734   -1.109)   13.513
   2.891   (   5.188    2.995   25.364)   26.062
   3.625   (  18.723   10.809  -10.605)   24.080
   4.509   (  30.295   17.491   24.956)   42.971
   5.498   (   4.246    2.451  -24.535)   25.020
  11.057   ( -11.549   -6.668    9.753)   16.522
  11.419   ( -10.109   -5.836   -8.180)   14.253
  11.689   (  -8.640   -4.988    2.111)   10.197
  11.823   (  -8.465   -4.887   -4.376)   10.709
  12.037   (  -3.429   -1.980    1.546)    4.251
  12.128   (  -1.639   -0.946   -2.832)    3.406
======================= Grid point 727 (86/135) =======================
q-point: ( 0.27  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.471   (   8.350    4.821    8.390)   12.781
   2.647   (   5.545    3.201   -0.261)    6.408
   2.911   (  -1.459   -0.842   20.663)   20.731
   3.773   (  -3.902   -2.253  -19.689)   20.198
   5.247   (  32.174   18.576   16.771)   40.762
   5.788   (  19.366   11.181  -11.334)   25.070
  10.850   (  -6.163   -3.558   10.692)   12.843
  11.228   (  -6.310   -3.643   -8.149)   10.931
  11.503   (  -7.294   -4.211    2.077)    8.674
  11.642   (  -6.976   -4.027   -4.660)    9.306
  11.931   (  -5.676   -3.277    3.407)    7.386
  12.079   (  -2.720   -1.570   -3.693)    4.848
======================= Grid point 728 (87/135) =======================
q-point: ( 0.33  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.576   (   0.964    0.556    6.723)    6.814
   2.679   (  -2.707   -1.563    1.855)    3.635
   2.875   (  -1.186   -0.685   10.816)   10.903
   3.584   ( -10.266   -5.927  -21.918)   24.918
   5.916   (  25.817   14.905   13.629)   32.778
   6.257   (  19.679   11.362   -3.953)   23.065
  10.781   (  -0.038   -0.022   10.993)   10.993
  11.137   (  -1.745   -1.007   -6.889)    7.178
  11.369   (  -4.311   -2.489    2.135)    5.417
  11.518   (  -3.779   -2.182   -5.063)    6.684
  11.785   (  -6.955   -4.015    5.894)    9.961
  12.008   (  -3.339   -1.928   -4.853)    6.198
======================= Grid point 729 (88/135) =======================
q-point: ( 0.40  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.536   (  -3.880   -2.240    5.397)    7.015
   2.553   (  -7.062   -4.078    4.898)    9.513
   2.863   (   0.028    0.016    2.182)    2.182
   3.373   (  -7.486   -4.322  -20.900)   22.617
   6.403   (  16.586    9.576   11.866)   22.529
   6.630   (  12.451    7.189    0.186)   14.379
  10.829   (   3.673    2.120   10.237)   11.080
  11.133   (   0.848    0.490   -4.893)    4.990
  11.304   (  -1.587   -0.916    2.276)    2.922
  11.466   (  -1.038   -0.599   -5.532)    5.660
  11.627   (  -6.580   -3.799    8.695)   11.547
  11.935   (  -2.836   -1.637   -6.099)    6.922
======================= Grid point 730 (89/135) =======================
q-point: ( 0.47  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.405   (  -4.004   -2.312    7.640)    8.930
   2.460   (  -2.279   -1.316    4.410)    5.136
   2.868   (   0.209    0.121   -2.515)    2.527
   3.258   (  -2.526   -1.458  -19.813)   20.026
   6.662   (   5.881    3.395   10.679)   12.656
   6.816   (   3.942    2.276    2.680)    5.283
  10.912   (   2.616    1.510    8.472)    8.994
  11.153   (   0.546    0.315   -2.790)    2.861
  11.285   (  -0.275   -0.159    2.396)    2.417
  11.457   (  -0.040   -0.023   -5.861)    5.861
  11.509   (  -3.103   -1.792   10.946)   11.517
  11.888   (  -1.143   -0.660   -7.022)    7.145
======================= Grid point 739 (90/135) =======================
q-point: ( 0.07  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.835   (   9.528   16.503   13.403)   23.298
   2.113   (   8.421   14.586   -0.709)   16.857
   2.375   (  15.707   27.205   16.417)   35.445
   2.745   (  19.395   33.594   -2.855)   38.896
   3.929   (   4.612    7.989   47.343)   48.233
   5.581   (  -3.600   -6.236  -35.256)   35.984
  11.423   (  -7.060  -12.229    8.078)   16.268
  11.718   (  -6.009  -10.408   -6.109)   13.481
  11.952   (  -2.755   -4.773    1.134)    5.626
  12.025   (  -3.097   -5.363   -0.892)    6.257
  12.106   (  -0.767   -1.329   -1.994)    2.516
  12.146   (  -1.636   -2.834   -2.679)    4.229
======================= Grid point 740 (91/135) =======================
q-point: ( 0.13  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.201   (   8.866   22.686   10.821)   26.652
   2.441   (   8.007   20.476   -1.065)   22.012
   2.791   (  10.323    7.152   23.255)   26.429
   3.391   (  21.737   18.951   -7.310)   29.750
   4.256   (  19.128   19.185   31.260)   41.366
   5.480   (  -1.037   -1.971  -29.023)   29.108
  11.158   (  -9.543  -10.594    9.300)   17.024
  11.492   (  -7.610   -9.969   -6.894)   14.311
  11.806   (  -8.871   -3.560    0.734)    9.587
  11.892   (  -5.212   -7.315   -2.929)    9.447
  12.050   (  -2.982   -3.387   -0.066)    4.513
  12.102   (   0.768   -4.719   -2.404)    5.351
======================= Grid point 741 (92/135) =======================
q-point: ( 0.20  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.531   (   2.296   22.282    9.362)   24.278
   2.739   (   1.142   21.614   -0.798)   21.659
   2.907   (   2.189    0.355   22.417)   22.526
   3.738   (   6.310    0.080  -16.202)   17.387
   4.896   (  29.456   20.632   19.604)   40.959
   5.599   (  12.750    7.911  -16.712)   22.460
  10.923   (  -7.453   -6.252   10.511)   14.322
  11.288   (  -6.140   -6.476   -7.386)   11.584
  11.605   (  -9.484   -4.164    0.981)   10.404
  11.714   (  -6.963   -6.272   -4.083)   10.222
  11.949   (  -2.910   -6.764    2.150)    7.670
  12.061   (   0.148   -5.534   -3.248)    6.419
======================= Grid point 742 (93/135) =======================
q-point: ( 0.27  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.682   (  -3.327   10.623    9.188)   14.434
   2.876   (  -7.115   17.114    1.374)   18.585
   2.952   (  -4.069   11.426   11.253)   16.545
   3.674   (  -6.541   -6.594  -21.089)   23.044
   5.587   (  28.989   15.397   14.974)   36.078
   6.007   (  21.012   10.349   -6.805)   24.391
  10.795   (  -2.684   -1.933   11.412)   11.881
  11.166   (  -3.082   -2.159   -6.931)    7.886
  11.421   (  -6.470   -4.764    1.460)    8.166
  11.553   (  -4.668   -5.327   -4.722)    8.513
  11.819   (  -3.902   -7.594    4.311)    9.564
  11.998   (  -0.587   -6.174   -4.329)    7.564
======================= Grid point 743 (94/135) =======================
q-point: ( 0.33  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.681   (  -5.040    5.768    5.339)    9.337
   2.844   ( -15.682   16.691    1.036)   22.926
   2.966   (  -8.299   13.372    6.125)   16.888
   3.475   (  -9.111   -3.886  -20.368)   22.649
   6.157   (  22.656    9.310   12.725)   27.603
   6.434   (  17.777    6.387   -1.603)   18.957
  10.786   (   2.239    0.077   11.460)   11.677
  11.132   (  -0.356    0.420   -5.663)    5.689
  11.296   (  -2.118   -5.840    2.439)    6.674
  11.459   (  -0.977   -4.105   -5.117)    6.632
  11.675   (  -4.578   -6.126    6.244)    9.873
  11.926   (  -0.681   -5.876   -5.495)    8.074
======================= Grid point 744 (95/135) =======================
q-point: ( 0.40  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.640   (  -7.945    9.088    2.838)   12.401
   2.715   ( -18.247   18.946    3.266)   26.506
   2.950   (  -7.279   11.259    1.024)   13.446
   3.317   (  -6.280   -0.629  -18.863)   19.891
   6.530   (  13.298    3.304   11.181)   17.686
   6.712   (  10.032    1.182    1.732)   10.249
  10.851   (   4.747   -0.389   10.062)   11.132
  11.144   (   0.818    0.312   -4.265)    4.354
  11.221   (   1.221   -8.326    4.848)    9.712
  11.431   (   1.335   -2.648   -5.383)    6.146
  11.557   (  -3.840   -2.503    7.161)    8.503
  11.867   (   0.173   -4.701   -6.521)    8.040
======================= Grid point 745 (96/135) =======================
q-point: ( 0.47  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.612   (  -7.274   12.599    2.478)   14.758
   2.645   ( -13.661   23.661    4.565)   27.701
   2.944   (  -5.865   10.159   -1.388)   11.812
   3.259   (  -1.254    2.172  -18.294)   18.465
   6.663   (   1.876   -3.249   10.420)   11.075
   6.802   (   2.077   -3.598    3.158)    5.218
  10.901   (   1.642   -2.845    8.426)    9.044
  11.154   (   0.635   -1.100   -3.878)    4.080
  11.183   (   4.779   -8.278    7.543)   12.176
  11.435   (   1.288   -2.231   -6.362)    6.864
  11.508   (  -0.618    1.071    7.870)    7.967
  11.844   (   1.912   -3.311   -6.959)    7.940
======================= Grid point 756 (97/135) =======================
q-point: ( 0.13  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.689   (  13.960   24.179    8.428)   29.164
   2.820   (   2.034    3.522   17.875)   18.332
   2.910   (  10.690   18.515    4.438)   21.835
   3.649   (   3.479    6.026  -14.378)   15.973
   4.729   (  15.877   27.499   21.598)   38.402
   5.527   (   4.645    8.046  -19.399)   21.509
  10.960   (  -5.119   -8.866   10.457)   14.634
  11.308   (  -4.670   -8.088   -6.424)   11.336
  11.711   (  -4.297   -7.442   -1.536)    8.730
  11.720   (  -5.218   -9.038   -1.548)   10.550
  11.937   (  -4.320   -7.483    0.744)    8.672
  12.012   (  -2.292   -3.969   -2.973)    5.463
======================= Grid point 757 (98/135) =======================
q-point: ( 0.20  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.780   (  -2.610   -0.867   15.815)   16.052
   3.105   (   4.221   23.695    6.538)   24.941
   3.259   (  -1.271   26.364   -0.710)   26.404
   3.655   (   0.450   -6.847  -19.249)   20.436
   5.298   (  22.254   19.858   16.590)   34.129
   5.798   (  14.199   11.692   -9.544)   20.722
  10.823   (  -3.407   -3.842   11.989)   13.042
  11.188   (  -2.567   -3.510   -5.522)    7.029
  11.501   (  -9.028   -6.773   -0.988)   11.329
  11.562   (  -4.249   -9.504   -3.655)   11.034
  11.792   (  -3.141   -8.157    3.015)    9.246
  11.937   (  -0.736   -6.241   -4.072)    7.488
======================= Grid point 758 (99/135) =======================
q-point: ( 0.27  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.743   (  -0.911    0.009    7.473)    7.528
   3.246   (  -9.972   13.191   -6.804)   17.882
   3.380   ( -10.864   24.724    3.116)   27.184
   3.616   (  -5.678    7.335  -11.933)   15.114
   5.849   (  23.252   11.247   13.858)   29.312
   6.181   (  17.964    7.265   -3.447)   19.682
  10.763   (  -0.575   -1.945   13.272)   13.427
  11.145   (  -0.586   -0.518   -3.831)    3.910
  11.288   (  -6.398   -9.078   -0.180)   11.107
  11.420   (  -1.413   -7.953   -4.831)    9.412
  11.671   (  -2.648   -6.369    4.086)    8.017
  11.865   (  -0.221   -6.391   -5.030)    8.136
======================= Grid point 759 (100/135) =======================
q-point: ( 0.33  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.748   (  -0.343    2.680    0.538)    2.755
   3.149   ( -10.042    6.877  -13.294)   18.024
   3.393   ( -16.728   26.169    5.057)   31.468
   3.558   ( -15.362   17.467   -6.430)   24.133
   6.288   (  18.410    4.226   11.958)   22.355
   6.509   (  14.397    1.334    0.476)   14.466
  10.759   (   3.235   -3.503   12.790)   13.650
  11.095   (  -1.249  -12.333    2.758)   12.699
  11.157   (  -0.371   -0.186   -2.352)    2.389
  11.364   (   1.984   -4.948   -5.790)    7.871
  11.579   (  -2.326   -3.373    4.379)    5.997
  11.808   (   0.369   -5.204   -5.835)    7.827
======================= Grid point 760 (101/135) =======================
q-point: ( 0.40  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.769   (  -1.821    4.612   -2.987)    5.789
   3.074   (  -7.480    8.966  -13.172)   17.602
   3.358   ( -17.154   27.070    4.121)   32.312
   3.473   ( -15.988   21.932   -4.439)   27.502
   6.539   (   9.115   -2.107   10.488)   14.054
   6.681   (   7.228   -4.093    3.067)    8.855
  10.797   (   5.637   -5.241   10.029)   12.642
  10.967   (   3.400  -15.113    6.841)   16.934
  11.161   (   0.033    0.415   -1.243)    1.311
  11.368   (   2.797   -3.836   -6.712)    8.222
  11.528   (  -1.232   -0.673    4.749)    4.953
  11.778   (   1.286   -3.750   -6.328)    7.468
======================= Grid point 772 (102/135) =======================
q-point: ( 0.20  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.748   (  -0.503   -0.872    8.652)    8.710
   3.386   (  -1.806   -3.128  -16.149)   16.548
   3.631   (   7.402   12.821    0.036)   14.804
   3.734   (  10.487   18.164   -0.851)   20.991
   5.682   (  11.012   19.073   14.620)   26.434
   6.045   (   7.552   13.080   -4.420)   15.737
  10.760   (  -1.753   -3.037   14.287)   14.711
  11.142   (  -1.094   -1.895   -1.401)    2.598
  11.327   (  -6.524  -11.299   -4.735)   13.880
  11.366   (  -4.464   -7.731   -3.313)    9.522
  11.659   (  -3.017   -5.226    3.442)    6.947
  11.832   (  -2.455   -4.253   -4.886)    6.927
======================= Grid point 773 (103/135) =======================
q-point: ( 0.27  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.762   (   1.148    2.369    1.117)    2.860
   3.259   (  -4.810   -4.435  -20.813)   21.817
   3.844   (  -4.500   20.590    5.313)   21.736
   3.981   (  -4.322   20.750   -1.246)   21.232
   6.046   (  14.750    9.241   12.807)   21.610
   6.301   (  11.177    5.195   -0.528)   12.336
  10.691   (  -1.060   -5.737   15.097)   16.185
  11.015   (  -5.173  -13.639    4.721)   15.332
  11.187   (  -2.555    0.606   -6.155)    6.692
  11.276   (   1.589   -5.678   -5.494)    8.059
  11.579   (  -1.666   -3.181    3.410)    4.952
  11.767   (  -0.768   -3.442   -5.469)    6.507
======================= Grid point 774 (104/135) =======================
q-point: ( 0.33  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.819   (   0.831    4.671   -4.421)    6.485
   3.161   (  -2.908   -1.975  -19.458)   19.773
   3.899   ( -13.148   23.353    5.534)   27.365
   4.030   ( -14.521   24.619   -0.930)   28.598
   6.340   (  12.302    1.552   10.936)   16.533
   6.500   (   9.243   -1.862    2.562)    9.770
  10.635   (   3.633   -8.579   12.734)   15.778
  10.801   (  -1.921  -15.301    6.040)   16.562
  11.171   (  -0.718    0.793   -2.411)    2.637
  11.279   (   2.972   -3.463   -5.876)    7.440
  11.536   (  -0.875   -1.268    3.455)    3.783
  11.735   (   0.217   -2.232   -5.894)    6.306
======================= Grid point 775 (105/135) =======================
q-point: ( 0.40  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.855   (  -2.457    4.255   -6.638)    8.258
   3.125   (   0.128   -0.221  -18.466)   18.468
   3.896   ( -14.463   25.050    5.463)   29.437
   4.020   ( -15.789   27.348   -0.626)   31.585
   6.461   (   3.013   -5.219    9.780)   11.488
   6.568   (   3.889   -6.736    4.151)    8.816
  10.658   (   4.996   -8.654    8.560)   13.158
  10.671   (   7.421  -12.853    9.867)   17.822
  11.170   (  -0.305    0.529   -1.363)    1.493
  11.289   (   2.189   -3.791   -5.812)    7.276
  11.522   (   0.043   -0.074    3.525)    3.526
  11.727   (   0.878   -1.521   -6.058)    6.307
======================= Grid point 789 (106/135) =======================
q-point: ( 0.27  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.836   (   2.748    4.759   -4.734)    7.253
   3.175   (  -2.206   -3.820  -19.871)   20.355
   4.149   (   5.865   10.158    5.730)   13.054
   4.299   (   6.291   10.895   -1.512)   12.672
   6.225   (   5.028    8.708   11.275)   15.108
   6.395   (   2.448    4.240    2.299)    5.409
  10.555   (  -3.926   -6.799   13.753)   15.836
  10.740   (  -7.471  -12.941    5.476)   15.914
  11.202   (  -0.627   -1.086   -3.291)    3.522
  11.213   (   0.238    0.413   -4.888)    4.911
  11.536   (  -0.751   -1.300    3.145)    3.485
  11.723   (  -0.717   -1.242   -5.820)    5.994
======================= Grid point 790 (107/135) =======================
q-point: ( 0.33  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.926   (   0.647    5.779   -9.142)   10.835
   3.111   (  -0.941   -3.122  -17.371)   17.674
   4.265   (  -4.689   12.953    5.949)   15.005
   4.425   (  -5.151   14.141   -1.788)   15.156
   6.367   (   4.573    1.565    9.333)   10.510
   6.453   (   2.857   -2.544    4.810)    6.145
  10.469   (   2.503   -6.774   10.958)   13.124
  10.538   (  -1.811  -10.515    7.436)   13.005
  11.185   (  -0.803    0.633   -1.913)    2.169
  11.227   (   1.567   -1.892   -3.784)    4.511
  11.522   (  -0.171   -0.305    3.009)    3.029
  11.710   (   0.056   -0.573   -6.062)    6.090
======================= Grid point 805 (108/135) =======================
q-point: ( 0.33  0.33  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.033   (  -2.390   -4.140  -13.597)   14.413
   3.033   (   2.390    4.140  -13.597)   14.413
   4.394   (   0.000    0.000    6.109)    6.109
   4.566   (   0.000    0.000   -2.231)    2.231
   6.405   (  -1.107   -1.918    7.164)    7.499
   6.405   (   1.107    1.918    7.164)    7.499
  10.413   (  -0.851   -1.473    9.053)    9.212
  10.413   (   0.851    1.473    9.053)    9.212
  11.200   (  -0.486   -0.842   -2.282)    2.481
  11.200   (   0.486    0.842   -2.282)    2.481
  11.519   (  -0.000   -0.000    2.861)    2.861
  11.705   (  -0.000   -0.000   -6.158)    6.158
======================= Grid point 964 (109/135) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 81
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.746   (   0.000   -0.000   -8.146)    8.146
   1.746   (   0.000    0.000   -8.146)    8.146
   1.746   (  -0.000    0.000    8.146)    8.146
   1.746   (   0.000   -0.000    8.146)    8.146
   4.881   (   0.000    0.000  -48.880)   48.880
   4.881   (   0.000   -0.000   48.880)   48.880
  11.834   (   0.000   -0.000   -6.201)    6.201
  11.834   (   0.000    0.000    6.201)    6.201
  12.118   (  -0.000    0.000   -3.330)    3.330
  12.118   (   0.000    0.000   -3.330)    3.330
  12.118   (   0.000    0.000    3.330)    3.330
  12.118   (   0.000   -0.000    3.330)    3.330
======================= Grid point 966 (110/135) =======================
q-point: ( 0.07  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.839   (   7.840    4.527   -7.849)   11.982
   1.839   (   7.840    4.527    7.849)   11.982
   2.119   (  28.702   16.571   -8.106)   34.119
   2.119   (  28.702   16.571    8.106)   34.119
   4.850   (  -2.321   -1.340  -47.168)   47.244
   4.850   (  -2.321   -1.340   47.168)   47.244
  11.741   (  -7.708   -4.450   -6.976)   11.309
  11.741   (  -7.708   -4.450    6.976)   11.309
  12.069   (  -4.241   -2.448   -3.372)    5.945
  12.069   (  -4.241   -2.448    3.372)    5.945
  12.120   (   0.006    0.003   -2.786)    2.786
  12.120   (   0.006    0.003    2.786)    2.786
======================= Grid point 968 (111/135) =======================
q-point: ( 0.13  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.080   (  12.610    7.281   -7.102)   16.201
   2.080   (  12.610    7.281    7.102)   16.201
   2.843   (  30.545   17.635  -10.668)   36.848
   2.843   (  30.545   17.635   10.668)   36.848
   4.816   (   0.988    0.570  -40.554)   40.570
   4.816   (   0.988    0.570   40.554)   40.570
  11.511   ( -11.713   -6.763   -8.359)   15.900
  11.511   ( -11.713   -6.763    8.359)   15.900
  11.932   (  -7.460   -4.307   -3.462)    9.284
  11.932   (  -7.460   -4.307    3.462)    9.284
  12.111   (  -0.853   -0.493   -2.025)    2.252
  12.111   (  -0.853   -0.493    2.025)    2.252
======================= Grid point 970 (112/135) =======================
q-point: ( 0.20  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.374   (  12.274    7.087   -6.062)   15.415
   2.374   (  12.274    7.087    6.062)   15.415
   3.334   (  10.076    5.818  -19.345)   22.575
   3.334   (  10.076    5.818   19.345)   22.575
   5.005   (  17.288    9.981  -25.792)   32.615
   5.005   (  17.288    9.981   25.792)   32.615
  11.246   ( -10.709   -6.183   -9.513)   15.602
  11.246   ( -10.709   -6.183    9.513)   15.602
  11.744   (  -8.620   -4.977   -3.560)   10.571
  11.744   (  -8.620   -4.977    3.560)   10.571
  12.076   (  -2.354   -1.359   -2.435)    3.649
  12.076   (  -2.354   -1.359    2.435)    3.649
======================= Grid point 972 (113/135) =======================
q-point: ( 0.27  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.600   (   6.872    3.968   -4.657)    9.201
   2.600   (   6.872    3.968    4.657)    9.201
   3.348   (  -5.874   -3.392  -23.212)   24.183
   3.348   (  -5.874   -3.392   23.212)   24.183
   5.545   (  27.073   15.630  -13.831)   34.184
   5.545   (  27.073   15.630   13.831)   34.184
  11.052   (  -5.931   -3.424   -9.873)   12.016
  11.052   (  -5.931   -3.424    9.873)   12.016
  11.559   (  -7.210   -4.163   -3.704)    9.112
  11.559   (  -7.210   -4.163    3.704)    9.112
  12.003   (  -3.932   -2.270   -3.911)    5.992
  12.003   (  -3.932   -2.270    3.911)    5.992
======================= Grid point 974 (114/135) =======================
q-point: ( 0.33  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.671   (  -0.606   -0.350   -2.743)    2.831
   2.671   (  -0.606   -0.350    2.743)    2.831
   3.174   (  -7.890   -4.556  -19.191)   21.244
   3.174   (  -7.890   -4.556   19.191)   21.244
   6.135   (  23.386   13.502   -8.748)   28.385
   6.135   (  23.386   13.502    8.748)   28.385
  10.980   (  -0.546   -0.315   -9.275)    9.296
  10.980   (  -0.546   -0.315    9.275)    9.296
  11.429   (  -4.124   -2.381   -3.967)    6.198
  11.429   (  -4.124   -2.381    3.967)    6.198
  11.902   (  -4.810   -2.777   -5.913)    8.112
  11.902   (  -4.810   -2.777    5.913)    8.112
======================= Grid point 976 (115/135) =======================
q-point: ( 0.40  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.596   (  -5.024   -2.901   -0.532)    5.825
   2.596   (  -5.024   -2.901    0.532)    5.825
   3.024   (  -4.979   -2.874  -13.937)   15.076
   3.024   (  -4.979   -2.874   13.937)   15.076
   6.577   (  14.912    8.610   -5.827)   18.179
   6.577   (  14.912    8.610    5.827)   18.179
  11.008   (   2.397    1.384   -7.915)    8.385
  11.008   (   2.397    1.384    7.915)    8.385
  11.368   (  -1.378   -0.795   -4.314)    4.598
  11.368   (  -1.378   -0.795    4.314)    4.598
  11.795   (  -4.312   -2.489   -8.205)    9.597
  11.795   (  -4.312   -2.489    8.205)    9.597
======================= Grid point 978 (116/135) =======================
q-point: ( 0.47  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.494   (  -2.898   -1.673   -1.277)    3.582
   2.494   (  -2.898   -1.673    1.277)    3.582
   2.950   (  -1.553   -0.897  -10.759)   10.907
   2.950   (  -1.553   -0.897   10.759)   10.907
   6.806   (   5.059    2.921   -3.938)    7.045
   6.806   (   5.059    2.921    3.938)    7.045
  11.060   (   1.526    0.881   -6.437)    6.674
  11.060   (   1.526    0.881    6.437)    6.674
  11.353   (  -0.183   -0.106   -4.569)    4.574
  11.353   (  -0.183   -0.106    4.569)    4.574
  11.720   (  -1.871   -1.080  -10.105)   10.333
  11.720   (  -1.871   -1.080   10.105)   10.333
======================= Grid point 996 (117/135) =======================
q-point: ( 0.07  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.038   (   8.526   14.768   -7.207)   18.513
   2.038   (   8.526   14.768    7.207)   18.513
   2.627   (  17.777   30.791   -9.476)   36.796
   2.627   (  17.777   30.791    9.476)   36.796
   4.815   (  -0.676   -1.172  -43.085)   43.106
   4.815   (  -0.676   -1.172   43.085)   43.106
  11.580   (  -6.442  -11.158   -7.823)   15.073
  11.580   (  -6.442  -11.158    7.823)   15.073
  11.989   (  -3.352   -5.805   -2.377)    7.112
  11.989   (  -3.352   -5.805    2.377)    7.112
  12.103   (  -0.883   -1.529   -1.429)    2.272
  12.103   (  -0.883   -1.529    1.429)    2.272
======================= Grid point 998 (118/135) =======================
q-point: ( 0.13  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.370   (   8.091   20.947   -6.285)   23.319
   2.370   (   8.091   20.947    6.285)   23.319
   3.172   (  16.102   12.844  -15.524)   25.792
   3.172   (  16.102   12.844   15.524)   25.792
   4.879   (   7.582    7.240  -32.155)   33.821
   4.879   (   7.582    7.240   32.155)   33.821
  11.337   (  -8.658   -9.907   -8.849)   15.856
  11.337   (  -8.658   -9.907    8.849)   15.856
  11.838   (  -7.435   -5.519   -2.385)    9.562
  11.838   (  -7.435   -5.519    2.385)    9.562
  12.063   (  -0.526   -4.217   -1.421)    4.480
  12.063   (  -0.526   -4.217    1.421)    4.480
======================= Grid point 1000 (119/135) =======================
q-point: ( 0.20  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.677   (   1.668   22.187   -5.447)   22.907
   2.677   (   1.668   22.187    5.447)   22.907
   3.356   (   1.407   -3.125  -22.414)   22.674
   3.356   (   1.407   -3.125   22.414)   22.674
   5.262   (  22.307   15.391  -18.013)   32.542
   5.262   (  22.307   15.391   18.013)   32.542
  11.121   (  -6.844   -5.728   -9.624)   13.126
  11.121   (  -6.844   -5.728    9.624)   13.126
  11.644   (  -8.180   -5.521   -2.961)   10.304
  11.644   (  -8.180   -5.521    2.961)   10.304
  11.999   (  -1.047   -6.134   -2.950)    6.887
  11.999   (  -1.047   -6.134    2.950)    6.887
======================= Grid point 1002 (120/135) =======================
q-point: ( 0.27  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.831   (  -6.442   16.485   -5.257)   18.463
   2.831   (  -6.442   16.485    5.257)   18.463
   3.265   (  -6.110   -2.478  -19.754)   20.825
   3.265   (  -6.110   -2.478   19.754)   20.825
   5.839   (  25.770   13.483  -10.736)   31.003
   5.839   (  25.770   13.483   10.736)   31.003
  11.003   (  -2.885   -1.271   -9.740)   10.237
  11.003   (  -2.885   -1.271    9.740)   10.237
  11.470   (  -5.404   -5.282   -3.570)    8.357
  11.470   (  -5.404   -5.282    3.570)    8.357
  11.908   (  -1.909   -6.920   -4.757)    8.611
  11.908   (  -1.909   -6.920    4.757)    8.611
======================= Grid point 1004 (121/135) =======================
q-point: ( 0.33  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.803   ( -10.430    9.573   -5.399)   15.152
   2.803   ( -10.430    9.573    5.399)   15.152
   3.142   (  -9.519    4.811  -12.881)   16.723
   3.142   (  -9.519    4.811   12.881)   16.723
   6.351   (  20.664    8.186   -7.131)   23.342
   6.351   (  20.664    8.186    7.131)   23.342
  10.988   (   0.851    0.787   -9.069)    9.142
  10.988   (   0.851    0.787    9.069)    9.142
  11.363   (  -1.417   -4.652   -4.367)    6.536
  11.363   (  -1.417   -4.652    4.367)    6.536
  11.805   (  -2.206   -6.361   -6.719)    9.512
  11.805   (  -2.206   -6.361    6.719)    9.512
======================= Grid point 1006 (122/135) =======================
q-point: ( 0.40  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.711   ( -11.031   10.646   -3.449)   15.714
   2.711   ( -11.031   10.646    3.449)   15.714
   3.042   (  -8.798    7.520   -8.789)   14.533
   3.042   (  -8.798    7.520    8.789)   14.533
   6.685   (  11.922    2.426   -4.681)   13.036
   6.685   (  11.922    2.426    4.681)   13.036
  11.026   (   2.532   -0.178   -7.670)    8.079
  11.026   (   2.532   -0.178    7.670)    8.079
  11.321   (   1.205   -4.125   -5.464)    6.952
  11.321   (   1.205   -4.125    5.464)    6.952
  11.719   (  -1.279   -4.470   -8.500)    9.689
  11.719   (  -1.279   -4.470    8.500)    9.689
======================= Grid point 1008 (123/135) =======================
q-point: ( 0.47  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.664   (  -8.325   14.420   -2.222)   16.798
   2.664   (  -8.325   14.420    2.222)   16.798
   3.002   (  -5.148    8.917   -7.524)   12.753
   3.002   (  -5.148    8.917    7.524)   12.753
   6.800   (   1.968   -3.408   -3.548)    5.298
   6.800   (   1.968   -3.408    3.548)    5.298
  11.050   (   1.259   -2.181   -6.675)    7.135
  11.050   (   1.259   -2.181    6.675)    7.135
  11.313   (   2.203   -3.816   -6.225)    7.626
  11.313   (   2.203   -3.816    6.225)    7.626
  11.684   (   1.265   -2.191   -9.290)    9.628
  11.684   (   1.265   -2.191    9.290)    9.628
======================= Grid point 1030 (124/135) =======================
q-point: ( 0.13  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.824   (  12.966   22.458   -5.059)   26.421
   2.824   (  12.966   22.458    5.059)   26.421
   3.295   (   0.077    0.133  -21.455)   21.455
   3.295   (   0.077    0.133   21.455)   21.455
   5.139   (  10.924   18.921  -20.527)   29.978
   5.139   (  10.924   18.921   20.527)   29.978
  11.154   (  -4.594   -7.957   -9.222)   13.018
  11.154   (  -4.594   -7.957    9.222)   13.018
  11.699   (  -5.004   -8.668   -0.353)   10.015
  11.699   (  -5.004   -8.668    0.353)   10.015
  11.963   (  -3.134   -5.428   -1.957)    6.567
  11.963   (  -3.134   -5.428    1.957)    6.567
======================= Grid point 1032 (125/135) =======================
q-point: ( 0.20  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.094   (  -5.769    8.463  -11.882)   15.687
   3.094   (  -5.769    8.463   11.882)   15.687
   3.319   (   4.147    9.672  -11.236)   15.395
   3.319   (   4.147    9.672   11.236)   15.395
   5.584   (  18.960   16.718  -12.734)   28.304
   5.584   (  18.960   16.718   12.734)   28.304
  11.038   (  -2.837   -2.520   -9.713)   10.428
  11.038   (  -2.837   -2.520    9.713)   10.428
  11.508   (  -6.610   -8.698   -1.707)   11.057
  11.508   (  -6.610   -8.698    1.707)   11.057
  11.859   (  -1.713   -7.197   -3.800)    8.317
  11.859   (  -1.713   -7.197    3.800)    8.317
======================= Grid point 1034 (126/135) =======================
q-point: ( 0.27  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.979   (  -6.329   -2.481  -15.613)   17.029
   2.979   (  -6.329   -2.481   15.613)   17.029
   3.473   (  -8.570   22.939   -5.226)   25.039
   3.473   (  -8.570   22.939    5.226)   25.039
   6.067   (  21.034    9.754   -8.503)   24.696
   6.067   (  21.034    9.754    8.503)   24.696
  10.997   (  -1.178    0.182  -10.419)   10.487
  10.997   (  -1.178    0.182   10.419)   10.487
  11.332   (  -3.010   -8.916   -3.876)   10.178
  11.332   (  -3.010   -8.916    3.876)   10.178
  11.764   (  -1.287   -6.731   -5.212)    8.610
  11.764   (  -1.287   -6.731    5.212)    8.610
======================= Grid point 1036 (127/135) =======================
q-point: ( 0.33  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.875   (  -4.071   -0.318  -11.893)   12.574
   2.875   (  -4.071   -0.318   11.893)   12.574
   3.479   ( -17.465   25.374   -3.627)   31.016
   3.479   ( -17.465   25.374    3.627)   31.016
   6.460   (  16.681    3.068   -5.667)   17.883
   6.460   (  16.681    3.068    5.667)   17.883
  10.981   (   0.564   -3.021   -9.655)   10.133
  10.981   (   0.564   -3.021    9.655)   10.133
  11.243   (   1.362   -6.409   -6.320)    9.103
  11.243   (   1.362   -6.409    6.320)    9.103
  11.687   (  -0.823   -4.781   -6.275)    7.931
  11.687   (  -0.823   -4.781    6.275)    7.931
======================= Grid point 1038 (128/135) =======================
q-point: ( 0.40  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.831   (  -2.480    2.327   -9.124)    9.737
   2.831   (  -2.480    2.327    9.124)    9.737
   3.424   ( -18.383   27.807   -2.460)   33.425
   3.424   ( -18.383   27.807    2.460)   33.425
   6.678   (   8.275   -2.988   -3.620)    9.514
   6.678   (   8.275   -2.988    3.620)    9.514
  10.962   (   3.477   -7.264   -6.193)   10.159
  10.962   (   3.477   -7.264    6.193)   10.159
  11.229   (   2.197   -3.815   -7.012)    8.280
  11.229   (   2.197   -3.815    7.012)    8.280
  11.646   (   0.280   -2.641   -6.899)    7.393
  11.646   (   0.280   -2.641    6.899)    7.393
======================= Grid point 1062 (129/135) =======================
q-point: ( 0.20  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.017   (  -3.906   -6.766  -18.049)   19.667
   3.017   (  -3.906   -6.766   18.049)   19.667
   3.692   (  10.332   17.896   -3.322)   20.930
   3.692   (  10.332   17.896    3.322)   20.930
   5.915   (   9.645   16.706   -9.251)   21.394
   5.915   (   9.645   16.706    9.251)   21.394
  11.014   (  -0.391   -0.676  -11.065)   11.093
  11.014   (  -0.391   -0.676   11.065)   11.093
  11.308   (  -6.016  -10.420   -2.090)   12.212
  11.308   (  -6.016  -10.420    2.090)   12.212
  11.738   (  -2.826   -4.894   -4.611)    7.294
  11.738   (  -2.826   -4.894    4.611)    7.294
======================= Grid point 1064 (130/135) =======================
q-point: ( 0.27  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.911   (  -2.551   -2.853  -13.745)   14.268
   2.911   (  -2.551   -2.853   13.745)   14.268
   3.934   (  -4.376   20.867   -3.636)   21.629
   3.934   (  -4.376   20.867    3.636)   21.629
   6.235   (  13.186    7.607   -6.489)   16.548
   6.235   (  13.186    7.607    6.489)   16.548
  10.955   (  -3.306   -6.950  -11.158)   13.554
  10.955   (  -3.306   -6.950   11.158)   13.554
  11.171   (   0.096   -4.389   -4.818)    6.518
  11.171   (   0.096   -4.389    4.818)    6.518
  11.663   (  -1.256   -3.531   -5.058)    6.295
  11.663   (  -1.256   -3.531    5.058)    6.295
======================= Grid point 1066 (131/135) =======================
q-point: ( 0.33  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.870   (  -1.181    0.292   -9.461)    9.539
   2.870   (  -1.181    0.292    9.461)    9.539
   3.988   ( -13.710   23.962   -3.456)   27.822
   3.988   ( -13.710   23.962    3.456)   27.822
   6.487   (  10.888    0.045   -4.069)   11.624
   6.487   (  10.888    0.045    4.069)   11.624
  10.822   (   0.441  -11.719   -5.099)   12.788
  10.822   (   0.441  -11.719    5.099)   12.788
  11.174   (   1.615   -0.846   -3.643)    4.073
  11.174   (   1.615   -0.846    3.643)    4.073
  11.623   (  -0.288   -1.956   -5.383)    5.735
  11.623   (  -0.288   -1.956    5.383)    5.735
======================= Grid point 1068 (132/135) =======================
q-point: ( 0.40  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.863   (  -0.798    1.381   -7.471)    7.640
   2.863   (  -0.798    1.381    7.471)    7.640
   3.983   ( -15.050   26.067   -3.273)   30.277
   3.983   ( -15.050   26.067    3.273)   30.277
   6.583   (   3.421   -5.925   -2.697)    7.354
   6.583   (   3.421   -5.925    2.697)    7.354
  10.766   (   6.179  -10.703   -1.061)   12.404
  10.766   (   6.179  -10.703    1.061)   12.404
  11.186   (   0.669   -1.159   -3.711)    3.945
  11.186   (   0.669   -1.159    3.711)    3.945
  11.612   (   0.577   -0.999   -5.521)    5.640
  11.612   (   0.577   -0.999    5.521)    5.640
======================= Grid point 1096 (133/135) =======================
q-point: ( 0.27  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.881   (  -0.189   -0.328   -9.337)    9.345
   2.881   (  -0.189   -0.328    9.337)    9.345
   4.246   (   6.050   10.479   -3.988)   12.740
   4.246   (   6.050   10.479    3.988)   12.740
   6.377   (   3.861    6.687   -4.328)    8.852
   6.377   (   3.861    6.687    4.328)    8.852
  10.755   (  -6.031  -10.446   -5.301)   13.176
  10.755   (  -6.031  -10.446    5.301)   13.176
  11.156   (   0.498    0.863   -0.203)    1.017
  11.156   (   0.498    0.863    0.203)    1.017
  11.616   (  -0.816   -1.413   -4.977)    5.237
  11.616   (  -0.816   -1.413    4.977)    5.237
======================= Grid point 1098 (134/135) =======================
q-point: ( 0.33  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.885   (  -0.195    0.959   -5.093)    5.186
   2.885   (  -0.195    0.959    5.093)    5.186
   4.366   (  -4.839   13.407   -4.219)   14.865
   4.366   (  -4.839   13.407    4.219)   14.865
   6.480   (   3.733   -0.405   -2.161)    4.332
   6.480   (   3.733   -0.405    2.161)    4.332
  10.603   (   0.275   -8.923   -1.788)    9.104
  10.603   (   0.275   -8.923    1.788)    9.104
  11.171   (   0.465    0.079   -1.187)    1.278
  11.171   (   0.465    0.079    1.187)    1.278
  11.600   (  -0.049   -0.566   -5.012)    5.044
  11.600   (  -0.049   -0.566    5.012)    5.044
======================= Grid point 1128 (135/135) =======================
q-point: ( 0.33  0.33 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 2.01e-04 2.01e-04 5.87e-04 5.87e-04 
Number of triplets: 119
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.896   (  -0.000   -0.000    0.000)    0.000
   2.896   (  -0.000   -0.000    0.000)    0.000
   4.500   (   0.000    0.000   -4.539)    4.539
   4.500   (   0.000    0.000    4.539)    4.539
   6.476   (  -0.000   -0.000    0.000)    0.000
   6.476   (  -0.000   -0.000    0.000)    0.000
  10.510   (  -0.000   -0.000    0.000)    0.000
  10.510   (  -0.000   -0.000    0.000)    0.000
  11.172   (  -0.000   -0.000    0.000)    0.000
  11.172   (  -0.000   -0.000    0.000)    0.000
  11.595   (  -0.000   -0.000   -4.936)    4.936
  11.595   (  -0.000   -0.000    4.936)    4.936
=================== 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   2120.525   2120.525   1011.589      0.000     -0.000     -0.000 3/21600
   20.0   1134.901   1134.901    707.280      0.000     -0.000     -0.000 3/21600
   30.0    632.560    632.560    455.439      0.000     -0.000     -0.000 3/21600
   40.0    483.794    483.794    364.710      0.000     -0.000     -0.000 3/21600
   50.0    408.477    408.477    311.637      0.000     -0.000     -0.000 3/21600
   60.0    358.472    358.472    273.549      0.000     -0.000     -0.000 3/21600
   70.0    320.235    320.235    243.116      0.000     -0.000      0.000 3/21600
   80.0    288.964    288.964    217.602      0.000     -0.000      0.000 3/21600
   90.0    262.711    262.711    195.961      0.000     -0.000      0.000 3/21600
  100.0    240.480    240.480    177.648      0.000     -0.000      0.000 3/21600
  110.0    221.584    221.584    162.199      0.000     -0.000      0.000 3/21600
  120.0    205.462    205.462    149.158      0.000     -0.000      0.000 3/21600
  130.0    191.632    191.632    138.100      0.000     -0.000      0.000 3/21600
  140.0    179.688    179.688    128.660      0.000     -0.000      0.000 3/21600
  150.0    169.295    169.295    120.535      0.000     -0.000      0.000 3/21600
  160.0    160.180    160.180    113.481      0.000     -0.000      0.000 3/21600
  170.0    152.125    152.125    107.304      0.000     -0.000      0.000 3/21600
  180.0    144.953    144.953    101.850      0.000     -0.000      0.000 3/21600
  190.0    138.524    138.524     96.999      0.000     -0.000      0.000 3/21600
  200.0    132.724    132.724     92.653      0.000     -0.000      0.000 3/21600
  210.0    127.460    127.460     88.734      0.000     -0.000      0.000 3/21600
  220.0    122.658    122.658     85.180      0.000     -0.000      0.000 3/21600
  230.0    118.255    118.255     81.938      0.000     -0.000      0.000 3/21600
  240.0    114.200    114.200     78.969      0.000     -0.000      0.000 3/21600
  250.0    110.451    110.451     76.235      0.000     -0.000      0.000 3/21600
  260.0    106.971    106.971     73.710      0.000     -0.000      0.000 3/21600
  270.0    103.731    103.731     71.367      0.000     -0.000      0.000 3/21600
  280.0    100.705    100.705     69.187      0.000     -0.000      0.000 3/21600
  290.0     97.870     97.870     67.152      0.000     -0.000      0.000 3/21600
  300.0     95.208     95.208     65.247      0.000     -0.000      0.000 3/21600
  310.0     92.701     92.701     63.459      0.000     -0.000      0.000 3/21600
  320.0     90.337     90.337     61.777      0.000     -0.000      0.000 3/21600
  330.0     88.102     88.102     60.191      0.000     -0.000      0.000 3/21600
  340.0     85.985     85.985     58.692      0.000     -0.000      0.000 3/21600
  350.0     83.976     83.976     57.273      0.000     -0.000      0.000 3/21600
  360.0     82.066     82.066     55.927      0.000     -0.000      0.000 3/21600
  370.0     80.249     80.249     54.649      0.000     -0.000      0.000 3/21600
  380.0     78.517     78.517     53.433      0.000     -0.000      0.000 3/21600
  390.0     76.863     76.863     52.274      0.000     -0.000      0.000 3/21600
  400.0     75.283     75.283     51.169      0.000     -0.000      0.000 3/21600
  410.0     73.770     73.770     50.112      0.000     -0.000      0.000 3/21600
  420.0     72.322     72.322     49.102      0.000     -0.000      0.000 3/21600
  430.0     70.933     70.933     48.134      0.000     -0.000      0.000 3/21600
  440.0     69.599     69.599     47.207      0.000     -0.000      0.000 3/21600
  450.0     68.318     68.318     46.317      0.000     -0.000      0.000 3/21600
  460.0     67.085     67.085     45.462      0.000     -0.000      0.000 3/21600
  470.0     65.899     65.899     44.639      0.000     -0.000      0.000 3/21600
  480.0     64.756     64.756     43.848      0.000     -0.000      0.000 3/21600
  490.0     63.654     63.654     43.086      0.000     -0.000      0.000 3/21600
  500.0     62.591     62.591     42.352      0.000     -0.000      0.000 3/21600
  510.0     61.565     61.565     41.643      0.000     -0.000      0.000 3/21600
  520.0     60.573     60.573     40.959      0.000     -0.000      0.000 3/21600
  530.0     59.614     59.614     40.299      0.000     -0.000      0.000 3/21600
  540.0     58.686     58.686     39.660      0.000     -0.000      0.000 3/21600
  550.0     57.788     57.788     39.042      0.000     -0.000      0.000 3/21600
  560.0     56.918     56.918     38.444      0.000     -0.000      0.000 3/21600
  570.0     56.075     56.075     37.865      0.000     -0.000      0.000 3/21600
  580.0     55.258     55.258     37.304      0.000     -0.000      0.000 3/21600
  590.0     54.465     54.465     36.760      0.000     -0.000      0.000 3/21600
  600.0     53.695     53.695     36.232      0.000     -0.000      0.000 3/21600
  610.0     52.947     52.947     35.720      0.000     -0.000      0.000 3/21600
  620.0     52.221     52.221     35.223      0.000     -0.000      0.000 3/21600
  630.0     51.515     51.515     34.739      0.000     -0.000      0.000 3/21600
  640.0     50.828     50.828     34.270      0.000     -0.000      0.000 3/21600
  650.0     50.161     50.161     33.813      0.000     -0.000      0.000 3/21600
  660.0     49.511     49.511     33.369      0.000     -0.000      0.000 3/21600
  670.0     48.878     48.878     32.937      0.000     -0.000      0.000 3/21600
  680.0     48.261     48.261     32.516      0.000     -0.000      0.000 3/21600
  690.0     47.661     47.661     32.106      0.000     -0.000      0.000 3/21600
  700.0     47.075     47.075     31.707      0.000     -0.000      0.000 3/21600
  710.0     46.504     46.504     31.318      0.000     -0.000      0.000 3/21600
  720.0     45.947     45.947     30.938      0.000     -0.000      0.000 3/21600
  730.0     45.404     45.404     30.568      0.000     -0.000      0.000 3/21600
  740.0     44.874     44.874     30.207      0.000     -0.000      0.000 3/21600
  750.0     44.356     44.356     29.854      0.000     -0.000      0.000 3/21600
  760.0     43.850     43.850     29.510      0.000     -0.000      0.000 3/21600
  770.0     43.356     43.356     29.174      0.000     -0.000      0.000 3/21600
  780.0     42.874     42.874     28.846      0.000     -0.000      0.000 3/21600
  790.0     42.402     42.402     28.525      0.000     -0.000      0.000 3/21600
  800.0     41.941     41.941     28.212      0.000     -0.000      0.000 3/21600
  810.0     41.489     41.489     27.905      0.000     -0.000      0.000 3/21600
  820.0     41.048     41.048     27.606      0.000     -0.000      0.000 3/21600
  830.0     40.616     40.616     27.312      0.000     -0.000      0.000 3/21600
  840.0     40.194     40.194     27.025      0.000     -0.000      0.000 3/21600
  850.0     39.780     39.780     26.745      0.000     -0.000      0.000 3/21600
  860.0     39.375     39.375     26.470      0.000     -0.000      0.000 3/21600
  870.0     38.978     38.978     26.200      0.000     -0.000      0.000 3/21600
  880.0     38.589     38.589     25.937      0.000     -0.000      0.000 3/21600
  890.0     38.208     38.208     25.678      0.000     -0.000      0.000 3/21600
  900.0     37.834     37.834     25.425      0.000     -0.000      0.000 3/21600
  910.0     37.468     37.468     25.177      0.000     -0.000      0.000 3/21600
  920.0     37.110     37.110     24.934      0.000     -0.000      0.000 3/21600
  930.0     36.758     36.758     24.696      0.000     -0.000      0.000 3/21600
  940.0     36.412     36.412     24.462      0.000     -0.000      0.000 3/21600
  950.0     36.074     36.074     24.233      0.000     -0.000      0.000 3/21600
  960.0     35.741     35.741     24.008      0.000     -0.000      0.000 3/21600
  970.0     35.415     35.415     23.787      0.000     -0.000      0.000 3/21600
  980.0     35.095     35.095     23.570      0.000     -0.000      0.000 3/21600
  990.0     34.780     34.780     23.357      0.000     -0.000      0.000 3/21600
 1000.0     34.472     34.472     23.149      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-09 04:43:54]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

