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

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

-------------------------[time 2026-01-08 10:50:19]-------------------------
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/mmm (123)
Number of symmetry operations in supercell: 288
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.877260440000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.877260440000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.474208570000000
Atomic positions (fractional):
   *1 Bi  0.50000000000000  0.50000000000000  0.25903075860376 208.980
    2 Bi  0.50000000000000  0.50000000000000  0.74096924139624 208.980
   *3 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259
   *4 O   0.50000000000000  0.00000000000000  0.85883740963116  15.999
    5 O   0.50000000000000  0.00000000000000  0.14116259036884  15.999
    6 O   0.00000000000000  0.50000000000000  0.14116259036884  15.999
    7 O   0.00000000000000  0.50000000000000  0.85883740963116  15.999
   *8 I   0.00000000000000  0.00000000000000  0.50000000000000 126.904
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.877260440000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.877260440000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.474208570000000
Atomic positions (fractional):
   *1 Bi  0.50000000000000  0.50000000000000  0.25903075860376 208.980 > 1
    2 Bi  0.50000000000000  0.50000000000000  0.74096924139624 208.980 > 2
   *3 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259 > 3
   *4 O   0.50000000000000  0.00000000000000  0.85883740963116  15.999 > 4
    5 O   0.50000000000000  0.00000000000000  0.14116259036884  15.999 > 5
    6 O   0.00000000000000  0.50000000000000  0.14116259036884  15.999 > 6
    7 O   0.00000000000000  0.50000000000000  0.85883740963116  15.999 > 7
   *8 I   0.00000000000000  0.00000000000000  0.50000000000000 126.904 > 8
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   11.631781320000000    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.631781320000000    0.000000000000000
  c    0.000000000000000    0.000000000000000   18.948417140000000
Atomic positions (fractional):
   *1 Bi  0.16666666666667  0.16666666666667  0.12951537930188 208.980 > 1
    2 Bi  0.50000000000000  0.16666666666667  0.12951537930188 208.980 > 1
    3 Bi  0.83333333333333  0.16666666666667  0.12951537930188 208.980 > 1
    4 Bi  0.16666666666667  0.50000000000000  0.12951537930188 208.980 > 1
    5 Bi  0.50000000000000  0.50000000000000  0.12951537930188 208.980 > 1
    6 Bi  0.83333333333333  0.50000000000000  0.12951537930188 208.980 > 1
    7 Bi  0.16666666666667  0.83333333333333  0.12951537930188 208.980 > 1
    8 Bi  0.50000000000000  0.83333333333333  0.12951537930188 208.980 > 1
    9 Bi  0.83333333333333  0.83333333333333  0.12951537930188 208.980 > 1
   10 Bi  0.16666666666667  0.16666666666667  0.62951537930188 208.980 > 1
   11 Bi  0.50000000000000  0.16666666666667  0.62951537930188 208.980 > 1
   12 Bi  0.83333333333333  0.16666666666667  0.62951537930188 208.980 > 1
   13 Bi  0.16666666666667  0.50000000000000  0.62951537930188 208.980 > 1
   14 Bi  0.50000000000000  0.50000000000000  0.62951537930188 208.980 > 1
   15 Bi  0.83333333333333  0.50000000000000  0.62951537930188 208.980 > 1
   16 Bi  0.16666666666667  0.83333333333333  0.62951537930188 208.980 > 1
   17 Bi  0.50000000000000  0.83333333333333  0.62951537930188 208.980 > 1
   18 Bi  0.83333333333333  0.83333333333333  0.62951537930188 208.980 > 1
   19 Bi  0.16666666666667  0.16666666666667  0.37048462069812 208.980 > 2
   20 Bi  0.50000000000000  0.16666666666667  0.37048462069812 208.980 > 2
   21 Bi  0.83333333333333  0.16666666666667  0.37048462069812 208.980 > 2
   22 Bi  0.16666666666667  0.50000000000000  0.37048462069812 208.980 > 2
   23 Bi  0.50000000000000  0.50000000000000  0.37048462069812 208.980 > 2
   24 Bi  0.83333333333333  0.50000000000000  0.37048462069812 208.980 > 2
   25 Bi  0.16666666666667  0.83333333333333  0.37048462069812 208.980 > 2
   26 Bi  0.50000000000000  0.83333333333333  0.37048462069812 208.980 > 2
   27 Bi  0.83333333333333  0.83333333333333  0.37048462069812 208.980 > 2
   28 Bi  0.16666666666667  0.16666666666667  0.87048462069812 208.980 > 2
   29 Bi  0.50000000000000  0.16666666666667  0.87048462069812 208.980 > 2
   30 Bi  0.83333333333333  0.16666666666667  0.87048462069812 208.980 > 2
   31 Bi  0.16666666666667  0.50000000000000  0.87048462069812 208.980 > 2
   32 Bi  0.50000000000000  0.50000000000000  0.87048462069812 208.980 > 2
   33 Bi  0.83333333333333  0.50000000000000  0.87048462069812 208.980 > 2
   34 Bi  0.16666666666667  0.83333333333333  0.87048462069812 208.980 > 2
   35 Bi  0.50000000000000  0.83333333333333  0.87048462069812 208.980 > 2
   36 Bi  0.83333333333333  0.83333333333333  0.87048462069812 208.980 > 2
  *37 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259 > 3
   38 Er  0.33333333333333  0.00000000000000  0.00000000000000 167.259 > 3
   39 Er  0.66666666666667  0.00000000000000  0.00000000000000 167.259 > 3
   40 Er  0.00000000000000  0.33333333333333  0.00000000000000 167.259 > 3
   41 Er  0.33333333333333  0.33333333333333  0.00000000000000 167.259 > 3
   42 Er  0.66666666666667  0.33333333333333  0.00000000000000 167.259 > 3
   43 Er  0.00000000000000  0.66666666666667  0.00000000000000 167.259 > 3
   44 Er  0.33333333333333  0.66666666666667  0.00000000000000 167.259 > 3
   45 Er  0.66666666666667  0.66666666666667  0.00000000000000 167.259 > 3
   46 Er  0.00000000000000  0.00000000000000  0.50000000000000 167.259 > 3
   47 Er  0.33333333333333  0.00000000000000  0.50000000000000 167.259 > 3
   48 Er  0.66666666666667  0.00000000000000  0.50000000000000 167.259 > 3
   49 Er  0.00000000000000  0.33333333333333  0.50000000000000 167.259 > 3
   50 Er  0.33333333333333  0.33333333333333  0.50000000000000 167.259 > 3
   51 Er  0.66666666666667  0.33333333333333  0.50000000000000 167.259 > 3
   52 Er  0.00000000000000  0.66666666666667  0.50000000000000 167.259 > 3
   53 Er  0.33333333333333  0.66666666666667  0.50000000000000 167.259 > 3
   54 Er  0.66666666666667  0.66666666666667  0.50000000000000 167.259 > 3
  *55 O   0.16666666666667  0.00000000000000  0.42941870481558  15.999 > 4
   56 O   0.50000000000000  0.00000000000000  0.42941870481558  15.999 > 4
   57 O   0.83333333333333  0.00000000000000  0.42941870481558  15.999 > 4
   58 O   0.16666666666667  0.33333333333333  0.42941870481558  15.999 > 4
   59 O   0.50000000000000  0.33333333333333  0.42941870481558  15.999 > 4
   60 O   0.83333333333333  0.33333333333333  0.42941870481558  15.999 > 4
   61 O   0.16666666666667  0.66666666666667  0.42941870481558  15.999 > 4
   62 O   0.50000000000000  0.66666666666667  0.42941870481558  15.999 > 4
   63 O   0.83333333333333  0.66666666666667  0.42941870481558  15.999 > 4
   64 O   0.16666666666667  0.00000000000000  0.92941870481558  15.999 > 4
   65 O   0.50000000000000  0.00000000000000  0.92941870481558  15.999 > 4
   66 O   0.83333333333333  0.00000000000000  0.92941870481558  15.999 > 4
   67 O   0.16666666666667  0.33333333333333  0.92941870481558  15.999 > 4
   68 O   0.50000000000000  0.33333333333333  0.92941870481558  15.999 > 4
   69 O   0.83333333333333  0.33333333333333  0.92941870481558  15.999 > 4
   70 O   0.16666666666667  0.66666666666667  0.92941870481558  15.999 > 4
   71 O   0.50000000000000  0.66666666666667  0.92941870481558  15.999 > 4
   72 O   0.83333333333333  0.66666666666667  0.92941870481558  15.999 > 4
   73 O   0.16666666666667  0.00000000000000  0.07058129518442  15.999 > 5
   74 O   0.50000000000000  0.00000000000000  0.07058129518442  15.999 > 5
   75 O   0.83333333333333  0.00000000000000  0.07058129518442  15.999 > 5
   76 O   0.16666666666667  0.33333333333333  0.07058129518442  15.999 > 5
   77 O   0.50000000000000  0.33333333333333  0.07058129518442  15.999 > 5
   78 O   0.83333333333333  0.33333333333333  0.07058129518442  15.999 > 5
   79 O   0.16666666666667  0.66666666666667  0.07058129518442  15.999 > 5
   80 O   0.50000000000000  0.66666666666667  0.07058129518442  15.999 > 5
   81 O   0.83333333333333  0.66666666666667  0.07058129518442  15.999 > 5
   82 O   0.16666666666667  0.00000000000000  0.57058129518442  15.999 > 5
   83 O   0.50000000000000  0.00000000000000  0.57058129518442  15.999 > 5
   84 O   0.83333333333333  0.00000000000000  0.57058129518442  15.999 > 5
   85 O   0.16666666666667  0.33333333333333  0.57058129518442  15.999 > 5
   86 O   0.50000000000000  0.33333333333333  0.57058129518442  15.999 > 5
   87 O   0.83333333333333  0.33333333333333  0.57058129518442  15.999 > 5
   88 O   0.16666666666667  0.66666666666667  0.57058129518442  15.999 > 5
   89 O   0.50000000000000  0.66666666666667  0.57058129518442  15.999 > 5
   90 O   0.83333333333333  0.66666666666667  0.57058129518442  15.999 > 5
   91 O   0.00000000000000  0.16666666666667  0.07058129518442  15.999 > 6
   92 O   0.33333333333333  0.16666666666667  0.07058129518442  15.999 > 6
   93 O   0.66666666666667  0.16666666666667  0.07058129518442  15.999 > 6
   94 O   0.00000000000000  0.50000000000000  0.07058129518442  15.999 > 6
   95 O   0.33333333333333  0.50000000000000  0.07058129518442  15.999 > 6
   96 O   0.66666666666667  0.50000000000000  0.07058129518442  15.999 > 6
   97 O   0.00000000000000  0.83333333333333  0.07058129518442  15.999 > 6
   98 O   0.33333333333333  0.83333333333333  0.07058129518442  15.999 > 6
   99 O   0.66666666666667  0.83333333333333  0.07058129518442  15.999 > 6
  100 O   0.00000000000000  0.16666666666667  0.57058129518442  15.999 > 6
  101 O   0.33333333333333  0.16666666666667  0.57058129518442  15.999 > 6
  102 O   0.66666666666667  0.16666666666667  0.57058129518442  15.999 > 6
  103 O   0.00000000000000  0.50000000000000  0.57058129518442  15.999 > 6
  104 O   0.33333333333333  0.50000000000000  0.57058129518442  15.999 > 6
  105 O   0.66666666666667  0.50000000000000  0.57058129518442  15.999 > 6
  106 O   0.00000000000000  0.83333333333333  0.57058129518442  15.999 > 6
  107 O   0.33333333333333  0.83333333333333  0.57058129518442  15.999 > 6
  108 O   0.66666666666667  0.83333333333333  0.57058129518442  15.999 > 6
  109 O   0.00000000000000  0.16666666666667  0.42941870481558  15.999 > 7
  110 O   0.33333333333333  0.16666666666667  0.42941870481558  15.999 > 7
  111 O   0.66666666666667  0.16666666666667  0.42941870481558  15.999 > 7
  112 O   0.00000000000000  0.50000000000000  0.42941870481558  15.999 > 7
  113 O   0.33333333333333  0.50000000000000  0.42941870481558  15.999 > 7
  114 O   0.66666666666667  0.50000000000000  0.42941870481558  15.999 > 7
  115 O   0.00000000000000  0.83333333333333  0.42941870481558  15.999 > 7
  116 O   0.33333333333333  0.83333333333333  0.42941870481558  15.999 > 7
  117 O   0.66666666666667  0.83333333333333  0.42941870481558  15.999 > 7
  118 O   0.00000000000000  0.16666666666667  0.92941870481558  15.999 > 7
  119 O   0.33333333333333  0.16666666666667  0.92941870481558  15.999 > 7
  120 O   0.66666666666667  0.16666666666667  0.92941870481558  15.999 > 7
  121 O   0.00000000000000  0.50000000000000  0.92941870481558  15.999 > 7
  122 O   0.33333333333333  0.50000000000000  0.92941870481558  15.999 > 7
  123 O   0.66666666666667  0.50000000000000  0.92941870481558  15.999 > 7
  124 O   0.00000000000000  0.83333333333333  0.92941870481558  15.999 > 7
  125 O   0.33333333333333  0.83333333333333  0.92941870481558  15.999 > 7
  126 O   0.66666666666667  0.83333333333333  0.92941870481558  15.999 > 7
 *127 I   0.00000000000000  0.00000000000000  0.25000000000000 126.904 > 8
  128 I   0.33333333333333  0.00000000000000  0.25000000000000 126.904 > 8
  129 I   0.66666666666667  0.00000000000000  0.25000000000000 126.904 > 8
  130 I   0.00000000000000  0.33333333333333  0.25000000000000 126.904 > 8
  131 I   0.33333333333333  0.33333333333333  0.25000000000000 126.904 > 8
  132 I   0.66666666666667  0.33333333333333  0.25000000000000 126.904 > 8
  133 I   0.00000000000000  0.66666666666667  0.25000000000000 126.904 > 8
  134 I   0.33333333333333  0.66666666666667  0.25000000000000 126.904 > 8
  135 I   0.66666666666667  0.66666666666667  0.25000000000000 126.904 > 8
  136 I   0.00000000000000  0.00000000000000  0.75000000000000 126.904 > 8
  137 I   0.33333333333333  0.00000000000000  0.75000000000000 126.904 > 8
  138 I   0.66666666666667  0.00000000000000  0.75000000000000 126.904 > 8
  139 I   0.00000000000000  0.33333333333333  0.75000000000000 126.904 > 8
  140 I   0.33333333333333  0.33333333333333  0.75000000000000 126.904 > 8
  141 I   0.66666666666667  0.33333333333333  0.75000000000000 126.904 > 8
  142 I   0.00000000000000  0.66666666666667  0.75000000000000 126.904 > 8
  143 I   0.33333333333333  0.66666666666667  0.75000000000000 126.904 > 8
  144 I   0.66666666666667  0.66666666666667  0.75000000000000 126.904 > 8
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            7.6300485    0.0000000    0.0000000
            0.0000000    7.6300485    0.0000000
            0.0000000    0.0000000    6.1975336
-------------------------- Born effective charges --------------------------
    1 Bi    5.3567548    0.0000000    0.0000000
            0.0000000    5.3567548    0.0000000
            0.0000000    0.0000000    2.9289379
    2 Bi    5.3567548    0.0000000    0.0000000
            0.0000000    5.3567548    0.0000000
            0.0000000    0.0000000    2.9289379
    3 Er    3.7416256    0.0000000    0.0000000
            0.0000000    3.7416256    0.0000000
            0.0000000    0.0000000    5.8899948
    4 O    -2.1740292    0.0000000    0.0000000
            0.0000000   -4.0009984    0.0000000
            0.0000000    0.0000000   -2.4759477
    5 O    -2.1740292    0.0000000    0.0000000
            0.0000000   -4.0009984    0.0000000
            0.0000000    0.0000000   -2.4759477
    6 O    -4.0009984    0.0000000    0.0000000
            0.0000000   -2.1740292    0.0000000
            0.0000000    0.0000000   -2.4759477
    7 O    -4.0009984    0.0000000    0.0000000
            0.0000000   -2.1740292    0.0000000
            0.0000000    0.0000000   -2.4759477
    8 I    -2.1050801    0.0000000    0.0000000
            0.0000000   -2.1050801    0.0000000
            0.0000000    0.0000000   -1.8440798
----------------------------------------------------------------------------
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: 428
Number of blocks in projector: 428
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 8
--- Eigsh_solver_block: 1 / 8 ---
Block_size: 164
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 8 ---
Block_size: 124
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 8 ---
Block_size: 99
Use standard eigh solver.
--- Eigsh_solver_block: 4 / 8 ---
Block_size: 37
Use standard eigh solver.
--- Eigsh_solver_block: 5 / 8 ---
Block_size: 1
Use standard eigh solver.
--- Eigsh_solver_block: 6 / 8 ---
Block_size: 1
Use standard eigh solver.
--- Eigsh_solver_block: 7 / 8 ---
Block_size: 1
Use standard eigh solver.
--- Eigsh_solver_block: 8 / 8 ---
Block_size: 1
Use standard eigh solver.
Tree of FC basis block matrices:
- (428, 419), data: False
|-- (1, 1), data: True
|-- (1, 1), data: True
|-- (1, 1), data: True
|-- (1, 1), data: True
|-- (37, 37), data: True
|-- (99, 95), data: True
|-- (124, 124), data: True
|-- (164, 159), data: True
-----
Solver_atoms: 1 -- 144 / 144
Time (Solver_compr_matrix_reshape): 0.008
Solver_block: 80 / 80
 - Time: 0.378
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.436
--------------------------------- Symfc end --------------------------------
Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) 
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: 428
Number of blocks in projector: 428
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 8
--- Eigsh_solver_block: 1 / 8 ---
Block_size: 164
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 8 ---
Block_size: 124
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 8 ---
Block_size: 99
Use standard eigh solver.
--- Eigsh_solver_block: 4 / 8 ---
Block_size: 37
Use standard eigh solver.
--- Eigsh_solver_block: 5 / 8 ---
Block_size: 1
Use standard eigh solver.
--- Eigsh_solver_block: 6 / 8 ---
Block_size: 1
Use standard eigh solver.
--- Eigsh_solver_block: 7 / 8 ---
Block_size: 1
Use standard eigh solver.
--- Eigsh_solver_block: 8 / 8 ---
Block_size: 1
Use standard eigh solver.
Tree of FC basis block matrices:
- (428, 419), data: False
|-- (1, 1), data: True
|-- (1, 1), data: True
|-- (1, 1), data: True
|-- (1, 1), data: True
|-- (37, 37), data: True
|-- (99, 95), data: True
|-- (124, 124), data: True
|-- (164, 159), data: True
Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 10:50:26]-------------------------
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/mmm (123)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.877260440000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.877260440000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.474208570000000
Atomic positions (fractional):
    1 Bi  0.50000000000000  0.50000000000000  0.25903075860376 208.980
    2 Bi  0.50000000000000  0.50000000000000  0.74096924139624 208.980
    3 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259
    4 O   0.50000000000000  0.00000000000000  0.85883740963116  15.999
    5 O   0.50000000000000  0.00000000000000  0.14116259036884  15.999
    6 O   0.00000000000000  0.50000000000000  0.14116259036884  15.999
    7 O   0.00000000000000  0.50000000000000  0.85883740963116  15.999
    8 I   0.00000000000000  0.00000000000000  0.50000000000000 126.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.631781320000000    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.631781320000000    0.000000000000000
  c    0.000000000000000    0.000000000000000   18.948417140000000
Atomic positions (fractional):
    1 Bi  0.16666666666667  0.16666666666667  0.12951537930188 208.980 > 1
    2 Bi  0.50000000000000  0.16666666666667  0.12951537930188 208.980 > 1
    3 Bi  0.83333333333333  0.16666666666667  0.12951537930188 208.980 > 1
    4 Bi  0.16666666666667  0.50000000000000  0.12951537930188 208.980 > 1
    5 Bi  0.50000000000000  0.50000000000000  0.12951537930188 208.980 > 1
    6 Bi  0.83333333333333  0.50000000000000  0.12951537930188 208.980 > 1
    7 Bi  0.16666666666667  0.83333333333333  0.12951537930188 208.980 > 1
    8 Bi  0.50000000000000  0.83333333333333  0.12951537930188 208.980 > 1
    9 Bi  0.83333333333333  0.83333333333333  0.12951537930188 208.980 > 1
   10 Bi  0.16666666666667  0.16666666666667  0.62951537930188 208.980 > 1
   11 Bi  0.50000000000000  0.16666666666667  0.62951537930188 208.980 > 1
   12 Bi  0.83333333333333  0.16666666666667  0.62951537930188 208.980 > 1
   13 Bi  0.16666666666667  0.50000000000000  0.62951537930188 208.980 > 1
   14 Bi  0.50000000000000  0.50000000000000  0.62951537930188 208.980 > 1
   15 Bi  0.83333333333333  0.50000000000000  0.62951537930188 208.980 > 1
   16 Bi  0.16666666666667  0.83333333333333  0.62951537930188 208.980 > 1
   17 Bi  0.50000000000000  0.83333333333333  0.62951537930188 208.980 > 1
   18 Bi  0.83333333333333  0.83333333333333  0.62951537930188 208.980 > 1
   19 Bi  0.16666666666667  0.16666666666667  0.37048462069812 208.980 > 19
   20 Bi  0.50000000000000  0.16666666666667  0.37048462069812 208.980 > 19
   21 Bi  0.83333333333333  0.16666666666667  0.37048462069812 208.980 > 19
   22 Bi  0.16666666666667  0.50000000000000  0.37048462069812 208.980 > 19
   23 Bi  0.50000000000000  0.50000000000000  0.37048462069812 208.980 > 19
   24 Bi  0.83333333333333  0.50000000000000  0.37048462069812 208.980 > 19
   25 Bi  0.16666666666667  0.83333333333333  0.37048462069812 208.980 > 19
   26 Bi  0.50000000000000  0.83333333333333  0.37048462069812 208.980 > 19
   27 Bi  0.83333333333333  0.83333333333333  0.37048462069812 208.980 > 19
   28 Bi  0.16666666666667  0.16666666666667  0.87048462069812 208.980 > 19
   29 Bi  0.50000000000000  0.16666666666667  0.87048462069812 208.980 > 19
   30 Bi  0.83333333333333  0.16666666666667  0.87048462069812 208.980 > 19
   31 Bi  0.16666666666667  0.50000000000000  0.87048462069812 208.980 > 19
   32 Bi  0.50000000000000  0.50000000000000  0.87048462069812 208.980 > 19
   33 Bi  0.83333333333333  0.50000000000000  0.87048462069812 208.980 > 19
   34 Bi  0.16666666666667  0.83333333333333  0.87048462069812 208.980 > 19
   35 Bi  0.50000000000000  0.83333333333333  0.87048462069812 208.980 > 19
   36 Bi  0.83333333333333  0.83333333333333  0.87048462069812 208.980 > 19
   37 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259 > 37
   38 Er  0.33333333333333  0.00000000000000  0.00000000000000 167.259 > 37
   39 Er  0.66666666666667  0.00000000000000  0.00000000000000 167.259 > 37
   40 Er  0.00000000000000  0.33333333333333  0.00000000000000 167.259 > 37
   41 Er  0.33333333333333  0.33333333333333  0.00000000000000 167.259 > 37
   42 Er  0.66666666666667  0.33333333333333  0.00000000000000 167.259 > 37
   43 Er  0.00000000000000  0.66666666666667  0.00000000000000 167.259 > 37
   44 Er  0.33333333333333  0.66666666666667  0.00000000000000 167.259 > 37
   45 Er  0.66666666666667  0.66666666666667  0.00000000000000 167.259 > 37
   46 Er  0.00000000000000  0.00000000000000  0.50000000000000 167.259 > 37
   47 Er  0.33333333333333  0.00000000000000  0.50000000000000 167.259 > 37
   48 Er  0.66666666666667  0.00000000000000  0.50000000000000 167.259 > 37
   49 Er  0.00000000000000  0.33333333333333  0.50000000000000 167.259 > 37
   50 Er  0.33333333333333  0.33333333333333  0.50000000000000 167.259 > 37
   51 Er  0.66666666666667  0.33333333333333  0.50000000000000 167.259 > 37
   52 Er  0.00000000000000  0.66666666666667  0.50000000000000 167.259 > 37
   53 Er  0.33333333333333  0.66666666666667  0.50000000000000 167.259 > 37
   54 Er  0.66666666666667  0.66666666666667  0.50000000000000 167.259 > 37
   55 O   0.16666666666667  0.00000000000000  0.42941870481558  15.999 > 55
   56 O   0.50000000000000  0.00000000000000  0.42941870481558  15.999 > 55
   57 O   0.83333333333333  0.00000000000000  0.42941870481558  15.999 > 55
   58 O   0.16666666666667  0.33333333333333  0.42941870481558  15.999 > 55
   59 O   0.50000000000000  0.33333333333333  0.42941870481558  15.999 > 55
   60 O   0.83333333333333  0.33333333333333  0.42941870481558  15.999 > 55
   61 O   0.16666666666667  0.66666666666667  0.42941870481558  15.999 > 55
   62 O   0.50000000000000  0.66666666666667  0.42941870481558  15.999 > 55
   63 O   0.83333333333333  0.66666666666667  0.42941870481558  15.999 > 55
   64 O   0.16666666666667  0.00000000000000  0.92941870481558  15.999 > 55
   65 O   0.50000000000000  0.00000000000000  0.92941870481558  15.999 > 55
   66 O   0.83333333333333  0.00000000000000  0.92941870481558  15.999 > 55
   67 O   0.16666666666667  0.33333333333333  0.92941870481558  15.999 > 55
   68 O   0.50000000000000  0.33333333333333  0.92941870481558  15.999 > 55
   69 O   0.83333333333333  0.33333333333333  0.92941870481558  15.999 > 55
   70 O   0.16666666666667  0.66666666666667  0.92941870481558  15.999 > 55
   71 O   0.50000000000000  0.66666666666667  0.92941870481558  15.999 > 55
   72 O   0.83333333333333  0.66666666666667  0.92941870481558  15.999 > 55
   73 O   0.16666666666667  0.00000000000000  0.07058129518442  15.999 > 73
   74 O   0.50000000000000  0.00000000000000  0.07058129518442  15.999 > 73
   75 O   0.83333333333333  0.00000000000000  0.07058129518442  15.999 > 73
   76 O   0.16666666666667  0.33333333333333  0.07058129518442  15.999 > 73
   77 O   0.50000000000000  0.33333333333333  0.07058129518442  15.999 > 73
   78 O   0.83333333333333  0.33333333333333  0.07058129518442  15.999 > 73
   79 O   0.16666666666667  0.66666666666667  0.07058129518442  15.999 > 73
   80 O   0.50000000000000  0.66666666666667  0.07058129518442  15.999 > 73
   81 O   0.83333333333333  0.66666666666667  0.07058129518442  15.999 > 73
   82 O   0.16666666666667  0.00000000000000  0.57058129518442  15.999 > 73
   83 O   0.50000000000000  0.00000000000000  0.57058129518442  15.999 > 73
   84 O   0.83333333333333  0.00000000000000  0.57058129518442  15.999 > 73
   85 O   0.16666666666667  0.33333333333333  0.57058129518442  15.999 > 73
   86 O   0.50000000000000  0.33333333333333  0.57058129518442  15.999 > 73
   87 O   0.83333333333333  0.33333333333333  0.57058129518442  15.999 > 73
   88 O   0.16666666666667  0.66666666666667  0.57058129518442  15.999 > 73
   89 O   0.50000000000000  0.66666666666667  0.57058129518442  15.999 > 73
   90 O   0.83333333333333  0.66666666666667  0.57058129518442  15.999 > 73
   91 O   0.00000000000000  0.16666666666667  0.07058129518442  15.999 > 91
   92 O   0.33333333333333  0.16666666666667  0.07058129518442  15.999 > 91
   93 O   0.66666666666667  0.16666666666667  0.07058129518442  15.999 > 91
   94 O   0.00000000000000  0.50000000000000  0.07058129518442  15.999 > 91
   95 O   0.33333333333333  0.50000000000000  0.07058129518442  15.999 > 91
   96 O   0.66666666666667  0.50000000000000  0.07058129518442  15.999 > 91
   97 O   0.00000000000000  0.83333333333333  0.07058129518442  15.999 > 91
   98 O   0.33333333333333  0.83333333333333  0.07058129518442  15.999 > 91
   99 O   0.66666666666667  0.83333333333333  0.07058129518442  15.999 > 91
  100 O   0.00000000000000  0.16666666666667  0.57058129518442  15.999 > 91
  101 O   0.33333333333333  0.16666666666667  0.57058129518442  15.999 > 91
  102 O   0.66666666666667  0.16666666666667  0.57058129518442  15.999 > 91
  103 O   0.00000000000000  0.50000000000000  0.57058129518442  15.999 > 91
  104 O   0.33333333333333  0.50000000000000  0.57058129518442  15.999 > 91
  105 O   0.66666666666667  0.50000000000000  0.57058129518442  15.999 > 91
  106 O   0.00000000000000  0.83333333333333  0.57058129518442  15.999 > 91
  107 O   0.33333333333333  0.83333333333333  0.57058129518442  15.999 > 91
  108 O   0.66666666666667  0.83333333333333  0.57058129518442  15.999 > 91
  109 O   0.00000000000000  0.16666666666667  0.42941870481558  15.999 > 109
  110 O   0.33333333333333  0.16666666666667  0.42941870481558  15.999 > 109
  111 O   0.66666666666667  0.16666666666667  0.42941870481558  15.999 > 109
  112 O   0.00000000000000  0.50000000000000  0.42941870481558  15.999 > 109
  113 O   0.33333333333333  0.50000000000000  0.42941870481558  15.999 > 109
  114 O   0.66666666666667  0.50000000000000  0.42941870481558  15.999 > 109
  115 O   0.00000000000000  0.83333333333333  0.42941870481558  15.999 > 109
  116 O   0.33333333333333  0.83333333333333  0.42941870481558  15.999 > 109
  117 O   0.66666666666667  0.83333333333333  0.42941870481558  15.999 > 109
  118 O   0.00000000000000  0.16666666666667  0.92941870481558  15.999 > 109
  119 O   0.33333333333333  0.16666666666667  0.92941870481558  15.999 > 109
  120 O   0.66666666666667  0.16666666666667  0.92941870481558  15.999 > 109
  121 O   0.00000000000000  0.50000000000000  0.92941870481558  15.999 > 109
  122 O   0.33333333333333  0.50000000000000  0.92941870481558  15.999 > 109
  123 O   0.66666666666667  0.50000000000000  0.92941870481558  15.999 > 109
  124 O   0.00000000000000  0.83333333333333  0.92941870481558  15.999 > 109
  125 O   0.33333333333333  0.83333333333333  0.92941870481558  15.999 > 109
  126 O   0.66666666666667  0.83333333333333  0.92941870481558  15.999 > 109
  127 I   0.00000000000000  0.00000000000000  0.25000000000000 126.904 > 127
  128 I   0.33333333333333  0.00000000000000  0.25000000000000 126.904 > 127
  129 I   0.66666666666667  0.00000000000000  0.25000000000000 126.904 > 127
  130 I   0.00000000000000  0.33333333333333  0.25000000000000 126.904 > 127
  131 I   0.33333333333333  0.33333333333333  0.25000000000000 126.904 > 127
  132 I   0.66666666666667  0.33333333333333  0.25000000000000 126.904 > 127
  133 I   0.00000000000000  0.66666666666667  0.25000000000000 126.904 > 127
  134 I   0.33333333333333  0.66666666666667  0.25000000000000 126.904 > 127
  135 I   0.66666666666667  0.66666666666667  0.25000000000000 126.904 > 127
  136 I   0.00000000000000  0.00000000000000  0.75000000000000 126.904 > 127
  137 I   0.33333333333333  0.00000000000000  0.75000000000000 126.904 > 127
  138 I   0.66666666666667  0.00000000000000  0.75000000000000 126.904 > 127
  139 I   0.00000000000000  0.33333333333333  0.75000000000000 126.904 > 127
  140 I   0.33333333333333  0.33333333333333  0.75000000000000 126.904 > 127
  141 I   0.66666666666667  0.33333333333333  0.75000000000000 126.904 > 127
  142 I   0.00000000000000  0.66666666666667  0.75000000000000 126.904 > 127
  143 I   0.33333333333333  0.66666666666667  0.75000000000000 126.904 > 127
  144 I   0.66666666666667  0.66666666666667  0.75000000000000 126.904 > 127
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            7.6300485    0.0000000    0.0000000
            0.0000000    7.6300485    0.0000000
            0.0000000    0.0000000    6.1975336
-------------------------- Born effective charges --------------------------
    1 Bi    5.3567548    0.0000000    0.0000000
            0.0000000    5.3567548    0.0000000
            0.0000000    0.0000000    2.9289379
    2 Bi    5.3567548    0.0000000    0.0000000
            0.0000000    5.3567548    0.0000000
            0.0000000    0.0000000    2.9289379
    3 Er    3.7416256    0.0000000    0.0000000
            0.0000000    3.7416256    0.0000000
            0.0000000    0.0000000    5.8899948
    4 O    -2.1740292    0.0000000    0.0000000
            0.0000000   -4.0009984    0.0000000
            0.0000000    0.0000000   -2.4759477
    5 O    -2.1740292    0.0000000    0.0000000
            0.0000000   -4.0009984    0.0000000
            0.0000000    0.0000000   -2.4759477
    6 O    -4.0009984    0.0000000    0.0000000
            0.0000000   -2.1740292    0.0000000
            0.0000000    0.0000000   -2.4759477
    7 O    -4.0009984    0.0000000    0.0000000
            0.0000000   -2.1740292    0.0000000
            0.0000000    0.0000000   -2.4759477
    8 I    -2.1050801    0.0000000    0.0000000
            0.0000000   -2.1050801    0.0000000
            0.0000000    0.0000000   -1.8440798
----------------------------------------------------------------------------
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[ 55, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Computing fc3[ 127, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000000 (xxz) 0.00000000 (xxz) 0.00000000 (xzx)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 10:50:35]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 10:50:36]-------------------------
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/mmm (123)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.877260440000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.877260440000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.474208570000000
Atomic positions (fractional):
    1 Bi  0.50000000000000  0.50000000000000  0.25903075860376 208.980
    2 Bi  0.50000000000000  0.50000000000000  0.74096924139624 208.980
    3 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259
    4 O   0.50000000000000  0.00000000000000  0.85883740963116  15.999
    5 O   0.50000000000000  0.00000000000000  0.14116259036884  15.999
    6 O   0.00000000000000  0.50000000000000  0.14116259036884  15.999
    7 O   0.00000000000000  0.50000000000000  0.85883740963116  15.999
    8 I   0.00000000000000  0.00000000000000  0.50000000000000 126.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.631781320000000    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.631781320000000    0.000000000000000
  c    0.000000000000000    0.000000000000000   18.948417140000000
Atomic positions (fractional):
    1 Bi  0.16666666666667  0.16666666666667  0.12951537930188 208.980 > 1
    2 Bi  0.50000000000000  0.16666666666667  0.12951537930188 208.980 > 1
    3 Bi  0.83333333333333  0.16666666666667  0.12951537930188 208.980 > 1
    4 Bi  0.16666666666667  0.50000000000000  0.12951537930188 208.980 > 1
    5 Bi  0.50000000000000  0.50000000000000  0.12951537930188 208.980 > 1
    6 Bi  0.83333333333333  0.50000000000000  0.12951537930188 208.980 > 1
    7 Bi  0.16666666666667  0.83333333333333  0.12951537930188 208.980 > 1
    8 Bi  0.50000000000000  0.83333333333333  0.12951537930188 208.980 > 1
    9 Bi  0.83333333333333  0.83333333333333  0.12951537930188 208.980 > 1
   10 Bi  0.16666666666667  0.16666666666667  0.62951537930188 208.980 > 1
   11 Bi  0.50000000000000  0.16666666666667  0.62951537930188 208.980 > 1
   12 Bi  0.83333333333333  0.16666666666667  0.62951537930188 208.980 > 1
   13 Bi  0.16666666666667  0.50000000000000  0.62951537930188 208.980 > 1
   14 Bi  0.50000000000000  0.50000000000000  0.62951537930188 208.980 > 1
   15 Bi  0.83333333333333  0.50000000000000  0.62951537930188 208.980 > 1
   16 Bi  0.16666666666667  0.83333333333333  0.62951537930188 208.980 > 1
   17 Bi  0.50000000000000  0.83333333333333  0.62951537930188 208.980 > 1
   18 Bi  0.83333333333333  0.83333333333333  0.62951537930188 208.980 > 1
   19 Bi  0.16666666666667  0.16666666666667  0.37048462069812 208.980 > 19
   20 Bi  0.50000000000000  0.16666666666667  0.37048462069812 208.980 > 19
   21 Bi  0.83333333333333  0.16666666666667  0.37048462069812 208.980 > 19
   22 Bi  0.16666666666667  0.50000000000000  0.37048462069812 208.980 > 19
   23 Bi  0.50000000000000  0.50000000000000  0.37048462069812 208.980 > 19
   24 Bi  0.83333333333333  0.50000000000000  0.37048462069812 208.980 > 19
   25 Bi  0.16666666666667  0.83333333333333  0.37048462069812 208.980 > 19
   26 Bi  0.50000000000000  0.83333333333333  0.37048462069812 208.980 > 19
   27 Bi  0.83333333333333  0.83333333333333  0.37048462069812 208.980 > 19
   28 Bi  0.16666666666667  0.16666666666667  0.87048462069812 208.980 > 19
   29 Bi  0.50000000000000  0.16666666666667  0.87048462069812 208.980 > 19
   30 Bi  0.83333333333333  0.16666666666667  0.87048462069812 208.980 > 19
   31 Bi  0.16666666666667  0.50000000000000  0.87048462069812 208.980 > 19
   32 Bi  0.50000000000000  0.50000000000000  0.87048462069812 208.980 > 19
   33 Bi  0.83333333333333  0.50000000000000  0.87048462069812 208.980 > 19
   34 Bi  0.16666666666667  0.83333333333333  0.87048462069812 208.980 > 19
   35 Bi  0.50000000000000  0.83333333333333  0.87048462069812 208.980 > 19
   36 Bi  0.83333333333333  0.83333333333333  0.87048462069812 208.980 > 19
   37 Er  0.00000000000000  0.00000000000000  0.00000000000000 167.259 > 37
   38 Er  0.33333333333333  0.00000000000000  0.00000000000000 167.259 > 37
   39 Er  0.66666666666667  0.00000000000000  0.00000000000000 167.259 > 37
   40 Er  0.00000000000000  0.33333333333333  0.00000000000000 167.259 > 37
   41 Er  0.33333333333333  0.33333333333333  0.00000000000000 167.259 > 37
   42 Er  0.66666666666667  0.33333333333333  0.00000000000000 167.259 > 37
   43 Er  0.00000000000000  0.66666666666667  0.00000000000000 167.259 > 37
   44 Er  0.33333333333333  0.66666666666667  0.00000000000000 167.259 > 37
   45 Er  0.66666666666667  0.66666666666667  0.00000000000000 167.259 > 37
   46 Er  0.00000000000000  0.00000000000000  0.50000000000000 167.259 > 37
   47 Er  0.33333333333333  0.00000000000000  0.50000000000000 167.259 > 37
   48 Er  0.66666666666667  0.00000000000000  0.50000000000000 167.259 > 37
   49 Er  0.00000000000000  0.33333333333333  0.50000000000000 167.259 > 37
   50 Er  0.33333333333333  0.33333333333333  0.50000000000000 167.259 > 37
   51 Er  0.66666666666667  0.33333333333333  0.50000000000000 167.259 > 37
   52 Er  0.00000000000000  0.66666666666667  0.50000000000000 167.259 > 37
   53 Er  0.33333333333333  0.66666666666667  0.50000000000000 167.259 > 37
   54 Er  0.66666666666667  0.66666666666667  0.50000000000000 167.259 > 37
   55 O   0.16666666666667  0.00000000000000  0.42941870481558  15.999 > 55
   56 O   0.50000000000000  0.00000000000000  0.42941870481558  15.999 > 55
   57 O   0.83333333333333  0.00000000000000  0.42941870481558  15.999 > 55
   58 O   0.16666666666667  0.33333333333333  0.42941870481558  15.999 > 55
   59 O   0.50000000000000  0.33333333333333  0.42941870481558  15.999 > 55
   60 O   0.83333333333333  0.33333333333333  0.42941870481558  15.999 > 55
   61 O   0.16666666666667  0.66666666666667  0.42941870481558  15.999 > 55
   62 O   0.50000000000000  0.66666666666667  0.42941870481558  15.999 > 55
   63 O   0.83333333333333  0.66666666666667  0.42941870481558  15.999 > 55
   64 O   0.16666666666667  0.00000000000000  0.92941870481558  15.999 > 55
   65 O   0.50000000000000  0.00000000000000  0.92941870481558  15.999 > 55
   66 O   0.83333333333333  0.00000000000000  0.92941870481558  15.999 > 55
   67 O   0.16666666666667  0.33333333333333  0.92941870481558  15.999 > 55
   68 O   0.50000000000000  0.33333333333333  0.92941870481558  15.999 > 55
   69 O   0.83333333333333  0.33333333333333  0.92941870481558  15.999 > 55
   70 O   0.16666666666667  0.66666666666667  0.92941870481558  15.999 > 55
   71 O   0.50000000000000  0.66666666666667  0.92941870481558  15.999 > 55
   72 O   0.83333333333333  0.66666666666667  0.92941870481558  15.999 > 55
   73 O   0.16666666666667  0.00000000000000  0.07058129518442  15.999 > 73
   74 O   0.50000000000000  0.00000000000000  0.07058129518442  15.999 > 73
   75 O   0.83333333333333  0.00000000000000  0.07058129518442  15.999 > 73
   76 O   0.16666666666667  0.33333333333333  0.07058129518442  15.999 > 73
   77 O   0.50000000000000  0.33333333333333  0.07058129518442  15.999 > 73
   78 O   0.83333333333333  0.33333333333333  0.07058129518442  15.999 > 73
   79 O   0.16666666666667  0.66666666666667  0.07058129518442  15.999 > 73
   80 O   0.50000000000000  0.66666666666667  0.07058129518442  15.999 > 73
   81 O   0.83333333333333  0.66666666666667  0.07058129518442  15.999 > 73
   82 O   0.16666666666667  0.00000000000000  0.57058129518442  15.999 > 73
   83 O   0.50000000000000  0.00000000000000  0.57058129518442  15.999 > 73
   84 O   0.83333333333333  0.00000000000000  0.57058129518442  15.999 > 73
   85 O   0.16666666666667  0.33333333333333  0.57058129518442  15.999 > 73
   86 O   0.50000000000000  0.33333333333333  0.57058129518442  15.999 > 73
   87 O   0.83333333333333  0.33333333333333  0.57058129518442  15.999 > 73
   88 O   0.16666666666667  0.66666666666667  0.57058129518442  15.999 > 73
   89 O   0.50000000000000  0.66666666666667  0.57058129518442  15.999 > 73
   90 O   0.83333333333333  0.66666666666667  0.57058129518442  15.999 > 73
   91 O   0.00000000000000  0.16666666666667  0.07058129518442  15.999 > 91
   92 O   0.33333333333333  0.16666666666667  0.07058129518442  15.999 > 91
   93 O   0.66666666666667  0.16666666666667  0.07058129518442  15.999 > 91
   94 O   0.00000000000000  0.50000000000000  0.07058129518442  15.999 > 91
   95 O   0.33333333333333  0.50000000000000  0.07058129518442  15.999 > 91
   96 O   0.66666666666667  0.50000000000000  0.07058129518442  15.999 > 91
   97 O   0.00000000000000  0.83333333333333  0.07058129518442  15.999 > 91
   98 O   0.33333333333333  0.83333333333333  0.07058129518442  15.999 > 91
   99 O   0.66666666666667  0.83333333333333  0.07058129518442  15.999 > 91
  100 O   0.00000000000000  0.16666666666667  0.57058129518442  15.999 > 91
  101 O   0.33333333333333  0.16666666666667  0.57058129518442  15.999 > 91
  102 O   0.66666666666667  0.16666666666667  0.57058129518442  15.999 > 91
  103 O   0.00000000000000  0.50000000000000  0.57058129518442  15.999 > 91
  104 O   0.33333333333333  0.50000000000000  0.57058129518442  15.999 > 91
  105 O   0.66666666666667  0.50000000000000  0.57058129518442  15.999 > 91
  106 O   0.00000000000000  0.83333333333333  0.57058129518442  15.999 > 91
  107 O   0.33333333333333  0.83333333333333  0.57058129518442  15.999 > 91
  108 O   0.66666666666667  0.83333333333333  0.57058129518442  15.999 > 91
  109 O   0.00000000000000  0.16666666666667  0.42941870481558  15.999 > 109
  110 O   0.33333333333333  0.16666666666667  0.42941870481558  15.999 > 109
  111 O   0.66666666666667  0.16666666666667  0.42941870481558  15.999 > 109
  112 O   0.00000000000000  0.50000000000000  0.42941870481558  15.999 > 109
  113 O   0.33333333333333  0.50000000000000  0.42941870481558  15.999 > 109
  114 O   0.66666666666667  0.50000000000000  0.42941870481558  15.999 > 109
  115 O   0.00000000000000  0.83333333333333  0.42941870481558  15.999 > 109
  116 O   0.33333333333333  0.83333333333333  0.42941870481558  15.999 > 109
  117 O   0.66666666666667  0.83333333333333  0.42941870481558  15.999 > 109
  118 O   0.00000000000000  0.16666666666667  0.92941870481558  15.999 > 109
  119 O   0.33333333333333  0.16666666666667  0.92941870481558  15.999 > 109
  120 O   0.66666666666667  0.16666666666667  0.92941870481558  15.999 > 109
  121 O   0.00000000000000  0.50000000000000  0.92941870481558  15.999 > 109
  122 O   0.33333333333333  0.50000000000000  0.92941870481558  15.999 > 109
  123 O   0.66666666666667  0.50000000000000  0.92941870481558  15.999 > 109
  124 O   0.00000000000000  0.83333333333333  0.92941870481558  15.999 > 109
  125 O   0.33333333333333  0.83333333333333  0.92941870481558  15.999 > 109
  126 O   0.66666666666667  0.83333333333333  0.92941870481558  15.999 > 109
  127 I   0.00000000000000  0.00000000000000  0.25000000000000 126.904 > 127
  128 I   0.33333333333333  0.00000000000000  0.25000000000000 126.904 > 127
  129 I   0.66666666666667  0.00000000000000  0.25000000000000 126.904 > 127
  130 I   0.00000000000000  0.33333333333333  0.25000000000000 126.904 > 127
  131 I   0.33333333333333  0.33333333333333  0.25000000000000 126.904 > 127
  132 I   0.66666666666667  0.33333333333333  0.25000000000000 126.904 > 127
  133 I   0.00000000000000  0.66666666666667  0.25000000000000 126.904 > 127
  134 I   0.33333333333333  0.66666666666667  0.25000000000000 126.904 > 127
  135 I   0.66666666666667  0.66666666666667  0.25000000000000 126.904 > 127
  136 I   0.00000000000000  0.00000000000000  0.75000000000000 126.904 > 127
  137 I   0.33333333333333  0.00000000000000  0.75000000000000 126.904 > 127
  138 I   0.66666666666667  0.00000000000000  0.75000000000000 126.904 > 127
  139 I   0.00000000000000  0.33333333333333  0.75000000000000 126.904 > 127
  140 I   0.33333333333333  0.33333333333333  0.75000000000000 126.904 > 127
  141 I   0.66666666666667  0.33333333333333  0.75000000000000 126.904 > 127
  142 I   0.00000000000000  0.66666666666667  0.75000000000000 126.904 > 127
  143 I   0.33333333333333  0.66666666666667  0.75000000000000 126.904 > 127
  144 I   0.66666666666667  0.66666666666667  0.75000000000000 126.904 > 127
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            7.6300485    0.0000000    0.0000000
            0.0000000    7.6300485    0.0000000
            0.0000000    0.0000000    6.1975336
-------------------------- Born effective charges --------------------------
    1 Bi    5.3567548    0.0000000    0.0000000
            0.0000000    5.3567548    0.0000000
            0.0000000    0.0000000    2.9289379
    2 Bi    5.3567548    0.0000000    0.0000000
            0.0000000    5.3567548    0.0000000
            0.0000000    0.0000000    2.9289379
    3 Er    3.7416256    0.0000000    0.0000000
            0.0000000    3.7416256    0.0000000
            0.0000000    0.0000000    5.8899948
    4 O    -2.1740292    0.0000000    0.0000000
            0.0000000   -4.0009984    0.0000000
            0.0000000    0.0000000   -2.4759477
    5 O    -2.1740292    0.0000000    0.0000000
            0.0000000   -4.0009984    0.0000000
            0.0000000    0.0000000   -2.4759477
    6 O    -4.0009984    0.0000000    0.0000000
            0.0000000   -2.1740292    0.0000000
            0.0000000    0.0000000   -2.4759477
    7 O    -4.0009984    0.0000000    0.0000000
            0.0000000   -2.1740292    0.0000000
            0.0000000    0.0000000   -2.4759477
    8 I    -2.1050801    0.0000000    0.0000000
            0.0000000   -2.1050801    0.0000000
            0.0000000    0.0000000   -1.8440798
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (xxz) 0.00000000 (xxz) 0.00000000 (xzx)
Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 13 13 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.80, Number of G-points: 291, Lambda: 0.22
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/84) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   1.784   (   0.000    0.000    0.000)    0.000
   1.784   (   0.000    0.000    0.000)    0.000
   1.882   (   0.000    0.000    0.000)    0.000
   1.882   (   0.000    0.000    0.000)    0.000
   2.120   (   0.000    0.000    0.000)    0.000
   3.110   (   0.000    0.000    0.000)    0.000
   3.110   (   0.000    0.000    0.000)    0.000
   3.693   (   0.000    0.000    0.000)    0.000
   5.549   (   0.000    0.000    0.000)    0.000
   8.111   (   0.000    0.000    0.000)    0.000
   8.111   (   0.000    0.000    0.000)    0.000
   8.612   (   0.000    0.000    0.000)    0.000
   8.612   (   0.000    0.000    0.000)    0.000
   9.182   (   0.000    0.000    0.000)    0.000
  11.833   (   0.000    0.000    0.000)    0.000
  13.060   (   0.000    0.000    0.000)    0.000
  13.692   (   0.000    0.000    0.000)    0.000
  13.692   (   0.000    0.000    0.000)    0.000
  13.832   (   0.000    0.000    0.000)    0.000
  13.832   (   0.000    0.000    0.000)    0.000
  16.042   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/84) =======================
q-point: ( 0.08  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.321   (  17.377    0.000    0.000)   17.377
   0.585   (  28.054    0.000    0.000)   28.054
   0.851   (  41.653    0.000    0.000)   41.653
   1.825   (   3.893    0.000    0.000)    3.893
   1.827   (   3.057    0.000    0.000)    3.057
   1.887   (   0.534    0.000    0.000)    0.534
   2.193   (   7.947    0.000    0.000)    7.947
   2.252   (   6.067    0.000    0.000)    6.067
   3.102   (  -0.842    0.000    0.000)    0.842
   3.191   (   3.827    0.000    0.000)    3.827
   3.735   (   4.687    0.000    0.000)    4.687
   5.403   ( -13.243    0.000    0.000)   13.243
   8.239   (  12.652    0.000    0.000)   12.652
   8.738   (  12.588    0.000    0.000)   12.588
   8.748   (  10.715    0.000    0.000)   10.715
   9.157   (  -0.806    0.000    0.000)    0.806
  11.414   (   6.446    0.000    0.000)    6.446
  11.819   (  -1.457    0.000    0.000)    1.457
  13.037   (  -2.167    0.000    0.000)    2.167
  13.649   (  -4.358    0.000    0.000)    4.358
  13.767   (   7.227    0.000    0.000)    7.227
  13.807   (  -2.474    0.000    0.000)    2.474
  15.698   (  -5.705    0.000    0.000)    5.705
  16.017   (  -2.844    0.000    0.000)    2.844
======================= Grid point 2 (3/84) =======================
q-point: ( 0.15  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.693   (  19.915    0.000    0.000)   19.915
   1.090   (  22.420    0.000    0.000)   22.420
   1.628   (  35.955    0.000    0.000)   35.955
   1.866   (   0.757    0.000    0.000)    0.757
   1.905   (   1.448    0.000    0.000)    1.448
   1.923   (   5.578    0.000    0.000)    5.578
   2.400   (   7.844    0.000    0.000)    7.844
   2.432   (  15.017    0.000    0.000)   15.017
   3.077   (  -1.613    0.000    0.000)    1.613
   3.278   (   5.049    0.000    0.000)    5.049
   3.896   (  11.491    0.000    0.000)   11.491
   5.097   ( -15.921    0.000    0.000)   15.921
   8.596   (  22.660    0.000    0.000)   22.660
   8.791   (  -8.367    0.000    0.000)    8.367
   9.096   (  22.937    0.000    0.000)   22.937
   9.327   (  18.595    0.000    0.000)   18.595
  11.589   (  10.772    0.000    0.000)   10.772
  11.774   (  -3.108    0.000    0.000)    3.108
  12.987   (  -2.458    0.000    0.000)    2.458
  13.520   (  -8.674    0.000    0.000)    8.674
  13.734   (  -4.851    0.000    0.000)    4.851
  13.959   (  11.437    0.000    0.000)   11.437
  15.552   (  -8.585    0.000    0.000)    8.585
  15.919   (  -7.167    0.000    0.000)    7.167
======================= Grid point 3 (4/84) =======================
q-point: ( 0.23  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.099   (  20.657    0.000    0.000)   20.657
   1.459   (  14.483    0.000    0.000)   14.483
   1.869   (  -0.134    0.000    0.000)    0.134
   1.952   (   3.607    0.000    0.000)    3.607
   2.032   (   5.123    0.000    0.000)    5.123
   2.248   (  25.865    0.000    0.000)   25.865
   2.524   (   3.888    0.000    0.000)    3.888
   2.732   (  13.891    0.000    0.000)   13.891
   3.039   (  -2.219    0.000    0.000)    2.219
   3.428   (  11.207    0.000    0.000)   11.207
   4.167   (  14.972    0.000    0.000)   14.972
   4.824   ( -10.697    0.000    0.000)   10.697
   8.534   ( -15.682    0.000    0.000)   15.682
   9.107   (  27.937    0.000    0.000)   27.937
   9.624   (  29.579    0.000    0.000)   29.579
   9.746   (  21.441    0.000    0.000)   21.441
  11.695   (  -4.710    0.000    0.000)    4.710
  11.814   (  11.127    0.000    0.000)   11.127
  12.962   (   0.509    0.000    0.000)    0.509
  13.305   ( -12.916    0.000    0.000)   12.916
  13.616   (  -6.925    0.000    0.000)    6.925
  14.194   (  11.674    0.000    0.000)   11.674
  15.366   ( -10.154    0.000    0.000)   10.154
  15.739   ( -10.512    0.000    0.000)   10.512
======================= Grid point 4 (5/84) =======================
q-point: ( 0.31  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.504   (  19.950    0.000    0.000)   19.950
   1.660   (   6.175    0.000    0.000)    6.175
   1.869   (   0.301    0.000    0.000)    0.301
   2.055   (   6.483    0.000    0.000)    6.483
   2.120   (   3.697    0.000    0.000)    3.697
   2.548   (  -1.031    0.000    0.000)    1.031
   2.645   (  14.672    0.000    0.000)   14.672
   2.939   (   6.287    0.000    0.000)    6.287
   2.991   (  -2.507    0.000    0.000)    2.507
   3.728   (  17.880    0.000    0.000)   17.880
   4.458   (  13.685    0.000    0.000)   13.685
   4.694   (  -2.262    0.000    0.000)    2.262
   8.214   ( -15.791    0.000    0.000)   15.791
   9.666   (  27.425    0.000    0.000)   27.425
  10.134   (  16.998    0.000    0.000)   16.998
  10.243   (  32.159    0.000    0.000)   32.159
  11.595   (  -5.131    0.000    0.000)    5.131
  12.003   (   7.387    0.000    0.000)    7.387
  13.008   ( -16.998    0.000    0.000)   16.998
  13.025   (   5.952    0.000    0.000)    5.952
  13.465   (  -8.170    0.000    0.000)    8.170
  14.404   (   9.159    0.000    0.000)    9.159
  15.152   ( -11.268    0.000    0.000)   11.268
  15.521   ( -10.915    0.000    0.000)   10.915
======================= Grid point 5 (6/84) =======================
q-point: ( 0.38  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.734   (   2.004    0.000    0.000)    2.004
   1.874   (  11.600    0.000    0.000)   11.600
   1.884   (   6.353    0.000    0.000)    6.353
   2.178   (   2.129    0.000    0.000)    2.129
   2.183   (   5.723    0.000    0.000)    5.723
   2.509   (  -2.293    0.000    0.000)    2.293
   2.858   (   7.351    0.000    0.000)    7.351
   2.943   (  -2.194    0.000    0.000)    2.194
   2.979   (  -1.691    0.000    0.000)    1.691
   4.073   (  15.344    0.000    0.000)   15.344
   4.689   (   9.207    0.000    0.000)    9.207
   4.711   (   2.926    0.000    0.000)    2.926
   7.937   ( -11.518    0.000    0.000)   11.518
  10.154   (  20.741    0.000    0.000)   20.741
  10.405   (  10.170    0.000    0.000)   10.170
  10.874   (  30.827    0.000    0.000)   30.827
  11.504   (  -3.763    0.000    0.000)    3.763
  12.105   (   3.249    0.000    0.000)    3.249
  12.635   ( -20.377    0.000    0.000)   20.377
  13.179   (   8.612    0.000    0.000)    8.612
  13.305   (  -7.413    0.000    0.000)    7.413
  14.551   (   5.583    0.000    0.000)    5.583
  14.932   ( -10.333    0.000    0.000)   10.333
  15.328   (  -8.143    0.000    0.000)    8.143
======================= Grid point 6 (7/84) =======================
q-point: ( 0.46  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.757   (   0.482    0.000    0.000)    0.482
   1.894   (   0.409    0.000    0.000)    0.409
   2.146   (   8.250    0.000    0.000)    8.250
   2.206   (   0.686    0.000    0.000)    0.686
   2.264   (   2.196    0.000    0.000)    2.196
   2.473   (  -1.021    0.000    0.000)    1.021
   2.911   (  -0.932    0.000    0.000)    0.932
   2.915   (  -3.080    0.000    0.000)    3.080
   2.952   (   2.288    0.000    0.000)    2.288
   4.291   (   5.915    0.000    0.000)    5.915
   4.767   (   1.856    0.000    0.000)    1.856
   4.813   (   3.170    0.000    0.000)    3.170
   7.777   (  -4.218    0.000    0.000)    4.218
  10.448   (   7.947    0.000    0.000)    7.947
  10.537   (   3.238    0.000    0.000)    3.238
  11.419   (  21.078    0.000    0.000)   21.078
  11.453   (  -1.309    0.000    0.000)    1.309
  12.144   (   0.866    0.000    0.000)    0.866
  12.234   ( -17.366    0.000    0.000)   17.366
  13.194   (  -3.204    0.000    0.000)    3.204
  13.319   (   4.295    0.000    0.000)    4.295
  14.624   (   1.853    0.000    0.000)    1.853
  14.774   (  -4.597    0.000    0.000)    4.597
  15.214   (  -3.042    0.000    0.000)    3.042
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.459   (  12.543   12.543    0.000)   17.739
   0.541   (  13.358   13.358    0.000)   18.891
   1.296   (  28.207   28.207    0.000)   39.891
   1.821   (   1.481    1.481    0.000)    2.094
   1.931   (   6.905    6.905    0.000)    9.765
   1.943   (   3.152    3.152    0.000)    4.457
   2.240   (   5.505    5.505    0.000)    7.785
   2.319   (   6.734    6.734    0.000)    9.523
   3.088   (  -2.370   -2.370    0.000)    3.351
   3.260   (   7.297    7.297    0.000)   10.319
   3.759   (   2.971    2.971    0.000)    4.201
   5.300   (  -9.842   -9.842    0.000)   13.918
   8.128   (   0.676    0.676    0.000)    0.956
   8.514   (  -4.998   -4.998    0.000)    7.068
   9.108   (  -3.778   -3.778    0.000)    5.343
   9.168   (  22.867   22.867    0.000)   32.339
  11.570   (  10.685   10.685    0.000)   15.112
  11.819   (  -0.029   -0.029    0.000)    0.041
  12.976   (  -4.655   -4.655    0.000)    6.584
  13.537   (  -7.769   -7.769    0.000)   10.988
  13.707   (  -5.972   -5.972    0.000)    8.446
  13.907   (  10.317   10.317    0.000)   14.591
  15.729   (  -3.156   -3.156    0.000)    4.464
  16.043   (   1.811    1.811    0.000)    2.561
======================= Grid point 15 (9/84) =======================
q-point: ( 0.15  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.771   (  18.025    7.772    0.000)   19.629
   0.977   (  24.878   -4.396    0.000)   25.263
   1.796   (  18.173   13.216    0.000)   22.470
   1.857   (   1.645    0.623    0.000)    1.759
   2.028   (   9.064    7.754    0.000)   11.928
   2.053   (   5.130   10.604    0.000)   11.780
   2.436   (  13.593    0.636    0.000)   13.608
   2.485   (   8.704    9.015    0.000)   12.531
   3.046   (  -1.447   -1.349    0.000)    1.979
   3.407   (   6.724   11.033    0.000)   12.921
   3.887   (  10.491   -0.458    0.000)   10.501
   5.054   ( -13.636   -4.352    0.000)   14.314
   8.340   (  18.638  -17.704    0.000)   25.706
   8.525   (   3.854  -19.342    0.000)   19.722
   8.989   (  -5.622   -4.521    0.000)    7.214
   9.651   (  24.399   19.193    0.000)   31.043
  11.789   (  11.444   16.248    0.000)   19.873
  11.813   (  -0.653    3.947    0.000)    4.001
  12.879   (  -4.838  -10.194    0.000)   11.284
  13.360   ( -10.280  -12.596    0.000)   16.258
  13.594   (  -5.800  -10.173    0.000)   11.711
  14.140   (  12.771   14.502    0.000)   19.324
  15.590   (  -9.335   -1.007    0.000)    9.389
  16.033   (  -3.219   10.658    0.000)   11.134
======================= Grid point 16 (10/84) =======================
q-point: ( 0.23  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.150   (  19.802    5.315    0.000)   20.503
   1.425   (  19.228   -2.788    0.000)   19.429
   1.876   (   0.378    1.005    0.000)    1.074
   1.953   (   3.560    4.398    0.000)    5.659
   2.149   (   4.808   10.767    0.000)   11.792
   2.354   (  17.533    4.116    0.000)   18.010
   2.621   (   4.589    9.644    0.000)   10.680
   2.721   (  13.692   -0.979    0.000)   13.727
   3.031   (  -0.521    1.295    0.000)    1.396
   3.561   (  10.283   11.921    0.000)   15.744
   4.153   (  15.171   -1.145    0.000)   15.214
   4.807   ( -10.208   -1.688    0.000)   10.346
   8.484   ( -10.247   -6.759    0.000)   12.275
   8.800   (  26.463  -25.739    0.000)   36.916
   9.047   (  13.326  -30.343    0.000)   33.140
  10.117   (  22.227   12.118    0.000)   25.315
  11.790   (  -1.694    9.472    0.000)    9.622
  12.007   (   9.716   16.196    0.000)   18.887
  12.802   (  -2.534  -15.648    0.000)   15.852
  13.123   ( -13.613  -15.499    0.000)   20.629
  13.472   (  -6.234  -10.679    0.000)   12.365
  14.395   (  12.462   17.078    0.000)   21.141
  15.386   ( -10.943   -0.810    0.000)   10.973
  15.925   (  -7.256   16.834    0.000)   18.331
======================= Grid point 17 (11/84) =======================
q-point: ( 0.31  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.544   (  19.659    4.254    0.000)   20.114
   1.711   (   9.280    4.050    0.000)   10.125
   1.881   (   0.305    1.171    0.000)    1.210
   2.039   (   5.704    0.942    0.000)    5.781
   2.241   (   4.266   11.284    0.000)   12.064
   2.579   (   3.351    0.872    0.000)    3.462
   2.727   (   8.411    6.998    0.000)   10.941
   2.931   (   6.617   -0.729    0.000)    6.657
   3.007   (  -1.939    3.191    0.000)    3.734
   3.832   (  16.088    9.776    0.000)   18.825
   4.450   (  14.005   -0.592    0.000)   14.018
   4.677   (  -2.656   -1.710    0.000)    3.159
   8.213   ( -14.872   -0.321    0.000)   14.875
   9.346   (  27.491  -28.524    0.000)   39.616
   9.407   (  20.366  -44.344    0.000)   48.797
  10.521   (  18.314    0.915    0.000)   18.336
  11.754   (  -1.559   15.458    0.000)   15.536
  12.149   (   4.449   10.648    0.000)   11.540
  12.773   (  -1.468  -22.093    0.000)   22.142
  12.821   ( -16.742  -16.375    0.000)   23.419
  13.373   (  -2.750   -6.718    0.000)    7.260
  14.618   (   9.770   18.739    0.000)   21.133
  15.163   ( -11.257   -0.132    0.000)   11.258
  15.765   (  -8.336   21.228    0.000)   22.805
======================= Grid point 18 (12/84) =======================
q-point: ( 0.38  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.818   (   2.647    7.876    0.000)    8.309
   1.889   (   1.561    1.078    0.000)    1.897
   1.918   (  16.376    4.006    0.000)   16.858
   2.162   (   5.720   -0.549    0.000)    5.747
   2.313   (   2.842   12.613    0.000)   12.929
   2.587   (  -0.628    5.696    0.000)    5.730
   2.880   (   5.289    1.054    0.000)    5.393
   2.967   (  -1.277    4.019    0.000)    4.217
   2.976   (  -1.587   -0.186    0.000)    1.598
   4.143   (  13.990    6.896    0.000)   15.597
   4.686   (   2.611   -2.518    0.000)    3.627
   4.687   (   9.453   -0.070    0.000)    9.453
   7.944   ( -11.431    0.488    0.000)   11.441
   9.789   (  16.862  -45.252    0.000)   48.292
   9.849   (  22.063  -27.830    0.000)   35.514
  10.835   (  13.083  -16.053    0.000)   20.709
  11.747   (   1.405   20.660    0.000)   20.708
  12.204   (   1.902    4.630    0.000)    5.005
  12.470   ( -18.162  -13.896    0.000)   22.868
  12.691   (  -6.987  -31.944    0.000)   32.700
  13.392   (   4.161   -0.547    0.000)    4.197
  14.775   (   5.991   19.822    0.000)   20.708
  14.953   (  -9.448    1.795    0.000)    9.617
  15.613   (  -6.482   24.210    0.000)   25.063
======================= Grid point 19 (13/84) =======================
q-point: ( 0.46  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.846   (   0.579    8.581    0.000)    8.600
   1.902   (   0.356    0.796    0.000)    0.872
   2.192   (   8.833    4.715    0.000)   10.012
   2.241   (   2.009   -0.538    0.000)    2.080
   2.348   (   0.622   13.437    0.000)   13.451
   2.578   (  -0.181    8.128    0.000)    8.130
   2.905   (  -0.635   -0.514    0.000)    0.817
   2.911   (  -3.243   -0.340    0.000)    3.261
   2.989   (   1.671    3.951    0.000)    4.290
   4.344   (   5.446    5.212    0.000)    7.538
   4.739   (   1.788   -2.883    0.000)    3.393
   4.814   (   3.276    0.257    0.000)    3.286
   7.784   (  -4.250    0.378    0.000)    4.267
  10.030   (   6.556  -39.816    0.000)   40.353
  10.169   (   8.932  -25.512    0.000)   27.030
  11.023   (   5.199  -33.576    0.000)   33.976
  11.832   (   7.174   21.941    0.000)   23.084
  12.131   ( -14.478   -7.795    0.000)   16.443
  12.242   (   1.767    6.033    0.000)    6.287
  12.546   (  -5.557  -42.857    0.000)   43.216
  13.478   (   3.004    0.761    0.000)    3.099
  14.813   (  -3.955    3.959    0.000)    5.596
  14.855   (   1.994   20.365    0.000)   20.462
  15.522   (  -2.448   25.798    0.000)   25.914
======================= Grid point 28 (14/84) =======================
q-point: ( 0.15  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.993   (  14.150   14.150    0.000)   20.011
   1.068   (  13.429   13.429    0.000)   18.992
   1.884   (   1.579    1.579    0.000)    2.233
   2.075   (  11.345   11.345    0.000)   16.045
   2.130   (   6.067    6.067    0.000)    8.580
   2.281   (  10.205   10.205    0.000)   14.432
   2.487   (   6.269    6.269    0.000)    8.866
   2.681   (   7.343    7.343    0.000)   10.385
   3.072   (   4.488    4.488    0.000)    6.346
   3.653   (  11.898   11.898    0.000)   16.826
   3.918   (   5.486    5.486    0.000)    7.758
   4.927   (  -8.181   -8.181    0.000)   11.570
   8.140   (  -0.492   -0.492    0.000)    0.696
   8.221   (  -9.470   -9.470    0.000)   13.392
   8.866   (  -8.746   -8.746    0.000)   12.369
   9.973   (  13.272   13.272    0.000)   18.770
  11.920   (   6.311    6.311    0.000)    8.925
  12.119   (  15.606   15.606    0.000)   22.070
  12.629   ( -13.684  -13.684    0.000)   19.352
  13.082   ( -15.128  -15.128    0.000)   21.394
  13.413   (  -7.929   -7.929    0.000)   11.213
  14.456   (  16.585   16.585    0.000)   23.455
  15.492   (  -7.542   -7.542    0.000)   10.665
  16.217   (   5.983    5.983    0.000)    8.461
======================= Grid point 29 (15/84) =======================
q-point: ( 0.23  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.312   (  17.622   10.918    0.000)   20.730
   1.424   (  19.996    4.368    0.000)   20.468
   1.907   (   0.622    1.790    0.000)    1.895
   2.134   (   1.016   12.314    0.000)   12.356
   2.385   (  13.083    0.393    0.000)   13.089
   2.415   (   5.036   15.065    0.000)   15.884
   2.703   (  12.937   -0.509    0.000)   12.947
   2.753   (   1.543    2.035    0.000)    2.554
   3.168   (   3.024   11.251    0.000)   11.650
   3.849   (   9.827   16.182    0.000)   18.932
   4.144   (  15.371    1.168    0.000)   15.415
   4.754   (  -8.387   -3.754    0.000)    9.189
   8.164   (   2.833  -22.897    0.000)   23.072
   8.304   (  15.660  -20.435    0.000)   25.745
   8.642   ( -12.799  -11.873    0.000)   17.458
  10.193   (   9.143   -2.455    0.000)    9.467
  12.039   (   5.038   14.188    0.000)   15.056
  12.275   (  -8.205   -3.253    0.000)    8.826
  12.469   (   5.389   -3.600    0.000)    6.481
  12.765   ( -16.613  -20.030    0.000)   26.024
  13.277   (  -5.494   -8.969    0.000)   10.519
  14.775   (  15.052   19.852    0.000)   24.913
  15.323   (  -9.039   -4.720    0.000)   10.198
  16.261   (  -1.012   14.488    0.000)   14.523
======================= Grid point 30 (16/84) =======================
q-point: ( 0.31  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.677   (  18.774    9.068    0.000)   20.849
   1.789   (  15.605    4.459    0.000)   16.229
   1.913   (   0.181    1.849    0.000)    1.858
   2.158   (   1.613   10.424    0.000)   10.548
   2.520   (   5.561   15.594    0.000)   16.556
   2.590   (   6.964    0.652    0.000)    6.995
   2.797   (   3.577    0.174    0.000)    3.581
   2.920   (   7.709    0.036    0.000)    7.709
   3.171   (  -2.011   11.478    0.000)   11.653
   4.081   (  13.022   14.542    0.000)   19.521
   4.452   (  14.604    1.397    0.000)   14.671
   4.628   (  -3.608   -3.142    0.000)    4.784
   8.175   (  -8.629   -5.143    0.000)   10.046
   8.520   (   7.496  -42.133    0.000)   42.795
   8.705   (  23.452  -32.456    0.000)   40.042
  10.329   (   4.093  -16.570    0.000)   17.068
  12.028   ( -13.315  -27.366    0.000)   30.433
  12.088   (  -1.667   12.755    0.000)   12.864
  12.456   ( -12.522  -16.829    0.000)   20.976
  12.571   (   0.642    7.108    0.000)    7.137
  13.227   (   1.211   -8.184    0.000)    8.273
  15.039   (  11.401   21.964    0.000)   24.747
  15.149   (  -8.180   -0.770    0.000)    8.217
  16.202   (  -4.295   20.223    0.000)   20.674
======================= Grid point 31 (17/84) =======================
q-point: ( 0.38  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.917   (   0.463    1.660    0.000)    1.724
   2.008   (   6.226   10.088    0.000)   11.855
   2.042   (  17.284    8.333    0.000)   19.188
   2.211   (   3.784    6.231    0.000)    7.290
   2.625   (   4.639   17.225    0.000)   17.839
   2.679   (   3.071    3.531    0.000)    4.680
   2.873   (   2.854   -1.121    0.000)    3.066
   2.984   (  -1.077    1.746    0.000)    2.051
   3.125   (  -1.782   10.287    0.000)   10.440
   4.331   (  11.171   11.462    0.000)   16.005
   4.614   (   1.607   -4.624    0.000)    4.895
   4.699   (   9.874    1.697    0.000)   10.018
   7.928   ( -12.275   -3.548    0.000)   12.778
   8.831   (  19.606  -48.521    0.000)   52.332
   9.182   (  23.100  -36.555    0.000)   43.242
  10.333   (  -4.097  -30.668    0.000)   30.941
  11.771   ( -12.317  -39.595    0.000)   41.467
  11.934   ( -11.116  -13.279    0.000)   17.317
  12.340   (  -0.729   10.645    0.000)   10.670
  12.498   (  -6.182    4.499    0.000)    7.645
  13.304   (   5.077   -7.713    0.000)    9.234
  15.010   (  -5.575    4.259    0.000)    7.015
  15.221   (   6.882   23.122    0.000)   24.124
  16.114   (  -4.121   23.653    0.000)   24.009
======================= Grid point 32 (18/84) =======================
q-point: ( 0.46  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.924   (   0.236    1.430    0.000)    1.450
   2.070   (   1.063   13.125    0.000)   13.168
   2.277   (   2.013    4.337    0.000)    4.782
   2.326   (   9.388    8.511    0.000)   12.671
   2.691   (   1.692   19.184    0.000)   19.258
   2.728   (   1.662    6.222    0.000)    6.441
   2.891   (  -0.401   -0.867    0.000)    0.955
   2.915   (  -3.962    1.737    0.000)    4.326
   3.113   (   0.114    8.155    0.000)    8.156
   4.492   (   4.403    9.400    0.000)   10.380
   4.655   (   1.557   -5.465    0.000)    5.683
   4.833   (   3.452    1.879    0.000)    3.930
   7.751   (  -4.825   -5.829    0.000)    7.567
   9.187   (  12.746  -41.678    0.000)   43.584
   9.544   (  10.933  -35.433    0.000)   37.082
  10.196   (  -7.002  -46.404    0.000)   46.929
  11.570   (  -6.432  -48.437    0.000)   48.862
  11.765   (  -4.655  -20.123    0.000)   20.654
  12.345   (   0.417   19.663    0.000)   19.667
  12.394   (  -3.161    3.355    0.000)    4.610
  13.383   (   2.245   -9.365    0.000)    9.630
  14.934   (  -1.938    8.209    0.000)    8.434
  15.312   (   2.274   23.611    0.000)   23.720
  16.054   (  -1.642   25.273    0.000)   25.326
======================= Grid point 42 (19/84) =======================
q-point: ( 0.23  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.573   (  15.089   15.089    0.000)   21.340
   1.614   (  13.968   13.968    0.000)   19.753
   1.936   (   1.023    1.023    0.000)    1.446
   2.387   (   6.098    6.098    0.000)    8.624
   2.395   (   6.888    6.888    0.000)    9.741
   2.714   (  11.455   11.455    0.000)   16.200
   2.720   (   5.516    5.516    0.000)    7.800
   2.762   (   0.203    0.203    0.000)    0.287
   3.345   (   4.787    4.787    0.000)    6.770
   4.157   (  12.917   12.917    0.000)   18.268
   4.251   (  11.144   11.144    0.000)   15.760
   4.660   (  -5.500   -5.500    0.000)    7.779
   7.790   ( -12.075  -12.075    0.000)   17.077
   8.065   (  -3.265   -3.265    0.000)    4.617
   8.399   ( -14.493  -14.493    0.000)   20.497
  10.101   (  -5.517   -5.517    0.000)    7.802
  11.877   ( -23.627  -23.627    0.000)   33.413
  12.304   (  11.322   11.322    0.000)   16.012
  12.338   ( -22.371  -22.371    0.000)   31.637
  12.666   (   9.831    9.831    0.000)   13.903
  13.111   (  -7.469   -7.469    0.000)   10.563
  15.149   (  17.303   17.303    0.000)   24.469
  15.224   (  -4.686   -4.686    0.000)    6.627
  16.464   (   6.061    6.061    0.000)    8.571
======================= Grid point 43 (20/84) =======================
q-point: ( 0.31  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.897   (  17.023   12.863    0.000)   21.336
   1.924   (  15.895    9.498    0.000)   18.517
   1.948   (   0.474    1.635    0.000)    1.702
   2.397   (   0.302   12.571    0.000)   12.574
   2.611   (   8.889    1.379    0.000)    8.995
   2.786   (   2.482   -0.449    0.000)    2.522
   2.798   (   5.168   10.257    0.000)   11.486
   2.971   (  10.675    7.025    0.000)   12.779
   3.364   (  -1.766    6.749    0.000)    6.976
   4.378   (  10.759   14.704    0.000)   18.220
   4.528   (  14.076    6.562    0.000)   15.531
   4.555   (  -4.351   -4.171    0.000)    6.028
   7.708   (   2.302  -32.407    0.000)   32.489
   8.079   ( -15.392  -10.749    0.000)   18.774
   8.160   (  12.025  -20.756    0.000)   23.987
   9.979   (  -7.672  -16.839    0.000)   18.505
  11.419   ( -22.385  -30.986    0.000)   38.226
  11.871   ( -23.218  -25.589    0.000)   34.552
  12.500   (   7.880   12.047    0.000)   14.395
  12.720   (  -5.473    5.340    0.000)    7.646
  13.039   (   1.288  -11.041    0.000)   11.116
  15.152   (  -2.444    1.386    0.000)    2.810
  15.447   (  12.608   18.542    0.000)   22.423
  16.522   (   0.407   11.786    0.000)   11.793
======================= Grid point 44 (21/84) =======================
q-point: ( 0.38  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.952   (   0.268    1.759    0.000)    1.779
   2.197   (  10.741    9.423    0.000)   14.288
   2.240   (  16.623   11.391    0.000)   20.152
   2.402   (   0.273   11.771    0.000)   11.774
   2.742   (   4.898    3.059    0.000)    5.775
   2.850   (   3.025   -1.048    0.000)    3.201
   2.891   (   3.205    7.403    0.000)    8.067
   3.120   (   4.250   13.561    0.000)   14.211
   3.315   (  -2.491    7.913    0.000)    8.296
   4.509   (  -0.196   -5.770    0.000)    5.773
   4.577   (   8.708   12.814    0.000)   15.493
   4.768   (   9.631    5.376    0.000)   11.030
   7.617   ( -12.923  -34.076    0.000)   36.444
   8.081   (  17.449  -20.174    0.000)   26.673
   8.490   (  19.415  -31.331    0.000)   36.859
   9.744   ( -17.171  -26.105    0.000)   31.246
  10.996   ( -19.860  -36.707    0.000)   41.735
  11.465   ( -16.321  -29.635    0.000)   33.832
  12.541   ( -10.045   -0.839    0.000)   10.079
  12.621   (   4.401   13.217    0.000)   13.930
  13.091   (   2.553  -13.778    0.000)   14.012
  15.127   (  -0.294    7.370    0.000)    7.376
  15.646   (   7.399   19.039    0.000)   20.426
  16.505   (  -1.518   15.477    0.000)   15.552
======================= Grid point 45 (22/84) =======================
q-point: ( 0.46  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.957   (   0.157    1.760    0.000)    1.767
   2.335   (   3.149   12.715    0.000)   13.099
   2.410   (   0.486    9.002    0.000)    9.015
   2.519   (   9.690   10.731    0.000)   14.459
   2.822   (   3.125    3.419    0.000)    4.632
   2.876   (  -0.421   -0.361    0.000)    0.555
   2.898   (  -1.992    3.460    0.000)    3.993
   3.161   (   0.701   20.795    0.000)   20.807
   3.281   (  -0.798    8.037    0.000)    8.076
   4.526   (   1.048   -7.321    0.000)    7.395
   4.702   (   3.429   11.379    0.000)   11.885
   4.899   (   3.389    4.856    0.000)    5.922
   7.410   (  -5.996  -31.206    0.000)   31.777
   8.537   (  23.628  -21.813    0.000)   32.157
   8.830   (  11.394  -35.172    0.000)   36.971
   9.313   ( -22.813  -38.886    0.000)   45.084
  10.675   ( -10.255  -40.108    0.000)   41.398
  11.254   (  -5.064  -31.016    0.000)   31.427
  12.386   (  -4.507   -4.039    0.000)    6.052
  12.678   (   1.419   13.621    0.000)   13.695
  13.124   (   0.783  -16.473    0.000)   16.491
  15.128   (   0.202   10.889    0.000)   10.891
  15.742   (   2.407   19.173    0.000)   19.323
  16.479   (  -0.813   17.232    0.000)   17.251
======================= Grid point 56 (23/84) =======================
q-point: ( 0.31  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.971   (   0.843    0.843    0.000)    1.192
   2.154   (  12.833   12.833    0.000)   18.149
   2.181   (  15.257   15.257    0.000)   21.577
   2.612   (   5.308    5.308    0.000)    7.506
   2.650   (   5.699    5.699    0.000)    8.060
   2.795   (   1.492    1.492    0.000)    2.111
   2.929   (   4.952    4.952    0.000)    7.003
   3.189   (  12.502   12.502    0.000)   17.681
   3.434   (   0.897    0.897    0.000)    1.269
   4.468   (  -4.420   -4.420    0.000)    6.252
   4.637   (  10.930   10.930    0.000)   15.457
   4.706   (  10.681   10.681    0.000)   15.106
   7.191   ( -20.091  -20.091    0.000)   28.412
   7.851   ( -10.038  -10.038    0.000)   14.196
   7.919   (  -3.497   -3.497    0.000)    4.946
   9.683   ( -13.919  -13.919    0.000)   19.684
  10.838   ( -27.197  -27.197    0.000)   38.463
  11.317   ( -28.722  -28.722    0.000)   40.620
  12.706   (   8.340    8.340    0.000)   11.795
  12.762   (  -5.308   -5.308    0.000)    7.507
  12.793   (  -8.963   -8.963    0.000)   12.676
  15.206   (   3.824    3.824    0.000)    5.409
  15.762   (  13.072   13.072    0.000)   18.487
  16.681   (   4.767    4.767    0.000)    6.741
======================= Grid point 57 (24/84) =======================
q-point: ( 0.38  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.984   (   0.530    1.418    0.000)    1.514
   2.396   (  10.657   10.444    0.000)   14.921
   2.486   (  14.975   13.220    0.000)   19.976
   2.631   (  -0.220   10.502    0.000)   10.504
   2.800   (   5.938    2.367    0.000)    6.393
   2.838   (   2.510    0.120    0.000)    2.513
   3.003   (   1.897    5.488    0.000)    5.806
   3.421   (  10.133   14.802    0.000)   17.938
   3.429   (  -0.804    3.831    0.000)    3.915
   4.395   (  -2.474   -5.481    0.000)    6.013
   4.819   (   7.500   11.165    0.000)   13.450
   4.904   (   8.353    7.861    0.000)   11.470
   6.854   ( -14.790  -36.990    0.000)   39.837
   7.808   (   8.053  -10.678    0.000)   13.374
   7.993   (   9.436  -17.997    0.000)   20.320
   9.351   ( -21.603  -13.879    0.000)   25.677
  10.352   ( -21.371  -27.645    0.000)   34.942
  10.807   ( -20.523  -35.829    0.000)   41.290
  12.478   ( -14.620   -4.892    0.000)   15.417
  12.762   (  -1.224  -18.991    0.000)   19.031
  12.833   (   4.510    8.199    0.000)    9.358
  15.287   (   3.918    8.303    0.000)    9.181
  15.965   (   7.481   13.117    0.000)   15.100
  16.733   (   1.072    7.913    0.000)    7.986
======================= Grid point 58 (25/84) =======================
q-point: ( 0.46  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.991   (   0.217    1.683    0.000)    1.697
   2.555   (   5.201    9.817    0.000)   11.110
   2.617   (  -1.363   10.522    0.000)   10.610
   2.743   (   9.443   11.624    0.000)   14.976
   2.860   (   0.987    0.762    0.000)    1.247
   2.891   (   1.893    1.590    0.000)    2.472
   2.989   (  -2.859    6.101    0.000)    6.738
   3.416   (  -0.358    5.559    0.000)    5.571
   3.561   (   3.540   18.532    0.000)   18.867
   4.374   (  -0.054   -7.727    0.000)    7.727
   4.926   (   2.900   10.750    0.000)   11.134
   5.019   (   2.992    6.836    0.000)    7.462
   6.629   (  -6.522  -42.641    0.000)   43.137
   8.166   (  25.699  -18.585    0.000)   31.715
   8.193   (   7.493  -28.106    0.000)   29.088
   8.811   ( -30.622  -10.836    0.000)   32.482
  10.027   (  -9.704  -24.459    0.000)   26.314
  10.547   (  -5.955  -39.658    0.000)   40.103
  12.255   (  -6.429   -8.384    0.000)   10.565
  12.741   (  -0.630  -21.623    0.000)   21.632
  12.890   (   1.415    8.030    0.000)    8.154
  15.345   (   1.604   10.398    0.000)   10.521
  16.063   (   2.400   13.057    0.000)   13.276
  16.741   (   0.032    9.513    0.000)    9.513
======================= Grid point 70 (26/84) =======================
q-point: ( 0.38  0.38  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.007   (   0.958    0.958    0.000)    1.355
   2.592   (   8.758    8.758    0.000)   12.386
   2.750   (  12.858   12.858    0.000)   18.184
   2.786   (   3.040    3.040    0.000)    4.299
   2.833   (   3.437    3.437    0.000)    4.861
   2.865   (   1.835    1.835    0.000)    2.595
   3.104   (   3.737    3.737    0.000)    5.285
   3.482   (   1.845    1.845    0.000)    2.610
   3.689   (  12.097   12.097    0.000)   17.108
   4.299   (  -4.073   -4.073    0.000)    5.761
   5.007   (   7.604    7.604    0.000)   10.754
   5.055   (   6.733    6.733    0.000)    9.522
   6.218   ( -26.394  -26.394    0.000)   37.327
   7.709   (   2.655    2.655    0.000)    3.755
   7.773   (  -5.089   -5.089    0.000)    7.196
   9.112   ( -13.501  -13.501    0.000)   19.093
   9.916   ( -16.080  -16.080    0.000)   22.741
  10.114   ( -30.844  -30.844    0.000)   43.620
  12.364   ( -12.332  -12.332    0.000)   17.440
  12.384   ( -11.398  -11.398    0.000)   16.119
  12.954   (   4.229    4.229    0.000)    5.980
  15.436   (   6.274    6.274    0.000)    8.873
  16.168   (   7.386    7.386    0.000)   10.445
  16.838   (   3.104    3.104    0.000)    4.390
======================= Grid point 71 (27/84) =======================
q-point: ( 0.46  0.38  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 255
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.021   (   0.353    1.208    0.000)    1.258
   2.726   (   4.420    7.079    0.000)    8.346
   2.789   (  -1.007    6.205    0.000)    6.286
   2.881   (   0.451    1.591    0.000)    1.654
   2.915   (   2.212    0.735    0.000)    2.331
   2.968   (   8.002   10.594    0.000)   13.276
   3.125   (  -1.480    6.735    0.000)    6.895
   3.507   (   0.668    3.695    0.000)    3.755
   3.879   (   5.907   13.253    0.000)   14.510
   4.237   (  -1.981   -5.611    0.000)    5.950
   5.114   (   2.905    7.888    0.000)    8.406
   5.150   (   2.496    5.875    0.000)    6.383
   5.845   ( -10.164  -34.277    0.000)   35.752
   7.737   (  -0.041  -17.625    0.000)   17.625
   7.915   (  14.630   -3.500    0.000)   15.043
   8.736   ( -20.694   -0.235    0.000)   20.696
   9.711   (  -4.853   -7.972    0.000)    9.333
   9.717   (  -8.624  -42.593    0.000)   43.457
  12.085   (  -8.380   -7.940    0.000)   11.544
  12.302   (  -2.046  -21.301    0.000)   21.399
  13.007   (   1.285    3.994    0.000)    4.196
  15.523   (   2.302    7.237    0.000)    7.595
  16.263   (   2.349    7.300    0.000)    7.669
  16.873   (   0.724    4.199    0.000)    4.261
======================= Grid point 84 (28/84) =======================
q-point: ( 0.46  0.46  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 147
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.037   (   0.438    0.438    0.000)    0.619
   2.825   (   2.931    2.931    0.000)    4.146
   2.850   (   0.531    0.531    0.000)    0.751
   2.922   (   1.108    1.108    0.000)    1.566
   2.924   (   0.877    0.877    0.000)    1.240
   3.134   (   5.216    5.216    0.000)    7.376
   3.218   (   2.101    2.101    0.000)    2.971
   3.560   (   1.419    1.419    0.000)    2.007
   4.073   (   6.042    6.042    0.000)    8.544
   4.161   (  -2.360   -2.360    0.000)    3.337
   5.225   (   3.058    3.058    0.000)    4.325
   5.234   (   2.259    2.259    0.000)    3.195
   5.348   ( -13.833  -13.833    0.000)   19.563
   7.501   (  -6.021   -6.021    0.000)    8.516
   8.045   (  13.530   13.530    0.000)   19.134
   8.711   (  -5.642   -5.642    0.000)    7.979
   8.978   ( -24.398  -24.398    0.000)   34.504
   9.652   (  -0.006   -0.006    0.000)    0.009
  11.964   (  -6.030   -6.030    0.000)    8.528
  11.976   (  -7.047   -7.047    0.000)    9.966
  13.057   (   1.181    1.181    0.000)    1.671
  15.622   (   2.570    2.570    0.000)    3.634
  16.358   (   2.315    2.315    0.000)    3.274
  16.922   (   1.109    1.109    0.000)    1.569
======================= Grid point 169 (29/84) =======================
q-point: ( 0.00  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.386   (   0.000   -0.000   17.498)   17.498
   0.386   (  -0.000   -0.000   17.498)   17.498
   0.632   (   0.000    0.000   28.125)   28.125
   1.608   (   0.000    0.000  -11.190)   11.190
   1.608   (   0.000    0.000  -11.190)   11.190
   1.990   (   0.000    0.000    5.024)    5.024
   1.990   (  -0.000   -0.000    5.024)    5.024
   2.646   (   0.000    0.000  -12.525)   12.525
   3.115   (   0.000   -0.000    0.377)    0.377
   3.115   (  -0.000    0.000    0.377)    0.377
   3.749   (   0.000    0.000    4.144)    4.144
   5.644   (   0.000    0.000   -0.003)    0.003
   8.127   (  -0.000   -0.000    1.299)    1.299
   8.127   (   0.000   -0.000    1.299)    1.299
   8.592   (  -0.000    0.000   -1.654)    1.654
   8.592   (   0.000    0.000   -1.654)    1.654
   9.192   (   0.000   -0.000    0.807)    0.807
  13.055   (   0.000   -0.000   -0.417)    0.417
  13.684   (   0.000    0.000   -0.653)    0.653
  13.684   (   0.000   -0.000   -0.653)    0.653
  13.842   (   0.000    0.000    0.836)    0.836
  13.842   (   0.000    0.000    0.836)    0.836
  15.142   (  -0.000    0.000    2.754)    2.754
  16.039   (  -0.000   -0.000   -0.318)    0.318
======================= Grid point 170 (30/84) =======================
q-point: ( 0.08  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.508   (   7.291    0.000   14.492)   16.223
   0.698   (  23.233    0.000    9.538)   25.115
   1.003   (  32.047    0.000   11.019)   33.888
   1.648   (   3.736    0.000  -10.696)   11.330
   1.753   (   9.827    0.000   -8.867)   13.236
   2.000   (   1.082    0.000    4.873)    4.992
   2.031   (   4.565    0.000    0.892)    4.651
   2.449   (  -2.049    0.000    1.461)    2.517
   3.106   (  -0.881    0.000    0.345)    0.946
   3.192   (   4.967    0.000   -0.041)    4.967
   3.731   (   2.371    0.000    0.596)    2.445
   5.452   ( -14.354    0.000    2.980)   14.660
   8.254   (  12.592    0.000    1.245)   12.653
   8.720   (  12.688    0.000   -1.567)   12.784
   8.754   (  13.438    0.000    0.377)   13.444
   9.166   (  -1.320    0.000    0.756)    1.521
   9.804   (  84.739    0.000  -52.380)   99.621
  13.003   (  -8.481    0.000    4.011)    9.382
  13.181   ( -35.449    0.000   33.958)   49.090
  13.642   (  -4.171    0.000   -0.513)    4.203
  13.816   (  -2.664    0.000    0.694)    2.753
  13.826   (   8.301    0.000    7.037)   10.883
  15.260   (   5.527    0.000  -26.133)   26.711
  15.988   (  -4.678    0.000   -0.807)    4.747
======================= Grid point 171 (31/84) =======================
q-point: ( 0.15  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.734   (  15.395    0.000    3.697)   15.833
   1.149   (  20.784    0.000    5.390)   21.472
   1.649   (  28.714    0.000   -0.448)   28.718
   1.742   (   5.495    0.000   -9.086)   10.618
   1.883   (   5.621    0.000    3.070)    6.405
   2.033   (   2.236    0.000    4.188)    4.747
   2.225   (  13.488    0.000   -6.992)   15.193
   2.567   (  10.799    0.000    3.850)   11.465
   3.080   (  -1.676    0.000    0.259)    1.696
   3.286   (   4.639    0.000    0.663)    4.686
   3.865   (  10.978    0.000   -2.345)   11.226
   5.137   ( -16.129    0.000    2.998)   16.406
   8.609   (  22.551    0.000    1.093)   22.577
   8.850   (  -7.316    0.000    4.735)    8.715
   9.080   (  23.122    0.000   -1.321)   23.159
   9.328   (  19.381    0.000    0.067)   19.381
  10.946   (  38.458    0.000  -27.105)   47.050
  12.461   ( -27.314    0.000   30.622)   41.034
  12.971   (  -3.618    0.000   -1.268)    3.834
  13.518   (  -8.394    0.000   -0.145)    8.395
  13.738   (  -5.146    0.000    0.312)    5.155
  14.034   (  11.717    0.000    7.606)   13.969
  15.319   (   0.628    0.000  -21.158)   21.167
  15.858   (  -8.425    0.000   -2.902)    8.911
======================= Grid point 172 (32/84) =======================
q-point: ( 0.23  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.087   (  19.482    0.000   -0.965)   19.506
   1.490   (  13.057    0.000    2.864)   13.367
   1.859   (   1.491    0.000   -0.945)    1.766
   1.864   (   7.010    0.000   -5.929)    9.181
   2.086   (   2.966    0.000    2.800)    4.079
   2.252   (  21.817    0.000    0.206)   21.818
   2.459   (   8.660    0.000   -3.402)    9.304
   2.813   (  12.959    0.000    4.997)   13.889
   3.041   (  -2.284    0.000    0.151)    2.289
   3.417   (   9.843    0.000   -0.820)    9.877
   4.136   (  15.219    0.000   -2.573)   15.435
   4.856   ( -11.265    0.000    2.541)   11.548
   8.595   ( -16.494    0.000    4.920)   17.212
   9.117   (  27.802    0.000    0.874)   27.815
   9.612   (  29.822    0.000   -0.953)   29.837
   9.758   (  21.559    0.000    0.939)   21.579
  11.514   (  19.887    0.000   -8.266)   21.536
  12.037   ( -15.002    0.000   11.843)   19.113
  12.928   (  -0.128    0.000   -2.816)    2.819
  13.309   ( -12.629    0.000    0.325)   12.633
  13.614   (  -7.249    0.000   -0.199)    7.252
  14.267   (  11.393    0.000    6.671)   13.202
  15.284   (  -4.330    0.000   -8.345)    9.402
  15.662   ( -10.850    0.000   -5.280)   12.066
======================= Grid point 173 (33/84) =======================
q-point: ( 0.31  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.487   (  20.403    0.000   -1.451)   20.454
   1.666   (   5.289    0.000    0.476)    5.310
   1.877   (   0.708    0.000    0.687)    0.986
   2.022   (   8.358    0.000   -2.270)    8.661
   2.144   (   2.716    0.000    1.426)    3.068
   2.511   (   2.949    0.000   -3.751)    4.772
   2.622   (   9.870    0.000   -0.075)    9.870
   2.992   (  -2.552    0.000    0.059)    2.553
   3.029   (   7.583    0.000    6.260)    9.833
   3.686   (  16.181    0.000   -3.460)   16.547
   4.434   (  14.102    0.000   -1.977)   14.240
   4.716   (  -2.599    0.000    1.684)    3.097
   8.239   ( -18.361    0.000    2.057)   18.476
   9.674   (  27.300    0.000    0.645)   27.308
  10.143   (  16.805    0.000    0.689)   16.819
  10.236   (  32.411    0.000   -0.531)   32.415
  11.632   (  -6.639    0.000    3.539)    7.523
  11.980   (   7.737    0.000   -2.210)    8.046
  12.987   (   6.275    0.000   -3.099)    6.999
  13.017   ( -16.756    0.000    0.773)   16.773
  13.456   (  -8.475    0.000   -0.730)    8.506
  14.474   (   9.221    0.000    6.217)   11.121
  15.145   (  -9.485    0.000   -0.760)    9.515
  15.451   (  -9.886    0.000   -6.038)   11.584
======================= Grid point 174 (34/84) =======================
q-point: ( 0.38  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.732   (   1.904    0.000   -0.225)    1.918
   1.879   (  16.188    0.000    0.282)   16.190
   1.893   (   3.276    0.000    0.708)    3.352
   2.168   (   5.418    0.000   -0.537)    5.444
   2.192   (   2.563    0.000    0.435)    2.600
   2.513   (  -1.232    0.000    0.077)    1.235
   2.754   (   2.721    0.000   -4.833)    5.547
   2.943   (  -2.210    0.000    0.010)    2.210
   3.107   (   1.116    0.000    7.286)    7.371
   3.994   (  13.440    0.000   -6.391)   14.882
   4.672   (   9.500    0.000   -1.196)    9.575
   4.737   (   3.699    0.000    1.707)    4.074
   7.905   ( -14.383    0.000   -2.614)   14.619
  10.160   (  20.667    0.000    0.474)   20.673
  10.413   (  10.259    0.000    0.579)   10.276
  10.872   (  30.981    0.000   -0.176)   30.981
  11.490   (  -5.990    0.000   -1.002)    6.073
  12.113   (   4.586    0.000    0.619)    4.628
  12.649   ( -20.175    0.000    1.139)   20.208
  13.158   (   9.833    0.000   -1.685)    9.976
  13.291   (  -7.649    0.000   -1.188)    7.741
  14.625   (   5.793    0.000    6.499)    8.706
  14.934   ( -10.766    0.000    0.157)   10.768
  15.283   (  -6.804    0.000   -4.059)    7.923
======================= Grid point 175 (35/84) =======================
q-point: ( 0.46  0.00  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.754   (   0.499    0.000   -0.244)    0.556
   1.901   (   0.325    0.000    0.617)    0.698
   2.193   (  10.739    0.000    3.759)   11.378
   2.209   (   0.644    0.000    0.274)    0.699
   2.261   (   2.230    0.000   -0.324)    2.254
   2.482   (  -1.675    0.000    0.902)    1.902
   2.742   (  -2.365    0.000   -8.102)    8.441
   2.911   (  -0.932    0.000   -0.000)    0.932
   3.113   (   0.034    0.000    7.102)    7.102
   4.182   (   4.950    0.000   -8.485)    9.824
   4.800   (   3.366    0.000    0.736)    3.445
   4.810   (   2.387    0.000    0.954)    2.571
   7.702   (  -5.466    0.000   -6.090)    8.183
  10.453   (   7.929    0.000    0.401)    7.939
  10.548   (   3.371    0.000    0.902)    3.490
  11.410   (  -2.006    0.000   -3.437)    3.980
  11.416   (  20.830    0.000   -0.194)   20.831
  12.168   (   1.233    0.000    1.922)    2.284
  12.252   ( -16.989    0.000    1.565)   17.061
  13.176   (  -3.300    0.000   -1.472)    3.613
  13.320   (   5.016    0.000    0.085)    5.017
  14.702   (   2.030    0.000    6.595)    6.900
  14.764   (  -5.186    0.000   -0.687)    5.231
  15.189   (  -2.515    0.000   -2.457)    3.516
======================= Grid point 183 (36/84) =======================
q-point: ( 0.08  0.08  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.596   (   6.159    6.159   10.637)   13.748
   0.665   (  10.858   10.858   10.249)   18.462
   1.373   (  25.337   25.337    5.723)   36.286
   1.732   (   5.717    5.717  -10.897)   13.569
   1.795   (   4.528    4.528   -3.782)    7.437
   2.070   (   4.232    4.232    4.540)    7.512
   2.072   (   6.149    6.149   -3.867)    9.517
   2.451   (   2.900    2.900    2.134)    4.623
   3.089   (  -2.361   -2.361   -0.107)    3.341
   3.267   (   7.379    7.379    0.543)   10.450
   3.746   (   2.263    2.263   -0.565)    3.250
   5.344   ( -10.099  -10.099    3.072)   14.609
   8.146   (   0.788    0.788    1.586)    1.939
   8.507   (  -4.427   -4.427   -0.692)    6.299
   9.120   (  -3.685   -3.685    1.103)    5.326
   9.177   (  24.035   24.035    0.760)   33.999
  10.503   (  46.007   46.007  -40.690)   76.740
  12.876   ( -23.688  -23.688   40.674)   52.694
  12.963   (  -5.057   -5.057   -1.137)    7.241
  13.534   (  -7.517   -7.517   -0.213)   10.632
  13.716   (  -6.012   -6.012    0.717)    8.533
  13.971   (  10.891   10.891    7.249)   17.022
  15.395   (   7.508    7.508  -27.810)   29.768
  15.959   (  -2.520   -2.520   -2.591)    4.406
======================= Grid point 184 (37/84) =======================
q-point: ( 0.15  0.08  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.796   (  14.581    6.571    2.290)   16.156
   1.048   (  22.479   -3.867    6.101)   23.611
   1.779   (   8.646   10.903   -4.207)   14.537
   1.847   (   4.555    4.382   -2.195)    6.691
   1.940   (  16.210    7.353    0.744)   17.815
   2.150   (   4.646    9.015    3.792)   10.827
   2.282   (  13.634    6.834   -7.462)   16.978
   2.592   (  10.087    3.249    3.727)   11.234
   3.041   (  -1.952   -2.345   -0.431)    3.082
   3.414   (   6.503   11.050    0.635)   12.837
   3.862   (  10.167    0.006   -1.960)   10.355
   5.094   ( -13.881   -4.358    3.039)   14.863
   8.355   (  18.370  -17.304    1.324)   25.272
   8.544   (   6.168  -20.454    1.317)   21.404
   9.001   (  -6.539   -3.979    1.100)    7.733
   9.674   (  24.604   20.768    1.897)   32.253
  11.259   (  30.165   25.008  -20.091)   44.034
  12.414   ( -20.089   -4.026   25.682)   32.853
  12.857   (  -5.546  -10.649   -1.744)   12.133
  13.364   (  -9.913  -12.440    0.360)   15.911
  13.600   (  -5.969   -9.890    0.521)   11.563
  14.210   (  12.700   14.195    6.999)   20.292
  15.454   (  -1.539    6.859  -14.655)   16.254
  15.898   (  -3.671    5.525   -7.872)   10.294
======================= Grid point 185 (38/84) =======================
q-point: ( 0.23  0.08  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.137   (  19.018    5.239   -1.075)   19.756
   1.462   (  17.949   -2.406    3.342)   18.416
   1.866   (   2.030    3.771   -3.329)    5.425
   1.903   (   3.427    3.372   -1.617)    5.073
   2.176   (   5.009    9.115    2.128)   10.616
   2.376   (  16.346    4.991    1.592)   17.165
   2.523   (   8.751    8.622   -5.700)   13.543
   2.806   (  10.702   -0.677    4.697)   11.707
   3.020   (  -0.122   -0.007   -1.031)    1.038
   3.554   (   9.114   12.213   -0.442)   15.245
   4.128   (  15.384   -0.516   -2.025)   15.526
   4.839   ( -10.774   -1.737    2.505)   11.197
   8.542   ( -10.528   -7.194    4.736)   13.602
   8.811   (  26.279  -25.636    0.933)   36.724
   9.037   (  13.643  -30.496   -0.924)   33.421
  10.132   (  21.310   12.965    1.267)   24.976
  11.717   (  14.997   16.890   -4.175)   22.970
  12.119   (  -8.434    8.212    7.504)   13.960
  12.764   (  -3.329  -15.859   -3.146)   16.508
  13.134   ( -13.318  -15.193    0.824)   20.220
  13.475   (  -6.470  -10.285    0.204)   12.152
  14.460   (  12.166   16.494    5.936)   21.338
  15.356   (  -7.790    1.940   -3.154)    8.626
  15.808   (  -5.438   15.324   -9.376)   18.770
======================= Grid point 186 (39/84) =======================
q-point: ( 0.31  0.08  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.529   (  20.138    4.495   -1.242)   20.671
   1.725   (   8.230    4.694    1.292)    9.562
   1.888   (   0.705    1.283    0.376)    1.511
   2.006   (   6.844    0.593   -2.579)    7.338
   2.256   (   3.512   10.787    1.274)   11.415
   2.592   (   2.711    2.576    0.664)    3.799
   2.668   (   9.110    7.344   -3.523)   12.220
   2.950   (   1.660   -2.583    0.390)    3.095
   3.050   (   3.847    1.861    3.719)    5.665
   3.797   (  14.556   10.396   -2.765)   18.099
   4.431   (  14.337   -0.091   -1.544)   14.420
   4.698   (  -2.984   -1.713    1.675)    3.827
   8.239   ( -17.336   -0.321    2.048)   17.459
   9.353   (  27.329  -28.828    0.670)   39.729
   9.410   (  21.002  -44.247    0.197)   48.979
  10.515   (  17.172    0.741   -0.565)   17.197
  11.796   (  -2.932   15.542    3.815)   16.270
  12.128   (   5.301   11.008   -1.910)   12.366
  12.728   (  -1.479  -22.079   -3.329)   22.378
  12.840   ( -15.893  -16.193    1.113)   22.716
  13.368   (  -3.239   -6.510   -0.376)    7.281
  14.681   (   9.807   18.074    5.528)   21.293
  15.166   ( -10.728    0.196    0.254)   10.732
  15.688   (  -6.374   21.312   -6.684)   23.227
======================= Grid point 187 (40/84) =======================
q-point: ( 0.38  0.08  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.816   (   2.484    8.042   -0.102)    8.417
   1.898   (   1.774    1.079    0.746)    2.206
   1.924   (  17.365    4.035    0.341)   17.831
   2.145   (   6.112   -1.044   -1.389)    6.354
   2.314   (   2.244   11.883    0.192)   12.094
   2.604   (  -0.253    6.548    1.428)    6.707
   2.797   (   3.011    4.557   -4.652)    7.175
   2.931   (  -1.943   -0.561   -0.860)    2.198
   3.101   (   0.666   -0.360    6.107)    6.154
   4.077   (  12.330    7.961   -5.321)   15.612
   4.673   (   9.534    0.220   -0.881)    9.577
   4.711   (   3.457   -2.533    1.626)    4.584
   7.916   ( -14.101    0.786   -2.356)   14.318
   9.797   (  16.971  -45.250    0.662)   48.332
   9.856   (  21.962  -27.831    0.595)   35.457
  10.811   (  12.516  -16.637   -1.996)   20.915
  11.758   (   0.380   21.637    1.011)   21.664
  12.207   (   2.789    3.911    0.217)    4.809
  12.484   ( -17.958  -13.964    1.200)   22.780
  12.675   (  -6.128  -31.096   -1.381)   31.724
  13.382   (   4.254   -0.342   -0.874)    4.356
  14.838   (   5.631   18.276    5.449)   19.885
  14.957   (  -9.305    2.582    0.441)    9.667
  15.569   (  -5.267   24.268   -4.008)   25.154
======================= Grid point 188 (41/84) =======================
q-point: ( 0.46  0.08  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.843   (   0.565    8.559   -0.251)    8.581
   1.909   (   0.275    0.764    0.573)    0.994
   2.227   (   2.363   -0.779   -0.948)    2.662
   2.233   (  10.451    3.816    3.207)   11.579
   2.344   (   0.723   12.397   -0.399)   12.424
   2.588   (  -1.389    8.107    1.024)    8.288
   2.793   (  -1.914    5.295   -6.021)    8.244
   2.899   (  -0.989   -0.998   -0.706)    1.572
   3.103   (   0.039   -0.704    5.810)    5.853
   4.250   (   4.600    6.647   -7.117)   10.770
   4.775   (   2.374   -2.923    2.246)    4.385
   4.808   (   3.346    0.356   -0.480)    3.399
   7.715   (  -5.405    1.008   -5.587)    7.839
  10.037   (   6.433  -39.973    0.549)   40.491
  10.176   (   8.907  -25.466    0.515)   26.983
  10.993   (   5.066  -33.061   -2.539)   33.543
  11.836   (   7.207   22.625    0.325)   23.747
  12.144   ( -14.548   -8.322    1.120)   16.797
  12.259   (   2.174    5.125    1.332)    5.724
  12.539   (  -5.508  -41.966   -0.614)   42.330
  13.474   (   3.309    0.729   -0.361)    3.408
  14.801   (  -4.068    4.097   -0.934)    5.849
  14.929   (   1.876   19.663    6.355)   20.749
  15.493   (  -2.053   25.617   -2.708)   25.842
======================= Grid point 197 (42/84) =======================
q-point: ( 0.15  0.15  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.994   (  13.029   13.029    0.033)   18.426
   1.134   (  12.527   12.527    5.854)   18.658
   1.895   (   1.772    1.772    1.168)    2.764
   2.007   (   7.587    7.587   -9.024)   14.020
   2.098   (  10.582   10.582    1.899)   15.085
   2.333   (   8.828    8.828    2.712)   12.776
   2.460   (  10.040   10.040   -2.190)   14.367
   2.707   (   6.384    6.384    2.027)    9.253
   3.031   (   2.275    2.275   -3.288)    4.600
   3.661   (  11.904   11.904    0.721)   16.850
   3.901   (   5.779    5.779   -1.386)    8.290
   4.965   (  -8.412   -8.412    2.942)   12.254
   8.158   (  -0.784   -0.784    1.620)    1.963
   8.236   (  -9.000   -9.000    0.997)   12.766
   8.881   (  -8.723   -8.723    1.320)   12.407
  10.017   (  13.437   13.437    3.561)   19.333
  11.747   (  21.859   21.859   -8.747)   32.127
  12.344   (  -2.036   -2.036   12.929)   13.246
  12.600   ( -13.882  -13.882   -2.340)   19.772
  13.088   ( -14.974  -14.974    0.500)   21.183
  13.423   (  -7.883   -7.883    0.812)   11.177
  14.519   (  16.100   16.100    5.884)   23.516
  15.459   (  -4.851   -4.851   -3.465)    7.685
  16.057   (   7.766    7.766  -12.882)   16.928
======================= Grid point 198 (43/84) =======================
q-point: ( 0.23  0.15  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.300   (  17.444   11.122   -0.954)   20.710
   1.467   (  18.779    4.518    3.869)   19.698
   1.918   (   0.539    1.799    1.019)    2.137
   2.078   (   2.115   12.307   -5.262)   13.550
   2.346   (  11.080    3.618   -0.808)   11.684
   2.481   (   7.519   10.944    3.300)   13.682
   2.682   (   9.846    4.683   -1.264)   10.976
   2.788   (   3.841   -1.174    2.698)    4.838
   3.113   (   3.869    8.929   -4.922)   10.905
   3.848   (   8.907   16.507   -0.032)   18.757
   4.134   (  15.618    1.790   -0.845)   15.743
   4.784   (  -8.882   -3.864    2.395)    9.978
   8.190   (   3.988  -24.695    1.895)   25.087
   8.318   (  14.975  -20.173    1.286)   25.156
   8.651   ( -13.623  -10.913    0.760)   17.471
  10.220   (   7.521   -2.407    2.216)    8.201
  12.062   (   8.802   16.264    1.861)   18.587
  12.268   ( -11.180   -7.713   -0.306)   13.586
  12.438   (   3.522   -1.747   -2.649)    4.740
  12.776   ( -16.182  -20.141    0.902)   25.852
  13.285   (  -5.764   -8.813    0.656)   10.551
  14.828   (  14.643   19.226    4.769)   24.634
  15.321   (  -8.222   -4.064   -0.126)    9.173
  16.150   (   1.579   15.891   -9.441)   18.552
======================= Grid point 199 (44/84) =======================
q-point: ( 0.31  0.15  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.668   (  19.214    9.419   -0.702)   21.410
   1.811   (  14.674    4.748    2.039)   15.557
   1.923   (   0.160    1.798    0.906)    2.020
   2.119   (   2.289   10.459   -3.521)   11.271
   2.501   (   5.441   10.089   -0.592)   11.478
   2.638   (   6.612    5.684    3.194)    9.286
   2.782   (   3.809    1.509   -1.280)    4.292
   2.929   (   6.609    0.693    0.850)    6.700
   3.146   (  -0.174    8.258   -2.342)    8.586
   4.060   (  11.968   15.179   -1.648)   19.400
   4.445   (  14.734    1.975   -0.579)   14.877
   4.650   (  -3.899   -3.112    1.676)    5.262
   8.203   ( -11.537   -4.645    2.245)   12.638
   8.509   (   9.491  -43.209   -0.789)   44.246
   8.718   (  22.890  -32.333    0.966)   39.627
  10.322   (   2.620  -16.424   -0.531)   16.640
  12.006   ( -12.804  -27.043   -1.667)   29.967
  12.124   (  -3.320   11.108    3.075)   11.995
  12.467   ( -10.271  -13.959    0.758)   17.347
  12.551   (   0.626    5.927   -1.748)    6.211
  13.225   (   0.576   -8.128   -0.164)    8.150
  15.080   (   8.627   18.252    3.516)   20.492
  15.162   (  -5.523    2.186    1.181)    6.056
  16.135   (  -2.523   20.991   -5.762)   21.913
======================= Grid point 200 (45/84) =======================
q-point: ( 0.38  0.15  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.925   (   0.359    1.614    0.669)    1.784
   2.014   (   5.419   10.511    0.588)   11.840
   2.046   (  18.030    8.168    0.221)   19.795
   2.185   (   4.429    5.994   -2.220)    7.776
   2.583   (   2.943   11.643   -1.869)   12.154
   2.730   (   3.589    7.898    2.734)    9.096
   2.862   (   2.821    0.814   -0.937)    3.082
   2.982   (  -0.891    4.865   -0.245)    4.952
   3.129   (  -1.134    5.056    0.274)    5.189
   4.288   (  10.050   12.607   -3.251)   16.447
   4.637   (   2.095   -4.580    1.618)    5.290
   4.694   (  10.003    2.035   -0.363)   10.214
   7.907   ( -14.246   -3.251   -1.790)   14.721
   8.839   (  19.746  -48.438    0.634)   52.312
   9.190   (  22.960  -36.592    0.632)   43.204
  10.308   (  -4.424  -30.372   -2.023)   30.759
  11.759   ( -11.824  -39.076   -0.907)   40.836
  11.950   ( -11.411  -13.802    1.237)   17.951
  12.361   (  -1.113   10.927    1.749)   11.122
  12.497   (  -5.386    4.598   -0.101)    7.082
  13.296   (   5.109   -7.676   -0.683)    9.246
  15.003   (  -5.811    4.239   -0.433)    7.205
  15.279   (   6.898   22.156    4.756)   23.687
  16.073   (  -3.199   23.998   -3.543)   24.468
======================= Grid point 201 (46/84) =======================
q-point: ( 0.46  0.15  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.930   (   0.169    1.412    0.498)    1.507
   2.067   (   0.903   13.082   -0.264)   13.115
   2.257   (   2.159    4.223   -1.610)    5.009
   2.350   (  10.429    7.632    1.964)   13.072
   2.619   (   0.611   11.393   -3.712)   11.998
   2.783   (   1.061   13.465    2.461)   13.729
   2.872   (  -0.757    0.165   -1.620)    1.796
   2.931   (  -2.803    5.020    1.145)    5.863
   3.117   (  -0.035    3.811    0.324)    3.825
   4.429   (   3.772   10.807   -4.308)   12.230
   4.690   (   2.035   -5.383    2.043)    6.107
   4.830   (   3.516    2.074   -0.194)    4.087
   7.702   (  -5.598   -4.509   -4.051)    8.251
   9.191   (  12.474  -41.828    0.359)   43.650
   9.549   (  10.861  -35.550    0.463)   37.175
  10.173   (  -6.702  -46.227   -1.890)   46.748
  11.567   (  -6.139  -47.482   -0.154)   47.877
  11.778   (  -4.722  -20.518    0.971)   21.077
  12.358   (   0.115   19.343    1.072)   19.373
  12.402   (  -2.963    3.521    0.637)    4.646
  13.377   (   2.344   -9.491   -0.496)    9.789
  14.922   (  -2.087    8.164   -0.963)    8.481
  15.371   (   2.326   22.707    5.003)   23.367
  16.025   (  -1.366   25.423   -2.597)   25.592
======================= Grid point 211 (47/84) =======================
q-point: ( 0.23  0.23  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.567   (  15.375   15.375   -0.529)   21.751
   1.653   (  13.488   13.488    3.565)   19.405
   1.946   (   0.894    0.894    0.765)    1.478
   2.308   (   7.420    7.420   -6.996)   12.611
   2.397   (   5.791    5.791    0.796)    8.228
   2.742   (   2.578    2.578    0.320)    3.660
   2.743   (  11.377   11.377    1.924)   16.204
   2.780   (   2.918    2.918    3.526)    5.428
   3.272   (   5.220    5.220   -6.467)    9.814
   4.163   (  12.812   12.812    0.546)   18.127
   4.241   (  11.112   11.112   -0.714)   15.731
   4.687   (  -5.727   -5.727    2.126)    8.374
   7.796   ( -12.589  -12.589    0.392)   17.808
   8.078   (  -3.501   -3.501    1.236)    5.103
   8.414   ( -14.513  -14.513    1.225)   20.561
  10.118   (  -6.393   -6.393    1.395)    9.147
  11.854   ( -23.195  -23.195   -1.889)   32.857
  12.344   ( -22.518  -22.518    0.517)   31.849
  12.345   (  11.173   11.173    3.689)   16.226
  12.634   (   9.896    9.896   -2.913)   14.295
  13.120   (  -7.549   -7.549    0.754)   10.702
  15.179   (  11.730   11.730    2.638)   16.797
  15.241   (   0.465    0.465    1.493)    1.631
  16.385   (   7.712    7.712   -6.679)   12.788
======================= Grid point 212 (48/84) =======================
q-point: ( 0.31  0.23  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.896   (  17.271   13.246   -0.053)   21.766
   1.952   (  13.818    9.017    2.108)   16.634
   1.956   (   1.812    2.154    0.970)    2.977
   2.360   (   1.030   12.617   -3.610)   13.164
   2.579   (   8.834    1.284   -1.783)    9.104
   2.766   (   1.830    0.030   -1.467)    2.345
   2.853   (   5.174    8.388    4.280)   10.745
   2.993   (  10.363    7.815    1.714)   13.092
   3.311   (  -0.586    6.664   -4.766)    8.214
   4.369   (   9.901   15.308   -0.654)   18.243
   4.524   (  13.887    5.952   -0.147)   15.110
   4.579   (  -4.070   -3.676    1.683)    5.737
   7.705   (   1.647  -33.177   -0.328)   33.220
   8.085   ( -14.066  -10.724    0.759)   17.704
   8.175   (  10.378  -20.805    0.985)   23.270
   9.978   (  -8.467  -16.361   -0.043)   18.422
  11.406   ( -21.803  -30.437   -1.055)   37.455
  11.877   ( -23.194  -26.083    0.465)   34.907
  12.531   (   7.066   11.602    2.692)   13.848
  12.707   (  -3.926    5.924   -1.234)    7.213
  13.037   (   0.484  -11.020   -0.127)   11.032
  15.147   (  -2.695    1.202   -0.192)    2.957
  15.489   (  12.372   17.870    3.441)   22.005
  16.472   (   1.583   12.672   -4.245)   13.457
======================= Grid point 213 (49/84) =======================
q-point: ( 0.38  0.23  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.959   (   0.221    1.749    0.606)    1.864
   2.203   (  10.497    9.680    0.736)   14.298
   2.249   (  16.285   11.571    0.679)   19.989
   2.377   (   0.809   11.913   -2.275)   12.155
   2.711   (   4.721    3.694   -1.965)    6.308
   2.818   (   2.750   -0.074   -2.318)    3.598
   2.941   (   2.931    4.374    3.683)    6.425
   3.141   (   4.312   14.156    1.581)   14.882
   3.281   (  -1.703    7.860   -2.973)    8.574
   4.513   (   1.482   -0.855    0.400)    1.758
   4.572   (   6.664    8.693   -0.414)   10.961
   4.767   (   9.663    5.400   -0.024)   11.069
   7.604   ( -13.104  -33.160   -1.059)   35.670
   8.085   (  16.850  -21.058    0.340)   26.971
   8.496   (  19.203  -31.469    0.493)   36.869
   9.731   ( -17.649  -25.161   -1.085)   30.753
  10.995   ( -19.227  -36.132   -0.056)   40.929
  11.474   ( -16.097  -29.934    0.703)   33.995
  12.531   (  -9.034   -0.581   -0.550)    9.069
  12.653   (   3.244   12.900    2.380)   13.513
  13.083   (   2.503  -13.795   -0.704)   14.038
  15.116   (  -0.586    7.107   -0.811)    7.177
  15.686   (   7.376   18.228    3.341)   19.946
  16.473   (  -0.879   15.963   -2.742)   16.220
======================= Grid point 214 (50/84) =======================
q-point: ( 0.46  0.23  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.963   (   0.122    1.813    0.523)    1.891
   2.330   (   2.673   12.635   -0.340)   12.919
   2.394   (   0.834    9.435   -1.328)    9.564
   2.530   (   9.936   10.380    0.893)   14.397
   2.773   (   1.181    5.853   -3.134)    6.744
   2.849   (   0.336   -0.118   -2.013)    2.044
   2.952   (  -1.190    0.548    3.471)    3.710
   3.181   (   0.571   20.339    1.487)   20.401
   3.259   (  -0.458    8.485   -1.921)    8.712
   4.540   (   1.179   -6.258    1.203)    6.480
   4.685   (   3.260   11.158   -1.498)   11.721
   4.899   (   3.396    4.891    0.006)    5.955
   7.392   (  -6.153  -29.586   -1.425)   30.253
   8.534   (  23.575  -22.337   -0.247)   32.478
   8.832   (  11.258  -35.400    0.186)   37.148
   9.293   ( -23.130  -38.536   -1.636)   44.974
  10.685   (  -9.871  -39.512    0.841)   40.735
  11.266   (  -4.958  -30.925    1.003)   31.336
  12.383   (  -4.431   -4.410   -0.220)    6.255
  12.700   (   1.213   13.920    1.760)   14.083
  13.116   (   0.780  -16.598   -0.716)   16.632
  15.114   (   0.101   10.666   -1.191)   10.733
  15.782   (   2.414   18.231    3.384)   18.699
  16.455   (  -0.617   17.563   -2.077)   17.696
======================= Grid point 225 (51/84) =======================
q-point: ( 0.31  0.31  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.978   (   0.876    0.876    0.612)    1.382
   2.180   (  12.613   12.613    2.393)   17.998
   2.186   (  15.456   15.456    0.422)   21.863
   2.585   (   6.354    6.354   -5.397)   10.482
   2.608   (   5.027    5.027   -0.311)    7.116
   2.776   (   1.405    1.405   -1.662)    2.591
   2.964   (   4.351    4.351    3.108)    6.894
   3.218   (  12.363   12.363    2.190)   17.621
   3.390   (   1.702    1.702   -3.809)    4.506
   4.486   (  -4.456   -4.456    1.485)    6.474
   4.639   (  10.797   10.797    0.132)   15.269
   4.699   (  10.665   10.665   -0.593)   15.095
   7.183   ( -19.938  -19.938   -0.659)   28.204
   7.852   ( -10.776  -10.776    0.117)   15.240
   7.926   (  -3.682   -3.682    0.612)    5.243
   9.693   ( -13.442  -13.442    0.846)   19.029
  10.836   ( -26.613  -26.613   -0.130)   37.637
  11.316   ( -28.874  -28.874   -0.071)   40.834
  12.703   (   1.603    1.603   -0.411)    2.304
  12.778   (   1.471    1.471    1.398)    2.507
  12.796   (  -9.162   -9.162    0.299)   12.960
  15.195   (   3.478    3.478   -0.769)    4.978
  15.794   (  12.657   12.657    2.594)   18.087
  16.647   (   5.475    5.475   -2.852)    8.251
======================= Grid point 226 (52/84) =======================
q-point: ( 0.38  0.31  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.992   (   0.542    1.495    0.668)    1.725
   2.411   (   9.966   10.910    1.485)   14.851
   2.495   (  14.994   13.082    0.715)   19.912
   2.609   (   0.441   10.819   -1.996)   11.010
   2.769   (   6.569    2.043   -2.385)    7.281
   2.818   (   2.526    0.737   -1.640)    3.101
   3.027   (   1.398    5.053    2.026)    5.620
   3.399   (   0.007    4.562   -2.389)    5.150
   3.445   (   9.502   14.192    1.690)   17.163
   4.410   (  -2.655   -5.617    1.218)    6.331
   4.815   (   7.262   11.526   -0.388)   13.629
   4.901   (   8.535    7.712   -0.204)   11.505
   6.853   ( -14.402  -36.581   -0.095)   39.314
   7.795   (   7.517  -11.566   -1.054)   13.834
   7.997   (   9.162  -18.125    0.252)   20.310
   9.363   ( -21.955  -12.406    1.005)   25.237
  10.362   ( -20.782  -27.144    0.822)   34.196
  10.810   ( -19.983  -36.104    0.229)   41.266
  12.468   ( -14.450   -4.994   -0.785)   15.309
  12.750   (  -1.042  -18.434   -0.808)   18.482
  12.866   (   3.913    7.611    2.538)    8.926
  15.271   (   3.658    7.979   -1.335)    8.879
  15.991   (   7.251   12.453    2.085)   14.561
  16.711   (   1.499    8.366   -1.893)    8.708
======================= Grid point 227 (53/84) =======================
q-point: ( 0.46  0.31  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.999   (   0.207    1.791    0.664)    1.921
   2.556   (   4.547   10.463    0.163)   11.410
   2.604   (  -0.845   10.541   -0.965)   10.619
   2.752   (   9.205   11.712    0.750)   14.916
   2.848   (   1.559    1.718   -1.084)    2.561
   2.865   (   1.752    1.519   -2.079)    3.114
   3.008   (  -2.778    5.826    1.579)    6.645
   3.396   (  -0.132    5.685   -1.697)    5.935
   3.576   (   3.249   18.017    1.139)   18.343
   4.387   (  -0.099   -8.103    1.031)    8.169
   4.919   (   2.856   11.222   -0.582)   11.594
   5.018   (   3.031    6.776   -0.045)    7.423
   6.633   (  -6.404  -42.006    0.318)   42.493
   8.146   (  25.286  -19.703   -1.594)   32.096
   8.191   (   7.312  -28.247   -0.209)   29.179
   8.814   ( -30.727   -9.087    0.298)   32.044
  10.047   (  -9.334  -24.050    1.657)   25.851
  10.562   (  -5.550  -39.549    1.189)   39.954
  12.246   (  -6.412   -8.576   -0.692)   10.730
  12.730   (  -0.617  -21.654   -0.886)   21.681
  12.917   (   1.263    8.201    2.202)    8.585
  15.325   (   1.514   10.084   -1.657)   10.331
  16.085   (   2.326   12.235    1.861)   12.593
  16.724   (   0.171    9.867   -1.429)    9.971
======================= Grid point 239 (54/84) =======================
q-point: ( 0.38  0.38  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.017   (   1.011    1.011    0.784)    1.631
   2.610   (   8.554    8.554    1.688)   12.214
   2.747   (   9.690    9.690   -0.394)   13.709
   2.791   (   6.273    6.273    0.508)    8.886
   2.799   (   4.273    4.273   -2.916)    6.710
   2.855   (   2.408    2.408   -0.885)    3.518
   3.123   (   3.689    3.689    1.666)    5.477
   3.459   (   2.007    2.007   -1.971)    3.456
   3.702   (  11.527   11.527    1.077)   16.338
   4.309   (  -4.423   -4.423    0.793)    6.305
   5.005   (   7.598    7.598   -0.129)   10.746
   5.053   (   6.867    6.867   -0.199)    9.714
   6.225   ( -26.027  -26.027    0.540)   36.812
   7.681   (   2.074    2.074   -2.131)    3.626
   7.773   (  -5.250   -5.250   -0.160)    7.426
   9.148   ( -12.914  -12.914    2.938)   18.498
   9.934   ( -15.909  -15.909    1.524)   22.551
  10.116   ( -30.338  -30.338    0.086)   42.904
  12.347   ( -12.162  -12.162   -1.367)   17.254
  12.379   ( -11.597  -11.597   -0.400)   16.405
  12.985   (   4.155    4.155    2.515)    6.392
  15.414   (   6.023    6.023   -1.845)    8.715
  16.182   (   6.951    6.951    1.168)    9.899
  16.823   (   3.411    3.411   -1.233)    4.979
======================= Grid point 240 (55/84) =======================
q-point: ( 0.46  0.38  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.031   (   0.364    1.301    0.834)    1.587
   2.737   (   4.064    7.263    0.862)    8.367
   2.780   (  -0.686    6.650   -0.547)    6.707
   2.880   (   1.376    1.634   -0.235)    2.149
   2.892   (   2.076    1.709   -1.702)    3.182
   2.982   (   7.538   10.719    1.078)   13.148
   3.147   (  -1.129    6.964    1.877)    7.300
   3.486   (   0.691    3.566   -1.839)    4.071
   3.883   (   5.586   12.660    0.276)   13.841
   4.240   (  -2.204   -6.132    0.259)    6.522
   5.113   (   2.964    8.035   -0.086)    8.565
   5.149   (   2.518    5.947   -0.044)    6.458
   5.857   ( -10.004  -34.025    0.999)   35.480
   7.728   (  -0.161  -17.880   -0.700)   17.894
   7.883   (  14.092   -3.884   -2.622)   14.851
   8.769   ( -20.842    0.973    2.727)   21.042
   9.711   (  -6.437  -28.242    1.298)   28.996
   9.759   (  -6.008  -22.075    2.192)   22.983
  12.074   (  -8.339   -7.997   -0.903)   11.589
  12.291   (  -2.092  -21.287   -0.922)   21.409
  13.037   (   1.252    4.087    2.455)    4.929
  15.497   (   2.210    6.981   -2.149)    7.631
  16.272   (   2.198    6.742    0.691)    7.125
  16.862   (   0.831    4.456   -0.886)    4.619
======================= Grid point 253 (56/84) =======================
q-point: ( 0.46  0.46  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.049   (   0.473    0.473    0.937)    1.151
   2.835   (   2.735    2.735    0.710)    3.933
   2.853   (   0.929    0.929    0.174)    1.325
   2.916   (   1.561    1.561   -0.427)    2.248
   2.922   (   0.957    0.957   -0.072)    1.356
   3.151   (   5.274    5.274    1.373)    7.584
   3.244   (   2.266    2.266    2.280)    3.933
   3.536   (   1.383    1.383   -2.058)    2.840
   4.067   (   5.730    5.730   -0.560)    8.122
   4.155   (  -2.683   -2.683   -0.513)    3.829
   5.226   (   3.133    3.133    0.067)    4.431
   5.235   (   2.298    2.298    0.099)    3.251
   5.363   ( -13.779  -13.779    1.233)   19.525
   7.488   (  -6.152   -6.152   -1.058)    8.764
   8.005   (  12.657   12.657   -3.424)   18.225
   8.761   (  -5.510   -5.510    4.080)    8.795
   8.997   ( -23.224  -23.224    1.810)   32.893
   9.684   (   0.011    0.011    2.445)    2.446
  11.952   (  -5.947   -5.947   -0.988)    8.469
  11.966   (  -7.116   -7.116   -0.860)   10.101
  13.088   (   1.206    1.206    2.560)    3.076
  15.592   (   2.473    2.473   -2.453)    4.272
  16.358   (   2.118    2.118    0.060)    2.995
  16.915   (   1.202    1.202   -0.587)    1.798
======================= Grid point 338 (57/84) =======================
q-point: ( 0.00  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.701   (   0.000    0.000   10.400)   10.400
   0.701   (   0.000    0.000   10.400)   10.400
   1.119   (   0.000    0.000   15.036)   15.036
   1.381   (   0.000    0.000   -8.339)    8.339
   1.381   (   0.000    0.000   -8.339)    8.339
   2.061   (  -0.000   -0.000    1.676)    1.676
   2.061   (   0.000    0.000    1.676)    1.676
   2.368   (   0.000   -0.000  -10.289)   10.289
   3.122   (  -0.000    0.000    0.233)    0.233
   3.122   (   0.000    0.000    0.233)    0.233
   3.821   (   0.000    0.000    2.005)    2.005
   5.644   (   0.000    0.000   -0.001)    0.001
   8.154   (   0.000    0.000    0.939)    0.939
   8.154   (  -0.000   -0.000    0.939)    0.939
   8.557   (  -0.000    0.000   -1.158)    1.158
   8.557   (   0.000    0.000   -1.158)    1.158
   9.208   (   0.000    0.000    0.500)    0.500
  13.047   (   0.000    0.000   -0.259)    0.259
  13.672   (   0.000    0.000   -0.336)    0.336
  13.672   (  -0.000   -0.000   -0.336)    0.336
  13.858   (  -0.000    0.000    0.449)    0.449
  13.858   (  -0.000   -0.000    0.449)    0.449
  15.197   (   0.000    0.000    1.767)    1.767
  16.031   (   0.000    0.000   -0.265)    0.265
======================= Grid point 339 (58/84) =======================
q-point: ( 0.08  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.835   (  -1.902    0.000   14.714)   14.837
   0.914   (  18.054    0.000    8.504)   19.956
   1.112   (  17.385    0.000   -4.218)   17.889
   1.429   (   4.517    0.000   -8.292)    9.443
   1.628   (  18.453    0.000   -0.645)   18.464
   2.069   (   0.793    0.000    1.605)    1.790
   2.123   (  -0.604    0.000    6.640)    6.667
   2.332   (   6.288    0.000   -9.540)   11.426
   3.113   (  -0.944    0.000    0.213)    0.968
   3.192   (   5.827    0.000    0.068)    5.827
   3.756   (  -2.156    0.000    1.023)    2.386
   5.494   ( -13.282    0.000    0.971)   13.317
   8.280   (  12.467    0.000    0.892)   12.498
   8.687   (  12.877    0.000   -1.091)   12.923
   8.759   (  17.486    0.000    0.092)   17.486
   9.182   (  -1.584    0.000    0.557)    1.679
   9.214   (  78.591    0.000  -10.640)   79.308
  13.004   (  -3.863    0.000   -0.396)    3.883
  13.531   ( -17.622    0.000    4.254)   18.128
  13.633   (  -3.938    0.000   -0.274)    3.948
  13.829   (  -2.903    0.000    0.386)    2.928
  14.056   (  16.670    0.000    9.583)   19.228
  14.871   ( -21.070    0.000  -10.100)   23.366
  15.978   (  -5.313    0.000   -0.241)    5.318
======================= Grid point 340 (59/84) =======================
q-point: ( 0.15  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.823   (   5.886    0.000    3.449)    6.822
   1.297   (  18.425    0.000    7.536)   19.906
   1.546   (   7.174    0.000   -8.624)   11.218
   1.574   (  21.524    0.000   -4.048)   21.901
   2.005   (  15.972    0.000    6.051)   17.079
   2.092   (   1.515    0.000    1.372)    2.044
   2.131   (   7.789    0.000   -2.731)    8.254
   2.585   (  15.124    0.000   -0.920)   15.152
   3.085   (  -1.779    0.000    0.160)    1.786
   3.300   (   4.436    0.000    0.446)    4.458
   3.824   (   9.302    0.000   -1.163)    9.375
   5.189   ( -16.008    0.000    1.470)   16.076
   8.632   (  22.329    0.000    0.762)   22.342
   8.942   (  -5.802    0.000    2.867)    6.472
   9.053   (  23.465    0.000   -0.903)   23.483
   9.330   (  20.830    0.000    0.051)   20.830
  10.569   (  53.919    0.000   -8.710)   54.618
  12.887   ( -20.007    0.000    7.497)   21.366
  12.960   ( -22.107    0.000    2.103)   22.206
  13.515   (  -7.972    0.000   -0.084)    7.973
  13.744   (  -5.591    0.000    0.187)    5.594
  14.300   (   6.721    0.000   17.137)   18.408
  14.803   (  12.324    0.000  -24.164)   27.125
  15.824   ( -10.085    0.000   -0.639)   10.106
======================= Grid point 341 (60/84) =======================
q-point: ( 0.23  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.068   (  17.216    0.000   -0.615)   17.227
   1.570   (   8.193    0.000    4.018)    9.125
   1.739   (  12.592    0.000   -5.162)   13.609
   1.835   (   5.705    0.000   -0.921)    5.778
   2.127   (   1.958    0.000    0.981)    2.190
   2.255   (  12.471    0.000    0.178)   12.472
   2.420   (  14.459    0.000   -0.602)   14.471
   2.888   (  14.948    0.000    1.774)   15.053
   3.044   (  -2.389    0.000    0.093)    2.391
   3.404   (   7.536    0.000   -0.291)    7.541
   4.085   (  15.567    0.000   -1.560)   15.645
   4.903   ( -11.950    0.000    1.375)   12.029
   8.691   ( -17.736    0.000    2.976)   17.984
   9.135   (  27.540    0.000    0.587)   27.546
   9.593   (  30.257    0.000   -0.636)   30.263
   9.776   (  21.581    0.000    0.574)   21.589
  11.400   (  31.791    0.000   -2.598)   31.897
  12.229   ( -29.526    0.000    5.254)   29.990
  12.868   (  -0.769    0.000   -2.058)    2.197
  13.316   ( -12.138    0.000    0.207)   12.139
  13.610   (  -7.801    0.000   -0.129)    7.802
  14.441   (   8.413    0.000    7.770)   11.452
  15.049   (   8.815    0.000  -10.355)   13.598
  15.587   ( -13.358    0.000   -1.599)   13.454
======================= Grid point 342 (61/84) =======================
q-point: ( 0.31  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.457   (  21.143    0.000   -0.971)   21.165
   1.676   (   3.566    0.000    0.343)    3.583
   1.892   (   1.265    0.000    0.496)    1.359
   1.982   (  10.923    0.000   -1.120)   10.980
   2.166   (   1.886    0.000    0.549)    1.964
   2.440   (   5.637    0.000   -2.031)    5.992
   2.635   (   7.196    0.000    0.602)    7.221
   2.993   (  -2.625    0.000    0.036)    2.625
   3.137   (   8.899    0.000    3.122)    9.431
   3.616   (  12.982    0.000   -2.282)   13.181
   4.395   (  14.857    0.000   -1.176)   14.904
   4.747   (  -3.198    0.000    0.948)    3.335
   8.279   ( -22.590    0.000    1.195)   22.622
   9.687   (  27.075    0.000    0.417)   27.078
  10.154   (  16.362    0.000    0.286)   16.365
  10.225   (  32.840    0.000   -0.347)   32.842
  11.730   ( -13.592    0.000    4.663)   14.369
  11.913   (  12.796    0.000   -3.550)   13.280
  12.925   (   6.936    0.000   -1.956)    7.207
  13.033   ( -16.323    0.000    0.515)   16.331
  13.441   (  -9.007    0.000   -0.488)    9.020
  14.621   (   9.065    0.000    5.711)   10.714
  15.112   (  -3.072    0.000   -2.156)    3.753
  15.325   ( -11.504    0.000   -3.968)   12.169
======================= Grid point 343 (62/84) =======================
q-point: ( 0.38  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.727   (   1.746    0.000   -0.140)    1.752
   1.881   (  20.877    0.000   -0.049)   20.877
   1.907   (   0.653    0.000    0.455)    0.796
   2.159   (   6.795    0.000   -0.218)    6.798
   2.199   (   1.392    0.000    0.156)    1.401
   2.503   (   1.467    0.000   -0.688)    1.620
   2.704   (  -0.471    0.000   -0.405)    0.621
   2.944   (  -2.236    0.000    0.006)    2.236
   3.226   (   0.839    0.000    3.315)    3.419
   3.866   (  10.997    0.000   -4.162)   11.758
   4.650   (  10.332    0.000   -0.675)   10.354
   4.765   (   3.718    0.000    0.700)    3.783
   7.852   ( -18.842    0.000   -1.667)   18.916
  10.169   (  20.541    0.000    0.299)   20.543
  10.422   (  10.514    0.000    0.194)   10.515
  10.869   (  31.231    0.000   -0.118)   31.231
  11.475   (  -9.807    0.000   -0.356)    9.814
  12.125   (   6.760    0.000    0.347)    6.769
  12.672   ( -19.823    0.000    0.769)   19.838
  13.125   (  11.795    0.000   -1.030)   11.840
  13.267   (  -8.051    0.000   -0.797)    8.091
  14.779   (   6.359    0.000    6.177)    8.865
  14.936   ( -11.654    0.000   -0.025)   11.654
  15.178   (  -4.667    0.000   -4.595)    6.549
======================= Grid point 344 (63/84) =======================
q-point: ( 0.46  0.00  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.749   (   0.528    0.000   -0.151)    0.549
   1.913   (   0.199    0.000    0.379)    0.428
   2.215   (   0.492    0.000    0.234)    0.545
   2.253   (   2.368    0.000   -0.265)    2.383
   2.261   (  15.787    0.000    1.891)   15.900
   2.511   (  -0.983    0.000    1.408)    1.717
   2.616   (  -7.084    0.000   -3.329)    7.827
   2.911   (  -0.933    0.000   -0.000)    0.933
   3.220   (  -0.542    0.000    2.716)    2.769
   4.020   (   4.106    0.000   -4.986)    6.459
   4.791   (   3.709    0.000   -0.380)    3.728
   4.837   (   2.380    0.000    0.546)    2.441
   7.584   (  -7.287    0.000   -3.629)    8.141
  10.461   (   7.900    0.000    0.252)    7.904
  10.566   (   3.778    0.000    0.544)    3.817
  11.344   (  -3.359    0.000   -2.041)    3.930
  11.412   (  20.437    0.000   -0.131)   20.437
  12.205   (   1.803    0.000    1.118)    2.121
  12.284   ( -16.380    0.000    1.046)   16.413
  13.146   (  -3.459    0.000   -0.982)    3.595
  13.322   (   6.166    0.000    0.043)    6.166
  14.743   (  -5.758    0.000   -0.532)    5.783
  14.867   (   2.062    0.000    6.774)    7.081
  15.112   (  -1.808    0.000   -4.163)    4.538
======================= Grid point 352 (64/84) =======================
q-point: ( 0.08  0.08  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.810   (  -3.085   -3.085    7.352)    8.549
   0.889   (   8.194    8.194    8.349)   14.282
   1.424   (  22.016   22.016   -1.659)   31.180
   1.513   (   6.184    6.184   -7.955)   11.822
   1.733   (  11.123   11.123   -0.050)   15.730
   2.103   (  -0.973   -0.973    4.192)    4.412
   2.132   (   3.751    3.751    1.402)    5.487
   2.395   (   8.272    8.272   -5.090)   12.758
   3.085   (  -2.364   -2.364   -0.094)    3.345
   3.278   (   7.507    7.507    0.331)   10.622
   3.745   (   0.173    0.173    0.201)    0.317
   5.393   (  -9.659   -9.659    1.254)   13.717
   8.181   (   1.052    1.052    1.206)    1.915
   8.488   (  -3.705   -3.705   -0.769)    5.296
   9.144   (  -3.403   -3.403    0.809)    4.880
   9.192   (  25.852   25.852    0.454)   36.562
   9.980   (  55.470   55.470  -10.970)   79.210
  12.940   (  -5.700   -5.700   -0.708)    8.092
  13.378   ( -19.595  -19.595    9.197)   29.197
  13.530   (  -7.180   -7.180   -0.131)   10.155
  13.730   (  -6.026   -6.026    0.434)    8.533
  14.263   (  12.596   12.596   19.114)   26.128
  14.813   (   1.610    1.610  -24.591)   24.696
  15.933   (  -4.447   -4.447   -0.441)    6.305
======================= Grid point 353 (65/84) =======================
q-point: ( 0.15  0.08  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.847   (   8.275    3.431    1.779)    9.133
   1.192   (  19.023   -2.993    5.856)   20.128
   1.627   (   5.679    8.929   -6.716)   12.533
   1.805   (  11.911   12.074   -1.340)   17.014
   2.020   (  15.863    2.812    3.610)   16.510
   2.170   (   9.879    6.990   -2.668)   12.393
   2.204   (   4.690    7.885    1.050)    9.234
   2.626   (  13.201    4.222    0.116)   13.861
   3.032   (  -2.719   -3.830   -0.280)    4.706
   3.428   (   6.332   11.011    0.426)   12.709
   3.826   (   9.204    0.498   -1.047)    9.277
   5.148   ( -13.913   -4.267    1.522)   14.632
   8.385   (  17.890  -16.448    1.072)   24.326
   8.560   (   9.526  -21.967    0.202)   23.945
   9.027   (  -7.641   -3.455    1.009)    8.446
   9.712   (  24.973   23.100    1.159)   34.038
  10.978   (  42.349   32.609   -6.544)   53.848
  12.785   ( -19.884   -9.651    6.248)   22.968
  12.836   ( -19.031  -11.806    2.175)   22.501
  13.371   (  -9.449  -12.292    0.222)   15.506
  13.611   (  -6.235   -9.431    0.334)   11.310
  14.438   (   7.416   10.914   13.016)   18.535
  15.036   (  15.234   18.712  -19.515)   31.033
  15.807   (  -7.986   -1.220   -1.554)    8.227
======================= Grid point 354 (66/84) =======================
q-point: ( 0.23  0.08  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.116   (  17.552    5.171   -0.676)   18.310
   1.557   (  16.196   -2.695    5.102)   17.193
   1.737   (   5.684    6.721   -6.382)   10.873
   1.912   (   2.178    3.056    0.684)    3.815
   2.221   (   4.135    8.399    1.616)    9.500
   2.398   (  13.643   10.029    0.009)   16.933
   2.432   (  12.371    3.961   -2.073)   13.154
   2.885   (  12.055   -0.405    2.512)   12.320
   2.990   (  -0.909   -3.335   -1.494)    3.766
   3.548   (   7.356   12.608   -0.080)   14.597
   4.089   (  15.720    0.524   -1.208)   15.775
   4.885   ( -11.478   -1.788    1.364)   11.696
   8.636   ( -10.876   -7.807    2.975)   13.715
   8.832   (  25.890  -25.173    0.733)   36.117
   9.015   (  14.123  -31.045   -0.814)   34.116
  10.157   (  19.816   14.290    0.784)   24.444
  11.657   (  26.496   21.566   -1.353)   34.190
  12.250   ( -22.327    2.642    3.841)   22.808
  12.698   (  -4.220  -16.217   -2.257)   16.909
  13.149   ( -12.996  -14.842    0.401)   19.732
  13.479   (  -6.870   -9.643    0.134)   11.841
  14.610   (   9.821   14.622    6.267)   18.695
  15.246   (   3.912   11.969   -5.800)   13.864
  15.640   (  -7.744    8.174   -4.426)   12.098
======================= Grid point 355 (67/84) =======================
q-point: ( 0.31  0.08  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.504   (  20.881    4.933   -0.819)   21.471
   1.760   (   4.662    7.392    1.675)    8.898
   1.878   (   5.629   -0.252   -1.351)    5.795
   1.965   (   5.801   -0.462   -0.719)    5.864
   2.283   (   2.449   10.786    0.864)   11.094
   2.573   (   4.692   10.537   -1.506)   11.632
   2.643   (   8.131    1.273    0.397)    8.240
   2.948   (  -1.335   -3.124   -0.262)    3.408
   3.117   (   7.920   -1.892    1.888)    8.359
   3.744   (  11.823   11.699   -1.603)   16.710
   4.401   (  14.920    0.728   -0.897)   14.964
   4.730   (  -3.574   -1.734    0.937)    4.082
   8.278   ( -21.412   -0.317    1.194)   21.448
   9.368   (  27.143  -28.501    0.530)   39.362
   9.412   (  21.946  -44.817   -0.021)   49.902
  10.503   (  15.327    0.535   -0.358)   15.341
  11.892   (  -8.504   15.010    4.135)   17.740
  12.074   (   9.909   12.377   -2.709)   16.085
  12.664   (   0.310  -22.604   -2.029)   22.697
  12.858   ( -16.179  -15.770    0.498)   22.599
  13.361   (  -4.015   -6.116   -0.227)    7.319
  14.810   (   9.758   16.659    4.877)   19.913
  15.171   (  -9.656    0.938    0.126)    9.703
  15.537   (  -3.327   21.651   -5.570)   22.602
======================= Grid point 356 (68/84) =======================
q-point: ( 0.38  0.08  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.815   (   2.066    8.402   -0.034)    8.652
   1.914   (   1.844    1.032    0.546)    2.183
   1.926   (  18.988    4.123   -0.140)   19.431
   2.118   (   6.777   -2.046   -0.823)    7.126
   2.321   (   1.398   11.353    0.275)   11.442
   2.632   (   1.803    9.724    0.749)    9.918
   2.733   (   0.709    3.526   -1.226)    3.800
   2.914   (  -1.912   -1.930   -0.656)    2.795
   3.198   (   0.741   -2.723    2.572)    3.818
   3.974   (  10.225   10.128   -3.188)   14.741
   4.656   (  10.328    0.753   -0.514)   10.368
   4.738   (   3.343   -2.626    0.694)    4.307
   7.868   ( -18.262    1.276   -1.505)   18.369
   9.810   (  17.067  -45.349    0.417)   48.456
   9.867   (  21.884  -27.663    0.325)   35.274
  10.772   (  11.581  -17.456   -1.233)   20.985
  11.780   (  -1.462   23.028    0.751)   23.086
  12.211   (   4.280    2.841    0.119)    5.138
  12.511   ( -17.921  -13.823    0.999)   22.655
  12.644   (  -4.344  -30.141   -1.134)   30.473
  13.364   (   4.361    0.001   -0.555)    4.396
  14.942   (  -7.875    4.001    1.791)    9.012
  14.987   (   3.840   15.054    3.207)   15.864
  15.472   (  -3.297   24.993   -3.849)   25.502
======================= Grid point 357 (69/84) =======================
q-point: ( 0.46  0.08  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.838   (   0.543    8.560   -0.151)    8.579
   1.920   (   0.147    0.671    0.343)    0.768
   2.207   (   2.211   -1.888   -0.593)    2.967
   2.290   (  13.957    3.418    1.129)   14.414
   2.341   (   1.347   10.885    0.430)   10.977
   2.624   (  -4.917    7.032    2.213)    8.861
   2.694   (  -1.434    9.667   -3.289)   10.311
   2.881   (  -1.068   -2.376   -0.729)    2.705
   3.192   (  -0.473   -2.717    2.243)    3.555
   4.119   (   3.873    9.294   -3.926)   10.807
   4.797   (   3.692    0.582   -0.205)    3.743
   4.806   (   2.328   -2.910    0.497)    3.759
   7.607   (  -7.107    1.974   -3.343)    8.098
  10.048   (   6.267  -40.274    0.355)   40.760
  10.186   (   8.862  -25.378    0.289)   26.882
  10.941   (   4.763  -32.234   -1.623)   32.624
  11.843   (   7.153   23.784    0.260)   24.837
  12.169   ( -14.701   -8.894    0.923)   17.207
  12.284   (   3.037    3.308    0.720)    4.548
  12.525   (  -5.507  -40.493   -0.512)   40.869
  13.466   (   3.822    0.669   -0.239)    3.888
  14.783   (  -4.834    4.022   -0.567)    6.314
  15.074   (   2.366   18.161    5.542)   19.135
  15.420   (  -1.537   26.029   -3.328)   26.286
======================= Grid point 366 (70/84) =======================
q-point: ( 0.15  0.15  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.992   (  11.329   11.329   -0.133)   16.022
   1.275   (  10.891   10.891    5.651)   16.405
   1.820   (   8.639    8.639   -6.776)   13.970
   1.925   (   2.042    2.042    1.119)    3.098
   2.136   (   8.973    8.973    1.210)   12.747
   2.370   (   8.270    8.270    0.848)   11.726
   2.416   (  15.038   15.038   -1.461)   21.318
   2.747   (   7.222    7.222    1.346)   10.302
   2.966   (  -2.242   -2.242   -2.142)    3.827
   3.676   (  11.913   11.913    0.436)   16.854
   3.874   (   6.183    6.183   -0.832)    8.784
   5.018   (  -8.616   -8.616    1.549)   12.283
   8.198   (  -1.063   -1.063    1.639)    2.224
   8.245   (  -8.364   -8.364   -0.213)   11.830
   8.910   (  -8.756   -8.756    1.030)   12.426
  10.086   (  13.786   13.786    2.110)   19.611
  11.621   (  29.710   29.710   -2.986)   42.122
  12.551   ( -12.995  -12.995    5.477)   19.176
  12.555   ( -14.169  -14.169   -1.414)   20.087
  13.097   ( -14.749  -14.749    0.238)   20.860
  13.439   (  -7.823   -7.823    0.505)   11.075
  14.678   (  12.956   12.956    7.135)   19.662
  15.334   (   7.238    7.238   -6.936)   12.365
  15.823   (   4.368    4.368   -6.276)    8.807
======================= Grid point 367 (71/84) =======================
q-point: ( 0.23  0.15  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.282   (  17.244   11.441   -0.594)   20.703
   1.563   (  16.439    4.678    4.018)   17.558
   1.934   (   1.964    6.973   -1.338)    7.367
   1.957   (   3.060    6.152   -2.960)    7.481
   2.348   (  10.569    2.725    0.391)   10.921
   2.532   (   8.505   11.138    1.377)   14.081
   2.662   (   8.088    9.667   -0.531)   12.616
   2.864   (   4.513   -0.925    4.701)    6.582
   2.985   (   4.516    3.686   -6.457)    8.699
   3.849   (   7.701   16.831    0.061)   18.509
   4.118   (  16.020    2.793   -0.481)   16.269
   4.829   (  -9.548   -4.011    1.320)   10.440
   8.220   (   5.766  -27.254    0.721)   27.866
   8.349   (  14.471  -19.994    1.140)   24.707
   8.668   ( -15.509   -9.438    0.578)   18.164
  10.264   (   4.946   -2.310    1.360)    5.626
  12.094   (  13.730   18.183    0.761)   22.797
  12.281   ( -10.069  -15.051    1.544)   18.175
  12.373   (  -4.036    2.722   -2.920)    5.677
  12.792   ( -15.682  -20.261    0.454)   25.625
  13.298   (  -6.184   -8.570    0.407)   10.576
  14.941   (  13.214   17.538    4.321)   22.380
  15.318   (  -6.226   -2.407   -0.191)    6.678
  15.947   (   5.721   18.176   -7.035)   20.312
======================= Grid point 368 (72/84) =======================
q-point: ( 0.31  0.15  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.654   (  19.837   10.097   -0.437)   22.263
   1.866   (  12.322    4.328    2.468)   13.291
   1.948   (   0.957    2.337    1.630)    3.006
   2.026   (   3.514   10.345   -4.554)   11.837
   2.503   (   5.041    7.561    0.406)    9.097
   2.700   (   6.567    7.481    2.476)   10.258
   2.748   (   3.671    3.541   -1.772)    5.399
   2.954   (   3.249    5.146    1.310)    6.225
   3.086   (   4.561    0.197   -2.528)    5.218
   4.031   (  10.328   16.215   -0.817)   19.242
   4.434   (  14.941    2.854   -0.321)   15.214
   4.681   (  -4.455   -3.124    0.917)    5.518
   8.245   ( -16.922   -3.875    1.269)   17.407
   8.497   (  12.620  -44.718   -0.244)   46.465
   8.735   (  22.697  -32.270    0.430)   39.455
  10.312   (   0.262  -16.199   -0.323)   16.204
  11.975   ( -12.173  -27.183   -0.962)   29.800
  12.192   (  -7.620    7.660    2.494)   11.088
  12.468   (   3.707    4.556   -0.759)    5.923
  12.520   (  -7.891   -8.480   -0.553)   11.597
  13.222   (  -0.457   -8.022   -0.073)    8.035
  15.150   (  -5.651    2.226    1.851)    6.350
  15.197   (   8.226   16.511    1.922)   18.546
  16.009   (   0.674   22.732   -4.409)   23.165
======================= Grid point 369 (73/84) =======================
q-point: ( 0.38  0.15  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.938   (   0.185    1.526    0.420)    1.594
   2.031   (   4.199   12.321    0.683)   13.035
   2.050   (  17.775    8.079    0.307)   19.527
   2.134   (   6.832    4.435   -2.023)    8.393
   2.566   (   1.632   10.324   -0.031)   10.452
   2.765   (   2.455    2.345    0.497)    3.431
   2.847   (   3.295    8.480   -0.218)    9.100
   2.976   (  -0.595    9.576   -0.254)    9.598
   3.133   (   0.037   -2.848    0.082)    2.849
   4.229   (   8.871   14.720   -1.672)   17.267
   4.664   (   2.123   -4.728    0.680)    5.227
   4.687   (  10.203    2.594   -0.201)   10.530
   7.871   ( -17.279   -2.712   -1.137)   17.528
   8.852   (  19.830  -48.325    0.410)   52.237
   9.202   (  22.781  -36.642    0.355)   43.148
  10.268   (  -4.976  -29.910   -1.233)   30.346
  11.743   ( -10.922  -38.626   -0.453)   40.143
  11.973   ( -12.060  -14.471    0.723)   18.852
  12.398   (  -1.943   11.914    1.316)   12.142
  12.492   (  -3.769    4.420   -0.351)    5.819
  13.282   (   5.128   -7.619   -0.434)    9.194
  14.997   (  -6.750    3.777   -0.100)    7.735
  15.379   (   7.441   20.596    3.482)   22.174
  15.993   (  -1.691   24.991   -2.890)   25.214
======================= Grid point 370 (74/84) =======================
q-point: ( 0.46  0.15  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.940   (   0.054    1.371    0.307)    1.406
   2.061   (   0.617   13.074   -0.164)   13.090
   2.226   (   2.349    4.112   -0.904)    4.821
   2.386   (  13.116    6.428    0.952)   14.637
   2.569   (  -1.801    9.612   -0.851)    9.816
   2.811   (   1.584    3.320   -0.083)    3.679
   2.841   (  -2.775   11.293   -0.532)   11.641
   2.951   (  -1.204   11.774    0.606)   11.851
   3.123   (  -0.285   -3.489    0.159)    3.504
   4.356   (   3.420   13.780   -2.036)   14.343
   4.718   (   2.077   -5.781    0.574)    6.170
   4.827   (   3.611    2.386   -0.106)    4.330
   7.622   (  -6.776   -2.446   -2.446)    7.608
   9.198   (  12.061  -42.079    0.216)   43.774
   9.558   (  10.741  -35.723    0.263)   37.303
  10.135   (  -6.300  -45.876   -1.201)   46.322
  11.567   (  -5.584  -46.890    0.095)   47.221
  11.795   (  -4.907  -20.301    0.462)   20.891
  12.381   (  -0.280   20.124    0.837)   20.144
  12.412   (  -2.709    2.544    0.211)    3.722
  13.367   (   2.513   -9.689   -0.319)   10.014
  14.904   (  -2.449    8.065   -0.589)    8.449
  15.480   (   2.565   20.788    3.852)   21.297
  15.963   (  -0.972   26.155   -2.384)   26.282
======================= Grid point 380 (75/84) =======================
q-point: ( 0.23  0.23  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.556   (  15.867   15.867   -0.328)   22.441
   1.742   (  12.507   12.507    3.656)   18.061
   1.961   (   0.586    0.586    0.489)    0.962
   2.160   (   8.239    8.239   -5.335)   12.815
   2.413   (   5.468    5.468    0.535)    7.752
   2.736   (  -0.107   -0.107   -0.350)    0.382
   2.773   (  11.637   11.637    0.775)   16.475
   2.905   (   4.574    4.574    7.332)    9.778
   3.103   (   6.251    6.251   -8.591)   12.327
   4.174   (  12.646   12.646    0.334)   17.887
   4.228   (  11.157   11.157   -0.378)   15.783
   4.726   (  -6.108   -6.108    1.170)    8.718
   7.801   ( -13.219  -13.219    0.072)   18.694
   8.106   (  -3.914   -3.914    0.962)    5.618
   8.438   ( -14.730  -14.730    0.738)   20.844
  10.146   (  -7.810   -7.810    0.862)   11.078
  11.817   ( -22.529  -22.529   -1.139)   31.881
  12.354   ( -22.704  -22.704    0.278)   32.109
  12.443   (  10.832   10.832    5.184)   16.172
  12.551   (  10.101   10.101   -4.736)   15.049
  13.135   (  -7.675   -7.675    0.461)   10.865
  15.230   (  -2.998   -2.998    1.373)    4.456
  15.283   (  14.045   14.045    1.977)   19.961
  16.243   (  10.834   10.834   -4.889)   16.083
======================= Grid point 381 (76/84) =======================
q-point: ( 0.31  0.23  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.896   (  17.905   13.934    0.018)   22.688
   1.966   (   0.732    1.877    0.344)    2.044
   2.017   (  13.213    9.217    2.787)   16.350
   2.266   (   3.387   12.034   -3.880)   13.090
   2.561   (   8.017    0.959   -0.106)    8.075
   2.735   (   0.713   -0.601   -1.054)    1.408
   2.950   (   5.575   10.400    3.672)   12.358
   3.034   (   8.280    5.535    2.342)   10.232
   3.188   (   2.302    7.462   -5.724)    9.682
   4.358   (   9.039   16.216   -0.275)   18.567
   4.523   (  14.449    6.158   -0.017)   15.707
   4.608   (  -4.944   -4.125    0.833)    6.492
   7.697   (   0.736  -33.962   -0.274)   33.971
   8.106   ( -15.144  -11.002    0.844)   18.737
   8.189   (  10.766  -21.003    0.297)   23.603
   9.978   (  -9.754  -15.586   -0.029)   18.387
  11.385   ( -20.871  -29.598   -0.636)   36.222
  11.886   ( -23.122  -26.861    0.282)   35.443
  12.598   (   5.005   10.616    3.259)   12.180
  12.667   (  -0.582    7.154   -2.434)    7.580
  13.035   (  -0.901  -10.977   -0.039)   11.014
  15.148   (  -3.489    0.560    0.141)    3.536
  15.560   (  12.150   16.804    2.317)   20.866
  16.382   (   3.700   14.421   -3.047)   15.196
======================= Grid point 382 (77/84) =======================
q-point: ( 0.38  0.23  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.971   (   0.142    1.735    0.380)    1.782
   2.232   (  10.867   10.235    1.963)   15.056
   2.259   (  14.638   11.387    0.017)   18.546
   2.319   (   2.264   12.499   -2.627)   12.971
   2.688   (   4.342    3.226   -0.143)    5.411
   2.768   (   2.276   -0.031   -1.890)    2.958
   3.022   (   1.791    1.783    3.203)    4.080
   3.169   (   4.392   14.871    0.809)   15.527
   3.208   (   0.165    8.785   -3.109)    9.321
   4.524   (   6.004   11.222    0.464)   12.736
   4.562   (   1.529   -2.173   -0.469)    2.698
   4.767   (   9.718    5.439   -0.010)   11.137
   7.584   ( -13.420  -31.723   -0.645)   34.451
   8.091   (  15.868  -22.461    0.194)   27.501
   8.506   (  18.903  -31.642    0.296)   36.860
   9.709   ( -18.413  -23.648   -0.669)   29.979
  10.994   ( -18.214  -35.261   -0.023)   39.688
  11.488   ( -15.737  -30.378    0.438)   34.215
  12.524   (  -8.871   -1.199   -0.146)    8.953
  12.696   (   2.875   13.535    1.297)   13.898
  13.069   (   2.427  -13.902   -0.472)   14.120
  15.101   (  -1.246    6.598   -0.416)    6.728
  15.755   (   7.466   16.705    2.310)   18.442
  16.414   (   0.183   17.027   -2.020)   17.147
======================= Grid point 383 (78/84) =======================
q-point: ( 0.46  0.23  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.973   (   0.059    1.908    0.321)    1.936
   2.326   (   1.180   12.865   -0.023)   12.919
   2.368   (   2.111    9.793   -0.845)   10.054
   2.547   (  10.929    9.778    0.470)   14.673
   2.725   (  -0.887    6.753   -1.061)    6.893
   2.805   (   1.091    0.396   -1.586)    1.965
   3.020   (  -1.325   -2.361    2.308)    3.558
   3.200   (  -0.067   13.144   -0.475)   13.152
   3.221   (   0.771   16.826   -0.574)   16.853
   4.572   (   1.377   -6.271    1.530)    6.600
   4.647   (   2.988   12.536   -1.724)   13.002
   4.899   (   3.409    4.941    0.009)    6.003
   7.364   (  -6.408  -26.994   -0.878)   27.758
   8.529   (  23.546  -23.232   -0.157)   33.078
   8.836   (  10.999  -35.725    0.115)   37.380
   9.261   ( -23.627  -37.902   -1.000)   44.675
  10.702   (  -9.270  -38.654    0.516)   39.753
  11.286   (  -4.788  -30.728    0.621)   31.105
  12.379   (  -4.359   -5.040   -0.087)    6.664
  12.734   (   0.926   14.469    1.058)   14.537
  13.101   (   0.797  -16.814   -0.456)   16.839
  15.090   (  -0.129   10.285   -0.742)   10.313
  15.854   (   2.493   16.409    2.418)   16.773
  16.409   (  -0.318   18.406   -1.620)   18.480
======================= Grid point 394 (79/84) =======================
q-point: ( 0.31  0.31  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.990   (   0.937    0.937    0.365)    1.374
   2.195   (  15.823   15.823    0.297)   22.379
   2.243   (  12.265   12.265    2.875)   17.582
   2.467   (   7.200    7.200   -4.465)   11.119
   2.604   (   4.453    4.453   -0.038)    6.297
   2.741   (   1.321    1.321   -1.163)    2.200
   3.039   (   2.673    2.673    3.079)    4.875
   3.257   (  12.336   12.336    1.080)   17.479
   3.301   (   3.710    3.710   -3.519)    6.317
   4.516   (  -4.725   -4.725    0.987)    6.755
   4.641   (  10.586   10.586    0.083)   14.971
   4.686   (  10.862   10.862   -0.472)   15.369
   7.170   ( -19.718  -19.718   -0.395)   27.888
   7.855   ( -11.930  -11.930    0.072)   16.871
   7.938   (  -3.993   -3.993    0.365)    5.659
   9.710   ( -12.668  -12.668    0.525)   17.922
  10.833   ( -25.669  -25.669   -0.081)   36.302
  11.314   ( -29.121  -29.121   -0.042)   41.184
  12.693   (  -3.661   -3.661   -0.328)    5.188
  12.802   (  -9.479   -9.479    0.184)   13.407
  12.807   (   6.829    6.829    0.921)    9.701
  15.182   (   2.770    2.770   -0.343)    3.933
  15.846   (  11.953   11.953    1.688)   16.988
  16.587   (   6.774    6.774   -1.967)    9.779
======================= Grid point 395 (80/84) =======================
q-point: ( 0.38  0.31  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.005   (   0.576    1.608    0.413)    1.757
   2.456   (   8.638   11.687    2.528)   14.751
   2.505   (  13.262   12.952   -0.410)   18.542
   2.557   (   4.286   11.086   -2.115)   12.072
   2.728   (   6.344    1.396   -1.081)    6.585
   2.784   (   2.789    1.818   -1.100)    3.506
   3.072   (   0.174    4.213    1.635)    4.523
   3.346   (   1.250    5.046   -1.888)    5.531
   3.476   (   9.135   14.209    0.900)   16.916
   4.434   (  -2.954   -5.929    0.759)    6.667
   4.806   (   6.913   12.226   -0.278)   14.048
   4.897   (   8.829    7.464   -0.108)   11.561
   6.851   ( -13.772  -35.942   -0.066)   38.490
   7.775   (   6.628  -12.798   -0.615)   14.425
   8.001   (   8.732  -18.508    0.122)   20.464
   9.383   ( -22.517  -10.019    0.620)   24.653
  10.378   ( -19.853  -26.300    0.506)   32.955
  10.814   ( -19.082  -36.590    0.151)   41.267
  12.453   ( -14.260   -5.316   -0.437)   15.225
  12.736   (  -1.458  -18.908   -0.370)   18.967
  12.913   (   3.767    8.204    1.394)    9.135
  15.245   (   3.123    7.389   -0.789)    8.061
  16.033   (   6.922   11.271    1.364)   13.297
  16.671   (   2.230    9.247   -1.294)    9.599
======================= Grid point 396 (81/84) =======================
q-point: ( 0.46  0.31  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.013   (   0.194    1.967    0.416)    2.020
   2.561   (   1.965   10.847    0.159)   11.024
   2.587   (   1.600   11.375   -0.461)   11.497
   2.766   (   8.840   11.934    0.393)   14.857
   2.818   (   2.418    2.149   -1.196)    3.448
   2.832   (   1.460    2.112   -0.770)    2.680
   3.040   (  -2.784    5.447    1.068)    6.210
   3.361   (   0.327    5.817   -1.141)    5.936
   3.596   (   2.847   17.357    0.591)   17.599
   4.407   (  -0.172   -8.783    0.641)    8.808
   4.907   (   2.799   12.069   -0.374)   12.395
   5.017   (   3.086    6.696   -0.017)    7.373
   6.639   (  -6.210  -40.985    0.194)   41.453
   8.114   (  24.607  -21.146   -1.005)   32.460
   8.187   (   6.920  -28.804   -0.138)   29.624
   8.821   ( -30.734   -6.300    0.240)   31.374
  10.079   (  -8.794  -23.367    0.982)   24.986
  10.586   (  -4.892  -39.372    0.759)   39.682
  12.233   (  -6.395   -8.929   -0.399)   10.990
  12.713   (  -0.701  -21.785   -0.549)   21.803
  12.961   (   1.139    8.598    1.348)    8.777
  15.292   (   1.324    9.549   -1.048)    9.697
  16.123   (   2.241   10.752    1.257)   11.055
  16.694   (   0.396   10.598   -0.996)   10.652
======================= Grid point 408 (82/84) =======================
q-point: ( 0.38  0.38  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.032   (   1.101    1.101    0.482)    1.630
   2.659   (   7.836    7.836    2.817)   11.434
   2.727   (   5.939    5.939   -3.495)    9.097
   2.737   (   6.633    6.633   -0.333)    9.386
   2.805   (   9.563    9.563    0.642)   13.540
   2.835   (   3.210    3.210   -0.744)    4.600
   3.160   (   3.559    3.559    1.306)    5.200
   3.416   (   2.294    2.294   -1.527)    3.586
   3.722   (  10.669   10.669    0.566)   15.100
   4.324   (  -4.982   -4.982    0.478)    7.062
   5.003   (   7.610    7.610   -0.055)   10.763
   5.049   (   7.096    7.096   -0.120)   10.037
   6.236   ( -25.424  -25.424    0.320)   35.956
   7.642   (   1.504    1.504   -1.166)    2.426
   7.767   (  -5.921   -5.921   -0.247)    8.377
   9.206   ( -11.952  -11.952    1.829)   17.002
   9.966   ( -15.072  -15.072    1.122)   21.345
  10.115   ( -30.050  -30.050   -0.124)   42.498
  12.320   ( -12.001  -12.001   -0.811)   16.991
  12.371   ( -11.929  -11.929   -0.252)   16.873
  13.034   (   4.161    4.161    1.535)    6.081
  15.377   (   5.555    5.555   -1.152)    7.940
  16.205   (   6.237    6.237    0.741)    8.851
  16.798   (   3.954    3.954   -0.804)    5.649
======================= Grid point 409 (83/84) =======================
q-point: ( 0.46  0.38  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 423
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.047   (   0.385    1.441    0.517)    1.578
   2.745   (   1.441    7.466   -0.126)    7.605
   2.787   (   1.952    7.677    1.236)    8.017
   2.844   (   2.785    1.796   -1.970)    3.855
   2.883   (   1.590    2.662    0.040)    3.101
   3.002   (   7.032   11.190    0.563)   13.228
   3.188   (  -0.546    7.703    1.429)    7.853
   3.446   (   0.744    3.022   -1.425)    3.423
   3.887   (   5.077   11.700    0.088)   12.754
   4.245   (  -2.582   -7.001    0.143)    7.463
   5.111   (   3.038    8.306   -0.033)    8.844
   5.149   (   2.583    6.070   -0.005)    6.597
   5.876   (  -9.735  -33.609    0.616)   34.996
   7.715   (  -0.514  -18.367   -0.408)   18.378
   7.831   (  13.353   -4.550   -1.688)   14.208
   8.823   ( -21.029    2.903    1.707)   21.297
   9.753   (  -6.028  -35.720    1.923)   36.276
   9.787   (  -4.770  -14.046    0.272)   14.837
  12.056   (  -8.266   -8.140   -0.566)   11.615
  12.272   (  -2.213  -21.280   -0.588)   21.403
  13.086   (   1.235    4.288    1.551)    4.725
  15.454   (   2.036    6.538   -1.375)    6.985
  16.286   (   1.965    5.791    0.448)    6.132
  16.844   (   1.008    4.924   -0.572)    5.059
======================= Grid point 422 (84/84) =======================
q-point: ( 0.46  0.46  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 7.23e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.067   (   0.522    0.522    0.573)    0.934
   2.845   (   2.146    2.146    0.167)    3.039
   2.853   (   1.342    1.342   -0.085)    1.900
   2.912   (   2.572    2.572    0.069)    3.638
   2.925   (   1.289    1.289    0.219)    1.836
   3.178   (   5.420    5.420    0.809)    7.708
   3.298   (   2.778    2.778    2.118)    4.464
   3.487   (   1.110    1.110   -1.984)    2.530
   4.055   (   5.218    5.218   -0.438)    7.392
   4.144   (  -3.276   -3.276   -0.376)    4.648
   5.228   (   3.243    3.243    0.076)    4.587
   5.237   (   2.374    2.374    0.083)    3.358
   5.388   ( -13.677  -13.677    0.772)   19.357
   7.467   (  -6.376   -6.376   -0.666)    9.042
   7.935   (  11.145   11.145   -2.259)   15.923
   8.842   (  -5.290   -5.290    2.550)    7.904
   9.039   ( -21.273  -21.273    1.494)   30.122
   9.729   (   0.109    0.109    1.290)    1.299
  11.932   (  -5.828   -5.828   -0.645)    8.267
  11.948   (  -7.245   -7.245   -0.564)   10.261
  13.140   (   1.265    1.265    1.644)    2.430
  15.543   (   2.304    2.304   -1.580)    3.622
  16.359   (   1.794    1.794    0.041)    2.537
  16.904   (   1.358    1.358   -0.368)    1.956
=================== 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/20280
   10.0    905.362    905.362    337.696      0.000      0.000      0.000 3/20280
   20.0    287.941    287.941     43.556      0.000     -0.000      0.000 3/20280
   30.0    155.438    155.438     19.124      0.000     -0.000      0.000 3/20280
   40.0    101.778    101.778     12.137      0.000     -0.000      0.000 3/20280
   50.0     74.590     74.590      8.867      0.000     -0.000      0.000 3/20280
   60.0     58.329     58.329      6.946      0.000     -0.000      0.000 3/20280
   70.0     47.496     47.496      5.673      0.000     -0.000      0.000 3/20280
   80.0     39.778     39.778      4.768      0.000     -0.000      0.000 3/20280
   90.0     34.035     34.035      4.094      0.000     -0.000      0.000 3/20280
  100.0     29.631     29.631      3.577      0.000     -0.000      0.000 3/20280
  110.0     26.171     26.171      3.169      0.000     -0.000      0.000 3/20280
  120.0     23.399     23.399      2.841      0.000     -0.000      0.000 3/20280
  130.0     21.140     21.140      2.573      0.000     -0.000      0.000 3/20280
  140.0     19.270     19.270      2.350      0.000     -0.000      0.000 3/20280
  150.0     17.701     17.701      2.162      0.000     -0.000      0.000 3/20280
  160.0     16.368     16.368      2.003      0.000     -0.000      0.000 3/20280
  170.0     15.223     15.223      1.865      0.000     -0.000      0.000 3/20280
  180.0     14.230     14.230      1.746      0.000     -0.000      0.000 3/20280
  190.0     13.360     13.360      1.641      0.000     -0.000      0.000 3/20280
  200.0     12.593     12.593      1.548      0.000     -0.000      0.000 3/20280
  210.0     11.911     11.911      1.466      0.000     -0.000      0.000 3/20280
  220.0     11.301     11.301      1.392      0.000     -0.000      0.000 3/20280
  230.0     10.753     10.753      1.325      0.000     -0.000      0.000 3/20280
  240.0     10.256     10.256      1.265      0.000     -0.000      0.000 3/20280
  250.0      9.804      9.804      1.210      0.000     -0.000      0.000 3/20280
  260.0      9.392      9.392      1.160      0.000     -0.000      0.000 3/20280
  270.0      9.014      9.014      1.114      0.000     -0.000      0.000 3/20280
  280.0      8.666      8.666      1.071      0.000     -0.000      0.000 3/20280
  290.0      8.344      8.344      1.032      0.000     -0.000      0.000 3/20280
  300.0      8.046      8.046      0.995      0.000     -0.000      0.000 3/20280
  310.0      7.769      7.769      0.961      0.000     -0.000      0.000 3/20280
  320.0      7.511      7.511      0.930      0.000     -0.000      0.000 3/20280
  330.0      7.270      7.270      0.900      0.000     -0.000      0.000 3/20280
  340.0      7.044      7.044      0.873      0.000     -0.000      0.000 3/20280
  350.0      6.832      6.832      0.847      0.000     -0.000      0.000 3/20280
  360.0      6.633      6.633      0.822      0.000     -0.000      0.000 3/20280
  370.0      6.445      6.445      0.799      0.000     -0.000      0.000 3/20280
  380.0      6.268      6.268      0.777      0.000     -0.000      0.000 3/20280
  390.0      6.101      6.101      0.757      0.000     -0.000      0.000 3/20280
  400.0      5.942      5.942      0.737      0.000     -0.000      0.000 3/20280
  410.0      5.791      5.791      0.719      0.000     -0.000      0.000 3/20280
  420.0      5.648      5.648      0.701      0.000     -0.000      0.000 3/20280
  430.0      5.513      5.513      0.684      0.000     -0.000      0.000 3/20280
  440.0      5.383      5.383      0.668      0.000     -0.000      0.000 3/20280
  450.0      5.260      5.260      0.653      0.000     -0.000      0.000 3/20280
  460.0      5.142      5.142      0.638      0.000     -0.000      0.000 3/20280
  470.0      5.029      5.029      0.625      0.000     -0.000      0.000 3/20280
  480.0      4.922      4.922      0.611      0.000     -0.000      0.000 3/20280
  490.0      4.818      4.818      0.599      0.000     -0.000      0.000 3/20280
  500.0      4.720      4.720      0.586      0.000     -0.000      0.000 3/20280
  510.0      4.625      4.625      0.575      0.000     -0.000      0.000 3/20280
  520.0      4.534      4.534      0.563      0.000     -0.000      0.000 3/20280
  530.0      4.446      4.446      0.553      0.000     -0.000      0.000 3/20280
  540.0      4.362      4.362      0.542      0.000     -0.000      0.000 3/20280
  550.0      4.281      4.281      0.532      0.000     -0.000      0.000 3/20280
  560.0      4.203      4.203      0.522      0.000     -0.000      0.000 3/20280
  570.0      4.128      4.128      0.513      0.000     -0.000      0.000 3/20280
  580.0      4.055      4.055      0.504      0.000     -0.000      0.000 3/20280
  590.0      3.985      3.985      0.495      0.000     -0.000      0.000 3/20280
  600.0      3.918      3.918      0.487      0.000     -0.000      0.000 3/20280
  610.0      3.853      3.853      0.479      0.000     -0.000      0.000 3/20280
  620.0      3.789      3.789      0.471      0.000     -0.000      0.000 3/20280
  630.0      3.728      3.728      0.464      0.000     -0.000      0.000 3/20280
  640.0      3.669      3.669      0.456      0.000     -0.000      0.000 3/20280
  650.0      3.612      3.612      0.449      0.000     -0.000      0.000 3/20280
  660.0      3.556      3.556      0.442      0.000     -0.000      0.000 3/20280
  670.0      3.502      3.502      0.436      0.000     -0.000      0.000 3/20280
  680.0      3.450      3.450      0.429      0.000     -0.000      0.000 3/20280
  690.0      3.400      3.400      0.423      0.000     -0.000      0.000 3/20280
  700.0      3.350      3.350      0.417      0.000     -0.000      0.000 3/20280
  710.0      3.303      3.303      0.411      0.000     -0.000      0.000 3/20280
  720.0      3.256      3.256      0.405      0.000     -0.000      0.000 3/20280
  730.0      3.211      3.211      0.399      0.000     -0.000      0.000 3/20280
  740.0      3.167      3.167      0.394      0.000     -0.000      0.000 3/20280
  750.0      3.124      3.124      0.389      0.000     -0.000      0.000 3/20280
  760.0      3.083      3.083      0.384      0.000     -0.000      0.000 3/20280
  770.0      3.042      3.042      0.379      0.000     -0.000      0.000 3/20280
  780.0      3.003      3.003      0.374      0.000     -0.000      0.000 3/20280
  790.0      2.965      2.965      0.369      0.000     -0.000      0.000 3/20280
  800.0      2.927      2.927      0.364      0.000     -0.000      0.000 3/20280
  810.0      2.891      2.891      0.360      0.000     -0.000      0.000 3/20280
  820.0      2.855      2.855      0.355      0.000     -0.000      0.000 3/20280
  830.0      2.820      2.820      0.351      0.000     -0.000      0.000 3/20280
  840.0      2.787      2.787      0.347      0.000     -0.000      0.000 3/20280
  850.0      2.754      2.754      0.343      0.000     -0.000      0.000 3/20280
  860.0      2.721      2.721      0.339      0.000     -0.000      0.000 3/20280
  870.0      2.690      2.690      0.335      0.000     -0.000      0.000 3/20280
  880.0      2.659      2.659      0.331      0.000     -0.000      0.000 3/20280
  890.0      2.629      2.629      0.327      0.000     -0.000      0.000 3/20280
  900.0      2.599      2.599      0.324      0.000     -0.000      0.000 3/20280
  910.0      2.571      2.571      0.320      0.000     -0.000      0.000 3/20280
  920.0      2.542      2.542      0.316      0.000     -0.000      0.000 3/20280
  930.0      2.515      2.515      0.313      0.000     -0.000      0.000 3/20280
  940.0      2.488      2.488      0.310      0.000     -0.000      0.000 3/20280
  950.0      2.462      2.462      0.306      0.000     -0.000      0.000 3/20280
  960.0      2.436      2.436      0.303      0.000     -0.000      0.000 3/20280
  970.0      2.410      2.410      0.300      0.000     -0.000      0.000 3/20280
  980.0      2.386      2.386      0.297      0.000     -0.000      0.000 3/20280
  990.0      2.361      2.361      0.294      0.000     -0.000      0.000 3/20280
 1000.0      2.338      2.338      0.291      0.000     -0.000      0.000 3/20280

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

