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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [5 5 1]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
Number of symmetry operations in supercell: 600
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.269505289205699    0.000000000000000    0.000000000000000
  b   -1.634752644602850    2.831474638259723    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
   *1 Se  0.33333333333333  0.66666666666667  0.12482814326795  78.960
    2 Se  0.66666666666667  0.33333333333333  0.62482814326795  78.960
    3 Se  0.33333333333333  0.66666666666667  0.37517185673205  78.960
    4 Se  0.66666666666667  0.33333333333333  0.87517185673205  78.960
   *5 W   0.66666666666667  0.33333333333333  0.25000000000000 183.840
    6 W   0.33333333333333  0.66666666666667  0.75000000000000 183.840
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.269505289205699    0.000000000000000    0.000000000000000
  b   -1.634752644602850    2.831474638259723    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
   *1 Se  0.33333333333333  0.66666666666667  0.12482814326795  78.960 > 1
    2 Se  0.66666666666667  0.33333333333333  0.62482814326795  78.960 > 2
    3 Se  0.33333333333333  0.66666666666667  0.37517185673205  78.960 > 3
    4 Se  0.66666666666667  0.33333333333333  0.87517185673205  78.960 > 4
   *5 W   0.66666666666667  0.33333333333333  0.25000000000000 183.840 > 5
    6 W   0.33333333333333  0.66666666666667  0.75000000000000 183.840 > 6
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   16.347526446028496    0.000000000000000    0.000000000000000
  b   -8.173763223014250   14.157373191298614    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
   *1 Se  0.06666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    2 Se  0.26666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    3 Se  0.46666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    4 Se  0.66666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    5 Se  0.86666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    6 Se  0.06666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    7 Se  0.26666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    8 Se  0.46666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    9 Se  0.66666666666667  0.33333333333333  0.12482814326795  78.960 > 1
   10 Se  0.86666666666667  0.33333333333333  0.12482814326795  78.960 > 1
   11 Se  0.06666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   12 Se  0.26666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   13 Se  0.46666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   14 Se  0.66666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   15 Se  0.86666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   16 Se  0.06666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   17 Se  0.26666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   18 Se  0.46666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   19 Se  0.66666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   20 Se  0.86666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   21 Se  0.06666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   22 Se  0.26666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   23 Se  0.46666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   24 Se  0.66666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   25 Se  0.86666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   26 Se  0.13333333333333  0.06666666666667  0.62482814326795  78.960 > 2
   27 Se  0.33333333333333  0.06666666666667  0.62482814326795  78.960 > 2
   28 Se  0.53333333333333  0.06666666666667  0.62482814326795  78.960 > 2
   29 Se  0.73333333333333  0.06666666666667  0.62482814326795  78.960 > 2
   30 Se  0.93333333333333  0.06666666666667  0.62482814326795  78.960 > 2
   31 Se  0.13333333333333  0.26666666666667  0.62482814326795  78.960 > 2
   32 Se  0.33333333333333  0.26666666666667  0.62482814326795  78.960 > 2
   33 Se  0.53333333333333  0.26666666666667  0.62482814326795  78.960 > 2
   34 Se  0.73333333333333  0.26666666666667  0.62482814326795  78.960 > 2
   35 Se  0.93333333333333  0.26666666666667  0.62482814326795  78.960 > 2
   36 Se  0.13333333333333  0.46666666666667  0.62482814326795  78.960 > 2
   37 Se  0.33333333333333  0.46666666666667  0.62482814326795  78.960 > 2
   38 Se  0.53333333333333  0.46666666666667  0.62482814326795  78.960 > 2
   39 Se  0.73333333333333  0.46666666666667  0.62482814326795  78.960 > 2
   40 Se  0.93333333333333  0.46666666666667  0.62482814326795  78.960 > 2
   41 Se  0.13333333333333  0.66666666666667  0.62482814326795  78.960 > 2
   42 Se  0.33333333333333  0.66666666666667  0.62482814326795  78.960 > 2
   43 Se  0.53333333333333  0.66666666666667  0.62482814326795  78.960 > 2
   44 Se  0.73333333333333  0.66666666666667  0.62482814326795  78.960 > 2
   45 Se  0.93333333333333  0.66666666666667  0.62482814326795  78.960 > 2
   46 Se  0.13333333333333  0.86666666666667  0.62482814326795  78.960 > 2
   47 Se  0.33333333333333  0.86666666666667  0.62482814326795  78.960 > 2
   48 Se  0.53333333333333  0.86666666666667  0.62482814326795  78.960 > 2
   49 Se  0.73333333333333  0.86666666666667  0.62482814326795  78.960 > 2
   50 Se  0.93333333333333  0.86666666666667  0.62482814326795  78.960 > 2
   51 Se  0.06666666666667  0.13333333333333  0.37517185673205  78.960 > 3
   52 Se  0.26666666666667  0.13333333333333  0.37517185673205  78.960 > 3
   53 Se  0.46666666666667  0.13333333333333  0.37517185673205  78.960 > 3
   54 Se  0.66666666666667  0.13333333333333  0.37517185673205  78.960 > 3
   55 Se  0.86666666666667  0.13333333333333  0.37517185673205  78.960 > 3
   56 Se  0.06666666666667  0.33333333333333  0.37517185673205  78.960 > 3
   57 Se  0.26666666666667  0.33333333333333  0.37517185673205  78.960 > 3
   58 Se  0.46666666666667  0.33333333333333  0.37517185673205  78.960 > 3
   59 Se  0.66666666666667  0.33333333333333  0.37517185673205  78.960 > 3
   60 Se  0.86666666666667  0.33333333333333  0.37517185673205  78.960 > 3
   61 Se  0.06666666666667  0.53333333333333  0.37517185673205  78.960 > 3
   62 Se  0.26666666666667  0.53333333333333  0.37517185673205  78.960 > 3
   63 Se  0.46666666666667  0.53333333333333  0.37517185673205  78.960 > 3
   64 Se  0.66666666666667  0.53333333333333  0.37517185673205  78.960 > 3
   65 Se  0.86666666666667  0.53333333333333  0.37517185673205  78.960 > 3
   66 Se  0.06666666666667  0.73333333333333  0.37517185673205  78.960 > 3
   67 Se  0.26666666666667  0.73333333333333  0.37517185673205  78.960 > 3
   68 Se  0.46666666666667  0.73333333333333  0.37517185673205  78.960 > 3
   69 Se  0.66666666666667  0.73333333333333  0.37517185673205  78.960 > 3
   70 Se  0.86666666666667  0.73333333333333  0.37517185673205  78.960 > 3
   71 Se  0.06666666666667  0.93333333333333  0.37517185673205  78.960 > 3
   72 Se  0.26666666666667  0.93333333333333  0.37517185673205  78.960 > 3
   73 Se  0.46666666666667  0.93333333333333  0.37517185673205  78.960 > 3
   74 Se  0.66666666666667  0.93333333333333  0.37517185673205  78.960 > 3
   75 Se  0.86666666666667  0.93333333333333  0.37517185673205  78.960 > 3
   76 Se  0.13333333333333  0.06666666666667  0.87517185673205  78.960 > 4
   77 Se  0.33333333333333  0.06666666666667  0.87517185673205  78.960 > 4
   78 Se  0.53333333333333  0.06666666666667  0.87517185673205  78.960 > 4
   79 Se  0.73333333333333  0.06666666666667  0.87517185673205  78.960 > 4
   80 Se  0.93333333333333  0.06666666666667  0.87517185673205  78.960 > 4
   81 Se  0.13333333333333  0.26666666666667  0.87517185673205  78.960 > 4
   82 Se  0.33333333333333  0.26666666666667  0.87517185673205  78.960 > 4
   83 Se  0.53333333333333  0.26666666666667  0.87517185673205  78.960 > 4
   84 Se  0.73333333333333  0.26666666666667  0.87517185673205  78.960 > 4
   85 Se  0.93333333333333  0.26666666666667  0.87517185673205  78.960 > 4
   86 Se  0.13333333333333  0.46666666666667  0.87517185673205  78.960 > 4
   87 Se  0.33333333333333  0.46666666666667  0.87517185673205  78.960 > 4
   88 Se  0.53333333333333  0.46666666666667  0.87517185673205  78.960 > 4
   89 Se  0.73333333333333  0.46666666666667  0.87517185673205  78.960 > 4
   90 Se  0.93333333333333  0.46666666666667  0.87517185673205  78.960 > 4
   91 Se  0.13333333333333  0.66666666666667  0.87517185673205  78.960 > 4
   92 Se  0.33333333333333  0.66666666666667  0.87517185673205  78.960 > 4
   93 Se  0.53333333333333  0.66666666666667  0.87517185673205  78.960 > 4
   94 Se  0.73333333333333  0.66666666666667  0.87517185673205  78.960 > 4
   95 Se  0.93333333333333  0.66666666666667  0.87517185673205  78.960 > 4
   96 Se  0.13333333333333  0.86666666666667  0.87517185673205  78.960 > 4
   97 Se  0.33333333333333  0.86666666666667  0.87517185673205  78.960 > 4
   98 Se  0.53333333333333  0.86666666666667  0.87517185673205  78.960 > 4
   99 Se  0.73333333333333  0.86666666666667  0.87517185673205  78.960 > 4
  100 Se  0.93333333333333  0.86666666666667  0.87517185673205  78.960 > 4
 *101 W   0.13333333333333  0.06666666666667  0.25000000000000 183.840 > 5
  102 W   0.33333333333333  0.06666666666667  0.25000000000000 183.840 > 5
  103 W   0.53333333333333  0.06666666666667  0.25000000000000 183.840 > 5
  104 W   0.73333333333333  0.06666666666667  0.25000000000000 183.840 > 5
  105 W   0.93333333333333  0.06666666666667  0.25000000000000 183.840 > 5
  106 W   0.13333333333333  0.26666666666667  0.25000000000000 183.840 > 5
  107 W   0.33333333333333  0.26666666666667  0.25000000000000 183.840 > 5
  108 W   0.53333333333333  0.26666666666667  0.25000000000000 183.840 > 5
  109 W   0.73333333333333  0.26666666666667  0.25000000000000 183.840 > 5
  110 W   0.93333333333333  0.26666666666667  0.25000000000000 183.840 > 5
  111 W   0.13333333333333  0.46666666666667  0.25000000000000 183.840 > 5
  112 W   0.33333333333333  0.46666666666667  0.25000000000000 183.840 > 5
  113 W   0.53333333333333  0.46666666666667  0.25000000000000 183.840 > 5
  114 W   0.73333333333333  0.46666666666667  0.25000000000000 183.840 > 5
  115 W   0.93333333333333  0.46666666666667  0.25000000000000 183.840 > 5
  116 W   0.13333333333333  0.66666666666667  0.25000000000000 183.840 > 5
  117 W   0.33333333333333  0.66666666666667  0.25000000000000 183.840 > 5
  118 W   0.53333333333333  0.66666666666667  0.25000000000000 183.840 > 5
  119 W   0.73333333333333  0.66666666666667  0.25000000000000 183.840 > 5
  120 W   0.93333333333333  0.66666666666667  0.25000000000000 183.840 > 5
  121 W   0.13333333333333  0.86666666666667  0.25000000000000 183.840 > 5
  122 W   0.33333333333333  0.86666666666667  0.25000000000000 183.840 > 5
  123 W   0.53333333333333  0.86666666666667  0.25000000000000 183.840 > 5
  124 W   0.73333333333333  0.86666666666667  0.25000000000000 183.840 > 5
  125 W   0.93333333333333  0.86666666666667  0.25000000000000 183.840 > 5
  126 W   0.06666666666667  0.13333333333333  0.75000000000000 183.840 > 6
  127 W   0.26666666666667  0.13333333333333  0.75000000000000 183.840 > 6
  128 W   0.46666666666667  0.13333333333333  0.75000000000000 183.840 > 6
  129 W   0.66666666666667  0.13333333333333  0.75000000000000 183.840 > 6
  130 W   0.86666666666667  0.13333333333333  0.75000000000000 183.840 > 6
  131 W   0.06666666666667  0.33333333333333  0.75000000000000 183.840 > 6
  132 W   0.26666666666667  0.33333333333333  0.75000000000000 183.840 > 6
  133 W   0.46666666666667  0.33333333333333  0.75000000000000 183.840 > 6
  134 W   0.66666666666667  0.33333333333333  0.75000000000000 183.840 > 6
  135 W   0.86666666666667  0.33333333333333  0.75000000000000 183.840 > 6
  136 W   0.06666666666667  0.53333333333333  0.75000000000000 183.840 > 6
  137 W   0.26666666666667  0.53333333333333  0.75000000000000 183.840 > 6
  138 W   0.46666666666667  0.53333333333333  0.75000000000000 183.840 > 6
  139 W   0.66666666666667  0.53333333333333  0.75000000000000 183.840 > 6
  140 W   0.86666666666667  0.53333333333333  0.75000000000000 183.840 > 6
  141 W   0.06666666666667  0.73333333333333  0.75000000000000 183.840 > 6
  142 W   0.26666666666667  0.73333333333333  0.75000000000000 183.840 > 6
  143 W   0.46666666666667  0.73333333333333  0.75000000000000 183.840 > 6
  144 W   0.66666666666667  0.73333333333333  0.75000000000000 183.840 > 6
  145 W   0.86666666666667  0.73333333333333  0.75000000000000 183.840 > 6
  146 W   0.06666666666667  0.93333333333333  0.75000000000000 183.840 > 6
  147 W   0.26666666666667  0.93333333333333  0.75000000000000 183.840 > 6
  148 W   0.46666666666667  0.93333333333333  0.75000000000000 183.840 > 6
  149 W   0.66666666666667  0.93333333333333  0.75000000000000 183.840 > 6
  150 W   0.86666666666667  0.93333333333333  0.75000000000000 183.840 > 6
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           14.4530973   -0.0000000    0.0000000
           -0.0000000   14.4530973    0.0000000
            0.0000000    0.0000000    7.3460639
-------------------------- Born effective charges --------------------------
    1 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    2 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    3 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    4 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    5 W    -1.2844985    0.0000000    0.0000000
            0.0000000   -1.2844985    0.0000000
            0.0000000    0.0000000   -0.4074682
    6 W    -1.2844985    0.0000000    0.0000000
            0.0000000   -1.2844985    0.0000000
            0.0000000    0.0000000   -0.4074682
----------------------------------------------------------------------------
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: 450/450
Permutation basis: 4689/4689
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 226
Number of blocks in projector: 182
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: 112
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 64
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 50
Use standard eigh solver.
Tree of FC basis block matrices:
- (226, 222), data: False
|-- (50, 48), data: True
|-- (64, 64), data: True
|-- (112, 110), data: True
-----
Solver_atoms: 1 -- 150 / 150
Time (Solver_compr_matrix_reshape): 0.013
Solver_block: 80 / 80
 - Time: 0.188
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.202
--------------------------------- Symfc end --------------------------------
Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) 
Permutation basis: 450/450
Permutation basis: 4689/4689
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 226
Number of blocks in projector: 182
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: 112
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 64
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 50
Use standard eigh solver.
Tree of FC basis block matrices:
- (226, 222), data: False
|-- (50, 48), data: True
|-- (64, 64), data: True
|-- (112, 110), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-07 18:59:09]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [5 5 1]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.269505289205699    0.000000000000000    0.000000000000000
  b   -1.634752644602850    2.831474638259723    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
    1 Se  0.33333333333333  0.66666666666667  0.12482814326795  78.960
    2 Se  0.66666666666667  0.33333333333333  0.62482814326795  78.960
    3 Se  0.33333333333333  0.66666666666667  0.37517185673205  78.960
    4 Se  0.66666666666667  0.33333333333333  0.87517185673205  78.960
    5 W   0.66666666666667  0.33333333333333  0.25000000000000 183.840
    6 W   0.33333333333333  0.66666666666667  0.75000000000000 183.840
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.347526446028496    0.000000000000000    0.000000000000000
  b   -8.173763223014250   14.157373191298614    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
    1 Se  0.06666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    2 Se  0.26666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    3 Se  0.46666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    4 Se  0.66666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    5 Se  0.86666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    6 Se  0.06666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    7 Se  0.26666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    8 Se  0.46666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    9 Se  0.66666666666667  0.33333333333333  0.12482814326795  78.960 > 1
   10 Se  0.86666666666667  0.33333333333333  0.12482814326795  78.960 > 1
   11 Se  0.06666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   12 Se  0.26666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   13 Se  0.46666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   14 Se  0.66666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   15 Se  0.86666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   16 Se  0.06666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   17 Se  0.26666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   18 Se  0.46666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   19 Se  0.66666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   20 Se  0.86666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   21 Se  0.06666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   22 Se  0.26666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   23 Se  0.46666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   24 Se  0.66666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   25 Se  0.86666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   26 Se  0.13333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   27 Se  0.33333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   28 Se  0.53333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   29 Se  0.73333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   30 Se  0.93333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   31 Se  0.13333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   32 Se  0.33333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   33 Se  0.53333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   34 Se  0.73333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   35 Se  0.93333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   36 Se  0.13333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   37 Se  0.33333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   38 Se  0.53333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   39 Se  0.73333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   40 Se  0.93333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   41 Se  0.13333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   42 Se  0.33333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   43 Se  0.53333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   44 Se  0.73333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   45 Se  0.93333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   46 Se  0.13333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   47 Se  0.33333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   48 Se  0.53333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   49 Se  0.73333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   50 Se  0.93333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   51 Se  0.06666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   52 Se  0.26666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   53 Se  0.46666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   54 Se  0.66666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   55 Se  0.86666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   56 Se  0.06666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   57 Se  0.26666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   58 Se  0.46666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   59 Se  0.66666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   60 Se  0.86666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   61 Se  0.06666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   62 Se  0.26666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   63 Se  0.46666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   64 Se  0.66666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   65 Se  0.86666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   66 Se  0.06666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   67 Se  0.26666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   68 Se  0.46666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   69 Se  0.66666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   70 Se  0.86666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   71 Se  0.06666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   72 Se  0.26666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   73 Se  0.46666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   74 Se  0.66666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   75 Se  0.86666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   76 Se  0.13333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   77 Se  0.33333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   78 Se  0.53333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   79 Se  0.73333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   80 Se  0.93333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   81 Se  0.13333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   82 Se  0.33333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   83 Se  0.53333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   84 Se  0.73333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   85 Se  0.93333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   86 Se  0.13333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   87 Se  0.33333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   88 Se  0.53333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   89 Se  0.73333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   90 Se  0.93333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   91 Se  0.13333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   92 Se  0.33333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   93 Se  0.53333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   94 Se  0.73333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   95 Se  0.93333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   96 Se  0.13333333333333  0.86666666666667  0.87517185673205  78.960 > 76
   97 Se  0.33333333333333  0.86666666666667  0.87517185673205  78.960 > 76
   98 Se  0.53333333333333  0.86666666666667  0.87517185673205  78.960 > 76
   99 Se  0.73333333333333  0.86666666666667  0.87517185673205  78.960 > 76
  100 Se  0.93333333333333  0.86666666666667  0.87517185673205  78.960 > 76
  101 W   0.13333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  102 W   0.33333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  103 W   0.53333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  104 W   0.73333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  105 W   0.93333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  106 W   0.13333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  107 W   0.33333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  108 W   0.53333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  109 W   0.73333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  110 W   0.93333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  111 W   0.13333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  112 W   0.33333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  113 W   0.53333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  114 W   0.73333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  115 W   0.93333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  116 W   0.13333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  117 W   0.33333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  118 W   0.53333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  119 W   0.73333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  120 W   0.93333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  121 W   0.13333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  122 W   0.33333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  123 W   0.53333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  124 W   0.73333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  125 W   0.93333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  126 W   0.06666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  127 W   0.26666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  128 W   0.46666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  129 W   0.66666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  130 W   0.86666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  131 W   0.06666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  132 W   0.26666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  133 W   0.46666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  134 W   0.66666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  135 W   0.86666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  136 W   0.06666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  137 W   0.26666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  138 W   0.46666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  139 W   0.66666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  140 W   0.86666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  141 W   0.06666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  142 W   0.26666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  143 W   0.46666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  144 W   0.66666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  145 W   0.86666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  146 W   0.06666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  147 W   0.26666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  148 W   0.46666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  149 W   0.66666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  150 W   0.86666666666667  0.93333333333333  0.75000000000000 183.840 > 126
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
           14.4530973   -0.0000000    0.0000000
           -0.0000000   14.4530973    0.0000000
            0.0000000    0.0000000    7.3460639
-------------------------- Born effective charges --------------------------
    1 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    2 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    3 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    4 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    5 W    -1.2844985    0.0000000    0.0000000
            0.0000000   -1.2844985    0.0000000
            0.0000000    0.0000000   -0.4074682
    6 W    -1.2844985    0.0000000    0.0000000
            0.0000000   -1.2844985    0.0000000
            0.0000000    0.0000000   -0.4074682
----------------------------------------------------------------------------
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[ 101, 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.00000000 (xxz) 0.00000000 (xxz) 0.00000000 (xzx)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-07 18:59:16]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 18:59:16]-------------------------
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): [5 5 1]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.269505289205699    0.000000000000000    0.000000000000000
  b   -1.634752644602850    2.831474638259723    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
    1 Se  0.33333333333333  0.66666666666667  0.12482814326795  78.960
    2 Se  0.66666666666667  0.33333333333333  0.62482814326795  78.960
    3 Se  0.33333333333333  0.66666666666667  0.37517185673205  78.960
    4 Se  0.66666666666667  0.33333333333333  0.87517185673205  78.960
    5 W   0.66666666666667  0.33333333333333  0.25000000000000 183.840
    6 W   0.33333333333333  0.66666666666667  0.75000000000000 183.840
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   16.347526446028496    0.000000000000000    0.000000000000000
  b   -8.173763223014250   14.157373191298614    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.404458419999999
Atomic positions (fractional):
    1 Se  0.06666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    2 Se  0.26666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    3 Se  0.46666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    4 Se  0.66666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    5 Se  0.86666666666667  0.13333333333333  0.12482814326795  78.960 > 1
    6 Se  0.06666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    7 Se  0.26666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    8 Se  0.46666666666667  0.33333333333333  0.12482814326795  78.960 > 1
    9 Se  0.66666666666667  0.33333333333333  0.12482814326795  78.960 > 1
   10 Se  0.86666666666667  0.33333333333333  0.12482814326795  78.960 > 1
   11 Se  0.06666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   12 Se  0.26666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   13 Se  0.46666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   14 Se  0.66666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   15 Se  0.86666666666667  0.53333333333333  0.12482814326795  78.960 > 1
   16 Se  0.06666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   17 Se  0.26666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   18 Se  0.46666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   19 Se  0.66666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   20 Se  0.86666666666667  0.73333333333333  0.12482814326795  78.960 > 1
   21 Se  0.06666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   22 Se  0.26666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   23 Se  0.46666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   24 Se  0.66666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   25 Se  0.86666666666667  0.93333333333333  0.12482814326795  78.960 > 1
   26 Se  0.13333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   27 Se  0.33333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   28 Se  0.53333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   29 Se  0.73333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   30 Se  0.93333333333333  0.06666666666667  0.62482814326795  78.960 > 26
   31 Se  0.13333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   32 Se  0.33333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   33 Se  0.53333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   34 Se  0.73333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   35 Se  0.93333333333333  0.26666666666667  0.62482814326795  78.960 > 26
   36 Se  0.13333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   37 Se  0.33333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   38 Se  0.53333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   39 Se  0.73333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   40 Se  0.93333333333333  0.46666666666667  0.62482814326795  78.960 > 26
   41 Se  0.13333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   42 Se  0.33333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   43 Se  0.53333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   44 Se  0.73333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   45 Se  0.93333333333333  0.66666666666667  0.62482814326795  78.960 > 26
   46 Se  0.13333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   47 Se  0.33333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   48 Se  0.53333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   49 Se  0.73333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   50 Se  0.93333333333333  0.86666666666667  0.62482814326795  78.960 > 26
   51 Se  0.06666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   52 Se  0.26666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   53 Se  0.46666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   54 Se  0.66666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   55 Se  0.86666666666667  0.13333333333333  0.37517185673205  78.960 > 51
   56 Se  0.06666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   57 Se  0.26666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   58 Se  0.46666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   59 Se  0.66666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   60 Se  0.86666666666667  0.33333333333333  0.37517185673205  78.960 > 51
   61 Se  0.06666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   62 Se  0.26666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   63 Se  0.46666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   64 Se  0.66666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   65 Se  0.86666666666667  0.53333333333333  0.37517185673205  78.960 > 51
   66 Se  0.06666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   67 Se  0.26666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   68 Se  0.46666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   69 Se  0.66666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   70 Se  0.86666666666667  0.73333333333333  0.37517185673205  78.960 > 51
   71 Se  0.06666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   72 Se  0.26666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   73 Se  0.46666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   74 Se  0.66666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   75 Se  0.86666666666667  0.93333333333333  0.37517185673205  78.960 > 51
   76 Se  0.13333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   77 Se  0.33333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   78 Se  0.53333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   79 Se  0.73333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   80 Se  0.93333333333333  0.06666666666667  0.87517185673205  78.960 > 76
   81 Se  0.13333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   82 Se  0.33333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   83 Se  0.53333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   84 Se  0.73333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   85 Se  0.93333333333333  0.26666666666667  0.87517185673205  78.960 > 76
   86 Se  0.13333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   87 Se  0.33333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   88 Se  0.53333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   89 Se  0.73333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   90 Se  0.93333333333333  0.46666666666667  0.87517185673205  78.960 > 76
   91 Se  0.13333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   92 Se  0.33333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   93 Se  0.53333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   94 Se  0.73333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   95 Se  0.93333333333333  0.66666666666667  0.87517185673205  78.960 > 76
   96 Se  0.13333333333333  0.86666666666667  0.87517185673205  78.960 > 76
   97 Se  0.33333333333333  0.86666666666667  0.87517185673205  78.960 > 76
   98 Se  0.53333333333333  0.86666666666667  0.87517185673205  78.960 > 76
   99 Se  0.73333333333333  0.86666666666667  0.87517185673205  78.960 > 76
  100 Se  0.93333333333333  0.86666666666667  0.87517185673205  78.960 > 76
  101 W   0.13333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  102 W   0.33333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  103 W   0.53333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  104 W   0.73333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  105 W   0.93333333333333  0.06666666666667  0.25000000000000 183.840 > 101
  106 W   0.13333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  107 W   0.33333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  108 W   0.53333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  109 W   0.73333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  110 W   0.93333333333333  0.26666666666667  0.25000000000000 183.840 > 101
  111 W   0.13333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  112 W   0.33333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  113 W   0.53333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  114 W   0.73333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  115 W   0.93333333333333  0.46666666666667  0.25000000000000 183.840 > 101
  116 W   0.13333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  117 W   0.33333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  118 W   0.53333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  119 W   0.73333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  120 W   0.93333333333333  0.66666666666667  0.25000000000000 183.840 > 101
  121 W   0.13333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  122 W   0.33333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  123 W   0.53333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  124 W   0.73333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  125 W   0.93333333333333  0.86666666666667  0.25000000000000 183.840 > 101
  126 W   0.06666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  127 W   0.26666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  128 W   0.46666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  129 W   0.66666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  130 W   0.86666666666667  0.13333333333333  0.75000000000000 183.840 > 126
  131 W   0.06666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  132 W   0.26666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  133 W   0.46666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  134 W   0.66666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  135 W   0.86666666666667  0.33333333333333  0.75000000000000 183.840 > 126
  136 W   0.06666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  137 W   0.26666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  138 W   0.46666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  139 W   0.66666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  140 W   0.86666666666667  0.53333333333333  0.75000000000000 183.840 > 126
  141 W   0.06666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  142 W   0.26666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  143 W   0.46666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  144 W   0.66666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  145 W   0.86666666666667  0.73333333333333  0.75000000000000 183.840 > 126
  146 W   0.06666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  147 W   0.26666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  148 W   0.46666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  149 W   0.66666666666667  0.93333333333333  0.75000000000000 183.840 > 126
  150 W   0.86666666666667  0.93333333333333  0.75000000000000 183.840 > 126
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
           14.4530973   -0.0000000    0.0000000
           -0.0000000   14.4530973    0.0000000
            0.0000000    0.0000000    7.3460639
-------------------------- Born effective charges --------------------------
    1 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    2 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    3 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    4 Se    0.6422492   -0.0000000    0.0000000
            0.0000000    0.6422492    0.0000000
            0.0000000    0.0000000    0.2037341
    5 W    -1.2844985    0.0000000    0.0000000
            0.0000000   -1.2844985    0.0000000
            0.0000000    0.0000000   -0.4074682
    6 W    -1.2844985    0.0000000    0.0000000
            0.0000000   -1.2844985    0.0000000
            0.0000000    0.0000000   -0.4074682
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (xxz) 0.00000000 (xxz) 0.00000000 (xzx)
Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 18 18 4 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.83, Number of G-points: 305, Lambda: 0.30
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/111) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 111
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   0.503   (   0.000    0.000    0.000)    0.000
   0.503   (   0.000    0.000    0.000)    0.000
   0.651   (   0.000    0.000    0.000)    0.000
   5.181   (   0.000    0.000    0.000)    0.000
   5.181   (   0.000    0.000    0.000)    0.000
   5.186   (   0.000    0.000    0.000)    0.000
   5.186   (   0.000    0.000    0.000)    0.000
   7.324   (   0.000    0.000    0.000)    0.000
   7.324   (   0.000    0.000    0.000)    0.000
   7.342   (   0.000    0.000    0.000)    0.000
   7.342   (   0.000    0.000    0.000)    0.000
   7.410   (   0.000    0.000    0.000)    0.000
   7.445   (   0.000    0.000    0.000)    0.000
   9.182   (   0.000    0.000    0.000)    0.000
   9.207   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/111) =======================
q-point: ( 0.06  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.224   (  12.681    7.321    0.000)   14.642
   0.583   (  25.261   14.584    0.000)   29.169
   0.633   (  -0.357   -0.206    0.000)    0.412
   0.772   (  19.253   11.116    0.000)   22.232
   0.892   (  38.697   22.342    0.000)   44.684
   1.017   (  33.311   19.232    0.000)   38.465
   5.201   (   1.703    0.983    0.000)    1.966
   5.205   (   1.642    0.948    0.000)    1.896
   5.254   (   6.218    3.590    0.000)    7.180
   5.269   (   7.026    4.057    0.000)    8.113
   7.220   (  -7.067   -4.080    0.000)    8.160
   7.251   (  -7.730   -4.463    0.000)    8.926
   7.324   (   0.020    0.012    0.000)    0.023
   7.341   (  -0.127   -0.073    0.000)    0.146
   7.513   (   5.328    3.076    0.000)    6.152
   7.554   (   5.262    3.038    0.000)    6.076
   9.139   (  -3.827   -2.209    0.000)    4.419
   9.159   (  -4.109   -2.372    0.000)    4.744
======================= Grid point 2 (3/111) =======================
q-point: ( 0.11  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.592   (  19.624   11.330    0.000)   22.660
   0.723   (   9.822    5.671    0.000)   11.341
   1.135   (  23.288   13.445    0.000)   26.891
   1.245   (  21.276   12.284    0.000)   24.568
   1.737   (  35.538   20.518    0.000)   41.035
   1.794   (  33.906   19.576    0.000)   39.151
   5.256   (   3.091    1.785    0.000)    3.570
   5.258   (   2.985    1.724    0.000)    3.447
   5.447   (  10.436    6.025    0.000)   12.050
   5.483   (  11.334    6.544    0.000)   13.088
   7.047   (  -8.004   -4.621    0.000)    9.243
   7.064   (  -8.562   -4.943    0.000)    9.887
   7.323   (  -0.139   -0.080    0.000)    0.160
   7.336   (  -0.326   -0.188    0.000)    0.376
   7.627   (   4.693    2.709    0.000)    5.419
   7.665   (   4.508    2.602    0.000)    5.205
   9.011   (  -7.330   -4.232    0.000)    8.464
   9.025   (  -7.641   -4.412    0.000)    8.823
======================= Grid point 3 (4/111) =======================
q-point: ( 0.17  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.076   (  20.392   11.773    0.000)   23.547
   1.095   (  24.165   13.952    0.000)   27.904
   1.630   (  20.282   11.710    0.000)   23.420
   1.708   (  19.243   11.110    0.000)   22.220
   2.485   (  30.056   17.353    0.000)   34.706
   2.522   (  29.949   17.291    0.000)   34.582
   5.334   (   3.744    2.162    0.000)    4.323
   5.338   (   3.984    2.300    0.000)    4.601
   5.707   (  12.085    6.977    0.000)   13.954
   5.761   (  12.788    7.383    0.000)   14.766
   6.865   (  -7.907   -4.565    0.000)    9.130
   6.873   (  -7.946   -4.588    0.000)    9.176
   7.316   (  -0.540   -0.312    0.000)    0.623
   7.325   (  -0.611   -0.353    0.000)    0.705
   7.725   (   3.971    2.293    0.000)    4.585
   7.760   (   3.923    2.265    0.000)    4.529
   8.811   ( -10.214   -5.897    0.000)   11.795
   8.819   ( -10.344   -5.972    0.000)   11.944
======================= Grid point 4 (5/111) =======================
q-point: ( 0.22  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.594   (  24.471   14.129    0.000)   28.257
   1.667   (  25.888   14.946    0.000)   29.892
   2.049   (  16.622    9.597    0.000)   19.194
   2.108   (  15.985    9.229    0.000)   18.458
   3.085   (  22.629   13.065    0.000)   26.130
   3.130   (  23.243   13.419    0.000)   26.839
   5.422   (   3.897    2.250    0.000)    4.500
   5.432   (   4.204    2.427    0.000)    4.854
   5.978   (  11.462    6.618    0.000)   13.235
   6.046   (  12.028    6.944    0.000)   13.889
   6.701   (  -6.238   -3.601    0.000)    7.203
   6.718   (  -5.590   -3.227    0.000)    6.454
   7.298   (  -1.033   -0.597    0.000)    1.193
   7.307   (  -0.970   -0.560    0.000)    1.120
   7.807   (   3.246    1.874    0.000)    3.748
   7.842   (   3.236    1.868    0.000)    3.737
   8.556   ( -12.090   -6.980    0.000)   13.960
   8.564   ( -11.996   -6.926    0.000)   13.851
======================= Grid point 5 (6/111) =======================
q-point: ( 0.28  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.154   (  24.454   14.118    0.000)   28.237
   2.249   (  25.038   14.456    0.000)   28.911
   2.381   (  12.725    7.347    0.000)   14.693
   2.430   (  12.381    7.148    0.000)   14.297
   3.500   (  13.835    7.988    0.000)   15.975
   3.558   (  14.266    8.236    0.000)   16.473
   5.507   (   3.513    2.028    0.000)    4.057
   5.524   (   3.799    2.193    0.000)    4.387
   6.211   (   8.815    5.089    0.000)   10.178
   6.291   (   9.352    5.399    0.000)   10.799
   6.598   (  -2.518   -1.454    0.000)    2.907
   6.630   (  -2.019   -1.166    0.000)    2.332
   7.270   (  -1.416   -0.817    0.000)    1.634
   7.282   (  -1.280   -0.739    0.000)    1.478
   7.871   (   2.336    1.349    0.000)    2.697
   7.904   (   2.234    1.290    0.000)    2.580
   8.276   ( -12.267   -7.082    0.000)   14.165
   8.287   ( -12.186   -7.035    0.000)   14.071
======================= Grid point 6 (7/111) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.627   (   8.962    5.174    0.000)   10.349
   2.670   (   8.868    5.120    0.000)   10.240
   2.683   (  21.795   12.583    0.000)   25.167
   2.786   (  22.100   12.759    0.000)   25.519
   3.713   (   5.288    3.053    0.000)    6.106
   3.775   (   5.180    2.991    0.000)    5.981
   5.578   (   2.764    1.596    0.000)    3.192
   5.601   (   2.968    1.714    0.000)    3.427
   6.363   (   4.317    2.492    0.000)    4.984
   6.455   (   4.826    2.786    0.000)    5.573
   6.591   (   1.710    0.987    0.000)    1.974
   6.628   (   1.711    0.988    0.000)    1.975
   7.236   (  -1.542   -0.890    0.000)    1.780
   7.251   (  -1.355   -0.782    0.000)    1.564
   7.912   (   1.256    0.725    0.000)    1.451
   7.941   (   1.000    0.577    0.000)    1.154
   8.018   ( -10.043   -5.799    0.000)   11.597
   8.029   ( -10.128   -5.848    0.000)   11.695
======================= Grid point 7 (8/111) =======================
q-point: ( 0.39  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.791   (   5.583    3.223    0.000)    6.447
   2.834   (   5.640    3.256    0.000)    6.512
   3.125   (  16.929    9.774    0.000)   19.548
   3.235   (  17.133    9.892    0.000)   19.783
   3.762   (  -0.269   -0.155    0.000)    0.310
   3.817   (  -0.715   -0.413    0.000)    0.825
   5.631   (   1.849    1.068    0.000)    2.135
   5.657   (   1.957    1.130    0.000)    2.260
   6.403   (  -0.578   -0.334    0.000)    0.668
   6.504   (  -0.328   -0.189    0.000)    0.378
   6.660   (   3.927    2.267    0.000)    4.535
   6.692   (   3.529    2.037    0.000)    4.075
   7.203   (  -1.308   -0.755    0.000)    1.510
   7.222   (  -1.132   -0.654    0.000)    1.308
   7.837   (  -5.827   -3.364    0.000)    6.729
   7.845   (  -5.939   -3.429    0.000)    6.857
   7.929   (   0.339    0.196    0.000)    0.391
   7.951   (  -0.061   -0.035    0.000)    0.071
======================= Grid point 8 (9/111) =======================
q-point: ( 0.44  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.883   (   2.646    1.528    0.000)    3.056
   2.928   (   2.720    1.571    0.000)    3.141
   3.431   (   9.591    5.538    0.000)   11.075
   3.544   (   9.716    5.609    0.000)   11.219
   3.734   (  -1.592   -0.919    0.000)    1.839
   3.779   (  -1.946   -1.124    0.000)    2.247
   5.662   (   0.912    0.527    0.000)    1.054
   5.690   (   0.951    0.549    0.000)    1.098
   6.361   (  -2.343   -1.353    0.000)    2.706
   6.464   (  -2.405   -1.388    0.000)    2.776
   6.745   (   3.039    1.754    0.000)    3.509
   6.767   (   2.625    1.515    0.000)    3.031
   7.179   (  -0.746   -0.431    0.000)    0.862
   7.202   (  -0.645   -0.372    0.000)    0.744
   7.750   (  -2.091   -1.207    0.000)    2.414
   7.757   (  -2.121   -1.225    0.000)    2.450
   7.931   (  -0.025   -0.014    0.000)    0.029
   7.943   (  -0.467   -0.269    0.000)    0.539
======================= Grid point 9 (10/111) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.913   (   0.000    0.000    0.000)    0.000
   2.959   (   0.000    0.000    0.000)    0.000
   3.542   (   0.000    0.000    0.000)    0.000
   3.657   (   0.000    0.000    0.000)    0.000
   3.713   (   0.000    0.000    0.000)    0.000
   3.754   (   0.000    0.000    0.000)    0.000
   5.672   (   0.000    0.000    0.000)    0.000
   5.701   (   0.000    0.000    0.000)    0.000
   6.329   (   0.000    0.000    0.000)    0.000
   6.431   (   0.000    0.000    0.000)    0.000
   6.781   (  -0.000   -0.000    0.000)    0.000
   6.798   (   0.000    0.000    0.000)    0.000
   7.170   (   0.000    0.000    0.000)    0.000
   7.194   (   0.000    0.000    0.000)    0.000
   7.729   (   0.000    0.000    0.000)    0.000
   7.735   (   0.000    0.000    0.000)    0.000
   7.931   (   0.000    0.000    0.000)    0.000
   7.936   (   0.000    0.000    0.000)    0.000
======================= Grid point 20 (11/111) =======================
q-point: ( 0.06  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.477   (  10.274   17.795    0.000)   20.548
   0.673   (   3.593    6.223    0.000)    7.186
   1.007   (  14.446   25.021    0.000)   28.891
   1.129   (  13.005   22.526    0.000)   26.010
   1.513   (  20.849   36.112    0.000)   41.699
   1.574   (  19.351   33.517    0.000)   38.702
   5.240   (   1.691    2.929    0.000)    3.382
   5.242   (   1.612    2.792    0.000)    3.224
   5.385   (   5.454    9.446    0.000)   10.907
   5.415   (   5.983   10.363    0.000)   11.966
   7.095   (  -4.568   -7.912    0.000)    9.136
   7.115   (  -4.905   -8.496    0.000)    9.810
   7.325   (  -0.001   -0.002    0.000)    0.003
   7.338   (  -0.114   -0.198    0.000)    0.229
   7.598   (   2.807    4.862    0.000)    5.614
   7.637   (   2.695    4.667    0.000)    5.389
   9.053   (  -3.714   -6.433    0.000)    7.429
   9.069   (  -3.904   -6.762    0.000)    7.808
======================= Grid point 21 (12/111) =======================
q-point: ( 0.11  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.901   (  17.352   19.578    0.000)   26.160
   0.921   (  13.181   14.743    0.000)   19.776
   1.498   (  15.425   21.901    0.000)   26.788
   1.586   (  14.543   20.950    0.000)   25.503
   2.220   (  24.579   26.870    0.000)   36.416
   2.247   (  24.798   25.071    0.000)   35.263
   5.310   (   2.639    3.674    0.000)    4.523
   5.312   (   2.790    3.804    0.000)    4.717
   5.607   (   8.955    9.896    0.000)   13.346
   5.655   (   9.549   10.543    0.000)   14.225
   6.927   (  -6.035   -7.204    0.000)    9.398
   6.937   (  -6.389   -7.646    0.000)    9.964
   7.322   (  -0.299   -0.049    0.000)    0.303
   7.332   (  -0.394   -0.174    0.000)    0.430
   7.691   (   3.189    3.554    0.000)    4.775
   7.727   (   3.118    3.502    0.000)    4.689
   8.890   (  -7.010   -8.053    0.000)   10.676
   8.899   (  -7.176   -8.254    0.000)   10.937
======================= Grid point 22 (13/111) =======================
q-point: ( 0.17  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.371   (  19.773   17.969    0.000)   26.718
   1.429   (  21.519   19.641    0.000)   29.134
   1.945   (  13.031   19.426    0.000)   23.392
   2.013   (  12.287   19.072    0.000)   22.687
   2.843   (  22.370   18.064    0.000)   28.753
   2.877   (  23.472   18.066    0.000)   29.619
   5.398   (   2.937    4.327    0.000)    5.229
   5.405   (   3.200    4.470    0.000)    5.497
   5.866   (  10.189    8.908    0.000)   13.534
   5.928   (  10.713    9.348    0.000)   14.218
   6.757   (  -5.946   -6.168    0.000)    8.568
   6.766   (  -5.334   -6.041    0.000)    8.059
   7.312   (  -0.843   -0.082    0.000)    0.847
   7.321   (  -0.825   -0.095    0.000)    0.830
   7.774   (   3.008    2.508    0.000)    3.916
   7.809   (   2.985    2.577    0.000)    3.943
   8.663   (  -9.599   -9.070    0.000)   13.206
   8.670   (  -9.572   -9.119    0.000)   13.220
======================= Grid point 23 (14/111) =======================
q-point: ( 0.22  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.909   (  22.170   17.625    0.000)   28.322
   1.997   (  22.889   18.222    0.000)   29.257
   2.316   (   9.528   16.988    0.000)   19.477
   2.370   (   8.940   16.850    0.000)   19.075
   3.320   (  17.001   10.229    0.000)   19.841
   3.372   (  17.682   10.683    0.000)   20.658
   5.490   (   2.567    4.728    0.000)    5.380
   5.503   (   2.889    4.833    0.000)    5.630
   6.113   (   9.291    6.861    0.000)   11.550
   6.187   (   9.787    7.208    0.000)   12.155
   6.623   (  -3.494   -4.049    0.000)    5.348
   6.646   (  -2.774   -3.881    0.000)    4.770
   7.292   (  -1.452   -0.008    0.000)    1.452
   7.302   (  -1.340    0.062    0.000)    1.341
   7.841   (   2.581    1.479    0.000)    2.974
   7.876   (   2.486    1.528    0.000)    2.918
   8.395   ( -10.961   -9.195    0.000)   14.307
   8.404   ( -10.814   -9.206    0.000)   14.202
======================= Grid point 24 (15/111) =======================
q-point: ( 0.28  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.448   (  21.539   15.480    0.000)   26.524
   2.548   (  21.900   15.683    0.000)   26.936
   2.600   (   5.812   14.477    0.000)   15.600
   2.646   (   5.489   14.345    0.000)   15.359
   3.613   (   9.741    3.259    0.000)   10.272
   3.676   (   9.829    3.405    0.000)   10.402
   5.575   (   1.698    4.908    0.000)    5.194
   5.594   (   1.993    4.955    0.000)    5.341
   6.302   (   6.315    3.934    0.000)    7.440
   6.388   (   6.828    4.231    0.000)    8.033
   6.570   (   0.448   -1.369    0.000)    1.440
   6.603   (   0.843   -1.449    0.000)    1.676
   7.264   (  -1.924    0.249    0.000)    1.940
   7.277   (  -1.757    0.400    0.000)    1.802
   7.889   (   1.903    0.388    0.000)    1.942
   7.921   (   1.677    0.310    0.000)    1.706
   8.125   ( -10.326   -8.092    0.000)   13.119
   8.135   ( -10.275   -8.197    0.000)   13.145
======================= Grid point 25 (16/111) =======================
q-point: ( 0.33  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.799   (   2.409   11.894    0.000)   12.136
   2.840   (   2.402   11.599    0.000)   11.845
   2.928   (  18.541   12.219    0.000)   22.205
   3.034   (  18.800   12.221    0.000)   22.423
   3.727   (   3.358   -1.580    0.000)    3.712
   3.785   (   2.993   -1.975    0.000)    3.586
   5.642   (   0.543    4.961    0.000)    4.991
   5.666   (   0.759    4.954    0.000)    5.012
   6.394   (   1.767    0.550    0.000)    1.850
   6.490   (   2.181    0.680    0.000)    2.284
   6.605   (   3.804    0.339    0.000)    3.819
   6.638   (   3.713   -0.014    0.000)    3.713
   7.234   (  -2.101    0.819    0.000)    2.255
   7.252   (  -1.934    1.013    0.000)    2.183
   7.904   (  -7.253   -5.720    0.000)    9.237
   7.913   (  -7.312   -5.883    0.000)    9.385
   7.913   (   1.149   -0.593    0.000)    1.293
   7.939   (   0.802   -0.823    0.000)    1.149
======================= Grid point 26 (17/111) =======================
q-point: ( 0.39  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.919   (  -0.381    9.321    0.000)    9.328
   2.958   (  -0.108    8.740    0.000)    8.741
   3.301   (  13.127    7.931    0.000)   15.337
   3.410   (  13.374    7.767    0.000)   15.466
   3.726   (   0.160   -3.182    0.000)    3.187
   3.771   (  -0.157   -3.984    0.000)    3.987
   5.690   (  -0.718    4.947    0.000)    4.999
   5.716   (  -0.588    4.920    0.000)    4.955
   6.382   (  -1.756   -1.736    0.000)    2.470
   6.483   (  -1.630   -1.960    0.000)    2.549
   6.684   (   4.668   -0.085    0.000)    4.669
   6.709   (   4.339   -0.508    0.000)    4.369
   7.211   (  -1.997    1.736    0.000)    2.646
   7.234   (  -1.896    1.909    0.000)    2.690
   7.773   (  -3.160   -3.118    0.000)    4.440
   7.780   (  -3.175   -3.176    0.000)    4.491
   7.919   (   0.687   -1.193    0.000)    1.376
   7.936   (   0.274   -1.508    0.000)    1.533
======================= Grid point 27 (18/111) =======================
q-point: ( 0.44  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.973   (  -2.602    7.139    0.000)    7.598
   3.012   (  -2.146    6.318    0.000)    6.672
   3.511   (   5.001    2.575    0.000)    5.624
   3.622   (   5.232    2.267    0.000)    5.702
   3.697   (   0.232   -2.377    0.000)    2.389
   3.732   (   0.398   -3.266    0.000)    3.291
   5.714   (  -1.991    4.745    0.000)    5.146
   5.742   (  -1.938    4.726    0.000)    5.108
   6.333   (  -1.195   -1.416    0.000)    1.853
   6.432   (  -1.101   -1.828    0.000)    2.134
   6.743   (   3.074   -1.968    0.000)    3.650
   6.760   (   2.905   -2.214    0.000)    3.652
   7.201   (  -1.782    2.598    0.000)    3.151
   7.226   (  -1.780    2.714    0.000)    3.246
   7.724   (  -0.352   -1.492    0.000)    1.533
   7.731   (  -0.357   -1.494    0.000)    1.536
   7.918   (   0.739   -1.299    0.000)    1.495
   7.925   (   0.402   -1.574    0.000)    1.625
======================= Grid point 39 (19/111) =======================
q-point: ( 0.11  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.289   (  12.827   22.217    0.000)   25.654
   1.340   (  14.227   24.642    0.000)   28.454
   1.941   (  12.927   22.391    0.000)   25.855
   2.015   (  12.609   21.840    0.000)   25.219
   2.734   (  13.938   24.141    0.000)   27.876
   2.748   (  14.863   25.743    0.000)   29.726
   5.397   (   2.875    4.979    0.000)    5.749
   5.402   (   2.984    5.169    0.000)    5.968
   5.821   (   6.640   11.501    0.000)   13.280
   5.880   (   6.983   12.095    0.000)   13.966
   6.779   (  -4.622   -8.005    0.000)    9.244
   6.780   (  -4.311   -7.467    0.000)    8.622
   7.320   (  -0.145   -0.251    0.000)    0.290
   7.328   (  -0.143   -0.248    0.000)    0.286
   7.760   (   1.918    3.323    0.000)    3.837
   7.796   (   1.953    3.383    0.000)    3.906
   8.700   (  -6.405  -11.094    0.000)   12.810
   8.707   (  -6.458  -11.186    0.000)   12.916
======================= Grid point 40 (20/111) =======================
q-point: ( 0.17  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.765   (  17.053   21.716    0.000)   27.611
   1.848   (  17.731   22.613    0.000)   28.735
   2.353   (   9.866   21.276    0.000)   23.452
   2.417   (   9.358   21.139    0.000)   23.118
   3.163   (  13.948   14.061    0.000)   19.805
   3.208   (  14.920   15.100    0.000)   21.228
   5.500   (   2.851    5.891    0.000)    6.544
   5.509   (   3.063    6.022    0.000)    6.756
   6.044   (   7.466    8.974    0.000)   11.673
   6.115   (   7.842    9.419    0.000)   12.257
   6.637   (  -3.750   -5.765    0.000)    6.878
   6.650   (  -3.182   -5.444    0.000)    6.306
   7.312   (  -0.670    0.002    0.000)    0.670
   7.322   (  -0.614    0.176    0.000)    0.638
   7.819   (   1.901    1.962    0.000)    2.732
   7.856   (   1.873    2.043    0.000)    2.771
   8.462   (  -8.508  -11.211    0.000)   14.074
   8.468   (  -8.429  -11.242    0.000)   14.051
======================= Grid point 41 (21/111) =======================
q-point: ( 0.22  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 492
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.272   (  18.534   19.067    0.000)   26.591
   2.369   (  18.859   19.352    0.000)   27.021
   2.685   (   5.345   19.868    0.000)   20.575
   2.739   (   4.703   19.862    0.000)   20.411
   3.468   (  10.894    4.575    0.000)   11.816
   3.531   (  11.086    5.139    0.000)   12.219
   5.604   (   2.009    6.708    0.000)    7.003
   5.619   (   2.282    6.768    0.000)    7.142
   6.237   (   6.400    5.571    0.000)    8.485
   6.318   (   6.811    5.933    0.000)    9.033
   6.552   (  -0.931   -3.020    0.000)    3.160
   6.576   (  -0.430   -2.972    0.000)    3.003
   7.297   (  -1.286    0.572    0.000)    1.408
   7.311   (  -1.199    0.855    0.000)    1.473
   7.861   (   1.652    0.459    0.000)    1.714
   7.896   (   1.458    0.418    0.000)    1.517
   8.203   (  -9.184  -10.205    0.000)   13.729
   8.211   (  -9.082  -10.320    0.000)   13.747
======================= Grid point 42 (22/111) =======================
q-point: ( 0.28  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.752   (  17.599   15.331    0.000)   23.340
   2.853   (  17.797   15.282    0.000)   23.458
   2.921   (   0.946   17.354    0.000)   17.380
   2.963   (   0.395   17.089    0.000)   17.093
   3.625   (   6.361   -1.811    0.000)    6.614
   3.686   (   5.867   -2.249    0.000)    6.283
   5.697   (   0.676    7.374    0.000)    7.405
   5.716   (   0.932    7.363    0.000)    7.421
   6.359   (   3.407    1.720    0.000)    3.817
   6.450   (   3.786    1.909    0.000)    4.240
   6.546   (   2.426   -1.049    0.000)    2.643
   6.576   (   2.585   -1.343    0.000)    2.913
   7.281   (  -1.771    1.566    0.000)    2.364
   7.299   (  -1.662    1.860    0.000)    2.494
   7.883   (   1.293   -0.982    0.000)    1.624
   7.913   (   0.950   -1.207    0.000)    1.536
   7.966   (  -7.730   -7.990    0.000)   11.117
   7.973   (  -7.711   -8.210    0.000)   11.263
======================= Grid point 43 (23/111) =======================
q-point: ( 0.33  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 492
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.059   (  -2.248   13.112    0.000)   13.304
   3.088   (  -2.517   12.253    0.000)   12.509
   3.156   (  14.317   10.897    0.000)   17.992
   3.259   (  14.508   10.559    0.000)   17.944
   3.671   (   2.560   -2.934    0.000)    3.894
   3.711   (   2.039   -4.613    0.000)    5.043
   5.770   (  -0.870    7.873    0.000)    7.921
   5.793   (  -0.658    7.810    0.000)    7.838
   6.384   (  -0.327   -1.504    0.000)    1.539
   6.480   (  -0.178   -1.755    0.000)    1.764
   6.598   (   4.350   -1.269    0.000)    4.532
   6.622   (   4.190   -1.728    0.000)    4.532
   7.271   (  -2.006    2.918    0.000)    3.541
   7.293   (  -1.905    3.145    0.000)    3.677
   7.797   (  -4.363   -5.221    0.000)    6.804
   7.802   (  -4.310   -5.312    0.000)    6.840
   7.888   (   1.040   -2.006    0.000)    2.259
   7.908   (   0.638   -2.359    0.000)    2.444
======================= Grid point 44 (24/111) =======================
q-point: ( 0.39  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.112   (  -3.383    8.305    0.000)    8.967
   3.129   (  -3.223    6.835    0.000)    7.557
   3.435   (   8.342    5.771    0.000)   10.144
   3.537   (   8.615    5.171    0.000)   10.048
   3.675   (   0.155   -0.275    0.000)    0.316
   3.692   (   0.314   -2.507    0.000)    2.527
   5.820   (  -2.538    8.076    0.000)    8.465
   5.845   (  -2.368    8.007    0.000)    8.350
   6.339   (  -1.691   -2.445    0.000)    2.973
   6.431   (  -1.667   -3.169    0.000)    3.580
   6.652   (   4.257   -3.124    0.000)    5.281
   6.669   (   4.055   -3.485    0.000)    5.347
   7.270   (  -2.194    4.054    0.000)    4.609
   7.296   (  -2.151    4.201    0.000)    4.719
   7.713   (  -1.002   -3.111    0.000)    3.268
   7.719   (  -0.973   -3.114    0.000)    3.263
   7.883   (   1.121   -2.394    0.000)    2.643
   7.893   (   0.747   -2.742    0.000)    2.842
======================= Grid point 45 (25/111) =======================
q-point: ( 0.44  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.122   (  -3.719    6.441    0.000)    7.437
   3.131   (  -2.491    4.314    0.000)    4.981
   3.536   (  -0.154    0.266    0.000)    0.307
   3.638   (   0.241   -0.417    0.000)    0.481
   3.675   (  -0.809    1.402    0.000)    1.618
   3.681   (   0.362   -0.628    0.000)    0.725
   5.839   (  -4.358    7.548    0.000)    8.716
   5.866   (  -4.302    7.451    0.000)    8.604
   6.308   (   0.618   -1.070    0.000)    1.236
   6.396   (   1.032   -1.787    0.000)    2.064
   6.674   (   2.773   -4.802    0.000)    5.545
   6.687   (   2.843   -4.924    0.000)    5.686
   7.272   (  -2.506    4.341    0.000)    5.012
   7.299   (  -2.555    4.426    0.000)    5.111
   7.691   (   1.109   -1.921    0.000)    2.219
   7.697   (   1.112   -1.925    0.000)    2.223
   7.882   (   1.359   -2.354    0.000)    2.718
   7.884   (   1.418   -2.456    0.000)    2.836
======================= Grid point 59 (26/111) =======================
q-point: ( 0.17  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.206   (  13.191   22.848    0.000)   26.382
   2.301   (  13.408   23.223    0.000)   26.815
   2.751   (  10.791   18.690    0.000)   21.581
   2.812   (  10.704   18.540    0.000)   21.409
   3.384   (   4.821    8.350    0.000)    9.642
   3.450   (   5.316    9.208    0.000)   10.632
   5.623   (   3.716    6.436    0.000)    7.432
   5.635   (   3.817    6.611    0.000)    7.634
   6.207   (   4.258    7.376    0.000)    8.517
   6.287   (   4.524    7.836    0.000)    9.049
   6.544   (  -1.943   -3.365    0.000)    3.885
   6.564   (  -1.750   -3.031    0.000)    3.500
   7.312   (   0.043    0.075    0.000)    0.087
   7.327   (   0.195    0.338    0.000)    0.391
   7.848   (   0.585    1.013    0.000)    1.169
   7.885   (   0.527    0.914    0.000)    1.055
   8.232   (  -6.918  -11.982    0.000)   13.836
   8.238   (  -6.942  -12.023    0.000)   13.883
======================= Grid point 60 (27/111) =======================
q-point: ( 0.22  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.647   (  14.160   18.823    0.000)   23.555
   2.745   (  14.231   18.794    0.000)   23.574
   3.065   (   4.341   17.994    0.000)   18.511
   3.121   (   3.502   18.371    0.000)   18.702
   3.499   (   5.318   -0.930    0.000)    5.399
   3.569   (   4.510   -1.089    0.000)    4.640
   5.746   (   2.606    7.596    0.000)    8.031
   5.762   (   2.769    7.705    0.000)    8.188
   6.324   (   3.087    3.098    0.000)    4.374
   6.412   (   3.344    3.400    0.000)    4.769
   6.509   (   0.422   -1.324    0.000)    1.389
   6.533   (   0.554   -1.410    0.000)    1.515
   7.316   (  -0.379    1.345    0.000)    1.397
   7.335   (  -0.279    1.641    0.000)    1.664
   7.858   (   0.452   -0.748    0.000)    0.874
   7.891   (   0.184   -1.023    0.000)    1.040
   8.004   (  -6.897   -9.837    0.000)   12.014
   8.008   (  -6.903  -10.058    0.000)   12.199
======================= Grid point 61 (28/111) =======================
q-point: ( 0.28  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.043   (  13.045   14.233    0.000)   19.307
   3.141   (  13.056   13.819    0.000)   19.012
   3.220   (  -2.711    9.522    0.000)    9.900
   3.263   (  -3.854   11.110    0.000)   11.760
   3.589   (   5.805    1.091    0.000)    5.907
   3.620   (   4.607   -1.812    0.000)    4.951
   5.856   (   0.832    8.628    0.000)    8.668
   5.875   (   0.971    8.668    0.000)    8.722
   6.367   (   0.699   -0.839    0.000)    1.092
   6.459   (   0.739   -0.974    0.000)    1.223
   6.523   (   2.459   -1.511    0.000)    2.886
   6.546   (   2.400   -1.873    0.000)    3.045
   7.327   (  -0.854    3.091    0.000)    3.207
   7.351   (  -0.793    3.318    0.000)    3.412
   7.817   (  -5.045   -7.107    0.000)    8.715
   7.820   (  -4.917   -7.176    0.000)    8.699
   7.851   (   0.559   -2.209    0.000)    2.279
   7.874   (   0.166   -2.617    0.000)    2.622
======================= Grid point 62 (29/111) =======================
q-point: ( 0.33  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.197   (  -2.774   -1.313    0.000)    3.069
   3.238   (  -4.222    0.325    0.000)    4.235
   3.354   (   9.474    9.302    0.000)   13.277
   3.446   (   9.552    8.477    0.000)   12.771
   3.689   (   1.381    5.915    0.000)    6.074
   3.707   (   1.037    7.917    0.000)    7.985
   5.941   (  -1.197    9.253    0.000)    9.330
   5.963   (  -1.057    9.178    0.000)    9.238
   6.339   (  -1.302   -2.816    0.000)    3.103
   6.425   (  -1.541   -3.655    0.000)    3.966
   6.551   (   3.427   -3.493    0.000)    4.893
   6.567   (   3.270   -3.866    0.000)    5.064
   7.346   (  -1.481    4.525    0.000)    4.761
   7.372   (  -1.430    4.658    0.000)    4.872
   7.699   (  -1.997   -4.690    0.000)    5.097
   7.705   (  -1.922   -4.667    0.000)    5.047
   7.838   (   0.960   -2.957    0.000)    3.109
   7.850   (   0.594   -3.362    0.000)    3.414
======================= Grid point 63 (30/111) =======================
q-point: ( 0.39  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.145   (   0.457   -3.959    0.000)    3.986
   3.174   (  -0.599   -3.163    0.000)    3.219
   3.531   (   3.092    4.110    0.000)    5.143
   3.616   (   3.470    3.051    0.000)    4.621
   3.748   (  -3.137    8.632    0.000)    9.185
   3.775   (  -3.541    9.733    0.000)   10.357
   5.993   (  -3.557    9.080    0.000)    9.751
   6.016   (  -3.351    8.858    0.000)    9.470
   6.292   (  -0.619   -2.191    0.000)    2.276
   6.364   (  -0.475   -3.471    0.000)    3.504
   6.567   (   3.536   -5.306    0.000)    6.377
   6.579   (   3.489   -5.521    0.000)    6.531
   7.361   (  -2.326    4.918    0.000)    5.440
   7.388   (  -2.309    4.980    0.000)    5.489
   7.651   (   0.562   -3.211    0.000)    3.260
   7.656   (   0.586   -3.237    0.000)    3.290
   7.829   (   1.541   -3.007    0.000)    3.379
   7.833   (   1.255   -3.273    0.000)    3.506
======================= Grid point 78 (31/111) =======================
q-point: ( 0.22  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 273
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.001   (   9.844   17.051    0.000)   19.689
   3.096   (   9.637   16.691    0.000)   19.273
   3.301   (   0.742    1.286    0.000)    1.485
   3.385   (   1.351    2.340    0.000)    2.702
   3.516   (   4.578    7.929    0.000)    9.155
   3.559   (   3.561    6.169    0.000)    7.123
   5.894   (   4.188    7.254    0.000)    8.376
   5.913   (   4.280    7.414    0.000)    8.561
   6.355   (   0.031    0.053    0.000)    0.061
   6.446   (   0.008    0.014    0.000)    0.017
   6.494   (  -0.291   -0.504    0.000)    0.582
   6.515   (  -0.428   -0.742    0.000)    0.856
   7.351   (   1.277    2.213    0.000)    2.555
   7.375   (   1.381    2.391    0.000)    2.761
   7.826   (  -4.778   -8.275    0.000)    9.556
   7.826   (  -4.645   -8.045    0.000)    9.289
   7.834   (  -0.914   -1.584    0.000)    1.829
   7.859   (  -1.201   -2.080    0.000)    2.402
======================= Grid point 79 (32/111) =======================
q-point: ( 0.28  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.206   (  -4.001   -7.240    0.000)    8.272
   3.275   (  -5.087   -7.726    0.000)    9.251
   3.305   (   8.608   12.217    0.000)   14.945
   3.389   (   8.328   11.344    0.000)   14.073
   3.734   (   5.141   12.428    0.000)   13.450
   3.767   (   4.901   12.665    0.000)   13.580
   6.024   (   2.002    8.137    0.000)    8.380
   6.045   (   1.982    8.311    0.000)    8.544
   6.331   (  -1.027   -2.584    0.000)    2.780
   6.415   (  -1.486   -3.323    0.000)    3.640
   6.482   (   0.698   -2.678    0.000)    2.768
   6.498   (   0.550   -3.006    0.000)    3.056
   7.398   (   0.511    3.968    0.000)    4.001
   7.425   (   0.542    4.070    0.000)    4.106
   7.691   (  -2.980   -5.651    0.000)    6.389
   7.695   (  -2.897   -5.569    0.000)    6.277
   7.802   (  -0.265   -2.694    0.000)    2.707
   7.816   (  -0.607   -3.173    0.000)    3.231
======================= Grid point 80 (33/111) =======================
q-point: ( 0.33  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 492
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.082   (  -1.003   -7.855    0.000)    7.919
   3.127   (  -1.960   -8.613    0.000)    8.833
   3.519   (   4.819    7.373    0.000)    8.808
   3.589   (   4.777    6.000    0.000)    7.669
   3.897   (  -0.498   12.397    0.000)   12.407
   3.930   (  -0.499   12.312    0.000)   12.322
   6.116   (  -0.526    7.971    0.000)    7.989
   6.136   (  -0.847    7.823    0.000)    7.868
   6.286   (  -1.019   -2.111    0.000)    2.343
   6.347   (  -1.762   -3.921    0.000)    4.298
   6.466   (   1.918   -4.816    0.000)    5.184
   6.477   (   1.825   -5.084    0.000)    5.402
   7.438   (  -0.975    4.643    0.000)    4.744
   7.466   (  -0.952    4.671    0.000)    4.767
   7.614   (  -0.579   -3.967    0.000)    4.009
   7.619   (  -0.561   -4.017    0.000)    4.056
   7.777   (   0.820   -3.055    0.000)    3.163
   7.783   (   0.586   -3.359    0.000)    3.410
======================= Grid point 81 (34/111) =======================
q-point: ( 0.39  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.034   (   3.434   -5.948    0.000)    6.868
   3.064   (   3.628   -6.284    0.000)    7.257
   3.598   (  -1.611    2.790    0.000)    3.222
   3.660   (  -0.849    1.471    0.000)    1.699
   3.954   (  -6.034   10.451    0.000)   12.068
   3.987   (  -5.941   10.289    0.000)   11.881
   6.156   (  -3.967    6.871    0.000)    7.934
   6.165   (  -2.937    5.087    0.000)    5.873
   6.268   (  -0.344    0.595    0.000)    0.687
   6.304   (   1.321   -2.288    0.000)    2.642
   6.458   (   3.074   -5.324    0.000)    6.147
   6.466   (   3.164   -5.481    0.000)    6.329
   7.453   (  -2.458    4.258    0.000)    4.916
   7.481   (  -2.445    4.235    0.000)    4.890
   7.590   (   1.620   -2.805    0.000)    3.239
   7.595   (   1.666   -2.885    0.000)    3.331
   7.770   (   1.687   -2.922    0.000)    3.374
   7.771   (   1.666   -2.885    0.000)    3.331
======================= Grid point 98 (35/111) =======================
q-point: ( 0.28  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 270
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.049   (  -4.274   -7.403    0.000)    8.549
   3.099   (  -5.012   -8.682    0.000)   10.025
   3.514   (   5.097    8.828    0.000)   10.194
   3.577   (   4.352    7.538    0.000)    8.704
   3.967   (   5.866   10.161    0.000)   11.733
   3.999   (   5.743    9.948    0.000)   11.486
   6.158   (   2.845    4.928    0.000)    5.690
   6.186   (   3.114    5.393    0.000)    6.227
   6.283   (  -0.889   -1.539    0.000)    1.778
   6.344   (  -1.871   -3.240    0.000)    3.741
   6.424   (  -1.723   -2.985    0.000)    3.447
   6.434   (  -1.877   -3.250    0.000)    3.753
   7.472   (   1.885    3.266    0.000)    3.771
   7.500   (   1.893    3.279    0.000)    3.787
   7.598   (  -2.129   -3.688    0.000)    4.259
   7.603   (  -2.156   -3.735    0.000)    4.313
   7.753   (  -1.219   -2.111    0.000)    2.437
   7.758   (  -1.422   -2.463    0.000)    2.844
======================= Grid point 99 (36/111) =======================
q-point: ( 0.33  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.946   (  -0.114   -5.182    0.000)    5.183
   2.971   (  -0.629   -6.385    0.000)    6.416
   3.638   (   1.345    4.636    0.000)    4.828
   3.676   (   1.266    2.807    0.000)    3.079
   4.108   (  -0.331    8.096    0.000)    8.102
   4.136   (  -0.303    7.877    0.000)    7.883
   6.197   (  -0.008    0.100    0.000)    0.100
   6.227   (  -0.440   -0.125    0.000)    0.457
   6.294   (  -0.043    2.857    0.000)    2.858
   6.321   (  -0.941    2.935    0.000)    3.082
   6.378   (   0.572   -3.646    0.000)    3.690
   6.384   (   0.537   -3.905    0.000)    3.942
   7.517   (  -0.436    3.167    0.000)    3.196
   7.545   (  -0.443    3.145    0.000)    3.176
   7.549   (  -0.068   -2.491    0.000)    2.492
   7.552   (  -0.069   -2.631    0.000)    2.632
   7.724   (   0.409   -2.217    0.000)    2.254
   7.725   (   0.360   -2.327    0.000)    2.354
======================= Grid point 117 (37/111) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 111
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.899   (  -0.319   -0.553    0.000)    0.638
   2.899   (   0.319    0.553    0.000)    0.638
   3.695   (  -0.503   -0.871    0.000)    1.005
   3.695   (   0.503    0.871    0.000)    1.005
   4.190   (  -0.000   -0.000    0.000)    0.000
   4.216   (  -0.000   -0.000    0.000)    0.000
   6.201   (  -0.471   -0.816    0.000)    0.943
   6.201   (   0.471    0.816    0.000)    0.943
   6.335   (  -0.000   -0.000    0.000)    0.000
   6.341   (  -0.071   -0.123    0.000)    0.142
   6.341   (   0.071    0.123    0.000)    0.142
   6.362   (  -0.000   -0.000    0.000)    0.000
   7.524   (  -0.044   -0.076    0.000)    0.087
   7.524   (   0.044    0.076    0.000)    0.087
   7.551   (  -0.000   -0.000    0.000)    0.000
   7.579   (  -0.000   -0.000    0.000)    0.000
   7.700   (  -0.003   -0.005    0.000)    0.006
   7.700   (   0.003    0.005    0.000)    0.006
======================= Grid point 343 (38/111) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 74
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.192   (   0.000   -0.000    9.783)    9.783
   0.192   (   0.000    0.000    9.783)    9.783
   0.249   (   0.000    0.000   12.649)   12.649
   0.465   (   0.000    0.000   -4.056)    4.056
   0.465   (   0.000    0.000   -4.056)    4.056
   0.601   (  -0.000   -0.000   -5.270)    5.270
   5.182   (   0.000   -0.000    0.069)    0.069
   5.182   (   0.000   -0.000    0.069)    0.069
   5.185   (   0.000    0.000   -0.069)    0.069
   5.185   (   0.000   -0.000   -0.069)    0.069
   7.326   (   0.000    0.000    0.272)    0.272
   7.326   (  -0.000    0.000    0.272)    0.272
   7.339   (   0.000    0.000   -0.271)    0.271
   7.339   (   0.000   -0.000   -0.271)    0.271
   7.415   (   0.000   -0.000    0.516)    0.516
   7.440   (   0.000    0.000   -0.515)    0.515
   9.195   (   0.000   -0.000    0.212)    0.212
   9.205   (   0.000    0.000   -0.209)    0.209
======================= Grid point 344 (39/111) =======================
q-point: ( 0.06  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.314   (   7.660    4.422    7.917)   11.871
   0.588   (   0.448    0.258   -4.654)    4.683
   0.614   (  24.004   13.859    3.114)   27.891
   0.748   (  19.861   11.467   -2.562)   23.077
   0.913   (  37.702   21.767    2.071)   43.584
   1.001   (  33.900   19.572   -1.650)   39.180
   5.201   (   1.694    0.978    0.059)    1.957
   5.204   (   1.651    0.953   -0.059)    1.908
   5.257   (   6.334    3.657    0.276)    7.319
   5.267   (   6.922    3.997   -0.187)    7.996
   7.221   (  -6.986   -4.033    0.170)    8.068
   7.239   (  -6.930   -4.001   -0.854)    8.047
   7.327   (  -0.001   -0.001    0.246)    0.246
   7.338   (  -0.105   -0.061   -0.245)    0.273
   7.516   (   5.385    3.109    0.360)    6.229
   7.542   (   5.685    3.282   -0.957)    6.634
   9.144   (  -4.047   -2.337    0.427)    4.693
   9.157   (  -4.107   -2.371   -0.226)    4.748
======================= Grid point 345 (40/111) =======================
q-point: ( 0.11  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.611   (  17.942   10.359    1.940)   20.808
   0.704   (  11.048    6.379   -1.955)   12.906
   1.152   (  22.956   13.254    1.699)   26.562
   1.230   (  21.537   12.435   -1.595)   24.920
   1.747   (  35.274   20.366    0.908)   40.741
   1.786   (  34.120   19.699   -0.772)   39.406
   5.256   (   3.076    1.776    0.028)    3.552
   5.258   (   3.001    1.733   -0.027)    3.465
   5.453   (  10.544    6.087    0.563)   12.188
   5.478   (  11.190    6.460   -0.498)   12.930
   7.049   (  -8.039   -4.641    0.193)    9.285
   7.061   (  -8.406   -4.853   -0.313)    9.711
   7.325   (  -0.166   -0.096    0.185)    0.266
   7.334   (  -0.298   -0.172   -0.184)    0.391
   7.632   (   4.723    2.727    0.476)    5.474
   7.658   (   4.627    2.671   -0.654)    5.382
   9.014   (  -7.401   -4.273    0.229)    8.549
   9.023   (  -7.611   -4.394   -0.182)    8.791
======================= Grid point 346 (41/111) =======================
q-point: ( 0.17  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.078   (  20.975   12.110    0.203)   24.220
   1.091   (  23.645   13.651   -0.348)   27.305
   1.642   (  20.121   11.617    1.179)   23.264
   1.697   (  19.387   11.193   -1.146)   22.415
   2.491   (  30.031   17.339    0.598)   34.682
   2.517   (  29.956   17.295   -0.517)   34.594
   5.335   (   3.779    2.182    0.057)    4.365
   5.338   (   3.949    2.280   -0.056)    4.561
   5.715   (  12.177    7.030    0.805)   14.084
   5.753   (  12.676    7.318   -0.795)   14.658
   6.867   (  -7.904   -4.563    0.105)    9.127
   6.872   (  -7.930   -4.578   -0.122)    9.157
   7.317   (  -0.550   -0.318    0.137)    0.650
   7.324   (  -0.600   -0.347   -0.137)    0.707
   7.730   (   3.980    2.298    0.504)    4.623
   7.755   (   3.948    2.279   -0.535)    4.590
   8.812   ( -10.237   -5.910    0.130)   11.821
   8.818   ( -10.328   -5.963   -0.119)   11.926
======================= Grid point 347 (42/111) =======================
q-point: ( 0.22  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.605   (  24.696   14.258    1.067)   28.536
   1.656   (  25.697   14.836   -1.117)   29.693
   2.058   (  16.525    9.541    0.892)   19.103
   2.100   (  16.075    9.281   -0.880)   18.582
   3.092   (  22.711   13.112    0.701)   26.233
   3.124   (  23.145   13.363   -0.650)   26.734
   5.423   (   3.942    2.276    0.153)    4.554
   5.431   (   4.159    2.401   -0.150)    4.804
   5.988   (  11.543    6.664    1.005)   13.366
   6.036   (  11.943    6.895   -1.016)   13.828
   6.703   (  -6.140   -3.545    0.258)    7.095
   6.715   (  -5.682   -3.281   -0.256)    6.566
   7.299   (  -1.024   -0.591    0.139)    1.191
   7.306   (  -0.979   -0.565   -0.138)    1.139
   7.812   (   3.249    1.876    0.516)    3.786
   7.837   (   3.242    1.872   -0.512)    3.778
   8.557   ( -12.076   -6.972    0.127)   13.945
   8.563   ( -12.009   -6.934   -0.120)   13.868
======================= Grid point 348 (43/111) =======================
q-point: ( 0.28  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.168   (  24.546   14.172    1.400)   28.378
   2.235   (  24.959   14.410   -1.408)   28.854
   2.389   (  12.673    7.317    0.723)   14.652
   2.423   (  12.431    7.177   -0.719)   14.372
   3.508   (  13.889    8.019    0.877)   16.061
   3.549   (  14.193    8.194   -0.861)   16.412
   5.509   (   3.555    2.052    0.255)    4.113
   5.521   (   3.757    2.169   -0.253)    4.345
   6.223   (   8.894    5.135    1.188)   10.339
   6.280   (   9.274    5.354   -1.201)   10.775
   6.603   (  -2.442   -1.410    0.470)    2.858
   6.625   (  -2.090   -1.207   -0.459)    2.456
   7.272   (  -1.395   -0.806    0.173)    1.620
   7.280   (  -1.299   -0.750   -0.172)    1.510
   7.876   (   2.322    1.340    0.504)    2.728
   7.899   (   2.249    1.299   -0.492)    2.643
   8.278   ( -12.253   -7.075    0.165)   14.150
   8.285   ( -12.196   -7.041   -0.154)   14.084
======================= Grid point 349 (44/111) =======================
q-point: ( 0.33  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.633   (   8.948    5.166    0.650)   10.353
   2.664   (   8.881    5.128   -0.649)   10.276
   2.698   (  21.841   12.610    1.549)   25.267
   2.771   (  22.057   12.734   -1.543)   25.516
   3.722   (   5.267    3.041    0.918)    6.150
   3.766   (   5.190    2.997   -0.933)    6.065
   5.582   (   2.794    1.613    0.338)    3.243
   5.598   (   2.938    1.696   -0.337)    3.409
   6.376   (   4.382    2.530    1.362)    5.240
   6.441   (   4.737    2.735   -1.388)    5.643
   6.596   (   1.725    0.996    0.572)    2.072
   6.623   (   1.726    0.997   -0.534)    2.064
   7.238   (  -1.514   -0.874    0.230)    1.763
   7.249   (  -1.382   -0.798   -0.227)    1.611
   7.916   (   1.213    0.700    0.441)    1.468
   7.937   (   1.025    0.592   -0.437)    1.261
   8.020   ( -10.043   -5.798    0.172)   11.598
   8.028   ( -10.110   -5.837   -0.149)   11.675
======================= Grid point 350 (45/111) =======================
q-point: ( 0.39  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.797   (   5.591    3.228    0.647)    6.489
   2.828   (   5.631    3.251   -0.646)    6.535
   3.141   (  16.957    9.790    1.631)   19.648
   3.219   (  17.100    9.873   -1.627)   19.813
   3.770   (  -0.332   -0.192    0.809)    0.895
   3.809   (  -0.647   -0.374   -0.832)    1.118
   5.635   (   1.865    1.077    0.389)    2.188
   5.653   (   1.941    1.121   -0.390)    2.275
   6.417   (  -0.553   -0.319    1.471)    1.604
   6.489   (  -0.377   -0.217   -1.555)    1.615
   6.666   (   3.881    2.241    0.529)    4.513
   6.688   (   3.603    2.080   -0.422)    4.182
   7.206   (  -1.282   -0.740    0.294)    1.509
   7.220   (  -1.158   -0.668   -0.289)    1.367
   7.838   (  -5.841   -3.373    0.120)    6.746
   7.844   (  -5.919   -3.417   -0.119)    6.836
   7.932   (   0.274    0.158    0.338)    0.463
   7.948   (  -0.007   -0.004   -0.314)    0.314
======================= Grid point 351 (46/111) =======================
q-point: ( 0.44  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.890   (   2.657    1.534    0.672)    3.141
   2.922   (   2.710    1.564   -0.672)    3.200
   3.447   (   9.598    5.541    1.674)   11.208
   3.527   (   9.681    5.589   -1.697)   11.307
   3.740   (  -1.631   -0.942    0.677)    2.002
   3.772   (  -1.877   -1.084   -0.670)    2.269
   5.666   (   0.918    0.530    0.412)    1.137
   5.686   (   0.945    0.546   -0.414)    1.167
   6.376   (  -2.354   -1.359    1.485)    3.097
   6.448   (  -2.395   -1.383   -1.590)    3.190
   6.749   (   2.981    1.721    0.394)    3.465
   6.764   (   2.689    1.552   -0.258)    3.115
   7.182   (  -0.731   -0.422    0.343)    0.911
   7.199   (  -0.659   -0.381   -0.336)    0.832
   7.751   (  -2.098   -1.211    0.094)    2.425
   7.756   (  -2.120   -1.224   -0.100)    2.450
   7.933   (  -0.090   -0.052    0.187)    0.214
   7.942   (  -0.403   -0.232   -0.172)    0.496
======================= Grid point 352 (47/111) =======================
q-point: (-0.50  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 182
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.919   (   0.000    0.000    0.686)    0.686
   2.952   (   0.000    0.000   -0.686)    0.686
   3.559   (   0.000    0.000    1.671)    1.671
   3.640   (   0.000    0.000   -1.751)    1.751
   3.719   (   0.000    0.000    0.640)    0.640
   3.748   (   0.000    0.000   -0.574)    0.574
   5.676   (   0.000    0.000    0.418)    0.418
   5.696   (   0.000    0.000   -0.421)    0.421
   6.344   (   0.000    0.000    1.465)    1.465
   6.416   (   0.000    0.000   -1.570)    1.570
   6.785   (  -0.000   -0.000    0.316)    0.316
   6.796   (   0.000    0.000   -0.175)    0.175
   7.174   (   0.000    0.000    0.360)    0.360
   7.191   (   0.000    0.000   -0.353)    0.353
   7.730   (   0.000    0.000    0.086)    0.086
   7.734   (   0.000    0.000   -0.098)    0.098
   7.932   (   0.000    0.000    0.079)    0.079
   7.935   (   0.000    0.000   -0.065)    0.065
======================= Grid point 363 (48/111) =======================
q-point: ( 0.06  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.508   (   8.940   15.485    3.105)   18.148
   0.646   (   4.318    7.479   -2.759)    9.066
   1.026   (  14.199   24.594    1.910)   28.463
   1.112   (  13.184   22.836   -1.733)   26.425
   1.522   (  20.611   35.699    0.977)   41.233
   1.566   (  19.549   33.860   -0.836)   39.107
   5.240   (   1.679    2.909    0.027)    3.359
   5.241   (   1.623    2.812   -0.025)    3.247
   5.390   (   5.516    9.555    0.487)   11.043
   5.411   (   5.900   10.219   -0.406)   11.807
   7.097   (  -4.577   -7.927    0.215)    9.156
   7.111   (  -4.780   -8.279   -0.411)    9.569
   7.327   (  -0.018   -0.031    0.199)    0.202
   7.336   (  -0.098   -0.170   -0.200)    0.280
   7.602   (   2.833    4.906    0.457)    5.684
   7.629   (   2.794    4.839   -0.721)    5.634
   9.056   (  -3.767   -6.524    0.272)    7.538
   9.067   (  -3.889   -6.736   -0.200)    7.781
======================= Grid point 364 (49/111) =======================
q-point: ( 0.11  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.902   (  16.757   18.883    0.203)   25.247
   0.917   (  13.805   15.461   -0.391)   20.731
   1.511   (  15.278   21.757    1.344)   26.619
   1.573   (  14.654   21.085   -1.275)   25.709
   2.224   (  24.610   26.591    0.445)   36.235
   2.243   (  24.762   25.319   -0.365)   35.417
   5.311   (   2.661    3.693    0.026)    4.552
   5.312   (   2.768    3.785   -0.025)    4.689
   5.614   (   9.031    9.978    0.720)   13.477
   5.648   (   9.453   10.437   -0.694)   14.099
   6.929   (  -6.075   -7.254    0.120)    9.462
   6.935   (  -6.321   -7.563   -0.154)    9.858
   7.324   (  -0.312   -0.068    0.145)    0.351
   7.331   (  -0.379   -0.156   -0.147)    0.436
   7.696   (   3.198    3.570    0.502)    4.819
   7.721   (   3.152    3.538   -0.563)    4.772
   8.891   (  -7.040   -8.089    0.153)   10.724
   8.898   (  -7.156   -8.229   -0.136)   10.906
======================= Grid point 365 (50/111) =======================
q-point: ( 0.17  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.379   (  20.050   18.231    0.819)   27.112
   1.420   (  21.284   19.413   -0.903)   28.822
   1.955   (  12.913   19.377    1.033)   23.308
   2.003   (  12.387   19.126   -1.001)   22.809
   2.848   (  22.528   18.054    0.530)   28.875
   2.872   (  23.307   18.056   -0.476)   29.487
   5.399   (   2.975    4.348    0.100)    5.269
   5.404   (   3.161    4.449   -0.098)    5.458
   5.875   (  10.262    8.969    0.923)   13.660
   5.919   (  10.632    9.280   -0.929)   14.143
   6.758   (  -5.854   -6.147    0.127)    8.489
   6.765   (  -5.421   -6.057   -0.129)    8.129
   7.313   (  -0.840   -0.085    0.125)    0.853
   7.319   (  -0.827   -0.094   -0.126)    0.842
   7.779   (   3.011    2.524    0.523)    3.963
   7.804   (   2.995    2.573   -0.527)    3.983
   8.664   (  -9.595   -9.077    0.110)   13.209
   8.669   (  -9.576   -9.112   -0.103)   13.219
======================= Grid point 366 (51/111) =======================
q-point: ( 0.22  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.922   (  22.285   17.719    1.291)   28.500
   1.984   (  22.793   18.141   -1.315)   29.161
   2.324   (   9.436   16.972    0.822)   19.436
   2.362   (   9.021   16.874   -0.805)   19.151
   3.328   (  17.094   10.285    0.794)   19.966
   3.365   (  17.576   10.606   -0.770)   20.542
   5.492   (   2.614    4.743    0.194)    5.420
   5.502   (   2.842    4.817   -0.193)    5.596
   6.124   (   9.364    6.913    1.103)   11.691
   6.176   (   9.714    7.159   -1.115)   12.118
   6.626   (  -3.386   -4.023    0.339)    5.269
   6.642   (  -2.878   -3.904   -0.332)    4.861
   7.293   (  -1.435    0.002    0.147)    1.443
   7.300   (  -1.356    0.052   -0.148)    1.365
   7.846   (   2.569    1.487    0.528)    3.015
   7.871   (   2.502    1.522   -0.517)    2.973
   8.397   ( -10.939   -9.195    0.135)   14.291
   8.403   ( -10.835   -9.204   -0.127)   14.217
======================= Grid point 367 (52/111) =======================
q-point: ( 0.28  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.462   (  21.595   15.511    1.486)   26.630
   2.533   (  21.847   15.654   -1.487)   26.917
   2.607   (   5.764   14.464    0.686)   15.585
   2.639   (   5.533   14.370   -0.672)   15.413
   3.622   (   9.749    3.269    0.924)   10.324
   3.666   (   9.811    3.373   -0.936)   10.417
   5.578   (   1.741    4.915    0.285)    5.222
   5.591   (   1.950    4.948   -0.284)    5.326
   6.314   (   6.387    3.977    1.279)    7.632
   6.375   (   6.748    4.186   -1.291)    8.045
   6.575   (   0.513   -1.376    0.503)    1.552
   6.598   (   0.792   -1.432   -0.482)    1.706
   7.266   (  -1.899    0.271    0.200)    1.928
   7.275   (  -1.781    0.378   -0.199)    1.831
   7.894   (   1.869    0.376    0.487)    1.967
   7.916   (   1.708    0.320   -0.476)    1.802
   8.127   ( -10.316   -8.106    0.159)   13.121
   8.134   ( -10.281   -8.180   -0.144)   13.139
======================= Grid point 368 (53/111) =======================
q-point: ( 0.33  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.805   (   2.405   11.858    0.613)   12.115
   2.834   (   2.400   11.649   -0.603)   11.909
   2.944   (  18.577   12.219    1.579)   22.291
   3.019   (  18.761   12.220   -1.572)   22.445
   3.735   (   3.307   -1.646    0.853)    3.791
   3.776   (   3.048   -1.925   -0.885)    3.712
   5.646   (   0.575    4.960    0.352)    5.006
   5.663   (   0.727    4.955   -0.352)    5.021
   6.407   (   1.812    0.555    1.412)    2.363
   6.476   (   2.101    0.642   -1.474)    2.645
   6.610   (   3.807    0.303    0.536)    3.857
   6.634   (   3.747    0.059   -0.454)    3.774
   7.237   (  -2.076    0.848    0.271)    2.259
   7.250   (  -1.958    0.985   -0.268)    2.208
   7.905   (  -7.197   -5.702    0.110)    9.183
   7.911   (  -6.390   -5.280   -0.129)    8.290
   7.918   (   0.180   -1.209    0.401)    1.287
   7.936   (   0.791   -0.830   -0.357)    1.201
======================= Grid point 369 (54/111) =======================
q-point: ( 0.39  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.925   (  -0.345    9.244    0.586)    9.269
   2.952   (  -0.152    8.834   -0.573)    8.854
   3.317   (  13.156    7.900    1.622)   15.431
   3.394   (  13.330    7.782   -1.631)   15.521
   3.732   (   0.125   -3.301    0.670)    3.370
   3.764   (  -0.098   -3.865   -0.691)    3.928
   5.694   (  -0.700    4.943    0.392)    5.008
   5.712   (  -0.607    4.924   -0.393)    4.977
   6.396   (  -1.740   -1.779    1.451)    2.880
   6.467   (  -1.648   -1.938   -1.562)    2.985
   6.688   (   4.623   -0.137    0.442)    4.646
   6.706   (   4.389   -0.434   -0.304)    4.421
   7.215   (  -1.982    1.761    0.335)    2.672
   7.231   (  -1.910    1.884   -0.330)    2.703
   7.774   (  -3.165   -3.128    0.098)    4.451
   7.779   (  -3.175   -3.169   -0.099)    4.487
   7.922   (   0.625   -1.240    0.254)    1.412
   7.933   (   0.333   -1.463   -0.238)    1.520
======================= Grid point 370 (55/111) =======================
q-point: ( 0.44  0.06  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.979   (  -2.540    7.028    0.578)    7.495
   3.006   (  -2.218    6.448   -0.563)    6.842
   3.527   (   5.020    2.506    1.616)    5.839
   3.605   (   5.169    2.264   -1.698)    5.893
   3.702   (   0.276   -2.491    0.551)    2.566
   3.727   (   0.408   -3.096   -0.496)    3.161
   5.719   (  -1.984    4.743    0.412)    5.157
   5.738   (  -1.946    4.729   -0.414)    5.130
   6.347   (  -1.177   -1.484    1.415)    2.364
   6.417   (  -1.110   -1.775   -1.531)    2.594
   6.746   (   3.047   -1.999    0.323)    3.659
   6.758   (   2.926   -2.172   -0.175)    3.649
   7.205   (  -1.781    2.616    0.371)    3.186
   7.222   (  -1.780    2.698   -0.363)    3.252
   7.725   (  -0.354   -1.494    0.088)    1.538
   7.730   (  -0.358   -1.495   -0.099)    1.541
   7.919   (   0.690   -1.340    0.107)    1.511
   7.924   (   0.451   -1.534   -0.094)    1.602
======================= Grid point 382 (56/111) =======================
q-point: ( 0.11  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.296   (  13.048   22.599    0.701)   26.105
   1.332   (  14.037   24.313   -0.803)   28.086
   1.953   (  12.878   22.305    1.126)   25.780
   2.005   (  12.653   21.916   -1.073)   25.329
   2.736   (  14.069   24.368    0.224)   28.139
   2.746   (  14.723   25.501   -0.180)   29.447
   5.397   (   2.891    5.007    0.075)    5.782
   5.401   (   2.968    5.141   -0.073)    5.937
   5.829   (   6.687   11.583    0.888)   13.404
   5.872   (   6.930   12.003   -0.890)   13.889
   6.779   (  -4.574   -7.922    0.011)    9.147
   6.780   (  -4.354   -7.542   -0.014)    8.708
   7.321   (  -0.145   -0.251    0.125)    0.316
   7.327   (  -0.143   -0.248   -0.127)    0.314
   7.765   (   1.928    3.340    0.531)    3.893
   7.790   (   1.953    3.383   -0.539)    3.944
   8.701   (  -6.414  -11.109    0.103)   12.828
   8.706   (  -6.451  -11.173   -0.095)   12.902
======================= Grid point 383 (57/111) =======================
q-point: ( 0.17  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.777   (  17.162   21.858    1.213)   27.817
   1.836   (  17.641   22.493   -1.250)   28.612
   2.362   (   9.787   21.257    0.969)   23.422
   2.408   (   9.428   21.160   -0.931)   23.184
   3.170   (  14.086   14.203    0.677)   20.015
   3.201   (  14.774   14.937   -0.658)   21.019
   5.501   (   2.882    5.909    0.139)    6.576
   5.508   (   3.032    6.002   -0.138)    6.726
   6.055   (   7.521    9.040    1.050)   11.806
   6.105   (   7.787    9.355   -1.059)   12.218
   6.639   (  -3.666   -5.717    0.197)    6.794
   6.648   (  -3.264   -5.490   -0.192)    6.390
   7.313   (  -0.662    0.028    0.147)    0.679
   7.320   (  -0.622    0.151   -0.149)    0.657
   7.824   (   1.899    1.976    0.555)    2.796
   7.850   (   1.879    2.033   -0.547)    2.822
   8.463   (  -8.496  -11.215    0.099)   14.070
   8.467   (  -8.440  -11.237   -0.091)   14.054
======================= Grid point 384 (58/111) =======================
q-point: ( 0.22  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.286   (  18.586   19.113    1.430)   26.698
   2.354   (  18.816   19.315   -1.436)   27.003
   2.693   (   5.249   19.873    0.811)   20.570
   2.731   (   4.795   19.869   -0.777)   20.454
   3.477   (  10.916    4.643    0.933)   11.899
   3.522   (  11.051    5.042   -0.951)   12.185
   5.606   (   2.049    6.716    0.216)    7.025
   5.617   (   2.243    6.758   -0.215)    7.124
   6.249   (   6.459    5.626    1.211)    8.651
   6.306   (   6.750    5.882   -1.215)    9.035
   6.555   (  -0.854   -3.010    0.372)    3.151
   6.573   (  -0.500   -2.976   -0.359)    3.039
   7.299   (  -1.273    0.614    0.201)    1.428
   7.309   (  -1.212    0.814   -0.201)    1.474
   7.866   (   1.623    0.452    0.534)    1.768
   7.891   (   1.486    0.423   -0.522)    1.631
   8.204   (  -9.168  -10.220    0.119)   13.730
   8.210   (  -9.096  -10.302   -0.107)   13.743
======================= Grid point 385 (59/111) =======================
q-point: ( 0.28  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.767   (  17.628   15.325    1.514)   23.407
   2.839   (  17.767   15.292   -1.509)   23.490
   2.928   (   0.863   17.330    0.646)   17.363
   2.958   (   0.472   17.144   -0.602)   17.161
   3.633   (   6.289   -1.896    0.881)    6.627
   3.676   (   5.939   -2.204   -0.940)    6.404
   5.699   (   0.714    7.372    0.287)    7.412
   5.713   (   0.895    7.364   -0.287)    7.424
   6.372   (   3.447    1.730    1.340)    4.083
   6.436   (   3.708    1.855   -1.381)    4.370
   6.551   (   2.467   -1.072    0.463)    2.730
   6.572   (   2.584   -1.270   -0.402)    2.908
   7.284   (  -1.755    1.610    0.270)    2.397
   7.297   (  -1.677    1.818   -0.268)    2.488
   7.888   (   1.237   -1.020    0.440)    1.662
   7.908   (   0.985   -1.186   -0.440)    1.603
   7.967   (  -7.712   -8.012    0.114)   11.121
   7.972   (  -7.707   -8.174   -0.090)   11.235
======================= Grid point 386 (60/111) =======================
q-point: ( 0.33  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.063   (  -2.297   13.008    0.465)   13.217
   3.084   (  -2.486   12.399   -0.408)   12.652
   3.171   (  14.338   10.849    1.530)   18.044
   3.244   (  14.474   10.604   -1.531)   18.008
   3.676   (   2.499   -3.199    0.572)    4.099
   3.705   (   2.131   -4.383   -0.644)    4.916
   5.773   (  -0.839    7.864    0.340)    7.916
   5.790   (  -0.690    7.820   -0.340)    7.858
   6.397   (  -0.312   -1.564    1.367)    2.101
   6.465   (  -0.202   -1.748   -1.499)    2.312
   6.602   (   4.333   -1.315    0.445)    4.550
   6.620   (   4.216   -1.630   -0.290)    4.529
   7.274   (  -1.991    2.952    0.338)    3.576
   7.290   (  -1.919    3.112   -0.333)    3.672
   7.798   (  -4.355   -5.234    0.084)    6.809
   7.801   (  -4.317   -5.298   -0.083)    6.834
   7.891   (   0.978   -2.061    0.308)    2.302
   7.905   (   0.695   -2.311   -0.291)    2.431
======================= Grid point 387 (61/111) =======================
q-point: ( 0.39  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.115   (  -3.366    8.094    0.275)    8.771
   3.127   (  -3.253    7.054   -0.222)    7.771
   3.449   (   8.358    5.678    1.493)   10.214
   3.521   (   8.535    5.229   -1.551)   10.129
   3.677   (   0.210   -0.593    0.248)    0.676
   3.689   (   0.339   -2.156   -0.253)    2.197
   5.824   (  -2.514    8.066    0.379)    8.457
   5.842   (  -2.394    8.018   -0.378)    8.376
   6.352   (  -1.680   -2.568    1.292)    3.330
   6.417   (  -1.658   -3.081   -1.448)    3.786
   6.656   (   4.221   -3.165    0.340)    5.287
   6.668   (   4.074   -3.418   -0.164)    5.320
   7.274   (  -2.187    4.076    0.385)    4.641
   7.292   (  -2.157    4.180   -0.378)    4.719
   7.714   (  -1.000   -3.113    0.088)    3.271
   7.718   (  -0.980   -3.116   -0.094)    3.267
   7.885   (   1.067   -2.446    0.153)    2.673
   7.892   (   0.802   -2.692   -0.141)    2.813
======================= Grid point 388 (62/111) =======================
q-point: ( 0.44  0.11  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.124   (  -3.538    6.127    0.158)    7.077
   3.130   (  -2.666    4.618   -0.114)    5.334
   3.551   (  -0.108    0.188    1.433)    1.449
   3.621   (   0.181   -0.313   -1.611)    1.651
   3.676   (  -0.630    1.092    0.146)    1.269
   3.682   (   0.190   -0.329   -0.017)    0.381
   5.843   (  -4.350    7.535    0.398)    8.710
   5.862   (  -4.311    7.467   -0.398)    8.631
   6.320   (   0.688   -1.191    1.227)    1.843
   6.382   (   0.981   -1.700   -1.382)    2.400
   6.677   (   2.777   -4.810    0.282)    5.561
   6.686   (   2.826   -4.894   -0.108)    5.652
   7.276   (  -2.513    4.353    0.402)    5.043
   7.295   (  -2.548    4.414   -0.394)    5.112
   7.692   (   1.110   -1.922    0.088)    2.221
   7.696   (   1.111   -1.925   -0.100)    2.225
   7.882   (   1.368   -2.370    0.041)    2.737
   7.884   (   1.410   -2.442   -0.030)    2.820
======================= Grid point 402 (63/111) =======================
q-point: ( 0.17  0.17  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.220   (  13.227   22.909    1.407)   26.490
   2.287   (  13.380   23.174   -1.415)   26.796
   2.760   (  10.778   18.668    0.933)   21.576
   2.803   (  10.717   18.562   -0.896)   21.453
   3.393   (   4.887    8.465    0.971)    9.822
   3.440   (   5.237    9.071   -0.987)   10.520
   5.624   (   3.730    6.461    0.182)    7.463
   5.633   (   3.802    6.585   -0.183)    7.606
   6.219   (   4.299    7.445    1.187)    8.679
   6.276   (   4.486    7.771   -1.188)    9.051
   6.547   (  -1.913   -3.313    0.302)    3.837
   6.561   (  -1.776   -3.077   -0.292)    3.565
   7.314   (   0.066    0.114    0.216)    0.253
   7.325   (   0.174    0.301   -0.216)    0.409
   7.854   (   0.576    0.998    0.558)    1.281
   7.880   (   0.536    0.928   -0.547)    1.203
   8.233   (  -6.921  -11.987    0.092)   13.842
   8.237   (  -6.938  -12.016   -0.080)   13.875
======================= Grid point 403 (64/111) =======================
q-point: ( 0.22  0.17  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.661   (  14.172   18.821    1.473)   23.606
   2.731   (  14.222   18.801   -1.467)   23.620
   3.074   (   4.235   18.063    0.859)   18.572
   3.113   (   3.645   18.331   -0.796)   18.707
   3.509   (   5.175   -0.975    0.994)    5.359
   3.558   (   4.607   -1.085   -1.068)    4.852
   5.749   (   2.631    7.612    0.237)    8.057
   5.760   (   2.746    7.689   -0.239)    8.168
   6.336   (   3.111    3.122    1.302)    4.596
   6.399   (   3.286    3.325   -1.330)    4.860
   6.513   (   0.457   -1.313    0.385)    1.443
   6.530   (   0.556   -1.363   -0.334)    1.510
   7.319   (  -0.364    1.389    0.292)    1.466
   7.333   (  -0.293    1.598   -0.289)    1.651
   7.863   (   0.411   -0.791    0.487)    1.016
   7.886   (   0.220   -0.987   -0.481)    1.120
   8.005   (  -6.895   -9.867    0.072)   12.038
   8.008   (  -6.901  -10.024   -0.055)   12.170
======================= Grid point 404 (65/111) =======================
q-point: ( 0.28  0.17  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.058   (  13.035   14.156    1.444)   19.297
   3.126   (  13.039   13.860   -1.450)   19.085
   3.227   (  -2.883    9.786    0.730)   10.228
   3.258   (  -3.699   10.917   -0.570)   11.541
   3.592   (   5.654    0.647    0.372)    5.703
   3.614   (   4.802   -1.403   -0.540)    5.032
   5.859   (   0.853    8.634    0.288)    8.681
   5.873   (   0.952    8.662   -0.289)    8.719
   6.380   (   0.688   -0.902    1.298)    1.724
   6.444   (   0.717   -1.023   -1.465)    1.925
   6.527   (   2.466   -1.522    0.423)    2.929
   6.544   (   2.424   -1.752   -0.231)    3.000
   7.331   (  -0.844    3.125    0.352)    3.256
   7.347   (  -0.802    3.286   -0.347)    3.400
   7.817   (  -5.018   -7.109    0.046)    8.702
   7.820   (  -4.926   -7.155   -0.049)    8.687
   7.855   (   0.489   -2.284    0.356)    2.362
   7.871   (   0.215   -2.567   -0.336)    2.598
======================= Grid point 405 (66/111) =======================
q-point: ( 0.33  0.17  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.203   (  -2.955   -1.143    0.602)    3.225
   3.232   (  -3.988    0.011   -0.622)    4.036
   3.368   (   9.442    9.209    1.396)   13.263
   3.433   (   9.503    8.635   -1.339)   12.910
   3.692   (   1.350    6.254    0.243)    6.403
   3.704   (   1.099    7.665   -0.290)    7.749
   5.945   (  -1.176    9.243    0.321)    9.323
   5.960   (  -1.077    9.190   -0.318)    9.258
   6.351   (  -1.327   -2.971    1.165)    3.456
   6.411   (  -1.480   -3.569   -1.397)    4.109
   6.555   (   3.392   -3.522    0.362)    4.904
   6.567   (   3.265   -3.780   -0.118)    4.996
   7.349   (  -1.473    4.545    0.393)    4.794
   7.368   (  -1.437    4.639   -0.385)    4.872
   7.700   (  -1.988   -4.688    0.077)    5.093
   7.704   (  -1.935   -4.672   -0.079)    5.058
   7.840   (   0.906   -3.018    0.190)    3.157
   7.848   (   0.648   -3.304   -0.181)    3.372
======================= Grid point 406 (67/111) =======================
q-point: ( 0.39  0.17  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.149   (   0.320   -3.882    0.410)    3.917
   3.170   (  -0.426   -3.322   -0.467)    3.381
   3.544   (   3.120    3.996    1.294)    5.232
   3.604   (   3.366    3.272   -1.264)    4.862
   3.752   (  -3.165    8.775    0.418)    9.338
   3.771   (  -3.472    9.574   -0.388)   10.192
   5.996   (  -3.528    9.053    0.348)    9.723
   6.012   (  -3.384    8.897   -0.339)    9.525
   6.301   (  -0.581   -2.411    0.963)    2.661
   6.352   (  -0.474   -3.318   -1.198)    3.559
   6.570   (   3.516   -5.319    0.283)    6.382
   6.578   (   3.476   -5.468   -0.054)    6.480
   7.365   (  -2.323    4.927    0.417)    5.463
   7.384   (  -2.311    4.971   -0.409)    5.497
   7.651   (   0.564   -3.216    0.083)    3.266
   7.656   (   0.581   -3.234   -0.093)    3.287
   7.830   (   1.500   -3.047    0.058)    3.397
   7.832   (   1.297   -3.235   -0.050)    3.486
======================= Grid point 421 (68/111) =======================
q-point: ( 0.22  0.22  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.015   (   9.806   16.984    1.412)   19.662
   3.082   (   9.658   16.729   -1.415)   19.369
   3.313   (   0.809    1.401    1.259)    2.050
   3.373   (   1.235    2.139   -1.237)    2.763
   3.522   (   4.459    7.724    0.635)    8.941
   3.553   (   3.746    6.487   -0.667)    7.521
   5.897   (   4.201    7.277    0.278)    8.407
   5.910   (   4.267    7.390   -0.281)    8.538
   6.367   (  -0.002   -0.004    1.278)    1.278
   6.431   (  -0.039   -0.068   -1.467)    1.469
   6.498   (  -0.282   -0.488    0.410)    0.697
   6.513   (  -0.358   -0.620   -0.195)    0.742
   7.355   (   1.293    2.240    0.364)    2.612
   7.372   (   1.366    2.366   -0.358)    2.755
   7.826   (  -4.690   -8.123    0.001)    9.380
   7.826   (  -4.649   -8.052   -0.016)    9.297
   7.838   (  -1.026   -1.778    0.385)    2.089
   7.855   (  -1.176   -2.037   -0.354)    2.379
======================= Grid point 422 (69/111) =======================
q-point: ( 0.28  0.22  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.214   (  -4.032   -7.168    0.922)    8.275
   3.264   (  -4.813   -7.541   -1.142)    9.019
   3.318   (   8.441   11.928    1.371)   14.677
   3.378   (   8.261   11.361   -1.153)   14.094
   3.738   (   5.104   12.475    0.501)   13.488
   3.762   (   4.933   12.641   -0.507)   13.579
   6.027   (   2.001    8.160    0.316)    8.407
   6.042   (   1.987    8.282   -0.320)    8.524
   6.342   (  -1.099   -2.731    1.099)    3.143
   6.400   (  -1.401   -3.264   -1.433)    3.830
   6.486   (   0.675   -2.690    0.393)    2.801
   6.498   (   0.549   -2.911   -0.043)    2.963
   7.402   (   0.516    3.984    0.401)    4.037
   7.421   (   0.538    4.056   -0.392)    4.110
   7.691   (  -2.970   -5.642    0.064)    6.376
   7.694   (  -2.911   -5.583   -0.065)    6.296
   7.804   (  -0.316   -2.766    0.215)    2.792
   7.814   (  -0.558   -3.104   -0.208)    3.161
======================= Grid point 423 (70/111) =======================
q-point: ( 0.33  0.22  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.088   (  -1.137   -7.970    0.608)    8.074
   3.120   (  -1.813   -8.508   -0.732)    8.729
   3.530   (   4.806    7.182    1.107)    8.712
   3.579   (   4.774    6.216   -0.979)    7.899
   3.902   (  -0.494   12.382    0.487)   12.402
   3.925   (  -0.497   12.324   -0.483)   12.344
   6.119   (  -0.557    7.956    0.296)    7.981
   6.133   (  -0.783    7.860   -0.278)    7.904
   6.292   (  -1.123   -2.428    0.732)    2.774
   6.335   (  -1.625   -3.707   -1.099)    4.194
   6.470   (   1.881   -4.824    0.317)    5.187
   6.478   (   1.794   -5.011    0.024)    5.323
   7.442   (  -0.972    4.647    0.417)    4.766
   7.462   (  -0.957    4.667   -0.408)    4.781
   7.615   (  -0.577   -3.976    0.072)    4.018
   7.618   (  -0.564   -4.011   -0.081)    4.052
   7.778   (   0.786   -3.100    0.080)    3.199
   7.782   (   0.620   -3.315   -0.075)    3.374
======================= Grid point 424 (71/111) =======================
q-point: ( 0.39  0.22  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.037   (   3.466   -6.003    0.389)    6.942
   3.059   (   3.604   -6.242   -0.509)    7.225
   3.608   (  -1.507    2.611    0.984)    3.171
   3.652   (  -0.971    1.681   -0.853)    2.120
   3.959   (  -6.017   10.421    0.483)   12.043
   3.982   (  -5.952   10.309   -0.475)   11.914
   6.158   (  -3.859    6.685    0.186)    7.721
   6.164   (  -3.118    5.401   -0.073)    6.237
   6.271   (  -0.046    0.080    0.314)    0.327
   6.296   (   1.146   -1.985   -0.763)    2.415
   6.461   (   3.068   -5.313    0.274)    6.141
   6.467   (   3.129   -5.420    0.044)    6.259
   7.457   (  -2.456    4.254    0.421)    4.930
   7.477   (  -2.446    4.237   -0.414)    4.910
   7.591   (   1.626   -2.816    0.068)    3.252
   7.594   (   1.659   -2.873   -0.078)    3.318
   7.770   (   1.684   -2.917    0.014)    3.369
   7.771   (   1.669   -2.891   -0.009)    3.338
======================= Grid point 441 (72/111) =======================
q-point: ( 0.28  0.28  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 342
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.056   (  -4.380   -7.586    0.670)    8.785
   3.091   (  -4.902   -8.490   -0.804)    9.836
   3.524   (   4.989    8.641    1.012)   10.029
   3.569   (   4.463    7.730   -0.866)    8.967
   3.972   (   5.848   10.130    0.469)   11.706
   3.994   (   5.761    9.979   -0.470)   11.532
   6.162   (   2.859    4.952    0.379)    5.731
   6.181   (   3.054    5.289   -0.429)    6.122
   6.290   (  -1.023   -1.773    0.726)    2.172
   6.332   (  -1.714   -2.969   -1.127)    3.609
   6.428   (  -1.735   -3.004    0.352)    3.487
   6.435   (  -1.851   -3.205    0.088)    3.702
   7.476   (   1.887    3.268    0.418)    3.796
   7.496   (   1.892    3.277   -0.407)    3.806
   7.599   (  -2.134   -3.697    0.065)    4.269
   7.602   (  -2.154   -3.730   -0.074)    4.308
   7.754   (  -1.249   -2.163    0.088)    2.499
   7.758   (  -1.392   -2.411   -0.083)    2.786
======================= Grid point 442 (73/111) =======================
q-point: ( 0.33  0.28  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 648
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.949   (  -0.189   -5.358    0.303)    5.370
   2.966   (  -0.553   -6.210   -0.431)    6.249
   3.645   (   1.334    4.375    0.649)    4.620
   3.672   (   1.277    3.082   -0.490)    3.372
   4.112   (  -0.326    8.063    0.423)    8.080
   4.132   (  -0.307    7.908   -0.423)    7.925
   6.200   (  -0.052   -0.108    0.266)    0.292
   6.221   (  -0.333   -0.349   -0.618)    0.784
   6.297   (  -0.211    3.052    0.330)    3.077
   6.316   (  -0.779    3.093   -0.463)    3.223
   6.382   (   0.546   -3.647    0.297)    3.699
   6.386   (   0.503   -3.839    0.152)    3.875
   7.521   (  -0.440    3.158    0.412)    3.215
   7.541   (  -0.448    3.062   -0.409)    3.122
   7.549   (  -0.058   -2.440    0.044)    2.441
   7.551   (  -0.070   -2.591   -0.051)    2.593
   7.724   (   0.401   -2.233    0.020)    2.269
   7.725   (   0.367   -2.311   -0.016)    2.340
======================= Grid point 460 (74/111) =======================
q-point: ( 0.33  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 128
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.898   (  -0.226   -0.391   -0.064)    0.456
   2.898   (   0.226    0.391   -0.064)    0.456
   3.696   (  -0.356   -0.616    0.086)    0.716
   3.696   (   0.356    0.616    0.086)    0.716
   4.194   (  -0.000   -0.000    0.387)    0.387
   4.213   (  -0.000   -0.000   -0.389)    0.389
   6.198   (  -0.328   -0.568   -0.274)    0.711
   6.198   (   0.328    0.568   -0.274)    0.711
   6.339   (  -0.000   -0.000    0.400)    0.400
   6.344   (  -0.056   -0.096    0.254)    0.277
   6.344   (   0.056    0.096    0.254)    0.277
   6.358   (  -0.000   -0.000   -0.409)    0.409
   7.524   (  -0.031   -0.053   -0.009)    0.062
   7.524   (   0.031    0.053   -0.009)    0.062
   7.555   (  -0.000   -0.000    0.415)    0.415
   7.575   (  -0.000   -0.000   -0.404)    0.404
   7.700   (  -0.002   -0.003    0.003)    0.005
   7.700   (   0.002    0.003    0.003)    0.005
======================= Grid point 686 (75/111) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 74
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.356   (   0.000    0.000   -7.492)    7.492
   0.356   (   0.000    0.000   -7.492)    7.492
   0.356   (   0.000    0.000    7.492)    7.492
   0.356   (   0.000    0.000    7.492)    7.492
   0.460   (  -0.000    0.000   -9.710)    9.710
   0.460   (   0.000   -0.000    9.710)    9.710
   5.183   (   0.000   -0.000   -0.097)    0.097
   5.183   (   0.000    0.000   -0.097)    0.097
   5.183   (   0.000    0.000    0.097)    0.097
   5.183   (   0.000    0.000    0.097)    0.097
   7.333   (   0.000    0.000   -0.384)    0.384
   7.333   (  -0.000   -0.000   -0.384)    0.384
   7.333   (   0.000   -0.000    0.384)    0.384
   7.333   (  -0.000   -0.000    0.384)    0.384
   7.428   (   0.000    0.000   -0.729)    0.729
   7.428   (   0.000    0.000    0.729)    0.729
   9.200   (   0.000   -0.000   -0.298)    0.298
   9.200   (  -0.000    0.000    0.298)    0.298
======================= Grid point 688 (76/111) =======================
q-point: ( 0.06  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.468   (   2.898    1.673   -7.848)    8.531
   0.468   (   2.898    1.673    7.848)    8.531
   0.684   (  21.626   12.486   -3.957)   25.283
   0.684   (  21.626   12.486    3.957)   25.283
   0.960   (  35.605   20.556   -2.625)   41.197
   0.960   (  35.605   20.556    2.625)   41.197
   5.203   (   1.673    0.966   -0.084)    1.933
   5.203   (   1.673    0.966    0.084)    1.933
   5.262   (   6.651    3.840   -0.290)    7.685
   5.262   (   6.651    3.840    0.290)    7.685
   7.227   (  -6.787   -3.919   -0.441)    7.849
   7.227   (  -6.787   -3.919    0.441)    7.849
   7.332   (  -0.054   -0.031   -0.347)    0.352
   7.332   (  -0.054   -0.031    0.347)    0.352
   7.526   (   5.563    3.212   -0.695)    6.461
   7.526   (   5.563    3.212    0.695)    6.461
   9.151   (  -4.109   -2.373   -0.380)    4.760
   9.151   (  -4.109   -2.373    0.380)    4.760
======================= Grid point 690 (77/111) =======================
q-point: ( 0.11  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.657   (  14.270    8.239   -2.750)   16.706
   0.657   (  14.270    8.239    2.750)   16.706
   1.191   (  22.211   12.824   -2.326)   25.752
   1.191   (  22.211   12.824    2.326)   25.752
   1.767   (  34.673   20.018   -1.186)   40.054
   1.767   (  34.673   20.018    1.186)   40.054
   5.257   (   3.039    1.754   -0.039)    3.509
   5.257   (   3.039    1.754    0.039)    3.509
   5.466   (  10.850    6.264   -0.740)   12.551
   5.466   (  10.850    6.264    0.740)   12.551
   7.054   (  -8.166   -4.715   -0.339)    9.436
   7.054   (  -8.166   -4.715    0.339)    9.436
   7.329   (  -0.232   -0.134   -0.261)    0.374
   7.329   (  -0.232   -0.134    0.261)    0.374
   7.644   (   4.742    2.738   -0.774)    5.531
   7.644   (   4.742    2.738    0.774)    5.531
   9.019   (  -7.525   -4.344   -0.284)    8.693
   9.019   (  -7.525   -4.344    0.284)    8.693
======================= Grid point 692 (78/111) =======================
q-point: ( 0.17  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.084   (  22.341   12.899   -0.389)   25.800
   1.084   (  22.341   12.899    0.389)   25.800
   1.669   (  19.745   11.400   -1.643)   22.859
   1.669   (  19.745   11.400    1.643)   22.859
   2.505   (  29.984   17.312   -0.786)   34.632
   2.505   (  29.984   17.312    0.786)   34.632
   5.336   (   3.865    2.231   -0.080)    4.463
   5.336   (   3.865    2.231    0.080)    4.463
   5.734   (  12.417    7.169   -1.130)   14.382
   5.734   (  12.417    7.169    1.130)   14.382
   6.869   (  -7.907   -4.565   -0.160)    9.132
   6.869   (  -7.907   -4.565    0.160)    9.132
   7.321   (  -0.575   -0.332   -0.194)    0.692
   7.321   (  -0.575   -0.332    0.194)    0.692
   7.742   (   3.981    2.299   -0.733)    4.655
   7.742   (   3.981    2.299    0.733)    4.655
   8.815   ( -10.286   -5.938   -0.176)   11.878
   8.815   ( -10.286   -5.938    0.176)   11.878
======================= Grid point 694 (79/111) =======================
q-point: ( 0.22  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.630   (  25.213   14.557   -1.544)   29.154
   1.630   (  25.213   14.557    1.544)   29.154
   2.079   (  16.297    9.409   -1.253)   18.860
   2.079   (  16.297    9.409    1.253)   18.860
   3.108   (  22.919   13.233   -0.955)   26.482
   3.108   (  22.919   13.233    0.955)   26.482
   5.427   (   4.050    2.339   -0.214)    4.682
   5.427   (   4.050    2.339    0.214)    4.682
   6.012   (  11.741    6.778   -1.428)   13.632
   6.012   (  11.741    6.778    1.428)   13.632
   6.709   (  -5.909   -3.411   -0.364)    6.833
   6.709   (  -5.909   -3.411    0.364)    6.833
   7.303   (  -1.002   -0.578   -0.196)    1.173
   7.303   (  -1.002   -0.578    0.196)    1.173
   7.824   (   3.250    1.876   -0.726)    3.822
   7.824   (   3.250    1.876    0.726)    3.822
   8.560   ( -12.042   -6.953   -0.174)   13.906
   8.560   ( -12.042   -6.953    0.174)   13.906
======================= Grid point 696 (80/111) =======================
q-point: ( 0.28  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.201   (  24.759   14.294   -1.985)   28.658
   2.201   (  24.759   14.294    1.985)   28.658
   2.406   (  12.551    7.246   -1.020)   14.528
   2.406   (  12.551    7.246    1.020)   14.528
   3.529   (  14.032    8.101   -1.229)   16.249
   3.529   (  14.032    8.101    1.229)   16.249
   5.515   (   3.656    2.111   -0.359)    4.236
   5.515   (   3.656    2.111    0.359)    4.236
   6.251   (   9.085    5.245   -1.688)   10.625
   6.251   (   9.085    5.245    1.688)   10.625
   6.614   (  -2.263   -1.307   -0.656)    2.694
   6.614   (  -2.263   -1.307    0.656)    2.694
   7.276   (  -1.347   -0.778   -0.244)    1.574
   7.276   (  -1.347   -0.778    0.244)    1.574
   7.888   (   2.286    1.320   -0.703)    2.732
   7.888   (   2.286    1.320    0.703)    2.732
   8.282   ( -12.223   -7.057   -0.225)   14.116
   8.282   ( -12.223   -7.057    0.225)   14.116
======================= Grid point 698 (81/111) =======================
q-point: ( 0.33  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.648   (   8.914    5.147   -0.918)   10.334
   2.648   (   8.914    5.147    0.918)   10.334
   2.734   (  21.950   12.673   -2.185)   25.440
   2.734   (  21.950   12.673    2.185)   25.440
   3.744   (   5.223    3.015   -1.308)    6.171
   3.744   (   5.223    3.015    1.308)    6.171
   5.590   (   2.865    1.654   -0.477)    3.343
   5.590   (   2.865    1.654    0.477)    3.343
   6.409   (   4.548    2.626   -1.939)    5.598
   6.409   (   4.548    2.626    1.939)    5.598
   6.610   (   1.740    1.004   -0.780)    2.155
   6.610   (   1.740    1.004    0.780)    2.155
   7.244   (  -1.447   -0.836   -0.323)    1.702
   7.244   (  -1.447   -0.836    0.323)    1.702
   7.926   (   1.111    0.641   -0.617)    1.423
   7.926   (   1.111    0.641    0.617)    1.423
   8.024   ( -10.068   -5.813   -0.224)   11.628
   8.024   ( -10.068   -5.813    0.224)   11.628
======================= Grid point 700 (82/111) =======================
q-point: ( 0.39  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.812   (   5.611    3.240   -0.914)    6.544
   2.812   (   5.611    3.240    0.914)    6.544
   3.180   (  17.026    9.830   -2.303)   19.794
   3.180   (  17.026    9.830    2.303)   19.794
   3.790   (  -0.487   -0.281   -1.160)    1.289
   3.790   (  -0.487   -0.281    1.160)    1.289
   5.644   (   1.902    1.098   -0.551)    2.265
   5.644   (   1.902    1.098    0.551)    2.265
   6.453   (  -0.477   -0.275   -2.127)    2.197
   6.453   (  -0.477   -0.275    2.127)    2.197
   6.678   (   3.757    2.169   -0.678)    4.391
   6.678   (   3.757    2.169    0.678)    4.391
   7.213   (  -1.219   -0.704   -0.412)    1.467
   7.213   (  -1.219   -0.704    0.412)    1.467
   7.841   (  -5.878   -3.394   -0.169)    6.789
   7.841   (  -5.878   -3.394    0.169)    6.789
   7.940   (   0.128    0.074   -0.460)    0.484
   7.940   (   0.128    0.074    0.460)    0.484
======================= Grid point 702 (83/111) =======================
q-point: ( 0.44  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.906   (   2.683    1.549   -0.951)    3.241
   2.906   (   2.683    1.549    0.951)    3.241
   3.487   (   9.626    5.558   -2.381)   11.367
   3.487   (   9.626    5.558    2.381)   11.367
   3.756   (  -1.740   -1.004   -0.956)    2.225
   3.756   (  -1.740   -1.004    0.956)    2.225
   5.676   (   0.931    0.538   -0.584)    1.224
   5.676   (   0.931    0.538    0.584)    1.224
   6.411   (  -2.376   -1.372   -2.163)    3.494
   6.411   (  -2.376   -1.372    2.163)    3.494
   6.757   (   2.838    1.639   -0.470)    3.311
   6.757   (   2.838    1.639    0.470)    3.311
   7.191   (  -0.695   -0.401   -0.480)    0.935
   7.191   (  -0.695   -0.401    0.480)    0.935
   7.754   (  -2.112   -1.219   -0.137)    2.442
   7.754   (  -2.112   -1.219    0.137)    2.442
   7.937   (  -0.247   -0.143   -0.254)    0.382
   7.937   (  -0.247   -0.143    0.254)    0.382
======================= Grid point 704 (84/111) =======================
q-point: (-0.50  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 141
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.936   (   0.000    0.000   -0.970)    0.970
   2.936   (   0.000    0.000    0.970)    0.970
   3.599   (   0.000    0.000   -2.402)    2.402
   3.599   (   0.000    0.000    2.402)    2.402
   3.734   (   0.000    0.000   -0.876)    0.876
   3.734   (   0.000    0.000    0.876)    0.876
   5.686   (   0.000    0.000   -0.593)    0.593
   5.686   (   0.000    0.000    0.593)    0.593
   6.379   (   0.000    0.000   -2.136)    2.136
   6.379   (   0.000    0.000    2.136)    2.136
   6.791   (   0.000    0.000   -0.356)    0.356
   6.791   (   0.000    0.000    0.356)    0.356
   7.183   (   0.000    0.000   -0.505)    0.505
   7.183   (   0.000    0.000    0.505)    0.505
   7.732   (   0.000    0.000   -0.131)    0.131
   7.732   (   0.000    0.000    0.131)    0.131
   7.934   (   0.000    0.000   -0.102)    0.102
   7.934   (   0.000    0.000    0.102)    0.102
======================= Grid point 726 (85/111) =======================
q-point: ( 0.06  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.579   (   6.342   10.985   -4.093)   13.328
   0.579   (   6.342   10.985    4.093)   13.328
   1.070   (  13.659   23.658   -2.571)   27.439
   1.070   (  13.659   23.658    2.571)   27.439
   1.545   (  20.059   34.743   -1.287)   40.138
   1.545   (  20.059   34.743    1.287)   40.138
   5.241   (   1.651    2.860   -0.037)    3.303
   5.241   (   1.651    2.860    0.037)    3.303
   5.401   (   5.701    9.874   -0.614)   11.418
   5.401   (   5.701    9.874    0.614)   11.418
   7.103   (  -4.628   -8.016   -0.401)    9.265
   7.103   (  -4.628   -8.016    0.401)    9.265
   7.331   (  -0.058   -0.101   -0.282)    0.305
   7.331   (  -0.058   -0.101    0.282)    0.305
   7.614   (   2.864    4.961   -0.783)    5.781
   7.614   (   2.864    4.961    0.783)    5.781
   9.062   (  -3.843   -6.657   -0.320)    7.694
   9.062   (  -3.843   -6.657    0.320)    7.694
======================= Grid point 728 (86/111) =======================
q-point: ( 0.11  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.908   (  15.295   17.183   -0.420)   23.008
   0.908   (  15.295   17.183    0.420)   23.008
   1.543   (  14.948   21.416   -1.851)   26.183
   1.543   (  14.948   21.416    1.851)   26.183
   2.234   (  24.684   25.940   -0.572)   35.812
   2.234   (  24.684   25.940    0.572)   35.812
   5.311   (   2.715    3.739   -0.036)    4.621
   5.311   (   2.715    3.739    0.036)    4.621
   5.631   (   9.232   10.196   -0.996)   13.791
   5.631   (   9.232   10.196    0.996)   13.791
   6.932   (  -6.184   -7.392   -0.190)    9.640
   6.932   (  -6.184   -7.392    0.190)    9.640
   7.327   (  -0.346   -0.113   -0.206)    0.418
   7.327   (  -0.346   -0.113    0.206)    0.418
   7.708   (   3.197    3.579   -0.749)    4.857
   7.708   (   3.197    3.579    0.749)    4.857
   8.895   (  -7.103   -8.164   -0.203)   10.823
   8.895   (  -7.103   -8.164    0.203)   10.823
======================= Grid point 730 (87/111) =======================
q-point: ( 0.17  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.399   (  20.688   18.839   -1.218)   28.007
   1.399   (  20.688   18.839    1.218)   28.007
   1.979   (  12.641   19.254   -1.438)   23.077
   1.979   (  12.641   19.254    1.438)   23.077
   2.861   (  22.915   18.045   -0.710)   29.175
   2.861   (  22.915   18.045    0.710)   29.175
   5.402   (   3.068    4.398   -0.140)    5.365
   5.402   (   3.068    4.398    0.140)    5.365
   5.897   (  10.443    9.122   -1.308)   13.928
   5.897   (  10.443    9.122    1.308)   13.928
   6.762   (  -5.634   -6.098   -0.181)    8.305
   6.762   (  -5.634   -6.098    0.181)    8.305
   7.316   (  -0.833   -0.090   -0.177)    0.856
   7.316   (  -0.833   -0.090    0.177)    0.856
   7.791   (   3.009    2.554   -0.742)    4.016
   7.791   (   3.009    2.554    0.742)    4.016
   8.666   (  -9.586   -9.095   -0.150)   13.215
   8.666   (  -9.586   -9.095    0.150)   13.215
======================= Grid point 732 (88/111) =======================
q-point: ( 0.22  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.953   (  22.548   17.936   -1.843)   28.871
   1.953   (  22.548   17.936    1.843)   28.871
   2.343   (   9.224   16.927   -1.150)   19.311
   2.343   (   9.224   16.927    1.150)   19.311
   3.346   (  17.329   10.435   -1.105)   20.259
   3.346   (  17.329   10.435    1.105)   20.259
   5.497   (   2.728    4.780   -0.274)    5.510
   5.497   (   2.728    4.780    0.274)    5.510
   6.150   (   9.539    7.037   -1.568)   11.957
   6.150   (   9.539    7.037    1.568)   11.957
   6.634   (  -3.130   -3.962   -0.474)    5.072
   6.634   (  -3.130   -3.962    0.474)    5.072
   7.297   (  -1.395    0.027   -0.209)    1.411
   7.297   (  -1.395    0.027    0.209)    1.411
   7.859   (   2.537    1.506   -0.738)    3.041
   7.859   (   2.537    1.506    0.738)    3.041
   8.400   ( -10.886   -9.199   -0.185)   14.253
   8.400   ( -10.886   -9.199    0.185)   14.253
======================= Grid point 734 (89/111) =======================
q-point: ( 0.28  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.498   (  21.724   15.584   -2.100)   26.818
   2.498   (  21.724   15.584    2.100)   26.818
   2.623   (   5.645   14.422   -0.960)   15.518
   2.623   (   5.645   14.422    0.960)   15.518
   3.644   (   9.775    3.311   -1.315)   10.404
   3.644   (   9.775    3.311    1.315)   10.404
   5.584   (   1.845    4.931   -0.402)    5.281
   5.584   (   1.845    4.931    0.402)    5.281
   6.345   (   6.564    4.080   -1.815)    7.939
   6.345   (   6.564    4.080    1.815)    7.939
   6.587   (   0.659   -1.399   -0.696)    1.696
   6.587   (   0.659   -1.399    0.696)    1.696
   7.271   (  -1.839    0.324   -0.283)    1.889
   7.271   (  -1.839    0.324    0.283)    1.889
   7.905   (   1.787    0.346   -0.679)    1.943
   7.905   (   1.787    0.346    0.679)    1.943
   8.131   ( -10.296   -8.141   -0.214)   13.128
   8.131   ( -10.296   -8.141    0.214)   13.128
======================= Grid point 736 (90/111) =======================
q-point: ( 0.33  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.819   (   2.400   11.760   -0.859)   12.033
   2.819   (   2.400   11.760    0.859)   12.033
   2.981   (  18.668   12.219   -2.227)   22.422
   2.981   (  18.668   12.219    2.227)   22.422
   3.756   (   3.179   -1.793   -1.228)    3.851
   3.756   (   3.179   -1.793    1.228)    3.851
   5.654   (   0.651    4.958   -0.498)    5.025
   5.654   (   0.651    4.958    0.498)    5.025
   6.441   (   1.940    0.583   -2.029)    2.867
   6.441   (   1.940    0.583    2.029)    2.867
   6.623   (   3.795    0.199   -0.704)    3.865
   6.623   (   3.795    0.199    0.704)    3.865
   7.243   (  -2.016    0.916   -0.381)    2.247
   7.243   (  -2.016    0.916    0.381)    2.247
   7.908   (  -7.000   -5.620   -0.158)    8.979
   7.908   (  -7.000   -5.620    0.158)    8.979
   7.927   (   0.692   -0.892   -0.524)    1.245
   7.927   (   0.692   -0.892    0.524)    1.245
======================= Grid point 738 (91/111) =======================
q-point: ( 0.39  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.939   (  -0.253    9.048   -0.820)    9.088
   2.939   (  -0.253    9.048    0.820)    9.088
   3.355   (  13.236    7.833   -2.299)   15.551
   3.355   (  13.236    7.833    2.299)   15.551
   3.748   (   0.026   -3.583   -0.963)    3.710
   3.748   (   0.026   -3.583    0.963)    3.710
   5.703   (  -0.654    4.934   -0.556)    5.008
   5.703   (  -0.654    4.934    0.556)    5.008
   6.431   (  -1.695   -1.869   -2.116)    3.293
   6.431   (  -1.695   -1.869    2.116)    3.293
   6.698   (   4.510   -0.275   -0.538)    4.550
   6.698   (   4.510   -0.275    0.538)    4.550
   7.223   (  -1.945    1.823   -0.470)    2.707
   7.223   (  -1.945    1.823    0.470)    2.707
   7.777   (  -3.172   -3.150   -0.140)    4.472
   7.777   (  -3.172   -3.150    0.140)    4.472
   7.928   (   0.478   -1.354   -0.348)    1.477
   7.928   (   0.478   -1.354    0.348)    1.477
======================= Grid point 740 (92/111) =======================
q-point: ( 0.44  0.06 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.993   (  -2.384    6.747   -0.807)    7.201
   2.993   (  -2.384    6.747    0.807)    7.201
   3.565   (   5.077    2.355   -2.325)    6.060
   3.565   (   5.077    2.355    2.325)    6.060
   3.715   (   0.365   -2.771   -0.758)    2.896
   3.715   (   0.365   -2.771    0.758)    2.896
   5.728   (  -1.965    4.736   -0.584)    5.161
   5.728   (  -1.965    4.736    0.584)    5.161
   6.381   (  -1.139   -1.637   -2.073)    2.876
   6.381   (  -1.139   -1.637    2.073)    2.876
   6.753   (   2.984   -2.081   -0.363)    3.656
   6.753   (   2.984   -2.081    0.363)    3.656
   7.213   (  -1.780    2.657   -0.519)    3.240
   7.213   (  -1.780    2.657    0.519)    3.240
   7.727   (  -0.358   -1.496   -0.133)    1.544
   7.727   (  -0.358   -1.496    0.133)    1.544
   7.922   (   0.571   -1.438   -0.142)    1.553
   7.922   (   0.571   -1.438    0.142)    1.553
======================= Grid point 764 (93/111) =======================
q-point: ( 0.11  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.313   (  13.558   23.483   -1.064)   27.137
   1.313   (  13.558   23.483    1.064)   27.137
   1.979   (  12.763   22.105   -1.555)   25.572
   1.979   (  12.763   22.105    1.555)   25.572
   2.741   (  14.392   24.927   -0.285)   28.785
   2.741   (  14.392   24.927    0.285)   28.785
   5.399   (   2.929    5.074   -0.105)    5.860
   5.399   (   2.929    5.074    0.105)    5.860
   5.851   (   6.806   11.788   -1.256)   13.669
   5.851   (   6.806   11.788    1.256)   13.669
   6.780   (  -4.461   -7.727   -0.018)    8.923
   6.780   (  -4.461   -7.727    0.018)    8.923
   7.324   (  -0.144   -0.250   -0.179)    0.339
   7.324   (  -0.144   -0.250    0.179)    0.339
   7.778   (   1.946    3.371   -0.756)    3.965
   7.778   (   1.946    3.371    0.756)    3.965
   8.704   (  -6.433  -11.142   -0.140)   12.866
   8.704   (  -6.433  -11.142    0.140)   12.866
======================= Grid point 766 (94/111) =======================
q-point: ( 0.17  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.806   (  17.411   22.187   -1.742)   28.256
   1.806   (  17.411   22.187    1.742)   28.256
   2.385   (   9.603   21.209   -1.344)   23.321
   2.385   (   9.603   21.209    1.344)   23.321
   3.186   (  14.426   14.559   -0.945)   20.517
   3.186   (  14.426   14.559    0.945)   20.517
   5.505   (   2.957    5.955   -0.196)    6.652
   5.505   (   2.957    5.955    0.196)    6.652
   6.080   (   7.654    9.199   -1.491)   12.059
   6.080   (   7.654    9.199    1.491)   12.059
   6.643   (  -3.464   -5.603   -0.275)    6.593
   6.643   (  -3.464   -5.603    0.275)    6.593
   7.317   (  -0.641    0.090   -0.209)    0.681
   7.317   (  -0.641    0.090    0.209)    0.681
   7.837   (   1.890    2.007   -0.779)    2.865
   7.837   (   1.890    2.007    0.779)    2.865
   8.465   (  -8.468  -11.225   -0.134)   14.062
   8.465   (  -8.468  -11.225    0.134)   14.062
======================= Grid point 768 (95/111) =======================
q-point: ( 0.22  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.320   (  18.705   19.219   -2.025)   26.895
   2.320   (  18.705   19.219    2.025)   26.895
   2.712   (   5.019   19.877   -1.123)   20.532
   2.712   (   5.019   19.877    1.123)   20.532
   3.499   (  10.977    4.828   -1.332)   12.066
   3.499   (  10.977    4.828    1.332)   12.066
   5.611   (   2.146    6.737   -0.305)    7.077
   5.611   (   2.146    6.737    0.305)    7.077
   6.278   (   6.604    5.756   -1.715)    8.926
   6.278   (   6.604    5.756    1.715)    8.926
   6.564   (  -0.674   -2.990   -0.516)    3.108
   6.564   (  -0.674   -2.990    0.516)    3.108
   7.304   (  -1.242    0.715   -0.284)    1.461
   7.304   (  -1.242    0.715    0.284)    1.461
   7.879   (   1.554    0.437   -0.746)    1.778
   7.879   (   1.554    0.437    0.746)    1.778
   8.207   (  -9.130  -10.260   -0.159)   13.735
   8.207   (  -9.130  -10.260    0.159)   13.735
======================= Grid point 770 (96/111) =======================
q-point: ( 0.28  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.803   (  17.698   15.310   -2.136)   23.499
   2.803   (  17.698   15.310    2.136)   23.499
   2.943   (   0.663   17.253   -0.883)   17.288
   2.943   (   0.663   17.253    0.883)   17.288
   3.654   (   6.114   -2.070   -1.288)    6.582
   3.654   (   6.114   -2.070    1.288)    6.582
   5.706   (   0.804    7.368   -0.406)    7.423
   5.706   (   0.804    7.368    0.406)    7.423
   6.404   (   3.560    1.772   -1.916)    4.414
   6.404   (   3.560    1.772    1.916)    4.414
   6.562   (   2.545   -1.148   -0.616)    2.859
   6.562   (   2.545   -1.148    0.616)    2.859
   7.290   (  -1.716    1.715   -0.381)    2.455
   7.290   (  -1.716    1.715    0.381)    2.455
   7.898   (   1.103   -1.109   -0.618)    1.682
   7.898   (   1.103   -1.109    0.618)    1.682
   7.970   (  -7.701   -8.088   -0.141)   11.169
   7.970   (  -7.701   -8.088    0.141)   11.169
======================= Grid point 772 (97/111) =======================
q-point: ( 0.33  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.074   (  -2.402   12.725   -0.619)   12.965
   3.074   (  -2.402   12.725    0.619)   12.965
   3.207   (  14.401   10.723   -2.162)   18.084
   3.207   (  14.401   10.723    2.162)   18.084
   3.690   (   2.331   -3.809   -0.861)    4.547
   3.690   (   2.331   -3.809    0.861)    4.547
   5.782   (  -0.765    7.843   -0.481)    7.894
   5.782   (  -0.765    7.843    0.481)    7.894
   6.430   (  -0.263   -1.682   -2.005)    2.630
   6.430   (  -0.263   -1.682    2.005)    2.630
   6.612   (   4.280   -1.447   -0.537)    4.550
   6.612   (   4.280   -1.447    0.537)    4.550
   7.282   (  -1.955    3.033   -0.474)    3.639
   7.282   (  -1.955    3.033    0.474)    3.639
   7.800   (  -4.335   -5.265   -0.118)    6.821
   7.800   (  -4.335   -5.265    0.118)    6.821
   7.898   (   0.834   -2.190   -0.423)    2.381
   7.898   (   0.834   -2.190    0.423)    2.381
======================= Grid point 774 (98/111) =======================
q-point: ( 0.39  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.121   (  -3.317    7.579   -0.351)    8.280
   3.121   (  -3.317    7.579    0.351)    8.280
   3.485   (   8.418    5.440   -2.134)   10.248
   3.485   (   8.418    5.440    2.134)   10.248
   3.683   (   0.311   -1.362   -0.369)    1.445
   3.683   (   0.311   -1.362    0.369)    1.445
   5.833   (  -2.455    8.043   -0.536)    8.426
   5.833   (  -2.455    8.043    0.536)    8.426
   6.383   (  -1.660   -2.841   -1.922)    3.811
   6.383   (  -1.660   -2.841    1.922)    3.811
   6.663   (   4.140   -3.278   -0.371)    5.294
   6.663   (   4.140   -3.278    0.371)    5.294
   7.283   (  -2.172    4.128   -0.539)    4.696
   7.283   (  -2.172    4.128    0.539)    4.696
   7.716   (  -0.993   -3.116   -0.129)    3.273
   7.716   (  -0.993   -3.116    0.129)    3.273
   7.888   (   0.935   -2.570   -0.207)    2.743
   7.888   (   0.935   -2.570    0.207)    2.743
======================= Grid point 776 (99/111) =======================
q-point: ( 0.44  0.11 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.127   (  -3.100    5.369   -0.190)    6.202
   3.127   (  -3.100    5.369    0.190)    6.202
   3.585   (   0.024   -0.041   -2.068)    2.069
   3.585   (   0.024   -0.041    2.068)    2.069
   3.680   (  -0.212    0.367   -0.197)    0.467
   3.680   (  -0.212    0.367    0.197)    0.467
   5.852   (  -4.331    7.502   -0.563)    8.681
   5.852   (  -4.331    7.502    0.563)    8.681
   6.350   (   0.844   -1.462   -1.832)    2.491
   6.350   (   0.844   -1.462    1.832)    2.491
   6.683   (   2.795   -4.841   -0.288)    5.597
   6.683   (   2.795   -4.841    0.288)    5.597
   7.285   (  -2.531    4.383   -0.563)    5.093
   7.285   (  -2.531    4.383    0.563)    5.093
   7.694   (   1.110   -1.923   -0.133)    2.225
   7.694   (   1.110   -1.923    0.133)    2.225
   7.883   (   1.390   -2.407   -0.050)    2.780
   7.883   (   1.390   -2.407    0.050)    2.780
======================= Grid point 804 (100/111) =======================
q-point: ( 0.17  0.17 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.254   (  13.307   23.048   -1.995)   26.688
   2.254   (  13.307   23.048    1.995)   26.688
   2.782   (  10.747   18.615   -1.293)   21.534
   2.782   (  10.747   18.615    1.293)   21.534
   3.416   (   5.056    8.757   -1.384)   10.205
   3.416   (   5.056    8.757    1.384)   10.205
   5.629   (   3.766    6.523   -0.258)    7.536
   5.629   (   3.766    6.523    0.258)    7.536
   6.247   (   4.394    7.610   -1.678)    8.946
   6.247   (   4.394    7.610    1.678)    8.946
   6.554   (  -1.843   -3.192   -0.420)    3.710
   6.554   (  -1.843   -3.192    0.420)    3.710
   7.319   (   0.120    0.208   -0.305)    0.388
   7.319   (   0.120    0.208    0.305)    0.388
   7.867   (   0.556    0.962   -0.781)    1.358
   7.867   (   0.556    0.962    0.781)    1.358
   8.235   (  -6.928  -12.000   -0.121)   13.857
   8.235   (  -6.928  -12.000    0.121)   13.857
======================= Grid point 806 (101/111) =======================
q-point: ( 0.22  0.17 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.696   (  14.198   18.813   -2.078)   23.661
   2.696   (  14.198   18.813    2.078)   23.661
   3.094   (   3.958   18.211   -1.172)   18.673
   3.094   (   3.958   18.211    1.172)   18.673
   3.533   (   4.868   -1.050   -1.460)    5.190
   3.533   (   4.868   -1.050    1.460)    5.190
   5.754   (   2.689    7.650   -0.337)    8.115
   5.754   (   2.689    7.650    0.337)    8.115
   6.367   (   3.183    3.200   -1.854)    4.879
   6.367   (   3.183    3.200    1.854)    4.879
   6.522   (   0.523   -1.312   -0.515)    1.503
   6.522   (   0.523   -1.312    0.515)    1.503
   7.326   (  -0.328    1.495   -0.411)    1.585
   7.326   (  -0.328    1.495    0.411)    1.585
   7.874   (   0.313   -0.892   -0.683)    1.166
   7.874   (   0.313   -0.892    0.683)    1.166
   8.006   (  -6.896   -9.944   -0.088)   12.101
   8.006   (  -6.896   -9.944    0.088)   12.101
======================= Grid point 808 (102/111) =======================
q-point: ( 0.28  0.17 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.092   (  13.025   13.991   -2.045)   19.224
   3.092   (  13.025   13.991    2.045)   19.224
   3.243   (  -3.301   10.386   -0.927)   10.938
   3.243   (  -3.301   10.386    0.927)   10.938
   3.602   (   5.249   -0.395   -0.650)    5.304
   3.602   (   5.249   -0.395    0.650)    5.304
   5.866   (   0.903    8.648   -0.408)    8.704
   5.866   (   0.903    8.648    0.408)    8.704
   6.411   (   0.685   -1.011   -1.917)    2.273
   6.411   (   0.685   -1.011    1.917)    2.273
   6.537   (   2.461   -1.590   -0.497)    2.971
   6.537   (   2.461   -1.590    0.497)    2.971
   7.339   (  -0.823    3.207   -0.494)    3.347
   7.339   (  -0.823    3.207    0.494)    3.347
   7.818   (  -4.964   -7.124   -0.067)    8.683
   7.818   (  -4.964   -7.124    0.067)    8.683
   7.863   (   0.343   -2.436   -0.488)    2.508
   7.863   (   0.343   -2.436    0.488)    2.508
======================= Grid point 810 (103/111) =======================
q-point: ( 0.33  0.17 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.217   (  -3.447   -0.641   -0.877)    3.614
   3.217   (  -3.447   -0.641    0.877)    3.614
   3.401   (   9.436    8.957   -1.915)   13.151
   3.401   (   9.436    8.957    1.915)   13.151
   3.698   (   1.240    7.003   -0.374)    7.121
   3.698   (   1.240    7.003    0.374)    7.121
   5.952   (  -1.126    9.218   -0.453)    9.297
   5.952   (  -1.126    9.218    0.453)    9.297
   6.379   (  -1.391   -3.303   -1.780)    4.002
   6.379   (  -1.391   -3.303    1.780)    4.002
   6.562   (   3.314   -3.624   -0.370)    4.924
   6.562   (   3.314   -3.624    0.370)    4.924
   7.359   (  -1.455    4.593   -0.550)    4.849
   7.359   (  -1.455    4.593    0.550)    4.849
   7.702   (  -1.963   -4.682   -0.111)    5.079
   7.702   (  -1.963   -4.682    0.111)    5.079
   7.844   (   0.777   -3.162   -0.262)    3.267
   7.844   (   0.777   -3.162    0.262)    3.267
======================= Grid point 812 (104/111) =======================
q-point: ( 0.39  0.17 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.159   (  -0.035   -3.643   -0.626)    3.696
   3.159   (  -0.035   -3.643    0.626)    3.696
   3.574   (   3.210    3.682   -1.786)    5.201
   3.574   (   3.210    3.682    1.786)    5.201
   3.762   (  -3.304    9.173   -0.553)    9.765
   3.762   (  -3.304    9.173    0.553)    9.765
   6.004   (  -3.458    8.982   -0.486)    9.637
   6.004   (  -3.458    8.982    0.486)    9.637
   6.325   (  -0.509   -2.898   -1.511)    3.307
   6.325   (  -0.509   -2.898    1.511)    3.307
   6.576   (   3.481   -5.373   -0.257)    6.407
   6.576   (   3.481   -5.373    0.257)    6.407
   7.375   (  -2.317    4.949   -0.584)    5.495
   7.375   (  -2.317    4.949    0.584)    5.495
   7.653   (   0.572   -3.225   -0.125)    3.278
   7.653   (   0.572   -3.225    0.125)    3.278
   7.831   (   1.399   -3.142   -0.076)    3.440
   7.831   (   1.399   -3.142    0.076)    3.440
======================= Grid point 842 (105/111) =======================
q-point: ( 0.22  0.22 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 262
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.049   (   9.724   16.842   -1.999)   19.550
   3.049   (   9.724   16.842    1.999)   19.550
   3.343   (   0.997    1.727   -1.768)    2.665
   3.343   (   0.997    1.727    1.768)    2.665
   3.537   (   4.135    7.162   -0.918)    8.321
   3.537   (   4.135    7.162    0.918)    8.321
   5.903   (   4.234    7.333   -0.395)    8.477
   5.903   (   4.234    7.333    0.395)    8.477
   6.398   (  -0.054   -0.093   -1.893)    1.896
   6.398   (  -0.054   -0.093    1.893)    1.896
   6.507   (  -0.287   -0.498   -0.475)    0.746
   6.507   (  -0.287   -0.498    0.475)    0.746
   7.363   (   1.330    2.304   -0.510)    2.709
   7.363   (   1.330    2.304    0.510)    2.709
   7.826   (  -4.657   -8.067   -0.018)    9.315
   7.826   (  -4.657   -8.067    0.018)    9.315
   7.847   (  -1.115   -1.931   -0.516)    2.288
   7.847   (  -1.115   -1.931    0.516)    2.288
======================= Grid point 844 (106/111) =======================
q-point: ( 0.28  0.22 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.238   (  -4.314   -7.247   -1.467)    8.560
   3.238   (  -4.314   -7.247    1.467)    8.560
   3.350   (   8.247   11.529   -1.775)   14.285
   3.350   (   8.247   11.529    1.775)   14.285
   3.750   (   5.016   12.569   -0.713)   13.552
   3.750   (   5.016   12.569    0.713)   13.552
   6.034   (   1.997    8.218   -0.450)    8.469
   6.034   (   1.997    8.218    0.450)    8.469
   6.369   (  -1.251   -3.038   -1.725)    3.711
   6.369   (  -1.251   -3.038    1.725)    3.711
   6.494   (   0.609   -2.762   -0.374)    2.853
   6.494   (   0.609   -2.762    0.374)    2.853
   7.412   (   0.527    4.021   -0.561)    4.094
   7.412   (   0.527    4.021    0.561)    4.094
   7.693   (  -2.942   -5.615   -0.091)    6.339
   7.693   (  -2.942   -5.615    0.091)    6.339
   7.809   (  -0.438   -2.936   -0.299)    2.984
   7.809   (  -0.438   -2.936    0.299)    2.984
======================= Grid point 846 (107/111) =======================
q-point: ( 0.33  0.22 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.103   (  -1.469   -8.244   -0.949)    8.427
   3.103   (  -1.469   -8.244    0.949)    8.427
   3.555   (   4.783    6.711   -1.472)    8.372
   3.555   (   4.783    6.711    1.472)    8.372
   3.914   (  -0.492   12.353   -0.685)   12.381
   3.914   (  -0.492   12.353    0.685)   12.381
   6.126   (  -0.654    7.919   -0.412)    7.956
   6.126   (  -0.654    7.919    0.412)    7.956
   6.312   (  -1.365   -3.117   -1.260)    3.629
   6.312   (  -1.365   -3.117    1.260)    3.629
   6.476   (   1.810   -4.887   -0.245)    5.217
   6.476   (   1.810   -4.887    0.245)    5.217
   7.452   (  -0.966    4.657   -0.583)    4.791
   7.452   (  -0.966    4.657    0.583)    4.791
   7.616   (  -0.571   -3.995   -0.108)    4.037
   7.616   (  -0.571   -3.995    0.108)    4.037
   7.780   (   0.703   -3.208   -0.110)    3.286
   7.780   (   0.703   -3.208    0.110)    3.286
======================= Grid point 848 (108/111) =======================
q-point: ( 0.39  0.22 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.047   (   3.538   -6.128   -0.637)    7.105
   3.047   (   3.538   -6.128    0.637)    7.105
   3.631   (  -1.248    2.161   -1.296)    2.812
   3.631   (  -1.248    2.161    1.296)    2.812
   3.971   (  -5.982   10.361   -0.676)   11.982
   3.971   (  -5.982   10.361    0.676)   11.982
   6.162   (  -3.530    6.115   -0.184)    7.063
   6.162   (  -3.530    6.115    0.184)    7.063
   6.281   (   0.615   -1.065   -0.740)    1.435
   6.281   (   0.615   -1.065    0.740)    1.435
   6.466   (   3.079   -5.333   -0.187)    6.161
   6.466   (   3.079   -5.333    0.187)    6.161
   7.467   (  -2.451    4.244   -0.590)    4.936
   7.467   (  -2.451    4.244    0.590)    4.936
   7.592   (   1.642   -2.843   -0.103)    3.285
   7.592   (   1.642   -2.843    0.103)    3.285
   7.770   (   1.677   -2.905   -0.016)    3.354
   7.770   (   1.677   -2.905    0.016)    3.354
======================= Grid point 882 (109/111) =======================
q-point: ( 0.28  0.28 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 261
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.072   (  -4.638   -8.034   -1.042)    9.335
   3.072   (  -4.638   -8.034    1.042)    9.335
   3.547   (   4.727    8.188   -1.328)    9.548
   3.547   (   4.727    8.188    1.328)    9.548
   3.983   (   5.805   10.054   -0.664)   11.628
   3.983   (   5.805   10.054    0.664)   11.628
   6.171   (   2.933    5.081   -0.576)    5.895
   6.171   (   2.933    5.081    0.576)    5.895
   6.309   (  -1.360   -2.355   -1.244)    2.991
   6.309   (  -1.360   -2.355    1.244)    2.991
   6.434   (  -1.783   -3.089   -0.249)    3.575
   6.434   (  -1.783   -3.089    0.249)    3.575
   7.486   (   1.889    3.272   -0.584)    3.823
   7.486   (   1.889    3.272    0.584)    3.823
   7.601   (  -2.145   -3.716   -0.099)    4.292
   7.601   (  -2.145   -3.716    0.099)    4.292
   7.756   (  -1.321   -2.287   -0.121)    2.644
   7.756   (  -1.321   -2.287    0.121)    2.644
======================= Grid point 884 (110/111) =======================
q-point: ( 0.33  0.28 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 486
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.957   (  -0.370   -5.784   -0.519)    5.819
   2.957   (  -0.370   -5.784    0.519)    5.819
   3.659   (   1.305    3.735   -0.806)    4.038
   3.659   (   1.305    3.735    0.806)    4.038
   4.122   (  -0.316    7.985   -0.598)    8.013
   4.122   (  -0.316    7.985    0.598)    8.013
   6.208   (  -0.167   -0.418   -0.629)    0.773
   6.208   (  -0.167   -0.418    0.629)    0.773
   6.306   (  -0.498    3.224   -0.561)    3.310
   6.306   (  -0.498    3.224    0.561)    3.310
   6.386   (   0.502   -3.712   -0.136)    3.749
   6.386   (   0.502   -3.712    0.136)    3.749
   7.531   (  -0.447    3.124   -0.576)    3.208
   7.531   (  -0.447    3.124    0.576)    3.208
   7.550   (  -0.061   -2.531   -0.065)    2.532
   7.550   (  -0.061   -2.531    0.065)    2.532
   7.725   (   0.384   -2.272   -0.026)    2.304
   7.725   (   0.384   -2.272    0.026)    2.304
======================= Grid point 920 (111/111) =======================
q-point: ( 0.33  0.33 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 4.63e-04 4.63e-04 4.63e-04 4.63e-04 6.97e-05 6.97e-05 
Number of triplets: 101
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.897   (  -0.000   -0.000    0.000)    0.000
   2.897   (  -0.000   -0.000    0.000)    0.000
   3.697   (  -0.000   -0.000    0.000)    0.000
   3.697   (  -0.000   -0.000    0.000)    0.000
   4.203   (  -0.000   -0.000   -0.549)    0.549
   4.203   (  -0.000   -0.000    0.549)    0.549
   6.194   (  -0.000   -0.000    0.000)    0.000
   6.194   (  -0.000   -0.000    0.000)    0.000
   6.347   (  -0.000   -0.000    0.000)    0.000
   6.347   (  -0.000   -0.000   -0.000)    0.000
   6.348   (   0.000    0.000   -0.572)    0.572
   6.348   (  -0.000   -0.000    0.572)    0.572
   7.524   (  -0.000   -0.000   -0.000)    0.000
   7.524   (  -0.000   -0.000    0.000)    0.000
   7.565   (  -0.000   -0.000   -0.579)    0.579
   7.565   (  -0.000   -0.000    0.579)    0.579
   7.700   (  -0.000   -0.000    0.000)    0.000
   7.700   (  -0.000   -0.000    0.000)    0.000
=================== 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/23328
   10.0    291.999    291.999     14.899     -0.000     -0.000      0.000 3/23328
   20.0    380.072    380.072      7.031     -0.000     -0.000      0.000 3/23328
   30.0    376.054    376.054      4.880     -0.000     -0.000      0.000 3/23328
   40.0    338.669    338.669      3.802     -0.000     -0.000      0.000 3/23328
   50.0    288.294    288.294      3.092     -0.000     -0.000      0.000 3/23328
   60.0    241.238    241.238      2.579     -0.000     -0.000      0.000 3/23328
   70.0    203.104    203.104      2.199     -0.000     -0.000      0.000 3/23328
   80.0    173.570    173.570      1.911     -0.000     -0.000      0.000 3/23328
   90.0    150.796    150.796      1.688     -0.000     -0.000      0.000 3/23328
  100.0    133.027    133.027      1.512     -0.000     -0.000      0.000 3/23328
  110.0    118.922    118.922      1.369     -0.000     -0.000      0.000 3/23328
  120.0    107.519    107.519      1.252     -0.000     -0.000      0.000 3/23328
  130.0     98.139     98.139      1.153     -0.000     -0.000      0.000 3/23328
  140.0     90.301     90.301      1.070     -0.000     -0.000      0.000 3/23328
  150.0     83.659     83.659      0.998     -0.000     -0.000      0.000 3/23328
  160.0     77.960     77.960      0.935     -0.000     -0.000      0.000 3/23328
  170.0     73.015     73.015      0.880     -0.000     -0.000      0.000 3/23328
  180.0     68.684     68.684      0.832     -0.000     -0.000      0.000 3/23328
  190.0     64.857     64.857      0.788     -0.000     -0.000      0.000 3/23328
  200.0     61.450     61.450      0.749     -0.000     -0.000      0.000 3/23328
  210.0     58.397     58.397      0.714     -0.000     -0.000      0.000 3/23328
  220.0     55.643     55.643      0.682     -0.000     -0.000      0.000 3/23328
  230.0     53.147     53.147      0.653     -0.000     -0.000      0.000 3/23328
  240.0     50.873     50.873      0.626     -0.000     -0.000      0.000 3/23328
  250.0     48.792     48.792      0.602     -0.000     -0.000      0.000 3/23328
  260.0     46.880     46.880      0.579     -0.000     -0.000      0.000 3/23328
  270.0     45.117     45.117      0.558     -0.000     -0.000      0.000 3/23328
  280.0     43.485     43.485      0.539     -0.000     -0.000      0.000 3/23328
  290.0     41.971     41.971      0.521     -0.000     -0.000      0.000 3/23328
  300.0     40.561     40.561      0.504     -0.000     -0.000      0.000 3/23328
  310.0     39.246     39.246      0.488     -0.000     -0.000      0.000 3/23328
  320.0     38.015     38.015      0.473     -0.000     -0.000      0.000 3/23328
  330.0     36.861     36.861      0.459     -0.000     -0.000      0.000 3/23328
  340.0     35.777     35.777      0.446     -0.000     -0.000      0.000 3/23328
  350.0     34.755     34.755      0.434     -0.000     -0.000      0.000 3/23328
  360.0     33.792     33.792      0.422     -0.000     -0.000      0.000 3/23328
  370.0     32.882     32.882      0.411     -0.000     -0.000      0.000 3/23328
  380.0     32.021     32.021      0.400     -0.000     -0.000      0.000 3/23328
  390.0     31.204     31.204      0.390     -0.000     -0.000      0.000 3/23328
  400.0     30.429     30.429      0.381     -0.000     -0.000      0.000 3/23328
  410.0     29.692     29.692      0.372     -0.000     -0.000      0.000 3/23328
  420.0     28.990     28.990      0.363     -0.000     -0.000      0.000 3/23328
  430.0     28.321     28.321      0.355     -0.000     -0.000      0.000 3/23328
  440.0     27.683     27.683      0.347     -0.000     -0.000      0.000 3/23328
  450.0     27.074     27.074      0.340     -0.000     -0.000      0.000 3/23328
  460.0     26.491     26.491      0.333     -0.000     -0.000      0.000 3/23328
  470.0     25.933     25.933      0.326     -0.000     -0.000      0.000 3/23328
  480.0     25.398     25.398      0.319     -0.000     -0.000      0.000 3/23328
  490.0     24.885     24.885      0.313     -0.000     -0.000      0.000 3/23328
  500.0     24.393     24.393      0.307     -0.000     -0.000      0.000 3/23328
  510.0     23.920     23.920      0.301     -0.000     -0.000      0.000 3/23328
  520.0     23.466     23.466      0.295     -0.000     -0.000      0.000 3/23328
  530.0     23.028     23.028      0.290     -0.000     -0.000      0.000 3/23328
  540.0     22.607     22.607      0.285     -0.000     -0.000      0.000 3/23328
  550.0     22.201     22.201      0.280     -0.000     -0.000      0.000 3/23328
  560.0     21.809     21.809      0.275     -0.000     -0.000      0.000 3/23328
  570.0     21.431     21.431      0.270     -0.000     -0.000      0.000 3/23328
  580.0     21.067     21.067      0.266     -0.000     -0.000      0.000 3/23328
  590.0     20.714     20.714      0.261     -0.000     -0.000      0.000 3/23328
  600.0     20.373     20.373      0.257     -0.000     -0.000      0.000 3/23328
  610.0     20.044     20.044      0.253     -0.000     -0.000      0.000 3/23328
  620.0     19.725     19.725      0.249     -0.000     -0.000      0.000 3/23328
  630.0     19.416     19.416      0.245     -0.000     -0.000      0.000 3/23328
  640.0     19.116     19.116      0.241     -0.000     -0.000      0.000 3/23328
  650.0     18.826     18.826      0.238     -0.000     -0.000      0.000 3/23328
  660.0     18.544     18.544      0.234     -0.000     -0.000      0.000 3/23328
  670.0     18.271     18.271      0.231     -0.000     -0.000      0.000 3/23328
  680.0     18.006     18.006      0.228     -0.000     -0.000      0.000 3/23328
  690.0     17.749     17.749      0.224     -0.000     -0.000      0.000 3/23328
  700.0     17.499     17.499      0.221     -0.000     -0.000      0.000 3/23328
  710.0     17.255     17.255      0.218     -0.000     -0.000      0.000 3/23328
  720.0     17.019     17.019      0.215     -0.000     -0.000      0.000 3/23328
  730.0     16.789     16.789      0.212     -0.000     -0.000      0.000 3/23328
  740.0     16.565     16.565      0.210     -0.000     -0.000      0.000 3/23328
  750.0     16.347     16.347      0.207     -0.000     -0.000      0.000 3/23328
  760.0     16.135     16.135      0.204     -0.000     -0.000      0.000 3/23328
  770.0     15.928     15.928      0.202     -0.000     -0.000      0.000 3/23328
  780.0     15.727     15.727      0.199     -0.000     -0.000      0.000 3/23328
  790.0     15.530     15.530      0.197     -0.000     -0.000      0.000 3/23328
  800.0     15.339     15.339      0.194     -0.000     -0.000      0.000 3/23328
  810.0     15.152     15.152      0.192     -0.000     -0.000      0.000 3/23328
  820.0     14.969     14.969      0.190     -0.000     -0.000      0.000 3/23328
  830.0     14.791     14.791      0.187     -0.000     -0.000      0.000 3/23328
  840.0     14.618     14.618      0.185     -0.000     -0.000      0.000 3/23328
  850.0     14.448     14.448      0.183     -0.000     -0.000      0.000 3/23328
  860.0     14.282     14.282      0.181     -0.000     -0.000      0.000 3/23328
  870.0     14.120     14.120      0.179     -0.000     -0.000      0.000 3/23328
  880.0     13.962     13.962      0.177     -0.000     -0.000      0.000 3/23328
  890.0     13.807     13.807      0.175     -0.000     -0.000      0.000 3/23328
  900.0     13.655     13.655      0.173     -0.000     -0.000      0.000 3/23328
  910.0     13.507     13.507      0.171     -0.000     -0.000      0.000 3/23328
  920.0     13.362     13.362      0.170     -0.000     -0.000      0.000 3/23328
  930.0     13.220     13.220      0.168     -0.000     -0.000      0.000 3/23328
  940.0     13.082     13.082      0.166     -0.000     -0.000      0.000 3/23328
  950.0     12.946     12.946      0.164     -0.000     -0.000      0.000 3/23328
  960.0     12.813     12.813      0.163     -0.000     -0.000      0.000 3/23328
  970.0     12.682     12.682      0.161     -0.000     -0.000      0.000 3/23328
  980.0     12.554     12.554      0.159     -0.000     -0.000      0.000 3/23328
  990.0     12.429     12.429      0.158     -0.000     -0.000      0.000 3/23328
 1000.0     12.306     12.306      0.156     -0.000     -0.000      0.000 3/23328

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

