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

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

-------------------------[time 2026-01-08 05:17:49]-------------------------
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: [3 3 2]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P4/nmm (129)
Number of symmetry operations in supercell: 288
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    5.108724409999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.108724409999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.813903529999999
Atomic positions (fractional):
   *1 Li  0.50000000000000  0.00000000000000  0.08382461632719   6.941
    2 Li  0.00000000000000  0.50000000000000  0.91617538367281   6.941
   *3 Br  0.00000000000000  0.50000000000000  0.66377417223600  79.904
    4 Br  0.50000000000000  0.00000000000000  0.33622582776400  79.904
   *5 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904
    6 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904
   *7 Cs  0.00000000000000  0.50000000000000  0.29968304843362 132.905
    8 Cs  0.50000000000000  0.00000000000000  0.70031695156638 132.905
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    5.108724409999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.108724409999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.813903529999999
Atomic positions (fractional):
   *1 Li  0.50000000000000  0.00000000000000  0.08382461632719   6.941 > 1
    2 Li  0.00000000000000  0.50000000000000  0.91617538367281   6.941 > 2
   *3 Br  0.00000000000000  0.50000000000000  0.66377417223600  79.904 > 3
    4 Br  0.50000000000000  0.00000000000000  0.33622582776400  79.904 > 4
   *5 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 5
    6 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 6
   *7 Cs  0.00000000000000  0.50000000000000  0.29968304843362 132.905 > 7
    8 Cs  0.50000000000000  0.00000000000000  0.70031695156638 132.905 > 8
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   15.326173229999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   15.326173229999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   19.627807059999999
Atomic positions (fractional):
   *1 Li  0.16666666666667  0.00000000000000  0.04191230816360   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.04191230816360   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.04191230816360   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.04191230816360   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.04191230816360   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.04191230816360   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.04191230816360   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.04191230816360   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.04191230816360   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.54191230816360   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.54191230816360   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.54191230816360   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.54191230816360   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.54191230816360   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.54191230816360   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.54191230816360   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.54191230816360   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.54191230816360   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.45808769183640   6.941 > 2
   20 Li  0.33333333333333  0.16666666666667  0.45808769183640   6.941 > 2
   21 Li  0.66666666666667  0.16666666666667  0.45808769183640   6.941 > 2
   22 Li  0.00000000000000  0.50000000000000  0.45808769183640   6.941 > 2
   23 Li  0.33333333333333  0.50000000000000  0.45808769183640   6.941 > 2
   24 Li  0.66666666666667  0.50000000000000  0.45808769183640   6.941 > 2
   25 Li  0.00000000000000  0.83333333333333  0.45808769183640   6.941 > 2
   26 Li  0.33333333333333  0.83333333333333  0.45808769183640   6.941 > 2
   27 Li  0.66666666666667  0.83333333333333  0.45808769183640   6.941 > 2
   28 Li  0.00000000000000  0.16666666666667  0.95808769183640   6.941 > 2
   29 Li  0.33333333333333  0.16666666666667  0.95808769183640   6.941 > 2
   30 Li  0.66666666666667  0.16666666666667  0.95808769183640   6.941 > 2
   31 Li  0.00000000000000  0.50000000000000  0.95808769183640   6.941 > 2
   32 Li  0.33333333333333  0.50000000000000  0.95808769183640   6.941 > 2
   33 Li  0.66666666666667  0.50000000000000  0.95808769183640   6.941 > 2
   34 Li  0.00000000000000  0.83333333333333  0.95808769183640   6.941 > 2
   35 Li  0.33333333333333  0.83333333333333  0.95808769183640   6.941 > 2
   36 Li  0.66666666666667  0.83333333333333  0.95808769183640   6.941 > 2
  *37 Br  0.00000000000000  0.16666666666667  0.33188708611800  79.904 > 3
   38 Br  0.33333333333333  0.16666666666667  0.33188708611800  79.904 > 3
   39 Br  0.66666666666667  0.16666666666667  0.33188708611800  79.904 > 3
   40 Br  0.00000000000000  0.50000000000000  0.33188708611800  79.904 > 3
   41 Br  0.33333333333333  0.50000000000000  0.33188708611800  79.904 > 3
   42 Br  0.66666666666667  0.50000000000000  0.33188708611800  79.904 > 3
   43 Br  0.00000000000000  0.83333333333333  0.33188708611800  79.904 > 3
   44 Br  0.33333333333333  0.83333333333333  0.33188708611800  79.904 > 3
   45 Br  0.66666666666667  0.83333333333333  0.33188708611800  79.904 > 3
   46 Br  0.00000000000000  0.16666666666667  0.83188708611800  79.904 > 3
   47 Br  0.33333333333333  0.16666666666667  0.83188708611800  79.904 > 3
   48 Br  0.66666666666667  0.16666666666667  0.83188708611800  79.904 > 3
   49 Br  0.00000000000000  0.50000000000000  0.83188708611800  79.904 > 3
   50 Br  0.33333333333333  0.50000000000000  0.83188708611800  79.904 > 3
   51 Br  0.66666666666667  0.50000000000000  0.83188708611800  79.904 > 3
   52 Br  0.00000000000000  0.83333333333333  0.83188708611800  79.904 > 3
   53 Br  0.33333333333333  0.83333333333333  0.83188708611800  79.904 > 3
   54 Br  0.66666666666667  0.83333333333333  0.83188708611800  79.904 > 3
   55 Br  0.16666666666667  0.00000000000000  0.16811291388200  79.904 > 4
   56 Br  0.50000000000000  0.00000000000000  0.16811291388200  79.904 > 4
   57 Br  0.83333333333333  0.00000000000000  0.16811291388200  79.904 > 4
   58 Br  0.16666666666667  0.33333333333333  0.16811291388200  79.904 > 4
   59 Br  0.50000000000000  0.33333333333333  0.16811291388200  79.904 > 4
   60 Br  0.83333333333333  0.33333333333333  0.16811291388200  79.904 > 4
   61 Br  0.16666666666667  0.66666666666667  0.16811291388200  79.904 > 4
   62 Br  0.50000000000000  0.66666666666667  0.16811291388200  79.904 > 4
   63 Br  0.83333333333333  0.66666666666667  0.16811291388200  79.904 > 4
   64 Br  0.16666666666667  0.00000000000000  0.66811291388200  79.904 > 4
   65 Br  0.50000000000000  0.00000000000000  0.66811291388200  79.904 > 4
   66 Br  0.83333333333333  0.00000000000000  0.66811291388200  79.904 > 4
   67 Br  0.16666666666667  0.33333333333333  0.66811291388200  79.904 > 4
   68 Br  0.50000000000000  0.33333333333333  0.66811291388200  79.904 > 4
   69 Br  0.83333333333333  0.33333333333333  0.66811291388200  79.904 > 4
   70 Br  0.16666666666667  0.66666666666667  0.66811291388200  79.904 > 4
   71 Br  0.50000000000000  0.66666666666667  0.66811291388200  79.904 > 4
   72 Br  0.83333333333333  0.66666666666667  0.66811291388200  79.904 > 4
  *73 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 5
   74 Br  0.33333333333333  0.00000000000000  0.00000000000000  79.904 > 5
   75 Br  0.66666666666667  0.00000000000000  0.00000000000000  79.904 > 5
   76 Br  0.00000000000000  0.33333333333333  0.00000000000000  79.904 > 5
   77 Br  0.33333333333333  0.33333333333333  0.00000000000000  79.904 > 5
   78 Br  0.66666666666667  0.33333333333333  0.00000000000000  79.904 > 5
   79 Br  0.00000000000000  0.66666666666667  0.00000000000000  79.904 > 5
   80 Br  0.33333333333333  0.66666666666667  0.00000000000000  79.904 > 5
   81 Br  0.66666666666667  0.66666666666667  0.00000000000000  79.904 > 5
   82 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 5
   83 Br  0.33333333333333  0.00000000000000  0.50000000000000  79.904 > 5
   84 Br  0.66666666666667  0.00000000000000  0.50000000000000  79.904 > 5
   85 Br  0.00000000000000  0.33333333333333  0.50000000000000  79.904 > 5
   86 Br  0.33333333333333  0.33333333333333  0.50000000000000  79.904 > 5
   87 Br  0.66666666666667  0.33333333333333  0.50000000000000  79.904 > 5
   88 Br  0.00000000000000  0.66666666666667  0.50000000000000  79.904 > 5
   89 Br  0.33333333333333  0.66666666666667  0.50000000000000  79.904 > 5
   90 Br  0.66666666666667  0.66666666666667  0.50000000000000  79.904 > 5
   91 Br  0.16666666666667  0.16666666666667  0.00000000000000  79.904 > 6
   92 Br  0.50000000000000  0.16666666666667  0.00000000000000  79.904 > 6
   93 Br  0.83333333333333  0.16666666666667  0.00000000000000  79.904 > 6
   94 Br  0.16666666666667  0.50000000000000  0.00000000000000  79.904 > 6
   95 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 6
   96 Br  0.83333333333333  0.50000000000000  0.00000000000000  79.904 > 6
   97 Br  0.16666666666667  0.83333333333333  0.00000000000000  79.904 > 6
   98 Br  0.50000000000000  0.83333333333333  0.00000000000000  79.904 > 6
   99 Br  0.83333333333333  0.83333333333333  0.00000000000000  79.904 > 6
  100 Br  0.16666666666667  0.16666666666667  0.50000000000000  79.904 > 6
  101 Br  0.50000000000000  0.16666666666667  0.50000000000000  79.904 > 6
  102 Br  0.83333333333333  0.16666666666667  0.50000000000000  79.904 > 6
  103 Br  0.16666666666667  0.50000000000000  0.50000000000000  79.904 > 6
  104 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 6
  105 Br  0.83333333333333  0.50000000000000  0.50000000000000  79.904 > 6
  106 Br  0.16666666666667  0.83333333333333  0.50000000000000  79.904 > 6
  107 Br  0.50000000000000  0.83333333333333  0.50000000000000  79.904 > 6
  108 Br  0.83333333333333  0.83333333333333  0.50000000000000  79.904 > 6
 *109 Cs  0.00000000000000  0.16666666666667  0.14984152421681 132.905 > 7
  110 Cs  0.33333333333333  0.16666666666667  0.14984152421681 132.905 > 7
  111 Cs  0.66666666666667  0.16666666666667  0.14984152421681 132.905 > 7
  112 Cs  0.00000000000000  0.50000000000000  0.14984152421681 132.905 > 7
  113 Cs  0.33333333333333  0.50000000000000  0.14984152421681 132.905 > 7
  114 Cs  0.66666666666667  0.50000000000000  0.14984152421681 132.905 > 7
  115 Cs  0.00000000000000  0.83333333333333  0.14984152421681 132.905 > 7
  116 Cs  0.33333333333333  0.83333333333333  0.14984152421681 132.905 > 7
  117 Cs  0.66666666666667  0.83333333333333  0.14984152421681 132.905 > 7
  118 Cs  0.00000000000000  0.16666666666667  0.64984152421681 132.905 > 7
  119 Cs  0.33333333333333  0.16666666666667  0.64984152421681 132.905 > 7
  120 Cs  0.66666666666667  0.16666666666667  0.64984152421681 132.905 > 7
  121 Cs  0.00000000000000  0.50000000000000  0.64984152421681 132.905 > 7
  122 Cs  0.33333333333333  0.50000000000000  0.64984152421681 132.905 > 7
  123 Cs  0.66666666666667  0.50000000000000  0.64984152421681 132.905 > 7
  124 Cs  0.00000000000000  0.83333333333333  0.64984152421681 132.905 > 7
  125 Cs  0.33333333333333  0.83333333333333  0.64984152421681 132.905 > 7
  126 Cs  0.66666666666667  0.83333333333333  0.64984152421681 132.905 > 7
  127 Cs  0.16666666666667  0.00000000000000  0.35015847578319 132.905 > 8
  128 Cs  0.50000000000000  0.00000000000000  0.35015847578319 132.905 > 8
  129 Cs  0.83333333333333  0.00000000000000  0.35015847578319 132.905 > 8
  130 Cs  0.16666666666667  0.33333333333333  0.35015847578319 132.905 > 8
  131 Cs  0.50000000000000  0.33333333333333  0.35015847578319 132.905 > 8
  132 Cs  0.83333333333333  0.33333333333333  0.35015847578319 132.905 > 8
  133 Cs  0.16666666666667  0.66666666666667  0.35015847578319 132.905 > 8
  134 Cs  0.50000000000000  0.66666666666667  0.35015847578319 132.905 > 8
  135 Cs  0.83333333333333  0.66666666666667  0.35015847578319 132.905 > 8
  136 Cs  0.16666666666667  0.00000000000000  0.85015847578319 132.905 > 8
  137 Cs  0.50000000000000  0.00000000000000  0.85015847578319 132.905 > 8
  138 Cs  0.83333333333333  0.00000000000000  0.85015847578319 132.905 > 8
  139 Cs  0.16666666666667  0.33333333333333  0.85015847578319 132.905 > 8
  140 Cs  0.50000000000000  0.33333333333333  0.85015847578319 132.905 > 8
  141 Cs  0.83333333333333  0.33333333333333  0.85015847578319 132.905 > 8
  142 Cs  0.16666666666667  0.66666666666667  0.85015847578319 132.905 > 8
  143 Cs  0.50000000000000  0.66666666666667  0.85015847578319 132.905 > 8
  144 Cs  0.83333333333333  0.66666666666667  0.85015847578319 132.905 > 8
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.0411543    0.0000000    0.0000000
            0.0000000    3.0411543    0.0000000
            0.0000000    0.0000000    2.9364170
-------------------------- Born effective charges --------------------------
    1 Li    1.1789278    0.0000000    0.0000000
            0.0000000    1.1789278    0.0000000
            0.0000000    0.0000000    0.9546349
    2 Li    1.1789278    0.0000000    0.0000000
            0.0000000    1.1789278    0.0000000
            0.0000000    0.0000000    0.9546349
    3 Br   -1.1320321    0.0000000    0.0000000
            0.0000000   -1.1320321    0.0000000
            0.0000000    0.0000000   -1.3349396
    4 Br   -1.1320321    0.0000000    0.0000000
            0.0000000   -1.1320321    0.0000000
            0.0000000    0.0000000   -1.3349396
    5 Br   -1.3515940    0.0000000    0.0000000
            0.0000000   -1.3515940    0.0000000
            0.0000000    0.0000000   -0.9134564
    6 Br   -1.3515940    0.0000000    0.0000000
            0.0000000   -1.3515940    0.0000000
            0.0000000    0.0000000   -0.9134564
    7 Cs    1.3046983    0.0000000    0.0000000
            0.0000000    1.3046983    0.0000000
            0.0000000    0.0000000    1.2937610
    8 Cs    1.3046983    0.0000000    0.0000000
            0.0000000    1.3046983    0.0000000
            0.0000000    0.0000000    1.2937610
----------------------------------------------------------------------------
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: 432/432
Permutation basis: 5784/5784
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 436
Number of blocks in projector: 436
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 4
--- Eigsh_solver_block: 1 / 4 ---
Block_size: 160
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 4 ---
Block_size: 128
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 4 ---
Block_size: 108
Use standard eigh solver.
--- Eigsh_solver_block: 4 / 4 ---
Block_size: 40
Use standard eigh solver.
Tree of FC basis block matrices:
- (436, 428), data: False
|-- (40, 40), data: True
|-- (108, 104), data: True
|-- (128, 128), data: True
|-- (160, 156), data: True
-----
Solver_atoms: 1 -- 144 / 144
Time (Solver_compr_matrix_reshape): 0.008
Solver_block: 80 / 80
 - Time: 0.356
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.367
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) 
Permutation basis: 432/432
Permutation basis: 5784/5784
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 436
Number of blocks in projector: 436
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 4
--- Eigsh_solver_block: 1 / 4 ---
Block_size: 160
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 4 ---
Block_size: 128
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 4 ---
Block_size: 108
Use standard eigh solver.
--- Eigsh_solver_block: 4 / 4 ---
Block_size: 40
Use standard eigh solver.
Tree of FC basis block matrices:
- (436, 428), data: False
|-- (40, 40), data: True
|-- (108, 104), data: True
|-- (128, 128), data: True
|-- (160, 156), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 05:17:54]-------------------------
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): [3 3 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P4/nmm (129)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    5.108724409999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.108724409999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.813903529999999
Atomic positions (fractional):
    1 Li  0.50000000000000  0.00000000000000  0.08382461632719   6.941
    2 Li  0.00000000000000  0.50000000000000  0.91617538367281   6.941
    3 Br  0.00000000000000  0.50000000000000  0.66377417223600  79.904
    4 Br  0.50000000000000  0.00000000000000  0.33622582776400  79.904
    5 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904
    6 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904
    7 Cs  0.00000000000000  0.50000000000000  0.29968304843362 132.905
    8 Cs  0.50000000000000  0.00000000000000  0.70031695156638 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.326173229999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   15.326173229999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   19.627807059999999
Atomic positions (fractional):
    1 Li  0.16666666666667  0.00000000000000  0.04191230816360   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.04191230816360   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.04191230816360   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.04191230816360   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.04191230816360   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.04191230816360   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.04191230816360   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.04191230816360   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.04191230816360   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.54191230816360   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.54191230816360   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.54191230816360   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.54191230816360   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.54191230816360   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.54191230816360   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.54191230816360   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.54191230816360   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.54191230816360   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.45808769183640   6.941 > 19
   20 Li  0.33333333333333  0.16666666666667  0.45808769183640   6.941 > 19
   21 Li  0.66666666666667  0.16666666666667  0.45808769183640   6.941 > 19
   22 Li  0.00000000000000  0.50000000000000  0.45808769183640   6.941 > 19
   23 Li  0.33333333333333  0.50000000000000  0.45808769183640   6.941 > 19
   24 Li  0.66666666666667  0.50000000000000  0.45808769183640   6.941 > 19
   25 Li  0.00000000000000  0.83333333333333  0.45808769183640   6.941 > 19
   26 Li  0.33333333333333  0.83333333333333  0.45808769183640   6.941 > 19
   27 Li  0.66666666666667  0.83333333333333  0.45808769183640   6.941 > 19
   28 Li  0.00000000000000  0.16666666666667  0.95808769183640   6.941 > 19
   29 Li  0.33333333333333  0.16666666666667  0.95808769183640   6.941 > 19
   30 Li  0.66666666666667  0.16666666666667  0.95808769183640   6.941 > 19
   31 Li  0.00000000000000  0.50000000000000  0.95808769183640   6.941 > 19
   32 Li  0.33333333333333  0.50000000000000  0.95808769183640   6.941 > 19
   33 Li  0.66666666666667  0.50000000000000  0.95808769183640   6.941 > 19
   34 Li  0.00000000000000  0.83333333333333  0.95808769183640   6.941 > 19
   35 Li  0.33333333333333  0.83333333333333  0.95808769183640   6.941 > 19
   36 Li  0.66666666666667  0.83333333333333  0.95808769183640   6.941 > 19
   37 Br  0.00000000000000  0.16666666666667  0.33188708611800  79.904 > 37
   38 Br  0.33333333333333  0.16666666666667  0.33188708611800  79.904 > 37
   39 Br  0.66666666666667  0.16666666666667  0.33188708611800  79.904 > 37
   40 Br  0.00000000000000  0.50000000000000  0.33188708611800  79.904 > 37
   41 Br  0.33333333333333  0.50000000000000  0.33188708611800  79.904 > 37
   42 Br  0.66666666666667  0.50000000000000  0.33188708611800  79.904 > 37
   43 Br  0.00000000000000  0.83333333333333  0.33188708611800  79.904 > 37
   44 Br  0.33333333333333  0.83333333333333  0.33188708611800  79.904 > 37
   45 Br  0.66666666666667  0.83333333333333  0.33188708611800  79.904 > 37
   46 Br  0.00000000000000  0.16666666666667  0.83188708611800  79.904 > 37
   47 Br  0.33333333333333  0.16666666666667  0.83188708611800  79.904 > 37
   48 Br  0.66666666666667  0.16666666666667  0.83188708611800  79.904 > 37
   49 Br  0.00000000000000  0.50000000000000  0.83188708611800  79.904 > 37
   50 Br  0.33333333333333  0.50000000000000  0.83188708611800  79.904 > 37
   51 Br  0.66666666666667  0.50000000000000  0.83188708611800  79.904 > 37
   52 Br  0.00000000000000  0.83333333333333  0.83188708611800  79.904 > 37
   53 Br  0.33333333333333  0.83333333333333  0.83188708611800  79.904 > 37
   54 Br  0.66666666666667  0.83333333333333  0.83188708611800  79.904 > 37
   55 Br  0.16666666666667  0.00000000000000  0.16811291388200  79.904 > 55
   56 Br  0.50000000000000  0.00000000000000  0.16811291388200  79.904 > 55
   57 Br  0.83333333333333  0.00000000000000  0.16811291388200  79.904 > 55
   58 Br  0.16666666666667  0.33333333333333  0.16811291388200  79.904 > 55
   59 Br  0.50000000000000  0.33333333333333  0.16811291388200  79.904 > 55
   60 Br  0.83333333333333  0.33333333333333  0.16811291388200  79.904 > 55
   61 Br  0.16666666666667  0.66666666666667  0.16811291388200  79.904 > 55
   62 Br  0.50000000000000  0.66666666666667  0.16811291388200  79.904 > 55
   63 Br  0.83333333333333  0.66666666666667  0.16811291388200  79.904 > 55
   64 Br  0.16666666666667  0.00000000000000  0.66811291388200  79.904 > 55
   65 Br  0.50000000000000  0.00000000000000  0.66811291388200  79.904 > 55
   66 Br  0.83333333333333  0.00000000000000  0.66811291388200  79.904 > 55
   67 Br  0.16666666666667  0.33333333333333  0.66811291388200  79.904 > 55
   68 Br  0.50000000000000  0.33333333333333  0.66811291388200  79.904 > 55
   69 Br  0.83333333333333  0.33333333333333  0.66811291388200  79.904 > 55
   70 Br  0.16666666666667  0.66666666666667  0.66811291388200  79.904 > 55
   71 Br  0.50000000000000  0.66666666666667  0.66811291388200  79.904 > 55
   72 Br  0.83333333333333  0.66666666666667  0.66811291388200  79.904 > 55
   73 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 73
   74 Br  0.33333333333333  0.00000000000000  0.00000000000000  79.904 > 73
   75 Br  0.66666666666667  0.00000000000000  0.00000000000000  79.904 > 73
   76 Br  0.00000000000000  0.33333333333333  0.00000000000000  79.904 > 73
   77 Br  0.33333333333333  0.33333333333333  0.00000000000000  79.904 > 73
   78 Br  0.66666666666667  0.33333333333333  0.00000000000000  79.904 > 73
   79 Br  0.00000000000000  0.66666666666667  0.00000000000000  79.904 > 73
   80 Br  0.33333333333333  0.66666666666667  0.00000000000000  79.904 > 73
   81 Br  0.66666666666667  0.66666666666667  0.00000000000000  79.904 > 73
   82 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 73
   83 Br  0.33333333333333  0.00000000000000  0.50000000000000  79.904 > 73
   84 Br  0.66666666666667  0.00000000000000  0.50000000000000  79.904 > 73
   85 Br  0.00000000000000  0.33333333333333  0.50000000000000  79.904 > 73
   86 Br  0.33333333333333  0.33333333333333  0.50000000000000  79.904 > 73
   87 Br  0.66666666666667  0.33333333333333  0.50000000000000  79.904 > 73
   88 Br  0.00000000000000  0.66666666666667  0.50000000000000  79.904 > 73
   89 Br  0.33333333333333  0.66666666666667  0.50000000000000  79.904 > 73
   90 Br  0.66666666666667  0.66666666666667  0.50000000000000  79.904 > 73
   91 Br  0.16666666666667  0.16666666666667  0.00000000000000  79.904 > 91
   92 Br  0.50000000000000  0.16666666666667  0.00000000000000  79.904 > 91
   93 Br  0.83333333333333  0.16666666666667  0.00000000000000  79.904 > 91
   94 Br  0.16666666666667  0.50000000000000  0.00000000000000  79.904 > 91
   95 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 91
   96 Br  0.83333333333333  0.50000000000000  0.00000000000000  79.904 > 91
   97 Br  0.16666666666667  0.83333333333333  0.00000000000000  79.904 > 91
   98 Br  0.50000000000000  0.83333333333333  0.00000000000000  79.904 > 91
   99 Br  0.83333333333333  0.83333333333333  0.00000000000000  79.904 > 91
  100 Br  0.16666666666667  0.16666666666667  0.50000000000000  79.904 > 91
  101 Br  0.50000000000000  0.16666666666667  0.50000000000000  79.904 > 91
  102 Br  0.83333333333333  0.16666666666667  0.50000000000000  79.904 > 91
  103 Br  0.16666666666667  0.50000000000000  0.50000000000000  79.904 > 91
  104 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 91
  105 Br  0.83333333333333  0.50000000000000  0.50000000000000  79.904 > 91
  106 Br  0.16666666666667  0.83333333333333  0.50000000000000  79.904 > 91
  107 Br  0.50000000000000  0.83333333333333  0.50000000000000  79.904 > 91
  108 Br  0.83333333333333  0.83333333333333  0.50000000000000  79.904 > 91
  109 Cs  0.00000000000000  0.16666666666667  0.14984152421681 132.905 > 109
  110 Cs  0.33333333333333  0.16666666666667  0.14984152421681 132.905 > 109
  111 Cs  0.66666666666667  0.16666666666667  0.14984152421681 132.905 > 109
  112 Cs  0.00000000000000  0.50000000000000  0.14984152421681 132.905 > 109
  113 Cs  0.33333333333333  0.50000000000000  0.14984152421681 132.905 > 109
  114 Cs  0.66666666666667  0.50000000000000  0.14984152421681 132.905 > 109
  115 Cs  0.00000000000000  0.83333333333333  0.14984152421681 132.905 > 109
  116 Cs  0.33333333333333  0.83333333333333  0.14984152421681 132.905 > 109
  117 Cs  0.66666666666667  0.83333333333333  0.14984152421681 132.905 > 109
  118 Cs  0.00000000000000  0.16666666666667  0.64984152421681 132.905 > 109
  119 Cs  0.33333333333333  0.16666666666667  0.64984152421681 132.905 > 109
  120 Cs  0.66666666666667  0.16666666666667  0.64984152421681 132.905 > 109
  121 Cs  0.00000000000000  0.50000000000000  0.64984152421681 132.905 > 109
  122 Cs  0.33333333333333  0.50000000000000  0.64984152421681 132.905 > 109
  123 Cs  0.66666666666667  0.50000000000000  0.64984152421681 132.905 > 109
  124 Cs  0.00000000000000  0.83333333333333  0.64984152421681 132.905 > 109
  125 Cs  0.33333333333333  0.83333333333333  0.64984152421681 132.905 > 109
  126 Cs  0.66666666666667  0.83333333333333  0.64984152421681 132.905 > 109
  127 Cs  0.16666666666667  0.00000000000000  0.35015847578319 132.905 > 127
  128 Cs  0.50000000000000  0.00000000000000  0.35015847578319 132.905 > 127
  129 Cs  0.83333333333333  0.00000000000000  0.35015847578319 132.905 > 127
  130 Cs  0.16666666666667  0.33333333333333  0.35015847578319 132.905 > 127
  131 Cs  0.50000000000000  0.33333333333333  0.35015847578319 132.905 > 127
  132 Cs  0.83333333333333  0.33333333333333  0.35015847578319 132.905 > 127
  133 Cs  0.16666666666667  0.66666666666667  0.35015847578319 132.905 > 127
  134 Cs  0.50000000000000  0.66666666666667  0.35015847578319 132.905 > 127
  135 Cs  0.83333333333333  0.66666666666667  0.35015847578319 132.905 > 127
  136 Cs  0.16666666666667  0.00000000000000  0.85015847578319 132.905 > 127
  137 Cs  0.50000000000000  0.00000000000000  0.85015847578319 132.905 > 127
  138 Cs  0.83333333333333  0.00000000000000  0.85015847578319 132.905 > 127
  139 Cs  0.16666666666667  0.33333333333333  0.85015847578319 132.905 > 127
  140 Cs  0.50000000000000  0.33333333333333  0.85015847578319 132.905 > 127
  141 Cs  0.83333333333333  0.33333333333333  0.85015847578319 132.905 > 127
  142 Cs  0.16666666666667  0.66666666666667  0.85015847578319 132.905 > 127
  143 Cs  0.50000000000000  0.66666666666667  0.85015847578319 132.905 > 127
  144 Cs  0.83333333333333  0.66666666666667  0.85015847578319 132.905 > 127
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.0411543    0.0000000    0.0000000
            0.0000000    3.0411543    0.0000000
            0.0000000    0.0000000    2.9364170
-------------------------- Born effective charges --------------------------
    1 Li    1.1789278    0.0000000    0.0000000
            0.0000000    1.1789278    0.0000000
            0.0000000    0.0000000    0.9546349
    2 Li    1.1789278    0.0000000    0.0000000
            0.0000000    1.1789278    0.0000000
            0.0000000    0.0000000    0.9546349
    3 Br   -1.1320321    0.0000000    0.0000000
            0.0000000   -1.1320321    0.0000000
            0.0000000    0.0000000   -1.3349396
    4 Br   -1.1320321    0.0000000    0.0000000
            0.0000000   -1.1320321    0.0000000
            0.0000000    0.0000000   -1.3349396
    5 Br   -1.3515940    0.0000000    0.0000000
            0.0000000   -1.3515940    0.0000000
            0.0000000    0.0000000   -0.9134564
    6 Br   -1.3515940    0.0000000    0.0000000
            0.0000000   -1.3515940    0.0000000
            0.0000000    0.0000000   -0.9134564
    7 Cs    1.3046983    0.0000000    0.0000000
            0.0000000    1.3046983    0.0000000
            0.0000000    0.0000000    1.2937610
    8 Cs    1.3046983    0.0000000    0.0000000
            0.0000000    1.3046983    0.0000000
            0.0000000    0.0000000    1.2937610
----------------------------------------------------------------------------
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[ 37, 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[ 73, 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[ 109, 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 (xzx) 0.00000000 (xzx) 0.00000000 (xxz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 05:18:03]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 05:18:03]-------------------------
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): [3 3 2]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P4/nmm (129)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    5.108724409999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.108724409999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.813903529999999
Atomic positions (fractional):
    1 Li  0.50000000000000  0.00000000000000  0.08382461632719   6.941
    2 Li  0.00000000000000  0.50000000000000  0.91617538367281   6.941
    3 Br  0.00000000000000  0.50000000000000  0.66377417223600  79.904
    4 Br  0.50000000000000  0.00000000000000  0.33622582776400  79.904
    5 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904
    6 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904
    7 Cs  0.00000000000000  0.50000000000000  0.29968304843362 132.905
    8 Cs  0.50000000000000  0.00000000000000  0.70031695156638 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.326173229999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   15.326173229999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   19.627807059999999
Atomic positions (fractional):
    1 Li  0.16666666666667  0.00000000000000  0.04191230816360   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.04191230816360   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.04191230816360   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.04191230816360   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.04191230816360   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.04191230816360   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.04191230816360   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.04191230816360   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.04191230816360   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.54191230816360   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.54191230816360   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.54191230816360   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.54191230816360   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.54191230816360   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.54191230816360   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.54191230816360   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.54191230816360   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.54191230816360   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.45808769183640   6.941 > 19
   20 Li  0.33333333333333  0.16666666666667  0.45808769183640   6.941 > 19
   21 Li  0.66666666666667  0.16666666666667  0.45808769183640   6.941 > 19
   22 Li  0.00000000000000  0.50000000000000  0.45808769183640   6.941 > 19
   23 Li  0.33333333333333  0.50000000000000  0.45808769183640   6.941 > 19
   24 Li  0.66666666666667  0.50000000000000  0.45808769183640   6.941 > 19
   25 Li  0.00000000000000  0.83333333333333  0.45808769183640   6.941 > 19
   26 Li  0.33333333333333  0.83333333333333  0.45808769183640   6.941 > 19
   27 Li  0.66666666666667  0.83333333333333  0.45808769183640   6.941 > 19
   28 Li  0.00000000000000  0.16666666666667  0.95808769183640   6.941 > 19
   29 Li  0.33333333333333  0.16666666666667  0.95808769183640   6.941 > 19
   30 Li  0.66666666666667  0.16666666666667  0.95808769183640   6.941 > 19
   31 Li  0.00000000000000  0.50000000000000  0.95808769183640   6.941 > 19
   32 Li  0.33333333333333  0.50000000000000  0.95808769183640   6.941 > 19
   33 Li  0.66666666666667  0.50000000000000  0.95808769183640   6.941 > 19
   34 Li  0.00000000000000  0.83333333333333  0.95808769183640   6.941 > 19
   35 Li  0.33333333333333  0.83333333333333  0.95808769183640   6.941 > 19
   36 Li  0.66666666666667  0.83333333333333  0.95808769183640   6.941 > 19
   37 Br  0.00000000000000  0.16666666666667  0.33188708611800  79.904 > 37
   38 Br  0.33333333333333  0.16666666666667  0.33188708611800  79.904 > 37
   39 Br  0.66666666666667  0.16666666666667  0.33188708611800  79.904 > 37
   40 Br  0.00000000000000  0.50000000000000  0.33188708611800  79.904 > 37
   41 Br  0.33333333333333  0.50000000000000  0.33188708611800  79.904 > 37
   42 Br  0.66666666666667  0.50000000000000  0.33188708611800  79.904 > 37
   43 Br  0.00000000000000  0.83333333333333  0.33188708611800  79.904 > 37
   44 Br  0.33333333333333  0.83333333333333  0.33188708611800  79.904 > 37
   45 Br  0.66666666666667  0.83333333333333  0.33188708611800  79.904 > 37
   46 Br  0.00000000000000  0.16666666666667  0.83188708611800  79.904 > 37
   47 Br  0.33333333333333  0.16666666666667  0.83188708611800  79.904 > 37
   48 Br  0.66666666666667  0.16666666666667  0.83188708611800  79.904 > 37
   49 Br  0.00000000000000  0.50000000000000  0.83188708611800  79.904 > 37
   50 Br  0.33333333333333  0.50000000000000  0.83188708611800  79.904 > 37
   51 Br  0.66666666666667  0.50000000000000  0.83188708611800  79.904 > 37
   52 Br  0.00000000000000  0.83333333333333  0.83188708611800  79.904 > 37
   53 Br  0.33333333333333  0.83333333333333  0.83188708611800  79.904 > 37
   54 Br  0.66666666666667  0.83333333333333  0.83188708611800  79.904 > 37
   55 Br  0.16666666666667  0.00000000000000  0.16811291388200  79.904 > 55
   56 Br  0.50000000000000  0.00000000000000  0.16811291388200  79.904 > 55
   57 Br  0.83333333333333  0.00000000000000  0.16811291388200  79.904 > 55
   58 Br  0.16666666666667  0.33333333333333  0.16811291388200  79.904 > 55
   59 Br  0.50000000000000  0.33333333333333  0.16811291388200  79.904 > 55
   60 Br  0.83333333333333  0.33333333333333  0.16811291388200  79.904 > 55
   61 Br  0.16666666666667  0.66666666666667  0.16811291388200  79.904 > 55
   62 Br  0.50000000000000  0.66666666666667  0.16811291388200  79.904 > 55
   63 Br  0.83333333333333  0.66666666666667  0.16811291388200  79.904 > 55
   64 Br  0.16666666666667  0.00000000000000  0.66811291388200  79.904 > 55
   65 Br  0.50000000000000  0.00000000000000  0.66811291388200  79.904 > 55
   66 Br  0.83333333333333  0.00000000000000  0.66811291388200  79.904 > 55
   67 Br  0.16666666666667  0.33333333333333  0.66811291388200  79.904 > 55
   68 Br  0.50000000000000  0.33333333333333  0.66811291388200  79.904 > 55
   69 Br  0.83333333333333  0.33333333333333  0.66811291388200  79.904 > 55
   70 Br  0.16666666666667  0.66666666666667  0.66811291388200  79.904 > 55
   71 Br  0.50000000000000  0.66666666666667  0.66811291388200  79.904 > 55
   72 Br  0.83333333333333  0.66666666666667  0.66811291388200  79.904 > 55
   73 Br  0.00000000000000  0.00000000000000  0.00000000000000  79.904 > 73
   74 Br  0.33333333333333  0.00000000000000  0.00000000000000  79.904 > 73
   75 Br  0.66666666666667  0.00000000000000  0.00000000000000  79.904 > 73
   76 Br  0.00000000000000  0.33333333333333  0.00000000000000  79.904 > 73
   77 Br  0.33333333333333  0.33333333333333  0.00000000000000  79.904 > 73
   78 Br  0.66666666666667  0.33333333333333  0.00000000000000  79.904 > 73
   79 Br  0.00000000000000  0.66666666666667  0.00000000000000  79.904 > 73
   80 Br  0.33333333333333  0.66666666666667  0.00000000000000  79.904 > 73
   81 Br  0.66666666666667  0.66666666666667  0.00000000000000  79.904 > 73
   82 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 73
   83 Br  0.33333333333333  0.00000000000000  0.50000000000000  79.904 > 73
   84 Br  0.66666666666667  0.00000000000000  0.50000000000000  79.904 > 73
   85 Br  0.00000000000000  0.33333333333333  0.50000000000000  79.904 > 73
   86 Br  0.33333333333333  0.33333333333333  0.50000000000000  79.904 > 73
   87 Br  0.66666666666667  0.33333333333333  0.50000000000000  79.904 > 73
   88 Br  0.00000000000000  0.66666666666667  0.50000000000000  79.904 > 73
   89 Br  0.33333333333333  0.66666666666667  0.50000000000000  79.904 > 73
   90 Br  0.66666666666667  0.66666666666667  0.50000000000000  79.904 > 73
   91 Br  0.16666666666667  0.16666666666667  0.00000000000000  79.904 > 91
   92 Br  0.50000000000000  0.16666666666667  0.00000000000000  79.904 > 91
   93 Br  0.83333333333333  0.16666666666667  0.00000000000000  79.904 > 91
   94 Br  0.16666666666667  0.50000000000000  0.00000000000000  79.904 > 91
   95 Br  0.50000000000000  0.50000000000000  0.00000000000000  79.904 > 91
   96 Br  0.83333333333333  0.50000000000000  0.00000000000000  79.904 > 91
   97 Br  0.16666666666667  0.83333333333333  0.00000000000000  79.904 > 91
   98 Br  0.50000000000000  0.83333333333333  0.00000000000000  79.904 > 91
   99 Br  0.83333333333333  0.83333333333333  0.00000000000000  79.904 > 91
  100 Br  0.16666666666667  0.16666666666667  0.50000000000000  79.904 > 91
  101 Br  0.50000000000000  0.16666666666667  0.50000000000000  79.904 > 91
  102 Br  0.83333333333333  0.16666666666667  0.50000000000000  79.904 > 91
  103 Br  0.16666666666667  0.50000000000000  0.50000000000000  79.904 > 91
  104 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 91
  105 Br  0.83333333333333  0.50000000000000  0.50000000000000  79.904 > 91
  106 Br  0.16666666666667  0.83333333333333  0.50000000000000  79.904 > 91
  107 Br  0.50000000000000  0.83333333333333  0.50000000000000  79.904 > 91
  108 Br  0.83333333333333  0.83333333333333  0.50000000000000  79.904 > 91
  109 Cs  0.00000000000000  0.16666666666667  0.14984152421681 132.905 > 109
  110 Cs  0.33333333333333  0.16666666666667  0.14984152421681 132.905 > 109
  111 Cs  0.66666666666667  0.16666666666667  0.14984152421681 132.905 > 109
  112 Cs  0.00000000000000  0.50000000000000  0.14984152421681 132.905 > 109
  113 Cs  0.33333333333333  0.50000000000000  0.14984152421681 132.905 > 109
  114 Cs  0.66666666666667  0.50000000000000  0.14984152421681 132.905 > 109
  115 Cs  0.00000000000000  0.83333333333333  0.14984152421681 132.905 > 109
  116 Cs  0.33333333333333  0.83333333333333  0.14984152421681 132.905 > 109
  117 Cs  0.66666666666667  0.83333333333333  0.14984152421681 132.905 > 109
  118 Cs  0.00000000000000  0.16666666666667  0.64984152421681 132.905 > 109
  119 Cs  0.33333333333333  0.16666666666667  0.64984152421681 132.905 > 109
  120 Cs  0.66666666666667  0.16666666666667  0.64984152421681 132.905 > 109
  121 Cs  0.00000000000000  0.50000000000000  0.64984152421681 132.905 > 109
  122 Cs  0.33333333333333  0.50000000000000  0.64984152421681 132.905 > 109
  123 Cs  0.66666666666667  0.50000000000000  0.64984152421681 132.905 > 109
  124 Cs  0.00000000000000  0.83333333333333  0.64984152421681 132.905 > 109
  125 Cs  0.33333333333333  0.83333333333333  0.64984152421681 132.905 > 109
  126 Cs  0.66666666666667  0.83333333333333  0.64984152421681 132.905 > 109
  127 Cs  0.16666666666667  0.00000000000000  0.35015847578319 132.905 > 127
  128 Cs  0.50000000000000  0.00000000000000  0.35015847578319 132.905 > 127
  129 Cs  0.83333333333333  0.00000000000000  0.35015847578319 132.905 > 127
  130 Cs  0.16666666666667  0.33333333333333  0.35015847578319 132.905 > 127
  131 Cs  0.50000000000000  0.33333333333333  0.35015847578319 132.905 > 127
  132 Cs  0.83333333333333  0.33333333333333  0.35015847578319 132.905 > 127
  133 Cs  0.16666666666667  0.66666666666667  0.35015847578319 132.905 > 127
  134 Cs  0.50000000000000  0.66666666666667  0.35015847578319 132.905 > 127
  135 Cs  0.83333333333333  0.66666666666667  0.35015847578319 132.905 > 127
  136 Cs  0.16666666666667  0.00000000000000  0.85015847578319 132.905 > 127
  137 Cs  0.50000000000000  0.00000000000000  0.85015847578319 132.905 > 127
  138 Cs  0.83333333333333  0.00000000000000  0.85015847578319 132.905 > 127
  139 Cs  0.16666666666667  0.33333333333333  0.85015847578319 132.905 > 127
  140 Cs  0.50000000000000  0.33333333333333  0.85015847578319 132.905 > 127
  141 Cs  0.83333333333333  0.33333333333333  0.85015847578319 132.905 > 127
  142 Cs  0.16666666666667  0.66666666666667  0.85015847578319 132.905 > 127
  143 Cs  0.50000000000000  0.66666666666667  0.85015847578319 132.905 > 127
  144 Cs  0.83333333333333  0.66666666666667  0.85015847578319 132.905 > 127
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.0411543    0.0000000    0.0000000
            0.0000000    3.0411543    0.0000000
            0.0000000    0.0000000    2.9364170
-------------------------- Born effective charges --------------------------
    1 Li    1.1789278    0.0000000    0.0000000
            0.0000000    1.1789278    0.0000000
            0.0000000    0.0000000    0.9546349
    2 Li    1.1789278    0.0000000    0.0000000
            0.0000000    1.1789278    0.0000000
            0.0000000    0.0000000    0.9546349
    3 Br   -1.1320321    0.0000000    0.0000000
            0.0000000   -1.1320321    0.0000000
            0.0000000    0.0000000   -1.3349396
    4 Br   -1.1320321    0.0000000    0.0000000
            0.0000000   -1.1320321    0.0000000
            0.0000000    0.0000000   -1.3349396
    5 Br   -1.3515940    0.0000000    0.0000000
            0.0000000   -1.3515940    0.0000000
            0.0000000    0.0000000   -0.9134564
    6 Br   -1.3515940    0.0000000    0.0000000
            0.0000000   -1.3515940    0.0000000
            0.0000000    0.0000000   -0.9134564
    7 Cs    1.3046983    0.0000000    0.0000000
            0.0000000    1.3046983    0.0000000
            0.0000000    0.0000000    1.2937610
    8 Cs    1.3046983    0.0000000    0.0000000
            0.0000000    1.3046983    0.0000000
            0.0000000    0.0000000    1.2937610
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (xzx) 0.00000000 (xzx) 0.00000000 (xxz)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 10 10 5 ]
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.65, Number of G-points: 313, Lambda: 0.12
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/63) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 63
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.807   (   0.000    0.000    0.000)    0.000
   0.807   (   0.000    0.000   -0.000)    0.000
   1.333   (  -0.000    0.000    0.000)    0.000
   1.333   (  -0.000    0.000   -0.000)    0.000
   1.377   (   0.000   -0.000    0.000)    0.000
   1.685   (   0.000    0.000   -0.000)    0.000
   1.763   (   0.000    0.000    0.000)    0.000
   2.160   (  -0.000    0.000   -0.000)    0.000
   2.160   (   0.000    0.000    0.000)    0.000
   2.161   (   0.000    0.000   -0.000)    0.000
   2.161   (   0.000   -0.000   -0.000)    0.000
   2.514   (   0.000    0.000    0.000)    0.000
   3.022   (   0.000    0.000    0.000)    0.000
   3.272   (   0.000    0.000    0.000)    0.000
   3.272   (  -0.000    0.000   -0.000)    0.000
   6.022   (   0.000    0.000   -0.000)    0.000
   6.022   (  -0.000    0.000    0.000)    0.000
   6.052   (   0.000    0.000   -0.000)    0.000
   6.052   (   0.000    0.000    0.000)    0.000
  10.097   (  -0.000    0.000   -0.000)    0.000
  10.586   (  -0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/63) =======================
q-point: ( 0.10  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.215   (  11.462    0.000    0.000)   11.462
   0.426   (  20.635    0.000    0.000)   20.635
   0.492   (  23.554    0.000    0.000)   23.554
   0.868   (   5.856    0.000    0.000)    5.856
   0.924   (  10.848    0.000    0.000)   10.848
   1.404   (   2.670    0.000    0.000)    2.670
   1.427   (   9.338    0.000    0.000)    9.338
   1.523   (  10.665    0.000    0.000)   10.665
   1.653   (  -2.986    0.000    0.000)    2.986
   1.757   (  -0.553    0.000    0.000)    0.553
   2.144   (  -1.711    0.000    0.000)    1.711
   2.147   (  -1.232    0.000    0.000)    1.232
   2.159   (  -0.172    0.000    0.000)    0.172
   2.522   (   0.833    0.000    0.000)    0.833
   2.532   (   0.884    0.000    0.000)    0.884
   3.006   (  -1.616    0.000    0.000)    1.616
   3.252   (  -2.071    0.000    0.000)    2.071
   3.323   (   4.582    0.000    0.000)    4.582
   6.026   (   0.425    0.000    0.000)    0.425
   6.053   (   0.053    0.000    0.000)    0.053
   6.138   (   8.792    0.000    0.000)    8.792
   7.951   (   1.537    0.000    0.000)    1.537
  10.081   (  -1.548    0.000    0.000)    1.548
  10.551   (  -3.552    0.000    0.000)    3.552
======================= Grid point 2 (3/63) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.446   (  12.000    0.000    0.000)   12.000
   0.794   (  16.924    0.000    0.000)   16.924
   0.898   (  17.566    0.000    0.000)   17.566
   1.017   (   8.866    0.000    0.000)    8.866
   1.183   (  14.592    0.000    0.000)   14.592
   1.474   (   4.279    0.000    0.000)    4.279
   1.577   (  -3.995    0.000    0.000)    3.995
   1.671   (  14.509    0.000    0.000)   14.509
   1.742   (  -0.906    0.000    0.000)    0.906
   1.814   (  17.713    0.000    0.000)   17.713
   2.095   (  -3.264    0.000    0.000)    3.264
   2.117   (  -1.476    0.000    0.000)    1.476
   2.152   (  -0.586    0.000    0.000)    0.586
   2.535   (  -0.144    0.000    0.000)    0.144
   2.560   (   2.640    0.000    0.000)    2.640
   2.961   (  -2.880    0.000    0.000)    2.880
   3.190   (  -4.351    0.000    0.000)    4.351
   3.416   (   3.880    0.000    0.000)    3.880
   6.037   (   0.656    0.000    0.000)    0.656
   6.054   (   0.127    0.000    0.000)    0.127
   6.394   (  17.143    0.000    0.000)   17.143
   7.981   (   1.025    0.000    0.000)    1.025
  10.041   (  -2.401    0.000    0.000)    2.401
  10.451   (  -6.627    0.000    0.000)    6.627
======================= Grid point 3 (4/63) =======================
q-point: ( 0.30  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.678   (  11.549    0.000    0.000)   11.549
   1.088   (  13.091    0.000    0.000)   13.091
   1.159   (   8.121    0.000    0.000)    8.121
   1.196   (   8.917    0.000    0.000)    8.917
   1.458   (  12.263    0.000    0.000)   12.263
   1.514   (  -1.468    0.000    0.000)    1.468
   1.569   (   5.419    0.000    0.000)    5.419
   1.730   (   0.144    0.000    0.000)    0.144
   1.908   (   5.626    0.000    0.000)    5.626
   2.019   (  -4.375    0.000    0.000)    4.375
   2.136   (  -1.067    0.000    0.000)    1.067
   2.149   (   8.179    0.000    0.000)    8.179
   2.161   (  15.467    0.000    0.000)   15.467
   2.518   (  -0.917    0.000    0.000)    0.917
   2.649   (   6.966    0.000    0.000)    6.966
   2.893   (  -4.089    0.000    0.000)    4.089
   3.079   (  -6.993    0.000    0.000)    6.993
   3.444   (  -1.449    0.000    0.000)    1.449
   6.050   (   0.604    0.000    0.000)    0.604
   6.058   (   0.191    0.000    0.000)    0.191
   6.791   (  22.684    0.000    0.000)   22.684
   7.974   (  -2.192    0.000    0.000)    2.192
   9.994   (  -2.075    0.000    0.000)    2.075
  10.299   (  -8.529    0.000    0.000)    8.529
======================= Grid point 4 (5/63) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.893   (  10.242    0.000    0.000)   10.242
   1.193   (  -3.543    0.000    0.000)    3.543
   1.303   (   8.574    0.000    0.000)    8.574
   1.348   (   6.040    0.000    0.000)    6.040
   1.534   (   2.566    0.000    0.000)    2.566
   1.596   (   1.416    0.000    0.000)    1.416
   1.692   (   6.812    0.000    0.000)    6.812
   1.784   (   4.652    0.000    0.000)    4.652
   1.908   (  -2.330    0.000    0.000)    2.330
   1.932   (  -4.102    0.000    0.000)    4.102
   2.124   (   0.820    0.000    0.000)    0.820
   2.299   (  -1.924    0.000    0.000)    1.924
   2.409   (  15.058    0.000    0.000)   15.058
   2.565   (   6.097    0.000    0.000)    6.097
   2.802   (  -5.205    0.000    0.000)    5.205
   2.856   (  13.857    0.000    0.000)   13.857
   2.914   (  -9.990    0.000    0.000)    9.990
   3.350   (  -8.100    0.000    0.000)    8.100
   6.059   (   0.338    0.000    0.000)    0.338
   6.061   (   0.160    0.000    0.000)    0.160
   7.244   (  22.523    0.000    0.000)   22.523
   7.876   (  -8.333    0.000    0.000)    8.333
   9.971   (   0.102    0.000    0.000)    0.102
  10.133   (  -7.876    0.000    0.000)    7.876
======================= Grid point 5 (6/63) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 54
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.073   (  -8.020    0.000    0.000)    8.020
   1.073   (   8.020    0.000    0.000)    8.020
   1.403   (  -1.099    0.000    0.000)    1.099
   1.403   (   1.099    0.000    0.000)    1.099
   1.581   (  -1.939    0.000    0.000)    1.939
   1.581   (   1.939    0.000    0.000)    1.939
   1.803   (  -3.367    0.000    0.000)    3.367
   1.803   (   3.367    0.000    0.000)    3.367
   1.881   (  -0.504    0.000    0.000)    0.504
   1.881   (   0.504    0.000    0.000)    0.504
   2.199   (  -6.104    0.000    0.000)    6.104
   2.199   (   6.104    0.000    0.000)    6.104
   2.688   ( -13.042    0.000    0.000)   13.042
   2.688   (  13.042    0.000    0.000)   13.042
   2.691   (  -6.130    0.000    0.000)    6.130
   2.691   (   6.130    0.000    0.000)    6.130
   3.136   ( -13.413    0.000    0.000)   13.413
   3.136   (  13.413    0.000    0.000)   13.413
   6.063   (  -0.031    0.000    0.000)    0.031
   6.063   (   0.031    0.000    0.000)    0.031
   7.634   ( -16.484    0.000    0.000)   16.484
   7.634   (  16.484    0.000    0.000)   16.484
  10.012   (  -4.206    0.000    0.000)    4.206
  10.012   (   4.206    0.000    0.000)    4.206
======================= Grid point 12 (7/63) =======================
q-point: ( 0.10  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.097   (   4.653    4.653    0.000)    6.580
   0.308   (   8.208    8.208    0.000)   11.608
   0.864   (   2.766    2.766    0.000)    3.912
   0.866   (  18.476   18.476    0.000)   26.129
   1.065   (  10.614   10.614    0.000)   15.011
   1.347   (   0.583    0.583    0.000)    0.824
   1.428   (   2.294    2.294    0.000)    3.244
   1.622   (  -2.544   -2.544    0.000)    3.597
   1.752   (  -0.524   -0.524    0.000)    0.741
   1.774   (  15.956   15.956    0.000)   22.566
   2.116   (  -2.219   -2.219    0.000)    3.138
   2.169   (   0.360    0.360    0.000)    0.509
   2.171   (   0.535    0.535    0.000)    0.757
   2.490   (  -1.963   -1.963    0.000)    2.777
   2.536   (   1.454    1.454    0.000)    2.056
   2.988   (  -1.694   -1.694    0.000)    2.396
   3.271   (  -1.152   -1.152    0.000)    1.629
   3.324   (   2.661    2.661    0.000)    3.763
   5.947   (  -3.079   -3.079    0.000)    4.354
   6.063   (   0.245    0.245    0.000)    0.346
   6.219   (   9.150    9.150    0.000)   12.940
   8.046   (   5.513    5.513    0.000)    7.796
  10.066   (  -1.461   -1.461    0.000)    2.067
  10.521   (  -3.127   -3.127    0.000)    4.423
======================= Grid point 13 (8/63) =======================
q-point: ( 0.20  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.481   (  21.144  -17.121    0.000)   27.206
   0.497   (  10.533    4.919    0.000)   11.625
   0.979   (   8.231   -1.400    0.000)    8.349
   1.172   (  12.466   14.204    0.000)   18.898
   1.293   (  12.007    8.492    0.000)   14.706
   1.450   (   8.149   -6.759    0.000)   10.587
   1.502   (   6.368   -1.704    0.000)    6.592
   1.562   (  -3.054   -1.301    0.000)    3.320
   1.737   (  -1.036   -0.619    0.000)    1.207
   2.063   (  -3.013   -2.477    0.000)    3.901
   2.081   (  10.551   14.441    0.000)   17.885
   2.163   (   2.667    4.850    0.000)    5.535
   2.169   (  -0.442    1.043    0.000)    1.132
   2.470   (  -0.525   -6.030    0.000)    6.053
   2.593   (   4.972    4.128    0.000)    6.462
   2.942   (  -2.966   -1.872    0.000)    3.507
   3.190   (  -6.304   -1.636    0.000)    6.513
   3.412   (   4.493    0.855    0.000)    4.573
   5.935   (   1.112   -6.851    0.000)    6.941
   6.059   (  -0.542    0.027    0.000)    0.543
   6.496   (  18.675   10.741    0.000)   21.543
   8.126   (   2.642   11.579    0.000)   11.877
  10.028   (  -2.265   -1.244    0.000)    2.584
  10.431   (  -6.001   -2.054    0.000)    6.343
======================= Grid point 14 (9/63) =======================
q-point: ( 0.30  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.707   (  10.687    2.886    0.000)   11.069
   0.852   (  16.397  -13.986    0.000)   21.551
   1.159   (   9.605   -2.746    0.000)    9.989
   1.322   (   2.083    6.404    0.000)    6.734
   1.507   (   8.995    4.024    0.000)    9.854
   1.524   (   0.656    0.904    0.000)    1.117
   1.587   (   6.353    2.242    0.000)    6.737
   1.713   (  -1.299   -1.878    0.000)    2.283
   1.734   (  14.595  -13.853    0.000)   20.123
   2.007   (  -2.231   -0.562    0.000)    2.301
   2.071   (  -2.308    1.608    0.000)    2.813
   2.151   (  -1.262    1.293    0.000)    1.807
   2.343   (   5.114    4.149    0.000)    6.585
   2.472   (   2.100   -3.660    0.000)    4.220
   2.760   (  12.252   11.581    0.000)   16.859
   2.873   (  -4.096   -1.987    0.000)    4.553
   3.030   (  -9.987   -5.927    0.000)   11.613
   3.466   (   0.584    2.974    0.000)    3.030
   5.969   (   1.930   -6.027    0.000)    6.328
   6.047   (  -0.600   -0.995    0.000)    1.161
   6.925   (  24.294   13.341    0.000)   27.716
   8.142   (  -1.354   14.729    0.000)   14.791
   9.985   (  -1.916   -0.963    0.000)    2.144
  10.291   (  -8.062   -0.923    0.000)    8.115
======================= Grid point 15 (10/63) =======================
q-point: ( 0.40  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.907   (   9.523    1.422    0.000)    9.629
   1.096   (   6.612   -8.227    0.000)   10.555
   1.312   (   0.134    0.234    0.000)    0.270
   1.328   (   6.794   -1.905    0.000)    7.056
   1.556   (   1.863    1.632    0.000)    2.477
   1.605   (   0.649    0.858    0.000)    1.076
   1.694   (  -0.442   -7.203    0.000)    7.217
   1.697   (   3.332    0.209    0.000)    3.338
   1.920   (   2.437   -0.218    0.000)    2.447
   1.995   (   0.488    4.445    0.000)    4.472
   2.131   (   0.128    0.670    0.000)    0.682
   2.152   (  11.666  -11.913    0.000)   16.674
   2.335   (  -1.937   -4.145    0.000)    4.575
   2.566   (   5.810   -0.898    0.000)    5.879
   2.781   (  -6.168   -2.846    0.000)    6.793
   2.802   ( -12.374  -10.016    0.000)   15.920
   3.030   (  13.954   13.897    0.000)   19.694
   3.424   (  -4.977    6.801    0.000)    8.428
   6.003   (   1.440   -4.237    0.000)    4.475
   6.036   (  -0.530   -1.985    0.000)    2.054
   7.406   (  23.684   15.550    0.000)   28.332
   8.055   (  -8.061   16.353    0.000)   18.231
   9.965   (   0.263   -0.665    0.000)    0.715
  10.130   (  -7.776   -0.355    0.000)    7.784
======================= Grid point 16 (11/63) =======================
q-point: (-0.50  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.069   (  -6.531   -0.487    0.000)    6.549
   1.069   (   6.531   -0.487    0.000)    6.549
   1.376   (  -2.696   -2.751    0.000)    3.852
   1.376   (   2.696   -2.751    0.000)    3.852
   1.585   (  -1.278    0.240    0.000)    1.300
   1.585   (   1.278    0.240    0.000)    1.300
   1.701   (  -1.272   -6.405    0.000)    6.530
   1.701   (   1.272   -6.405    0.000)    6.530
   1.970   (  -2.945    5.245    0.000)    6.016
   1.970   (   2.945    5.245    0.000)    6.016
   2.190   (  -4.756   -0.792    0.000)    4.822
   2.190   (   4.756   -0.792    0.000)    4.822
   2.509   ( -15.499  -14.767    0.000)   21.408
   2.509   (  15.499  -14.767    0.000)   21.408
   2.679   (  -5.655   -1.633    0.000)    5.886
   2.679   (   5.655   -1.633    0.000)    5.886
   3.273   ( -10.345   11.016    0.000)   15.112
   3.273   (  10.345   11.016    0.000)   15.112
   6.024   (  -0.783   -2.914    0.000)    3.018
   6.024   (   0.783   -2.914    0.000)    3.018
   7.811   ( -16.934   16.663    0.000)   23.757
   7.811   (  16.934   16.663    0.000)   23.757
  10.008   (  -4.323   -0.399    0.000)    4.341
  10.008   (   4.323   -0.399    0.000)    4.341
======================= Grid point 24 (12/63) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 93
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.328   (   6.597    6.597    0.000)    9.330
   0.624   (   7.702    7.702    0.000)   10.893
   1.002   (   3.881    3.881    0.000)    5.488
   1.360   (  -0.172   -0.172    0.000)    0.244
   1.381   (   5.104    5.104    0.000)    7.218
   1.461   (   7.865    7.865    0.000)   11.123
   1.526   (   2.384    2.384    0.000)    3.372
   1.558   (   2.679    2.679    0.000)    3.789
   1.713   (  -1.995   -1.995    0.000)    2.821
   2.019   (  -1.610   -1.610    0.000)    2.277
   2.182   (   0.218    0.218    0.000)    0.308
   2.197   (   0.725    0.725    0.000)    1.026
   2.323   (   0.590    0.590    0.000)    0.834
   2.354   (   1.146    1.146    0.000)    1.620
   2.758   (  12.807   12.807    0.000)   18.112
   2.892   (  -3.196   -3.196    0.000)    4.520
   3.076   (  -9.702   -9.702    0.000)   13.720
   3.477   (   4.991    4.991    0.000)    7.058
   5.861   (  -0.488   -0.488    0.000)    0.690
   6.048   (  -1.047   -1.047    0.000)    1.480
   6.813   (  21.081   21.081    0.000)   29.813
   8.335   (   8.699    8.699    0.000)   12.302
   9.996   (  -1.930   -1.930    0.000)    2.729
  10.369   (  -4.326   -4.326    0.000)    6.118
======================= Grid point 25 (13/63) =======================
q-point: ( 0.30  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.612   (  17.368   -8.846    0.000)   19.491
   0.786   (   8.530    4.983    0.000)    9.879
   1.120   (   7.705   -0.808    0.000)    7.747
   1.374   (  -2.594   -0.218    0.000)    2.603
   1.451   (   7.941   -9.587    0.000)   12.449
   1.577   (   3.522    2.200    0.000)    4.153
   1.578   (   3.246    0.335    0.000)    3.264
   1.651   (  -4.234   -4.054    0.000)    5.862
   1.666   (   5.497    4.759    0.000)    7.271
   2.019   (   1.917    1.988    0.000)    2.762
   2.125   (  -5.907    2.510    0.000)    6.418
   2.175   (  -0.972    0.903    0.000)    1.327
   2.286   (  -0.674   -6.580    0.000)    6.615
   2.449   (   4.502    0.591    0.000)    4.540
   2.819   (  -4.188   -3.392    0.000)    5.389
   2.833   ( -14.660  -13.780    0.000)   20.120
   3.044   (  14.980   15.568    0.000)   21.605
   3.565   (   3.219    6.428    0.000)    7.189
   5.893   (   2.976   -0.908    0.000)    3.111
   6.025   (  -1.077   -1.012    0.000)    1.478
   7.288   (  26.454   22.924    0.000)   35.005
   8.449   (   2.654   14.822    0.000)   15.058
   9.959   (  -1.557   -1.539    0.000)    2.189
  10.260   (  -6.718   -2.359    0.000)    7.121
======================= Grid point 26 (14/63) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 156
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.919   (  12.849   -8.951    0.000)   15.659
   0.946   (   7.464    2.465    0.000)    7.860
   1.279   (   7.332   -2.962    0.000)    7.908
   1.309   (  -3.176   -0.418    0.000)    3.203
   1.534   (  -5.234   -8.603    0.000)   10.070
   1.573   (   2.397   -1.053    0.000)    2.618
   1.626   (   0.396    0.964    0.000)    1.042
   1.670   (  -0.738   -4.422    0.000)    4.483
   1.771   (  12.483  -11.315    0.000)   16.848
   2.059   (   0.409    1.438    0.000)    1.495
   2.078   (   2.834    3.998    0.000)    4.901
   2.148   (  -1.185    0.786    0.000)    1.422
   2.260   (  -2.194   -2.987    0.000)    3.706
   2.514   ( -17.457  -17.420    0.000)   24.661
   2.539   (   4.775   -1.569    0.000)    5.026
   2.732   (  -4.759   -3.213    0.000)    5.742
   3.312   (  11.957   13.957    0.000)   18.378
   3.584   (  -1.544    8.785    0.000)    8.920
   5.952   (   2.621   -0.271    0.000)    2.635
   6.008   (  -0.772   -0.533    0.000)    0.938
   7.802   (  24.736   23.587    0.000)   34.179
   8.424   (  -5.626   19.485    0.000)   20.281
   9.946   (   0.579   -1.222    0.000)    1.352
  10.117   (  -7.408   -1.173    0.000)    7.500
======================= Grid point 27 (15/63) =======================
q-point: (-0.50  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.048   (  -1.103   -1.998    0.000)    2.282
   1.048   (   1.103   -1.998    0.000)    2.282
   1.306   (  -3.419   -3.839    0.000)    5.140
   1.306   (   3.419   -3.839    0.000)    5.140
   1.541   (  -5.741   -5.975    0.000)    8.287
   1.541   (   5.741   -5.975    0.000)    8.287
   1.624   (  -0.083   -0.341    0.000)    0.351
   1.624   (   0.083   -0.341    0.000)    0.351
   2.036   (  -8.112   -2.456    0.000)    8.476
   2.036   (   8.112   -2.456    0.000)    8.476
   2.148   (  -6.676   -6.543    0.000)    9.348
   2.148   (   6.676   -6.543    0.000)    9.348
   2.221   (  -3.155   -6.611    0.000)    7.325
   2.221   (   3.155   -6.611    0.000)    7.325
   2.636   (  -5.016   -2.642    0.000)    5.669
   2.636   (   5.016   -2.642    0.000)    5.669
   3.500   (  -7.077   11.414    0.000)   13.430
   3.500   (   7.077   11.414    0.000)   13.430
   5.989   (  -1.279   -0.075    0.000)    1.281
   5.989   (   1.279   -0.075    0.000)    1.281
   8.212   ( -16.272   22.537    0.000)   27.798
   8.212   (  16.272   22.537    0.000)   27.798
   9.995   (  -4.533   -1.019    0.000)    4.646
   9.995   (   4.533   -1.019    0.000)    4.646
======================= Grid point 36 (16/63) =======================
q-point: ( 0.30  0.30  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.577   (   5.835    5.835    0.000)    8.252
   0.894   (   5.759    5.759    0.000)    8.145
   1.138   (   2.783    2.783    0.000)    3.935
   1.324   (  -1.577   -1.577    0.000)    2.230
   1.351   (  -2.267   -2.267    0.000)    3.206
   1.512   (  -7.136   -7.136    0.000)   10.092
   1.601   (   1.400    1.400    0.000)    1.980
   1.634   (   0.911    0.911    0.000)    1.288
   1.716   (  -1.050   -1.050    0.000)    1.485
   2.079   (   3.763    3.763    0.000)    5.322
   2.109   (  -3.054   -3.054    0.000)    4.319
   2.181   (  -0.539   -0.539    0.000)    0.762
   2.226   (   0.735    0.735    0.000)    1.039
   2.468   (   1.481    1.481    0.000)    2.095
   2.514   ( -18.079  -18.079    0.000)   25.568
   2.746   (  -4.092   -4.092    0.000)    5.787
   3.332   (  13.411   13.411    0.000)   18.966
   3.684   (   5.000    5.000    0.000)    7.071
   5.923   (   3.190    3.190    0.000)    4.511
   6.014   (  -0.106   -0.106    0.000)    0.150
   7.784   (  26.622   26.622    0.000)   37.649
   8.683   (   8.580    8.580    0.000)   12.134
   9.930   (  -1.287   -1.287    0.000)    1.820
  10.192   (  -4.749   -4.749    0.000)    6.716
======================= Grid point 37 (17/63) =======================
q-point: ( 0.40  0.30  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.778   (  12.302   -4.506    0.000)   13.101
   0.995   (   3.827    2.124    0.000)    4.377
   1.225   (   5.308   -2.167    0.000)    5.734
   1.296   (  -3.028   -0.962    0.000)    3.177
   1.366   (  -6.361   -8.084    0.000)   10.286
   1.400   (   8.535  -11.907    0.000)   14.650
   1.614   (  -6.935   -6.563    0.000)    9.548
   1.629   (   0.869   -0.724    0.000)    1.132
   1.642   (  -0.117    0.609    0.000)    0.620
   2.083   ( -15.304  -15.012    0.000)   21.437
   2.099   (   1.789    2.186    0.000)    2.825
   2.152   (   3.280    3.183    0.000)    4.571
   2.213   (  -6.054   -4.823    0.000)    7.740
   2.237   (  -0.230    0.444    0.000)    0.500
   2.512   (   2.959   -1.038    0.000)    3.136
   2.663   (  -4.293   -3.732    0.000)    5.688
   3.558   (   9.546   10.967    0.000)   14.540
   3.741   (   0.514    6.676    0.000)    6.696
   5.987   (   2.765    3.261    0.000)    4.275
   6.019   (   0.451    1.538    0.000)    1.603
   8.283   (  23.095   24.507    0.000)   33.674
   8.767   (  -0.515   14.699    0.000)   14.708
   9.919   (   0.581   -1.413    0.000)    1.527
  10.078   (  -6.441   -2.928    0.000)    7.076
======================= Grid point 38 (18/63) =======================
q-point: (-0.50  0.30  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.981   (  -6.645   -4.605    0.000)    8.085
   0.981   (   6.645   -4.605    0.000)    8.085
   1.250   (  -0.548   -1.582    0.000)    1.674
   1.250   (   0.548   -1.582    0.000)    1.674
   1.407   (  -7.867   -7.051    0.000)   10.564
   1.407   (   7.867   -7.051    0.000)   10.564
   1.626   (  -2.432   -1.001    0.000)    2.630
   1.626   (   2.432   -1.001    0.000)    2.630
   1.734   ( -15.535  -18.056    0.000)   23.820
   1.734   (  15.535  -18.056    0.000)   23.820
   2.147   (  -2.538    2.076    0.000)    3.279
   2.147   (   2.538    2.076    0.000)    3.279
   2.213   (  -2.080    1.504    0.000)    2.567
   2.213   (   2.080    1.504    0.000)    2.567
   2.581   (  -3.970   -2.795    0.000)    4.855
   2.581   (   3.970   -2.795    0.000)    4.855
   3.699   (  -4.762    8.620    0.000)    9.848
   3.699   (   4.762    8.620    0.000)    9.848
   6.020   (  -0.615    2.813    0.000)    2.880
   6.020   (   0.615    2.813    0.000)    2.880
   8.642   ( -12.632   20.365    0.000)   23.965
   8.642   (  12.632   20.365    0.000)   23.965
   9.966   (  -4.331   -1.946    0.000)    4.748
   9.966   (   4.331   -1.946    0.000)    4.748
======================= Grid point 48 (19/63) =======================
q-point: ( 0.40  0.40  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 93
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.773   (   4.007    4.007    0.000)    5.667
   0.999   (  -2.718   -2.718    0.000)    3.844
   1.211   (   0.996    0.996    0.000)    1.409
   1.249   (  -2.666   -2.666    0.000)    3.771
   1.252   (  -1.838   -1.838    0.000)    2.599
   1.271   (  -1.609   -1.609    0.000)    2.275
   1.449   (  -8.605   -8.605    0.000)   12.169
   1.639   (   0.600    0.600    0.000)    0.848
   1.649   (   0.077    0.077    0.000)    0.109
   1.701   ( -19.765  -19.765    0.000)   27.951
   2.143   (   2.177    2.177    0.000)    3.078
   2.195   (  -0.313   -0.313    0.000)    0.442
   2.210   (   2.500    2.500    0.000)    3.535
   2.252   (   0.537    0.537    0.000)    0.760
   2.505   (   0.458    0.458    0.000)    0.648
   2.591   (  -3.477   -3.477    0.000)    4.917
   3.735   (   6.993    6.993    0.000)    9.890
   3.831   (   2.306    2.306    0.000)    3.261
   6.051   (   2.698    2.698    0.000)    3.816
   6.057   (   1.921    1.921    0.000)    2.716
   8.718   (  18.734   18.734    0.000)   26.493
   8.975   (   5.964    5.964    0.000)    8.435
   9.899   (  -0.344   -0.344    0.000)    0.487
  10.002   (  -4.587   -4.587    0.000)    6.487
======================= Grid point 49 (20/63) =======================
q-point: (-0.50  0.40  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.894   (  -6.896   -3.503    0.000)    7.735
   0.894   (   6.896   -3.503    0.000)    7.735
   1.236   (  -1.614   -0.271    0.000)    1.637
   1.236   (   1.614   -0.271    0.000)    1.637
   1.284   (  -3.718   -4.927    0.000)    6.173
   1.284   (   3.718   -4.927    0.000)    6.173
   1.378   ( -10.362  -13.844    0.000)   17.292
   1.378   (  10.362  -13.844    0.000)   17.292
   1.648   (  -0.295    0.194    0.000)    0.353
   1.648   (   0.295    0.194    0.000)    0.353
   2.179   (  -1.393    1.186    0.000)    1.830
   2.179   (   1.393    1.186    0.000)    1.830
   2.247   (  -1.124    1.537    0.000)    1.904
   2.247   (   1.124    1.537    0.000)    1.904
   2.533   (  -2.330   -1.966    0.000)    3.048
   2.533   (   2.330   -1.966    0.000)    3.048
   3.829   (  -2.518    4.472    0.000)    5.132
   3.829   (   2.518    4.472    0.000)    5.132
   6.079   (  -0.013    2.612    0.000)    2.612
   6.079   (   0.013    2.612    0.000)    2.612
   8.974   (  -6.658   12.723    0.000)   14.360
   8.974   (   6.658   12.723    0.000)   14.360
   9.921   (  -2.907   -2.314    0.000)    3.716
   9.921   (   2.907   -2.314    0.000)    3.716
======================= Grid point 65 (21/63) =======================
q-point: (-0.50 -0.50  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 36
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.860   (   0.000    0.000    0.000)    0.000
   0.860   (   0.000    0.000    0.000)    0.000
   1.208   (   0.000    0.000    0.000)    0.000
   1.208   (   0.000    0.000    0.000)    0.000
   1.239   (  -0.000   -0.000    0.000)    0.000
   1.239   (  -0.000   -0.000    0.000)    0.000
   1.249   (   0.000    0.000    0.000)    0.000
   1.249   (   0.000    0.000    0.000)    0.000
   1.650   (   0.000    0.000    0.000)    0.000
   1.650   (   0.000    0.000    0.000)    0.000
   2.191   (   0.000    0.000    0.000)    0.000
   2.191   (   0.000    0.000    0.000)    0.000
   2.263   (   0.000    0.000    0.000)    0.000
   2.263   (   0.000    0.000    0.000)    0.000
   2.512   (   0.000    0.000    0.000)    0.000
   2.512   (   0.000    0.000    0.000)    0.000
   3.873   (   0.000    0.000    0.000)    0.000
   3.873   (   0.000    0.000    0.000)    0.000
   6.106   (   0.000    0.000    0.000)    0.000
   6.106   (   0.000    0.000    0.000)    0.000
   9.105   (   0.000    0.000    0.000)    0.000
   9.105   (   0.000    0.000    0.000)    0.000
   9.894   (   0.000    0.000    0.000)    0.000
   9.894   (   0.000    0.000    0.000)    0.000
======================= Grid point 121 (22/63) =======================
q-point: ( 0.00  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.208   (   0.000    0.000    9.731)    9.731
   0.208   (   0.000    0.000    9.731)    9.731
   0.583   (   0.000   -0.000   27.857)   27.857
   0.768   (   0.000    0.000   -3.606)    3.606
   0.768   (   0.000    0.000   -3.606)    3.606
   1.341   (   0.000    0.000    0.704)    0.704
   1.341   (   0.000    0.000    0.704)    0.704
   1.382   (  -0.000   -0.000    0.453)    0.453
   1.711   (   0.000    0.000    1.273)    1.273
   2.099   (   0.000   -0.000  -14.440)   14.440
   2.151   (   0.000    0.000   -0.484)    0.484
   2.151   (  -0.000   -0.000   -0.484)    0.484
   2.167   (   0.000    0.000    0.209)    0.209
   2.167   (   0.000    0.000    0.209)    0.209
   2.663   (  -0.000   -0.000    8.452)    8.452
   3.058   (  -0.000   -0.000   -4.541)    4.541
   3.273   (   0.000    0.000    0.051)    0.051
   3.273   (   0.000    0.000    0.051)    0.051
   6.024   (  -0.000    0.000    0.172)    0.172
   6.024   (   0.000    0.000    0.172)    0.172
   6.051   (   0.000   -0.000   -0.104)    0.104
   6.051   (   0.000    0.000   -0.104)    0.104
  10.605   (   0.000    0.000    1.636)    1.636
  10.947   (  -0.000   -0.000   -0.351)    0.351
======================= Grid point 122 (23/63) =======================
q-point: ( 0.10  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.474   (  18.431    0.000    4.308)   18.927
   0.483   (  14.074    0.000    9.004)   16.708
   0.648   (  10.137    0.000   19.578)   22.047
   0.835   (   6.423    0.000   -3.063)    7.115
   0.892   (  11.030    0.000   -2.573)   11.326
   1.405   (   2.405    0.000    0.139)    2.409
   1.433   (   9.084    0.000    0.453)    9.095
   1.484   (  12.538    0.000   -1.246)   12.600
   1.676   (  -3.024    0.000    1.105)    3.219
   1.810   (  -8.988    0.000   -0.262)    8.992
   2.136   (  -1.520    0.000   -0.446)    1.584
   2.154   (  -1.295    0.000    0.286)    1.326
   2.165   (   0.290    0.000    0.524)    0.599
   2.381   (   6.633    0.000   -9.430)   11.529
   2.618   (  -2.115    0.000    6.815)    7.135
   2.988   (  -3.318    0.000   -2.641)    4.241
   3.252   (  -2.106    0.000    0.020)    2.106
   3.327   (   4.576    0.000    0.144)    4.578
   6.028   (   0.394    0.000    0.143)    0.419
   6.052   (   0.066    0.000   -0.091)    0.112
   6.138   (   8.866    0.000    0.009)    8.867
   6.979   (  51.617    0.000  -40.944)   65.884
  10.566   (  -4.214    0.000    1.833)    4.595
  10.587   ( -22.387    0.000   17.492)   28.410
======================= Grid point 123 (24/63) =======================
q-point: ( 0.20  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 156
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.652   (   7.312    0.000   15.200)   16.867
   0.817   (  16.044    0.000    2.095)   16.180
   0.955   (  15.903    0.000    6.123)   17.041
   0.997   (   9.479    0.000   -1.949)    9.677
   1.148   (  13.967    0.000   -3.247)   14.340
   1.470   (   4.185    0.000   -0.284)    4.195
   1.612   (  -3.270    0.000    2.607)    4.182
   1.673   (  14.350    0.000    0.099)   14.350
   1.707   (  -1.956    0.000   -3.901)    4.364
   1.798   (  17.762    0.000   -1.349)   17.813
   2.091   (  -3.057    0.000   -0.278)    3.070
   2.122   (  -1.509    0.000    0.380)    1.556
   2.170   (   0.158    0.000    1.619)    1.627
   2.463   (   2.529    0.000   -5.990)    6.502
   2.601   (   0.267    0.000    3.741)    3.751
   2.929   (  -3.094    0.000   -3.033)    4.333
   3.189   (  -4.393    0.000   -0.049)    4.393
   3.418   (   3.797    0.000    0.147)    3.800
   6.038   (   0.615    0.000    0.077)    0.620
   6.054   (   0.148    0.000   -0.058)    0.159
   6.394   (  17.049    0.000   -0.020)   17.049
   7.632   (  19.249    0.000  -24.079)   30.827
  10.262   ( -11.228    0.000   14.625)   18.438
  10.460   (  -6.821    0.000    0.817)    6.870
======================= Grid point 124 (25/63) =======================
q-point: ( 0.30  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.809   (   8.580    0.000   11.265)   14.160
   1.097   (  12.602    0.000    0.856)   12.631
   1.186   (   9.365    0.000   -0.945)    9.413
   1.195   (   7.719    0.000    3.707)    8.563
   1.392   (   9.620    0.000   -5.225)   10.948
   1.561   (  -0.627    0.000    2.514)    2.591
   1.572   (   5.588    0.000    0.573)    5.617
   1.698   (   1.866    0.000   -3.156)    3.666
   1.908   (   5.634    0.000   -0.070)    5.634
   2.018   (  -4.260    0.000   -0.054)    4.261
   2.134   (   9.353    0.000   -0.300)    9.357
   2.153   (   8.052    0.000    0.368)    8.060
   2.185   (   6.450    0.000    2.560)    6.939
   2.480   (  -0.769    0.000   -3.785)    3.862
   2.651   (   5.887    0.000    0.921)    5.959
   2.860   (  -3.966    0.000   -3.077)    5.020
   3.078   (  -7.004    0.000   -0.098)    7.005
   3.446   (  -1.440    0.000    0.125)    1.446
   6.050   (   0.578    0.000    0.017)    0.578
   6.057   (   0.210    0.000   -0.021)    0.211
   6.790   (  22.744    0.000   -0.116)   22.744
   7.853   (   4.916    0.000   -9.435)   10.639
  10.102   (  -5.778    0.000    8.551)   10.320
  10.308   (  -8.394    0.000    0.703)    8.423
======================= Grid point 125 (26/63) =======================
q-point: ( 0.40  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 156
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.979   (   8.626    0.000    7.883)   11.686
   1.238   (  -2.636    0.000    4.056)    4.838
   1.306   (   8.392    0.000    0.243)    8.396
   1.345   (   6.322    0.000   -0.268)    6.328
   1.503   (   2.914    0.000   -4.328)    5.218
   1.561   (  -0.088    0.000   -2.357)    2.358
   1.701   (   6.987    0.000    1.200)    7.089
   1.789   (   5.289    0.000    0.244)    5.295
   1.907   (  -2.379    0.000   -0.085)    2.381
   1.933   (  -4.101    0.000    0.044)    4.101
   2.157   (   0.912    0.000    3.081)    3.213
   2.321   (  -1.682    0.000    1.995)    2.610
   2.411   (  14.961    0.000    0.157)   14.962
   2.510   (   6.112    0.000   -4.891)    7.829
   2.751   (  -3.868    0.000   -4.037)    5.591
   2.869   (  12.284    0.000    0.636)   12.300
   2.913   (  -9.953    0.000   -0.077)    9.953
   3.353   (  -8.043    0.000    0.181)    8.046
   6.059   (   0.334    0.000   -0.008)    0.334
   6.061   (   0.166    0.000    0.001)    0.166
   7.247   (  22.775    0.000    0.223)   22.776
   7.847   (  -5.514    0.000   -2.337)    5.989
  10.026   (  -1.861    0.000    4.645)    5.004
  10.148   (  -7.463    0.000    1.120)    7.546
======================= Grid point 126 (27/63) =======================
q-point: (-0.50  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.136   (  -7.039    0.000    5.666)    9.036
   1.136   (   7.039    0.000    5.666)    9.036
   1.404   (  -1.012    0.000    0.050)    1.014
   1.404   (   1.012    0.000    0.050)    1.014
   1.544   (  -1.481    0.000   -3.671)    3.959
   1.544   (   1.481    0.000   -3.671)    3.959
   1.813   (  -3.230    0.000    1.147)    3.427
   1.813   (   3.230    0.000    1.147)    3.427
   1.881   (  -0.592    0.000   -0.062)    0.595
   1.881   (   0.592    0.000   -0.062)    0.595
   2.225   (  -5.688    0.000    2.329)    6.146
   2.225   (   5.688    0.000    2.329)    6.146
   2.642   (  -6.399    0.000   -4.283)    7.700
   2.642   (   6.399    0.000   -4.283)    7.700
   2.688   ( -12.968    0.000    0.019)   12.968
   2.688   (  12.968    0.000    0.019)   12.968
   3.140   ( -13.266    0.000    0.312)   13.270
   3.140   (  13.266    0.000    0.312)   13.270
   6.063   (  -0.039    0.000   -0.003)    0.039
   6.063   (   0.039    0.000   -0.003)    0.039
   7.636   ( -15.928    0.000    0.180)   15.929
   7.636   (  15.928    0.000    0.180)   15.929
  10.038   (  -3.245    0.000    2.266)    3.958
  10.038   (   3.245    0.000    2.266)    3.958
======================= Grid point 133 (28/63) =======================
q-point: ( 0.10  0.10  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 140
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.242   (   2.369    2.369    9.544)   10.115
   0.610   (   3.863    3.863   20.306)   21.028
   0.819   (   2.518    2.518   -4.135)    5.457
   0.904   (  16.834   16.834    4.426)   24.215
   1.046   (  11.056   11.056   -1.776)   15.736
   1.357   (   0.644    0.644    0.793)    1.207
   1.429   (   2.170    2.170    0.120)    3.071
   1.650   (  -2.270   -2.270    1.629)    3.600
   1.709   (   8.225    8.225   -1.881)   11.784
   1.774   (   3.990    3.990   -2.718)    6.263
   2.129   (  -1.429   -1.429    1.188)    2.345
   2.158   (   0.317    0.317   -0.537)    0.700
   2.181   (   0.638    0.638    0.420)    0.995
   2.398   (   1.721    1.721   -7.675)    8.051
   2.595   (  -1.136   -1.136    5.668)    5.891
   2.960   (  -2.242   -2.242   -2.874)    4.280
   3.270   (  -1.187   -1.187   -0.063)    1.680
   3.330   (   2.696    2.696    0.336)    3.827
   5.955   (  -2.838   -2.838    0.698)    4.074
   6.067   (   0.394    0.394    0.328)    0.647
   6.216   (   9.062    9.062   -0.320)   12.819
   7.408   (  28.906   28.906  -36.219)   54.616
  10.424   ( -12.300  -12.300   18.258)   25.218
  10.533   (  -3.369   -3.369    1.077)    4.885
======================= Grid point 134 (29/63) =======================
q-point: ( 0.20  0.10  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.523   (  18.141  -13.548    4.837)   23.153
   0.715   (   7.147    2.170   16.527)   18.137
   0.938   (   8.763   -2.460   -3.766)    9.850
   1.189   (  11.774   12.686    1.702)   17.391
   1.274   (  11.336    9.485   -1.556)   14.862
   1.457   (   7.787   -5.470    0.517)    9.531
   1.502   (   6.606   -2.302   -0.003)    6.996
   1.599   (  -2.500   -1.089    2.972)    4.033
   1.692   (  -1.424   -1.185   -4.685)    5.038
   2.013   (   6.221    6.019   -3.176)    9.221
   2.109   (  -0.398    4.563    1.024)    4.693
   2.167   (   4.821    5.467    0.805)    7.333
   2.182   (   0.031    1.068    1.173)    1.587
   2.432   (   1.119   -3.511   -4.202)    5.588
   2.609   (   3.358    2.522    2.499)    4.887
   2.908   (  -3.093   -2.027   -3.141)    4.852
   3.189   (  -6.311   -1.667   -0.056)    6.528
   3.416   (   4.403    1.024    0.318)    4.532
   5.944   (   0.988   -6.322    0.713)    6.438
   6.062   (  -0.627    0.292    0.307)    0.757
   6.493   (  18.699   10.544   -0.307)   21.469
   7.845   (  15.345   17.422  -20.261)   30.814
  10.214   (  -8.531   -4.251   13.026)   16.141
  10.439   (  -6.115   -2.111    0.731)    6.510
======================= Grid point 135 (30/63) =======================
q-point: ( 0.30  0.10  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.811   (  10.380   -2.141    7.429)   12.943
   0.915   (  13.229   -8.081    7.428)   17.189
   1.130   (  10.174   -3.881   -2.777)   11.238
   1.339   (   2.959    6.198    1.535)    7.038
   1.451   (   6.125    4.422   -3.461)    8.310
   1.565   (   1.240    0.385    0.966)    1.618
   1.594   (   5.446    2.377    0.889)    6.009
   1.687   (   1.083   -1.051   -2.182)    2.653
   1.737   (  14.528  -13.626    0.263)   19.920
   2.009   (  -1.527   -1.353    0.236)    2.053
   2.062   (  -2.582    1.987   -0.721)    3.337
   2.173   (  -0.711    0.936    1.902)    2.236
   2.355   (   7.276    5.727    0.715)    9.287
   2.431   (  -0.329   -4.844   -3.176)    5.802
   2.754   (  10.278   10.032   -0.491)   14.371
   2.843   (  -2.418   -0.782   -2.553)    3.602
   3.029   (  -9.979   -5.884   -0.057)   11.585
   3.470   (   0.555    3.100    0.256)    3.160
   5.974   (   1.773   -5.656    0.455)    5.944
   6.049   (  -0.662   -0.835    0.172)    1.079
   6.923   (  24.417   13.269   -0.217)   27.790
   8.034   (   4.448   15.907   -8.583)   18.615
  10.084   (  -5.001   -1.759    7.978)    9.579
  10.300   (  -7.921   -0.905    0.692)    8.002
======================= Grid point 136 (31/63) =======================
q-point: ( 0.40  0.10  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.987   (   8.226    0.545    7.041)   10.842
   1.137   (   6.976   -7.183    4.050)   10.801
   1.304   (   6.041   -2.440   -1.419)    6.668
   1.348   (   1.234    0.702    2.201)    2.619
   1.525   (   2.242    1.838   -3.305)    4.396
   1.572   (  -0.398    0.733   -2.949)    3.064
   1.692   (   2.343   -2.466   -0.237)    3.410
   1.710   (   0.343   -4.407    1.539)    4.681
   1.926   (   3.059    0.325    0.447)    3.108
   1.993   (  -0.129    4.477   -0.159)    4.481
   2.148   (   9.444   -8.969    0.966)   13.060
   2.167   (   3.662   -3.560    1.974)    5.475
   2.350   (  -2.814   -3.833    1.052)    4.870
   2.518   (   6.076   -0.522   -3.983)    7.284
   2.737   (  -5.331   -2.041   -3.930)    6.931
   2.801   ( -13.114  -10.693   -0.047)   16.921
   3.035   (  13.756   13.660    0.379)   19.390
   3.427   (  -4.973    6.864    0.247)    8.479
   6.005   (   1.339   -4.055    0.230)    4.277
   6.037   (  -0.544   -1.898    0.108)    1.977
   7.407   (  23.861   15.436    0.101)   28.419
   8.025   (  -5.526   16.290   -2.468)   17.378
  10.018   (  -1.537   -0.836    4.521)    4.848
  10.144   (  -7.374   -0.355    1.104)    7.465
======================= Grid point 137 (32/63) =======================
q-point: (-0.50  0.10  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.128   (  -5.583   -0.818    5.295)    7.738
   1.128   (   5.583   -0.818    5.295)    7.738
   1.381   (  -1.938   -2.123    0.441)    2.908
   1.381   (   1.938   -2.123    0.441)    2.908
   1.555   (  -1.079    0.582   -3.161)    3.390
   1.555   (   1.079    0.582   -3.161)    3.390
   1.701   (  -0.932   -6.783    0.351)    6.855
   1.701   (   0.932   -6.783    0.351)    6.855
   1.975   (  -2.399    5.308    0.433)    5.841
   1.975   (   2.399    5.308    0.433)    5.841
   2.213   (  -4.331   -1.097    2.113)    4.942
   2.213   (   4.331   -1.097    2.113)    4.942
   2.505   ( -15.265  -14.419   -0.439)   21.003
   2.505   (  15.265  -14.419   -0.439)   21.003
   2.636   (  -5.698   -1.848   -3.679)    7.030
   2.636   (   5.698   -1.848   -3.679)    7.030
   3.276   ( -10.312   10.990    0.267)   15.073
   3.276   (  10.312   10.990    0.267)   15.073
   6.025   (  -0.754   -2.829    0.121)    2.930
   6.025   (   0.754   -2.829    0.121)    2.930
   7.810   ( -16.330   16.418   -0.079)   23.157
   7.810   (  16.330   16.418   -0.079)   23.157
  10.034   (  -3.375   -0.422    2.247)    4.076
  10.034   (   3.375   -0.422    2.247)    4.076
======================= Grid point 145 (33/63) =======================
q-point: ( 0.20  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 141
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.409   (   5.655    5.655    6.957)   10.600
   0.796   (   5.479    5.479   14.502)   16.442
   0.948   (   3.698    3.698   -4.904)    7.169
   1.372   (   5.297    5.297   -0.777)    7.531
   1.373   (  -0.037   -0.037    1.067)    1.068
   1.465   (   8.218    8.218    0.264)   11.624
   1.526   (   2.468    2.468    0.035)    3.491
   1.589   (   0.977    0.977    1.945)    2.386
   1.673   (  -0.575   -0.575   -3.444)    3.539
   2.010   (  -1.654   -1.654   -0.654)    2.429
   2.169   (   0.167    0.167   -0.789)    0.824
   2.209   (   0.717    0.717    0.684)    1.223
   2.307   (   2.774    2.774   -3.595)    5.321
   2.355   (  -0.752   -0.752    2.237)    2.477
   2.754   (  11.542   11.542   -0.169)   16.324
   2.858   (  -2.340   -2.340   -2.943)    4.428
   3.075   (  -9.654   -9.654   -0.037)   13.653
   3.483   (   4.974    4.974    0.467)    7.049
   5.874   (  -0.536   -0.536    1.132)    1.363
   6.053   (  -1.092   -1.092    0.473)    1.615
   6.808   (  21.135   21.135   -0.399)   29.892
   8.171   (  14.059   14.059  -12.731)   23.609
  10.122   (  -4.726   -4.726    9.678)   11.762
  10.377   (  -4.297   -4.297    0.654)    6.111
======================= Grid point 146 (34/63) =======================
q-point: ( 0.30  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.660   (  15.850   -7.279    4.571)   18.030
   0.919   (   6.849    3.747   11.435)   13.845
   1.071   (   8.308   -1.632   -4.412)    9.547
   1.393   (  -1.630   -0.010    1.518)    2.228
   1.454   (   7.086   -7.900    0.301)   10.616
   1.561   (   3.254    1.030   -2.157)    4.038
   1.581   (   2.132    0.509   -0.883)    2.363
   1.628   (  -0.279   -1.965   -0.822)    2.148
   1.687   (   2.856    2.846    1.903)    4.459
   2.002   (   1.451    1.089   -0.819)    1.991
   2.122   (  -4.226    3.182   -0.495)    5.313
   2.192   (  -1.174    0.677    1.312)    1.886
   2.299   (  -0.504   -6.918    0.198)    6.939
   2.417   (   3.780    1.193   -1.858)    4.378
   2.773   (  -4.108   -3.457   -4.067)    6.736
   2.833   ( -14.653  -13.733    0.043)   20.082
   3.048   (  14.768   15.319    0.287)   21.280
   3.570   (   3.170    6.473    0.420)    7.220
   5.902   (   2.668   -0.957    0.766)    2.936
   6.029   (  -1.189   -1.018    0.312)    1.597
   7.286   (  26.577   22.959   -0.235)   35.121
   8.369   (   6.035   16.473   -6.421)   18.682
  10.040   (  -3.485   -2.600    6.606)    7.909
  10.270   (  -6.565   -2.270    0.767)    6.988
======================= Grid point 147 (35/63) =======================
q-point: ( 0.40  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 252
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.942   (  11.891   -6.475    3.195)   13.911
   1.053   (   6.330    0.375    8.528)   10.627
   1.242   (   7.808   -3.465   -3.134)    9.099
   1.347   (  -1.873   -0.779    2.828)    3.480
   1.527   (  -2.750   -5.606   -0.231)    6.249
   1.561   (  -1.042   -1.300   -3.218)    3.624
   1.604   (   0.162   -0.842   -0.675)    1.092
   1.680   (  -0.813   -3.133    1.084)    3.414
   1.774   (  12.980  -12.465    0.229)   17.997
   2.050   (   2.583    1.722   -0.288)    3.118
   2.081   (   0.607    4.082   -0.092)    4.128
   2.171   (  -0.572    0.310    1.980)    2.084
   2.273   (  -2.205   -3.319    1.010)    4.111
   2.495   (   1.488   -3.299   -3.132)    4.787
   2.517   ( -14.446  -15.322   -0.243)   21.060
   2.686   (  -4.711   -3.423   -4.084)    7.113
   3.314   (  11.958   13.893    0.208)   18.332
   3.589   (  -1.591    8.822    0.348)    8.971
   5.955   (   2.410   -0.389    0.307)    2.461
   6.010   (  -0.844   -0.586    0.149)    1.038
   7.801   (  24.766   23.511   -0.062)   34.148
   8.396   (  -3.773   19.701   -2.384)   20.200
   9.994   (  -0.822   -1.553    4.106)    4.466
  10.131   (  -7.050   -1.112    1.164)    7.232
======================= Grid point 148 (36/63) =======================
q-point: (-0.50  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.098   (  -3.451   -2.475    4.499)    6.187
   1.098   (   3.451   -2.475    4.499)    6.187
   1.330   (  -3.309   -2.821    1.826)    4.716
   1.330   (   3.309   -2.821    1.826)    4.716
   1.510   (  -3.945   -5.893   -2.464)    7.507
   1.510   (   3.945   -5.893   -2.464)    7.507
   1.624   (  -2.093   -0.269   -0.108)    2.113
   1.624   (   2.093   -0.269   -0.108)    2.113
   2.042   (  -7.471   -2.527    0.398)    7.897
   2.042   (   7.471   -2.527    0.398)    7.897
   2.152   (  -8.721   -8.445    0.361)   12.145
   2.152   (   8.721   -8.445    0.361)   12.145
   2.231   (  -1.546   -5.144    0.986)    5.461
   2.231   (   1.546   -5.144    0.986)    5.461
   2.591   (  -4.876   -2.683   -3.900)    6.797
   2.591   (   4.876   -2.683   -3.900)    6.797
   3.503   (  -7.124   11.411    0.265)   13.455
   3.503   (   7.124   11.411    0.265)   13.455
   5.991   (  -1.250   -0.193    0.106)    1.269
   5.991   (   1.250   -0.193    0.106)    1.269
   8.206   ( -15.662   22.387   -0.435)   27.325
   8.206   (  15.662   22.387   -0.435)   27.325
  10.021   (  -3.685   -1.060    2.196)    4.419
  10.021   (   3.685   -1.060    2.196)    4.419
======================= Grid point 157 (37/63) =======================
q-point: ( 0.30  0.30  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 140
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.636   (   5.489    5.489    5.356)    9.432
   1.004   (   4.553    4.553    9.659)   11.609
   1.080   (   2.751    2.751   -5.056)    6.380
   1.343   (  -1.450   -1.450    1.340)    2.450
   1.370   (  -2.233   -2.233    1.720)    3.596
   1.539   (  -6.068   -6.068    1.729)    8.753
   1.579   (  -0.203   -0.203   -4.721)    4.729
   1.606   (   1.518    1.518    0.449)    2.194
   1.724   (  -0.717   -0.717    0.732)    1.250
   2.047   (   2.884    2.884   -1.094)    4.223
   2.158   (  -1.613   -1.613    2.269)    3.218
   2.168   (  -0.407   -0.407   -0.924)    1.088
   2.235   (   0.592    0.592    0.602)    1.031
   2.437   (   1.221    1.221   -2.411)    2.966
   2.516   ( -18.062  -18.062    0.140)   25.544
   2.696   (  -4.200   -4.200   -4.336)    7.355
   3.334   (  13.371   13.371    0.136)   18.910
   3.690   (   4.977    4.977    0.439)    7.052
   5.928   (   2.869    2.869    0.394)    4.076
   6.016   (  -0.236   -0.236    0.185)    0.382
   7.782   (  26.660   26.660   -0.157)   37.703
   8.635   (  10.087   10.087   -3.936)   14.799
   9.990   (  -2.312   -2.312    5.024)    5.994
  10.204   (  -4.601   -4.601    0.972)    6.578
======================= Grid point 158 (38/63) =======================
q-point: ( 0.40  0.30  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.825   (  11.597   -4.005    4.314)   13.005
   1.077   (   2.020    0.904    7.117)    7.453
   1.178   (   6.727   -2.696   -4.110)    8.331
   1.321   (  -2.079   -1.969    2.331)    3.693
   1.381   (  -1.066   -7.470    1.195)    7.640
   1.438   (   2.159   -9.890    3.062)   10.576
   1.545   (  -3.470   -4.978   -4.767)    7.717
   1.630   (   0.984    0.851    0.195)    1.315
   1.670   (  -2.779   -1.898    0.821)    3.464
   2.082   ( -13.441  -13.013   -0.071)   18.708
   2.100   (   1.481    1.367    0.045)    2.016
   2.167   (   0.958    1.807    1.052)    2.300
   2.211   (  -5.295   -5.128    0.023)    7.371
   2.245   (  -0.311    0.204    0.599)    0.705
   2.473   (   2.555   -1.077   -3.314)    4.321
   2.613   (  -4.182   -3.860   -4.347)    7.161
   3.560   (   9.585   10.952    0.151)   14.555
   3.746   (   0.456    6.663    0.371)    6.689
   5.986   (   2.558    2.968   -0.081)    3.919
   6.018   (   0.307    1.350   -0.060)    1.386
   8.281   (  23.021   24.455   -0.163)   33.586
   8.745   (   0.574   15.142   -1.837)   15.264
   9.960   (  -0.346   -1.842    3.457)    3.932
  10.094   (  -6.156   -2.829    1.330)    6.905
======================= Grid point 159 (39/63) =======================
q-point: (-0.50  0.30  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.025   (  -7.653   -4.823    3.947)    9.870
   1.025   (   7.653   -4.823    3.947)    9.870
   1.284   (  -3.816   -1.880    1.810)    4.623
   1.284   (   3.816   -1.880    1.810)    4.623
   1.384   (  -6.965   -5.996   -0.863)    9.230
   1.384   (   6.965   -5.996   -0.863)    9.230
   1.630   (  -2.216   -0.587    0.173)    2.299
   1.630   (   2.216   -0.587    0.173)    2.299
   1.734   ( -15.819  -18.350    0.006)   24.227
   1.734   (  15.819  -18.350    0.006)   24.227
   2.145   (  -2.171    1.944   -0.107)    2.916
   2.145   (   2.171    1.944   -0.107)    2.916
   2.222   (  -1.867    1.173    0.845)    2.361
   2.222   (   1.867    1.173    0.845)    2.361
   2.535   (  -3.671   -2.880   -4.000)    6.146
   2.535   (   3.671   -2.880   -4.000)    6.146
   3.703   (  -4.833    8.608    0.256)    9.876
   3.703   (   4.833    8.608    0.256)    9.876
   6.018   (  -0.646    2.563   -0.225)    2.652
   6.018   (   0.646    2.563   -0.225)    2.652
   8.635   ( -12.158   20.362   -0.560)   23.722
   8.635   (  12.158   20.362   -0.560)   23.722
   9.991   (  -3.698   -2.008    2.112)    4.708
   9.991   (   3.698   -2.008    2.112)    4.708
======================= Grid point 169 (40/63) =======================
q-point: ( 0.40  0.40  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 141
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.820   (   3.723    3.723    4.337)    6.821
   1.051   (  -4.152   -4.152    4.569)    7.440
   1.157   (   1.325    1.325   -4.170)    4.571
   1.273   (  -1.956   -1.956    1.009)    2.945
   1.287   (  -0.490   -0.490    4.048)    4.107
   1.320   (  -3.218   -3.218    2.690)    5.286
   1.420   (  -6.898   -6.898   -1.682)    9.899
   1.646   (   0.554    0.554    0.605)    0.990
   1.653   (  -0.168   -0.168    0.102)    0.259
   1.703   ( -19.787  -19.787    0.164)   27.983
   2.144   (   2.044    2.044    0.097)    2.892
   2.188   (  -0.167   -0.167   -0.608)    0.653
   2.218   (   2.170    2.170    0.664)    3.141
   2.256   (   0.437    0.437    0.322)    0.697
   2.463   (   0.181    0.181   -3.670)    3.679
   2.540   (  -3.419   -3.419   -4.398)    6.536
   3.737   (   7.020    7.020    0.165)    9.929
   3.835   (   2.262    2.262    0.320)    3.215
   6.046   (   2.519    2.519   -0.506)    3.599
   6.052   (   1.728    1.728   -0.395)    2.475
   8.715   (  18.652   18.652   -0.264)   26.379
   8.962   (   6.415    6.415   -1.076)    9.136
   9.931   (  -0.800   -0.800    2.706)    2.933
  10.020   (  -4.418   -4.418    1.563)    6.440
======================= Grid point 170 (41/63) =======================
q-point: (-0.50  0.40  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.935   (  -6.784   -3.542    3.801)    8.545
   0.935   (   6.784   -3.542    3.801)    8.545
   1.233   (  -4.825   -3.241   -2.379)    6.281
   1.233   (   4.825   -3.241   -2.379)    6.281
   1.311   (  -2.969   -1.445    4.149)    5.303
   1.311   (   2.969   -1.445    4.149)    5.303
   1.378   ( -10.945  -14.083    0.121)   17.837
   1.378   (  10.945  -14.083    0.121)   17.837
   1.653   (  -0.116    0.358    0.337)    0.505
   1.653   (   0.116    0.358    0.337)    0.505
   2.176   (  -1.149    1.119   -0.273)    1.627
   2.176   (   1.149    1.119   -0.273)    1.627
   2.251   (  -1.052    1.358    0.343)    1.752
   2.251   (   1.052    1.358    0.343)    1.752
   2.485   (  -2.087   -2.063   -4.149)    5.082
   2.485   (   2.087   -2.063   -4.149)    5.082
   3.832   (  -2.568    4.461    0.235)    5.153
   3.832   (   2.568    4.461    0.235)    5.153
   6.071   (  -0.038    2.423   -0.615)    2.500
   6.071   (   0.038    2.423   -0.615)    2.500
   8.968   (  -6.423   12.777   -0.503)   14.309
   8.968   (   6.423   12.777   -0.503)   14.309
   9.945   (  -2.558   -2.374    2.006)    4.025
   9.945   (   2.558   -2.374    2.006)    4.025
======================= Grid point 186 (42/63) =======================
q-point: (-0.50 -0.50  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 54
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.902   (   0.000    0.000    3.780)    3.780
   0.902   (   0.000    0.000    3.780)    3.780
   1.209   (   0.000    0.000   -2.858)    2.858
   1.209   (   0.000    0.000   -2.858)    2.858
   1.209   (  -0.000   -0.000    0.066)    0.066
   1.209   (   0.000    0.000    0.066)    0.066
   1.305   (   0.000    0.000    4.863)    4.863
   1.305   (   0.000    0.000    4.863)    4.863
   1.657   (   0.000    0.000    0.540)    0.540
   1.657   (   0.000    0.000    0.540)    0.540
   2.187   (   0.000    0.000   -0.337)    0.337
   2.187   (   0.000    0.000   -0.337)    0.337
   2.266   (   0.000    0.000    0.185)    0.185
   2.266   (   0.000    0.000    0.185)    0.185
   2.463   (   0.000    0.000   -4.245)    4.245
   2.463   (   0.000    0.000   -4.245)    4.245
   3.875   (   0.000    0.000    0.225)    0.225
   3.875   (   0.000    0.000    0.225)    0.225
   6.097   (   0.000    0.000   -0.781)    0.781
   6.097   (   0.000    0.000   -0.781)    0.781
   9.100   (   0.000    0.000   -0.451)    0.451
   9.100   (   0.000    0.000   -0.451)    0.451
   9.917   (   0.000    0.000    1.953)    1.953
   9.917   (   0.000    0.000    1.953)    1.953
======================= Grid point 242 (43/63) =======================
q-point: ( 0.00  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.377   (  -0.000    0.000    5.791)    5.791
   0.377   (   0.000    0.000    5.791)    5.791
   0.682   (  -0.000   -0.000   -3.808)    3.808
   0.682   (   0.000    0.000   -3.808)    3.808
   1.118   (   0.000    0.000   23.833)   23.833
   1.354   (   0.000    0.000    0.398)    0.398
   1.354   (  -0.000    0.000    0.398)    0.398
   1.391   (   0.000    0.000    0.278)    0.278
   1.583   (  -0.000   -0.000  -16.138)   16.138
   1.916   (   0.000    0.000   -1.474)    1.474
   2.143   (   0.000    0.000   -0.216)    0.216
   2.143   (   0.000    0.000   -0.216)    0.216
   2.169   (   0.000   -0.000    0.045)    0.045
   2.169   (   0.000    0.000    0.045)    0.045
   2.810   (  -0.000   -0.000    6.045)    6.045
   2.944   (  -0.000    0.000   -5.862)    5.862
   3.274   (  -0.000    0.000    0.032)    0.032
   3.274   (  -0.000   -0.000    0.032)    0.032
   6.027   (   0.000   -0.000    0.132)    0.132
   6.027   (   0.000    0.000    0.132)    0.132
   6.049   (  -0.000    0.000   -0.089)    0.089
   6.049   (  -0.000    0.000   -0.089)    0.089
  10.638   (  -0.000    0.000    1.088)    1.088
  10.939   (   0.000    0.000   -0.296)    0.296
======================= Grid point 243 (44/63) =======================
q-point: ( 0.10  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.573   (  15.495    0.000    4.354)   16.096
   0.620   (  17.679    0.000    5.671)   18.566
   0.758   (   7.313    0.000   -3.684)    8.189
   0.782   (   9.793    0.000   -4.473)   10.767
   1.137   (   2.107    0.000   22.056)   22.156
   1.409   (   2.001    0.000    0.137)    2.006
   1.441   (   8.710    0.000    0.263)    8.714
   1.475   (  11.488    0.000   -0.224)   11.490
   1.551   (  -2.340    0.000  -14.748)   14.932
   1.794   (  -7.700    0.000    0.666)    7.729
   2.129   (  -1.440    0.000   -0.191)    1.453
   2.158   (  -1.167    0.000    0.091)    1.170
   2.175   (   0.837    0.000    0.231)    0.868
   2.261   (   7.146    0.000   -2.152)    7.463
   2.744   (  -4.787    0.000    4.718)    6.721
   2.909   (  -2.394    0.000   -3.999)    4.661
   3.253   (  -2.163    0.000    0.012)    2.163
   3.327   (   4.679    0.000   -0.063)    4.679
   6.031   (   0.335    0.000    0.108)    0.352
   6.050   (   0.097    0.000   -0.076)    0.123
   6.138   (   9.000    0.000    0.002)    9.000
   6.537   (  40.916    0.000   -8.151)   41.720
  10.592   (  -4.508    0.000    0.839)    4.585
  10.769   ( -14.635    0.000    3.182)   14.977
======================= Grid point 244 (45/63) =======================
q-point: ( 0.20  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 156
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.873   (  14.133    0.000    3.207)   14.492
   0.895   (  10.176    0.000    7.597)   12.699
   0.943   (  10.968    0.000   -3.119)   11.403
   1.044   (  14.504    0.000   -4.391)   15.154
   1.208   (   5.274    0.000   15.672)   16.535
   1.466   (   3.895    0.000   -0.196)    3.900
   1.525   (   0.451    0.000  -12.665)   12.673
   1.672   (  -4.544    0.000    1.272)    4.718
   1.675   (  14.100    0.000    0.060)   14.100
   1.778   (  18.107    0.000   -0.460)   18.113
   2.086   (  -2.884    0.000   -0.120)    2.886
   2.129   (  -1.404    0.000    0.183)    1.416
   2.204   (   2.172    0.000    1.201)    2.482
   2.363   (   3.176    0.000   -2.683)    4.158
   2.679   (  -1.601    0.000    3.054)    3.448
   2.856   (  -3.148    0.000   -3.161)    4.461
   3.188   (  -4.460    0.000   -0.031)    4.460
   3.419   (   3.815    0.000   -0.006)    3.815
   6.039   (   0.539    0.000    0.056)    0.542
   6.053   (   0.194    0.000   -0.044)    0.199
   6.393   (  16.855    0.000   -0.033)   16.855
   7.272   (  30.320    0.000   -9.205)   31.686
  10.463   ( -10.017    0.000    2.770)   10.393
  10.486   ( -10.925    0.000    3.079)   11.351
======================= Grid point 245 (46/63) =======================
q-point: ( 0.30  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.054   (   6.680    0.000   10.749)   12.655
   1.120   (  11.247    0.000    1.159)   11.306
   1.161   (  10.656    0.000   -1.213)   10.724
   1.260   (   7.262    0.000   -6.689)    9.873
   1.320   (   4.645    0.000    7.652)    8.951
   1.523   (  -0.436    0.000   -6.360)    6.375
   1.593   (   6.853    0.000   -0.559)    6.876
   1.649   (   4.145    0.000    0.306)    4.156
   1.906   (   5.626    0.000   -0.047)    5.626
   2.017   (  -4.096    0.000   -0.020)    4.096
   2.125   (  14.897    0.000   -0.428)   14.904
   2.160   (   7.888    0.000    0.216)    7.891
   2.256   (   3.512    0.000    3.645)    5.061
   2.391   (  -0.275    0.000   -4.001)    4.010
   2.691   (   3.281    0.000    2.522)    4.138
   2.784   (  -3.861    0.000   -3.662)    5.321
   3.076   (  -7.022    0.000   -0.060)    7.022
   3.448   (  -1.400    0.000    0.032)    1.400
   6.050   (   0.532    0.000    0.011)    0.532
   6.057   (   0.244    0.000   -0.013)    0.245
   6.787   (  22.859    0.000   -0.106)   22.860
   7.693   (  13.167    0.000   -4.650)   13.964
  10.257   (  -8.918    0.000    5.269)   10.358
  10.315   (  -8.550    0.000   -0.463)    8.563
======================= Grid point 246 (47/63) =======================
q-point: ( 0.40  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 156
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.173   (   5.685    0.000   10.002)   11.505
   1.311   (   7.944    0.000    0.226)    7.947
   1.335   (  -0.533    0.000    1.900)    1.973
   1.339   (   6.929    0.000   -0.239)    6.933
   1.379   (   3.117    0.000   -4.019)    5.086
   1.484   (  -2.319    0.000   -5.245)    5.735
   1.743   (   7.451    0.000    2.375)    7.821
   1.783   (   6.719    0.000   -1.040)    6.799
   1.905   (  -2.446    0.000   -0.061)    2.446
   1.933   (  -4.111    0.000    0.034)    4.112
   2.243   (   0.952    0.000    4.935)    5.026
   2.350   (  -3.192    0.000   -1.396)    3.484
   2.409   (   8.303    0.000   -2.742)    8.744
   2.414   (  14.807    0.000    0.097)   14.808
   2.670   (  -5.445    0.000   -2.919)    6.178
   2.876   (  12.762    0.000    0.141)   12.763
   2.911   (  -9.892    0.000   -0.047)    9.893
   3.356   (  -7.977    0.000    0.065)    7.977
   6.059   (   0.329    0.000   -0.004)    0.329
   6.061   (   0.172    0.000   -0.000)    0.172
   7.251   (  23.222    0.000    0.135)   23.222
   7.805   (  -1.490    0.000   -1.270)    1.958
  10.120   (  -4.839    0.000    3.685)    6.082
  10.161   (  -6.571    0.000   -0.278)    6.577
======================= Grid point 247 (48/63) =======================
q-point: (-0.50  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 90
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.275   (  -4.588    0.000    7.187)    8.527
   1.275   (   4.588    0.000    7.187)    8.527
   1.405   (  -0.845    0.000    0.031)    0.846
   1.405   (   0.845    0.000    0.031)    0.846
   1.438   (  -2.646    0.000   -6.214)    6.754
   1.438   (   2.646    0.000   -6.214)    6.754
   1.841   (  -1.509    0.000    1.110)    1.873
   1.841   (   1.509    0.000    1.110)    1.873
   1.879   (  -0.736    0.000   -0.039)    0.737
   1.879   (   0.736    0.000   -0.039)    0.737
   2.276   (  -3.206    0.000    2.018)    3.788
   2.276   (   3.206    0.000    2.018)    3.788
   2.551   (  -6.552    0.000   -3.379)    7.372
   2.551   (   6.552    0.000   -3.379)    7.372
   2.688   ( -12.850    0.000    0.011)   12.850
   2.688   (  12.850    0.000    0.011)   12.850
   3.145   ( -13.112    0.000    0.115)   13.113
   3.145   (  13.112    0.000    0.115)   13.113
   6.063   (  -0.047    0.000   -0.002)    0.047
   6.063   (   0.047    0.000   -0.002)    0.047
   7.639   ( -15.074    0.000    0.132)   15.075
   7.639   (  15.074    0.000    0.132)   15.075
  10.081   (  -1.144    0.000    1.377)    1.790
  10.081   (   1.144    0.000    1.377)    1.790
======================= Grid point 254 (49/63) =======================
q-point: ( 0.10  0.10  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 140
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.416   (   2.250    2.250    6.321)    7.077
   0.719   (   1.961    1.961   -4.487)    5.274
   0.878   (  10.480   10.480    5.342)   15.754
   0.984   (  14.477   14.477   -2.506)   20.626
   1.166   (   3.231    3.231   18.839)   19.385
   1.371   (   0.711    0.711    0.432)    1.095
   1.432   (   1.970    1.970    0.074)    2.787
   1.533   (  -1.239   -1.239  -14.144)   14.253
   1.701   (  13.581   13.581   -0.158)   19.207
   1.732   (  -3.635   -3.635    0.470)    5.161
   2.150   (   0.273    0.273   -0.221)    0.444
   2.154   (   0.292    0.292    0.829)    0.926
   2.187   (   0.779    0.779    0.149)    1.112
   2.283   (   2.499    2.499   -2.698)    4.446
   2.704   (  -2.942   -2.942    3.970)    5.751
   2.886   (  -2.008   -2.008   -3.364)    4.402
   3.269   (  -1.244   -1.244   -0.039)    1.760
   3.333   (   2.887    2.887    0.008)    4.082
   5.969   (  -2.446   -2.446    0.457)    3.489
   6.074   (   0.644    0.644    0.208)    0.934
   6.209   (   8.915    8.915   -0.224)   12.609
   6.942   (  31.735   31.735  -10.335)   46.055
  10.553   (  -3.869   -3.869    0.668)    5.512
  10.645   ( -10.910  -10.910    4.631)   16.110
======================= Grid point 255 (50/63) =======================
q-point: ( 0.20  0.10  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.641   (  16.008  -11.236    5.545)   20.328
   0.835   (   9.177   -3.884   -4.878)   11.095
   1.036   (   6.261    5.095   10.402)   13.167
   1.229   (   8.404    8.508   -0.604)   11.974
   1.282   (   8.712    8.780    7.380)   14.403
   1.465   (   6.689   -3.351    0.123)    7.483
   1.502   (   7.029   -3.636   -0.255)    7.918
   1.525   (   0.905   -0.012  -11.254)   11.290
   1.651   (  -2.455   -1.301    1.492)    3.153
   1.957   (   8.189    5.829   -1.662)   10.188
   2.118   (  -2.925    2.744    0.071)    4.012
   2.181   (   1.671    4.576    0.317)    4.881
   2.229   (   5.578    3.581    2.889)    7.231
   2.340   (   2.764   -1.656   -3.756)    4.949
   2.674   (   0.767    0.592    2.825)    2.987
   2.835   (  -3.253   -2.121   -3.138)    4.993
   3.188   (  -6.311   -1.702   -0.039)    6.536
   3.421   (   4.406    1.299    0.091)    4.595
   5.958   (   0.774   -5.471    0.466)    5.545
   6.068   (  -0.766    0.733    0.189)    1.077
   6.487   (  18.729   10.225   -0.212)   21.339
   7.527   (  25.072   21.529   -8.579)   34.142
  10.414   ( -10.811   -5.345    5.453)   13.236
  10.452   (  -6.799   -2.529    0.258)    7.259
======================= Grid point 256 (51/63) =======================
q-point: ( 0.30  0.10  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.935   (  13.024   -8.944    5.469)   16.719
   1.046   (  11.414   -6.237   -4.480)   13.757
   1.149   (   5.867    1.686   11.624)   13.130
   1.346   (   4.374    3.838   -4.072)    7.102
   1.407   (   2.885    6.265    2.181)    7.234
   1.511   (  -1.065   -0.935   -5.869)    6.038
   1.601   (   6.444    0.845   -0.256)    6.504
   1.673   (   4.391    2.374    0.726)    5.044
   1.742   (  14.435  -13.133    0.150)   19.516
   2.015   (   0.670   -2.163    0.181)    2.272
   2.049   (  -3.563    2.517   -0.451)    4.385
   2.216   (   1.692    0.632    1.670)    2.460
   2.337   (   2.120   -1.686   -1.708)    3.202
   2.401   (   3.544    3.473    0.093)    4.963
   2.740   (  -1.280    0.528   -0.768)    1.584
   2.797   (   7.276    7.103   -1.293)   10.250
   3.028   (  -9.967   -5.815   -0.036)   11.539
   3.474   (   0.548    3.309    0.111)    3.356
   5.983   (   1.507   -5.048    0.292)    5.276
   6.052   (  -0.770   -0.581    0.103)    0.970
   6.918   (  24.631   13.174   -0.156)   27.933
   7.885   (  11.494   17.129   -4.399)   21.092
  10.231   (  -7.832   -2.588    4.990)    9.640
  10.306   (  -7.968   -0.875   -0.384)    8.025
======================= Grid point 257 (52/63) =======================
q-point: ( 0.40  0.10  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.138   (   7.761   -3.880    6.562)   10.879
   1.237   (   6.507   -4.465   -1.026)    7.958
   1.276   (   6.909   -2.112    5.124)    8.857
   1.400   (  -0.002    1.327    2.046)    2.439
   1.428   (   1.892    2.279   -5.731)    6.451
   1.493   (  -0.541    1.997   -4.008)    4.511
   1.691   (   1.880   -5.330    0.079)    5.653
   1.744   (   1.030   -1.833    1.211)    2.426
   1.934   (   3.933    1.030    0.287)    4.076
   1.989   (  -0.933    4.440   -0.204)    4.542
   2.159   (  14.240  -14.061    0.269)   20.014
   2.236   (   0.584   -0.637    4.368)    4.452
   2.333   (  -1.601   -2.875   -3.636)    4.904
   2.466   (   3.888   -0.184   -0.057)    3.892
   2.651   (  -5.549   -2.091   -3.360)    6.816
   2.800   ( -13.195  -10.712   -0.020)   16.996
   3.041   (  13.623   13.473    0.151)   19.161
   3.431   (  -4.970    6.985    0.115)    8.573
   6.010   (   1.170   -3.752    0.149)    3.934
   6.039   (  -0.569   -1.770    0.057)    1.860
   7.409   (  24.177   15.256    0.061)   28.588
   7.980   (  -1.856   16.124   -1.378)   16.289
  10.110   (  -4.302   -1.139    3.516)    5.671
  10.158   (  -6.508   -0.288   -0.180)    6.517
======================= Grid point 258 (53/63) =======================
q-point: (-0.50  0.10  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.251   (  -3.757   -2.124    5.557)    7.036
   1.251   (   3.757   -2.124    5.557)    7.036
   1.383   (  -2.786   -1.580   -0.716)    3.282
   1.383   (   2.786   -1.580   -0.716)    3.282
   1.473   (  -2.067    2.159   -3.617)    4.692
   1.473   (   2.067    2.159   -3.617)    4.692
   1.715   (  -2.093   -7.140    0.656)    7.469
   1.715   (   2.093   -7.140    0.656)    7.469
   1.983   (  -1.231    5.222    0.276)    5.372
   1.983   (   1.231    5.222    0.276)    5.372
   2.261   (  -2.294   -1.464    1.934)    3.338
   2.261   (   2.294   -1.464    1.934)    3.338
   2.485   ( -13.771  -11.204   -1.421)   17.810
   2.485   (  13.771  -11.204   -1.421)   17.810
   2.569   (  -6.712   -4.867   -1.875)    8.500
   2.569   (   6.712   -4.867   -1.875)    8.500
   3.281   ( -10.282   10.993    0.119)   15.052
   3.281   (  10.282   10.993    0.119)   15.052
   6.028   (  -0.691   -2.693    0.072)    2.781
   6.028   (   0.691   -2.693    0.072)    2.781
   7.809   ( -15.396   16.018   -0.029)   22.218
   7.809   (  15.396   16.018   -0.029)   22.218
  10.077   (  -1.288   -0.456    1.366)    1.932
  10.077   (   1.288   -0.456    1.366)    1.932
======================= Grid point 266 (54/63) =======================
q-point: ( 0.20  0.20  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 141
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.562   (   4.926    4.926    6.464)    9.504
   0.829   (   3.399    3.399   -5.513)    7.314
   1.115   (   3.557    3.557   14.413)   15.265
   1.343   (   4.181    4.181   -2.984)    6.624
   1.391   (   0.080    0.080    0.529)    0.541
   1.461   (   4.240    4.240   -1.152)    6.107
   1.527   (   2.602    2.602    0.021)    3.680
   1.533   (   2.731    2.731   -4.898)    6.237
   1.662   (   2.882    2.882    1.213)    4.252
   1.991   (  -0.178   -0.178   -0.983)    1.014
   2.157   (   0.025    0.025   -0.329)    0.331
   2.219   (   0.763    0.763    0.263)    1.111
   2.255   (  -1.666   -1.666   -0.853)    2.505
   2.388   (   3.740    3.740    0.754)    5.342
   2.755   (   2.328    2.328    0.238)    3.301
   2.798   (   4.930    4.930   -2.167)    7.301
   3.075   (  -9.577   -9.577   -0.023)   13.543
   3.491   (   5.021    5.021    0.217)    7.103
   5.896   (  -0.627   -0.627    0.713)    1.138
   6.062   (  -1.168   -1.168    0.295)    1.678
   6.800   (  21.235   21.235   -0.261)   30.033
   7.955   (  19.733   19.733   -6.282)   28.605
  10.288   (  -6.818   -6.818    5.085)   10.901
  10.385   (  -4.398   -4.398   -0.012)    6.219
======================= Grid point 267 (55/63) =======================
q-point: ( 0.30  0.20  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.778   (  14.147   -5.670    6.029)   16.390
   0.957   (   9.050   -2.537   -5.792)   11.040
   1.186   (   3.756    2.134   13.201)   13.890
   1.404   (   2.380    2.207   -4.117)    5.242
   1.443   (  -1.078   -2.065   -2.699)    3.565
   1.471   (   5.342   -7.164   -1.343)    9.037
   1.572   (   2.521    1.471   -0.233)    2.928
   1.623   (   2.931   -0.338   -0.030)    2.950
   1.730   (   2.536    2.730    1.555)    4.037
   1.997   (   1.531    0.560    0.080)    1.632
   2.109   (  -3.365    3.145   -0.557)    4.640
   2.220   (  -0.145   -0.718    1.526)    1.692
   2.262   (   1.282   -3.850   -2.682)    4.864
   2.420   (   0.617   -0.757    1.171)    1.524
   2.689   (  -4.524   -3.822   -3.026)    6.651
   2.834   ( -14.592  -13.618    0.028)   19.960
   3.052   (  14.654   15.121    0.089)   21.057
   3.578   (   3.122    6.573    0.225)    7.280
   5.917   (   2.171   -1.047    0.480)    2.458
   6.035   (  -1.383   -1.025    0.188)    1.732
   7.281   (  26.783   23.027   -0.151)   35.322
   8.255   (  10.497   18.586   -3.498)   21.630
  10.164   (  -5.634   -3.950    4.152)    8.036
  10.279   (  -6.393   -2.039   -0.041)    6.711
======================= Grid point 268 (56/63) =======================
q-point: ( 0.40  0.20  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 252
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.035   (  11.295   -5.815    5.460)   13.827
   1.150   (   9.389   -4.514   -4.969)   11.542
   1.254   (   2.671   -0.459    9.545)    9.922
   1.395   (  -1.997   -1.456    0.047)    2.472
   1.450   (   1.603    0.694   -5.395)    5.671
   1.534   (  -1.525   -2.794    0.211)    3.190
   1.585   (  -0.640   -1.344   -1.057)    1.826
   1.713   (  -3.576   -2.251    1.479)    4.477
   1.778   (  13.997  -14.552    0.079)   20.191
   2.050   (   3.289    1.909    0.100)    3.805
   2.077   (   0.081    4.086   -0.191)    4.091
   2.219   (   0.147   -1.216    2.569)    2.846
   2.275   (  -0.357   -2.375   -1.427)    2.794
   2.443   (   2.036   -1.795   -0.947)    2.875
   2.515   ( -17.162  -16.766   -0.177)   23.993
   2.599   (  -4.838   -3.856   -3.197)    6.964
   3.318   (  11.990   13.812    0.087)   18.291
   3.595   (  -1.666    8.897    0.191)    9.054
   5.961   (   2.068   -0.584    0.200)    2.158
   6.012   (  -0.962   -0.669    0.078)    1.175
   7.800   (  24.825   23.388   -0.037)   34.107
   8.352   (  -1.029   19.987   -1.378)   20.061
  10.076   (  -3.016   -2.223    2.880)    4.726
  10.148   (  -6.302   -0.830    0.270)    6.362
======================= Grid point 269 (57/63) =======================
q-point: (-0.50  0.20  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.201   (  -4.982   -3.063    4.538)    7.403
   1.201   (   4.982   -3.063    4.538)    7.403
   1.336   (  -5.638   -3.602   -1.803)    6.929
   1.336   (   5.638   -3.602   -1.803)    6.929
   1.480   (  -2.271   -2.186   -0.012)    3.152
   1.480   (   2.271   -2.186   -0.012)    3.152
   1.619   (  -4.073   -1.536   -0.339)    4.366
   1.619   (   4.073   -1.536   -0.339)    4.366
   2.048   (  -6.431   -2.446    0.181)    6.882
   2.048   (   6.431   -2.446    0.181)    6.882
   2.159   ( -11.776  -11.484    0.220)   16.450
   2.159   (  11.776  -11.484    0.220)   16.450
   2.256   (  -0.787   -2.928    1.096)    3.224
   2.256   (   0.787   -2.928    1.096)    3.224
   2.510   (  -4.212   -2.786   -2.958)    5.853
   2.510   (   4.212   -2.786   -2.958)    5.853
   3.508   (  -7.211   11.420    0.136)   13.507
   3.508   (   7.211   11.420    0.136)   13.507
   5.993   (  -1.178   -0.384    0.063)    1.240
   5.993   (   1.178   -0.384    0.063)    1.240
   8.198   ( -14.706   22.136   -0.254)   26.577
   8.198   (  14.706   22.136   -0.254)   26.577
  10.062   (  -1.945   -1.121    1.338)    2.614
  10.062   (   1.945   -1.121    1.338)    2.614
======================= Grid point 278 (58/63) =======================
q-point: ( 0.30  0.30  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 140
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.766   (   5.050    5.050    6.322)    9.538
   0.959   (   2.920    2.920   -5.971)    7.260
   1.228   (   1.870    1.870   11.092)   11.402
   1.364   (  -1.479   -1.479    0.550)    2.163
   1.402   (  -2.213   -2.213   -0.974)    3.277
   1.443   (   1.474    1.474   -5.760)    6.125
   1.571   (  -4.925   -4.925    0.313)    6.973
   1.614   (   1.697    1.697    0.265)    2.415
   1.746   (  -1.755   -1.755    1.099)    2.714
   2.039   (   3.251    3.251    0.039)    4.597
   2.153   (  -0.365   -0.365   -0.457)    0.689
   2.175   (  -0.113   -0.113   -0.274)    0.317
   2.244   (   0.527    0.527    0.254)    0.787
   2.404   (  -0.521   -0.521   -0.520)    0.903
   2.518   ( -18.031  -18.031    0.087)   25.500
   2.607   (  -4.334   -4.334   -3.139)    6.887
   3.335   (  13.339   13.339    0.035)   18.864
   3.698   (   4.951    4.951    0.251)    7.006
   5.935   (   2.347    2.347    0.244)    3.329
   6.019   (  -0.447   -0.447    0.113)    0.642
   7.779   (  26.725   26.725   -0.098)   37.794
   8.563   (  12.230   12.230   -2.240)   17.441
  10.086   (  -3.774   -3.774    3.221)    6.234
  10.218   (  -4.274   -4.274    0.277)    6.051
======================= Grid point 279 (59/63) =======================
q-point: ( 0.40  0.30  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 250
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.937   (  10.472   -3.432    6.264)   12.676
   1.071   (   8.122   -3.124   -6.010)   10.576
   1.236   (  -2.013   -1.968    7.428)    7.944
   1.358   (  -3.066   -2.988   -2.660)    5.040
   1.396   (   1.470   -5.867    2.044)    6.384
   1.463   (   1.430   -2.309   -1.754)    3.233
   1.491   (   1.052   -8.952    1.092)    9.080
   1.633   (   0.224    1.661    0.118)    1.680
   1.675   (  -4.041   -2.447   -0.059)    4.724
   2.081   ( -12.840  -12.529   -0.034)   17.940
   2.102   (   2.651    2.553    0.137)    3.683
   2.181   (  -4.011   -2.282    0.255)    4.622
   2.219   (  -1.771   -3.483    0.522)    3.942
   2.255   (   0.180    0.042    0.236)    0.300
   2.409   (   1.191   -1.661   -2.108)    2.936
   2.524   (  -4.060   -3.873   -3.137)    6.428
   3.562   (   9.666   10.937    0.061)   14.596
   3.752   (   0.358    6.646    0.211)    6.659
   5.985   (   2.230    2.485   -0.046)    3.339
   6.017   (   0.067    1.058   -0.044)    1.061
   8.278   (  22.905   24.370   -0.099)   33.445
   8.711   (   2.213   15.810   -1.080)   16.001
  10.027   (  -1.842   -2.605    2.257)    3.908
  10.117   (  -5.572   -2.549    0.631)    6.160
======================= Grid point 280 (60/63) =======================
q-point: (-0.50  0.30  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.124   (  -7.990   -4.815    5.282)   10.720
   1.124   (   7.990   -4.815    5.282)   10.720
   1.242   (  -7.544   -5.290   -4.573)   10.286
   1.242   (   7.544   -5.290   -4.573)   10.286
   1.430   (  -2.882   -1.949    3.087)    4.651
   1.430   (   2.882   -1.949    3.087)    4.651
   1.625   (  -1.134    0.682   -0.552)    1.433
   1.625   (   1.134    0.682   -0.552)    1.433
   1.734   ( -16.636  -19.116   -0.023)   25.341
   1.734   (  16.636  -19.116   -0.023)   25.341
   2.144   (  -1.556    1.832   -0.045)    2.404
   2.144   (   1.556    1.832   -0.045)    2.404
   2.243   (  -1.357    0.395    0.866)    1.657
   2.243   (   1.357    0.395    0.866)    1.657
   2.453   (  -3.128   -2.965   -2.993)    5.247
   2.453   (   3.128   -2.965   -2.993)    5.247
   3.707   (  -4.959    8.591    0.135)    9.921
   3.707   (   4.959    8.591    0.135)    9.921
   6.013   (  -0.678    2.159   -0.140)    2.267
   6.013   (   0.678    2.159   -0.140)    2.267
   8.624   ( -11.415   20.350   -0.337)   23.335
   8.624   (  11.415   20.350   -0.337)   23.335
  10.031   (  -2.512   -2.103    1.291)    3.521
  10.031   (   2.512   -2.103    1.291)    3.521
======================= Grid point 290 (61/63) =======================
q-point: ( 0.40  0.40  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 141
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.931   (   3.097    3.097    6.120)    7.526
   1.055   (   2.254    2.254   -5.715)    6.544
   1.160   (  -5.836   -5.836    5.732)   10.049
   1.277   (  -5.732   -5.732   -5.012)    9.531
   1.286   (  -2.357   -2.357    0.292)    3.346
   1.393   (  -0.077   -0.077    4.370)    4.372
   1.445   (  -1.983   -1.983    1.876)    3.374
   1.647   (  -0.166   -0.166   -0.580)    0.626
   1.657   (   0.549    0.549    0.359)    0.856
   1.706   ( -19.882  -19.882    0.096)   28.117
   2.147   (   1.802    1.802    0.152)    2.553
   2.177   (  -0.002   -0.002   -0.334)    0.334
   2.233   (   1.685    1.685    0.600)    2.458
   2.262   (   0.346    0.346    0.152)    0.513
   2.387   (  -0.388   -0.388   -2.773)    2.827
   2.452   (  -3.329   -3.329   -3.001)    5.583
   3.740   (   7.069    7.069    0.074)    9.997
   3.841   (   2.186    2.186    0.179)    3.097
   6.036   (   2.229    2.229   -0.313)    3.168
   6.044   (   1.416    1.416   -0.246)    2.018
   8.710   (  18.519   18.519   -0.163)   26.190
   8.942   (   7.104    7.104   -0.643)   10.067
   9.983   (  -1.568   -1.568    1.717)    2.804
  10.049   (  -4.070   -4.070    0.886)    5.823
======================= Grid point 291 (62/63) =======================
q-point: (-0.50  0.40  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 150
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.032   (  -6.373   -3.794    5.014)    8.953
   1.032   (   6.373   -3.794    5.014)    8.953
   1.156   (  -6.410   -3.091   -4.514)    8.427
   1.156   (   6.410   -3.091   -4.514)    8.427
   1.373   ( -10.988  -12.856    0.333)   16.915
   1.373   (  10.988  -12.856    0.333)   16.915
   1.412   (  -2.281   -2.450    3.162)    4.605
   1.412   (   2.281   -2.450    3.162)    4.605
   1.656   (  -0.597    1.021   -0.062)    1.185
   1.656   (   0.597    1.021   -0.062)    1.185
   2.171   (  -0.707    0.966   -0.116)    1.203
   2.171   (   0.707    0.966   -0.116)    1.203
   2.258   (  -0.791    0.911    0.308)    1.246
   2.258   (   0.791    0.911    0.308)    1.246
   2.400   (  -1.807   -2.158   -2.992)    4.108
   2.400   (   1.807   -2.158   -2.992)    4.108
   3.836   (  -2.656    4.443    0.122)    5.178
   3.836   (   2.656    4.443    0.122)    5.178
   6.060   (  -0.112    2.117   -0.381)    2.154
   6.060   (   0.112    2.117   -0.381)    2.154
   8.958   (  -6.054   12.857   -0.308)   14.214
   8.958   (   6.054   12.857   -0.308)   14.214
   9.983   (  -1.935   -2.464    1.230)    3.366
   9.983   (   1.935   -2.464    1.230)    3.366
======================= Grid point 307 (63/63) =======================
q-point: (-0.50 -0.50  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.56e-04 1.56e-04 1.56e-04 1.56e-04 0.00e+00 0.00e+00 
Number of triplets: 54
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.995   (   0.000    0.000    4.606)    4.606
   0.995   (   0.000    0.000    4.606)    4.606
   1.133   (   0.000    0.000   -4.046)    4.046
   1.133   (   0.000    0.000   -4.046)    4.046
   1.210   (  -0.000   -0.000    0.041)    0.041
   1.210   (   0.000    0.000    0.041)    0.041
   1.403   (   0.000    0.000    3.324)    3.324
   1.403   (   0.000    0.000    3.324)    3.324
   1.667   (   0.000    0.000    0.315)    0.315
   1.667   (   0.000    0.000    0.315)    0.315
   2.181   (   0.000    0.000   -0.195)    0.195
   2.181   (   0.000    0.000   -0.195)    0.195
   2.269   (   0.000    0.000    0.099)    0.099
   2.269   (   0.000    0.000    0.099)    0.099
   2.377   (   0.000    0.000   -2.984)    2.984
   2.377   (   0.000    0.000   -2.984)    2.984
   3.879   (   0.000    0.000    0.116)    0.116
   3.879   (   0.000    0.000    0.116)    0.116
   6.082   (   0.000    0.000   -0.484)    0.484
   6.082   (   0.000    0.000   -0.484)    0.484
   9.091   (   0.000    0.000   -0.279)    0.279
   9.091   (   0.000    0.000   -0.279)    0.279
   9.954   (   0.000    0.000    1.201)    1.201
   9.954   (   0.000    0.000    1.201)    1.201
=================== 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/12000
   10.0     23.616     23.616     12.551      0.000      0.000      0.000 3/12000
   20.0     10.164     10.164      5.747      0.000      0.000      0.000 3/12000
   30.0      6.962      6.962      3.524      0.000      0.000      0.000 3/12000
   40.0      5.424      5.424      2.537      0.000      0.000      0.000 3/12000
   50.0      4.429      4.429      1.986      0.000      0.000      0.000 3/12000
   60.0      3.714      3.714      1.633      0.000      0.000      0.000 3/12000
   70.0      3.181      3.181      1.388      0.000      0.000      0.000 3/12000
   80.0      2.774      2.774      1.207      0.000      0.000      0.000 3/12000
   90.0      2.456      2.456      1.068      0.000      0.000      0.000 3/12000
  100.0      2.202      2.202      0.958      0.000      0.000      0.000 3/12000
  110.0      1.996      1.996      0.869      0.000      0.000      0.000 3/12000
  120.0      1.825      1.825      0.796      0.000      0.000      0.000 3/12000
  130.0      1.681      1.681      0.734      0.000      0.000      0.000 3/12000
  140.0      1.559      1.559      0.681      0.000      0.000      0.000 3/12000
  150.0      1.453      1.453      0.635      0.000      0.000      0.000 3/12000
  160.0      1.361      1.361      0.595      0.000      0.000      0.000 3/12000
  170.0      1.280      1.280      0.561      0.000      0.000      0.000 3/12000
  180.0      1.209      1.209      0.530      0.000      0.000      0.000 3/12000
  190.0      1.145      1.145      0.502      0.000      0.000      0.000 3/12000
  200.0      1.087      1.087      0.477      0.000      0.000      0.000 3/12000
  210.0      1.036      1.036      0.454      0.000      0.000      0.000 3/12000
  220.0      0.989      0.989      0.434      0.000      0.000      0.000 3/12000
  230.0      0.946      0.946      0.415      0.000      0.000      0.000 3/12000
  240.0      0.906      0.906      0.398      0.000      0.000      0.000 3/12000
  250.0      0.870      0.870      0.383      0.000      0.000      0.000 3/12000
  260.0      0.837      0.837      0.368      0.000      0.000      0.000 3/12000
  270.0      0.806      0.806      0.355      0.000      0.000      0.000 3/12000
  280.0      0.777      0.777      0.342      0.000      0.000      0.000 3/12000
  290.0      0.751      0.751      0.331      0.000      0.000      0.000 3/12000
  300.0      0.726      0.726      0.320      0.000      0.000      0.000 3/12000
  310.0      0.703      0.703      0.310      0.000      0.000      0.000 3/12000
  320.0      0.681      0.681      0.300      0.000      0.000      0.000 3/12000
  330.0      0.660      0.660      0.291      0.000      0.000      0.000 3/12000
  340.0      0.641      0.641      0.283      0.000      0.000      0.000 3/12000
  350.0      0.623      0.623      0.275      0.000      0.000      0.000 3/12000
  360.0      0.606      0.606      0.267      0.000      0.000      0.000 3/12000
  370.0      0.590      0.590      0.260      0.000      0.000      0.000 3/12000
  380.0      0.574      0.574      0.254      0.000      0.000      0.000 3/12000
  390.0      0.560      0.560      0.247      0.000      0.000      0.000 3/12000
  400.0      0.546      0.546      0.241      0.000      0.000      0.000 3/12000
  410.0      0.533      0.533      0.235      0.000      0.000      0.000 3/12000
  420.0      0.520      0.520      0.230      0.000      0.000      0.000 3/12000
  430.0      0.508      0.508      0.225      0.000      0.000      0.000 3/12000
  440.0      0.497      0.497      0.220      0.000      0.000      0.000 3/12000
  450.0      0.486      0.486      0.215      0.000      0.000      0.000 3/12000
  460.0      0.475      0.475      0.210      0.000      0.000      0.000 3/12000
  470.0      0.465      0.465      0.206      0.000      0.000      0.000 3/12000
  480.0      0.456      0.456      0.202      0.000      0.000      0.000 3/12000
  490.0      0.447      0.447      0.198      0.000      0.000      0.000 3/12000
  500.0      0.438      0.438      0.194      0.000      0.000      0.000 3/12000
  510.0      0.429      0.429      0.190      0.000      0.000      0.000 3/12000
  520.0      0.421      0.421      0.186      0.000      0.000      0.000 3/12000
  530.0      0.413      0.413      0.183      0.000      0.000      0.000 3/12000
  540.0      0.406      0.406      0.180      0.000      0.000      0.000 3/12000
  550.0      0.398      0.398      0.176      0.000      0.000      0.000 3/12000
  560.0      0.391      0.391      0.173      0.000      0.000      0.000 3/12000
  570.0      0.384      0.384      0.170      0.000      0.000      0.000 3/12000
  580.0      0.378      0.378      0.168      0.000      0.000      0.000 3/12000
  590.0      0.372      0.372      0.165      0.000      0.000      0.000 3/12000
  600.0      0.365      0.365      0.162      0.000      0.000      0.000 3/12000
  610.0      0.360      0.360      0.159      0.000      0.000      0.000 3/12000
  620.0      0.354      0.354      0.157      0.000      0.000      0.000 3/12000
  630.0      0.348      0.348      0.154      0.000      0.000      0.000 3/12000
  640.0      0.343      0.343      0.152      0.000      0.000      0.000 3/12000
  650.0      0.338      0.338      0.150      0.000      0.000      0.000 3/12000
  660.0      0.333      0.333      0.148      0.000      0.000      0.000 3/12000
  670.0      0.328      0.328      0.145      0.000      0.000      0.000 3/12000
  680.0      0.323      0.323      0.143      0.000      0.000      0.000 3/12000
  690.0      0.318      0.318      0.141      0.000      0.000      0.000 3/12000
  700.0      0.314      0.314      0.139      0.000      0.000      0.000 3/12000
  710.0      0.309      0.309      0.137      0.000      0.000      0.000 3/12000
  720.0      0.305      0.305      0.136      0.000      0.000      0.000 3/12000
  730.0      0.301      0.301      0.134      0.000      0.000      0.000 3/12000
  740.0      0.297      0.297      0.132      0.000      0.000      0.000 3/12000
  750.0      0.293      0.293      0.130      0.000      0.000      0.000 3/12000
  760.0      0.289      0.289      0.129      0.000      0.000      0.000 3/12000
  770.0      0.286      0.286      0.127      0.000      0.000      0.000 3/12000
  780.0      0.282      0.282      0.125      0.000      0.000      0.000 3/12000
  790.0      0.278      0.278      0.124      0.000      0.000      0.000 3/12000
  800.0      0.275      0.275      0.122      0.000      0.000      0.000 3/12000
  810.0      0.272      0.272      0.121      0.000      0.000      0.000 3/12000
  820.0      0.268      0.268      0.119      0.000      0.000      0.000 3/12000
  830.0      0.265      0.265      0.118      0.000      0.000      0.000 3/12000
  840.0      0.262      0.262      0.117      0.000      0.000      0.000 3/12000
  850.0      0.259      0.259      0.115      0.000      0.000      0.000 3/12000
  860.0      0.256      0.256      0.114      0.000      0.000      0.000 3/12000
  870.0      0.253      0.253      0.113      0.000      0.000      0.000 3/12000
  880.0      0.250      0.250      0.111      0.000      0.000      0.000 3/12000
  890.0      0.247      0.247      0.110      0.000      0.000      0.000 3/12000
  900.0      0.245      0.245      0.109      0.000      0.000      0.000 3/12000
  910.0      0.242      0.242      0.108      0.000      0.000      0.000 3/12000
  920.0      0.239      0.239      0.107      0.000      0.000      0.000 3/12000
  930.0      0.237      0.237      0.105      0.000      0.000      0.000 3/12000
  940.0      0.234      0.234      0.104      0.000      0.000      0.000 3/12000
  950.0      0.232      0.232      0.103      0.000      0.000      0.000 3/12000
  960.0      0.230      0.230      0.102      0.000      0.000      0.000 3/12000
  970.0      0.227      0.227      0.101      0.000      0.000      0.000 3/12000
  980.0      0.225      0.225      0.100      0.000      0.000      0.000 3/12000
  990.0      0.223      0.223      0.099      0.000      0.000      0.000 3/12000
 1000.0      0.220      0.220      0.098      0.000      0.000      0.000 3/12000

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

