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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [3 3 3]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: Pm-3m (221)
Number of symmetry operations in supercell: 1296
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.933007870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.933007870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.933007870000000
Atomic positions (fractional):
   *1 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941
    2 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941
    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941
   *4 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999
   *5 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.933007870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.933007870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.933007870000000
Atomic positions (fractional):
   *1 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1
    2 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 2
    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 3
   *4 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 4
   *5 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   11.799023609999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.799023609999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.799023609999999
Atomic positions (fractional):
   *1 Li  0.16666666666667  0.00000000000000  0.00000000000000   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.00000000000000   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.00000000000000   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.00000000000000   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.00000000000000   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.00000000000000   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.00000000000000   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.00000000000000   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.33333333333333   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.33333333333333   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.33333333333333   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.33333333333333   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.33333333333333   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.33333333333333   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.33333333333333   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.33333333333333   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.33333333333333   6.941 > 1
   19 Li  0.16666666666667  0.00000000000000  0.66666666666667   6.941 > 1
   20 Li  0.50000000000000  0.00000000000000  0.66666666666667   6.941 > 1
   21 Li  0.83333333333333  0.00000000000000  0.66666666666667   6.941 > 1
   22 Li  0.16666666666667  0.33333333333333  0.66666666666667   6.941 > 1
   23 Li  0.50000000000000  0.33333333333333  0.66666666666667   6.941 > 1
   24 Li  0.83333333333333  0.33333333333333  0.66666666666667   6.941 > 1
   25 Li  0.16666666666667  0.66666666666667  0.66666666666667   6.941 > 1
   26 Li  0.50000000000000  0.66666666666667  0.66666666666667   6.941 > 1
   27 Li  0.83333333333333  0.66666666666667  0.66666666666667   6.941 > 1
   28 Li  0.00000000000000  0.00000000000000  0.16666666666667   6.941 > 2
   29 Li  0.33333333333333  0.00000000000000  0.16666666666667   6.941 > 2
   30 Li  0.66666666666667  0.00000000000000  0.16666666666667   6.941 > 2
   31 Li  0.00000000000000  0.33333333333333  0.16666666666667   6.941 > 2
   32 Li  0.33333333333333  0.33333333333333  0.16666666666667   6.941 > 2
   33 Li  0.66666666666667  0.33333333333333  0.16666666666667   6.941 > 2
   34 Li  0.00000000000000  0.66666666666667  0.16666666666667   6.941 > 2
   35 Li  0.33333333333333  0.66666666666667  0.16666666666667   6.941 > 2
   36 Li  0.66666666666667  0.66666666666667  0.16666666666667   6.941 > 2
   37 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 2
   38 Li  0.33333333333333  0.00000000000000  0.50000000000000   6.941 > 2
   39 Li  0.66666666666667  0.00000000000000  0.50000000000000   6.941 > 2
   40 Li  0.00000000000000  0.33333333333333  0.50000000000000   6.941 > 2
   41 Li  0.33333333333333  0.33333333333333  0.50000000000000   6.941 > 2
   42 Li  0.66666666666667  0.33333333333333  0.50000000000000   6.941 > 2
   43 Li  0.00000000000000  0.66666666666667  0.50000000000000   6.941 > 2
   44 Li  0.33333333333333  0.66666666666667  0.50000000000000   6.941 > 2
   45 Li  0.66666666666667  0.66666666666667  0.50000000000000   6.941 > 2
   46 Li  0.00000000000000  0.00000000000000  0.83333333333333   6.941 > 2
   47 Li  0.33333333333333  0.00000000000000  0.83333333333333   6.941 > 2
   48 Li  0.66666666666667  0.00000000000000  0.83333333333333   6.941 > 2
   49 Li  0.00000000000000  0.33333333333333  0.83333333333333   6.941 > 2
   50 Li  0.33333333333333  0.33333333333333  0.83333333333333   6.941 > 2
   51 Li  0.66666666666667  0.33333333333333  0.83333333333333   6.941 > 2
   52 Li  0.00000000000000  0.66666666666667  0.83333333333333   6.941 > 2
   53 Li  0.33333333333333  0.66666666666667  0.83333333333333   6.941 > 2
   54 Li  0.66666666666667  0.66666666666667  0.83333333333333   6.941 > 2
   55 Li  0.00000000000000  0.16666666666667  0.00000000000000   6.941 > 3
   56 Li  0.33333333333333  0.16666666666667  0.00000000000000   6.941 > 3
   57 Li  0.66666666666667  0.16666666666667  0.00000000000000   6.941 > 3
   58 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 3
   59 Li  0.33333333333333  0.50000000000000  0.00000000000000   6.941 > 3
   60 Li  0.66666666666667  0.50000000000000  0.00000000000000   6.941 > 3
   61 Li  0.00000000000000  0.83333333333333  0.00000000000000   6.941 > 3
   62 Li  0.33333333333333  0.83333333333333  0.00000000000000   6.941 > 3
   63 Li  0.66666666666667  0.83333333333333  0.00000000000000   6.941 > 3
   64 Li  0.00000000000000  0.16666666666667  0.33333333333333   6.941 > 3
   65 Li  0.33333333333333  0.16666666666667  0.33333333333333   6.941 > 3
   66 Li  0.66666666666667  0.16666666666667  0.33333333333333   6.941 > 3
   67 Li  0.00000000000000  0.50000000000000  0.33333333333333   6.941 > 3
   68 Li  0.33333333333333  0.50000000000000  0.33333333333333   6.941 > 3
   69 Li  0.66666666666667  0.50000000000000  0.33333333333333   6.941 > 3
   70 Li  0.00000000000000  0.83333333333333  0.33333333333333   6.941 > 3
   71 Li  0.33333333333333  0.83333333333333  0.33333333333333   6.941 > 3
   72 Li  0.66666666666667  0.83333333333333  0.33333333333333   6.941 > 3
   73 Li  0.00000000000000  0.16666666666667  0.66666666666667   6.941 > 3
   74 Li  0.33333333333333  0.16666666666667  0.66666666666667   6.941 > 3
   75 Li  0.66666666666667  0.16666666666667  0.66666666666667   6.941 > 3
   76 Li  0.00000000000000  0.50000000000000  0.66666666666667   6.941 > 3
   77 Li  0.33333333333333  0.50000000000000  0.66666666666667   6.941 > 3
   78 Li  0.66666666666667  0.50000000000000  0.66666666666667   6.941 > 3
   79 Li  0.00000000000000  0.83333333333333  0.66666666666667   6.941 > 3
   80 Li  0.33333333333333  0.83333333333333  0.66666666666667   6.941 > 3
   81 Li  0.66666666666667  0.83333333333333  0.66666666666667   6.941 > 3
  *82 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 4
   83 O   0.33333333333333  0.00000000000000  0.00000000000000  15.999 > 4
   84 O   0.66666666666667  0.00000000000000  0.00000000000000  15.999 > 4
   85 O   0.00000000000000  0.33333333333333  0.00000000000000  15.999 > 4
   86 O   0.33333333333333  0.33333333333333  0.00000000000000  15.999 > 4
   87 O   0.66666666666667  0.33333333333333  0.00000000000000  15.999 > 4
   88 O   0.00000000000000  0.66666666666667  0.00000000000000  15.999 > 4
   89 O   0.33333333333333  0.66666666666667  0.00000000000000  15.999 > 4
   90 O   0.66666666666667  0.66666666666667  0.00000000000000  15.999 > 4
   91 O   0.00000000000000  0.00000000000000  0.33333333333333  15.999 > 4
   92 O   0.33333333333333  0.00000000000000  0.33333333333333  15.999 > 4
   93 O   0.66666666666667  0.00000000000000  0.33333333333333  15.999 > 4
   94 O   0.00000000000000  0.33333333333333  0.33333333333333  15.999 > 4
   95 O   0.33333333333333  0.33333333333333  0.33333333333333  15.999 > 4
   96 O   0.66666666666667  0.33333333333333  0.33333333333333  15.999 > 4
   97 O   0.00000000000000  0.66666666666667  0.33333333333333  15.999 > 4
   98 O   0.33333333333333  0.66666666666667  0.33333333333333  15.999 > 4
   99 O   0.66666666666667  0.66666666666667  0.33333333333333  15.999 > 4
  100 O   0.00000000000000  0.00000000000000  0.66666666666667  15.999 > 4
  101 O   0.33333333333333  0.00000000000000  0.66666666666667  15.999 > 4
  102 O   0.66666666666667  0.00000000000000  0.66666666666667  15.999 > 4
  103 O   0.00000000000000  0.33333333333333  0.66666666666667  15.999 > 4
  104 O   0.33333333333333  0.33333333333333  0.66666666666667  15.999 > 4
  105 O   0.66666666666667  0.33333333333333  0.66666666666667  15.999 > 4
  106 O   0.00000000000000  0.66666666666667  0.66666666666667  15.999 > 4
  107 O   0.33333333333333  0.66666666666667  0.66666666666667  15.999 > 4
  108 O   0.66666666666667  0.66666666666667  0.66666666666667  15.999 > 4
 *109 Br  0.16666666666667  0.16666666666667  0.16666666666667  79.904 > 5
  110 Br  0.50000000000000  0.16666666666667  0.16666666666667  79.904 > 5
  111 Br  0.83333333333333  0.16666666666667  0.16666666666667  79.904 > 5
  112 Br  0.16666666666667  0.50000000000000  0.16666666666667  79.904 > 5
  113 Br  0.50000000000000  0.50000000000000  0.16666666666667  79.904 > 5
  114 Br  0.83333333333333  0.50000000000000  0.16666666666667  79.904 > 5
  115 Br  0.16666666666667  0.83333333333333  0.16666666666667  79.904 > 5
  116 Br  0.50000000000000  0.83333333333333  0.16666666666667  79.904 > 5
  117 Br  0.83333333333333  0.83333333333333  0.16666666666667  79.904 > 5
  118 Br  0.16666666666667  0.16666666666667  0.50000000000000  79.904 > 5
  119 Br  0.50000000000000  0.16666666666667  0.50000000000000  79.904 > 5
  120 Br  0.83333333333333  0.16666666666667  0.50000000000000  79.904 > 5
  121 Br  0.16666666666667  0.50000000000000  0.50000000000000  79.904 > 5
  122 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 5
  123 Br  0.83333333333333  0.50000000000000  0.50000000000000  79.904 > 5
  124 Br  0.16666666666667  0.83333333333333  0.50000000000000  79.904 > 5
  125 Br  0.50000000000000  0.83333333333333  0.50000000000000  79.904 > 5
  126 Br  0.83333333333333  0.83333333333333  0.50000000000000  79.904 > 5
  127 Br  0.16666666666667  0.16666666666667  0.83333333333333  79.904 > 5
  128 Br  0.50000000000000  0.16666666666667  0.83333333333333  79.904 > 5
  129 Br  0.83333333333333  0.16666666666667  0.83333333333333  79.904 > 5
  130 Br  0.16666666666667  0.50000000000000  0.83333333333333  79.904 > 5
  131 Br  0.50000000000000  0.50000000000000  0.83333333333333  79.904 > 5
  132 Br  0.83333333333333  0.50000000000000  0.83333333333333  79.904 > 5
  133 Br  0.16666666666667  0.83333333333333  0.83333333333333  79.904 > 5
  134 Br  0.50000000000000  0.83333333333333  0.83333333333333  79.904 > 5
  135 Br  0.83333333333333  0.83333333333333  0.83333333333333  79.904 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.7131400    0.0000000    0.0000000
            0.0000000    3.7131400    0.0000000
            0.0000000    0.0000000    3.7131400
-------------------------- Born effective charges --------------------------
    1 Li    0.9769099    0.0000000    0.0000000
            0.0000000    1.2086006    0.0000000
            0.0000000    0.0000000    1.2086006
    2 Li    1.2086006    0.0000000    0.0000000
            0.0000000    1.2086006    0.0000000
            0.0000000    0.0000000    0.9769099
    3 Li    1.2086006    0.0000000    0.0000000
            0.0000000    0.9769099    0.0000000
            0.0000000    0.0000000    1.2086006
    4 O    -2.1770128    0.0000000    0.0000000
            0.0000000   -2.1770128    0.0000000
            0.0000000    0.0000000   -2.1770128
    5 Br   -1.2170983    0.0000000    0.0000000
            0.0000000   -1.2170983    0.0000000
            0.0000000    0.0000000   -1.2170983
----------------------------------------------------------------------------
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: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), data: True
-----
Solver_atoms: 1 -- 135 / 135
Time (Solver_compr_matrix_reshape): 0.003
Solver_block: 80 / 80
 - Time: 0.092
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.096
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) 
Permutation basis: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), 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 07:25:54]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 07:25:54]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [3 3 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: Pm-3m (221)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.933007870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.933007870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.933007870000000
Atomic positions (fractional):
    1 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941
    2 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941
    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941
    4 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999
    5 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.799023609999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.799023609999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.799023609999999
Atomic positions (fractional):
    1 Li  0.16666666666667  0.00000000000000  0.00000000000000   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.00000000000000   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.00000000000000   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.00000000000000   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.00000000000000   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.00000000000000   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.00000000000000   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.00000000000000   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.33333333333333   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.33333333333333   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.33333333333333   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.33333333333333   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.33333333333333   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.33333333333333   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.33333333333333   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.33333333333333   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.33333333333333   6.941 > 1
   19 Li  0.16666666666667  0.00000000000000  0.66666666666667   6.941 > 1
   20 Li  0.50000000000000  0.00000000000000  0.66666666666667   6.941 > 1
   21 Li  0.83333333333333  0.00000000000000  0.66666666666667   6.941 > 1
   22 Li  0.16666666666667  0.33333333333333  0.66666666666667   6.941 > 1
   23 Li  0.50000000000000  0.33333333333333  0.66666666666667   6.941 > 1
   24 Li  0.83333333333333  0.33333333333333  0.66666666666667   6.941 > 1
   25 Li  0.16666666666667  0.66666666666667  0.66666666666667   6.941 > 1
   26 Li  0.50000000000000  0.66666666666667  0.66666666666667   6.941 > 1
   27 Li  0.83333333333333  0.66666666666667  0.66666666666667   6.941 > 1
   28 Li  0.00000000000000  0.00000000000000  0.16666666666667   6.941 > 28
   29 Li  0.33333333333333  0.00000000000000  0.16666666666667   6.941 > 28
   30 Li  0.66666666666667  0.00000000000000  0.16666666666667   6.941 > 28
   31 Li  0.00000000000000  0.33333333333333  0.16666666666667   6.941 > 28
   32 Li  0.33333333333333  0.33333333333333  0.16666666666667   6.941 > 28
   33 Li  0.66666666666667  0.33333333333333  0.16666666666667   6.941 > 28
   34 Li  0.00000000000000  0.66666666666667  0.16666666666667   6.941 > 28
   35 Li  0.33333333333333  0.66666666666667  0.16666666666667   6.941 > 28
   36 Li  0.66666666666667  0.66666666666667  0.16666666666667   6.941 > 28
   37 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 28
   38 Li  0.33333333333333  0.00000000000000  0.50000000000000   6.941 > 28
   39 Li  0.66666666666667  0.00000000000000  0.50000000000000   6.941 > 28
   40 Li  0.00000000000000  0.33333333333333  0.50000000000000   6.941 > 28
   41 Li  0.33333333333333  0.33333333333333  0.50000000000000   6.941 > 28
   42 Li  0.66666666666667  0.33333333333333  0.50000000000000   6.941 > 28
   43 Li  0.00000000000000  0.66666666666667  0.50000000000000   6.941 > 28
   44 Li  0.33333333333333  0.66666666666667  0.50000000000000   6.941 > 28
   45 Li  0.66666666666667  0.66666666666667  0.50000000000000   6.941 > 28
   46 Li  0.00000000000000  0.00000000000000  0.83333333333333   6.941 > 28
   47 Li  0.33333333333333  0.00000000000000  0.83333333333333   6.941 > 28
   48 Li  0.66666666666667  0.00000000000000  0.83333333333333   6.941 > 28
   49 Li  0.00000000000000  0.33333333333333  0.83333333333333   6.941 > 28
   50 Li  0.33333333333333  0.33333333333333  0.83333333333333   6.941 > 28
   51 Li  0.66666666666667  0.33333333333333  0.83333333333333   6.941 > 28
   52 Li  0.00000000000000  0.66666666666667  0.83333333333333   6.941 > 28
   53 Li  0.33333333333333  0.66666666666667  0.83333333333333   6.941 > 28
   54 Li  0.66666666666667  0.66666666666667  0.83333333333333   6.941 > 28
   55 Li  0.00000000000000  0.16666666666667  0.00000000000000   6.941 > 55
   56 Li  0.33333333333333  0.16666666666667  0.00000000000000   6.941 > 55
   57 Li  0.66666666666667  0.16666666666667  0.00000000000000   6.941 > 55
   58 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 55
   59 Li  0.33333333333333  0.50000000000000  0.00000000000000   6.941 > 55
   60 Li  0.66666666666667  0.50000000000000  0.00000000000000   6.941 > 55
   61 Li  0.00000000000000  0.83333333333333  0.00000000000000   6.941 > 55
   62 Li  0.33333333333333  0.83333333333333  0.00000000000000   6.941 > 55
   63 Li  0.66666666666667  0.83333333333333  0.00000000000000   6.941 > 55
   64 Li  0.00000000000000  0.16666666666667  0.33333333333333   6.941 > 55
   65 Li  0.33333333333333  0.16666666666667  0.33333333333333   6.941 > 55
   66 Li  0.66666666666667  0.16666666666667  0.33333333333333   6.941 > 55
   67 Li  0.00000000000000  0.50000000000000  0.33333333333333   6.941 > 55
   68 Li  0.33333333333333  0.50000000000000  0.33333333333333   6.941 > 55
   69 Li  0.66666666666667  0.50000000000000  0.33333333333333   6.941 > 55
   70 Li  0.00000000000000  0.83333333333333  0.33333333333333   6.941 > 55
   71 Li  0.33333333333333  0.83333333333333  0.33333333333333   6.941 > 55
   72 Li  0.66666666666667  0.83333333333333  0.33333333333333   6.941 > 55
   73 Li  0.00000000000000  0.16666666666667  0.66666666666667   6.941 > 55
   74 Li  0.33333333333333  0.16666666666667  0.66666666666667   6.941 > 55
   75 Li  0.66666666666667  0.16666666666667  0.66666666666667   6.941 > 55
   76 Li  0.00000000000000  0.50000000000000  0.66666666666667   6.941 > 55
   77 Li  0.33333333333333  0.50000000000000  0.66666666666667   6.941 > 55
   78 Li  0.66666666666667  0.50000000000000  0.66666666666667   6.941 > 55
   79 Li  0.00000000000000  0.83333333333333  0.66666666666667   6.941 > 55
   80 Li  0.33333333333333  0.83333333333333  0.66666666666667   6.941 > 55
   81 Li  0.66666666666667  0.83333333333333  0.66666666666667   6.941 > 55
   82 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 82
   83 O   0.33333333333333  0.00000000000000  0.00000000000000  15.999 > 82
   84 O   0.66666666666667  0.00000000000000  0.00000000000000  15.999 > 82
   85 O   0.00000000000000  0.33333333333333  0.00000000000000  15.999 > 82
   86 O   0.33333333333333  0.33333333333333  0.00000000000000  15.999 > 82
   87 O   0.66666666666667  0.33333333333333  0.00000000000000  15.999 > 82
   88 O   0.00000000000000  0.66666666666667  0.00000000000000  15.999 > 82
   89 O   0.33333333333333  0.66666666666667  0.00000000000000  15.999 > 82
   90 O   0.66666666666667  0.66666666666667  0.00000000000000  15.999 > 82
   91 O   0.00000000000000  0.00000000000000  0.33333333333333  15.999 > 82
   92 O   0.33333333333333  0.00000000000000  0.33333333333333  15.999 > 82
   93 O   0.66666666666667  0.00000000000000  0.33333333333333  15.999 > 82
   94 O   0.00000000000000  0.33333333333333  0.33333333333333  15.999 > 82
   95 O   0.33333333333333  0.33333333333333  0.33333333333333  15.999 > 82
   96 O   0.66666666666667  0.33333333333333  0.33333333333333  15.999 > 82
   97 O   0.00000000000000  0.66666666666667  0.33333333333333  15.999 > 82
   98 O   0.33333333333333  0.66666666666667  0.33333333333333  15.999 > 82
   99 O   0.66666666666667  0.66666666666667  0.33333333333333  15.999 > 82
  100 O   0.00000000000000  0.00000000000000  0.66666666666667  15.999 > 82
  101 O   0.33333333333333  0.00000000000000  0.66666666666667  15.999 > 82
  102 O   0.66666666666667  0.00000000000000  0.66666666666667  15.999 > 82
  103 O   0.00000000000000  0.33333333333333  0.66666666666667  15.999 > 82
  104 O   0.33333333333333  0.33333333333333  0.66666666666667  15.999 > 82
  105 O   0.66666666666667  0.33333333333333  0.66666666666667  15.999 > 82
  106 O   0.00000000000000  0.66666666666667  0.66666666666667  15.999 > 82
  107 O   0.33333333333333  0.66666666666667  0.66666666666667  15.999 > 82
  108 O   0.66666666666667  0.66666666666667  0.66666666666667  15.999 > 82
  109 Br  0.16666666666667  0.16666666666667  0.16666666666667  79.904 > 109
  110 Br  0.50000000000000  0.16666666666667  0.16666666666667  79.904 > 109
  111 Br  0.83333333333333  0.16666666666667  0.16666666666667  79.904 > 109
  112 Br  0.16666666666667  0.50000000000000  0.16666666666667  79.904 > 109
  113 Br  0.50000000000000  0.50000000000000  0.16666666666667  79.904 > 109
  114 Br  0.83333333333333  0.50000000000000  0.16666666666667  79.904 > 109
  115 Br  0.16666666666667  0.83333333333333  0.16666666666667  79.904 > 109
  116 Br  0.50000000000000  0.83333333333333  0.16666666666667  79.904 > 109
  117 Br  0.83333333333333  0.83333333333333  0.16666666666667  79.904 > 109
  118 Br  0.16666666666667  0.16666666666667  0.50000000000000  79.904 > 109
  119 Br  0.50000000000000  0.16666666666667  0.50000000000000  79.904 > 109
  120 Br  0.83333333333333  0.16666666666667  0.50000000000000  79.904 > 109
  121 Br  0.16666666666667  0.50000000000000  0.50000000000000  79.904 > 109
  122 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 109
  123 Br  0.83333333333333  0.50000000000000  0.50000000000000  79.904 > 109
  124 Br  0.16666666666667  0.83333333333333  0.50000000000000  79.904 > 109
  125 Br  0.50000000000000  0.83333333333333  0.50000000000000  79.904 > 109
  126 Br  0.83333333333333  0.83333333333333  0.50000000000000  79.904 > 109
  127 Br  0.16666666666667  0.16666666666667  0.83333333333333  79.904 > 109
  128 Br  0.50000000000000  0.16666666666667  0.83333333333333  79.904 > 109
  129 Br  0.83333333333333  0.16666666666667  0.83333333333333  79.904 > 109
  130 Br  0.16666666666667  0.50000000000000  0.83333333333333  79.904 > 109
  131 Br  0.50000000000000  0.50000000000000  0.83333333333333  79.904 > 109
  132 Br  0.83333333333333  0.50000000000000  0.83333333333333  79.904 > 109
  133 Br  0.16666666666667  0.83333333333333  0.83333333333333  79.904 > 109
  134 Br  0.50000000000000  0.83333333333333  0.83333333333333  79.904 > 109
  135 Br  0.83333333333333  0.83333333333333  0.83333333333333  79.904 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.7131400    0.0000000    0.0000000
            0.0000000    3.7131400    0.0000000
            0.0000000    0.0000000    3.7131400
-------------------------- Born effective charges --------------------------
    1 Li    0.9769099    0.0000000    0.0000000
            0.0000000    1.2086006    0.0000000
            0.0000000    0.0000000    1.2086006
    2 Li    1.2086006    0.0000000    0.0000000
            0.0000000    1.2086006    0.0000000
            0.0000000    0.0000000    0.9769099
    3 Li    1.2086006    0.0000000    0.0000000
            0.0000000    0.9769099    0.0000000
            0.0000000    0.0000000    1.2086006
    4 O    -2.1770128    0.0000000    0.0000000
            0.0000000   -2.1770128    0.0000000
            0.0000000    0.0000000   -2.1770128
    5 Br   -1.2170983    0.0000000    0.0000000
            0.0000000   -1.2170983    0.0000000
            0.0000000    0.0000000   -1.2170983
----------------------------------------------------------------------------
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.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  0.0000]
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 82, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 109, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 07:25:58]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 07:25:58]-------------------------
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 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: Pm-3m (221)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.933007870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.933007870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.933007870000000
Atomic positions (fractional):
    1 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941
    2 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941
    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941
    4 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999
    5 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.799023609999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.799023609999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.799023609999999
Atomic positions (fractional):
    1 Li  0.16666666666667  0.00000000000000  0.00000000000000   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1
    3 Li  0.83333333333333  0.00000000000000  0.00000000000000   6.941 > 1
    4 Li  0.16666666666667  0.33333333333333  0.00000000000000   6.941 > 1
    5 Li  0.50000000000000  0.33333333333333  0.00000000000000   6.941 > 1
    6 Li  0.83333333333333  0.33333333333333  0.00000000000000   6.941 > 1
    7 Li  0.16666666666667  0.66666666666667  0.00000000000000   6.941 > 1
    8 Li  0.50000000000000  0.66666666666667  0.00000000000000   6.941 > 1
    9 Li  0.83333333333333  0.66666666666667  0.00000000000000   6.941 > 1
   10 Li  0.16666666666667  0.00000000000000  0.33333333333333   6.941 > 1
   11 Li  0.50000000000000  0.00000000000000  0.33333333333333   6.941 > 1
   12 Li  0.83333333333333  0.00000000000000  0.33333333333333   6.941 > 1
   13 Li  0.16666666666667  0.33333333333333  0.33333333333333   6.941 > 1
   14 Li  0.50000000000000  0.33333333333333  0.33333333333333   6.941 > 1
   15 Li  0.83333333333333  0.33333333333333  0.33333333333333   6.941 > 1
   16 Li  0.16666666666667  0.66666666666667  0.33333333333333   6.941 > 1
   17 Li  0.50000000000000  0.66666666666667  0.33333333333333   6.941 > 1
   18 Li  0.83333333333333  0.66666666666667  0.33333333333333   6.941 > 1
   19 Li  0.16666666666667  0.00000000000000  0.66666666666667   6.941 > 1
   20 Li  0.50000000000000  0.00000000000000  0.66666666666667   6.941 > 1
   21 Li  0.83333333333333  0.00000000000000  0.66666666666667   6.941 > 1
   22 Li  0.16666666666667  0.33333333333333  0.66666666666667   6.941 > 1
   23 Li  0.50000000000000  0.33333333333333  0.66666666666667   6.941 > 1
   24 Li  0.83333333333333  0.33333333333333  0.66666666666667   6.941 > 1
   25 Li  0.16666666666667  0.66666666666667  0.66666666666667   6.941 > 1
   26 Li  0.50000000000000  0.66666666666667  0.66666666666667   6.941 > 1
   27 Li  0.83333333333333  0.66666666666667  0.66666666666667   6.941 > 1
   28 Li  0.00000000000000  0.00000000000000  0.16666666666667   6.941 > 28
   29 Li  0.33333333333333  0.00000000000000  0.16666666666667   6.941 > 28
   30 Li  0.66666666666667  0.00000000000000  0.16666666666667   6.941 > 28
   31 Li  0.00000000000000  0.33333333333333  0.16666666666667   6.941 > 28
   32 Li  0.33333333333333  0.33333333333333  0.16666666666667   6.941 > 28
   33 Li  0.66666666666667  0.33333333333333  0.16666666666667   6.941 > 28
   34 Li  0.00000000000000  0.66666666666667  0.16666666666667   6.941 > 28
   35 Li  0.33333333333333  0.66666666666667  0.16666666666667   6.941 > 28
   36 Li  0.66666666666667  0.66666666666667  0.16666666666667   6.941 > 28
   37 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 28
   38 Li  0.33333333333333  0.00000000000000  0.50000000000000   6.941 > 28
   39 Li  0.66666666666667  0.00000000000000  0.50000000000000   6.941 > 28
   40 Li  0.00000000000000  0.33333333333333  0.50000000000000   6.941 > 28
   41 Li  0.33333333333333  0.33333333333333  0.50000000000000   6.941 > 28
   42 Li  0.66666666666667  0.33333333333333  0.50000000000000   6.941 > 28
   43 Li  0.00000000000000  0.66666666666667  0.50000000000000   6.941 > 28
   44 Li  0.33333333333333  0.66666666666667  0.50000000000000   6.941 > 28
   45 Li  0.66666666666667  0.66666666666667  0.50000000000000   6.941 > 28
   46 Li  0.00000000000000  0.00000000000000  0.83333333333333   6.941 > 28
   47 Li  0.33333333333333  0.00000000000000  0.83333333333333   6.941 > 28
   48 Li  0.66666666666667  0.00000000000000  0.83333333333333   6.941 > 28
   49 Li  0.00000000000000  0.33333333333333  0.83333333333333   6.941 > 28
   50 Li  0.33333333333333  0.33333333333333  0.83333333333333   6.941 > 28
   51 Li  0.66666666666667  0.33333333333333  0.83333333333333   6.941 > 28
   52 Li  0.00000000000000  0.66666666666667  0.83333333333333   6.941 > 28
   53 Li  0.33333333333333  0.66666666666667  0.83333333333333   6.941 > 28
   54 Li  0.66666666666667  0.66666666666667  0.83333333333333   6.941 > 28
   55 Li  0.00000000000000  0.16666666666667  0.00000000000000   6.941 > 55
   56 Li  0.33333333333333  0.16666666666667  0.00000000000000   6.941 > 55
   57 Li  0.66666666666667  0.16666666666667  0.00000000000000   6.941 > 55
   58 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 55
   59 Li  0.33333333333333  0.50000000000000  0.00000000000000   6.941 > 55
   60 Li  0.66666666666667  0.50000000000000  0.00000000000000   6.941 > 55
   61 Li  0.00000000000000  0.83333333333333  0.00000000000000   6.941 > 55
   62 Li  0.33333333333333  0.83333333333333  0.00000000000000   6.941 > 55
   63 Li  0.66666666666667  0.83333333333333  0.00000000000000   6.941 > 55
   64 Li  0.00000000000000  0.16666666666667  0.33333333333333   6.941 > 55
   65 Li  0.33333333333333  0.16666666666667  0.33333333333333   6.941 > 55
   66 Li  0.66666666666667  0.16666666666667  0.33333333333333   6.941 > 55
   67 Li  0.00000000000000  0.50000000000000  0.33333333333333   6.941 > 55
   68 Li  0.33333333333333  0.50000000000000  0.33333333333333   6.941 > 55
   69 Li  0.66666666666667  0.50000000000000  0.33333333333333   6.941 > 55
   70 Li  0.00000000000000  0.83333333333333  0.33333333333333   6.941 > 55
   71 Li  0.33333333333333  0.83333333333333  0.33333333333333   6.941 > 55
   72 Li  0.66666666666667  0.83333333333333  0.33333333333333   6.941 > 55
   73 Li  0.00000000000000  0.16666666666667  0.66666666666667   6.941 > 55
   74 Li  0.33333333333333  0.16666666666667  0.66666666666667   6.941 > 55
   75 Li  0.66666666666667  0.16666666666667  0.66666666666667   6.941 > 55
   76 Li  0.00000000000000  0.50000000000000  0.66666666666667   6.941 > 55
   77 Li  0.33333333333333  0.50000000000000  0.66666666666667   6.941 > 55
   78 Li  0.66666666666667  0.50000000000000  0.66666666666667   6.941 > 55
   79 Li  0.00000000000000  0.83333333333333  0.66666666666667   6.941 > 55
   80 Li  0.33333333333333  0.83333333333333  0.66666666666667   6.941 > 55
   81 Li  0.66666666666667  0.83333333333333  0.66666666666667   6.941 > 55
   82 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 82
   83 O   0.33333333333333  0.00000000000000  0.00000000000000  15.999 > 82
   84 O   0.66666666666667  0.00000000000000  0.00000000000000  15.999 > 82
   85 O   0.00000000000000  0.33333333333333  0.00000000000000  15.999 > 82
   86 O   0.33333333333333  0.33333333333333  0.00000000000000  15.999 > 82
   87 O   0.66666666666667  0.33333333333333  0.00000000000000  15.999 > 82
   88 O   0.00000000000000  0.66666666666667  0.00000000000000  15.999 > 82
   89 O   0.33333333333333  0.66666666666667  0.00000000000000  15.999 > 82
   90 O   0.66666666666667  0.66666666666667  0.00000000000000  15.999 > 82
   91 O   0.00000000000000  0.00000000000000  0.33333333333333  15.999 > 82
   92 O   0.33333333333333  0.00000000000000  0.33333333333333  15.999 > 82
   93 O   0.66666666666667  0.00000000000000  0.33333333333333  15.999 > 82
   94 O   0.00000000000000  0.33333333333333  0.33333333333333  15.999 > 82
   95 O   0.33333333333333  0.33333333333333  0.33333333333333  15.999 > 82
   96 O   0.66666666666667  0.33333333333333  0.33333333333333  15.999 > 82
   97 O   0.00000000000000  0.66666666666667  0.33333333333333  15.999 > 82
   98 O   0.33333333333333  0.66666666666667  0.33333333333333  15.999 > 82
   99 O   0.66666666666667  0.66666666666667  0.33333333333333  15.999 > 82
  100 O   0.00000000000000  0.00000000000000  0.66666666666667  15.999 > 82
  101 O   0.33333333333333  0.00000000000000  0.66666666666667  15.999 > 82
  102 O   0.66666666666667  0.00000000000000  0.66666666666667  15.999 > 82
  103 O   0.00000000000000  0.33333333333333  0.66666666666667  15.999 > 82
  104 O   0.33333333333333  0.33333333333333  0.66666666666667  15.999 > 82
  105 O   0.66666666666667  0.33333333333333  0.66666666666667  15.999 > 82
  106 O   0.00000000000000  0.66666666666667  0.66666666666667  15.999 > 82
  107 O   0.33333333333333  0.66666666666667  0.66666666666667  15.999 > 82
  108 O   0.66666666666667  0.66666666666667  0.66666666666667  15.999 > 82
  109 Br  0.16666666666667  0.16666666666667  0.16666666666667  79.904 > 109
  110 Br  0.50000000000000  0.16666666666667  0.16666666666667  79.904 > 109
  111 Br  0.83333333333333  0.16666666666667  0.16666666666667  79.904 > 109
  112 Br  0.16666666666667  0.50000000000000  0.16666666666667  79.904 > 109
  113 Br  0.50000000000000  0.50000000000000  0.16666666666667  79.904 > 109
  114 Br  0.83333333333333  0.50000000000000  0.16666666666667  79.904 > 109
  115 Br  0.16666666666667  0.83333333333333  0.16666666666667  79.904 > 109
  116 Br  0.50000000000000  0.83333333333333  0.16666666666667  79.904 > 109
  117 Br  0.83333333333333  0.83333333333333  0.16666666666667  79.904 > 109
  118 Br  0.16666666666667  0.16666666666667  0.50000000000000  79.904 > 109
  119 Br  0.50000000000000  0.16666666666667  0.50000000000000  79.904 > 109
  120 Br  0.83333333333333  0.16666666666667  0.50000000000000  79.904 > 109
  121 Br  0.16666666666667  0.50000000000000  0.50000000000000  79.904 > 109
  122 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 109
  123 Br  0.83333333333333  0.50000000000000  0.50000000000000  79.904 > 109
  124 Br  0.16666666666667  0.83333333333333  0.50000000000000  79.904 > 109
  125 Br  0.50000000000000  0.83333333333333  0.50000000000000  79.904 > 109
  126 Br  0.83333333333333  0.83333333333333  0.50000000000000  79.904 > 109
  127 Br  0.16666666666667  0.16666666666667  0.83333333333333  79.904 > 109
  128 Br  0.50000000000000  0.16666666666667  0.83333333333333  79.904 > 109
  129 Br  0.83333333333333  0.16666666666667  0.83333333333333  79.904 > 109
  130 Br  0.16666666666667  0.50000000000000  0.83333333333333  79.904 > 109
  131 Br  0.50000000000000  0.50000000000000  0.83333333333333  79.904 > 109
  132 Br  0.83333333333333  0.50000000000000  0.83333333333333  79.904 > 109
  133 Br  0.16666666666667  0.83333333333333  0.83333333333333  79.904 > 109
  134 Br  0.50000000000000  0.83333333333333  0.83333333333333  79.904 > 109
  135 Br  0.83333333333333  0.83333333333333  0.83333333333333  79.904 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.7131400    0.0000000    0.0000000
            0.0000000    3.7131400    0.0000000
            0.0000000    0.0000000    3.7131400
-------------------------- Born effective charges --------------------------
    1 Li    0.9769099    0.0000000    0.0000000
            0.0000000    1.2086006    0.0000000
            0.0000000    0.0000000    1.2086006
    2 Li    1.2086006    0.0000000    0.0000000
            0.0000000    1.2086006    0.0000000
            0.0000000    0.0000000    0.9769099
    3 Li    1.2086006    0.0000000    0.0000000
            0.0000000    0.9769099    0.0000000
            0.0000000    0.0000000    1.2086006
    4 O    -2.1770128    0.0000000    0.0000000
            0.0000000   -2.1770128    0.0000000
            0.0000000    0.0000000   -2.1770128
    5 Br   -1.2170983    0.0000000    0.0000000
            0.0000000   -1.2170983    0.0000000
            0.0000000    0.0000000   -1.2170983
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)
Max drift of fc2: -0.00000000 (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 13 ]
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: 1.06, Number of G-points: 305, Lambda: 0.21
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: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
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
   4.951   (   0.000    0.000    0.000)    0.000
   4.951   (   0.000    0.000    0.000)    0.000
   4.951   (   0.000    0.000    0.000)    0.000
   7.780   (   0.000    0.000    0.000)    0.000
   7.780   (   0.000    0.000    0.000)    0.000
   7.780   (   0.000    0.000    0.000)    0.000
  10.821   (   0.000    0.000    0.000)    0.000
  10.821   (   0.000    0.000    0.000)    0.000
  10.821   (   0.000    0.000    0.000)    0.000
  14.636   (   0.000    0.000    0.000)    0.000
  14.636   (   0.000    0.000    0.000)    0.000
  14.636   (   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: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.746   (  37.143    0.000    0.000)   37.143
   0.746   (  37.143    0.000    0.000)   37.143
   1.042   (  51.139    0.000    0.000)   51.139
   4.976   (   2.526    0.000    0.000)    2.526
   4.976   (   2.526    0.000    0.000)    2.526
   5.955   (  12.839    0.000    0.000)   12.839
   7.811   (   3.079    0.000    0.000)    3.079
   7.811   (   3.079    0.000    0.000)    3.079
   7.828   (   4.773    0.000    0.000)    4.773
  10.816   (  -0.516    0.000    0.000)    0.516
  10.816   (  -0.516    0.000    0.000)    0.516
  13.449   (   2.466    0.000    0.000)    2.466
  14.629   (  -0.676    0.000    0.000)    0.676
  14.629   (  -0.676    0.000    0.000)    0.676
  17.630   (  -5.116    0.000    0.000)    5.116
======================= Grid point 2 (3/84) =======================
q-point: ( 0.15  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.435   (  32.898    0.000    0.000)   32.898
   1.435   (  32.898    0.000    0.000)   32.898
   1.965   (  42.507    0.000    0.000)   42.507
   5.047   (   4.664    0.000    0.000)    4.664
   5.047   (   4.664    0.000    0.000)    4.664
   6.320   (  24.069    0.000    0.000)   24.069
   7.896   (   5.498    0.000    0.000)    5.498
   7.896   (   5.498    0.000    0.000)    5.498
   7.958   (   8.303    0.000    0.000)    8.303
  10.802   (  -0.938    0.000    0.000)    0.938
  10.802   (  -0.938    0.000    0.000)    0.938
  13.521   (   4.915    0.000    0.000)    4.915
  14.611   (  -1.191    0.000    0.000)    1.191
  14.611   (  -1.191    0.000    0.000)    1.191
  17.482   (  -9.922    0.000    0.000)    9.922
======================= Grid point 3 (4/84) =======================
q-point: ( 0.23  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.019   (  26.537    0.000    0.000)   26.537
   2.019   (  26.537    0.000    0.000)   26.537
   2.688   (  31.381    0.000    0.000)   31.381
   5.153   (   5.979    0.000    0.000)    5.979
   5.153   (   5.979    0.000    0.000)    5.979
   6.866   (  30.703    0.000    0.000)   30.703
   8.018   (   6.706    0.000    0.000)    6.706
   8.018   (   6.706    0.000    0.000)    6.706
   8.139   (   9.773    0.000    0.000)    9.773
  10.781   (  -1.180    0.000    0.000)    1.180
  10.781   (  -1.180    0.000    0.000)    1.180
  13.640   (   7.238    0.000    0.000)    7.238
  14.584   (  -1.425    0.000    0.000)    1.425
  14.584   (  -1.425    0.000    0.000)    1.425
  17.247   ( -13.852    0.000    0.000)   13.852
======================= Grid point 4 (5/84) =======================
q-point: ( 0.31  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.465   (  19.013    0.000    0.000)   19.013
   2.465   (  19.013    0.000    0.000)   19.013
   3.196   (  20.817    0.000    0.000)   20.817
   5.273   (   6.012    0.000    0.000)    6.012
   5.273   (   6.012    0.000    0.000)    6.012
   7.472   (  29.929    0.000    0.000)   29.929
   8.148   (   6.388    0.000    0.000)    6.388
   8.148   (   6.388    0.000    0.000)    6.388
   8.326   (   8.979    0.000    0.000)    8.979
  10.757   (  -1.168    0.000    0.000)    1.168
  10.757   (  -1.168    0.000    0.000)    1.168
  13.800   (   8.948    0.000    0.000)    8.948
  14.557   (  -1.328    0.000    0.000)    1.328
  14.557   (  -1.328    0.000    0.000)    1.328
  16.952   ( -15.800    0.000    0.000)   15.800
======================= Grid point 5 (6/84) =======================
q-point: ( 0.38  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.761   (  11.247    0.000    0.000)   11.247
   2.761   (  11.247    0.000    0.000)   11.247
   3.513   (  11.791    0.000    0.000)   11.791
   5.378   (   4.495    0.000    0.000)    4.495
   5.378   (   4.495    0.000    0.000)    4.495
   7.987   (  21.557    0.000    0.000)   21.557
   8.258   (   4.572    0.000    0.000)    4.572
   8.258   (   4.572    0.000    0.000)    4.572
   8.477   (   6.244    0.000    0.000)    6.244
  10.737   (  -0.867    0.000    0.000)    0.867
  10.737   (  -0.867    0.000    0.000)    0.867
  13.977   (   8.580    0.000    0.000)    8.580
  14.534   (  -0.932    0.000    0.000)    0.932
  14.534   (  -0.932    0.000    0.000)    0.932
  16.654   ( -13.823    0.000    0.000)   13.823
======================= Grid point 6 (7/84) =======================
q-point: ( 0.46  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.907   (   3.700    0.000    0.000)    3.700
   2.907   (   3.700    0.000    0.000)    3.700
   3.664   (   3.819    0.000    0.000)    3.819
   5.440   (   1.669    0.000    0.000)    1.669
   5.440   (   1.669    0.000    0.000)    1.669
   8.280   (   7.821    0.000    0.000)    7.821
   8.320   (   1.658    0.000    0.000)    1.658
   8.320   (   1.658    0.000    0.000)    1.658
   8.562   (   2.229    0.000    0.000)    2.229
  10.725   (  -0.322    0.000    0.000)    0.322
  10.725   (  -0.322    0.000    0.000)    0.322
  14.107   (   3.870    0.000    0.000)    3.870
  14.522   (  -0.334    0.000    0.000)    0.334
  14.522   (  -0.334    0.000    0.000)    0.334
  16.452   (  -5.840    0.000    0.000)    5.840
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.792   (  19.749   19.749    0.000)   27.929
   1.040   (  25.158   25.158    0.000)   35.579
   1.585   (  36.660   36.660    0.000)   51.845
   5.004   (   2.850    2.850    0.000)    4.031
   5.034   (   4.349    4.349    0.000)    6.150
   5.996   (   8.590    8.590    0.000)   12.148
   7.735   (  -2.537   -2.537    0.000)    3.588
   7.841   (   3.055    3.055    0.000)    4.320
   7.977   (  10.005   10.005    0.000)   14.149
  10.808   (  -0.785   -0.785    0.000)    1.110
  10.833   (   0.618    0.618    0.000)    0.875
  13.566   (   7.016    7.016    0.000)    9.922
  14.498   (  -6.948   -6.948    0.000)    9.826
  14.622   (  -0.674   -0.674    0.000)    0.953
  17.613   (  -3.550   -3.550    0.000)    5.020
======================= Grid point 15 (9/84) =======================
q-point: ( 0.15  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.347   (  31.749   -0.632    0.000)   31.755
   1.581   (  28.084   13.809    0.000)   31.296
   2.281   (  33.028   24.019    0.000)   40.839
   5.085   (   5.295    3.808    0.000)    6.522
   5.130   (   5.634    7.510    0.000)    9.389
   6.313   (  22.534    1.964    0.000)   22.620
   7.714   (   0.331   -9.095    0.000)    9.101
   7.926   (   5.456    2.989    0.000)    6.221
   8.207   (  13.276   15.807    0.000)   20.643
  10.786   (  -1.445   -1.563    0.000)    2.128
  10.833   (  -0.463    2.303    0.000)    2.349
  13.662   (   3.083   12.819    0.000)   13.185
  14.431   (  -0.373  -15.736    0.000)   15.740
  14.604   (  -1.186   -0.667    0.000)    1.361
  17.478   (  -9.963   -1.568    0.000)   10.086
======================= Grid point 16 (10/84) =======================
q-point: ( 0.23  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.940   (  27.926   -4.667    0.000)   28.313
   2.095   (  23.873    7.377    0.000)   24.987
   2.847   (  24.574   12.981    0.000)   27.792
   5.205   (   6.836    5.308    0.000)    8.654
   5.256   (   7.178    9.781    0.000)   12.132
   6.818   (  27.318   -3.000    0.000)   27.482
   7.749   (   3.461  -14.138    0.000)   14.556
   8.047   (   6.661    2.899    0.000)    7.265
   8.486   (  14.982   21.764    0.000)   26.422
  10.754   (  -1.848   -2.747    0.000)    3.311
  10.818   (  -0.997    3.084    0.000)    3.241
  13.700   (   1.009    5.208    0.000)    5.305
  14.469   (   4.147  -10.216    0.000)   11.026
  14.578   (  -1.418   -0.655    0.000)    1.562
  17.234   ( -14.717   -1.828    0.000)   14.830
======================= Grid point 17 (11/84) =======================
q-point: ( 0.31  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.418   (  20.669   -3.806    0.000)   21.016
   2.500   (  17.353    3.371    0.000)   17.678
   3.243   (  16.154    4.058    0.000)   16.656
   5.342   (   6.905    7.075    0.000)    9.886
   5.402   (   7.383   12.344    0.000)   14.383
   7.292   (  18.496  -15.882    0.000)   24.379
   7.878   (  10.886  -11.040    0.000)   15.504
   8.176   (   6.357    2.819    0.000)    6.954
   8.780   (  14.650   28.173    0.000)   31.754
  10.717   (  -1.863   -4.134    0.000)    4.534
  10.795   (  -1.386    3.243    0.000)    3.526
  13.707   (  -0.132   -8.804    0.000)    8.805
  14.551   (  -1.320   -0.640    0.000)    1.467
  14.588   (   7.884    2.719    0.000)    8.340
  16.915   ( -17.531   -3.759    0.000)   17.929
======================= Grid point 18 (12/84) =======================
q-point: ( 0.38  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.743   (  12.475   -1.711    0.000)   12.591
   2.770   (  10.299    0.899    0.000)   10.338
   3.488   (   9.059   -2.309    0.000)    9.349
   5.463   (   5.170    8.657    0.000)   10.084
   5.532   (   5.591   14.880    0.000)   15.896
   7.509   (   5.355  -32.703    0.000)   33.139
   8.155   (  14.742   -3.107    0.000)   15.066
   8.285   (   4.560    2.775    0.000)    5.338
   9.043   (  11.703   34.634    0.000)   36.558
  10.684   (  -1.404   -5.401    0.000)    5.580
  10.766   (  -1.516    2.594    0.000)    3.005
  13.699   (  -0.497  -20.972    0.000)   20.978
  14.528   (  -0.926   -0.625    0.000)    1.117
  14.766   (   9.669   16.007    0.000)   18.701
  16.573   ( -16.580   -7.708    0.000)   18.284
======================= Grid point 19 (13/84) =======================
q-point: ( 0.46  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.904   (   3.390   -0.297    0.000)    3.403
   2.905   (   4.147   -0.158    0.000)    4.150
   3.604   (   2.920   -5.690    0.000)    6.395
   5.535   (   1.920    9.590    0.000)    9.781
   5.609   (   2.088   16.507    0.000)   16.639
   7.564   (   1.087  -39.261    0.000)   39.276
   8.348   (   1.657    2.764    0.000)    3.223
   8.367   (   5.822   -1.842    0.000)    6.106
   9.211   (   4.761   39.097    0.000)   39.385
  10.664   (  -0.526   -6.171    0.000)    6.193
  10.741   (  -0.754    1.498    0.000)    1.677
  13.691   (  -0.244  -27.529    0.000)   27.530
  14.516   (  -0.332   -0.616    0.000)    0.700
  14.927   (   5.315   26.307    0.000)   26.839
  16.317   (  -7.844  -13.194    0.000)   15.349
======================= Grid point 28 (14/84) =======================
q-point: ( 0.15  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.521   (  17.144   17.144    0.000)   24.245
   1.916   (  19.041   19.041    0.000)   26.928
   2.743   (  21.621   21.621    0.000)   30.577
   5.193   (   7.143    7.143    0.000)   10.102
   5.298   (   9.089    9.089    0.000)   12.854
   6.483   (  15.559   15.559    0.000)   22.004
   7.564   (  -6.204   -6.204    0.000)    8.774
   8.009   (   5.337    5.337    0.000)    7.547
   8.547   (  18.740   18.740    0.000)   26.502
  10.742   (  -2.918   -2.918    0.000)    4.126
  10.866   (   0.868    0.868    0.000)    1.228
  13.943   (  11.743   11.743    0.000)   16.608
  14.108   ( -12.692  -12.692    0.000)   17.949
  14.585   (  -1.172   -1.172    0.000)    1.657
  17.383   (  -8.455   -8.455    0.000)   11.958
======================= Grid point 29 (15/84) =======================
q-point: ( 0.23  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.931   (  22.880    4.607    0.000)   23.339
   2.282   (  17.622   11.026    0.000)   20.787
   3.100   (  14.747   11.662    0.000)   18.801
   5.356   (   9.302    9.980    0.000)   13.643
   5.484   (  10.037   12.036    0.000)   15.672
   6.846   (  17.752    7.730    0.000)   19.362
   7.499   (   1.538  -11.751    0.000)   11.852
   8.127   (   6.523    5.164    0.000)    8.320
   8.931   (  20.228   23.352    0.000)   30.895
  10.675   (  -3.812   -5.191    0.000)    6.441
  10.869   (  -0.375    1.664    0.000)    1.705
  13.755   ( -14.931   -3.390    0.000)   15.311
  14.306   (  16.892   -2.114    0.000)   17.024
  14.560   (  -1.397   -1.148    0.000)    1.808
  17.152   ( -14.930   -7.026    0.000)   16.500
======================= Grid point 30 (16/84) =======================
q-point: ( 0.31  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.363   (  20.384   -0.954    0.000)   20.406
   2.586   (  13.188    5.103    0.000)   14.141
   3.323   (   8.243    3.459    0.000)    8.939
   5.544   (   9.436   13.260    0.000)   16.275
   5.690   (  10.849   14.702    0.000)   18.271
   7.016   (  -0.328   -8.645    0.000)    8.651
   7.714   (  18.822   -7.024    0.000)   20.090
   8.254   (   6.253    4.993    0.000)    8.002
   9.325   (  19.766   27.086    0.000)   33.531
  10.598   (  -3.942   -7.914    0.000)    8.841
  10.854   (  -1.247    2.134    0.000)    2.471
  13.486   ( -12.277  -13.214    0.000)   18.036
  14.533   (  -1.297   -1.119    0.000)    1.713
  14.647   (  18.052    3.169    0.000)   18.328
  16.811   ( -19.559   -6.668    0.000)   20.664
======================= Grid point 31 (17/84) =======================
q-point: ( 0.38  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.705   (  14.205   -1.733    0.000)   14.310
   2.793   (   7.898    1.276    0.000)    8.000
   3.433   (   3.302   -3.185    0.000)    4.588
   5.709   (   7.074   16.181    0.000)   17.659
   5.894   (   9.487   19.208    0.000)   21.423
   6.942   (  -5.337  -22.415    0.000)   23.041
   8.083   (  16.436   -4.453    0.000)   17.029
   8.362   (   4.519    4.889    0.000)    6.658
   9.685   (  16.359   30.300    0.000)   34.434
  10.528   (  -3.047  -10.479    0.000)   10.913
  10.819   (  -2.312    2.348    0.000)    3.295
  13.282   (  -8.397  -20.805    0.000)   22.435
  14.511   (  -0.908   -1.092    0.000)    1.420
  15.014   (  19.373    7.628    0.000)   20.820
  16.402   ( -21.801   -8.048    0.000)   23.239
======================= Grid point 32 (18/84) =======================
q-point: ( 0.46  0.15  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.895   (   2.607   -0.585    0.000)    2.671
   2.901   (   5.323   -0.257    0.000)    5.329
   3.466   (   0.550   -7.840    0.000)    7.859
   5.807   (   2.628   17.915    0.000)   18.106
   6.035   (   4.211   25.369    0.000)   25.716
   6.852   (  -2.921  -32.342    0.000)   32.474
   8.308   (   5.955   -4.230    0.000)    7.304
   8.424   (   1.652    4.860    0.000)    5.133
   9.927   (   7.170   33.340    0.000)   34.102
  10.485   (  -1.159  -12.083    0.000)   12.139
  10.773   (  -1.772    1.457    0.000)    2.294
  13.168   (  -3.044  -25.389    0.000)   25.571
  14.499   (  -0.325   -1.075    0.000)    1.124
  15.384   (  16.445   17.231    0.000)   23.819
  15.994   ( -17.550  -16.218    0.000)   23.896
======================= Grid point 42 (19/84) =======================
q-point: ( 0.23  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.111   (  12.911   12.911    0.000)   18.258
   2.500   (  10.726   10.726    0.000)   15.168
   3.282   (   6.753    6.753    0.000)    9.550
   5.585   (  13.025   13.025    0.000)   18.420
   5.697   (   9.851    9.851    0.000)   13.932
   7.115   (  15.386   15.386    0.000)   21.759
   7.294   (  -6.265   -6.265    0.000)    8.861
   8.241   (   6.302    6.302    0.000)    8.913
   9.389   (  23.180   23.180    0.000)   32.782
  10.555   (  -6.927   -6.927    0.000)    9.796
  10.881   (  -0.225   -0.225    0.000)    0.318
  13.535   ( -16.205  -16.205    0.000)   22.917
  14.432   (  12.343   12.343    0.000)   17.455
  14.535   (  -1.364   -1.364    0.000)    1.929
  16.956   ( -12.776  -12.776    0.000)   18.068
======================= Grid point 43 (20/84) =======================
q-point: ( 0.31  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.398   (  15.538    4.609    0.000)   16.207
   2.687   (   8.216    5.000    0.000)    9.618
   3.365   (   1.950    0.757    0.000)    2.092
   5.848   (  13.285   17.382    0.000)   21.878
   5.883   (   9.378    4.094    0.000)   10.233
   7.034   ( -10.862    9.822    0.000)   14.644
   7.593   (  22.284   -4.752    0.000)   22.785
   8.364   (   6.062    6.049    0.000)    8.564
   9.826   (  21.307   23.892    0.000)   32.013
  10.413   (  -7.382  -10.835    0.000)   13.111
  10.870   (  -0.816   -0.297    0.000)    0.868
  13.200   ( -17.023  -15.935    0.000)   23.318
  14.509   (  -1.263   -1.325    0.000)    1.830
  14.715   (  15.803    3.869    0.000)   16.269
  16.663   ( -16.861   -8.218    0.000)   18.757
======================= Grid point 44 (21/84) =======================
q-point: ( 0.38  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.689   (  13.567    0.368    0.000)   13.572
   2.817   (   4.976    1.127    0.000)    5.102
   3.365   (  -1.801   -3.811    0.000)    4.215
   6.066   (   9.182   -5.520    0.000)   10.714
   6.082   (  10.072   21.481    0.000)   23.725
   6.822   ( -10.486   12.943    0.000)   16.657
   7.982   (  16.262   -5.799    0.000)   17.265
   8.469   (   4.430    5.864    0.000)    7.349
  10.211   (  17.626   23.197    0.000)   29.134
  10.279   (  -5.909  -14.772    0.000)   15.910
  10.847   (  -1.709    0.632    0.000)    1.822
  12.896   ( -13.368  -18.727    0.000)   23.009
  14.487   (  -0.882   -1.291    0.000)    1.564
  15.034   (  16.526   -4.951    0.000)   17.252
  16.314   ( -18.232   -0.473    0.000)   18.238
======================= Grid point 45 (22/84) =======================
q-point: ( 0.46  0.23  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.882   (   1.651   -0.786    0.000)    1.828
   2.898   (   6.602    0.039    0.000)    6.602
   3.314   (  -2.467   -7.497    0.000)    7.892
   6.221   (   3.782   24.089    0.000)   24.385
   6.223   (   5.731  -20.308    0.000)   21.101
   6.647   (  -6.221   24.447    0.000)   25.226
   8.201   (   5.771   -6.648    0.000)    8.803
   8.530   (   1.641    5.804    0.000)    6.032
  10.196   (  -2.308  -17.405    0.000)   17.558
  10.489   (   9.466   23.434    0.000)   25.274
  10.796   (  -3.056    1.290    0.000)    3.317
  12.706   (  -5.334  -21.891    0.000)   22.532
  14.475   (  -0.316   -1.271    0.000)    1.310
  15.330   (  11.610  -19.361    0.000)   22.575
  15.991   ( -12.422   13.844    0.000)   18.600
======================= Grid point 56 (23/84) =======================
q-point: ( 0.31  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.534   (   8.867    8.867    0.000)   12.540
   2.775   (   3.840    3.840    0.000)    5.430
   3.351   (  -2.182   -2.182    0.000)    3.086
   5.894   (  -1.293   -1.293    0.000)    1.829
   6.202   (  18.153   18.153    0.000)   25.672
   7.229   (   4.239    4.239    0.000)    5.995
   7.606   (   9.115    9.115    0.000)   12.890
   8.482   (   5.789    5.789    0.000)    8.187
  10.185   ( -12.061  -12.061    0.000)   17.057
  10.253   (  19.687   19.687    0.000)   27.841
  10.857   (  -0.655   -0.655    0.000)    0.926
  12.871   ( -17.380  -17.380    0.000)   24.579
  14.483   (  -1.224   -1.224    0.000)    1.731
  14.798   (   4.502    4.502    0.000)    6.366
  16.500   (  -8.126   -8.126    0.000)   11.492
======================= Grid point 57 (24/84) =======================
q-point: ( 0.38  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.726   (  10.211    3.365    0.000)   10.751
   2.835   (   2.328    0.758    0.000)    2.448
   3.283   (  -4.688   -4.551    0.000)    6.534
   5.866   (  -1.277  -10.723    0.000)   10.799
   6.528   (  14.389   23.506    0.000)   27.560
   7.209   (  -3.230   20.931    0.000)   21.179
   7.866   (  13.304   -5.529    0.000)   14.407
   8.583   (   4.268    5.522    0.000)    6.979
   9.960   ( -10.263  -17.418    0.000)   20.217
  10.591   (  14.509   15.963    0.000)   21.572
  10.855   (   0.660    0.508    0.000)    0.833
  12.541   ( -15.621  -17.821    0.000)   23.698
  14.462   (  -0.854   -1.192    0.000)    1.466
  14.883   (   3.890   -8.741    0.000)    9.567
  16.357   (  -6.196    3.777    0.000)    7.256
======================= Grid point 58 (25/84) =======================
q-point: ( 0.46  0.31  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.866   (   0.773   -0.805    0.000)    1.116
   2.906   (   6.974    0.848    0.000)    7.026
   3.180   (  -4.832   -6.005    0.000)    7.708
   5.849   (  -0.446  -15.697    0.000)   15.703
   6.731   (   5.634   27.565    0.000)   28.135
   7.157   (  -1.514   24.383    0.000)   24.430
   8.052   (   4.975   -8.398    0.000)    9.761
   8.642   (   1.621    5.411    0.000)    5.648
   9.812   (  -4.236  -21.586    0.000)   21.998
  10.782   (   3.999    5.460    0.000)    6.768
  10.892   (   2.892   10.239    0.000)   10.640
  12.304   (  -7.186  -19.538    0.000)   20.817
  14.451   (  -0.305   -1.173    0.000)    1.212
  14.939   (   1.596  -17.901    0.000)   17.972
  16.271   (  -2.382   12.898    0.000)   13.116
======================= Grid point 70 (26/84) =======================
q-point: ( 0.38  0.38  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.816   (   5.590    5.590    0.000)    7.905
   2.847   (   0.400    0.400    0.000)    0.566
   3.187   (  -5.260   -5.260    0.000)    7.439
   5.697   (  -6.406   -6.406    0.000)    9.059
   6.969   (  20.528   20.528    0.000)   29.031
   7.516   (   7.993    7.993    0.000)   11.303
   7.809   (   1.443    1.443    0.000)    2.041
   8.678   (   4.014    4.014    0.000)    5.676
   9.621   ( -16.453  -16.453    0.000)   23.268
  10.854   (  10.641   10.641    0.000)   15.049
  10.874   (   2.189    2.189    0.000)    3.095
  12.201   ( -16.492  -16.492    0.000)   23.323
  14.442   (  -0.831   -0.831    0.000)    1.175
  14.738   (  -5.532   -5.532    0.000)    7.824
  16.431   (   3.281    3.281    0.000)    4.640
======================= Grid point 71 (27/84) =======================
q-point: ( 0.46  0.38  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.852   (   0.121   -0.603    0.000)    0.615
   2.931   (   5.537    1.623    0.000)    5.770
   3.081   (  -5.176   -4.196    0.000)    6.663
   5.615   (  -2.058   -8.487    0.000)    8.733
   7.275   (   9.061   27.196    0.000)   28.666
   7.552   (  -0.681   15.380    0.000)   15.396
   7.889   (   3.482   -7.610    0.000)    8.369
   8.734   (   1.579    3.810    0.000)    4.125
   9.368   (  -7.751  -23.132    0.000)   24.396
  10.890   (   0.530    7.305    0.000)    7.324
  11.046   (   6.523    4.748    0.000)    8.068
  11.936   (  -8.748  -17.979    0.000)   19.994
  14.431   (  -0.297   -0.818    0.000)    0.871
  14.669   (  -1.653   -9.562    0.000)    9.704
  16.476   (   1.193    7.950    0.000)    8.039
======================= Grid point 84 (28/84) =======================
q-point: ( 0.46  0.46  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.843   (  -0.224   -0.224    0.000)    0.316
   2.967   (   2.017    2.017    0.000)    2.852
   3.014   (  -2.719   -2.719    0.000)    3.846
   5.508   (  -2.637   -2.637    0.000)    3.729
   7.729   (  15.893   15.893    0.000)   22.476
   7.743   (   3.066    3.066    0.000)    4.336
   7.786   (  -1.243   -1.243    0.000)    1.758
   8.786   (   1.429    1.429    0.000)    2.021
   8.968   ( -14.519  -14.519    0.000)   20.533
  11.050   (   6.121    6.121    0.000)    8.657
  11.105   (   2.788    2.788    0.000)    3.943
  11.628   ( -11.523  -11.523    0.000)   16.296
  14.420   (  -0.293   -0.293    0.000)    0.414
  14.554   (  -2.689   -2.689    0.000)    3.803
  16.580   (   2.662    2.662    0.000)    3.764
======================= Grid point 183 (29/84) =======================
q-point: ( 0.08  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.075   (  17.442   17.442   17.442)   30.211
   1.075   (  17.442   17.442   17.442)   30.211
   1.941   (  27.855   27.855   27.855)   48.247
   5.056   (   3.581    3.581    3.581)    6.203
   5.056   (   3.581    3.581    3.581)    6.203
   6.051   (   8.095    8.095    8.095)   14.020
   7.652   (  -4.389   -4.389   -4.389)    7.602
   8.000   (   7.191    7.191    7.191)   12.455
   8.000   (   7.191    7.191    7.191)   12.455
  10.827   (   0.167    0.167    0.167)    0.289
  10.827   (   0.167    0.167    0.167)    0.289
  13.679   (   8.237    8.237    8.237)   14.266
  14.491   (  -4.857   -4.857   -4.857)    8.412
  14.491   (  -4.857   -4.857   -4.857)    8.412
  17.594   (  -3.234   -3.234   -3.234)    5.602
======================= Grid point 184 (30/84) =======================
q-point: ( 0.15  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.454   (  25.100    8.073    8.073)   27.575
   1.609   (  27.965    8.168    8.168)   30.257
   2.493   (  26.719   17.654   17.654)   36.568
   5.124   (   4.584    3.973    3.973)    7.251
   5.181   (   7.788    4.870    4.870)   10.396
   6.345   (  21.068    5.060    5.060)   22.251
   7.589   (  -2.145   -8.186   -8.186)   11.774
   8.074   (   5.275    8.806    8.806)   13.525
   8.272   (  15.399    8.202    8.202)   19.279
  10.808   (  -1.303    0.371    0.371)    1.404
  10.831   (  -0.444    0.816    0.816)    1.236
  13.781   (   1.064   11.115   11.115)   15.756
  14.396   (   0.523   -8.046   -8.046)   11.390
  14.473   (  -1.191   -6.936   -6.936)    9.881
  17.465   (  -9.856   -2.192   -2.192)   10.332
======================= Grid point 185 (31/84) =======================
q-point: ( 0.23  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.959   (  25.060    1.091    1.091)   25.107
   2.124   (  24.060    5.074    5.074)   25.107
   2.950   (  19.766    8.935    8.935)   23.460
   5.233   (   6.406    3.964    3.964)    8.512
   5.356   (   9.806    8.283    8.283)   15.277
   6.813   (  24.905    0.911    0.911)   24.938
   7.561   (  -0.932  -12.313  -12.313)   17.438
   8.217   (   9.675    9.116    9.116)   16.119
   8.589   (  16.665   11.689   11.689)   23.473
  10.776   (  -1.982   -0.112   -0.112)    1.988
  10.814   (  -1.200    1.198    1.198)    2.076
  13.736   (  -4.489    3.039    3.039)    6.215
  14.447   (  -1.421   -6.917   -6.917)    9.885
  14.497   (   8.731   -2.318   -2.318)    9.326
  17.217   ( -15.220   -2.145   -2.145)   15.520
======================= Grid point 186 (32/84) =======================
q-point: ( 0.31  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.399   (  19.411   -1.473   -1.473)   19.522
   2.534   (  17.635    3.345    3.345)   18.259
   3.267   (  12.832    2.138    2.138)   13.183
   5.367   (   6.972    4.438    4.438)    9.382
   5.550   (   9.695   12.286   12.286)   19.897
   7.248   (  18.188   -5.936   -5.936)   20.032
   7.547   (  -0.580  -16.413  -16.413)   23.219
   8.464   (  15.578   12.170   12.170)   23.214
   8.907   (  15.354   15.104   15.104)   26.306
  10.731   (  -2.509   -1.057   -1.057)    2.920
  10.786   (  -1.612    1.118    1.118)    2.259
  13.640   (  -4.892   -6.662   -6.662)   10.616
  14.419   (  -1.324   -6.894   -6.894)    9.839
  14.708   (  12.406    5.692    5.692)   14.789
  16.881   ( -18.775   -3.300   -3.300)   19.346
======================= Grid point 187 (33/84) =======================
q-point: ( 0.38  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.707   (  11.962   -2.024   -2.024)   12.300
   2.809   (  10.521    2.386    2.386)   11.049
   3.461   (   7.136   -2.525   -2.525)    7.979
   5.492   (   5.536    5.293    5.293)    9.310
   5.718   (   7.098   15.846   15.846)   23.507
   7.506   (   8.603  -12.885  -12.885)   20.151
   7.536   (  -0.498  -19.936  -19.936)   28.198
   8.791   (  16.190   16.106   16.106)   27.945
   9.171   (  11.180   18.296   18.296)   28.187
  10.681   (  -2.419   -2.495   -2.495)    4.278
  10.752   (  -1.805    0.581    0.581)    1.983
  13.557   (  -3.444  -13.491  -13.491)   19.387
  14.397   (  -0.933   -6.876   -6.876)    9.768
  14.971   (  14.028   12.857   12.857)   22.966
  16.503   ( -19.123   -6.289   -6.289)   21.090
======================= Grid point 188 (34/84) =======================
q-point: ( 0.46  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.864   (   4.014   -1.919   -1.919)    4.845
   2.946   (   3.471    1.952    1.952)    4.435
   3.552   (   2.292   -4.944   -4.944)    7.357
   5.570   (   2.126    6.005    6.005)    8.755
   5.816   (   2.604   17.978   17.978)   25.557
   7.529   (  -0.222  -22.058  -22.058)   31.196
   7.608   (   2.346  -16.624  -16.624)   23.627
   9.031   (   6.916   18.494   18.494)   27.053
   9.325   (   4.189   20.547   20.547)   29.358
  10.645   (  -1.081   -3.835   -3.835)    5.530
  10.723   (  -0.946   -0.230   -0.230)    1.001
  13.510   (  -1.227  -17.001  -17.001)   24.075
  14.384   (  -0.336   -6.867   -6.867)    9.717
  15.215   (   8.777   20.964   20.964)   30.919
  16.187   ( -10.745  -12.500  -12.500)   20.687
======================= Grid point 197 (35/84) =======================
q-point: ( 0.15  0.15  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.663   (  14.954   14.954   13.579)   25.132
   1.913   (  19.456   19.456    1.008)   27.533
   2.849   (  17.106   17.106    9.332)   25.930
   5.231   (   6.513    6.513    3.957)   10.025
   5.312   (   8.750    8.750    1.432)   12.456
   6.561   (  16.650   16.650    8.772)   25.128
   7.444   (  -6.563   -6.563   -9.141)   13.028
   8.197   (   5.386    5.386   15.951)   17.676
   8.545   (  17.951   17.951    0.008)   25.387
  10.781   (  -2.338   -2.338    3.151)    4.567
  10.858   (   0.842    0.842   -0.835)    1.454
  14.030   (   8.776    8.776    8.308)   14.935
  14.102   ( -12.649  -12.649   -0.602)   17.898
  14.444   (   0.646    0.646  -13.262)   13.294
  17.363   (  -8.544   -8.544   -2.491)   12.338
======================= Grid point 198 (36/84) =======================
q-point: ( 0.23  0.15  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.015   (  19.687    5.848    8.049)   22.058
   2.295   (  18.631   11.129    1.798)   21.776
   3.133   (  11.752    8.700    2.981)   14.923
   5.367   (   7.446    9.612    1.707)   12.278
   5.529   (  13.125    8.612    3.696)   16.127
   6.922   (  16.935   11.156    8.480)   21.981
   7.342   (  -2.766   -9.989  -12.182)   15.995
   8.333   (   9.267    4.183   16.953)   19.768
   8.926   (  20.356   20.891   -0.133)   29.169
  10.723   (  -3.481   -4.186    4.065)    6.795
  10.859   (  -0.577    1.861   -0.998)    2.189
  13.751   ( -16.304   -4.306   -0.532)   16.872
  14.261   (   6.201   -6.990   -5.883)   11.042
  14.522   (   9.219    3.170   -1.935)    9.939
  17.132   ( -14.819   -6.839   -2.124)   16.458
======================= Grid point 199 (37/84) =======================
q-point: ( 0.31  0.15  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.388   (  17.617    0.997    2.601)   17.836
   2.620   (  14.314    5.113    3.340)   15.563
   3.311   (   6.661    2.010   -1.055)    7.037
   5.522   (   8.171   11.190   -0.046)   13.855
   5.804   (  14.411   11.375    8.526)   20.243
   7.060   (  -2.477   -8.552    3.997)    9.759
   7.416   (   8.726    0.573  -18.885)   20.812
   8.578   (  15.652    2.229   22.790)   27.737
   9.323   (  19.853   24.216    0.141)   31.314
  10.647   (  -4.274   -6.448    4.221)    8.812
  10.838   (  -1.541    2.484   -1.425)    3.252
  13.458   ( -13.405  -11.610   -2.972)   17.981
  14.306   (   0.200   -4.547  -14.429)   15.130
  14.787   (  16.248    2.448    5.758)   17.411
  16.795   ( -19.411   -5.303   -1.412)   20.171
======================= Grid point 200 (38/84) =======================
q-point: ( 0.38  0.15  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.679   (  11.683   -0.614   -1.286)   11.770
   2.850   (   9.229    1.703    4.432)   10.379
   3.402   (   2.795   -3.261   -2.929)    5.199
   5.674   (   6.889   12.818   -0.607)   14.565
   6.071   (  12.276   16.806   13.983)   25.073
   6.959   (  -6.314  -21.160    1.083)   22.109
   7.549   (   4.140    2.253  -29.388)   29.763
   8.902   (  15.759    0.883   31.036)   34.819
   9.684   (  16.472   27.976    0.188)   32.465
  10.561   (  -4.249   -8.828    2.913)   10.221
  10.797   (  -2.665    2.704   -1.981)    4.282
  13.235   (  -9.111  -17.944   -4.902)   20.713
  14.302   (  -0.422   -3.354  -15.668)   16.028
  15.132   (  18.832    1.785    5.598)   19.728
  16.386   ( -21.968   -3.781   -0.329)   22.294
======================= Grid point 201 (39/84) =======================
q-point: ( 0.46  0.15  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.833   (   3.909   -1.154   -2.687)    4.881
   2.977   (   3.504    1.198    4.107)    5.530
   3.431   (   0.537   -6.891   -3.334)    7.674
   5.773   (   2.786   14.315   -0.407)   14.590
   6.256   (   5.654   25.124   18.828)   31.901
   6.851   (  -3.681  -31.661   -0.414)   31.877
   7.592   (   0.810    0.698  -37.359)   37.374
   9.132   (   6.567    0.551   36.854)   37.438
   9.931   (   7.430   32.008    0.575)   32.864
  10.494   (  -2.112  -10.881    0.507)   11.096
  10.745   (  -2.023    1.659   -2.590)    3.681
  13.112   (  -3.273  -21.910   -5.910)   22.928
  14.295   (  -0.193   -2.884  -15.908)   16.169
  15.522   (  21.027    1.134    4.339)   21.500
  15.949   ( -22.372   -2.606    0.871)   22.540
======================= Grid point 211 (40/84) =======================
q-point: ( 0.23  0.23  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.194   (  11.904   11.904    8.138)   18.698
   2.521   (  11.287   11.287    2.106)   16.101
   3.270   (   5.192    5.192   -1.072)    7.421
   5.588   (  11.816   11.816    0.301)   16.713
   5.696   (   9.556    9.556   -0.080)   13.514
   7.156   (  -7.138   -7.138  -11.864)   15.578
   7.228   (  15.955   15.955   11.957)   25.537
   8.426   (   6.272    6.272   16.533)   18.762
   9.357   (  22.513   22.513   -2.993)   31.978
  10.623   (  -5.954   -5.954    6.010)   10.345
  10.874   (  -0.133   -0.133   -0.712)    0.736
  13.532   ( -16.098  -16.098   -0.307)   22.768
  14.221   (   0.925    0.925  -15.078)   15.135
  14.626   (   7.620    7.620    3.038)   11.197
  16.948   ( -11.794  -11.794   -0.742)   16.695
======================= Grid point 212 (41/84) =======================
q-point: ( 0.31  0.23  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.450   (  13.488    5.376    5.172)   15.413
   2.721   (   8.992    4.920    3.301)   10.768
   3.336   (   1.616    0.398   -2.800)    3.257
   5.785   (   8.731   14.361   -3.547)   17.176
   5.937   (  13.909    3.113    2.377)   14.450
   6.998   (  -7.583    0.842   -7.156)   10.460
   7.517   (  10.536    9.607    1.035)   14.296
   8.609   (  12.776    1.696   20.181)   23.945
   9.788   (  21.121   22.505   -3.606)   31.074
  10.496   (  -6.892   -8.682    7.446)   13.353
  10.863   (  -0.921    0.092   -0.722)    1.174
  13.194   ( -17.319  -15.214   -0.627)   23.061
  14.241   (   0.785   -2.165  -19.252)   19.389
  14.811   (  10.818    0.227    1.804)   10.970
  16.684   ( -14.795   -5.767    2.477)   16.071
======================= Grid point 213 (42/84) =======================
q-point: ( 0.38  0.23  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.690   (  10.180    1.843    0.710)   10.370
   2.874   (   6.895    0.684    5.058)    8.579
   3.338   (  -1.329   -3.278   -2.577)    4.377
   5.952   (   8.427   11.249   -5.711)   15.172
   6.175   (   8.992   -2.712    3.531)   10.034
   6.878   (  -4.063   11.731    2.160)   12.601
   7.620   (   0.880    5.722  -18.203)   19.101
   8.897   (  15.003   -0.982   28.734)   32.430
  10.161   (  15.507   16.950   -3.476)   23.234
  10.370   (  -4.468   -7.201    7.259)   11.158
  10.836   (  -1.993    1.253   -1.097)    2.597
  12.885   ( -13.619  -17.687   -1.269)   22.359
  14.251   (   0.242   -2.040  -19.625)   19.732
  15.029   (  10.933  -10.797   -4.974)   16.151
  16.392   ( -14.302    3.849    8.283)   16.969
======================= Grid point 214 (43/84) =======================
q-point: ( 0.46  0.23  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.822   (   3.162    0.086   -2.268)    3.892
   2.991   (   4.248    0.175    5.614)    7.042
   3.298   (  -1.963   -6.447   -1.430)    6.889
   6.101   (   5.914   14.813   -6.796)   17.337
   6.252   (  -0.502  -17.977   -0.213)   17.985
   6.850   (   0.428   27.089   17.498)   32.252
   7.602   (  -1.369    0.373  -34.938)   34.967
   9.127   (   6.877   -0.921   35.619)   36.288
  10.224   (  -2.631  -15.464    2.639)   15.907
  10.485   (   9.209   22.686   -0.365)   24.486
  10.778   (  -3.557    2.035   -1.873)    4.506
  12.691   (  -5.412  -20.756   -1.640)   21.512
  14.253   (   0.051   -1.626  -19.316)   19.385
  15.205   (   5.705  -23.373  -13.992)   27.832
  16.168   (  -7.078   15.880   16.270)   23.812
======================= Grid point 225 (44/84) =======================
q-point: ( 0.31  0.31  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.589   (   8.376    8.376    5.465)   13.046
   2.807   (   3.839    3.839    3.134)    6.268
   3.324   (  -1.642   -1.642   -2.574)    3.467
   5.889   (  -1.424   -1.424   -0.579)    2.096
   6.145   (  16.321   16.321   -5.768)   23.791
   7.032   (   2.305    2.305  -17.404)   17.707
   7.730   (   9.227    9.227   13.182)   18.548
   8.679   (   6.481    6.481   17.315)   19.592
  10.151   (   8.991    8.991   -2.790)   13.018
  10.355   (  -0.098   -0.098    8.757)    8.758
  10.855   (  -0.518   -0.518   -0.202)    0.759
  12.873   ( -17.250  -17.250    0.188)   24.396
  14.218   (  -0.386   -0.386  -21.197)   21.204
  14.804   (  -0.727   -0.727   -4.853)    4.961
  16.577   (  -4.993   -4.993    7.784)   10.510
======================= Grid point 226 (45/84) =======================
q-point: ( 0.38  0.31  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.752   (   7.315    4.469    2.788)    9.014
   2.879   (   4.140   -0.101    4.182)    5.886
   3.272   (  -3.718   -3.528   -0.998)    5.222
   5.855   (  -1.511  -10.780   -1.104)   10.941
   6.429   (  11.154   20.560   -9.468)   25.234
   7.117   (   7.472    8.914  -14.200)   18.356
   7.790   (  -3.445   10.874    4.168)   12.145
   8.877   (  12.774   -0.448   22.845)   26.178
  10.050   (  -9.865  -14.286    8.343)   19.262
  10.583   (  13.984   15.522   -0.763)   20.906
  10.854   (   0.681    0.848   -0.110)    1.093
  12.542   ( -15.750  -17.436    0.121)   23.497
  14.220   (   0.394   -1.112  -21.066)   21.099
  14.788   (  -0.885  -12.242  -12.027)   17.184
  16.494   (  -3.340    5.667   12.934)   14.511
======================= Grid point 227 (46/84) =======================
q-point: ( 0.46  0.31  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.836   (   1.580    1.350   -0.676)    2.185
   2.986   (   5.586   -0.523    5.794)    8.065
   3.186   (  -4.322   -4.904    0.684)    6.572
   5.835   (  -0.519  -15.988   -1.409)   16.058
   6.560   (   2.777   23.396  -13.535)   27.171
   7.325   (  12.093   16.945    2.941)   21.024
   7.635   ( -10.536    5.731  -16.469)   20.373
   9.100   (   7.627   -1.688   31.123)   32.088
   9.881   (  -5.850  -18.978    6.746)   20.974
  10.771   (   3.994    5.876   -1.144)    7.196
  10.893   (   2.996   10.761    0.041)   11.170
  12.302   (  -7.305  -19.173   -0.203)   20.519
  14.227   (   0.263   -1.060  -20.272)   20.301
  14.774   (  -0.398  -19.125  -16.878)   25.511
  16.450   (  -1.171   12.377   16.124)   20.360
======================= Grid point 239 (47/84) =======================
q-point: ( 0.38  0.38  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.856   (   5.323    5.323    4.065)    8.555
   2.875   (   0.395    0.395    2.808)    2.863
   3.197   (  -4.227   -4.227    1.087)    6.075
   5.681   (  -6.681   -6.681   -1.596)    9.582
   6.819   (  17.202   17.202  -14.763)   28.456
   7.239   (   5.991    5.991  -23.701)   25.169
   7.946   (   2.399    2.399   14.642)   15.030
   8.917   (   5.376    5.376   19.994)   21.391
   9.761   ( -14.616  -14.616   13.030)   24.434
  10.850   (  11.134   11.134   -0.496)   15.753
  10.877   (   2.289    2.289    0.212)    3.245
  12.206   ( -16.508  -16.508    0.497)   23.352
  14.208   (  -0.097   -0.097  -21.068)   21.068
  14.588   (  -7.830   -7.830  -16.225)   19.644
  16.592   (   3.972    3.972   14.698)   15.734
======================= Grid point 240 (48/84) =======================
q-point: ( 0.46  0.38  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.867   (  -0.162    1.426    1.882)    2.367
   2.979   (   5.771    0.138    4.477)    7.305
   3.106   (  -4.883   -3.423    2.582)    6.498
   5.595   (  -2.162   -8.753   -1.983)    9.232
   7.016   (   3.373   22.179  -20.458)   30.361
   7.429   (  11.256   -1.262  -27.053)   29.328
   7.872   (  -7.364   12.985   16.713)   22.409
   9.068   (   7.693   -1.290   25.167)   26.348
   9.507   (  -9.737  -18.312   13.147)   24.555
  10.891   (   0.500    7.889   -0.005)    7.905
  11.056   (   7.036    5.214    0.923)    8.806
  11.938   (  -8.979  -18.015    0.168)   20.129
  14.211   (   0.251   -0.552  -20.330)   20.339
  14.485   (  -2.655  -10.359  -18.707)   21.548
  16.644   (   1.371    7.466   15.293)   17.073
======================= Grid point 253 (49/84) =======================
q-point: ( 0.46  0.46  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.883   (   0.280    0.280    3.953)    3.973
   3.000   (   1.925    1.925    3.377)    4.338
   3.047   (  -2.659   -2.659    3.356)    5.040
   5.484   (  -2.738   -2.738   -2.354)    4.532
   7.354   (   7.628    7.628  -27.087)   29.156
   7.405   (   2.180    2.180  -27.983)   28.153
   8.030   (   3.154    3.154   17.680)   18.234
   9.070   (   2.113    2.113   22.632)   22.829
   9.209   ( -10.961  -10.961   20.356)   25.586
  11.062   (   6.646    6.646    1.165)    9.471
  11.123   (   3.197    3.197    1.786)    4.862
  11.624   ( -12.049  -12.049   -0.350)   17.043
  14.205   (  -0.104   -0.104  -20.095)   20.096
  14.358   (  -3.032   -3.032  -19.782)   20.241
  16.742   (   2.480    2.480   14.754)   15.165
======================= Grid point 366 (50/84) =======================
q-point: ( 0.15  0.15  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.999   (  13.714   13.714   13.714)   23.754
   1.999   (  13.714   13.714   13.714)   23.754
   3.043   (   9.397    9.397    9.397)   16.277
   5.357   (   6.319    6.319    6.319)   10.945
   5.357   (   6.319    6.319    6.319)   10.945
   6.835   (  18.299   18.299   18.299)   31.695
   7.266   (  -8.680   -8.680   -8.680)   15.034
   8.560   (  11.201   11.201   11.201)   19.402
   8.560   (  11.201   11.201   11.201)   19.402
  10.832   (  -0.097   -0.097   -0.097)    0.168
  10.832   (  -0.097   -0.097   -0.097)    0.168
  14.086   (  -8.698   -8.698   -8.698)   15.065
  14.086   (  -8.698   -8.698   -8.698)   15.065
  14.277   (  10.407   10.407   10.407)   18.026
  17.266   (  -7.703   -7.703   -7.703)   13.342
======================= Grid point 367 (51/84) =======================
q-point: ( 0.23  0.15  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.210   (  12.764   10.706   10.706)   19.803
   2.379   (  18.689    7.463    7.463)   21.464
   3.196   (   6.139    3.030    3.030)    7.487
   5.445   (   5.601    6.539    6.539)   10.812
   5.600   (  15.304    3.516    3.516)   16.092
   7.105   (  -7.821  -11.235  -11.235)   17.709
   7.195   (  15.855   18.621   18.621)   30.738
   8.661   (   9.539   13.113   13.113)   20.853
   8.973   (  21.421    7.882    7.882)   24.148
  10.785   (  -3.096    0.617    0.617)    3.216
  10.841   (  -0.520    0.029    0.029)    0.522
  13.716   ( -19.137   -4.064   -4.064)   19.981
  14.063   (  -1.224  -12.344  -12.344)   17.500
  14.538   (  14.311    3.255    3.255)   15.033
  17.062   ( -13.015   -5.214   -5.214)   14.959
======================= Grid point 368 (52/84) =======================
q-point: ( 0.31  0.15  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.465   (  12.529    5.207    5.207)   14.533
   2.707   (  14.475    5.225    5.225)   16.252
   3.287   (   3.342   -1.267   -1.267)    3.792
   5.573   (   7.335    5.446    5.446)   10.636
   5.939   (  18.768    3.967    3.967)   19.589
   6.957   (  -7.327  -13.771  -13.771)   20.807
   7.375   (   1.880   15.264   15.264)   21.668
   8.960   (  20.324   12.504   12.504)   26.939
   9.383   (  20.175    8.919    8.919)   23.794
  10.708   (  -4.746    0.276    0.276)    4.762
  10.818   (  -1.819    0.366    0.366)    1.891
  13.363   ( -16.542   -6.996   -6.996)   19.275
  14.039   (  -1.132  -12.141  -12.141)   17.208
  14.819   (  14.193   -2.482   -2.482)   14.621
  16.770   ( -16.314   -0.761   -0.761)   16.350
======================= Grid point 369 (53/84) =======================
q-point: ( 0.38  0.15  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.678   (   8.837    1.293    1.293)    9.024
   2.937   (   8.867    3.939    3.939)   10.471
   3.333   (   1.484   -3.647   -3.647)    5.367
   5.719   (   7.076    5.183    5.183)   10.187
   6.315   (  19.243    7.029    7.029)   21.659
   6.822   (  -6.226  -16.305  -16.305)   23.884
   7.294   (  -8.945    8.467    8.467)   14.947
   9.389   (  21.486   14.072   14.072)   29.286
   9.748   (  16.775   10.285   10.285)   22.203
  10.598   (  -6.489   -0.741   -0.741)    6.573
  10.768   (  -3.284    0.223    0.223)    3.299
  13.085   ( -11.499  -10.418  -10.418)   18.690
  14.020   (  -0.797  -11.968  -11.968)   16.944
  15.088   (  13.014  -10.703  -10.703)   19.962
  16.446   ( -16.116    6.837    6.837)   18.795
======================= Grid point 370 (54/84) =======================
q-point: ( 0.46  0.15  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.798   (   3.189   -0.800   -0.800)    3.384
   3.052   (   2.956    3.358    3.358)    5.593
   3.350   (   0.372   -4.563   -4.563)    6.463
   5.826   (   3.207    5.837    5.837)    8.856
   6.668   (  15.506   17.152   17.152)   28.788
   6.729   (  -2.758  -18.473  -18.473)   26.270
   7.069   ( -12.397   -3.876   -3.876)   13.555
   9.718   (  10.050   16.322   16.322)   25.175
  10.013   (   8.922   13.605   13.605)   21.208
  10.471   (  -5.311   -3.927   -3.927)    7.684
  10.701   (  -2.705   -0.798   -0.798)    2.931
  12.929   (  -4.157  -12.750  -12.750)   18.503
  14.009   (  -0.287  -11.869  -11.869)   16.788
  15.299   (   7.061  -22.448  -22.448)   32.522
  16.187   (  -8.486   18.236   18.236)   27.149
======================= Grid point 380 (55/84) =======================
q-point: ( 0.23  0.23  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.407   (   9.613    9.613   12.947)   18.774
   2.581   (  12.028   12.028    4.042)   17.484
   3.243   (   1.741    1.741   -1.635)    2.956
   5.609   (   8.216    8.216    2.464)   11.878
   5.696   (   8.478    8.478    0.201)   11.991
   6.902   (  -9.447   -9.447  -13.106)   18.715
   7.580   (  17.678   17.678   23.094)   34.034
   8.808   (   6.309    6.309   21.310)   23.103
   9.283   (  20.608   20.608   -4.148)   29.438
  10.744   (  -3.697   -3.697    5.370)    7.494
  10.853   (   0.153    0.153   -1.516)    1.531
  13.523   ( -15.801  -15.801   -0.511)   22.352
  13.892   (  -4.180   -4.180  -18.276)   19.208
  14.623   (   4.805    4.805   -3.058)    7.452
  16.932   (  -7.845   -7.845   -0.903)   11.131
======================= Grid point 381 (56/84) =======================
q-point: ( 0.31  0.23  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.588   (   8.653    7.286    8.552)   14.181
   2.809   (  10.694    4.816    5.473)   12.943
   3.267   (   0.788   -0.704   -3.951)    4.090
   5.733   (   5.976   10.487   -1.338)   12.144
   5.953   (  15.364   -1.730   -0.719)   15.478
   6.738   (  -7.067   -7.915  -16.133)   19.310
   7.773   (   0.848   22.498   22.504)   31.832
   9.037   (  17.736   -0.817   22.287)   28.495
   9.697   (  20.856   19.882   -5.438)   29.322
  10.653   (  -5.556   -4.776    7.491)   10.479
  10.842   (  -1.158    0.979   -1.360)    2.037
  13.173   ( -18.481  -12.639   -1.722)   22.456
  13.852   (  -0.878   -6.853  -19.695)   20.872
  14.711   (   3.910   -7.778  -11.724)   14.603
  16.769   (  -8.380    0.495    6.015)   10.327
======================= Grid point 382 (57/84) =======================
q-point: ( 0.38  0.23  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.731   (   5.642    4.109    3.547)    7.829
   2.990   (   7.684    1.235    6.353)   10.047
   3.275   (  -0.037   -2.047   -3.537)    4.087
   5.860   (   6.952    7.821   -3.443)   11.016
   6.210   (   8.002  -12.795    0.351)   15.095
   6.669   (   2.391    0.120  -21.502)   21.635
   7.654   ( -11.218   24.075   19.656)   33.043
   9.438   (  20.960   -3.957   25.867)   33.527
  10.077   (  17.237   17.393   -5.400)   25.076
  10.530   (  -6.594   -4.367    8.397)   11.535
  10.806   (  -2.750    2.555   -1.924)    4.218
  12.840   ( -14.654  -14.145   -3.500)   20.666
  13.842   (  -0.214   -6.538  -20.439)   21.460
  14.769   (   2.015  -19.665  -20.835)   28.720
  16.621   (  -6.296    9.879   13.916)   18.190
======================= Grid point 383 (58/84) =======================
q-point: ( 0.46  0.23  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.804   (   1.839    1.503    0.553)    2.439
   3.102   (   3.377    1.554    5.515)    6.651
   3.266   (  -0.608   -3.732   -1.603)    4.107
   5.988   (   5.111    9.271   -4.415)   11.470
   6.229   (  -2.891  -22.691   -2.403)   23.001
   6.838   (  10.273    1.583  -27.714)   29.599
   7.411   ( -10.806   27.821   24.418)   38.563
   9.780   (  11.704   -3.017   30.369)   32.686
  10.275   (  -0.227   -6.232    1.473)    6.408
  10.482   (   4.862   17.165    0.540)   17.848
  10.727   (  -4.821    3.483   -3.075)    6.695
  12.633   (  -5.799  -16.833   -4.701)   18.415
  13.840   (  -0.023   -5.850  -20.852)   21.657
  14.792   (   0.502  -27.239  -27.043)   38.387
  16.535   (  -2.306   16.337   19.446)   25.502
======================= Grid point 394 (59/84) =======================
q-point: ( 0.31  0.31  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.736   (   7.167    7.167    9.173)   13.670
   2.890   (   3.788    3.788    5.206)    7.469
   3.260   (  -0.129   -0.129   -3.722)    3.726
   5.872   (  -1.707   -1.707   -1.174)    2.684
   5.995   (  11.256   11.256   -8.621)   18.103
   6.645   (  -1.634   -1.634  -20.625)   20.754
   8.126   (   9.885    9.885   26.627)   30.073
   9.066   (   6.922    6.922   20.898)   23.077
  10.072   (  17.008   17.008   -5.400)   24.652
  10.550   (  -4.905   -4.905   10.111)   12.262
  10.847   (  -0.156   -0.156   -0.581)    0.621
  12.878   ( -16.889  -16.889    0.362)   23.887
  13.767   (  -2.187   -2.187  -23.714)   23.915
  14.544   (  -8.604   -8.604  -20.995)   24.266
  16.784   (   0.913    0.913   12.387)   12.454
======================= Grid point 395 (60/84) =======================
q-point: ( 0.38  0.31  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.833   (   2.433    6.154    5.332)    8.498
   2.986   (   6.277   -1.292    6.437)    9.084
   3.252   (  -1.036   -0.662   -0.847)    1.493
   5.823   (  -2.195  -10.657   -2.198)   11.100
   6.199   (   7.370   11.631  -12.849)   18.833
   6.688   (   6.375    1.024  -25.537)   26.340
   8.121   (  -7.778   21.022   26.076)   34.386
   9.336   (  18.190   -4.779   22.938)   29.663
  10.231   (  -3.690   -3.965    8.119)    9.760
  10.572   (   9.310   11.030    0.347)   14.438
  10.849   (   0.623    1.755   -0.542)    1.939
  12.544   ( -16.320  -16.058   -0.047)   22.895
  13.751   (   0.212   -3.040  -25.171)   25.355
  14.392   (  -6.655  -17.580  -27.522)   33.329
  16.800   (   0.696    8.007   16.809)   18.631
======================= Grid point 396 (61/84) =======================
q-point: ( 0.46  0.31  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.855   (   0.306    3.509    2.543)    4.344
   3.109   (   5.436   -0.774    6.367)    8.408
   3.212   (  -2.976   -1.855    2.057)    4.066
   5.794   (  -0.741  -16.603   -2.799)   16.853
   6.274   (   1.278   14.496  -14.906)   20.831
   6.842   (   6.335   -0.736  -32.082)   32.710
   7.950   (  -6.593   25.324   28.761)   38.884
   9.679   (  14.599   -6.760   26.631)   31.113
  10.062   ( -10.299  -12.312   11.696)   19.861
  10.735   (   4.080    7.069   -2.555)    8.552
  10.888   (   3.037   12.085   -0.779)   12.485
  12.292   (  -7.773  -17.771   -1.026)   19.424
  13.761   (   0.476   -2.670  -25.501)   25.645
  14.298   (  -2.632  -22.074  -30.910)   38.074
  16.810   (   0.248   11.627   18.987)   22.266
======================= Grid point 408 (62/84) =======================
q-point: ( 0.38  0.38  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.950   (   0.551    0.551    4.657)    4.722
   2.968   (   4.653    4.653    7.055)    9.647
   3.232   (  -1.623   -1.623    2.593)    3.463
   5.636   (  -7.444   -7.444   -2.912)   10.922
   6.447   (  10.981   10.981  -21.566)   26.576
   6.719   (   3.068    3.068  -27.748)   28.085
   8.375   (   3.664    3.664   28.276)   28.747
   9.337   (   6.485    6.485   21.545)   23.416
  10.064   ( -11.003  -11.003   16.508)   22.685
  10.829   (  12.465   12.465   -1.861)   17.726
  10.880   (   2.524    2.524    0.024)    3.569
  12.221   ( -16.506  -16.506    0.978)   23.363
  13.718   (  -0.411   -0.411  -27.355)   27.361
  14.109   ( -11.099  -11.099  -31.741)   35.410
  16.926   (   4.831    4.831   17.932)   19.189
======================= Grid point 409 (63/84) =======================
q-point: ( 0.46  0.38  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.932   (  -0.600    3.996    4.746)    6.233
   3.088   (   5.721   -0.821    6.287)    8.541
   3.181   (  -3.705   -1.653    5.014)    6.450
   5.540   (  -2.456   -9.454   -3.609)   10.413
   6.564   (   1.807   14.013  -24.619)   28.385
   6.818   (   4.798   -1.416  -32.277)   32.662
   8.343   (  -3.136   13.840   29.474)   32.712
   9.551   (  13.053   -4.992   22.838)   26.775
   9.828   ( -12.671  -11.180   18.532)   25.079
  10.886   (   0.415    9.548   -0.569)    9.573
  11.077   (   8.618    6.665    1.027)   10.943
  11.942   (  -9.736  -17.828    0.286)   20.315
  13.726   (   0.724   -1.062  -27.596)   27.626
  13.958   (  -4.183  -12.426  -34.218)   36.644
  16.989   (   1.587    6.752   18.302)   19.572
======================= Grid point 422 (64/84) =======================
q-point: ( 0.46  0.46  0.15)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.990   (   1.515    1.515    6.778)    7.109
   3.094   (   1.699    1.699    6.016)    6.478
   3.141   (  -2.473   -2.473    5.992)    6.938
   5.419   (  -3.015   -3.015   -4.229)    6.005
   6.768   (   4.474    4.474  -31.285)   31.918
   6.802   (   1.022    1.022  -32.021)   32.054
   8.491   (   2.505    2.505   28.931)   29.148
   9.526   (   2.632    2.632   22.593)   22.897
   9.637   (  -8.119   -8.119   21.782)   24.623
  11.090   (   8.062    8.062    1.459)   11.494
  11.170   (   4.508    4.508    2.888)    6.999
  11.618   ( -13.427  -13.427   -0.148)   18.990
  13.717   (   0.027    0.027  -28.287)   28.287
  13.802   (  -3.966   -3.966  -36.141)   36.573
  17.076   (   2.194    2.194   17.823)   18.091
======================= Grid point 549 (65/84) =======================
q-point: ( 0.23  0.23  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.673   (   9.395    9.395    9.395)   16.272
   2.673   (   9.395    9.395    9.395)   16.272
   3.212   (  -1.394   -1.394   -1.394)    2.415
   5.709   (   4.586    4.586    4.586)    7.943
   5.709   (   4.586    4.586    4.586)    7.943
   6.659   ( -11.601  -11.601  -11.601)   20.094
   8.052   (  21.929   21.929   21.929)   37.981
   9.230   (  11.004   11.004   11.004)   19.059
   9.230   (  11.004   11.004   11.004)   19.059
  10.814   (  -0.401   -0.401   -0.401)    0.695
  10.814   (  -0.401   -0.401   -0.401)    0.695
  13.513   ( -10.444  -10.444  -10.444)   18.090
  13.513   ( -10.444  -10.444  -10.444)   18.090
  14.527   (  -6.167   -6.167   -6.167)   10.681
  16.915   (  -0.792   -0.792   -0.792)    1.371
======================= Grid point 550 (66/84) =======================
q-point: ( 0.31  0.23  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.768   (   4.118    9.339    9.339)   13.834
   2.921   (  11.283    5.495    5.495)   13.700
   3.197   (   0.082   -2.697   -2.697)    3.815
   5.747   (   2.898    3.056    3.056)    5.203
   5.915   (  12.673   -2.939   -2.939)   13.337
   6.448   ( -10.044  -12.583  -12.583)   20.434
   8.291   (   0.777   27.895   27.895)   39.457
   9.412   (  18.912   10.473   10.473)   24.022
   9.630   (  19.197    4.605    4.605)   20.271
  10.751   (  -3.920    1.498    1.498)    4.456
  10.819   (  -0.969   -0.223   -0.223)    1.019
  13.108   ( -20.548   -5.975   -5.975)   22.217
  13.502   (  -0.495  -14.931  -14.931)   21.122
  14.399   (  -6.639  -18.603  -18.603)   27.134
  16.901   (  -0.607    6.967    6.967)    9.872
======================= Grid point 551 (67/84) =======================
q-point: ( 0.38  0.23  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.825   (   1.765    5.888    5.888)    8.512
   3.104   (   7.166    4.426    4.426)    9.514
   3.215   (   1.524   -1.804   -1.804)    2.972
   5.824   (   4.862   -0.454   -0.454)    4.904
   6.177   (  13.141   -8.594   -8.594)   17.900
   6.264   (  -8.840  -12.358  -12.358)   19.586
   8.166   ( -10.355   29.457   29.457)   42.927
   9.871   (  24.677   10.419   10.419)   28.742
   9.972   (  15.566    2.658    2.658)   16.013
  10.658   (  -5.641    2.918    2.918)    6.989
  10.779   (  -3.306    0.422    0.422)    3.360
  12.735   ( -16.771   -7.633   -7.633)   19.945
  13.494   (  -0.332  -14.549  -14.549)   20.578
  14.278   (  -5.345  -27.460  -27.460)   39.200
  16.892   (  -0.349   12.955   12.955)   18.325
======================= Grid point 552 (68/84) =======================
q-point: ( 0.46  0.23  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.844   (   0.350    3.426    3.426)    4.858
   3.198   (   2.423    3.935    3.935)    6.069
   3.242   (   0.883   -0.651   -0.651)    1.275
   5.931   (   5.605   -2.114   -2.114)    6.352
   6.106   (  -6.852  -12.481  -12.481)   18.934
   6.381   (   6.312  -13.398  -13.398)   19.971
   7.985   (  -5.962   31.258   31.258)   44.605
  10.223   (   9.040   -0.764   -0.764)    9.104
  10.284   (  15.077   12.218   12.218)   22.932
  10.539   (  -5.601    5.977    5.977)   10.140
  10.683   (  -5.866    0.405    0.405)    5.893
  12.493   (  -6.906   -9.847   -9.847)   15.544
  13.489   (  -0.117  -14.328  -14.328)   20.263
  14.202   (  -2.129  -31.715  -31.715)   44.903
  16.887   (  -0.105   16.048   16.048)   22.696
======================= Grid point 563 (69/84) =======================
q-point: ( 0.31  0.31  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.931   (   5.783    5.783   10.295)   13.148
   3.003   (   3.607    3.607    6.131)    7.976
   3.188   (   2.031    2.031   -3.543)    4.562
   5.842   (  -1.774   -1.774   -1.880)    3.135
   5.849   (   4.310    4.310   -5.531)    8.230
   6.261   (  -6.396   -6.396  -17.945)   20.096
   8.727   (  11.723   11.723   33.393)   37.282
   9.465   (   6.474    6.474   19.444)   21.492
   9.937   (  19.174   19.174   -8.218)   28.334
  10.717   (  -3.610   -3.610    6.480)    8.250
  10.830   (   0.303    0.303   -1.262)    1.333
  12.887   ( -16.375  -16.375    0.486)   23.162
  13.316   (  -4.356   -4.356  -22.092)   22.934
  14.041   ( -16.849  -16.849  -28.735)   37.329
  17.025   (   5.476    5.476   11.588)   13.937
======================= Grid point 564 (70/84) =======================
q-point: ( 0.38  0.31  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.952   (  -1.074    6.786    6.572)    9.508
   3.119   (   6.480   -1.501    6.912)    9.593
   3.244   (   2.646    3.367    0.097)    4.283
   5.768   (  -3.300   -9.099   -3.539)   10.306
   5.968   (   4.271   -1.546   -9.940)   10.929
   6.205   (   0.838   -0.969  -22.868)   22.903
   8.716   (  -7.940   24.469   33.094)   41.917
   9.759   (  18.993   -7.382   20.244)   28.724
  10.282   (  14.038   14.790   -3.264)   20.651
  10.641   (  -2.949   -1.499    6.338)    7.149
  10.829   (   0.032    2.919   -1.545)    3.303
  12.535   ( -17.935  -12.806   -1.116)   22.066
  13.279   (  -0.393   -7.744  -22.307)   23.617
  13.765   ( -11.066  -23.437  -34.949)   43.511
  17.111   (   3.251    9.265   14.314)   17.358
======================= Grid point 565 (71/84) =======================
q-point: ( 0.46  0.31  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.933   (  -0.588    5.520    5.365)    7.720
   3.221   (   3.578   -0.346    4.789)    5.988
   3.270   (  -0.108    2.157    3.852)    4.416
   5.723   (  -1.174  -16.264   -4.673)   16.962
   6.002   (   0.235    2.901  -12.338)   12.677
   6.260   (   2.725   -2.974  -26.322)   26.629
   8.569   (  -4.838   27.000   33.259)   43.112
  10.086   (  12.928   -8.617   14.323)   21.131
  10.353   (  -5.317   -5.167   17.778)   19.262
  10.676   (   3.764    7.377   -3.094)    8.841
  10.852   (   2.297   13.844   -3.153)   14.383
  12.252   (  -8.914  -14.350   -3.360)   17.224
  13.282   (   0.353   -7.393  -22.501)   23.687
  13.615   (  -4.073  -26.827  -37.360)   46.175
  17.153   (   1.057   11.018   15.533)   19.073
======================= Grid point 577 (72/84) =======================
q-point: ( 0.38  0.38  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.053   (   1.263    1.263    5.786)    6.056
   3.121   (   3.868    3.868    8.262)    9.909
   3.293   (   1.330    1.330    3.408)    3.893
   5.570   (  -8.519   -8.519   -3.716)   12.608
   6.032   (   4.645    4.645  -19.852)   20.911
   6.195   (   0.168    0.168  -25.193)   25.194
   9.011   (   4.031    4.031   35.690)   36.142
   9.733   (   6.716    6.716   18.652)   20.930
  10.365   (  -5.714   -5.714   13.403)   15.650
  10.770   (  13.398   13.398   -4.254)   19.420
  10.874   (   2.739    2.739   -0.798)    3.956
  12.244   ( -16.382  -16.382    1.348)   23.207
  13.175   (  -3.133   -3.133  -28.019)   28.368
  13.402   ( -13.498  -13.498  -38.557)   43.023
  17.252   (   5.267    5.267   14.813)   16.581
======================= Grid point 578 (73/84) =======================
q-point: ( 0.46  0.38  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.048   (  -0.410    5.789    6.930)    9.039
   3.216   (   4.169   -0.057    6.655)    7.853
   3.292   (  -1.364    0.008    5.915)    6.070
   5.458   (  -2.887  -10.276   -4.586)   11.617
   6.091   (   1.181    5.191  -22.933)   23.542
   6.215   (   1.345   -1.785  -28.698)   28.785
   8.984   (  -2.665   14.194   35.005)   37.868
   9.944   (  11.581   -4.637   17.482)   21.477
  10.197   (  -9.799   -9.307   18.341)   22.782
  10.863   (   0.395   12.014   -1.904)   12.171
  11.086   (  11.291    9.258   -0.446)   14.608
  11.947   ( -11.262  -16.605    0.060)   20.063
  13.152   (   0.708   -9.944  -35.714)   37.080
  13.221   (  -5.520   -9.274  -34.137)   35.803
  17.320   (   1.677    6.155   14.924)   16.230
======================= Grid point 591 (74/84) =======================
q-point: ( 0.46  0.46  0.23)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.138   (   2.718    2.718    8.082)    8.950
   3.227   (   1.439    1.439    7.261)    7.540
   3.271   (  -1.968   -1.968    6.995)    7.528
   5.325   (  -3.388   -3.388   -5.211)    7.079
   6.171   (   1.997    1.997  -28.985)   29.123
   6.193   (  -0.178   -0.178  -29.562)   29.563
   9.129   (   2.272    2.272   35.322)   35.468
   9.929   (   2.667    2.667   18.365)   18.749
  10.029   (  -7.090   -7.090   18.031)   20.631
  11.111   (   9.725    9.725    0.519)   13.763
  11.228   (   7.027    7.027    2.728)   10.306
  11.623   ( -14.880  -14.880    0.771)   21.057
  12.977   (  -4.722   -4.722  -46.112)   46.593
  13.156   (  -1.095   -1.095  -28.882)   28.924
  17.398   (   1.944    1.944   14.491)   14.750
======================= Grid point 732 (75/84) =======================
q-point: ( 0.31  0.31  0.31)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.123   (   6.090    6.090    6.090)   10.547
   3.123   (   6.090    6.090    6.090)   10.547
   3.126   (   0.073    0.073    0.073)    0.126
   5.791   (  -2.175   -2.175   -2.175)    3.767
   5.791   (  -2.175   -2.175   -2.175)    3.767
   5.966   ( -11.405  -11.405  -11.405)   19.754
   9.328   (  21.218   21.218   21.218)   36.751
   9.815   (   9.060    9.060    9.060)   15.693
   9.815   (   9.060    9.060    9.060)   15.693
  10.796   (  -0.259   -0.259   -0.259)    0.449
  10.796   (  -0.259   -0.259   -0.259)    0.449
  12.897   ( -10.377  -10.377  -10.377)   17.974
  12.897   ( -10.377  -10.377  -10.377)   17.974
  13.493   ( -25.469  -25.469  -25.469)   44.113
  17.221   (   8.282    8.282    8.282)   14.345
======================= Grid point 733 (76/84) =======================
q-point: ( 0.38  0.31  0.31)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.080   (  -1.783    6.172    6.172)    8.909
   3.231   (   4.627    2.950    2.950)    6.230
   3.270   (   5.900    3.937    3.937)    8.113
   5.678   (  -6.433   -6.391   -6.391)   11.094
   5.816   (   2.421   -9.264   -9.264)   13.323
   5.837   (  -2.997   -9.866   -9.866)   14.271
   9.350   (  -8.935   29.675   29.675)   42.908
  10.092   (  11.848    4.448    4.448)   13.414
  10.203   (  27.201    6.610    6.610)   28.763
  10.742   (  -3.150    1.062    1.062)    3.490
  10.799   (  -0.811    0.417    0.417)    1.003
  12.479   ( -20.781   -5.550   -5.550)   22.214
  12.906   (   0.392  -15.352  -15.352)   21.715
  13.094   ( -14.918  -31.481  -31.481)   46.953
  17.347   (   4.670    9.781    9.781)   14.600
======================= Grid point 734 (77/84) =======================
q-point: ( 0.46  0.31  0.31)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.057   (  -0.617    6.981    6.981)    9.891
   3.295   (   1.741    2.981    2.981)    4.561
   3.348   (   2.027    3.606    3.606)    5.487
   5.587   (  -2.418  -10.106  -10.106)   14.495
   5.804   (  -0.736   -8.059   -8.059)   11.421
   5.860   (   1.496  -13.010  -13.010)   18.460
   9.199   (  -4.652   29.662   29.662)   42.205
  10.264   (   5.104    3.783    3.783)    7.394
  10.614   (  10.611    4.411    4.411)   12.309
  10.684   (  -2.156    1.246    1.246)    2.785
  10.804   (   2.571    5.996    5.996)    8.861
  12.144   ( -11.001   -8.035   -8.035)   15.816
  12.904   (  -4.793  -33.456  -33.456)   47.556
  12.911   (   0.148  -15.084  -15.084)   21.332
  17.407   (   1.476   10.391   10.391)   14.768
======================= Grid point 746 (78/84) =======================
q-point: ( 0.38  0.38  0.31)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.171   (   2.572    2.572    6.111)    7.112
   3.278   (   3.200    3.200    7.479)    8.742
   3.356   (   3.569    3.569    2.825)    5.784
   5.494   (  -9.568   -9.568   -4.106)   14.140
   5.697   (  -1.675   -1.675  -14.051)   14.249
   5.756   (  -2.620   -2.620  -19.090)   19.447
   9.729   (   3.607    3.607   36.627)   36.980
  10.062   (   6.460    6.460   15.023)   17.583
  10.531   (  14.096   14.096   -0.103)   19.935
  10.700   (   0.916    0.916    0.739)    1.491
  10.845   (   2.801    2.801   -2.241)    4.551
  12.271   ( -16.101  -16.101    1.436)   22.816
  12.559   ( -17.631  -17.631  -35.176)   43.117
  12.749   (  -5.169   -5.169  -25.995)   27.003
  17.494   (   5.341    5.341    9.838)   12.404
======================= Grid point 747 (79/84) =======================
q-point: ( 0.46  0.38  0.31)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.193   (   0.296    6.597    7.654)   10.109
   3.340   (   2.254    1.675    5.792)    6.437
   3.398   (   0.824    1.382    4.710)    4.978
   5.365   (  -3.382  -10.440   -4.925)   12.029
   5.691   (   0.369   -3.303  -17.733)   18.041
   5.719   (  -1.167   -3.944  -21.362)   21.755
   9.665   (  -3.688   16.425   33.085)   37.121
  10.252   (   7.728   -2.707   14.655)   16.787
  10.519   (  -4.862  -10.520   13.953)   18.138
  10.810   (   1.066   14.421   -3.370)   14.847
  11.046   (  14.317   12.287   -3.802)   19.246
  11.936   ( -13.785  -12.989   -1.643)   19.011
  12.353   (  -4.681  -21.930  -39.023)   45.007
  12.700   (  -0.882   -7.089  -23.170)   24.246
  17.561   (   1.651    5.570    9.741)   11.341
======================= Grid point 760 (80/84) =======================
q-point: ( 0.46  0.46  0.31)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.296   (   3.238    3.238    7.760)    9.011
   3.366   (   1.213    1.213    6.689)    6.905
   3.401   (  -0.976   -0.976    5.974)    6.132
   5.223   (  -3.742   -3.742   -5.036)    7.305
   5.658   (  -0.407   -0.407  -22.914)   22.922
   5.669   (  -1.413   -1.413  -23.504)   23.589
   9.831   (   2.305    2.305   35.401)   35.550
  10.240   (   2.091    2.091   13.418)   13.740
  10.337   (  -7.053   -7.053   13.497)   16.782
  11.107   (  10.916   10.916   -0.915)   15.465
  11.262   (  11.484   11.484   -0.033)   16.241
  11.648   ( -15.434  -15.434    1.668)   21.891
  12.063   (  -7.452   -7.452  -44.793)   46.016
  12.614   (  -2.040   -2.040  -26.291)   26.448
  17.631   (   1.697    1.697    9.292)    9.597
======================= Grid point 915 (81/84) =======================
q-point: ( 0.38  0.38  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.278   (   4.440    4.440    4.440)    7.690
   3.402   (   3.502    3.502    3.502)    6.065
   3.402   (   3.502    3.502    3.502)    6.065
   5.389   (  -7.839   -7.839   -7.839)   13.577
   5.487   (  -7.266   -7.266   -7.266)   12.584
   5.487   (  -7.266   -7.266   -7.266)   12.584
  10.322   (   8.006    8.006    8.006)   13.866
  10.322   (   8.006    8.006    8.006)   13.866
  10.470   (  16.338   16.338   16.338)   28.297
  10.783   (   0.360    0.360    0.360)    0.624
  10.783   (   0.360    0.360    0.360)    0.624
  11.958   ( -23.549  -23.549  -23.549)   40.788
  12.297   ( -10.136  -10.136  -10.136)   17.556
  12.297   ( -10.136  -10.136  -10.136)   17.556
  17.640   (   5.175    5.175    5.175)    8.964
======================= Grid point 916 (82/84) =======================
q-point: ( 0.46  0.38  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.334   (   1.288    6.420    6.420)    9.170
   3.433   (   1.023    3.541    3.541)    5.111
   3.469   (   1.742    2.475    2.475)    3.909
   5.258   (  -3.779   -6.615   -6.615)   10.090
   5.393   (  -2.252  -10.314  -10.314)   14.759
   5.422   (  -1.268  -10.328  -10.328)   14.661
  10.193   (  -5.298   18.961   18.961)   27.333
  10.516   (   7.894    7.375    7.375)   13.080
  10.687   (  -5.427    1.046    1.046)    5.625
  10.781   (   2.926    5.688    5.688)    8.560
  10.979   (  14.302    0.599    0.599)   14.327
  11.699   (  -1.696  -21.710  -21.710)   30.750
  11.904   ( -17.802   -7.112   -7.112)   20.447
  12.313   (   0.421  -15.529  -15.529)   21.965
  17.704   (   1.548    4.989    4.989)    7.223
======================= Grid point 929 (83/84) =======================
q-point: ( 0.46  0.46  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.431   (   3.034    3.034    5.760)    7.183
   3.477   (   1.018    1.018    4.453)    4.680
   3.495   (   0.277    0.277    3.505)    3.527
   5.135   (  -3.911   -3.911   -3.707)    6.659
   5.290   (  -2.560   -2.560  -14.831)   15.266
   5.291   (  -2.608   -2.608  -14.366)   14.832
  10.431   (   1.614    1.614   19.654)   19.786
  10.453   (   0.752    0.752    8.429)    8.496
  10.605   (  -5.959   -5.959   19.259)   21.022
  11.080   (   7.356    7.356  -24.022)   26.177
  11.080   (  12.003   12.003   -1.642)   17.054
  11.450   (  -1.702   -1.702  -10.151)   10.432
  11.682   ( -15.215  -15.215    1.626)   21.579
  12.132   (  -3.527   -3.527  -23.131)   23.663
  17.765   (   1.454    1.454    4.592)    5.032
======================= Grid point 1098 (84/84) =======================
q-point: ( 0.46  0.46  0.46)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.46e-03 1.46e-03 1.46e-03 3.36e-05 1.56e-04 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.513   (   2.523    2.523    2.523)    4.370
   3.536   (   1.073    1.073    1.073)    1.859
   3.536   (   1.073    1.073    1.073)    1.859
   5.078   (  -2.669   -2.669   -2.669)    4.622
   5.103   (  -4.250   -4.250   -4.250)    7.361
   5.103   (  -4.250   -4.250   -4.250)    7.361
  10.567   (   0.434    0.434    0.434)    0.751
  10.567   (   0.434    0.434    0.434)    0.751
  10.691   (  -7.435   -7.435   -7.435)   12.877
  11.053   (   8.570    8.570    8.570)   14.844
  11.053   (   8.570    8.570    8.570)   14.844
  11.420   (   4.451    4.451    4.451)    7.710
  11.705   (  -9.738   -9.738   -9.738)   16.867
  11.705   (  -9.738   -9.738   -9.738)   16.867
  17.820   (   1.288    1.288    1.288)    2.230
=================== 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/32955
   10.0   1835.590   1835.590   1835.590     -0.000     -0.000      0.000 3/32955
   20.0    581.817    581.817    581.817     -0.000     -0.000      0.000 3/32955
   30.0    340.234    340.234    340.234     -0.000     -0.000      0.000 3/32955
   40.0    205.350    205.350    205.350     -0.000     -0.000      0.000 3/32955
   50.0    120.130    120.130    120.130     -0.000     -0.000      0.000 3/32955
   60.0     74.006     74.006     74.006     -0.000     -0.000      0.000 3/32955
   70.0     49.382     49.382     49.382      0.000     -0.000      0.000 3/32955
   80.0     35.492     35.492     35.492      0.000     -0.000      0.000 3/32955
   90.0     27.117     27.117     27.117      0.000     -0.000      0.000 3/32955
  100.0     21.743     21.743     21.743      0.000     -0.000      0.000 3/32955
  110.0     18.103     18.103     18.103      0.000     -0.000      0.000 3/32955
  120.0     15.520     15.520     15.520      0.000     -0.000      0.000 3/32955
  130.0     13.611     13.611     13.611      0.000     -0.000      0.000 3/32955
  140.0     12.153     12.153     12.153      0.000     -0.000      0.000 3/32955
  150.0     11.005     11.005     11.005      0.000     -0.000      0.000 3/32955
  160.0     10.078     10.078     10.078      0.000     -0.000      0.000 3/32955
  170.0      9.314      9.314      9.314      0.000     -0.000      0.000 3/32955
  180.0      8.673      8.673      8.673      0.000     -0.000      0.000 3/32955
  190.0      8.125      8.125      8.125      0.000     -0.000      0.000 3/32955
  200.0      7.651      7.651      7.651      0.000     -0.000      0.000 3/32955
  210.0      7.237      7.237      7.237      0.000     -0.000      0.000 3/32955
  220.0      6.870      6.870      6.870      0.000     -0.000      0.000 3/32955
  230.0      6.543      6.543      6.543      0.000     -0.000      0.000 3/32955
  240.0      6.249      6.249      6.249      0.000     -0.000      0.000 3/32955
  250.0      5.983      5.983      5.983      0.000     -0.000      0.000 3/32955
  260.0      5.741      5.741      5.741      0.000     -0.000      0.000 3/32955
  270.0      5.519      5.519      5.519      0.000     -0.000      0.000 3/32955
  280.0      5.316      5.316      5.316      0.000     -0.000      0.000 3/32955
  290.0      5.127      5.127      5.127      0.000     -0.000      0.000 3/32955
  300.0      4.953      4.953      4.953      0.000     -0.000      0.000 3/32955
  310.0      4.791      4.791      4.791      0.000     -0.000      0.000 3/32955
  320.0      4.639      4.639      4.639      0.000     -0.000      0.000 3/32955
  330.0      4.498      4.498      4.498      0.000     -0.000      0.000 3/32955
  340.0      4.365      4.365      4.365      0.000     -0.000      0.000 3/32955
  350.0      4.240      4.240      4.240      0.000     -0.000      0.000 3/32955
  360.0      4.122      4.122      4.122      0.000     -0.000      0.000 3/32955
  370.0      4.011      4.011      4.011      0.000     -0.000      0.000 3/32955
  380.0      3.906      3.906      3.906      0.000     -0.000      0.000 3/32955
  390.0      3.807      3.807      3.807      0.000     -0.000      0.000 3/32955
  400.0      3.712      3.712      3.712      0.000     -0.000      0.000 3/32955
  410.0      3.623      3.623      3.623      0.000     -0.000      0.000 3/32955
  420.0      3.537      3.537      3.537      0.000     -0.000      0.000 3/32955
  430.0      3.456      3.456      3.456      0.000     -0.000      0.000 3/32955
  440.0      3.378      3.378      3.378      0.000     -0.000      0.000 3/32955
  450.0      3.304      3.304      3.304      0.000     -0.000      0.000 3/32955
  460.0      3.233      3.233      3.233      0.000     -0.000      0.000 3/32955
  470.0      3.165      3.165      3.165      0.000     -0.000      0.000 3/32955
  480.0      3.100      3.100      3.100      0.000     -0.000      0.000 3/32955
  490.0      3.038      3.038      3.038      0.000     -0.000      0.000 3/32955
  500.0      2.978      2.978      2.978      0.000     -0.000      0.000 3/32955
  510.0      2.920      2.920      2.920      0.000     -0.000      0.000 3/32955
  520.0      2.865      2.865      2.865      0.000     -0.000      0.000 3/32955
  530.0      2.812      2.812      2.812      0.000     -0.000      0.000 3/32955
  540.0      2.761      2.761      2.761      0.000     -0.000      0.000 3/32955
  550.0      2.711      2.711      2.711      0.000     -0.000      0.000 3/32955
  560.0      2.663      2.663      2.663      0.000     -0.000      0.000 3/32955
  570.0      2.617      2.617      2.617      0.000     -0.000      0.000 3/32955
  580.0      2.573      2.573      2.573      0.000     -0.000      0.000 3/32955
  590.0      2.530      2.530      2.530      0.000     -0.000      0.000 3/32955
  600.0      2.489      2.489      2.489      0.000     -0.000      0.000 3/32955
  610.0      2.448      2.448      2.448      0.000     -0.000      0.000 3/32955
  620.0      2.409      2.409      2.409      0.000     -0.000      0.000 3/32955
  630.0      2.372      2.372      2.372      0.000     -0.000      0.000 3/32955
  640.0      2.335      2.335      2.335      0.000     -0.000      0.000 3/32955
  650.0      2.300      2.300      2.300      0.000     -0.000      0.000 3/32955
  660.0      2.266      2.266      2.266      0.000     -0.000      0.000 3/32955
  670.0      2.232      2.232      2.232      0.000     -0.000      0.000 3/32955
  680.0      2.200      2.200      2.200      0.000     -0.000      0.000 3/32955
  690.0      2.169      2.169      2.169      0.000     -0.000      0.000 3/32955
  700.0      2.138      2.138      2.138      0.000     -0.000      0.000 3/32955
  710.0      2.108      2.108      2.108      0.000     -0.000      0.000 3/32955
  720.0      2.079      2.079      2.079      0.000     -0.000      0.000 3/32955
  730.0      2.051      2.051      2.051      0.000     -0.000      0.000 3/32955
  740.0      2.024      2.024      2.024      0.000     -0.000      0.000 3/32955
  750.0      1.997      1.997      1.997      0.000     -0.000      0.000 3/32955
  760.0      1.971      1.971      1.971      0.000     -0.000      0.000 3/32955
  770.0      1.946      1.946      1.946      0.000     -0.000      0.000 3/32955
  780.0      1.922      1.922      1.922      0.000     -0.000      0.000 3/32955
  790.0      1.898      1.898      1.898      0.000     -0.000      0.000 3/32955
  800.0      1.874      1.874      1.874      0.000     -0.000      0.000 3/32955
  810.0      1.851      1.851      1.851      0.000     -0.000      0.000 3/32955
  820.0      1.829      1.829      1.829      0.000     -0.000      0.000 3/32955
  830.0      1.807      1.807      1.807      0.000     -0.000      0.000 3/32955
  840.0      1.786      1.786      1.786      0.000     -0.000      0.000 3/32955
  850.0      1.765      1.765      1.765      0.000     -0.000      0.000 3/32955
  860.0      1.745      1.745      1.745      0.000     -0.000      0.000 3/32955
  870.0      1.725      1.725      1.725      0.000     -0.000      0.000 3/32955
  880.0      1.706      1.706      1.706      0.000     -0.000      0.000 3/32955
  890.0      1.687      1.687      1.687      0.000     -0.000      0.000 3/32955
  900.0      1.669      1.669      1.669      0.000     -0.000      0.000 3/32955
  910.0      1.650      1.650      1.650      0.000     -0.000      0.000 3/32955
  920.0      1.633      1.633      1.633      0.000     -0.000      0.000 3/32955
  930.0      1.615      1.615      1.615      0.000     -0.000      0.000 3/32955
  940.0      1.598      1.598      1.598      0.000     -0.000      0.000 3/32955
  950.0      1.582      1.582      1.582      0.000     -0.000      0.000 3/32955
  960.0      1.565      1.565      1.565      0.000     -0.000      0.000 3/32955
  970.0      1.550      1.550      1.550      0.000     -0.000      0.000 3/32955
  980.0      1.534      1.534      1.534      0.000     -0.000      0.000 3/32955
  990.0      1.519      1.519      1.519      0.000     -0.000      0.000 3/32955
 1000.0      1.504      1.504      1.504      0.000     -0.000      0.000 3/32955

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

