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

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

-------------------------[time 2026-01-08 05:23:48]-------------------------
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    4.842505820000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.842505820000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.364425330000000
Atomic positions (fractional):
   *1 Li  0.50000000000000  0.00000000000000  0.08114733798113   6.941
    2 Li  0.00000000000000  0.50000000000000  0.91885266201887   6.941
   *3 Cl  0.00000000000000  0.50000000000000  0.66980280217486  35.453
    4 Cl  0.50000000000000  0.00000000000000  0.33019719782514  35.453
   *5 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453
    6 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453
   *7 Cs  0.00000000000000  0.50000000000000  0.30233586163350 132.905
    8 Cs  0.50000000000000  0.00000000000000  0.69766413836650 132.905
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.842505820000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.842505820000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.364425330000000
Atomic positions (fractional):
   *1 Li  0.50000000000000  0.00000000000000  0.08114733798113   6.941 > 1
    2 Li  0.00000000000000  0.50000000000000  0.91885266201887   6.941 > 2
   *3 Cl  0.00000000000000  0.50000000000000  0.66980280217486  35.453 > 3
    4 Cl  0.50000000000000  0.00000000000000  0.33019719782514  35.453 > 4
   *5 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 5
    6 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 6
   *7 Cs  0.00000000000000  0.50000000000000  0.30233586163350 132.905 > 7
    8 Cs  0.50000000000000  0.00000000000000  0.69766413836650 132.905 > 8
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   14.527517460000002    0.000000000000000    0.000000000000000
  b    0.000000000000000   14.527517460000002    0.000000000000000
  c    0.000000000000000    0.000000000000000   18.728850659999999
Atomic positions (fractional):
   *1 Li  0.16666666666667  0.00000000000000  0.04057366899056   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.04057366899056   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.04057366899056   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.04057366899056   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.04057366899056   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.04057366899056   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.04057366899056   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.04057366899056   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.04057366899056   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.54057366899056   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.54057366899056   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.54057366899056   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.54057366899056   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.54057366899056   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.54057366899056   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.54057366899056   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.54057366899056   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.54057366899056   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.45942633100944   6.941 > 2
   20 Li  0.33333333333333  0.16666666666667  0.45942633100944   6.941 > 2
   21 Li  0.66666666666667  0.16666666666667  0.45942633100944   6.941 > 2
   22 Li  0.00000000000000  0.50000000000000  0.45942633100944   6.941 > 2
   23 Li  0.33333333333333  0.50000000000000  0.45942633100944   6.941 > 2
   24 Li  0.66666666666667  0.50000000000000  0.45942633100944   6.941 > 2
   25 Li  0.00000000000000  0.83333333333333  0.45942633100944   6.941 > 2
   26 Li  0.33333333333333  0.83333333333333  0.45942633100944   6.941 > 2
   27 Li  0.66666666666667  0.83333333333333  0.45942633100944   6.941 > 2
   28 Li  0.00000000000000  0.16666666666667  0.95942633100944   6.941 > 2
   29 Li  0.33333333333333  0.16666666666667  0.95942633100944   6.941 > 2
   30 Li  0.66666666666667  0.16666666666667  0.95942633100944   6.941 > 2
   31 Li  0.00000000000000  0.50000000000000  0.95942633100944   6.941 > 2
   32 Li  0.33333333333333  0.50000000000000  0.95942633100944   6.941 > 2
   33 Li  0.66666666666667  0.50000000000000  0.95942633100944   6.941 > 2
   34 Li  0.00000000000000  0.83333333333333  0.95942633100944   6.941 > 2
   35 Li  0.33333333333333  0.83333333333333  0.95942633100944   6.941 > 2
   36 Li  0.66666666666667  0.83333333333333  0.95942633100944   6.941 > 2
  *37 Cl  0.00000000000000  0.16666666666667  0.33490140108743  35.453 > 3
   38 Cl  0.33333333333333  0.16666666666667  0.33490140108743  35.453 > 3
   39 Cl  0.66666666666667  0.16666666666667  0.33490140108743  35.453 > 3
   40 Cl  0.00000000000000  0.50000000000000  0.33490140108743  35.453 > 3
   41 Cl  0.33333333333333  0.50000000000000  0.33490140108743  35.453 > 3
   42 Cl  0.66666666666667  0.50000000000000  0.33490140108743  35.453 > 3
   43 Cl  0.00000000000000  0.83333333333333  0.33490140108743  35.453 > 3
   44 Cl  0.33333333333333  0.83333333333333  0.33490140108743  35.453 > 3
   45 Cl  0.66666666666667  0.83333333333333  0.33490140108743  35.453 > 3
   46 Cl  0.00000000000000  0.16666666666667  0.83490140108743  35.453 > 3
   47 Cl  0.33333333333333  0.16666666666667  0.83490140108743  35.453 > 3
   48 Cl  0.66666666666667  0.16666666666667  0.83490140108743  35.453 > 3
   49 Cl  0.00000000000000  0.50000000000000  0.83490140108743  35.453 > 3
   50 Cl  0.33333333333333  0.50000000000000  0.83490140108743  35.453 > 3
   51 Cl  0.66666666666667  0.50000000000000  0.83490140108743  35.453 > 3
   52 Cl  0.00000000000000  0.83333333333333  0.83490140108743  35.453 > 3
   53 Cl  0.33333333333333  0.83333333333333  0.83490140108743  35.453 > 3
   54 Cl  0.66666666666667  0.83333333333333  0.83490140108743  35.453 > 3
   55 Cl  0.16666666666667  0.00000000000000  0.16509859891257  35.453 > 4
   56 Cl  0.50000000000000  0.00000000000000  0.16509859891257  35.453 > 4
   57 Cl  0.83333333333333  0.00000000000000  0.16509859891257  35.453 > 4
   58 Cl  0.16666666666667  0.33333333333333  0.16509859891257  35.453 > 4
   59 Cl  0.50000000000000  0.33333333333333  0.16509859891257  35.453 > 4
   60 Cl  0.83333333333333  0.33333333333333  0.16509859891257  35.453 > 4
   61 Cl  0.16666666666667  0.66666666666667  0.16509859891257  35.453 > 4
   62 Cl  0.50000000000000  0.66666666666667  0.16509859891257  35.453 > 4
   63 Cl  0.83333333333333  0.66666666666667  0.16509859891257  35.453 > 4
   64 Cl  0.16666666666667  0.00000000000000  0.66509859891257  35.453 > 4
   65 Cl  0.50000000000000  0.00000000000000  0.66509859891257  35.453 > 4
   66 Cl  0.83333333333333  0.00000000000000  0.66509859891257  35.453 > 4
   67 Cl  0.16666666666667  0.33333333333333  0.66509859891257  35.453 > 4
   68 Cl  0.50000000000000  0.33333333333333  0.66509859891257  35.453 > 4
   69 Cl  0.83333333333333  0.33333333333333  0.66509859891257  35.453 > 4
   70 Cl  0.16666666666667  0.66666666666667  0.66509859891257  35.453 > 4
   71 Cl  0.50000000000000  0.66666666666667  0.66509859891257  35.453 > 4
   72 Cl  0.83333333333333  0.66666666666667  0.66509859891257  35.453 > 4
  *73 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 5
   74 Cl  0.33333333333333  0.00000000000000  0.00000000000000  35.453 > 5
   75 Cl  0.66666666666667  0.00000000000000  0.00000000000000  35.453 > 5
   76 Cl  0.00000000000000  0.33333333333333  0.00000000000000  35.453 > 5
   77 Cl  0.33333333333333  0.33333333333333  0.00000000000000  35.453 > 5
   78 Cl  0.66666666666667  0.33333333333333  0.00000000000000  35.453 > 5
   79 Cl  0.00000000000000  0.66666666666667  0.00000000000000  35.453 > 5
   80 Cl  0.33333333333333  0.66666666666667  0.00000000000000  35.453 > 5
   81 Cl  0.66666666666667  0.66666666666667  0.00000000000000  35.453 > 5
   82 Cl  0.00000000000000  0.00000000000000  0.50000000000000  35.453 > 5
   83 Cl  0.33333333333333  0.00000000000000  0.50000000000000  35.453 > 5
   84 Cl  0.66666666666667  0.00000000000000  0.50000000000000  35.453 > 5
   85 Cl  0.00000000000000  0.33333333333333  0.50000000000000  35.453 > 5
   86 Cl  0.33333333333333  0.33333333333333  0.50000000000000  35.453 > 5
   87 Cl  0.66666666666667  0.33333333333333  0.50000000000000  35.453 > 5
   88 Cl  0.00000000000000  0.66666666666667  0.50000000000000  35.453 > 5
   89 Cl  0.33333333333333  0.66666666666667  0.50000000000000  35.453 > 5
   90 Cl  0.66666666666667  0.66666666666667  0.50000000000000  35.453 > 5
   91 Cl  0.16666666666667  0.16666666666667  0.00000000000000  35.453 > 6
   92 Cl  0.50000000000000  0.16666666666667  0.00000000000000  35.453 > 6
   93 Cl  0.83333333333333  0.16666666666667  0.00000000000000  35.453 > 6
   94 Cl  0.16666666666667  0.50000000000000  0.00000000000000  35.453 > 6
   95 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 6
   96 Cl  0.83333333333333  0.50000000000000  0.00000000000000  35.453 > 6
   97 Cl  0.16666666666667  0.83333333333333  0.00000000000000  35.453 > 6
   98 Cl  0.50000000000000  0.83333333333333  0.00000000000000  35.453 > 6
   99 Cl  0.83333333333333  0.83333333333333  0.00000000000000  35.453 > 6
  100 Cl  0.16666666666667  0.16666666666667  0.50000000000000  35.453 > 6
  101 Cl  0.50000000000000  0.16666666666667  0.50000000000000  35.453 > 6
  102 Cl  0.83333333333333  0.16666666666667  0.50000000000000  35.453 > 6
  103 Cl  0.16666666666667  0.50000000000000  0.50000000000000  35.453 > 6
  104 Cl  0.50000000000000  0.50000000000000  0.50000000000000  35.453 > 6
  105 Cl  0.83333333333333  0.50000000000000  0.50000000000000  35.453 > 6
  106 Cl  0.16666666666667  0.83333333333333  0.50000000000000  35.453 > 6
  107 Cl  0.50000000000000  0.83333333333333  0.50000000000000  35.453 > 6
  108 Cl  0.83333333333333  0.83333333333333  0.50000000000000  35.453 > 6
 *109 Cs  0.00000000000000  0.16666666666667  0.15116793081675 132.905 > 7
  110 Cs  0.33333333333333  0.16666666666667  0.15116793081675 132.905 > 7
  111 Cs  0.66666666666667  0.16666666666667  0.15116793081675 132.905 > 7
  112 Cs  0.00000000000000  0.50000000000000  0.15116793081675 132.905 > 7
  113 Cs  0.33333333333333  0.50000000000000  0.15116793081675 132.905 > 7
  114 Cs  0.66666666666667  0.50000000000000  0.15116793081675 132.905 > 7
  115 Cs  0.00000000000000  0.83333333333333  0.15116793081675 132.905 > 7
  116 Cs  0.33333333333333  0.83333333333333  0.15116793081675 132.905 > 7
  117 Cs  0.66666666666667  0.83333333333333  0.15116793081675 132.905 > 7
  118 Cs  0.00000000000000  0.16666666666667  0.65116793081675 132.905 > 7
  119 Cs  0.33333333333333  0.16666666666667  0.65116793081675 132.905 > 7
  120 Cs  0.66666666666667  0.16666666666667  0.65116793081675 132.905 > 7
  121 Cs  0.00000000000000  0.50000000000000  0.65116793081675 132.905 > 7
  122 Cs  0.33333333333333  0.50000000000000  0.65116793081675 132.905 > 7
  123 Cs  0.66666666666667  0.50000000000000  0.65116793081675 132.905 > 7
  124 Cs  0.00000000000000  0.83333333333333  0.65116793081675 132.905 > 7
  125 Cs  0.33333333333333  0.83333333333333  0.65116793081675 132.905 > 7
  126 Cs  0.66666666666667  0.83333333333333  0.65116793081675 132.905 > 7
  127 Cs  0.16666666666667  0.00000000000000  0.34883206918325 132.905 > 8
  128 Cs  0.50000000000000  0.00000000000000  0.34883206918325 132.905 > 8
  129 Cs  0.83333333333333  0.00000000000000  0.34883206918325 132.905 > 8
  130 Cs  0.16666666666667  0.33333333333333  0.34883206918325 132.905 > 8
  131 Cs  0.50000000000000  0.33333333333333  0.34883206918325 132.905 > 8
  132 Cs  0.83333333333333  0.33333333333333  0.34883206918325 132.905 > 8
  133 Cs  0.16666666666667  0.66666666666667  0.34883206918325 132.905 > 8
  134 Cs  0.50000000000000  0.66666666666667  0.34883206918325 132.905 > 8
  135 Cs  0.83333333333333  0.66666666666667  0.34883206918325 132.905 > 8
  136 Cs  0.16666666666667  0.00000000000000  0.84883206918325 132.905 > 8
  137 Cs  0.50000000000000  0.00000000000000  0.84883206918325 132.905 > 8
  138 Cs  0.83333333333333  0.00000000000000  0.84883206918325 132.905 > 8
  139 Cs  0.16666666666667  0.33333333333333  0.84883206918325 132.905 > 8
  140 Cs  0.50000000000000  0.33333333333333  0.84883206918325 132.905 > 8
  141 Cs  0.83333333333333  0.33333333333333  0.84883206918325 132.905 > 8
  142 Cs  0.16666666666667  0.66666666666667  0.84883206918325 132.905 > 8
  143 Cs  0.50000000000000  0.66666666666667  0.84883206918325 132.905 > 8
  144 Cs  0.83333333333333  0.66666666666667  0.84883206918325 132.905 > 8
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.7762129    0.0000000    0.0000000
            0.0000000    2.7762129    0.0000000
            0.0000000    0.0000000    2.7028699
-------------------------- Born effective charges --------------------------
    1 Li    1.1526049    0.0000000    0.0000000
            0.0000000    1.1526049    0.0000000
            0.0000000    0.0000000    0.9775342
    2 Li    1.1526049    0.0000000    0.0000000
            0.0000000    1.1526049    0.0000000
            0.0000000    0.0000000    0.9775342
    3 Cl   -1.1747772    0.0000000    0.0000000
            0.0000000   -1.1747772    0.0000000
            0.0000000    0.0000000   -1.2902662
    4 Cl   -1.1747772    0.0000000    0.0000000
            0.0000000   -1.1747772    0.0000000
            0.0000000    0.0000000   -1.2902662
    5 Cl   -1.3011168    0.0000000    0.0000000
            0.0000000   -1.3011168    0.0000000
            0.0000000    0.0000000   -0.9574103
    6 Cl   -1.3011168    0.0000000    0.0000000
            0.0000000   -1.3011168    0.0000000
            0.0000000    0.0000000   -0.9574103
    7 Cs    1.3232891    0.0000000    0.0000000
            0.0000000    1.3232891    0.0000000
            0.0000000    0.0000000    1.2701423
    8 Cs    1.3232891    0.0000000    0.0000000
            0.0000000    1.3232891    0.0000000
            0.0000000    0.0000000    1.2701423
----------------------------------------------------------------------------
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: 444
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.370
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.390
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
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: 444
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 (yy) -0.00000000 (yy) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 05:23: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    4.842505820000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.842505820000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.364425330000000
Atomic positions (fractional):
    1 Li  0.50000000000000  0.00000000000000  0.08114733798113   6.941
    2 Li  0.00000000000000  0.50000000000000  0.91885266201887   6.941
    3 Cl  0.00000000000000  0.50000000000000  0.66980280217486  35.453
    4 Cl  0.50000000000000  0.00000000000000  0.33019719782514  35.453
    5 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453
    6 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453
    7 Cs  0.00000000000000  0.50000000000000  0.30233586163350 132.905
    8 Cs  0.50000000000000  0.00000000000000  0.69766413836650 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   14.527517460000002    0.000000000000000    0.000000000000000
  b    0.000000000000000   14.527517460000002    0.000000000000000
  c    0.000000000000000    0.000000000000000   18.728850659999999
Atomic positions (fractional):
    1 Li  0.16666666666667  0.00000000000000  0.04057366899056   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.04057366899056   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.04057366899056   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.04057366899056   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.04057366899056   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.04057366899056   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.04057366899056   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.04057366899056   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.04057366899056   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.54057366899056   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.54057366899056   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.54057366899056   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.54057366899056   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.54057366899056   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.54057366899056   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.54057366899056   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.54057366899056   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.54057366899056   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.45942633100944   6.941 > 19
   20 Li  0.33333333333333  0.16666666666667  0.45942633100944   6.941 > 19
   21 Li  0.66666666666667  0.16666666666667  0.45942633100944   6.941 > 19
   22 Li  0.00000000000000  0.50000000000000  0.45942633100944   6.941 > 19
   23 Li  0.33333333333333  0.50000000000000  0.45942633100944   6.941 > 19
   24 Li  0.66666666666667  0.50000000000000  0.45942633100944   6.941 > 19
   25 Li  0.00000000000000  0.83333333333333  0.45942633100944   6.941 > 19
   26 Li  0.33333333333333  0.83333333333333  0.45942633100944   6.941 > 19
   27 Li  0.66666666666667  0.83333333333333  0.45942633100944   6.941 > 19
   28 Li  0.00000000000000  0.16666666666667  0.95942633100944   6.941 > 19
   29 Li  0.33333333333333  0.16666666666667  0.95942633100944   6.941 > 19
   30 Li  0.66666666666667  0.16666666666667  0.95942633100944   6.941 > 19
   31 Li  0.00000000000000  0.50000000000000  0.95942633100944   6.941 > 19
   32 Li  0.33333333333333  0.50000000000000  0.95942633100944   6.941 > 19
   33 Li  0.66666666666667  0.50000000000000  0.95942633100944   6.941 > 19
   34 Li  0.00000000000000  0.83333333333333  0.95942633100944   6.941 > 19
   35 Li  0.33333333333333  0.83333333333333  0.95942633100944   6.941 > 19
   36 Li  0.66666666666667  0.83333333333333  0.95942633100944   6.941 > 19
   37 Cl  0.00000000000000  0.16666666666667  0.33490140108743  35.453 > 37
   38 Cl  0.33333333333333  0.16666666666667  0.33490140108743  35.453 > 37
   39 Cl  0.66666666666667  0.16666666666667  0.33490140108743  35.453 > 37
   40 Cl  0.00000000000000  0.50000000000000  0.33490140108743  35.453 > 37
   41 Cl  0.33333333333333  0.50000000000000  0.33490140108743  35.453 > 37
   42 Cl  0.66666666666667  0.50000000000000  0.33490140108743  35.453 > 37
   43 Cl  0.00000000000000  0.83333333333333  0.33490140108743  35.453 > 37
   44 Cl  0.33333333333333  0.83333333333333  0.33490140108743  35.453 > 37
   45 Cl  0.66666666666667  0.83333333333333  0.33490140108743  35.453 > 37
   46 Cl  0.00000000000000  0.16666666666667  0.83490140108743  35.453 > 37
   47 Cl  0.33333333333333  0.16666666666667  0.83490140108743  35.453 > 37
   48 Cl  0.66666666666667  0.16666666666667  0.83490140108743  35.453 > 37
   49 Cl  0.00000000000000  0.50000000000000  0.83490140108743  35.453 > 37
   50 Cl  0.33333333333333  0.50000000000000  0.83490140108743  35.453 > 37
   51 Cl  0.66666666666667  0.50000000000000  0.83490140108743  35.453 > 37
   52 Cl  0.00000000000000  0.83333333333333  0.83490140108743  35.453 > 37
   53 Cl  0.33333333333333  0.83333333333333  0.83490140108743  35.453 > 37
   54 Cl  0.66666666666667  0.83333333333333  0.83490140108743  35.453 > 37
   55 Cl  0.16666666666667  0.00000000000000  0.16509859891257  35.453 > 55
   56 Cl  0.50000000000000  0.00000000000000  0.16509859891257  35.453 > 55
   57 Cl  0.83333333333333  0.00000000000000  0.16509859891257  35.453 > 55
   58 Cl  0.16666666666667  0.33333333333333  0.16509859891257  35.453 > 55
   59 Cl  0.50000000000000  0.33333333333333  0.16509859891257  35.453 > 55
   60 Cl  0.83333333333333  0.33333333333333  0.16509859891257  35.453 > 55
   61 Cl  0.16666666666667  0.66666666666667  0.16509859891257  35.453 > 55
   62 Cl  0.50000000000000  0.66666666666667  0.16509859891257  35.453 > 55
   63 Cl  0.83333333333333  0.66666666666667  0.16509859891257  35.453 > 55
   64 Cl  0.16666666666667  0.00000000000000  0.66509859891257  35.453 > 55
   65 Cl  0.50000000000000  0.00000000000000  0.66509859891257  35.453 > 55
   66 Cl  0.83333333333333  0.00000000000000  0.66509859891257  35.453 > 55
   67 Cl  0.16666666666667  0.33333333333333  0.66509859891257  35.453 > 55
   68 Cl  0.50000000000000  0.33333333333333  0.66509859891257  35.453 > 55
   69 Cl  0.83333333333333  0.33333333333333  0.66509859891257  35.453 > 55
   70 Cl  0.16666666666667  0.66666666666667  0.66509859891257  35.453 > 55
   71 Cl  0.50000000000000  0.66666666666667  0.66509859891257  35.453 > 55
   72 Cl  0.83333333333333  0.66666666666667  0.66509859891257  35.453 > 55
   73 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 73
   74 Cl  0.33333333333333  0.00000000000000  0.00000000000000  35.453 > 73
   75 Cl  0.66666666666667  0.00000000000000  0.00000000000000  35.453 > 73
   76 Cl  0.00000000000000  0.33333333333333  0.00000000000000  35.453 > 73
   77 Cl  0.33333333333333  0.33333333333333  0.00000000000000  35.453 > 73
   78 Cl  0.66666666666667  0.33333333333333  0.00000000000000  35.453 > 73
   79 Cl  0.00000000000000  0.66666666666667  0.00000000000000  35.453 > 73
   80 Cl  0.33333333333333  0.66666666666667  0.00000000000000  35.453 > 73
   81 Cl  0.66666666666667  0.66666666666667  0.00000000000000  35.453 > 73
   82 Cl  0.00000000000000  0.00000000000000  0.50000000000000  35.453 > 73
   83 Cl  0.33333333333333  0.00000000000000  0.50000000000000  35.453 > 73
   84 Cl  0.66666666666667  0.00000000000000  0.50000000000000  35.453 > 73
   85 Cl  0.00000000000000  0.33333333333333  0.50000000000000  35.453 > 73
   86 Cl  0.33333333333333  0.33333333333333  0.50000000000000  35.453 > 73
   87 Cl  0.66666666666667  0.33333333333333  0.50000000000000  35.453 > 73
   88 Cl  0.00000000000000  0.66666666666667  0.50000000000000  35.453 > 73
   89 Cl  0.33333333333333  0.66666666666667  0.50000000000000  35.453 > 73
   90 Cl  0.66666666666667  0.66666666666667  0.50000000000000  35.453 > 73
   91 Cl  0.16666666666667  0.16666666666667  0.00000000000000  35.453 > 91
   92 Cl  0.50000000000000  0.16666666666667  0.00000000000000  35.453 > 91
   93 Cl  0.83333333333333  0.16666666666667  0.00000000000000  35.453 > 91
   94 Cl  0.16666666666667  0.50000000000000  0.00000000000000  35.453 > 91
   95 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 91
   96 Cl  0.83333333333333  0.50000000000000  0.00000000000000  35.453 > 91
   97 Cl  0.16666666666667  0.83333333333333  0.00000000000000  35.453 > 91
   98 Cl  0.50000000000000  0.83333333333333  0.00000000000000  35.453 > 91
   99 Cl  0.83333333333333  0.83333333333333  0.00000000000000  35.453 > 91
  100 Cl  0.16666666666667  0.16666666666667  0.50000000000000  35.453 > 91
  101 Cl  0.50000000000000  0.16666666666667  0.50000000000000  35.453 > 91
  102 Cl  0.83333333333333  0.16666666666667  0.50000000000000  35.453 > 91
  103 Cl  0.16666666666667  0.50000000000000  0.50000000000000  35.453 > 91
  104 Cl  0.50000000000000  0.50000000000000  0.50000000000000  35.453 > 91
  105 Cl  0.83333333333333  0.50000000000000  0.50000000000000  35.453 > 91
  106 Cl  0.16666666666667  0.83333333333333  0.50000000000000  35.453 > 91
  107 Cl  0.50000000000000  0.83333333333333  0.50000000000000  35.453 > 91
  108 Cl  0.83333333333333  0.83333333333333  0.50000000000000  35.453 > 91
  109 Cs  0.00000000000000  0.16666666666667  0.15116793081675 132.905 > 109
  110 Cs  0.33333333333333  0.16666666666667  0.15116793081675 132.905 > 109
  111 Cs  0.66666666666667  0.16666666666667  0.15116793081675 132.905 > 109
  112 Cs  0.00000000000000  0.50000000000000  0.15116793081675 132.905 > 109
  113 Cs  0.33333333333333  0.50000000000000  0.15116793081675 132.905 > 109
  114 Cs  0.66666666666667  0.50000000000000  0.15116793081675 132.905 > 109
  115 Cs  0.00000000000000  0.83333333333333  0.15116793081675 132.905 > 109
  116 Cs  0.33333333333333  0.83333333333333  0.15116793081675 132.905 > 109
  117 Cs  0.66666666666667  0.83333333333333  0.15116793081675 132.905 > 109
  118 Cs  0.00000000000000  0.16666666666667  0.65116793081675 132.905 > 109
  119 Cs  0.33333333333333  0.16666666666667  0.65116793081675 132.905 > 109
  120 Cs  0.66666666666667  0.16666666666667  0.65116793081675 132.905 > 109
  121 Cs  0.00000000000000  0.50000000000000  0.65116793081675 132.905 > 109
  122 Cs  0.33333333333333  0.50000000000000  0.65116793081675 132.905 > 109
  123 Cs  0.66666666666667  0.50000000000000  0.65116793081675 132.905 > 109
  124 Cs  0.00000000000000  0.83333333333333  0.65116793081675 132.905 > 109
  125 Cs  0.33333333333333  0.83333333333333  0.65116793081675 132.905 > 109
  126 Cs  0.66666666666667  0.83333333333333  0.65116793081675 132.905 > 109
  127 Cs  0.16666666666667  0.00000000000000  0.34883206918325 132.905 > 127
  128 Cs  0.50000000000000  0.00000000000000  0.34883206918325 132.905 > 127
  129 Cs  0.83333333333333  0.00000000000000  0.34883206918325 132.905 > 127
  130 Cs  0.16666666666667  0.33333333333333  0.34883206918325 132.905 > 127
  131 Cs  0.50000000000000  0.33333333333333  0.34883206918325 132.905 > 127
  132 Cs  0.83333333333333  0.33333333333333  0.34883206918325 132.905 > 127
  133 Cs  0.16666666666667  0.66666666666667  0.34883206918325 132.905 > 127
  134 Cs  0.50000000000000  0.66666666666667  0.34883206918325 132.905 > 127
  135 Cs  0.83333333333333  0.66666666666667  0.34883206918325 132.905 > 127
  136 Cs  0.16666666666667  0.00000000000000  0.84883206918325 132.905 > 127
  137 Cs  0.50000000000000  0.00000000000000  0.84883206918325 132.905 > 127
  138 Cs  0.83333333333333  0.00000000000000  0.84883206918325 132.905 > 127
  139 Cs  0.16666666666667  0.33333333333333  0.84883206918325 132.905 > 127
  140 Cs  0.50000000000000  0.33333333333333  0.84883206918325 132.905 > 127
  141 Cs  0.83333333333333  0.33333333333333  0.84883206918325 132.905 > 127
  142 Cs  0.16666666666667  0.66666666666667  0.84883206918325 132.905 > 127
  143 Cs  0.50000000000000  0.66666666666667  0.84883206918325 132.905 > 127
  144 Cs  0.83333333333333  0.66666666666667  0.84883206918325 132.905 > 127
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.7762129    0.0000000    0.0000000
            0.0000000    2.7762129    0.0000000
            0.0000000    0.0000000    2.7028699
-------------------------- Born effective charges --------------------------
    1 Li    1.1526049    0.0000000    0.0000000
            0.0000000    1.1526049    0.0000000
            0.0000000    0.0000000    0.9775342
    2 Li    1.1526049    0.0000000    0.0000000
            0.0000000    1.1526049    0.0000000
            0.0000000    0.0000000    0.9775342
    3 Cl   -1.1747772    0.0000000    0.0000000
            0.0000000   -1.1747772    0.0000000
            0.0000000    0.0000000   -1.2902662
    4 Cl   -1.1747772    0.0000000    0.0000000
            0.0000000   -1.1747772    0.0000000
            0.0000000    0.0000000   -1.2902662
    5 Cl   -1.3011168    0.0000000    0.0000000
            0.0000000   -1.3011168    0.0000000
            0.0000000    0.0000000   -0.9574103
    6 Cl   -1.3011168    0.0000000    0.0000000
            0.0000000   -1.3011168    0.0000000
            0.0000000    0.0000000   -0.9574103
    7 Cs    1.3232891    0.0000000    0.0000000
            0.0000000    1.3232891    0.0000000
            0.0000000    0.0000000    1.2701423
    8 Cs    1.3232891    0.0000000    0.0000000
            0.0000000    1.3232891    0.0000000
            0.0000000    0.0000000    1.2701423
----------------------------------------------------------------------------
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 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 05:24:03]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 05:24:04]-------------------------
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    4.842505820000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.842505820000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.364425330000000
Atomic positions (fractional):
    1 Li  0.50000000000000  0.00000000000000  0.08114733798113   6.941
    2 Li  0.00000000000000  0.50000000000000  0.91885266201887   6.941
    3 Cl  0.00000000000000  0.50000000000000  0.66980280217486  35.453
    4 Cl  0.50000000000000  0.00000000000000  0.33019719782514  35.453
    5 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453
    6 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453
    7 Cs  0.00000000000000  0.50000000000000  0.30233586163350 132.905
    8 Cs  0.50000000000000  0.00000000000000  0.69766413836650 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   14.527517460000002    0.000000000000000    0.000000000000000
  b    0.000000000000000   14.527517460000002    0.000000000000000
  c    0.000000000000000    0.000000000000000   18.728850659999999
Atomic positions (fractional):
    1 Li  0.16666666666667  0.00000000000000  0.04057366899056   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.04057366899056   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.04057366899056   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.04057366899056   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.04057366899056   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.04057366899056   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.04057366899056   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.04057366899056   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.04057366899056   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.54057366899056   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.54057366899056   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.54057366899056   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.54057366899056   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.54057366899056   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.54057366899056   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.54057366899056   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.54057366899056   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.54057366899056   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.45942633100944   6.941 > 19
   20 Li  0.33333333333333  0.16666666666667  0.45942633100944   6.941 > 19
   21 Li  0.66666666666667  0.16666666666667  0.45942633100944   6.941 > 19
   22 Li  0.00000000000000  0.50000000000000  0.45942633100944   6.941 > 19
   23 Li  0.33333333333333  0.50000000000000  0.45942633100944   6.941 > 19
   24 Li  0.66666666666667  0.50000000000000  0.45942633100944   6.941 > 19
   25 Li  0.00000000000000  0.83333333333333  0.45942633100944   6.941 > 19
   26 Li  0.33333333333333  0.83333333333333  0.45942633100944   6.941 > 19
   27 Li  0.66666666666667  0.83333333333333  0.45942633100944   6.941 > 19
   28 Li  0.00000000000000  0.16666666666667  0.95942633100944   6.941 > 19
   29 Li  0.33333333333333  0.16666666666667  0.95942633100944   6.941 > 19
   30 Li  0.66666666666667  0.16666666666667  0.95942633100944   6.941 > 19
   31 Li  0.00000000000000  0.50000000000000  0.95942633100944   6.941 > 19
   32 Li  0.33333333333333  0.50000000000000  0.95942633100944   6.941 > 19
   33 Li  0.66666666666667  0.50000000000000  0.95942633100944   6.941 > 19
   34 Li  0.00000000000000  0.83333333333333  0.95942633100944   6.941 > 19
   35 Li  0.33333333333333  0.83333333333333  0.95942633100944   6.941 > 19
   36 Li  0.66666666666667  0.83333333333333  0.95942633100944   6.941 > 19
   37 Cl  0.00000000000000  0.16666666666667  0.33490140108743  35.453 > 37
   38 Cl  0.33333333333333  0.16666666666667  0.33490140108743  35.453 > 37
   39 Cl  0.66666666666667  0.16666666666667  0.33490140108743  35.453 > 37
   40 Cl  0.00000000000000  0.50000000000000  0.33490140108743  35.453 > 37
   41 Cl  0.33333333333333  0.50000000000000  0.33490140108743  35.453 > 37
   42 Cl  0.66666666666667  0.50000000000000  0.33490140108743  35.453 > 37
   43 Cl  0.00000000000000  0.83333333333333  0.33490140108743  35.453 > 37
   44 Cl  0.33333333333333  0.83333333333333  0.33490140108743  35.453 > 37
   45 Cl  0.66666666666667  0.83333333333333  0.33490140108743  35.453 > 37
   46 Cl  0.00000000000000  0.16666666666667  0.83490140108743  35.453 > 37
   47 Cl  0.33333333333333  0.16666666666667  0.83490140108743  35.453 > 37
   48 Cl  0.66666666666667  0.16666666666667  0.83490140108743  35.453 > 37
   49 Cl  0.00000000000000  0.50000000000000  0.83490140108743  35.453 > 37
   50 Cl  0.33333333333333  0.50000000000000  0.83490140108743  35.453 > 37
   51 Cl  0.66666666666667  0.50000000000000  0.83490140108743  35.453 > 37
   52 Cl  0.00000000000000  0.83333333333333  0.83490140108743  35.453 > 37
   53 Cl  0.33333333333333  0.83333333333333  0.83490140108743  35.453 > 37
   54 Cl  0.66666666666667  0.83333333333333  0.83490140108743  35.453 > 37
   55 Cl  0.16666666666667  0.00000000000000  0.16509859891257  35.453 > 55
   56 Cl  0.50000000000000  0.00000000000000  0.16509859891257  35.453 > 55
   57 Cl  0.83333333333333  0.00000000000000  0.16509859891257  35.453 > 55
   58 Cl  0.16666666666667  0.33333333333333  0.16509859891257  35.453 > 55
   59 Cl  0.50000000000000  0.33333333333333  0.16509859891257  35.453 > 55
   60 Cl  0.83333333333333  0.33333333333333  0.16509859891257  35.453 > 55
   61 Cl  0.16666666666667  0.66666666666667  0.16509859891257  35.453 > 55
   62 Cl  0.50000000000000  0.66666666666667  0.16509859891257  35.453 > 55
   63 Cl  0.83333333333333  0.66666666666667  0.16509859891257  35.453 > 55
   64 Cl  0.16666666666667  0.00000000000000  0.66509859891257  35.453 > 55
   65 Cl  0.50000000000000  0.00000000000000  0.66509859891257  35.453 > 55
   66 Cl  0.83333333333333  0.00000000000000  0.66509859891257  35.453 > 55
   67 Cl  0.16666666666667  0.33333333333333  0.66509859891257  35.453 > 55
   68 Cl  0.50000000000000  0.33333333333333  0.66509859891257  35.453 > 55
   69 Cl  0.83333333333333  0.33333333333333  0.66509859891257  35.453 > 55
   70 Cl  0.16666666666667  0.66666666666667  0.66509859891257  35.453 > 55
   71 Cl  0.50000000000000  0.66666666666667  0.66509859891257  35.453 > 55
   72 Cl  0.83333333333333  0.66666666666667  0.66509859891257  35.453 > 55
   73 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 73
   74 Cl  0.33333333333333  0.00000000000000  0.00000000000000  35.453 > 73
   75 Cl  0.66666666666667  0.00000000000000  0.00000000000000  35.453 > 73
   76 Cl  0.00000000000000  0.33333333333333  0.00000000000000  35.453 > 73
   77 Cl  0.33333333333333  0.33333333333333  0.00000000000000  35.453 > 73
   78 Cl  0.66666666666667  0.33333333333333  0.00000000000000  35.453 > 73
   79 Cl  0.00000000000000  0.66666666666667  0.00000000000000  35.453 > 73
   80 Cl  0.33333333333333  0.66666666666667  0.00000000000000  35.453 > 73
   81 Cl  0.66666666666667  0.66666666666667  0.00000000000000  35.453 > 73
   82 Cl  0.00000000000000  0.00000000000000  0.50000000000000  35.453 > 73
   83 Cl  0.33333333333333  0.00000000000000  0.50000000000000  35.453 > 73
   84 Cl  0.66666666666667  0.00000000000000  0.50000000000000  35.453 > 73
   85 Cl  0.00000000000000  0.33333333333333  0.50000000000000  35.453 > 73
   86 Cl  0.33333333333333  0.33333333333333  0.50000000000000  35.453 > 73
   87 Cl  0.66666666666667  0.33333333333333  0.50000000000000  35.453 > 73
   88 Cl  0.00000000000000  0.66666666666667  0.50000000000000  35.453 > 73
   89 Cl  0.33333333333333  0.66666666666667  0.50000000000000  35.453 > 73
   90 Cl  0.66666666666667  0.66666666666667  0.50000000000000  35.453 > 73
   91 Cl  0.16666666666667  0.16666666666667  0.00000000000000  35.453 > 91
   92 Cl  0.50000000000000  0.16666666666667  0.00000000000000  35.453 > 91
   93 Cl  0.83333333333333  0.16666666666667  0.00000000000000  35.453 > 91
   94 Cl  0.16666666666667  0.50000000000000  0.00000000000000  35.453 > 91
   95 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 91
   96 Cl  0.83333333333333  0.50000000000000  0.00000000000000  35.453 > 91
   97 Cl  0.16666666666667  0.83333333333333  0.00000000000000  35.453 > 91
   98 Cl  0.50000000000000  0.83333333333333  0.00000000000000  35.453 > 91
   99 Cl  0.83333333333333  0.83333333333333  0.00000000000000  35.453 > 91
  100 Cl  0.16666666666667  0.16666666666667  0.50000000000000  35.453 > 91
  101 Cl  0.50000000000000  0.16666666666667  0.50000000000000  35.453 > 91
  102 Cl  0.83333333333333  0.16666666666667  0.50000000000000  35.453 > 91
  103 Cl  0.16666666666667  0.50000000000000  0.50000000000000  35.453 > 91
  104 Cl  0.50000000000000  0.50000000000000  0.50000000000000  35.453 > 91
  105 Cl  0.83333333333333  0.50000000000000  0.50000000000000  35.453 > 91
  106 Cl  0.16666666666667  0.83333333333333  0.50000000000000  35.453 > 91
  107 Cl  0.50000000000000  0.83333333333333  0.50000000000000  35.453 > 91
  108 Cl  0.83333333333333  0.83333333333333  0.50000000000000  35.453 > 91
  109 Cs  0.00000000000000  0.16666666666667  0.15116793081675 132.905 > 109
  110 Cs  0.33333333333333  0.16666666666667  0.15116793081675 132.905 > 109
  111 Cs  0.66666666666667  0.16666666666667  0.15116793081675 132.905 > 109
  112 Cs  0.00000000000000  0.50000000000000  0.15116793081675 132.905 > 109
  113 Cs  0.33333333333333  0.50000000000000  0.15116793081675 132.905 > 109
  114 Cs  0.66666666666667  0.50000000000000  0.15116793081675 132.905 > 109
  115 Cs  0.00000000000000  0.83333333333333  0.15116793081675 132.905 > 109
  116 Cs  0.33333333333333  0.83333333333333  0.15116793081675 132.905 > 109
  117 Cs  0.66666666666667  0.83333333333333  0.15116793081675 132.905 > 109
  118 Cs  0.00000000000000  0.16666666666667  0.65116793081675 132.905 > 109
  119 Cs  0.33333333333333  0.16666666666667  0.65116793081675 132.905 > 109
  120 Cs  0.66666666666667  0.16666666666667  0.65116793081675 132.905 > 109
  121 Cs  0.00000000000000  0.50000000000000  0.65116793081675 132.905 > 109
  122 Cs  0.33333333333333  0.50000000000000  0.65116793081675 132.905 > 109
  123 Cs  0.66666666666667  0.50000000000000  0.65116793081675 132.905 > 109
  124 Cs  0.00000000000000  0.83333333333333  0.65116793081675 132.905 > 109
  125 Cs  0.33333333333333  0.83333333333333  0.65116793081675 132.905 > 109
  126 Cs  0.66666666666667  0.83333333333333  0.65116793081675 132.905 > 109
  127 Cs  0.16666666666667  0.00000000000000  0.34883206918325 132.905 > 127
  128 Cs  0.50000000000000  0.00000000000000  0.34883206918325 132.905 > 127
  129 Cs  0.83333333333333  0.00000000000000  0.34883206918325 132.905 > 127
  130 Cs  0.16666666666667  0.33333333333333  0.34883206918325 132.905 > 127
  131 Cs  0.50000000000000  0.33333333333333  0.34883206918325 132.905 > 127
  132 Cs  0.83333333333333  0.33333333333333  0.34883206918325 132.905 > 127
  133 Cs  0.16666666666667  0.66666666666667  0.34883206918325 132.905 > 127
  134 Cs  0.50000000000000  0.66666666666667  0.34883206918325 132.905 > 127
  135 Cs  0.83333333333333  0.66666666666667  0.34883206918325 132.905 > 127
  136 Cs  0.16666666666667  0.00000000000000  0.84883206918325 132.905 > 127
  137 Cs  0.50000000000000  0.00000000000000  0.84883206918325 132.905 > 127
  138 Cs  0.83333333333333  0.00000000000000  0.84883206918325 132.905 > 127
  139 Cs  0.16666666666667  0.33333333333333  0.84883206918325 132.905 > 127
  140 Cs  0.50000000000000  0.33333333333333  0.84883206918325 132.905 > 127
  141 Cs  0.83333333333333  0.33333333333333  0.84883206918325 132.905 > 127
  142 Cs  0.16666666666667  0.66666666666667  0.84883206918325 132.905 > 127
  143 Cs  0.50000000000000  0.66666666666667  0.84883206918325 132.905 > 127
  144 Cs  0.83333333333333  0.66666666666667  0.84883206918325 132.905 > 127
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.7762129    0.0000000    0.0000000
            0.0000000    2.7762129    0.0000000
            0.0000000    0.0000000    2.7028699
-------------------------- Born effective charges --------------------------
    1 Li    1.1526049    0.0000000    0.0000000
            0.0000000    1.1526049    0.0000000
            0.0000000    0.0000000    0.9775342
    2 Li    1.1526049    0.0000000    0.0000000
            0.0000000    1.1526049    0.0000000
            0.0000000    0.0000000    0.9775342
    3 Cl   -1.1747772    0.0000000    0.0000000
            0.0000000   -1.1747772    0.0000000
            0.0000000    0.0000000   -1.2902662
    4 Cl   -1.1747772    0.0000000    0.0000000
            0.0000000   -1.1747772    0.0000000
            0.0000000    0.0000000   -1.2902662
    5 Cl   -1.3011168    0.0000000    0.0000000
            0.0000000   -1.3011168    0.0000000
            0.0000000    0.0000000   -0.9574103
    6 Cl   -1.3011168    0.0000000    0.0000000
            0.0000000   -1.3011168    0.0000000
            0.0000000    0.0000000   -0.9574103
    7 Cs    1.3232891    0.0000000    0.0000000
            0.0000000    1.3232891    0.0000000
            0.0000000    0.0000000    1.2701423
    8 Cs    1.3232891    0.0000000    0.0000000
            0.0000000    1.3232891    0.0000000
            0.0000000    0.0000000    1.2701423
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 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.69, 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.974   (   0.000    0.000    0.000)    0.000
   0.974   (  -0.000    0.000    0.000)    0.000
   1.835   (   0.000    0.000    0.000)    0.000
   1.945   (   0.000    0.000   -0.000)    0.000
   1.945   (   0.000    0.000    0.000)    0.000
   2.291   (   0.000    0.000    0.000)    0.000
   2.325   (   0.000    0.000    0.000)    0.000
   3.340   (  -0.000    0.000   -0.000)    0.000
   3.340   (  -0.000    0.000    0.000)    0.000
   3.373   (  -0.000    0.000   -0.000)    0.000
   3.373   (   0.000    0.000    0.000)    0.000
   3.837   (   0.000    0.000   -0.000)    0.000
   4.808   (   0.000    0.000    0.000)    0.000
   5.164   (   0.000    0.000   -0.000)    0.000
   5.164   (   0.000    0.000   -0.000)    0.000
   6.199   (   0.000    0.000    0.000)    0.000
   6.199   (  -0.000    0.000   -0.000)    0.000
   6.408   (   0.000    0.000    0.000)    0.000
   6.408   (   0.000   -0.000   -0.000)    0.000
  10.981   (   0.000    0.000    0.000)    0.000
  11.397   (   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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.266   (  13.551    0.000    0.000)   13.551
   0.548   (  24.629    0.000    0.000)   24.629
   0.640   (  29.061    0.000    0.000)   29.061
   1.049   (   6.773    0.000    0.000)    6.773
   1.117   (  12.425    0.000    0.000)   12.425
   1.795   (  -3.391    0.000    0.000)    3.391
   2.090   (  13.606    0.000    0.000)   13.606
   2.295   (  14.535    0.000    0.000)   14.535
   2.297   (   0.568    0.000    0.000)    0.568
   2.355   (   2.880    0.000    0.000)    2.880
   3.334   (  -0.577    0.000    0.000)    0.577
   3.354   (   1.312    0.000    0.000)    1.312
   3.377   (   0.508    0.000    0.000)    0.508
   3.725   (   2.415    0.000    0.000)    2.415
   3.857   (   2.162    0.000    0.000)    2.162
   4.783   (  -2.382    0.000    0.000)    2.382
   5.132   (  -3.224    0.000    0.000)    3.224
   5.254   (   7.376    0.000    0.000)    7.376
   6.205   (   0.555    0.000    0.000)    0.555
   6.401   (  -0.712    0.000    0.000)    0.712
   6.460   (   5.741    0.000    0.000)    5.741
   8.815   (   0.252    0.000    0.000)    0.252
  10.954   (  -2.500    0.000    0.000)    2.500
  11.368   (  -2.900    0.000    0.000)    2.900
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.555   (  14.197    0.000    0.000)   14.197
   0.995   (  18.489    0.000    0.000)   18.489
   1.172   (  21.934    0.000    0.000)   21.934
   1.225   (   9.511    0.000    0.000)    9.511
   1.422   (  15.809    0.000    0.000)   15.809
   1.716   (  -3.608    0.000    0.000)    3.608
   2.314   (   1.024    0.000    0.000)    1.024
   2.409   (   2.497    0.000    0.000)    2.497
   2.468   (  21.467    0.000    0.000)   21.467
   2.719   (  24.664    0.000    0.000)   24.664
   3.317   (  -1.122    0.000    0.000)    1.122
   3.392   (   2.318    0.000    0.000)    2.318
   3.401   (   2.086    0.000    0.000)    2.086
   3.774   (   2.191    0.000    0.000)    2.191
   3.935   (   5.305    0.000    0.000)    5.305
   4.712   (  -4.420    0.000    0.000)    4.420
   5.028   (  -6.865    0.000    0.000)    6.865
   5.381   (   2.738    0.000    0.000)    2.738
   6.221   (   0.996    0.000    0.000)    0.996
   6.380   (  -1.180    0.000    0.000)    1.180
   6.692   (  17.604    0.000    0.000)   17.604
   8.809   (  -1.114    0.000    0.000)    1.114
  10.886   (  -3.756    0.000    0.000)    3.756
  11.277   (  -5.818    0.000    0.000)    5.818
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.841   (  13.230    0.000    0.000)   13.230
   1.314   (  12.503    0.000    0.000)   12.503
   1.415   (   8.331    0.000    0.000)    8.331
   1.475   (   3.878    0.000    0.000)    3.878
   1.724   (   8.033    0.000    0.000)    8.033
   1.731   (  13.568    0.000    0.000)   13.568
   2.336   (   0.952    0.000    0.000)    0.952
   2.436   (  -0.158    0.000    0.000)    0.158
   2.898   (  17.767    0.000    0.000)   17.767
   3.250   (  24.562    0.000    0.000)   24.562
   3.289   (  -1.554    0.000    0.000)    1.554
   3.449   (   3.131    0.000    0.000)    3.131
   3.497   (   8.694    0.000    0.000)    8.694
   3.825   (   3.409    0.000    0.000)    3.409
   4.075   (   8.469    0.000    0.000)    8.469
   4.602   (  -6.253    0.000    0.000)    6.253
   4.844   ( -11.079    0.000    0.000)   11.079
   5.325   (  -8.011    0.000    0.000)    8.011
   6.245   (   1.252    0.000    0.000)    1.252
   6.354   (  -1.340    0.000    0.000)    1.340
   7.174   (  27.581    0.000    0.000)   27.581
   8.753   (  -4.779    0.000    0.000)    4.779
  10.812   (  -3.099    0.000    0.000)    3.099
  11.131   (  -8.183    0.000    0.000)    8.183
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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|
   1.092   (  10.922    0.000    0.000)   10.922
   1.426   (  -5.325    0.000    0.000)    5.325
   1.513   (   6.716    0.000    0.000)    6.716
   1.551   (   4.515    0.000    0.000)    4.515
   1.946   (   9.772    0.000    0.000)    9.772
   1.964   (   8.641    0.000    0.000)    8.641
   2.346   (   0.008    0.000    0.000)    0.008
   2.402   (  -2.747    0.000    0.000)    2.747
   3.135   (   5.997    0.000    0.000)    5.997
   3.254   (  -1.833    0.000    0.000)    1.833
   3.520   (   3.716    0.000    0.000)    3.716
   3.588   (   5.695    0.000    0.000)    5.695
   3.797   (  19.095    0.000    0.000)   19.095
   3.992   (  13.485    0.000    0.000)   13.485
   4.308   (  15.031    0.000    0.000)   15.031
   4.453   (  -8.227    0.000    0.000)    8.227
   4.569   ( -15.558    0.000    0.000)   15.558
   5.071   ( -16.065    0.000    0.000)   16.065
   6.272   (   1.330    0.000    0.000)    1.330
   6.326   (  -1.332    0.000    0.000)    1.332
   7.760   (  27.662    0.000    0.000)   27.662
   8.592   ( -11.367    0.000    0.000)   11.367
  10.774   (  -0.142    0.000    0.000)    0.142
  10.956   (  -8.167    0.000    0.000)    8.167
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.288   (  -8.044    0.000    0.000)    8.044
   1.288   (   8.044    0.000    0.000)    8.044
   1.591   (  -0.837    0.000    0.000)    0.837
   1.591   (   0.837    0.000    0.000)    0.837
   2.061   (  -0.316    0.000    0.000)    0.316
   2.061   (   0.316    0.000    0.000)    0.316
   2.351   (  -1.307    0.000    0.000)    1.307
   2.351   (   1.307    0.000    0.000)    1.307
   3.211   (  -2.416    0.000    0.000)    2.416
   3.211   (   2.416    0.000    0.000)    2.416
   3.592   (  -2.739    0.000    0.000)    2.739
   3.592   (   2.739    0.000    0.000)    2.739
   4.206   ( -19.307    0.000    0.000)   19.307
   4.206   (  19.307    0.000    0.000)   19.307
   4.257   ( -10.960    0.000    0.000)   10.960
   4.257   (  10.960    0.000    0.000)   10.960
   4.688   ( -20.021    0.000    0.000)   20.021
   4.688   (  20.021    0.000    0.000)   20.021
   6.299   (  -1.322    0.000    0.000)    1.322
   6.299   (   1.322    0.000    0.000)    1.322
   8.265   ( -20.435    0.000    0.000)   20.435
   8.265   (  20.435    0.000    0.000)   20.435
  10.819   (  -4.584    0.000    0.000)    4.584
  10.819   (   4.584    0.000    0.000)    4.584
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.167   (   5.608    5.608    0.000)    7.931
   0.383   (   9.766    9.766    0.000)   13.811
   1.038   (   3.003    3.003    0.000)    4.247
   1.099   (  21.207   21.207    0.000)   29.992
   1.287   (  11.635   11.635    0.000)   16.454
   1.765   (  -2.602   -2.602    0.000)    3.680
   1.961   (   0.603    0.603    0.000)    0.853
   2.304   (   0.629    0.629    0.000)    0.890
   2.375   (   2.094    2.094    0.000)    2.961
   2.673   (  22.851   22.851    0.000)   32.316
   3.331   (  -0.386   -0.386    0.000)    0.547
   3.364   (   1.022    1.022    0.000)    1.445
   3.401   (   1.243    1.243    0.000)    1.758
   3.718   (   0.400    0.400    0.000)    0.566
   3.882   (   2.557    2.557    0.000)    3.616
   4.757   (  -2.435   -2.435    0.000)    3.443
   5.180   (  -2.354   -2.354    0.000)    3.329
   5.217   (   2.584    2.584    0.000)    3.654
   6.119   (  -3.276   -3.276    0.000)    4.632
   6.418   (   0.126    0.126    0.000)    0.178
   6.525   (   8.197    8.197    0.000)   11.593
   8.894   (   3.925    3.925    0.000)    5.550
  10.929   (  -2.323   -2.323    0.000)    3.285
  11.341   (  -2.623   -2.623    0.000)    3.710
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.619   (  12.408    5.856    0.000)   13.721
   0.638   (  25.242  -19.680    0.000)   32.007
   1.175   (   9.277   -1.853    0.000)    9.460
   1.449   (  11.710   12.667    0.000)   17.250
   1.538   (  12.055    8.260    0.000)   14.613
   1.714   (  -0.874    0.799    0.000)    1.184
   2.135   (  14.931  -14.499    0.000)   20.813
   2.321   (   0.923    0.600    0.000)    1.100
   2.434   (   3.428    0.884    0.000)    3.540
   3.139   (  19.502   24.121    0.000)   31.019
   3.319   (  -0.684    0.282    0.000)    0.740
   3.395   (   1.946    0.336    0.000)    1.975
   3.437   (   4.115    2.868    0.000)    5.016
   3.746   (   2.169   -2.520    0.000)    3.325
   3.974   (   6.786    4.751    0.000)    8.284
   4.684   (  -4.604   -2.658    0.000)    5.316
   5.010   ( -11.468   -4.648    0.000)   12.375
   5.327   (   4.165   -1.957    0.000)    4.602
   6.104   (   1.091   -7.823    0.000)    7.899
   6.398   (  -1.826    0.721    0.000)    1.963
   6.825   (  20.670   13.045    0.000)   24.442
   8.939   (   0.455    9.646    0.000)    9.657
  10.866   (  -3.493   -1.888    0.000)    3.971
  11.258   (  -5.403   -1.912    0.000)    5.732
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.876   (  12.097    3.282    0.000)   12.534
   1.099   (  18.818  -14.848    0.000)   23.970
   1.380   (   9.637   -2.608    0.000)    9.983
   1.523   (  -2.473    0.687    0.000)    2.567
   1.766   (   9.697    2.524    0.000)   10.020
   1.807   (   7.926    6.416    0.000)   10.197
   2.337   (   0.555    0.156    0.000)    0.576
   2.413   (   4.167   -4.522    0.000)    6.150
   2.611   (  19.619  -16.122    0.000)   25.393
   3.303   (  -0.915    1.371    0.000)    1.648
   3.335   (   2.967    5.151    0.000)    5.945
   3.444   (   2.769   -0.401    0.000)    2.798
   3.622   (   7.357    4.386    0.000)    8.565
   3.842   (   8.828    0.319    0.000)    8.833
   4.191   (  15.039   14.034    0.000)   20.570
   4.567   (  -7.008   -3.444    0.000)    7.809
   4.721   ( -16.195  -12.423    0.000)   20.411
   5.349   (  -2.291    3.515    0.000)    4.196
   6.141   (   2.201   -7.714    0.000)    8.022
   6.347   (  -2.923   -0.583    0.000)    2.980
   7.347   (  28.600   15.803    0.000)   32.675
   8.908   (  -3.764   12.837    0.000)   13.377
  10.797   (  -2.851   -1.395    0.000)    3.174
  11.119   (  -7.868   -1.147    0.000)    7.951
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.107   (  10.091    1.425    0.000)   10.191
   1.380   (   3.862   -5.827    0.000)    6.991
   1.468   (   2.025   -2.498    0.000)    3.216
   1.539   (   5.168   -1.115    0.000)    5.287
   1.931   (   6.084   -3.450    0.000)    6.994
   1.948   (   5.095   -0.050    0.000)    5.095
   2.343   (   0.249    0.107    0.000)    0.271
   2.412   (  -2.113    0.897    0.000)    2.295
   3.049   (  17.429  -12.571    0.000)   21.490
   3.280   (  -1.428    2.307    0.000)    2.713
   3.425   (   8.438   -7.722    0.000)   11.438
   3.508   (   3.301   -1.126    0.000)    3.487
   3.674   (   0.480   -4.062    0.000)    4.090
   4.055   (  10.228   -1.474    0.000)   10.333
   4.303   ( -20.051  -18.189    0.000)   27.071
   4.441   (  -9.025   -5.058    0.000)   10.346
   4.572   (  19.999   21.345    0.000)   29.250
   5.227   (  -9.600   11.783    0.000)   15.198
   6.187   (   2.207   -6.446    0.000)    6.813
   6.286   (  -2.844   -2.998    0.000)    4.132
   7.945   (  27.844   16.484    0.000)   32.357
   8.764   ( -10.661   14.917    0.000)   18.335
  10.764   (   0.029   -1.005    0.000)    1.005
  10.948   (  -8.097   -0.750    0.000)    8.132
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.288   (  -7.343   -0.020    0.000)    7.343
   1.288   (   7.343   -0.020    0.000)    7.343
   1.571   (  -2.418   -2.040    0.000)    3.163
   1.571   (   2.418   -2.040    0.000)    3.163
   2.001   (  -0.332   -5.542    0.000)    5.552
   2.001   (   0.332   -5.542    0.000)    5.552
   2.363   (  -1.952    1.326    0.000)    2.360
   2.363   (   1.952    1.326    0.000)    2.360
   3.235   (  -3.498    1.543    0.000)    3.823
   3.235   (   3.498    1.543    0.000)    3.823
   3.564   (  -0.919   -3.001    0.000)    3.138
   3.564   (   0.919   -3.001    0.000)    3.138
   3.860   ( -20.472  -22.152    0.000)   30.163
   3.860   (  20.472  -22.152    0.000)   30.163
   4.256   (  -9.307   -3.327    0.000)    9.883
   4.256   (   9.307   -3.327    0.000)    9.883
   4.956   ( -16.371   18.729    0.000)   24.875
   4.956   (  16.371   18.729    0.000)   24.875
   6.233   (  -2.300   -5.043    0.000)    5.543
   6.233   (   2.300   -5.043    0.000)    5.543
   8.448   ( -20.150   16.159    0.000)   25.829
   8.448   (  20.150   16.159    0.000)   25.829
  10.811   (  -4.672   -0.774    0.000)    4.736
  10.811   (   4.672   -0.774    0.000)    4.736
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.452   (   7.698    7.698    0.000)   10.887
   0.777   (   8.891    8.891    0.000)   12.573
   1.200   (   4.391    4.391    0.000)    6.210
   1.559   (  -0.823   -0.823    0.000)    1.165
   1.700   (   6.750    6.750    0.000)    9.546
   1.814   (   7.396    7.396    0.000)   10.460
   1.975   (  -0.154   -0.154    0.000)    0.218
   2.334   (   0.624    0.624    0.000)    0.882
   2.467   (   2.087    2.087    0.000)    2.952
   3.330   (   0.741    0.741    0.000)    1.048
   3.407   (   0.828    0.828    0.000)    1.171
   3.463   (   1.613    1.613    0.000)    2.282
   3.555   (  11.286   11.286    0.000)   15.961
   3.698   (  -0.624   -0.624    0.000)    0.883
   4.188   (  16.692   16.692    0.000)   23.605
   4.604   (  -5.042   -5.042    0.000)    7.130
   4.763   ( -17.702  -17.702    0.000)   25.035
   5.377   (   5.070    5.070    0.000)    7.170
   6.002   (  -1.576   -1.576    0.000)    2.228
   6.389   (  -1.646   -1.646    0.000)    2.328
   7.215   (  23.984   23.984    0.000)   33.918
   9.114   (   6.370    6.370    0.000)    9.009
  10.815   (  -2.859   -2.859    0.000)    4.043
  11.195   (  -4.254   -4.254    0.000)    6.016
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.811   (  20.855  -10.728    0.000)   23.453
   0.968   (   9.309    5.422    0.000)   10.773
   1.339   (   8.401   -0.728    0.000)    8.432
   1.508   (  -3.373   -1.633    0.000)    3.747
   1.800   (   2.696    0.160    0.000)    2.700
   1.921   (   2.890    3.764    0.000)    4.745
   2.136   (  14.285  -15.037    0.000)   20.741
   2.346   (   0.603    0.851    0.000)    1.043
   2.509   (   1.934    0.814    0.000)    2.098
   3.344   (   0.351    2.498    0.000)    2.522
   3.433   (   1.762   -0.702    0.000)    1.897
   3.446   (  -1.927    5.200    0.000)    5.545
   3.622   (  -0.940   -2.823    0.000)    2.976
   3.807   (   7.970   -2.633    0.000)    8.393
   4.331   ( -23.173  -22.402    0.000)   32.231
   4.483   (  -6.916   -6.255    0.000)    9.325
   4.608   (  21.878   23.080    0.000)   31.801
   5.471   (   3.182    7.366    0.000)    8.024
   6.019   (   2.437   -3.080    0.000)    3.928
   6.331   (  -3.912   -1.084    0.000)    4.060
   7.777   (  29.309   24.735    0.000)   38.352
   9.187   (   0.453   12.547    0.000)   12.555
  10.759   (  -2.292   -2.186    0.000)    3.167
  11.080   (  -6.904   -2.832    0.000)    7.463
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.148   (   7.898    2.472    0.000)    8.276
   1.197   (  15.046  -10.326    0.000)   18.248
   1.434   (  -3.125   -1.220    0.000)    3.355
   1.505   (   5.971   -2.196    0.000)    6.362
   1.815   (  -0.412   -7.702    0.000)    7.713
   1.930   (  -2.052   -2.062    0.000)    2.909
   2.367   (   1.615    2.345    0.000)    2.847
   2.412   (   4.711   -4.093    0.000)    6.241
   2.650   (  18.350  -18.338    0.000)   25.942
   3.332   (  -2.033    2.019    0.000)    2.865
   3.412   (  -0.909    3.165    0.000)    3.293
   3.475   (   1.909   -2.109    0.000)    2.845
   3.594   (  -1.637   -3.850    0.000)    4.184
   3.829   ( -24.231  -25.422    0.000)   35.120
   3.975   (   8.387   -5.356    0.000)    9.951
   4.327   (  -8.126   -6.509    0.000)   10.411
   5.033   (  18.447   21.912    0.000)   28.643
   5.483   (  -2.605   12.260    0.000)   12.534
   6.079   (   3.117   -3.060    0.000)    4.369
   6.237   (  -4.809   -1.413    0.000)    5.013
   8.374   (  27.208   23.675    0.000)   36.066
   9.118   (  -7.609   17.536    0.000)   19.115
  10.735   (   0.352   -1.765    0.000)    1.800
  10.922   (  -7.787   -1.893    0.000)    8.014
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.285   (  -4.850   -0.362    0.000)    4.863
   1.285   (   4.850   -0.362    0.000)    4.863
   1.498   (  -6.668   -5.103    0.000)    8.396
   1.498   (   6.668   -5.103    0.000)    8.396
   1.853   (  -4.312   -8.120    0.000)    9.194
   1.853   (   4.312   -8.120    0.000)    9.194
   2.409   (  -2.038    2.933    0.000)    3.572
   2.409   (   2.038    2.933    0.000)    3.572
   3.126   ( -21.728  -19.628    0.000)   29.281
   3.126   (  21.728  -19.628    0.000)   29.281
   3.430   (  -4.015   -3.679    0.000)    5.446
   3.430   (   4.015   -3.679    0.000)    5.446
   3.562   (  -0.565   -4.675    0.000)    4.709
   3.562   (   0.565   -4.675    0.000)    4.709
   4.152   (  -8.671   -6.389    0.000)   10.771
   4.152   (   8.671   -6.389    0.000)   10.771
   5.344   ( -11.109   17.935    0.000)   21.096
   5.344   (  11.109   17.935    0.000)   21.096
   6.147   (  -3.710   -2.509    0.000)    4.479
   6.147   (   3.710   -2.509    0.000)    4.479
   8.852   ( -18.330   21.310    0.000)   28.109
   8.852   (  18.330   21.310    0.000)   28.109
  10.786   (  -4.780   -1.659    0.000)    5.060
  10.786   (   4.780   -1.659    0.000)    5.060
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.759   (   6.728    6.728    0.000)    9.514
   1.089   (   5.971    5.971    0.000)    8.444
   1.363   (   3.179    3.179    0.000)    4.496
   1.464   (  -2.458   -2.458    0.000)    3.476
   1.755   (  -4.617   -4.617    0.000)    6.530
   1.938   (  -1.595   -1.595    0.000)    2.255
   1.956   (  -1.789   -1.789    0.000)    2.531
   2.378   (   2.304    2.304    0.000)    3.259
   2.537   (   1.316    1.316    0.000)    1.861
   3.399   (   2.253    2.253    0.000)    3.187
   3.408   (  -1.981   -1.981    0.000)    2.802
   3.527   (   1.444    1.444    0.000)    2.042
   3.546   (  -3.208   -3.208    0.000)    4.536
   3.781   (   1.078    1.078    0.000)    1.525
   3.815   ( -25.892  -25.892    0.000)   36.617
   4.341   (  -7.551   -7.551    0.000)   10.679
   5.070   (  20.881   20.881    0.000)   29.531
   5.619   (   6.369    6.369    0.000)    9.008
   6.012   (   1.780    1.780    0.000)    2.518
   6.296   (  -2.546   -2.546    0.000)    3.601
   8.329   (  27.734   27.734    0.000)   39.222
   9.389   (   6.559    6.559    0.000)    9.276
  10.715   (  -1.846   -1.846    0.000)    2.610
  10.996   (  -5.343   -5.343    0.000)    7.556
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.011   (  14.640   -6.280    0.000)   15.930
   1.204   (   4.861    2.684    0.000)    5.553
   1.409   (  -2.784   -1.286    0.000)    3.067
   1.458   (   4.264   -1.918    0.000)    4.675
   1.626   (  -6.061  -10.140    0.000)   11.813
   1.852   (  -6.871   -5.524    0.000)    8.816
   2.093   (  15.189  -18.049    0.000)   23.590
   2.435   (   2.843    3.890    0.000)    4.819
   2.550   (  -0.065    1.499    0.000)    1.501
   3.140   ( -25.986  -26.148    0.000)   36.864
   3.443   (   0.697    1.300    0.000)    1.475
   3.474   (  -2.347    0.909    0.000)    2.517
   3.489   (  -1.874   -3.628    0.000)    4.084
   3.539   (  -0.945    0.163    0.000)    0.959
   3.859   (   5.937   -5.307    0.000)    7.963
   4.176   (  -8.316   -7.931    0.000)   11.492
   5.448   (  15.233   17.892    0.000)   23.499
   5.713   (   2.152    9.698    0.000)    9.934
   6.064   (   2.878    1.134    0.000)    3.093
   6.222   (  -4.426   -0.459    0.000)    4.450
   8.876   (  24.110   24.052    0.000)   34.055
   9.440   (  -2.146   12.887    0.000)   13.065
  10.695   (   0.326   -1.925    0.000)    1.952
  10.865   (  -6.845   -3.701    0.000)    7.782
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.250   (  -4.119   -4.116    0.000)    5.823
   1.250   (   4.119   -4.116    0.000)    5.823
   1.388   (  -4.243   -3.931    0.000)    5.784
   1.388   (   4.243   -3.931    0.000)    5.784
   1.685   (  -8.341   -7.880    0.000)   11.474
   1.685   (   8.341   -7.880    0.000)   11.474
   2.450   (  -7.155   -4.099    0.000)    8.246
   2.450   (   7.155   -4.099    0.000)    8.246
   2.614   ( -16.934  -20.054    0.000)   26.247
   2.614   (  16.934  -20.054    0.000)   26.247
   3.446   (  -0.519    1.380    0.000)    1.474
   3.446   (   0.519    1.380    0.000)    1.474
   3.509   (  -1.665   -0.385    0.000)    1.709
   3.509   (   1.665   -0.385    0.000)    1.709
   4.006   (  -7.932   -7.468    0.000)   10.894
   4.006   (   7.932   -7.468    0.000)   10.894
   5.677   (  -6.375   14.010    0.000)   15.393
   5.677   (   6.375   14.010    0.000)   15.393
   6.131   (  -3.863    0.496    0.000)    3.894
   6.131   (   3.863    0.496    0.000)    3.894
   9.278   ( -13.899   18.945    0.000)   23.496
   9.278   (  13.899   18.945    0.000)   23.496
  10.742   (  -4.401   -2.585    0.000)    5.104
  10.742   (   4.401   -2.585    0.000)    5.104
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.985   (   4.007    4.007    0.000)    5.666
   1.241   (  -0.256   -0.256    0.000)    0.362
   1.381   (  -1.403   -1.403    0.000)    1.984
   1.420   (  -7.762   -7.762    0.000)   10.976
   1.452   (   1.176    1.176    0.000)    1.664
   1.708   (  -8.432   -8.432    0.000)   11.924
   1.866   (  -1.524   -1.524    0.000)    2.155
   2.515   (   3.474    3.474    0.000)    4.913
   2.550   ( -28.044  -28.044    0.000)   39.660
   2.580   (   0.792    0.792    0.000)    1.120
   3.451   (  -0.608   -0.608    0.000)    0.860
   3.455   (  -1.029   -1.029    0.000)    1.455
   3.495   (   2.573    2.573    0.000)    3.638
   3.578   (   0.983    0.983    0.000)    1.390
   3.788   (  -0.411   -0.411    0.000)    0.581
   4.008   (  -8.025   -8.025    0.000)   11.349
   5.759   (  11.786   11.786    0.000)   16.668
   5.873   (   5.394    5.394    0.000)    7.628
   6.100   (   1.905    1.905    0.000)    2.694
   6.199   (  -2.026   -2.026    0.000)    2.865
   9.327   (  18.567   18.567    0.000)   26.258
   9.626   (   4.566    4.566    0.000)    6.457
  10.665   (  -0.624   -0.624    0.000)    0.883
  10.772   (  -4.956   -4.956    0.000)    7.009
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.131   (  -8.487   -5.494    0.000)   10.110
   1.131   (   8.487   -5.494    0.000)   10.110
   1.356   (  -0.573   -0.049    0.000)    0.575
   1.356   (   0.573   -0.049    0.000)    0.575
   1.540   (  -6.807   -5.683    0.000)    8.867
   1.540   (   6.807   -5.683    0.000)    8.867
   2.050   ( -18.306  -19.356    0.000)   26.641
   2.050   (  18.306  -19.356    0.000)   26.641
   2.572   (  -1.676    2.272    0.000)    2.824
   2.572   (   1.676    2.272    0.000)    2.824
   3.445   (  -0.029   -0.502    0.000)    0.503
   3.445   (   0.029   -0.502    0.000)    0.503
   3.552   (  -2.647    3.257    0.000)    4.196
   3.552   (   2.647    3.257    0.000)    4.196
   3.857   (  -6.163   -6.484    0.000)    8.946
   3.857   (   6.163   -6.484    0.000)    8.946
   5.912   (  -2.468    8.486    0.000)    8.837
   5.912   (   2.468    8.486    0.000)    8.837
   6.144   (  -2.680    0.300    0.000)    2.696
   6.144   (   2.680    0.300    0.000)    2.696
   9.602   (  -7.354   11.722    0.000)   13.838
   9.602   (   7.354   11.722    0.000)   13.838
  10.685   (  -2.833   -2.607    0.000)    3.850
  10.685   (   2.833   -2.607    0.000)    3.850
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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|
   1.071   (   0.000    0.000    0.000)    0.000
   1.071   (   0.000    0.000    0.000)    0.000
   1.360   (   0.000    0.000    0.000)    0.000
   1.360   (   0.000    0.000    0.000)    0.000
   1.476   (   0.000    0.000    0.000)    0.000
   1.476   (   0.000    0.000    0.000)    0.000
   1.832   (   0.000    0.000    0.000)    0.000
   1.832   (   0.000    0.000    0.000)    0.000
   2.598   (   0.000    0.000    0.000)    0.000
   2.598   (   0.000    0.000    0.000)    0.000
   3.440   (   0.000    0.000    0.000)    0.000
   3.440   (   0.000    0.000    0.000)    0.000
   3.600   (   0.000    0.000    0.000)    0.000
   3.600   (   0.000    0.000    0.000)    0.000
   3.774   (   0.000    0.000    0.000)    0.000
   3.774   (   0.000    0.000    0.000)    0.000
   6.008   (  -0.000   -0.000    0.000)    0.000
   6.008   (   0.000    0.000    0.000)    0.000
   6.142   (   0.000    0.000    0.000)    0.000
   6.142   (   0.000    0.000    0.000)    0.000
   9.729   (   0.000    0.000    0.000)    0.000
   9.729   (   0.000    0.000    0.000)    0.000
  10.654   (   0.000    0.000    0.000)    0.000
  10.654   (   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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.262   (   0.000    0.000   11.840)   11.840
   0.262   (   0.000    0.000   11.840)   11.840
   0.719   (   0.000    0.000   32.660)   32.660
   0.925   (  -0.000    0.000   -4.411)    4.411
   0.925   (   0.000    0.000   -4.411)    4.411
   1.849   (   0.000    0.000    0.487)    0.487
   1.947   (   0.000    0.000    0.118)    0.118
   1.947   (   0.000    0.000    0.118)    0.118
   2.333   (  -0.000   -0.000    0.632)    0.632
   2.944   (  -0.000    0.000  -16.260)   16.260
   3.326   (  -0.000    0.000   -0.871)    0.871
   3.326   (  -0.000    0.000   -0.871)    0.871
   3.387   (  -0.000    0.000    0.874)    0.874
   3.387   (   0.000    0.000    0.874)    0.874
   3.982   (  -0.000   -0.000   10.308)   10.308
   4.781   (   0.000   -0.000   -5.833)    5.833
   5.166   (   0.000    0.000    0.157)    0.157
   5.166   (   0.000    0.000    0.157)    0.157
   6.199   (  -0.000    0.000    0.034)    0.034
   6.199   (   0.000    0.000    0.034)    0.034
   6.409   (   0.000    0.000    0.040)    0.040
   6.409   (   0.000    0.000    0.040)    0.040
  11.432   (   0.000    0.000    2.867)    2.867
  12.179   (   0.000    0.000   -1.106)    1.106
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.605   (  21.953    0.000    4.884)   22.490
   0.609   (  15.675    0.000   11.713)   19.568
   0.816   (  14.299    0.000   21.694)   25.982
   1.010   (   7.583    0.000   -3.573)    8.382
   1.079   (  12.729    0.000   -3.194)   13.123
   1.811   (  -3.073    0.000    0.600)    3.131
   2.092   (  13.598    0.000    0.094)   13.599
   2.139   (  14.524    0.000   -2.790)   14.789
   2.358   (   2.774    0.000    0.302)    2.790
   2.526   ( -11.755    0.000    3.396)   12.235
   3.321   (  -0.468    0.000   -0.844)    0.964
   3.372   (   2.025    0.000    1.588)    2.573
   3.390   (   0.397    0.000    0.844)    0.933
   3.643   (  10.337    0.000   -9.035)   13.729
   3.886   (  -3.061    0.000    6.788)    7.446
   4.726   (  -2.975    0.000   -4.972)    5.794
   5.133   (  -3.287    0.000    0.101)    3.289
   5.259   (   7.381    0.000    0.308)    7.387
   6.205   (   0.559    0.000    0.037)    0.560
   6.401   (  -0.722    0.000    0.030)    0.723
   6.461   (   5.814    0.000    0.086)    5.815
   7.486   (  63.384    0.000  -51.824)   81.874
  11.399   (  -3.175    0.000    2.608)    4.109
  11.685   ( -30.139    0.000   22.988)   37.905
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.814   (   8.756    0.000   18.065)   20.075
   1.018   (  17.420    0.000    2.042)   17.539
   1.204   (  10.392    0.000   -1.885)   10.561
   1.229   (  19.157    0.000    5.780)   20.010
   1.386   (  16.069    0.000   -2.917)   16.332
   1.750   (  -2.094    0.000    2.211)    3.045
   2.285   (   0.390    0.000   -1.866)    1.907
   2.407   (   0.453    0.000   -1.143)    1.229
   2.469   (  21.493    0.000    0.127)   21.494
   2.682   (  23.451    0.000   -3.140)   23.660
   3.307   (  -0.972    0.000   -0.733)    1.218
   3.411   (   1.948    0.000    0.740)    2.083
   3.433   (   3.796    0.000    3.453)    5.131
   3.798   (   4.748    0.000   -1.810)    5.081
   3.887   (   3.408    0.000   -0.006)    3.408
   4.651   (  -4.447    0.000   -5.125)    6.786
   5.028   (  -6.951    0.000   -0.032)    6.951
   5.384   (   2.631    0.000    0.239)    2.642
   6.222   (   1.001    0.000    0.045)    1.002
   6.381   (  -1.185    0.000    0.015)    1.186
   6.692   (  17.484    0.000    0.054)   17.484
   8.311   (  22.072    0.000  -31.767)   38.682
  11.215   ( -15.995    0.000   19.406)   25.148
  11.304   (  -5.956    0.000    2.412)    6.426
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.005   (   9.563    0.000   13.130)   16.243
   1.321   (  11.977    0.000    0.561)   11.990
   1.409   (   8.850    0.000   -0.518)    8.865
   1.484   (   3.941    0.000    0.307)    3.953
   1.712   (  14.787    0.000    0.539)   14.797
   1.783   (   6.824    0.000    3.181)    7.529
   2.293   (   0.786    0.000   -3.565)    3.650
   2.404   (  -1.144    0.000   -2.999)    3.210
   2.899   (  17.709    0.000    0.105)   17.709
   3.204   (  24.253    0.000   -3.017)   24.440
   3.282   (  -1.413    0.000   -0.547)    1.515
   3.504   (   8.569    0.000    0.587)    8.590
   3.523   (   4.772    0.000    5.660)    7.403
   3.858   (   2.516    0.000    1.839)    3.116
   4.023   (   9.224    0.000   -3.460)    9.852
   4.541   (  -6.255    0.000   -5.199)    8.134
   4.842   ( -11.117    0.000   -0.145)   11.118
   5.327   (  -8.032    0.000    0.146)    8.033
   6.245   (   1.254    0.000    0.051)    1.255
   6.354   (  -1.333    0.000    0.016)    1.333
   7.171   (  27.499    0.000   -0.184)   27.500
   8.563   (   4.402    0.000  -13.872)   14.554
  10.973   (  -8.218    0.000   11.717)   14.312
  11.158   (  -8.030    0.000    2.227)    8.334
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.196   (   8.607    0.000    8.594)   12.163
   1.462   (  -3.898    0.000    2.770)    4.782
   1.513   (   6.611    0.000    0.026)    6.612
   1.551   (   4.630    0.000    0.028)    4.630
   1.959   (   8.267    0.000    0.620)    8.290
   1.970   (   9.212    0.000    1.066)    9.273
   2.308   (   0.355    0.000   -3.319)    3.338
   2.358   (  -2.933    0.000   -3.678)    4.704
   3.135   (   5.925    0.000   -0.005)    5.925
   3.249   (  -1.710    0.000   -0.342)    1.744
   3.549   (   7.663    0.000    2.396)    8.029
   3.632   (   4.499    0.000    3.934)    5.977
   3.801   (  18.913    0.000    0.340)   18.916
   3.966   (   9.392    0.000   -2.643)    9.757
   4.274   (  13.525    0.000   -2.505)   13.755
   4.402   (  -5.157    0.000   -4.310)    6.721
   4.567   ( -15.496    0.000   -0.129)   15.497
   5.073   ( -16.006    0.000    0.169)   16.007
   6.272   (   1.328    0.000    0.052)    1.329
   6.327   (  -1.322    0.000    0.032)    1.322
   7.757   (  27.593    0.000   -0.215)   27.594
   8.530   (  -7.498    0.000   -4.771)    8.887
  10.861   (  -2.660    0.000    6.809)    7.311
  10.989   (  -7.711    0.000    2.621)    8.144
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.353   (  -6.535    0.000    5.360)    8.451
   1.353   (   6.535    0.000    5.360)    8.451
   1.591   (  -0.894    0.000    0.032)    0.894
   1.591   (   0.894    0.000    0.032)    0.894
   2.065   (  -0.839    0.000    0.380)    0.922
   2.065   (   0.839    0.000    0.380)    0.922
   2.312   (  -0.810    0.000   -3.301)    3.399
   2.312   (   0.810    0.000   -3.301)    3.399
   3.209   (  -2.313    0.000   -0.154)    2.319
   3.209   (   2.313    0.000   -0.154)    2.319
   3.641   (  -2.489    0.000    4.333)    4.997
   3.641   (   2.489    0.000    4.333)    4.997
   4.187   ( -10.266    0.000   -5.992)   11.886
   4.187   (  10.266    0.000   -5.992)   11.886
   4.207   ( -19.157    0.000    0.054)   19.157
   4.207   (  19.157    0.000    0.054)   19.157
   4.692   ( -19.996    0.000    0.252)   19.998
   4.692   (  19.996    0.000    0.252)   19.998
   6.300   (  -1.316    0.000    0.046)    1.317
   6.300   (   1.316    0.000    0.046)    1.317
   8.252   ( -19.272    0.000   -1.002)   19.298
   8.252   (  19.272    0.000   -1.002)   19.298
  10.868   (  -3.398    0.000    3.935)    5.199
  10.868   (   3.398    0.000    3.935)    5.199
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.328   (   3.590    3.590   11.047)   12.158
   0.755   (   4.642    4.642   23.616)   24.512
   0.985   (   2.836    2.836   -4.813)    6.265
   1.137   (  18.977   18.977    4.508)   27.214
   1.263   (  12.171   12.171   -2.082)   17.337
   1.786   (  -2.038   -2.038    1.025)    3.059
   1.965   (   0.688    0.688    0.324)    1.026
   2.327   (   0.769    0.769    1.313)    1.705
   2.376   (   1.906    1.906    0.124)    2.698
   2.638   (  17.082   17.082   -4.227)   24.525
   3.346   (   0.832    0.832   -1.095)    1.608
   3.363   (   0.960    0.960    2.683)    3.007
   3.419   (   1.420    1.420    1.115)    2.297
   3.693   (   4.092    4.092   -5.881)    8.251
   3.870   (  -0.062   -0.062    3.804)    3.805
   4.698   (  -2.573   -2.573   -5.074)    6.244
   5.180   (  -2.398   -2.398   -0.006)    3.391
   5.226   (   2.738    2.738    0.636)    3.924
   6.126   (  -3.014   -3.014    0.603)    4.305
   6.421   (   0.239    0.239    0.323)    0.468
   6.523   (   8.087    8.087   -0.199)   11.438
   8.013   (  33.548   33.548  -46.002)   66.084
  11.370   (  -2.817   -2.817    2.439)    4.672
  11.454   ( -16.875  -16.875   23.803)   33.706
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.685   (  21.169  -14.776    5.742)   26.447
   0.889   (   8.784    1.797   18.809)   20.837
   1.130   (  10.170   -3.030   -4.030)   11.352
   1.444   (   9.733    9.621   -0.640)   13.700
   1.531   (  13.016   10.411    0.110)   16.668
   1.756   (   0.296    1.192    2.813)    3.070
   2.132   (  14.228  -13.364   -0.212)   19.521
   2.302   (  -0.689   -0.734   -1.959)    2.203
   2.430   (   2.930    1.004   -0.500)    3.138
   3.051   (  18.164   18.894   -6.066)   26.901
   3.332   (  -1.118    2.614   -0.252)    2.854
   3.442   (   3.673    1.210    3.844)    5.452
   3.453   (   5.109    3.643    1.636)    6.484
   3.769   (   2.732   -2.474    0.024)    3.686
   3.933   (   6.906    5.568   -1.411)    8.982
   4.623   (  -4.621   -2.674   -5.141)    7.411
   5.011   ( -11.530   -4.698    0.013)   12.451
   5.335   (   4.064   -1.547    0.599)    4.389
   6.112   (   0.992   -7.292    0.666)    7.389
   6.402   (  -1.882    0.913    0.314)    2.115
   6.821   (  20.649   12.849   -0.254)   24.321
   8.537   (  16.726   17.430  -26.803)   36.083
  11.143   ( -12.304   -6.148   17.485)   22.247
  11.284   (  -5.470   -1.979    2.154)    6.203
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.013   (  10.615   -0.616    9.904)   14.531
   1.154   (  15.528  -10.125    5.839)   19.435
   1.357   (  10.683   -3.785   -2.081)   11.523
   1.526   (  -0.705    1.080   -0.040)    1.291
   1.767   (   9.616    3.229    0.510)   10.157
   1.838   (   6.702    4.914    2.445)    8.662
   2.293   (   0.600   -0.067   -3.793)    3.841
   2.391   (   2.856   -3.237   -2.075)    4.790
   2.611   (  19.832  -16.437    0.012)   25.759
   3.282   (   4.948    3.298   -2.089)    6.303
   3.309   (  -0.549    3.019   -0.949)    3.212
   3.498   (   3.309   -0.235    3.548)    4.858
   3.643   (   8.111    4.729    1.743)    9.550
   3.832   (   4.843   -2.882   -1.090)    5.740
   4.172   (  16.318   15.891   -0.957)   22.797
   4.507   (  -6.611   -3.132   -5.085)    8.909
   4.719   ( -16.451  -12.525   -0.132)   20.677
   5.355   (  -2.372    3.813    0.442)    4.512
   6.147   (   2.064   -7.363    0.445)    7.660
   6.350   (  -2.992   -0.460    0.189)    3.033
   7.344   (  28.611   15.754   -0.278)   32.663
   8.740   (   3.762   14.758  -12.409)   19.645
  10.945   (  -7.124   -2.617   10.893)   13.276
  11.146   (  -7.717   -1.177    2.138)    8.093
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.203   (   8.184    0.638    7.890)   11.386
   1.410   (   5.904   -5.704    2.359)    8.541
   1.487   (   0.265   -1.352    1.635)    2.138
   1.534   (   5.610   -1.584   -0.494)    5.850
   1.930   (   6.185   -3.788    0.277)    7.258
   1.965   (   4.633    0.004    1.220)    4.791
   2.304   (   0.558    0.031   -3.462)    3.507
   2.378   (  -2.341    1.790   -2.848)    4.098
   3.047   (  17.233  -12.599   -0.140)   21.348
   3.275   (  -1.553    2.263   -0.403)    2.774
   3.410   (   9.481   -7.886   -0.538)   12.344
   3.570   (   3.460   -1.948    4.774)    6.210
   3.711   (   0.537   -3.611    3.227)    4.872
   3.992   (   9.294   -2.164   -5.501)   11.014
   4.285   ( -17.164  -14.291   -2.289)   22.452
   4.386   ( -11.187   -8.098   -3.634)   14.281
   4.581   (  19.472   20.923    0.448)   28.585
   5.232   (  -9.607   11.959    0.366)   15.344
   6.190   (   2.113   -6.296    0.244)    6.646
   6.287   (  -2.860   -2.942    0.112)    4.105
   7.941   (  27.749   16.451   -0.266)   32.260
   8.706   (  -7.211   15.221   -4.564)   17.450
  10.846   (  -2.264   -1.378    6.526)    7.044
  10.981   (  -7.664   -0.818    2.541)    8.115
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.351   (  -6.000   -0.238    5.171)    7.924
   1.351   (   6.000   -0.238    5.171)    7.924
   1.572   (  -2.141   -1.911    0.027)    2.870
   1.572   (   2.141   -1.911    0.027)    2.870
   2.008   (  -0.853   -5.530    0.606)    5.628
   2.008   (   0.853   -5.530    0.606)    5.628
   2.327   (  -1.975    1.750   -3.039)    4.025
   2.327   (   1.975    1.750   -3.039)    4.025
   3.230   (  -3.446    1.268   -0.377)    3.691
   3.230   (   3.446    1.268   -0.377)    3.691
   3.607   (  -2.482   -3.831    3.753)    5.909
   3.607   (   2.482   -3.831    3.753)    5.909
   3.862   ( -19.168  -21.193    0.130)   28.576
   3.862   (  19.168  -21.193    0.130)   28.576
   4.189   (  -9.366   -3.567   -5.708)   11.533
   4.189   (   9.366   -3.567   -5.708)   11.533
   4.962   ( -16.302   18.772    0.402)   24.865
   4.962   (  16.302   18.772    0.402)   24.865
   6.235   (  -2.262   -4.990    0.132)    5.481
   6.235   (   2.262   -4.990    0.132)    5.481
   8.435   ( -19.023   16.104   -1.084)   24.947
   8.435   (  19.023   16.104   -1.084)   24.947
  10.859   (  -3.536   -0.913    3.818)    5.283
  10.859   (   3.536   -0.913    3.818)    5.283
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.553   (   6.721    6.721    8.364)   12.661
   0.979   (   5.903    5.903   15.992)   18.040
   1.142   (   4.398    4.398   -5.217)    8.118
   1.545   (   0.782    0.782   -1.256)    1.674
   1.718   (   6.569    6.569    1.210)    9.369
   1.836   (   6.057    6.057    1.948)    8.785
   1.981   (  -0.196   -0.196    0.409)    0.494
   2.294   (  -0.127   -0.127   -3.771)    3.775
   2.468   (   2.257    2.257    0.073)    3.192
   3.280   (   4.874    4.874   -4.195)    8.069
   3.380   (   0.654    0.654   -1.670)    1.908
   3.488   (   1.709    1.709    1.512)    2.850
   3.572   (   7.637    7.637    1.573)   10.914
   3.718   (  -1.876   -1.876    1.324)    2.966
   4.174   (  17.872   17.872   -0.699)   25.284
   4.543   (  -4.941   -4.941   -5.145)    8.678
   4.764   ( -17.646  -17.646    0.019)   24.956
   5.390   (   5.060    5.060    0.968)    7.220
   6.016   (  -1.607   -1.607    1.119)    2.533
   6.394   (  -1.684   -1.684    0.418)    2.418
   7.211   (  24.024   24.024   -0.357)   33.977
   8.875   (  13.462   13.462  -17.171)   25.638
  11.001   (  -6.896   -6.896   13.025)   16.271
  11.221   (  -4.235   -4.235    2.026)    6.323
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.864   (  18.474   -8.460    5.063)   20.940
   1.117   (   7.236    3.253   11.881)   14.286
   1.295   (   9.526   -1.541   -4.057)   10.468
   1.521   (  -2.291   -1.146    0.921)    2.723
   1.806   (   2.118    0.312    0.463)    2.190
   1.934   (   2.819    3.192    1.210)    4.427
   2.134   (  13.836  -14.772   -0.184)   20.240
   2.300   (   0.777    0.738   -4.156)    4.292
   2.509   (   1.729    1.496   -0.035)    2.287
   3.340   (   1.091    2.824   -0.389)    3.052
   3.378   (  -0.218    2.324   -2.448)    3.383
   3.506   (   0.820    1.511    2.960)    3.423
   3.656   (  -0.412   -2.333    2.769)    3.644
   3.755   (   6.358   -3.579   -4.043)    8.341
   4.330   ( -22.611  -21.268   -0.304)   31.043
   4.415   (  -6.798   -6.565   -5.219)   10.796
   4.616   (  21.335   22.449    0.370)   30.972
   5.481   (   3.034    7.442    0.800)    8.077
   6.028   (   2.129   -3.119    0.762)    3.853
   6.334   (  -4.011   -1.090    0.281)    4.166
   7.774   (  29.369   24.805   -0.247)   38.443
   9.066   (   4.924   14.998   -9.097)   18.219
  10.875   (  -4.967   -3.801    8.874)   10.857
  11.106   (  -6.778   -2.806    2.101)    7.630
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.194   (  11.384   -4.317    3.654)   12.712
   1.279   (   8.911   -3.210    6.926)   11.734
   1.462   (  -1.734   -1.321    1.803)    2.829
   1.490   (   6.337   -2.519   -1.042)    6.899
   1.814   (  -0.498   -7.328    0.033)    7.344
   1.939   (  -2.384   -2.761    0.810)    3.737
   2.331   (   2.305    2.691   -3.228)    4.793
   2.399   (   4.221   -3.053   -1.095)    5.323
   2.649   (  18.299  -18.315   -0.097)   25.890
   3.323   (  -2.506    1.512   -0.827)    3.042
   3.383   (   0.545    2.027   -1.207)    2.421
   3.534   (   1.397   -1.563    4.006)    4.521
   3.631   (  -1.695   -4.156    3.182)    5.502
   3.829   ( -23.873  -25.011   -0.086)   34.576
   3.910   (   7.888   -5.391   -5.445)   10.997
   4.259   (  -8.132   -6.500   -5.642)   11.841
   5.039   (  18.443   21.821    0.399)   28.574
   5.491   (  -2.704   12.312    0.589)   12.619
   6.083   (   2.900   -3.183    0.295)    4.316
   6.238   (  -4.887   -1.457    0.121)    5.101
   8.370   (  27.053   23.710   -0.267)   35.973
   9.069   (  -5.098   18.158   -3.843)   19.248
  10.806   (  -1.401   -2.396    5.709)    6.348
  10.953   (  -7.440   -1.941    2.447)    8.069
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.341   (  -4.286   -0.853    4.583)    6.333
   1.341   (   4.286   -0.853    4.583)    6.333
   1.508   (  -5.371   -4.225    0.723)    6.872
   1.508   (   5.371   -4.225    0.723)    6.872
   1.854   (  -4.626   -8.639    0.214)    9.802
   1.854   (   4.626   -8.639    0.214)    9.802
   2.386   (  -2.622    3.654   -1.974)    4.912
   2.386   (   2.622    3.654   -1.974)    4.912
   3.122   ( -21.239  -19.113   -0.387)   28.575
   3.122   (  21.239  -19.113   -0.387)   28.575
   3.426   (  -5.937   -5.254   -0.178)    7.930
   3.426   (   5.937   -5.254   -0.178)    7.930
   3.601   (  -1.585   -4.337    3.268)    5.657
   3.601   (   1.585   -4.337    3.268)    5.657
   4.085   (  -8.638   -6.232   -5.708)   12.084
   4.085   (   8.638   -6.232   -5.708)   12.084
   5.350   ( -11.165   17.934    0.457)   21.130
   5.350   (  11.165   17.934    0.457)   21.130
   6.148   (  -3.669   -2.635    0.078)    4.518
   6.148   (   3.669   -2.635    0.078)    4.518
   8.839   ( -17.337   21.381   -1.094)   27.548
   8.839   (  17.337   21.381   -1.094)   27.548
  10.830   (  -3.829   -1.884    3.513)    5.527
  10.830   (   3.829   -1.884    3.513)    5.527
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.830   (   6.144    6.144    6.167)   10.655
   1.201   (   4.449    4.449    9.261)   11.196
   1.312   (   3.498    3.498   -4.615)    6.765
   1.488   (  -1.995   -1.995    1.820)    3.357
   1.767   (  -4.381   -4.381    0.848)    6.253
   1.940   (  -1.653   -1.653    0.190)    2.345
   1.954   (  -2.183   -2.183   -0.019)    3.087
   2.337   (   2.851    2.851   -3.619)    5.418
   2.547   (   1.462    1.462    0.769)    2.205
   3.382   (  -1.589   -1.589   -1.870)    2.923
   3.398   (   2.025    2.025    0.088)    2.865
   3.550   (   1.192    1.192    1.586)    2.315
   3.584   (  -3.160   -3.160    2.942)    5.350
   3.722   (   0.915    0.915   -4.685)    4.860
   3.818   ( -26.000  -26.000    0.190)   36.770
   4.274   (  -7.502   -7.502   -5.615)   12.003
   5.075   (  20.829   20.829    0.298)   29.458
   5.628   (   6.246    6.246    0.733)    8.864
   6.017   (   1.451    1.451    0.389)    2.089
   6.298   (  -2.629   -2.629    0.201)    3.723
   8.327   (  27.780   27.780   -0.127)   39.287
   9.317   (   8.681    8.681   -5.580)   13.485
  10.799   (  -3.298   -3.298    6.621)    8.099
  11.024   (  -5.272   -5.272    2.184)    7.769
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.056   (  13.229   -5.459    4.058)   14.876
   1.288   (   3.604    1.573    6.734)    7.798
   1.425   (   4.022   -1.951   -1.769)    4.808
   1.446   (   0.035   -1.558    1.754)    2.347
   1.639   (  -6.557   -9.245    1.076)   11.386
   1.848   (  -7.154   -6.023   -0.269)    9.356
   2.090   (  15.190  -17.849   -0.173)   23.438
   2.408   (   3.678    4.492   -2.289)    6.241
   2.556   (  -0.530    1.956    0.447)    2.076
   3.139   ( -25.423  -25.358   -0.072)   35.908
   3.422   (   0.073    1.336   -0.778)    1.547
   3.452   (  -4.220   -2.932   -1.365)    5.317
   3.533   (  -0.382   -0.834    2.818)    2.964
   3.568   (  -0.809   -0.934    2.117)    2.452
   3.798   (   5.879   -5.145   -5.006)    9.279
   4.110   (  -8.191   -7.789   -5.488)   12.565
   5.452   (  15.276   17.852    0.305)   23.498
   5.720   (   2.050    9.577    0.520)    9.808
   6.063   (   2.667    0.833   -0.104)    2.796
   6.222   (  -4.582   -0.553   -0.002)    4.615
   8.874   (  23.945   24.099   -0.190)   33.973
   9.405   (  -0.625   13.675   -2.698)   13.953
  10.753   (  -0.823   -2.628    4.607)    5.367
  10.896   (  -6.631   -3.734    2.405)    7.981
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.290   (  -6.835   -4.962    3.238)    9.046
   1.290   (   6.835   -4.962    3.238)    9.046
   1.425   (  -3.458   -2.619    2.851)    5.191
   1.425   (   3.458   -2.619    2.851)    5.191
   1.675   (  -8.499   -8.260   -0.642)   11.869
   1.675   (   8.499   -8.260   -0.642)   11.869
   2.443   (  -6.065   -2.983   -0.648)    6.790
   2.443   (   6.065   -2.983   -0.648)    6.790
   2.612   ( -17.352  -20.442   -0.117)   26.813
   2.612   (  17.352  -20.442   -0.117)   26.813
   3.421   (  -0.318    0.957   -1.359)    1.692
   3.421   (   0.318    0.957   -1.359)    1.692
   3.547   (  -1.125   -0.951    2.540)    2.936
   3.547   (   1.125   -0.951    2.540)    2.936
   3.943   (  -7.787   -7.218   -5.240)   11.840
   3.943   (   7.787   -7.218   -5.240)   11.840
   5.682   (  -6.431   13.942    0.403)   15.359
   5.682   (   6.431   13.942    0.403)   15.359
   6.128   (  -3.929    0.286   -0.230)    3.946
   6.128   (   3.929    0.286   -0.230)    3.946
   9.267   ( -13.164   19.143   -0.870)   23.249
   9.267   (  13.164   19.143   -0.870)   23.249
  10.781   (  -3.747   -2.844    3.112)    5.640
  10.781   (   3.747   -2.844    3.112)    5.640
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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|
   1.033   (   3.461    3.461    4.068)    6.365
   1.298   (  -2.285   -2.285    4.413)    5.469
   1.414   (  -1.320   -1.320    3.053)    3.578
   1.420   (   1.740    1.740   -2.824)    3.746
   1.464   (  -5.814   -5.814    3.339)    8.874
   1.697   (  -8.455   -8.455   -0.803)   11.985
   1.868   (  -1.524   -1.524    0.115)    2.158
   2.503   (   4.188    4.188   -1.107)    6.025
   2.553   ( -27.965  -27.965    0.301)   39.549
   2.590   (   0.579    0.579    0.783)    1.133
   3.431   (  -0.412   -0.412   -1.521)    1.628
   3.432   (  -0.308   -0.308   -1.002)    1.092
   3.528   (   0.978    0.978    1.938)    2.381
   3.593   (   0.718    0.718    1.054)    1.464
   3.726   (  -0.442   -0.442   -5.027)    5.066
   3.946   (  -7.810   -7.810   -5.108)   12.169
   5.763   (  11.813   11.813    0.301)   16.708
   5.877   (   5.304    5.304    0.300)    7.507
   6.094   (   1.721    1.721   -0.545)    2.495
   6.197   (  -2.198   -2.198   -0.211)    3.116
   9.325   (  18.506   18.506   -0.179)   26.172
   9.608   (   5.289    5.289   -1.445)    7.618
  10.709   (  -1.242   -1.242    3.513)    3.927
  10.803   (  -4.910   -4.910    2.415)    7.352
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.166   (  -8.153   -5.476    2.904)   10.242
   1.166   (   8.153   -5.476    2.904)   10.242
   1.405   (  -0.930    0.133    3.602)    3.723
   1.405   (   0.930    0.133    3.602)    3.723
   1.527   (  -7.189   -5.723   -0.844)    9.227
   1.527   (   7.189   -5.723   -0.844)    9.227
   2.052   ( -18.415  -19.326    0.136)   26.695
   2.052   (  18.415  -19.326    0.136)   26.695
   2.574   (  -2.269    2.648    0.115)    3.489
   2.574   (   2.269    2.648    0.115)    3.489
   3.427   (  -0.085   -0.164   -1.275)    1.288
   3.427   (   0.085   -0.164   -1.275)    1.288
   3.565   (  -2.254    2.326    0.847)    3.348
   3.565   (   2.254    2.326    0.847)    3.348
   3.799   (  -6.135   -6.398   -4.743)   10.053
   3.799   (   6.135   -6.398   -4.743)   10.053
   5.916   (  -2.417    8.468    0.300)    8.811
   5.916   (   2.417    8.468    0.300)    8.811
   6.137   (  -2.863    0.107   -0.564)    2.920
   6.137   (   2.863    0.107   -0.564)    2.920
   9.596   (  -6.970   11.925   -0.517)   13.822
   9.596   (   6.970   11.925   -0.517)   13.822
  10.719   (  -2.506   -2.819    2.711)    4.645
  10.719   (   2.506   -2.819    2.711)    4.645
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.105   (   0.000    0.000    2.877)    2.877
   1.105   (   0.000    0.000    2.877)    2.877
   1.415   (   0.000    0.000    4.448)    4.448
   1.415   (   0.000    0.000    4.448)    4.448
   1.459   (   0.000    0.000   -1.493)    1.493
   1.459   (   0.000    0.000   -1.493)    1.493
   1.834   (   0.000    0.000    0.121)    0.121
   1.834   (   0.000    0.000    0.121)    0.121
   2.604   (   0.000    0.000    0.486)    0.486
   2.604   (   0.000    0.000    0.486)    0.486
   3.425   (   0.000    0.000   -1.199)    1.199
   3.425   (   0.000    0.000   -1.199)    1.199
   3.609   (   0.000    0.000    0.632)    0.632
   3.609   (   0.000    0.000    0.632)    0.632
   3.712   (   0.000    0.000   -5.082)    5.082
   3.712   (   0.000    0.000   -5.082)    5.082
   6.013   (  -0.000   -0.000    0.377)    0.377
   6.013   (   0.000    0.000    0.377)    0.377
   6.132   (   0.000    0.000   -0.828)    0.828
   6.132   (   0.000    0.000   -0.828)    0.828
   9.725   (   0.000    0.000   -0.327)    0.327
   9.725   (   0.000    0.000   -0.327)    0.327
  10.685   (   0.000    0.000    2.517)    2.517
  10.685   (   0.000    0.000    2.517)    2.517
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.486   (   0.000    0.000    8.011)    8.011
   0.486   (   0.000    0.000    8.011)    8.011
   0.806   (   0.000    0.000   -5.580)    5.580
   0.806   (  -0.000    0.000   -5.580)    5.580
   1.363   (   0.000    0.000   26.455)   26.455
   1.773   (   0.000    0.000  -10.304)   10.304
   1.949   (   0.000    0.000    0.070)    0.070
   1.949   (  -0.000    0.000    0.070)    0.070
   2.346   (  -0.000    0.000    0.388)    0.388
   2.622   (   0.000    0.000  -10.119)   10.119
   3.312   (   0.000    0.000   -0.365)    0.365
   3.312   (  -0.000    0.000   -0.365)    0.365
   3.401   (  -0.000    0.000    0.367)    0.367
   3.401   (   0.000    0.000    0.367)    0.367
   4.183   (   0.000    0.000    6.759)    6.759
   4.641   (   0.000    0.000   -5.564)    5.564
   5.170   (   0.000    0.000    0.097)    0.097
   5.170   (  -0.000    0.000    0.097)    0.097
   6.200   (   0.000    0.000    0.021)    0.021
   6.200   (  -0.000   -0.000    0.021)    0.021
   6.409   (   0.000    0.000    0.025)    0.025
   6.409   (   0.000    0.000    0.025)    0.025
  11.491   (   0.000    0.000    1.898)    1.898
  12.154   (   0.000    0.000   -0.819)    0.819
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.730   (  17.847    0.000    5.885)   18.792
   0.786   (  20.040    0.000    7.406)   21.364
   0.909   (   9.257    0.000   -5.169)   10.603
   0.951   (  13.731    0.000   -6.196)   15.064
   1.374   (   1.305    0.000   23.340)   23.376
   1.745   (  -1.479    0.000   -8.934)    9.055
   2.094   (  13.590    0.000    0.057)   13.590
   2.119   (  14.598    0.000    0.054)   14.598
   2.364   (   2.189    0.000    0.186)    2.197
   2.459   ( -10.490    0.000   -5.219)   11.716
   3.307   (  -0.451    0.000   -0.369)    0.583
   3.404   (   0.379    0.000    0.370)    0.530
   3.414   (   4.016    0.000    2.209)    4.583
   3.486   (   9.843    0.000   -4.772)   10.938
   4.044   (  -7.341    0.000    5.500)    9.173
   4.610   (  -2.504    0.000   -4.311)    4.986
   5.135   (  -3.390    0.000    0.062)    3.390
   5.263   (   7.511    0.000    0.075)    7.512
   6.206   (   0.565    0.000    0.023)    0.565
   6.402   (  -0.739    0.000    0.018)    0.739
   6.464   (   5.907    0.000    0.120)    5.908
   6.898   (  50.665    0.000  -10.457)   51.733
  11.452   (  -3.731    0.000    1.671)    4.088
  11.919   ( -19.361    0.000    3.597)   19.692
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.075   (  14.463    0.000    3.003)   14.771
   1.112   (  11.878    0.000    8.033)   14.339
   1.150   (  13.051    0.000   -2.908)   13.371
   1.296   (  13.511    0.000   -3.717)   14.013
   1.458   (   9.819    0.000   12.215)   15.673
   1.751   (   2.387    0.000   -3.153)    3.955
   2.255   (  -4.147    0.000   -0.784)    4.220
   2.347   (   0.152    0.000   -3.226)    3.230
   2.472   (  21.537    0.000    0.078)   21.538
   2.621   (  23.861    0.000   -1.645)   23.918
   3.294   (  -0.888    0.000   -0.363)    0.960
   3.424   (   1.882    0.000    0.367)    1.918
   3.519   (   6.542    0.000    4.019)    7.678
   3.678   (   7.742    0.000   -5.647)    9.583
   3.968   (  -0.397    0.000    3.650)    3.672
   4.539   (  -4.469    0.000   -3.906)    5.935
   5.027   (  -7.089    0.000   -0.020)    7.089
   5.388   (   2.536    0.000    0.097)    2.538
   6.223   (   1.008    0.000    0.028)    1.009
   6.381   (  -1.194    0.000    0.009)    1.194
   6.694   (  17.239    0.000    0.053)   17.239
   7.822   (  34.869    0.000  -11.850)   36.828
  11.343   (  -6.694    0.000    1.310)    6.821
  11.506   ( -18.292    0.000    6.638)   19.459
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.289   (   5.839    0.000   11.351)   12.765
   1.333   (  10.762    0.000    0.469)   10.772
   1.397   (  10.054    0.000   -0.443)   10.064
   1.452   (   4.734    0.000   -4.077)    6.248
   1.759   (  14.772    0.000    3.108)   15.096
   1.839   (   5.751    0.000    1.428)    5.926
   2.204   (  -0.194    0.000   -3.821)    3.825
   2.333   (  -1.699    0.000   -2.393)    2.935
   2.901   (  17.619    0.000    0.065)   17.619
   3.155   (  25.521    0.000   -1.249)   25.551
   3.271   (  -1.255    0.000   -0.311)    1.293
   3.515   (   8.436    0.000    0.335)    8.443
   3.659   (   5.848    0.000    6.775)    8.950
   3.805   (   4.716    0.000   -6.171)    7.767
   4.024   (   6.719    0.000    2.403)    7.136
   4.424   (  -6.692    0.000   -4.225)    7.914
   4.839   ( -11.178    0.000   -0.090)   11.178
   5.329   (  -8.067    0.000    0.062)    8.068
   6.246   (   1.258    0.000    0.031)    1.259
   6.355   (  -1.323    0.000    0.010)    1.323
   7.168   (  27.373    0.000   -0.121)   27.373
   8.320   (  14.286    0.000   -6.645)   15.756
  11.172   ( -10.072    0.000    4.544)   11.049
  11.209   ( -10.110    0.000    2.416)   10.394
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.388   (   4.177    0.000    7.397)    8.495
   1.502   (  -0.032    0.000   -0.080)    0.086
   1.513   (   6.447    0.000    0.015)    6.447
   1.552   (   4.812    0.000    0.018)    4.812
   1.969   (   6.372    0.000    0.443)    6.388
   2.006   (   8.517    0.000    1.584)    8.663
   2.227   (   1.473    0.000   -3.288)    3.603
   2.280   (  -3.072    0.000   -2.526)    3.977
   3.135   (   5.807    0.000    0.000)    5.807
   3.242   (  -1.524    0.000   -0.211)    1.538
   3.600   (  15.673    0.000    1.686)   15.764
   3.738   (   2.182    0.000    5.400)    5.824
   3.808   (  18.631    0.000    0.206)   18.632
   3.879   (   3.842    0.000   -5.023)    6.323
   4.234   (  -1.418    0.000   -0.910)    1.685
   4.309   (   9.935    0.000   -3.129)   10.416
   4.565   ( -15.397    0.000   -0.079)   15.398
   5.076   ( -15.943    0.000    0.063)   15.943
   6.273   (   1.326    0.000    0.032)    1.327
   6.327   (  -1.306    0.000    0.020)    1.306
   7.753   (  27.528    0.000   -0.116)   27.528
   8.441   (  -2.243    0.000   -2.552)    3.398
  10.996   (  -5.878    0.000    4.597)    7.462
  11.035   (  -6.961    0.000    0.827)    7.010
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.465   (  -3.104    0.000    3.627)    4.774
   1.465   (   3.104    0.000    3.627)    4.774
   1.592   (  -0.979    0.000    0.020)    0.979
   1.592   (   0.979    0.000    0.020)    0.979
   2.073   (  -2.382    0.000    0.277)    2.398
   2.073   (   2.382    0.000    0.277)    2.398
   2.241   (  -0.290    0.000   -2.428)    2.445
   2.241   (   0.290    0.000   -2.428)    2.445
   3.206   (  -2.143    0.000   -0.095)    2.145
   3.206   (   2.143    0.000   -0.095)    2.145
   3.752   (  -1.425    0.000    4.889)    5.093
   3.752   (   1.425    0.000    4.889)    5.093
   4.043   ( -10.713    0.000   -5.897)   12.229
   4.043   (  10.713    0.000   -5.897)   12.229
   4.208   ( -18.914    0.000    0.033)   18.914
   4.208   (  18.914    0.000    0.033)   18.914
   4.696   ( -20.088    0.000    0.082)   20.088
   4.696   (  20.088    0.000    0.082)   20.088
   6.300   (  -1.306    0.000    0.028)    1.307
   6.300   (   1.306    0.000    0.028)    1.307
   8.234   ( -17.540    0.000   -0.544)   17.549
   8.234   (  17.540    0.000   -0.544)   17.549
  10.945   (  -1.124    0.000    2.350)    2.604
  10.945   (   1.124    0.000    2.350)    2.604
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.556   (   3.434    3.434    8.805)   10.055
   0.852   (   2.292    2.292   -6.526)    7.287
   1.077   (  11.875   11.875    5.565)   17.691
   1.207   (  14.850   14.850   -2.671)   21.170
   1.403   (   3.618    3.618   18.348)   19.048
   1.741   (   0.186    0.186   -6.954)    6.959
   1.971   (   0.817    0.817    0.189)    1.170
   2.341   (  -3.860   -3.860    0.080)    5.459
   2.379   (   1.602    1.602    0.077)    2.268
   2.512   (  18.643   18.643   -5.087)   26.851
   3.329   (   0.659    0.659   -0.453)    1.037
   3.430   (   3.271    3.271    3.391)    5.735
   3.437   (   1.572    1.572    0.465)    2.271
   3.546   (   5.762    5.762   -5.677)    9.931
   3.991   (  -3.179   -3.179    4.575)    6.414
   4.584   (  -2.434   -2.434   -4.039)    5.307
   5.180   (  -2.468   -2.468   -0.003)    3.491
   5.236   (   3.061    3.061    0.251)    4.336
   6.138   (  -2.581   -2.581    0.382)    3.670
   6.429   (   0.410    0.410    0.255)    0.633
   6.518   (   7.895    7.895   -0.133)   11.166
   7.398   (  36.952   36.952  -12.939)   53.835
  11.418   (  -3.241   -3.241    1.501)    4.823
  11.744   ( -14.739  -14.739    5.506)   21.558
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.837   (  18.018  -11.427    7.679)   22.676
   1.010   (  12.130   -5.709   -6.745)   15.007
   1.258   (   5.811    4.620   11.954)   14.071
   1.414   (   5.940    5.573   -2.337)    8.474
   1.579   (  11.656    9.365    4.535)   15.624
   1.780   (   3.882    3.020   -1.371)    5.106
   2.131   (  13.303  -12.303    0.109)   18.120
   2.235   (  -3.707   -2.753   -3.313)    5.683
   2.417   (   2.303    1.891   -0.473)    3.018
   2.944   (  19.483   16.036   -2.971)   25.408
   3.322   (  -1.008    2.167   -0.392)    2.422
   3.467   (   2.156    2.864    0.458)    3.614
   3.556   (   7.695    3.606    5.628)   10.193
   3.685   (   6.814    1.603   -6.025)    9.236
   3.977   (   3.020    2.395    2.997)    4.882
   4.510   (  -4.707   -2.718   -3.917)    6.700
   5.011   ( -11.607   -4.764    0.003)   12.546
   5.346   (   3.977   -0.924    0.303)    4.094
   6.126   (   0.832   -6.432    0.430)    6.500
   6.408   (  -2.002    1.248    0.215)    2.369
   6.816   (  20.610   12.502   -0.162)   24.106
   8.105   (  28.091   22.677  -10.961)   37.729
  11.321   (  -6.100   -2.137    1.342)    6.602
  11.417   ( -15.052   -7.836    6.476)   18.164
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.169   (  13.570   -8.492    4.918)   16.747
   1.289   (  12.058   -6.263   -2.883)   13.890
   1.358   (   5.593    0.036    8.836)   10.458
   1.502   (   2.530    2.308   -2.619)    4.311
   1.797   (   9.043    2.297    1.975)    9.536
   1.886   (   5.610    3.962    1.402)    7.010
   2.194   (   0.116   -0.754   -4.295)    4.363
   2.344   (   1.735   -0.862   -1.521)    2.463
   2.611   (  20.081  -16.845    0.016)   26.211
   3.244   (   9.139    4.247   -1.064)   10.134
   3.296   (  -1.374    2.126   -0.347)    2.555
   3.557   (   7.058    1.059    1.529)    7.299
   3.691   (   5.028    0.473    2.528)    5.648
   3.790   (   2.701   -1.752   -2.456)    4.049
   4.179   (  16.355   16.184    1.292)   23.045
   4.390   (  -6.771   -3.018   -4.407)    8.624
   4.717   ( -16.651  -12.518   -0.057)   20.831
   5.363   (  -2.493    4.278    0.235)    4.957
   6.156   (   1.831   -6.789    0.287)    7.037
   6.353   (  -3.115   -0.258    0.119)    3.128
   7.338   (  28.642   15.677   -0.178)   32.653
   8.519   (  12.243   16.798   -6.136)   21.673
  11.141   ( -10.226   -3.398    5.839)   12.257
  11.185   (  -8.107   -1.617    0.820)    8.307
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.374   (   5.820   -1.815    6.537)    8.938
   1.446   (   6.329   -3.581   -0.083)    7.273
   1.520   (   1.863   -1.201    0.835)    2.369
   1.524   (   4.428   -1.508   -0.074)    4.678
   1.950   (   5.778   -4.675    1.354)    7.555
   1.988   (   3.862    0.512    0.647)    3.949
   2.218   (   1.687   -0.195   -3.487)    3.879
   2.319   (  -2.720    3.378   -1.883)    4.728
   3.045   (  16.961  -12.667   -0.082)   21.169
   3.266   (  -1.527    2.026   -0.290)    2.554
   3.408   (  10.896   -8.546    0.106)   13.848
   3.691   (   3.752   -2.690    6.248)    7.769
   3.787   (   2.019   -2.633    2.723)    4.292
   3.850   (   7.295   -2.947   -6.902)   10.466
   4.200   ( -11.586   -6.295   -3.816)   13.727
   4.344   ( -18.018  -16.495   -0.605)   24.436
   4.586   (  20.041   21.531    0.114)   29.415
   5.239   (  -9.638   12.245    0.183)   15.584
   6.195   (   1.955   -6.046    0.156)    6.356
   6.290   (  -2.886   -2.855    0.067)    4.060
   7.936   (  27.633   16.393   -0.150)   32.130
   8.620   (  -2.468   15.571   -2.483)   15.960
  10.977   (  -5.325   -1.924    4.392)    7.165
  11.025   (  -6.860   -0.880    0.847)    6.968
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.460   (  -2.998   -0.477    3.705)    4.790
   1.460   (   2.998   -0.477    3.705)    4.790
   1.569   (  -1.751   -2.091   -0.344)    2.749
   1.569   (   1.751   -2.091   -0.344)    2.749
   2.024   (  -0.913   -5.309    0.697)    5.432
   2.024   (   0.913   -5.309    0.697)    5.432
   2.259   (  -2.571    2.472   -2.396)    4.297
   2.259   (   2.571    2.472   -2.396)    4.297
   3.223   (  -3.312    0.914   -0.201)    3.441
   3.223   (   3.312    0.914   -0.201)    3.441
   3.699   (  -6.655   -7.089    3.692)   10.401
   3.699   (   6.655   -7.089    3.692)   10.401
   3.863   ( -13.849  -16.598   -0.077)   21.617
   3.863   (  13.849  -16.598   -0.077)   21.617
   4.060   ( -10.945   -5.168   -4.746)   13.001
   4.060   (  10.945   -5.168   -4.746)   13.001
   4.969   ( -16.248   18.884    0.183)   24.913
   4.969   (  16.248   18.884    0.183)   24.913
   6.237   (  -2.196   -4.904    0.083)    5.374
   6.237   (   2.196   -4.904    0.083)    5.374
   8.414   ( -17.329   15.990   -0.606)   23.587
   8.414   (  17.329   15.990   -0.606)   23.587
  10.934   (  -1.338   -1.122    2.283)    2.874
  10.934   (   1.338   -1.122    2.283)    2.874
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.759   (   5.947    5.947    9.740)   12.869
   0.991   (   4.252    4.252   -8.244)   10.204
   1.320   (   2.331    2.331   13.066)   13.476
   1.500   (   2.619    2.619   -3.482)    5.084
   1.748   (   6.073    6.073    1.598)    8.736
   1.873   (   5.350    5.350    0.812)    7.609
   1.988   (  -0.264   -0.264    0.231)    0.439
   2.190   (  -1.290   -1.290   -4.494)    4.850
   2.470   (   2.534    2.534    0.046)    3.583
   3.193   (   8.832    8.832   -2.758)   12.792
   3.353   (   0.398    0.398   -0.715)    0.910
   3.512   (   1.834    1.834    0.616)    2.666
   3.612   (   1.855    1.855    1.582)    3.063
   3.733   (   0.886    0.886   -0.023)    1.253
   4.180   (  17.790   17.790    1.040)   25.180
   4.428   (  -5.079   -5.079   -4.138)    8.289
   4.764   ( -17.556  -17.556    0.012)   24.828
   5.408   (   5.088    5.088    0.539)    7.216
   6.038   (  -1.666   -1.666    0.701)    2.458
   6.403   (  -1.767   -1.767    0.276)    2.513
   7.203   (  24.099   24.099   -0.231)   34.082
   8.577   (  20.441   20.441   -8.119)   30.027
  11.222   (  -9.249   -9.249    6.023)   14.400
  11.260   (  -4.701   -4.701    1.057)    6.732
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.007   (  14.531   -5.676    7.647)   17.373
   1.173   (  12.305   -3.530   -6.922)   14.552
   1.371   (   2.801    1.036    9.070)    9.549
   1.527   (   0.187    0.319   -0.735)    0.823
   1.820   (   1.340   -0.199    0.852)    1.601
   1.962   (   2.479    2.349    0.931)    3.540
   2.131   (  13.085  -14.281   -0.059)   19.369
   2.193   (   1.639    0.992   -4.434)    4.830
   2.508   (   1.392    2.646   -0.039)    2.990
   3.328   (   3.878    4.045   -0.548)    5.631
   3.341   (  -1.037    1.809   -0.831)    2.244
   3.563   (   3.501    0.092    2.115)    4.091
   3.661   (   3.259   -2.414   -2.774)    4.913
   3.727   (  -0.807   -3.860    1.858)    4.360
   4.284   ( -10.437   -7.637   -3.270)   13.340
   4.347   ( -19.828  -20.706   -0.470)   28.672
   4.620   (  21.934   22.959    0.080)   31.752
   5.497   (   2.809    7.583    0.458)    8.100
   6.043   (   1.630   -3.195    0.477)    3.618
   6.340   (  -4.185   -1.108    0.175)    4.332
   7.769   (  29.476   24.922   -0.155)   38.600
   8.900   (  10.447   18.003   -4.757)   21.351
  11.042   (  -7.620   -5.626    5.001)   10.711
  11.144   (  -6.666   -2.704    0.907)    7.251
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.282   (  11.500   -6.323    3.834)   13.672
   1.397   (   6.023   -1.897    0.350)    6.324
   1.460   (   6.538   -2.670    1.644)    7.251
   1.513   (  -1.176   -0.636    1.321)    1.880
   1.821   (  -0.478   -7.377    0.490)    7.409
   1.959   (  -2.846   -3.700    0.697)    4.720
   2.253   (   4.051    3.740   -2.937)    6.247
   2.377   (   3.316   -1.292   -0.739)    3.635
   2.647   (  18.251  -18.348   -0.074)   25.879
   3.305   (  -2.571    0.920   -0.554)    2.786
   3.372   (   1.028    1.886   -0.059)    2.148
   3.629   (   1.652   -3.243    4.666)    5.917
   3.706   (  -1.263   -4.872    2.814)    5.766
   3.772   (   6.492   -4.788   -5.801)    9.936
   3.830   ( -22.742  -24.018   -0.452)   33.080
   4.137   (  -8.595   -6.675   -4.235)   11.678
   5.046   (  18.478   21.740    0.189)   28.533
   5.502   (  -2.873   12.411    0.327)   12.743
   6.088   (   2.546   -3.389    0.187)    4.243
   6.241   (  -5.011   -1.533    0.070)    5.240
   8.365   (  26.824   23.757   -0.156)   35.833
   8.996   (  -1.563   19.014   -2.148)   19.198
  10.920   (  -3.831   -3.435    3.600)    6.279
  10.998   (  -6.777   -1.861    1.166)    7.124
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.441   (  -3.231   -1.691    3.779)    5.251
   1.441   (   3.231   -1.691    3.779)    5.251
   1.512   (  -3.068   -3.126   -0.698)    4.435
   1.512   (   3.068   -3.126   -0.698)    4.435
   1.865   (  -4.975   -9.278    0.587)   10.544
   1.865   (   4.975   -9.278    0.587)   10.544
   2.341   (  -3.962    4.964   -1.549)    6.538
   2.341   (   3.962    4.964   -1.549)    6.538
   3.114   ( -20.552  -18.406   -0.226)   27.590
   3.114   (  20.552  -18.406   -0.226)   27.590
   3.426   (  -7.217   -6.965    0.074)   10.030
   3.426   (   7.217   -6.965    0.074)   10.030
   3.680   (  -1.399   -4.825    3.341)    6.033
   3.680   (   1.399   -4.825    3.341)    6.033
   3.953   (  -9.039   -5.646   -5.005)   11.774
   3.953   (   9.039   -5.646   -5.005)   11.774
   5.358   ( -11.267   17.956    0.237)   21.200
   5.358   (  11.267   17.956    0.237)   21.200
   6.150   (  -3.596   -2.844    0.048)    4.585
   6.150   (   3.596   -2.844    0.048)    4.585
   8.818   ( -15.826   21.469   -0.632)   26.680
   8.818   (  15.826   21.469   -0.632)   26.680
  10.899   (  -2.087   -2.225    2.110)    3.709
  10.899   (   2.087   -2.225    2.110)    3.709
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.993   (   4.779    4.779    8.385)   10.770
   1.178   (   4.527    4.527   -7.537)    9.889
   1.398   (   1.526    1.526    6.681)    7.021
   1.521   (  -0.808   -0.808    0.864)    1.433
   1.781   (  -4.172   -4.172    0.395)    5.913
   1.944   (  -1.742   -1.742    0.102)    2.465
   1.958   (  -2.931   -2.931    0.200)    4.149
   2.250   (   4.321    4.321   -3.260)    6.926
   2.562   (   1.707    1.707    0.488)    2.463
   3.349   (  -1.382   -1.382   -0.910)    2.156
   3.403   (   2.512    2.512    0.275)    3.563
   3.577   (   1.210    1.210    0.720)    1.856
   3.615   (  -1.252   -1.252   -1.818)    2.538
   3.660   (  -2.001   -2.001    1.122)    3.044
   3.822   ( -26.182  -26.182    0.119)   37.027
   4.154   (  -7.523   -7.523   -4.029)   11.376
   5.080   (  20.791   20.791    0.137)   29.403
   5.643   (   6.050    6.050    0.426)    8.566
   6.025   (   0.917    0.917    0.240)    1.319
   6.302   (  -2.770   -2.770    0.117)    3.919
   8.324   (  27.856   27.856   -0.079)   39.395
   9.212   (  11.547   11.547   -3.054)   16.612
  10.927   (  -5.147   -5.147    3.880)    8.248
  11.064   (  -5.099   -5.099    1.119)    7.298
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.160   (   9.916   -4.509    4.614)   11.830
   1.344   (   9.428   -2.390   -3.620)   10.378
   1.430   (   1.199   -0.765    4.026)    4.270
   1.497   (  -1.309   -0.888    1.730)    2.344
   1.660   (  -6.074   -7.863    0.637)    9.956
   1.847   (  -7.459   -6.841    0.129)   10.122
   2.087   (  15.137  -17.525   -0.091)   23.158
   2.355   (   5.395    5.738   -1.923)    8.108
   2.564   (  -1.394    2.777    0.243)    3.117
   3.137   ( -24.403  -24.008   -0.070)   34.233
   3.417   (  -0.695    1.742    0.118)    1.879
   3.431   (  -5.266   -4.962   -0.483)    7.251
   3.586   (   0.306   -1.649    2.011)    2.619
   3.614   (  -0.662   -3.226    1.659)    3.688
   3.690   (   5.337   -3.964   -3.976)    7.746
   3.992   (  -8.034   -7.348   -4.032)   11.610
   5.458   (  15.352   17.809    0.154)   23.513
   5.730   (   1.876    9.376    0.290)    9.566
   6.061   (   2.322    0.344   -0.066)    2.348
   6.222   (  -4.830   -0.708   -0.010)    4.881
   8.870   (  23.690   24.169   -0.114)   33.844
   9.354   (   1.570   14.825   -1.537)   14.986
  10.844   (  -2.490   -3.723    2.814)    5.289
  10.942   (  -6.208   -3.692    1.344)    7.347
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.361   (  -8.290   -6.011    2.629)   10.572
   1.361   (   8.290   -6.011    2.629)   10.572
   1.466   (  -2.679   -1.282    0.396)    2.997
   1.466   (   2.679   -1.282    0.396)    2.997
   1.674   (  -7.915   -8.491    0.472)   11.618
   1.674   (   7.915   -8.491    0.472)   11.618
   2.428   (  -3.822   -0.868   -0.572)    3.961
   2.428   (   3.822   -0.868   -0.572)    3.961
   2.610   ( -18.244  -21.254   -0.104)   28.011
   2.610   (  18.244  -21.254   -0.104)   28.011
   3.403   (  -0.267    0.614   -0.354)    0.757
   3.403   (   0.267    0.614   -0.354)    0.757
   3.598   (  -0.724   -2.840    1.903)    3.494
   3.598   (   0.724   -2.840    1.903)    3.494
   3.829   (  -7.587   -6.237   -4.025)   10.614
   3.829   (   7.587   -6.237   -4.025)   10.614
   5.690   (  -6.521   13.839    0.215)   15.300
   5.690   (   6.521   13.839    0.215)   15.300
   6.124   (  -4.032   -0.059   -0.148)    4.035
   6.124   (   4.032   -0.059   -0.148)    4.035
   9.250   ( -12.040   19.442   -0.512)   22.874
   9.250   (  12.040   19.442   -0.512)   22.874
  10.842   (  -2.626   -3.238    1.880)    4.573
  10.842   (   2.626   -3.238    1.880)    4.573
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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|
   1.130   (   1.908    1.908    3.885)    4.730
   1.347   (   3.348    3.348   -3.158)    5.692
   1.382   (  -5.466   -5.466    2.668)    8.177
   1.479   (  -0.905   -0.905    1.580)    2.033
   1.527   (  -3.472   -3.472    2.113)    5.345
   1.688   (  -8.174   -8.174    0.057)   11.560
   1.870   (  -1.517   -1.517    0.064)    2.146
   2.476   (   5.567    5.567   -0.953)    7.930
   2.559   ( -28.069  -28.069    0.193)   39.697
   2.606   (   0.470    0.470    0.494)    0.828
   3.403   (  -0.269   -0.269   -0.804)    0.890
   3.425   (  -0.832   -0.832    0.133)    1.183
   3.562   (  -0.281   -0.281    1.137)    1.204
   3.611   (   0.463    0.463    0.498)    0.823
   3.622   (  -0.533   -0.533   -3.517)    3.597
   3.840   (  -7.148   -7.148   -3.478)   10.690
   5.768   (  11.863   11.863    0.162)   16.777
   5.882   (   5.132    5.132    0.136)    7.259
   6.083   (   1.423    1.423   -0.338)    2.040
   6.192   (  -2.462   -2.462   -0.135)    3.484
   9.321   (  18.406   18.406   -0.111)   26.031
   9.580   (   6.367    6.367   -0.840)    9.043
  10.778   (  -2.191   -2.191    2.146)    3.769
  10.850   (  -4.784   -4.784    1.429)    6.915
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.230   (  -7.257   -5.559    2.329)    9.434
   1.230   (   7.257   -5.559    2.329)    9.434
   1.441   (  -3.540   -1.344   -0.519)    3.822
   1.441   (   3.540   -1.344   -0.519)    3.822
   1.541   (  -4.721   -3.943    1.761)    6.399
   1.541   (   4.721   -3.943    1.761)    6.399
   2.054   ( -18.594  -19.277    0.083)   26.783
   2.054   (  18.594  -19.277    0.083)   26.783
   2.575   (  -3.411    3.406   -0.002)    4.821
   2.575   (   3.411    3.406   -0.002)    4.821
   3.406   (  -0.626   -0.271   -0.493)    0.842
   3.406   (   0.626   -0.271   -0.493)    0.842
   3.574   (  -1.325    0.509   -0.009)    1.420
   3.574   (   1.325    0.509   -0.009)    1.420
   3.706   (  -5.640   -5.411   -2.740)    8.282
   3.706   (   5.640   -5.411   -2.740)    8.282
   5.921   (  -2.287    8.419    0.144)    8.725
   5.921   (   2.287    8.419    0.144)    8.725
   6.126   (  -3.172   -0.186   -0.348)    3.196
   6.126   (   3.172   -0.186   -0.348)    3.196
   9.586   (  -6.379   12.234   -0.310)   13.801
   9.586   (   6.379   12.234   -0.310)   13.801
  10.773   (  -1.961   -3.138    1.651)    4.052
  10.773   (   1.961   -3.138    1.651)    4.052
======================= 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 5.83e-04 5.83e-04 5.83e-04 5.83e-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.168   (   0.000    0.000    2.180)    2.180
   1.168   (   0.000    0.000    2.180)    2.180
   1.423   (   0.000    0.000   -1.350)    1.350
   1.423   (   0.000    0.000   -1.350)    1.350
   1.505   (   0.000    0.000    2.828)    2.828
   1.505   (   0.000    0.000    2.828)    2.828
   1.836   (   0.000    0.000    0.075)    0.075
   1.836   (   0.000    0.000    0.075)    0.075
   2.614   (   0.000    0.000    0.313)    0.313
   2.614   (   0.000    0.000    0.313)    0.313
   3.402   (   0.000    0.000   -0.682)    0.682
   3.402   (   0.000    0.000   -0.682)    0.682
   3.609   (  -0.000   -0.000   -3.244)    3.244
   3.609   (   0.000    0.000   -3.244)    3.244
   3.620   (   0.000    0.000    0.318)    0.318
   3.620   (   0.000    0.000    0.318)    0.318
   6.019   (  -0.000   -0.000    0.193)    0.193
   6.019   (   0.000    0.000    0.193)    0.193
   6.115   (   0.000    0.000   -0.513)    0.513
   6.115   (   0.000    0.000   -0.513)    0.513
   9.719   (   0.000    0.000   -0.202)    0.202
   9.719   (   0.000    0.000   -0.202)    0.202
  10.735   (   0.000    0.000    1.541)    1.541
  10.735   (   0.000    0.000    1.541)    1.541
=================== 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     78.451     78.451     56.523      0.000      0.000      0.000 3/12000
   20.0     26.595     26.595     16.830      0.000      0.000      0.000 3/12000
   30.0     13.692     13.692      8.156      0.000     -0.000      0.000 3/12000
   40.0      9.044      9.044      5.135      0.000     -0.000      0.000 3/12000
   50.0      6.775      6.775      3.705      0.000     -0.000      0.000 3/12000
   60.0      5.437      5.437      2.892      0.000     -0.000      0.000 3/12000
   70.0      4.553      4.553      2.372      0.000     -0.000      0.000 3/12000
   80.0      3.924      3.924      2.013      0.000     -0.000      0.000 3/12000
   90.0      3.452      3.452      1.750      0.000     -0.000      0.000 3/12000
  100.0      3.085      3.085      1.550      0.000     -0.000      0.000 3/12000
  110.0      2.790      2.790      1.392      0.000     -0.000      0.000 3/12000
  120.0      2.549      2.549      1.264      0.000     -0.000      0.000 3/12000
  130.0      2.347      2.347      1.159      0.000     -0.000      0.000 3/12000
  140.0      2.175      2.175      1.070      0.000     -0.000      0.000 3/12000
  150.0      2.028      2.028      0.994      0.000     -0.000      0.000 3/12000
  160.0      1.899      1.899      0.929      0.000     -0.000      0.000 3/12000
  170.0      1.787      1.787      0.872      0.000     -0.000      0.000 3/12000
  180.0      1.687      1.687      0.822      0.000     -0.000      0.000 3/12000
  190.0      1.598      1.598      0.778      0.000     -0.000      0.000 3/12000
  200.0      1.518      1.518      0.738      0.000     -0.000      0.000 3/12000
  210.0      1.446      1.446      0.702      0.000     -0.000      0.000 3/12000
  220.0      1.380      1.380      0.670      0.000     -0.000      0.000 3/12000
  230.0      1.320      1.320      0.640      0.000     -0.000      0.000 3/12000
  240.0      1.265      1.265      0.613      0.000     -0.000      0.000 3/12000
  250.0      1.215      1.215      0.589      0.000     -0.000      0.000 3/12000
  260.0      1.169      1.169      0.566      0.000     -0.000      0.000 3/12000
  270.0      1.126      1.126      0.545      0.000     -0.000      0.000 3/12000
  280.0      1.086      1.086      0.525      0.000     -0.000      0.000 3/12000
  290.0      1.048      1.048      0.507      0.000     -0.000      0.000 3/12000
  300.0      1.014      1.014      0.490      0.000     -0.000      0.000 3/12000
  310.0      0.981      0.981      0.475      0.000     -0.000      0.000 3/12000
  320.0      0.951      0.951      0.460      0.000     -0.000      0.000 3/12000
  330.0      0.922      0.922      0.446      0.000     -0.000      0.000 3/12000
  340.0      0.895      0.895      0.433      0.000     -0.000      0.000 3/12000
  350.0      0.870      0.870      0.421      0.000     -0.000      0.000 3/12000
  360.0      0.846      0.846      0.409      0.000     -0.000      0.000 3/12000
  370.0      0.823      0.823      0.398      0.000     -0.000      0.000 3/12000
  380.0      0.802      0.802      0.388      0.000     -0.000      0.000 3/12000
  390.0      0.781      0.781      0.378      0.000     -0.000      0.000 3/12000
  400.0      0.762      0.762      0.368      0.000     -0.000      0.000 3/12000
  410.0      0.744      0.744      0.359      0.000     -0.000      0.000 3/12000
  420.0      0.726      0.726      0.351      0.000     -0.000      0.000 3/12000
  430.0      0.709      0.709      0.343      0.000     -0.000      0.000 3/12000
  440.0      0.693      0.693      0.335      0.000     -0.000      0.000 3/12000
  450.0      0.678      0.678      0.328      0.000     -0.000      0.000 3/12000
  460.0      0.663      0.663      0.321      0.000     -0.000      0.000 3/12000
  470.0      0.649      0.649      0.314      0.000     -0.000      0.000 3/12000
  480.0      0.636      0.636      0.307      0.000     -0.000      0.000 3/12000
  490.0      0.623      0.623      0.301      0.000     -0.000      0.000 3/12000
  500.0      0.611      0.611      0.295      0.000     -0.000      0.000 3/12000
  510.0      0.599      0.599      0.290      0.000     -0.000      0.000 3/12000
  520.0      0.587      0.587      0.284      0.000     -0.000      0.000 3/12000
  530.0      0.576      0.576      0.279      0.000     -0.000      0.000 3/12000
  540.0      0.566      0.566      0.274      0.000     -0.000      0.000 3/12000
  550.0      0.556      0.556      0.269      0.000     -0.000      0.000 3/12000
  560.0      0.546      0.546      0.264      0.000     -0.000      0.000 3/12000
  570.0      0.536      0.536      0.259      0.000     -0.000      0.000 3/12000
  580.0      0.527      0.527      0.255      0.000     -0.000      0.000 3/12000
  590.0      0.518      0.518      0.251      0.000     -0.000      0.000 3/12000
  600.0      0.510      0.510      0.247      0.000     -0.000      0.000 3/12000
  610.0      0.501      0.501      0.243      0.000     -0.000      0.000 3/12000
  620.0      0.493      0.493      0.239      0.000     -0.000      0.000 3/12000
  630.0      0.486      0.486      0.235      0.000     -0.000      0.000 3/12000
  640.0      0.478      0.478      0.231      0.000     -0.000      0.000 3/12000
  650.0      0.471      0.471      0.228      0.000     -0.000      0.000 3/12000
  660.0      0.464      0.464      0.224      0.000     -0.000      0.000 3/12000
  670.0      0.457      0.457      0.221      0.000     -0.000      0.000 3/12000
  680.0      0.450      0.450      0.218      0.000     -0.000      0.000 3/12000
  690.0      0.444      0.444      0.215      0.000     -0.000      0.000 3/12000
  700.0      0.437      0.437      0.212      0.000     -0.000      0.000 3/12000
  710.0      0.431      0.431      0.209      0.000     -0.000      0.000 3/12000
  720.0      0.425      0.425      0.206      0.000     -0.000      0.000 3/12000
  730.0      0.419      0.419      0.203      0.000     -0.000      0.000 3/12000
  740.0      0.414      0.414      0.200      0.000     -0.000      0.000 3/12000
  750.0      0.408      0.408      0.198      0.000     -0.000      0.000 3/12000
  760.0      0.403      0.403      0.195      0.000     -0.000      0.000 3/12000
  770.0      0.398      0.398      0.193      0.000     -0.000      0.000 3/12000
  780.0      0.393      0.393      0.190      0.000     -0.000      0.000 3/12000
  790.0      0.388      0.388      0.188      0.000     -0.000      0.000 3/12000
  800.0      0.383      0.383      0.185      0.000     -0.000      0.000 3/12000
  810.0      0.378      0.378      0.183      0.000     -0.000      0.000 3/12000
  820.0      0.374      0.374      0.181      0.000     -0.000      0.000 3/12000
  830.0      0.369      0.369      0.179      0.000     -0.000      0.000 3/12000
  840.0      0.365      0.365      0.177      0.000     -0.000      0.000 3/12000
  850.0      0.361      0.361      0.175      0.000     -0.000      0.000 3/12000
  860.0      0.356      0.356      0.173      0.000     -0.000      0.000 3/12000
  870.0      0.352      0.352      0.171      0.000     -0.000      0.000 3/12000
  880.0      0.348      0.348      0.169      0.000     -0.000      0.000 3/12000
  890.0      0.344      0.344      0.167      0.000     -0.000      0.000 3/12000
  900.0      0.341      0.341      0.165      0.000     -0.000      0.000 3/12000
  910.0      0.337      0.337      0.163      0.000     -0.000      0.000 3/12000
  920.0      0.333      0.333      0.161      0.000     -0.000      0.000 3/12000
  930.0      0.330      0.330      0.160      0.000     -0.000      0.000 3/12000
  940.0      0.326      0.326      0.158      0.000     -0.000      0.000 3/12000
  950.0      0.323      0.323      0.156      0.000     -0.000      0.000 3/12000
  960.0      0.319      0.319      0.155      0.000     -0.000      0.000 3/12000
  970.0      0.316      0.316      0.153      0.000     -0.000      0.000 3/12000
  980.0      0.313      0.313      0.152      0.000     -0.000      0.000 3/12000
  990.0      0.310      0.310      0.150      0.000     -0.000      0.000 3/12000
 1000.0      0.307      0.307      0.149      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:24:26]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

