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

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

-------------------------[time 2026-01-09 04:49: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.837396259999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.837396259999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.837396259999999
Atomic positions (fractional):
   *1 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941
    2 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941
    3 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941
   *4 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453
   *5 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.837396259999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.837396259999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.837396259999999
Atomic positions (fractional):
   *1 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1
    2 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 2
    3 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 3
   *4 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 4
   *5 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   11.512188779999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.512188779999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.512188779999997
Atomic positions (fractional):
   *1 Li  0.00000000000000  0.16666666666667  0.16666666666667   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.16666666666667   6.941 > 1
    3 Li  0.66666666666667  0.16666666666667  0.16666666666667   6.941 > 1
    4 Li  0.00000000000000  0.50000000000000  0.16666666666667   6.941 > 1
    5 Li  0.33333333333333  0.50000000000000  0.16666666666667   6.941 > 1
    6 Li  0.66666666666667  0.50000000000000  0.16666666666667   6.941 > 1
    7 Li  0.00000000000000  0.83333333333333  0.16666666666667   6.941 > 1
    8 Li  0.33333333333333  0.83333333333333  0.16666666666667   6.941 > 1
    9 Li  0.66666666666667  0.83333333333333  0.16666666666667   6.941 > 1
   10 Li  0.00000000000000  0.16666666666667  0.50000000000000   6.941 > 1
   11 Li  0.33333333333333  0.16666666666667  0.50000000000000   6.941 > 1
   12 Li  0.66666666666667  0.16666666666667  0.50000000000000   6.941 > 1
   13 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1
   14 Li  0.33333333333333  0.50000000000000  0.50000000000000   6.941 > 1
   15 Li  0.66666666666667  0.50000000000000  0.50000000000000   6.941 > 1
   16 Li  0.00000000000000  0.83333333333333  0.50000000000000   6.941 > 1
   17 Li  0.33333333333333  0.83333333333333  0.50000000000000   6.941 > 1
   18 Li  0.66666666666667  0.83333333333333  0.50000000000000   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.83333333333333   6.941 > 1
   20 Li  0.33333333333333  0.16666666666667  0.83333333333333   6.941 > 1
   21 Li  0.66666666666667  0.16666666666667  0.83333333333333   6.941 > 1
   22 Li  0.00000000000000  0.50000000000000  0.83333333333333   6.941 > 1
   23 Li  0.33333333333333  0.50000000000000  0.83333333333333   6.941 > 1
   24 Li  0.66666666666667  0.50000000000000  0.83333333333333   6.941 > 1
   25 Li  0.00000000000000  0.83333333333333  0.83333333333333   6.941 > 1
   26 Li  0.33333333333333  0.83333333333333  0.83333333333333   6.941 > 1
   27 Li  0.66666666666667  0.83333333333333  0.83333333333333   6.941 > 1
   28 Li  0.16666666666667  0.00000000000000  0.16666666666667   6.941 > 2
   29 Li  0.50000000000000  0.00000000000000  0.16666666666667   6.941 > 2
   30 Li  0.83333333333333  0.00000000000000  0.16666666666667   6.941 > 2
   31 Li  0.16666666666667  0.33333333333333  0.16666666666667   6.941 > 2
   32 Li  0.50000000000000  0.33333333333333  0.16666666666667   6.941 > 2
   33 Li  0.83333333333333  0.33333333333333  0.16666666666667   6.941 > 2
   34 Li  0.16666666666667  0.66666666666667  0.16666666666667   6.941 > 2
   35 Li  0.50000000000000  0.66666666666667  0.16666666666667   6.941 > 2
   36 Li  0.83333333333333  0.66666666666667  0.16666666666667   6.941 > 2
   37 Li  0.16666666666667  0.00000000000000  0.50000000000000   6.941 > 2
   38 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 2
   39 Li  0.83333333333333  0.00000000000000  0.50000000000000   6.941 > 2
   40 Li  0.16666666666667  0.33333333333333  0.50000000000000   6.941 > 2
   41 Li  0.50000000000000  0.33333333333333  0.50000000000000   6.941 > 2
   42 Li  0.83333333333333  0.33333333333333  0.50000000000000   6.941 > 2
   43 Li  0.16666666666667  0.66666666666667  0.50000000000000   6.941 > 2
   44 Li  0.50000000000000  0.66666666666667  0.50000000000000   6.941 > 2
   45 Li  0.83333333333333  0.66666666666667  0.50000000000000   6.941 > 2
   46 Li  0.16666666666667  0.00000000000000  0.83333333333333   6.941 > 2
   47 Li  0.50000000000000  0.00000000000000  0.83333333333333   6.941 > 2
   48 Li  0.83333333333333  0.00000000000000  0.83333333333333   6.941 > 2
   49 Li  0.16666666666667  0.33333333333333  0.83333333333333   6.941 > 2
   50 Li  0.50000000000000  0.33333333333333  0.83333333333333   6.941 > 2
   51 Li  0.83333333333333  0.33333333333333  0.83333333333333   6.941 > 2
   52 Li  0.16666666666667  0.66666666666667  0.83333333333333   6.941 > 2
   53 Li  0.50000000000000  0.66666666666667  0.83333333333333   6.941 > 2
   54 Li  0.83333333333333  0.66666666666667  0.83333333333333   6.941 > 2
   55 Li  0.16666666666667  0.16666666666667  0.00000000000000   6.941 > 3
   56 Li  0.50000000000000  0.16666666666667  0.00000000000000   6.941 > 3
   57 Li  0.83333333333333  0.16666666666667  0.00000000000000   6.941 > 3
   58 Li  0.16666666666667  0.50000000000000  0.00000000000000   6.941 > 3
   59 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 3
   60 Li  0.83333333333333  0.50000000000000  0.00000000000000   6.941 > 3
   61 Li  0.16666666666667  0.83333333333333  0.00000000000000   6.941 > 3
   62 Li  0.50000000000000  0.83333333333333  0.00000000000000   6.941 > 3
   63 Li  0.83333333333333  0.83333333333333  0.00000000000000   6.941 > 3
   64 Li  0.16666666666667  0.16666666666667  0.33333333333333   6.941 > 3
   65 Li  0.50000000000000  0.16666666666667  0.33333333333333   6.941 > 3
   66 Li  0.83333333333333  0.16666666666667  0.33333333333333   6.941 > 3
   67 Li  0.16666666666667  0.50000000000000  0.33333333333333   6.941 > 3
   68 Li  0.50000000000000  0.50000000000000  0.33333333333333   6.941 > 3
   69 Li  0.83333333333333  0.50000000000000  0.33333333333333   6.941 > 3
   70 Li  0.16666666666667  0.83333333333333  0.33333333333333   6.941 > 3
   71 Li  0.50000000000000  0.83333333333333  0.33333333333333   6.941 > 3
   72 Li  0.83333333333333  0.83333333333333  0.33333333333333   6.941 > 3
   73 Li  0.16666666666667  0.16666666666667  0.66666666666667   6.941 > 3
   74 Li  0.50000000000000  0.16666666666667  0.66666666666667   6.941 > 3
   75 Li  0.83333333333333  0.16666666666667  0.66666666666667   6.941 > 3
   76 Li  0.16666666666667  0.50000000000000  0.66666666666667   6.941 > 3
   77 Li  0.50000000000000  0.50000000000000  0.66666666666667   6.941 > 3
   78 Li  0.83333333333333  0.50000000000000  0.66666666666667   6.941 > 3
   79 Li  0.16666666666667  0.83333333333333  0.66666666666667   6.941 > 3
   80 Li  0.50000000000000  0.83333333333333  0.66666666666667   6.941 > 3
   81 Li  0.83333333333333  0.83333333333333  0.66666666666667   6.941 > 3
  *82 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 4
   83 Cl  0.33333333333333  0.00000000000000  0.00000000000000  35.453 > 4
   84 Cl  0.66666666666667  0.00000000000000  0.00000000000000  35.453 > 4
   85 Cl  0.00000000000000  0.33333333333333  0.00000000000000  35.453 > 4
   86 Cl  0.33333333333333  0.33333333333333  0.00000000000000  35.453 > 4
   87 Cl  0.66666666666667  0.33333333333333  0.00000000000000  35.453 > 4
   88 Cl  0.00000000000000  0.66666666666667  0.00000000000000  35.453 > 4
   89 Cl  0.33333333333333  0.66666666666667  0.00000000000000  35.453 > 4
   90 Cl  0.66666666666667  0.66666666666667  0.00000000000000  35.453 > 4
   91 Cl  0.00000000000000  0.00000000000000  0.33333333333333  35.453 > 4
   92 Cl  0.33333333333333  0.00000000000000  0.33333333333333  35.453 > 4
   93 Cl  0.66666666666667  0.00000000000000  0.33333333333333  35.453 > 4
   94 Cl  0.00000000000000  0.33333333333333  0.33333333333333  35.453 > 4
   95 Cl  0.33333333333333  0.33333333333333  0.33333333333333  35.453 > 4
   96 Cl  0.66666666666667  0.33333333333333  0.33333333333333  35.453 > 4
   97 Cl  0.00000000000000  0.66666666666667  0.33333333333333  35.453 > 4
   98 Cl  0.33333333333333  0.66666666666667  0.33333333333333  35.453 > 4
   99 Cl  0.66666666666667  0.66666666666667  0.33333333333333  35.453 > 4
  100 Cl  0.00000000000000  0.00000000000000  0.66666666666667  35.453 > 4
  101 Cl  0.33333333333333  0.00000000000000  0.66666666666667  35.453 > 4
  102 Cl  0.66666666666667  0.00000000000000  0.66666666666667  35.453 > 4
  103 Cl  0.00000000000000  0.33333333333333  0.66666666666667  35.453 > 4
  104 Cl  0.33333333333333  0.33333333333333  0.66666666666667  35.453 > 4
  105 Cl  0.66666666666667  0.33333333333333  0.66666666666667  35.453 > 4
  106 Cl  0.00000000000000  0.66666666666667  0.66666666666667  35.453 > 4
  107 Cl  0.33333333333333  0.66666666666667  0.66666666666667  35.453 > 4
  108 Cl  0.66666666666667  0.66666666666667  0.66666666666667  35.453 > 4
 *109 O   0.16666666666667  0.16666666666667  0.16666666666667  15.999 > 5
  110 O   0.50000000000000  0.16666666666667  0.16666666666667  15.999 > 5
  111 O   0.83333333333333  0.16666666666667  0.16666666666667  15.999 > 5
  112 O   0.16666666666667  0.50000000000000  0.16666666666667  15.999 > 5
  113 O   0.50000000000000  0.50000000000000  0.16666666666667  15.999 > 5
  114 O   0.83333333333333  0.50000000000000  0.16666666666667  15.999 > 5
  115 O   0.16666666666667  0.83333333333333  0.16666666666667  15.999 > 5
  116 O   0.50000000000000  0.83333333333333  0.16666666666667  15.999 > 5
  117 O   0.83333333333333  0.83333333333333  0.16666666666667  15.999 > 5
  118 O   0.16666666666667  0.16666666666667  0.50000000000000  15.999 > 5
  119 O   0.50000000000000  0.16666666666667  0.50000000000000  15.999 > 5
  120 O   0.83333333333333  0.16666666666667  0.50000000000000  15.999 > 5
  121 O   0.16666666666667  0.50000000000000  0.50000000000000  15.999 > 5
  122 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 5
  123 O   0.83333333333333  0.50000000000000  0.50000000000000  15.999 > 5
  124 O   0.16666666666667  0.83333333333333  0.50000000000000  15.999 > 5
  125 O   0.50000000000000  0.83333333333333  0.50000000000000  15.999 > 5
  126 O   0.83333333333333  0.83333333333333  0.50000000000000  15.999 > 5
  127 O   0.16666666666667  0.16666666666667  0.83333333333333  15.999 > 5
  128 O   0.50000000000000  0.16666666666667  0.83333333333333  15.999 > 5
  129 O   0.83333333333333  0.16666666666667  0.83333333333333  15.999 > 5
  130 O   0.16666666666667  0.50000000000000  0.83333333333333  15.999 > 5
  131 O   0.50000000000000  0.50000000000000  0.83333333333333  15.999 > 5
  132 O   0.83333333333333  0.50000000000000  0.83333333333333  15.999 > 5
  133 O   0.16666666666667  0.83333333333333  0.83333333333333  15.999 > 5
  134 O   0.50000000000000  0.83333333333333  0.83333333333333  15.999 > 5
  135 O   0.83333333333333  0.83333333333333  0.83333333333333  15.999 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.2441248    0.0000000    0.0000000
            0.0000000    3.2441248    0.0000000
            0.0000000    0.0000000    3.2441248
-------------------------- Born effective charges --------------------------
    1 Li    0.9906835    0.0000000    0.0000000
            0.0000000    1.1502198    0.0000000
            0.0000000    0.0000000    1.1502198
    2 Li    1.1502198    0.0000000    0.0000000
            0.0000000    0.9906835    0.0000000
            0.0000000    0.0000000    1.1502198
    3 Li    1.1502198    0.0000000    0.0000000
            0.0000000    1.1502198    0.0000000
            0.0000000    0.0000000    0.9906835
    4 Cl   -1.2699045    0.0000000    0.0000000
            0.0000000   -1.2699045    0.0000000
            0.0000000    0.0000000   -1.2699045
    5 O    -2.0212186    0.0000000    0.0000000
            0.0000000   -2.0212186    0.0000000
            0.0000000    0.0000000   -2.0212186
----------------------------------------------------------------------------
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.088
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.092
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
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 (yy) -0.00000000 (yy) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-09 04:49: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.837396259999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.837396259999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.837396259999999
Atomic positions (fractional):
    1 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941
    2 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941
    3 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941
    4 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453
    5 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.512188779999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.512188779999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.512188779999997
Atomic positions (fractional):
    1 Li  0.00000000000000  0.16666666666667  0.16666666666667   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.16666666666667   6.941 > 1
    3 Li  0.66666666666667  0.16666666666667  0.16666666666667   6.941 > 1
    4 Li  0.00000000000000  0.50000000000000  0.16666666666667   6.941 > 1
    5 Li  0.33333333333333  0.50000000000000  0.16666666666667   6.941 > 1
    6 Li  0.66666666666667  0.50000000000000  0.16666666666667   6.941 > 1
    7 Li  0.00000000000000  0.83333333333333  0.16666666666667   6.941 > 1
    8 Li  0.33333333333333  0.83333333333333  0.16666666666667   6.941 > 1
    9 Li  0.66666666666667  0.83333333333333  0.16666666666667   6.941 > 1
   10 Li  0.00000000000000  0.16666666666667  0.50000000000000   6.941 > 1
   11 Li  0.33333333333333  0.16666666666667  0.50000000000000   6.941 > 1
   12 Li  0.66666666666667  0.16666666666667  0.50000000000000   6.941 > 1
   13 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1
   14 Li  0.33333333333333  0.50000000000000  0.50000000000000   6.941 > 1
   15 Li  0.66666666666667  0.50000000000000  0.50000000000000   6.941 > 1
   16 Li  0.00000000000000  0.83333333333333  0.50000000000000   6.941 > 1
   17 Li  0.33333333333333  0.83333333333333  0.50000000000000   6.941 > 1
   18 Li  0.66666666666667  0.83333333333333  0.50000000000000   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.83333333333333   6.941 > 1
   20 Li  0.33333333333333  0.16666666666667  0.83333333333333   6.941 > 1
   21 Li  0.66666666666667  0.16666666666667  0.83333333333333   6.941 > 1
   22 Li  0.00000000000000  0.50000000000000  0.83333333333333   6.941 > 1
   23 Li  0.33333333333333  0.50000000000000  0.83333333333333   6.941 > 1
   24 Li  0.66666666666667  0.50000000000000  0.83333333333333   6.941 > 1
   25 Li  0.00000000000000  0.83333333333333  0.83333333333333   6.941 > 1
   26 Li  0.33333333333333  0.83333333333333  0.83333333333333   6.941 > 1
   27 Li  0.66666666666667  0.83333333333333  0.83333333333333   6.941 > 1
   28 Li  0.16666666666667  0.00000000000000  0.16666666666667   6.941 > 28
   29 Li  0.50000000000000  0.00000000000000  0.16666666666667   6.941 > 28
   30 Li  0.83333333333333  0.00000000000000  0.16666666666667   6.941 > 28
   31 Li  0.16666666666667  0.33333333333333  0.16666666666667   6.941 > 28
   32 Li  0.50000000000000  0.33333333333333  0.16666666666667   6.941 > 28
   33 Li  0.83333333333333  0.33333333333333  0.16666666666667   6.941 > 28
   34 Li  0.16666666666667  0.66666666666667  0.16666666666667   6.941 > 28
   35 Li  0.50000000000000  0.66666666666667  0.16666666666667   6.941 > 28
   36 Li  0.83333333333333  0.66666666666667  0.16666666666667   6.941 > 28
   37 Li  0.16666666666667  0.00000000000000  0.50000000000000   6.941 > 28
   38 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 28
   39 Li  0.83333333333333  0.00000000000000  0.50000000000000   6.941 > 28
   40 Li  0.16666666666667  0.33333333333333  0.50000000000000   6.941 > 28
   41 Li  0.50000000000000  0.33333333333333  0.50000000000000   6.941 > 28
   42 Li  0.83333333333333  0.33333333333333  0.50000000000000   6.941 > 28
   43 Li  0.16666666666667  0.66666666666667  0.50000000000000   6.941 > 28
   44 Li  0.50000000000000  0.66666666666667  0.50000000000000   6.941 > 28
   45 Li  0.83333333333333  0.66666666666667  0.50000000000000   6.941 > 28
   46 Li  0.16666666666667  0.00000000000000  0.83333333333333   6.941 > 28
   47 Li  0.50000000000000  0.00000000000000  0.83333333333333   6.941 > 28
   48 Li  0.83333333333333  0.00000000000000  0.83333333333333   6.941 > 28
   49 Li  0.16666666666667  0.33333333333333  0.83333333333333   6.941 > 28
   50 Li  0.50000000000000  0.33333333333333  0.83333333333333   6.941 > 28
   51 Li  0.83333333333333  0.33333333333333  0.83333333333333   6.941 > 28
   52 Li  0.16666666666667  0.66666666666667  0.83333333333333   6.941 > 28
   53 Li  0.50000000000000  0.66666666666667  0.83333333333333   6.941 > 28
   54 Li  0.83333333333333  0.66666666666667  0.83333333333333   6.941 > 28
   55 Li  0.16666666666667  0.16666666666667  0.00000000000000   6.941 > 55
   56 Li  0.50000000000000  0.16666666666667  0.00000000000000   6.941 > 55
   57 Li  0.83333333333333  0.16666666666667  0.00000000000000   6.941 > 55
   58 Li  0.16666666666667  0.50000000000000  0.00000000000000   6.941 > 55
   59 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 55
   60 Li  0.83333333333333  0.50000000000000  0.00000000000000   6.941 > 55
   61 Li  0.16666666666667  0.83333333333333  0.00000000000000   6.941 > 55
   62 Li  0.50000000000000  0.83333333333333  0.00000000000000   6.941 > 55
   63 Li  0.83333333333333  0.83333333333333  0.00000000000000   6.941 > 55
   64 Li  0.16666666666667  0.16666666666667  0.33333333333333   6.941 > 55
   65 Li  0.50000000000000  0.16666666666667  0.33333333333333   6.941 > 55
   66 Li  0.83333333333333  0.16666666666667  0.33333333333333   6.941 > 55
   67 Li  0.16666666666667  0.50000000000000  0.33333333333333   6.941 > 55
   68 Li  0.50000000000000  0.50000000000000  0.33333333333333   6.941 > 55
   69 Li  0.83333333333333  0.50000000000000  0.33333333333333   6.941 > 55
   70 Li  0.16666666666667  0.83333333333333  0.33333333333333   6.941 > 55
   71 Li  0.50000000000000  0.83333333333333  0.33333333333333   6.941 > 55
   72 Li  0.83333333333333  0.83333333333333  0.33333333333333   6.941 > 55
   73 Li  0.16666666666667  0.16666666666667  0.66666666666667   6.941 > 55
   74 Li  0.50000000000000  0.16666666666667  0.66666666666667   6.941 > 55
   75 Li  0.83333333333333  0.16666666666667  0.66666666666667   6.941 > 55
   76 Li  0.16666666666667  0.50000000000000  0.66666666666667   6.941 > 55
   77 Li  0.50000000000000  0.50000000000000  0.66666666666667   6.941 > 55
   78 Li  0.83333333333333  0.50000000000000  0.66666666666667   6.941 > 55
   79 Li  0.16666666666667  0.83333333333333  0.66666666666667   6.941 > 55
   80 Li  0.50000000000000  0.83333333333333  0.66666666666667   6.941 > 55
   81 Li  0.83333333333333  0.83333333333333  0.66666666666667   6.941 > 55
   82 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 82
   83 Cl  0.33333333333333  0.00000000000000  0.00000000000000  35.453 > 82
   84 Cl  0.66666666666667  0.00000000000000  0.00000000000000  35.453 > 82
   85 Cl  0.00000000000000  0.33333333333333  0.00000000000000  35.453 > 82
   86 Cl  0.33333333333333  0.33333333333333  0.00000000000000  35.453 > 82
   87 Cl  0.66666666666667  0.33333333333333  0.00000000000000  35.453 > 82
   88 Cl  0.00000000000000  0.66666666666667  0.00000000000000  35.453 > 82
   89 Cl  0.33333333333333  0.66666666666667  0.00000000000000  35.453 > 82
   90 Cl  0.66666666666667  0.66666666666667  0.00000000000000  35.453 > 82
   91 Cl  0.00000000000000  0.00000000000000  0.33333333333333  35.453 > 82
   92 Cl  0.33333333333333  0.00000000000000  0.33333333333333  35.453 > 82
   93 Cl  0.66666666666667  0.00000000000000  0.33333333333333  35.453 > 82
   94 Cl  0.00000000000000  0.33333333333333  0.33333333333333  35.453 > 82
   95 Cl  0.33333333333333  0.33333333333333  0.33333333333333  35.453 > 82
   96 Cl  0.66666666666667  0.33333333333333  0.33333333333333  35.453 > 82
   97 Cl  0.00000000000000  0.66666666666667  0.33333333333333  35.453 > 82
   98 Cl  0.33333333333333  0.66666666666667  0.33333333333333  35.453 > 82
   99 Cl  0.66666666666667  0.66666666666667  0.33333333333333  35.453 > 82
  100 Cl  0.00000000000000  0.00000000000000  0.66666666666667  35.453 > 82
  101 Cl  0.33333333333333  0.00000000000000  0.66666666666667  35.453 > 82
  102 Cl  0.66666666666667  0.00000000000000  0.66666666666667  35.453 > 82
  103 Cl  0.00000000000000  0.33333333333333  0.66666666666667  35.453 > 82
  104 Cl  0.33333333333333  0.33333333333333  0.66666666666667  35.453 > 82
  105 Cl  0.66666666666667  0.33333333333333  0.66666666666667  35.453 > 82
  106 Cl  0.00000000000000  0.66666666666667  0.66666666666667  35.453 > 82
  107 Cl  0.33333333333333  0.66666666666667  0.66666666666667  35.453 > 82
  108 Cl  0.66666666666667  0.66666666666667  0.66666666666667  35.453 > 82
  109 O   0.16666666666667  0.16666666666667  0.16666666666667  15.999 > 109
  110 O   0.50000000000000  0.16666666666667  0.16666666666667  15.999 > 109
  111 O   0.83333333333333  0.16666666666667  0.16666666666667  15.999 > 109
  112 O   0.16666666666667  0.50000000000000  0.16666666666667  15.999 > 109
  113 O   0.50000000000000  0.50000000000000  0.16666666666667  15.999 > 109
  114 O   0.83333333333333  0.50000000000000  0.16666666666667  15.999 > 109
  115 O   0.16666666666667  0.83333333333333  0.16666666666667  15.999 > 109
  116 O   0.50000000000000  0.83333333333333  0.16666666666667  15.999 > 109
  117 O   0.83333333333333  0.83333333333333  0.16666666666667  15.999 > 109
  118 O   0.16666666666667  0.16666666666667  0.50000000000000  15.999 > 109
  119 O   0.50000000000000  0.16666666666667  0.50000000000000  15.999 > 109
  120 O   0.83333333333333  0.16666666666667  0.50000000000000  15.999 > 109
  121 O   0.16666666666667  0.50000000000000  0.50000000000000  15.999 > 109
  122 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 109
  123 O   0.83333333333333  0.50000000000000  0.50000000000000  15.999 > 109
  124 O   0.16666666666667  0.83333333333333  0.50000000000000  15.999 > 109
  125 O   0.50000000000000  0.83333333333333  0.50000000000000  15.999 > 109
  126 O   0.83333333333333  0.83333333333333  0.50000000000000  15.999 > 109
  127 O   0.16666666666667  0.16666666666667  0.83333333333333  15.999 > 109
  128 O   0.50000000000000  0.16666666666667  0.83333333333333  15.999 > 109
  129 O   0.83333333333333  0.16666666666667  0.83333333333333  15.999 > 109
  130 O   0.16666666666667  0.50000000000000  0.83333333333333  15.999 > 109
  131 O   0.50000000000000  0.50000000000000  0.83333333333333  15.999 > 109
  132 O   0.83333333333333  0.50000000000000  0.83333333333333  15.999 > 109
  133 O   0.16666666666667  0.83333333333333  0.83333333333333  15.999 > 109
  134 O   0.50000000000000  0.83333333333333  0.83333333333333  15.999 > 109
  135 O   0.83333333333333  0.83333333333333  0.83333333333333  15.999 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.2441248    0.0000000    0.0000000
            0.0000000    3.2441248    0.0000000
            0.0000000    0.0000000    3.2441248
-------------------------- Born effective charges --------------------------
    1 Li    0.9906835    0.0000000    0.0000000
            0.0000000    1.1502198    0.0000000
            0.0000000    0.0000000    1.1502198
    2 Li    1.1502198    0.0000000    0.0000000
            0.0000000    0.9906835    0.0000000
            0.0000000    0.0000000    1.1502198
    3 Li    1.1502198    0.0000000    0.0000000
            0.0000000    1.1502198    0.0000000
            0.0000000    0.0000000    0.9906835
    4 Cl   -1.2699045    0.0000000    0.0000000
            0.0000000   -1.2699045    0.0000000
            0.0000000    0.0000000   -1.2699045
    5 O    -2.0212186    0.0000000    0.0000000
            0.0000000   -2.0212186    0.0000000
            0.0000000    0.0000000   -2.0212186
----------------------------------------------------------------------------
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 (yyy) -0.00000000 (yyy) -0.00000000 (yyy)
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-09 04:49:58]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-09 04:49: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.837396259999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    3.837396259999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    3.837396259999999
Atomic positions (fractional):
    1 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941
    2 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941
    3 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941
    4 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453
    5 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   11.512188779999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   11.512188779999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   11.512188779999997
Atomic positions (fractional):
    1 Li  0.00000000000000  0.16666666666667  0.16666666666667   6.941 > 1
    2 Li  0.33333333333333  0.16666666666667  0.16666666666667   6.941 > 1
    3 Li  0.66666666666667  0.16666666666667  0.16666666666667   6.941 > 1
    4 Li  0.00000000000000  0.50000000000000  0.16666666666667   6.941 > 1
    5 Li  0.33333333333333  0.50000000000000  0.16666666666667   6.941 > 1
    6 Li  0.66666666666667  0.50000000000000  0.16666666666667   6.941 > 1
    7 Li  0.00000000000000  0.83333333333333  0.16666666666667   6.941 > 1
    8 Li  0.33333333333333  0.83333333333333  0.16666666666667   6.941 > 1
    9 Li  0.66666666666667  0.83333333333333  0.16666666666667   6.941 > 1
   10 Li  0.00000000000000  0.16666666666667  0.50000000000000   6.941 > 1
   11 Li  0.33333333333333  0.16666666666667  0.50000000000000   6.941 > 1
   12 Li  0.66666666666667  0.16666666666667  0.50000000000000   6.941 > 1
   13 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1
   14 Li  0.33333333333333  0.50000000000000  0.50000000000000   6.941 > 1
   15 Li  0.66666666666667  0.50000000000000  0.50000000000000   6.941 > 1
   16 Li  0.00000000000000  0.83333333333333  0.50000000000000   6.941 > 1
   17 Li  0.33333333333333  0.83333333333333  0.50000000000000   6.941 > 1
   18 Li  0.66666666666667  0.83333333333333  0.50000000000000   6.941 > 1
   19 Li  0.00000000000000  0.16666666666667  0.83333333333333   6.941 > 1
   20 Li  0.33333333333333  0.16666666666667  0.83333333333333   6.941 > 1
   21 Li  0.66666666666667  0.16666666666667  0.83333333333333   6.941 > 1
   22 Li  0.00000000000000  0.50000000000000  0.83333333333333   6.941 > 1
   23 Li  0.33333333333333  0.50000000000000  0.83333333333333   6.941 > 1
   24 Li  0.66666666666667  0.50000000000000  0.83333333333333   6.941 > 1
   25 Li  0.00000000000000  0.83333333333333  0.83333333333333   6.941 > 1
   26 Li  0.33333333333333  0.83333333333333  0.83333333333333   6.941 > 1
   27 Li  0.66666666666667  0.83333333333333  0.83333333333333   6.941 > 1
   28 Li  0.16666666666667  0.00000000000000  0.16666666666667   6.941 > 28
   29 Li  0.50000000000000  0.00000000000000  0.16666666666667   6.941 > 28
   30 Li  0.83333333333333  0.00000000000000  0.16666666666667   6.941 > 28
   31 Li  0.16666666666667  0.33333333333333  0.16666666666667   6.941 > 28
   32 Li  0.50000000000000  0.33333333333333  0.16666666666667   6.941 > 28
   33 Li  0.83333333333333  0.33333333333333  0.16666666666667   6.941 > 28
   34 Li  0.16666666666667  0.66666666666667  0.16666666666667   6.941 > 28
   35 Li  0.50000000000000  0.66666666666667  0.16666666666667   6.941 > 28
   36 Li  0.83333333333333  0.66666666666667  0.16666666666667   6.941 > 28
   37 Li  0.16666666666667  0.00000000000000  0.50000000000000   6.941 > 28
   38 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 28
   39 Li  0.83333333333333  0.00000000000000  0.50000000000000   6.941 > 28
   40 Li  0.16666666666667  0.33333333333333  0.50000000000000   6.941 > 28
   41 Li  0.50000000000000  0.33333333333333  0.50000000000000   6.941 > 28
   42 Li  0.83333333333333  0.33333333333333  0.50000000000000   6.941 > 28
   43 Li  0.16666666666667  0.66666666666667  0.50000000000000   6.941 > 28
   44 Li  0.50000000000000  0.66666666666667  0.50000000000000   6.941 > 28
   45 Li  0.83333333333333  0.66666666666667  0.50000000000000   6.941 > 28
   46 Li  0.16666666666667  0.00000000000000  0.83333333333333   6.941 > 28
   47 Li  0.50000000000000  0.00000000000000  0.83333333333333   6.941 > 28
   48 Li  0.83333333333333  0.00000000000000  0.83333333333333   6.941 > 28
   49 Li  0.16666666666667  0.33333333333333  0.83333333333333   6.941 > 28
   50 Li  0.50000000000000  0.33333333333333  0.83333333333333   6.941 > 28
   51 Li  0.83333333333333  0.33333333333333  0.83333333333333   6.941 > 28
   52 Li  0.16666666666667  0.66666666666667  0.83333333333333   6.941 > 28
   53 Li  0.50000000000000  0.66666666666667  0.83333333333333   6.941 > 28
   54 Li  0.83333333333333  0.66666666666667  0.83333333333333   6.941 > 28
   55 Li  0.16666666666667  0.16666666666667  0.00000000000000   6.941 > 55
   56 Li  0.50000000000000  0.16666666666667  0.00000000000000   6.941 > 55
   57 Li  0.83333333333333  0.16666666666667  0.00000000000000   6.941 > 55
   58 Li  0.16666666666667  0.50000000000000  0.00000000000000   6.941 > 55
   59 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 55
   60 Li  0.83333333333333  0.50000000000000  0.00000000000000   6.941 > 55
   61 Li  0.16666666666667  0.83333333333333  0.00000000000000   6.941 > 55
   62 Li  0.50000000000000  0.83333333333333  0.00000000000000   6.941 > 55
   63 Li  0.83333333333333  0.83333333333333  0.00000000000000   6.941 > 55
   64 Li  0.16666666666667  0.16666666666667  0.33333333333333   6.941 > 55
   65 Li  0.50000000000000  0.16666666666667  0.33333333333333   6.941 > 55
   66 Li  0.83333333333333  0.16666666666667  0.33333333333333   6.941 > 55
   67 Li  0.16666666666667  0.50000000000000  0.33333333333333   6.941 > 55
   68 Li  0.50000000000000  0.50000000000000  0.33333333333333   6.941 > 55
   69 Li  0.83333333333333  0.50000000000000  0.33333333333333   6.941 > 55
   70 Li  0.16666666666667  0.83333333333333  0.33333333333333   6.941 > 55
   71 Li  0.50000000000000  0.83333333333333  0.33333333333333   6.941 > 55
   72 Li  0.83333333333333  0.83333333333333  0.33333333333333   6.941 > 55
   73 Li  0.16666666666667  0.16666666666667  0.66666666666667   6.941 > 55
   74 Li  0.50000000000000  0.16666666666667  0.66666666666667   6.941 > 55
   75 Li  0.83333333333333  0.16666666666667  0.66666666666667   6.941 > 55
   76 Li  0.16666666666667  0.50000000000000  0.66666666666667   6.941 > 55
   77 Li  0.50000000000000  0.50000000000000  0.66666666666667   6.941 > 55
   78 Li  0.83333333333333  0.50000000000000  0.66666666666667   6.941 > 55
   79 Li  0.16666666666667  0.83333333333333  0.66666666666667   6.941 > 55
   80 Li  0.50000000000000  0.83333333333333  0.66666666666667   6.941 > 55
   81 Li  0.83333333333333  0.83333333333333  0.66666666666667   6.941 > 55
   82 Cl  0.00000000000000  0.00000000000000  0.00000000000000  35.453 > 82
   83 Cl  0.33333333333333  0.00000000000000  0.00000000000000  35.453 > 82
   84 Cl  0.66666666666667  0.00000000000000  0.00000000000000  35.453 > 82
   85 Cl  0.00000000000000  0.33333333333333  0.00000000000000  35.453 > 82
   86 Cl  0.33333333333333  0.33333333333333  0.00000000000000  35.453 > 82
   87 Cl  0.66666666666667  0.33333333333333  0.00000000000000  35.453 > 82
   88 Cl  0.00000000000000  0.66666666666667  0.00000000000000  35.453 > 82
   89 Cl  0.33333333333333  0.66666666666667  0.00000000000000  35.453 > 82
   90 Cl  0.66666666666667  0.66666666666667  0.00000000000000  35.453 > 82
   91 Cl  0.00000000000000  0.00000000000000  0.33333333333333  35.453 > 82
   92 Cl  0.33333333333333  0.00000000000000  0.33333333333333  35.453 > 82
   93 Cl  0.66666666666667  0.00000000000000  0.33333333333333  35.453 > 82
   94 Cl  0.00000000000000  0.33333333333333  0.33333333333333  35.453 > 82
   95 Cl  0.33333333333333  0.33333333333333  0.33333333333333  35.453 > 82
   96 Cl  0.66666666666667  0.33333333333333  0.33333333333333  35.453 > 82
   97 Cl  0.00000000000000  0.66666666666667  0.33333333333333  35.453 > 82
   98 Cl  0.33333333333333  0.66666666666667  0.33333333333333  35.453 > 82
   99 Cl  0.66666666666667  0.66666666666667  0.33333333333333  35.453 > 82
  100 Cl  0.00000000000000  0.00000000000000  0.66666666666667  35.453 > 82
  101 Cl  0.33333333333333  0.00000000000000  0.66666666666667  35.453 > 82
  102 Cl  0.66666666666667  0.00000000000000  0.66666666666667  35.453 > 82
  103 Cl  0.00000000000000  0.33333333333333  0.66666666666667  35.453 > 82
  104 Cl  0.33333333333333  0.33333333333333  0.66666666666667  35.453 > 82
  105 Cl  0.66666666666667  0.33333333333333  0.66666666666667  35.453 > 82
  106 Cl  0.00000000000000  0.66666666666667  0.66666666666667  35.453 > 82
  107 Cl  0.33333333333333  0.66666666666667  0.66666666666667  35.453 > 82
  108 Cl  0.66666666666667  0.66666666666667  0.66666666666667  35.453 > 82
  109 O   0.16666666666667  0.16666666666667  0.16666666666667  15.999 > 109
  110 O   0.50000000000000  0.16666666666667  0.16666666666667  15.999 > 109
  111 O   0.83333333333333  0.16666666666667  0.16666666666667  15.999 > 109
  112 O   0.16666666666667  0.50000000000000  0.16666666666667  15.999 > 109
  113 O   0.50000000000000  0.50000000000000  0.16666666666667  15.999 > 109
  114 O   0.83333333333333  0.50000000000000  0.16666666666667  15.999 > 109
  115 O   0.16666666666667  0.83333333333333  0.16666666666667  15.999 > 109
  116 O   0.50000000000000  0.83333333333333  0.16666666666667  15.999 > 109
  117 O   0.83333333333333  0.83333333333333  0.16666666666667  15.999 > 109
  118 O   0.16666666666667  0.16666666666667  0.50000000000000  15.999 > 109
  119 O   0.50000000000000  0.16666666666667  0.50000000000000  15.999 > 109
  120 O   0.83333333333333  0.16666666666667  0.50000000000000  15.999 > 109
  121 O   0.16666666666667  0.50000000000000  0.50000000000000  15.999 > 109
  122 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 109
  123 O   0.83333333333333  0.50000000000000  0.50000000000000  15.999 > 109
  124 O   0.16666666666667  0.83333333333333  0.50000000000000  15.999 > 109
  125 O   0.50000000000000  0.83333333333333  0.50000000000000  15.999 > 109
  126 O   0.83333333333333  0.83333333333333  0.50000000000000  15.999 > 109
  127 O   0.16666666666667  0.16666666666667  0.83333333333333  15.999 > 109
  128 O   0.50000000000000  0.16666666666667  0.83333333333333  15.999 > 109
  129 O   0.83333333333333  0.16666666666667  0.83333333333333  15.999 > 109
  130 O   0.16666666666667  0.50000000000000  0.83333333333333  15.999 > 109
  131 O   0.50000000000000  0.50000000000000  0.83333333333333  15.999 > 109
  132 O   0.83333333333333  0.50000000000000  0.83333333333333  15.999 > 109
  133 O   0.16666666666667  0.83333333333333  0.83333333333333  15.999 > 109
  134 O   0.50000000000000  0.83333333333333  0.83333333333333  15.999 > 109
  135 O   0.83333333333333  0.83333333333333  0.83333333333333  15.999 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.2441248    0.0000000    0.0000000
            0.0000000    3.2441248    0.0000000
            0.0000000    0.0000000    3.2441248
-------------------------- Born effective charges --------------------------
    1 Li    0.9906835    0.0000000    0.0000000
            0.0000000    1.1502198    0.0000000
            0.0000000    0.0000000    1.1502198
    2 Li    1.1502198    0.0000000    0.0000000
            0.0000000    0.9906835    0.0000000
            0.0000000    0.0000000    1.1502198
    3 Li    1.1502198    0.0000000    0.0000000
            0.0000000    1.1502198    0.0000000
            0.0000000    0.0000000    0.9906835
    4 Cl   -1.2699045    0.0000000    0.0000000
            0.0000000   -1.2699045    0.0000000
            0.0000000    0.0000000   -1.2699045
    5 O    -2.0212186    0.0000000    0.0000000
            0.0000000   -2.0212186    0.0000000
            0.0000000    0.0000000   -2.0212186
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000000 (yyy) -0.00000000 (yyy) -0.00000000 (yyy)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 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.08, Number of G-points: 305, Lambda: 0.20
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 5.83e-04 3.36e-05 
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.541   (   0.000    0.000    0.000)    0.000
   4.541   (   0.000    0.000    0.000)    0.000
   4.541   (   0.000    0.000    0.000)    0.000
   6.594   (   0.000    0.000    0.000)    0.000
   6.594   (   0.000    0.000    0.000)    0.000
   6.594   (   0.000    0.000    0.000)    0.000
  11.103   (   0.000    0.000    0.000)    0.000
  11.103   (   0.000    0.000    0.000)    0.000
  11.103   (   0.000    0.000    0.000)    0.000
  17.330   (   0.000    0.000    0.000)    0.000
  17.330   (   0.000    0.000    0.000)    0.000
  17.330   (   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 5.83e-04 3.36e-05 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.941   (  45.939    0.000    0.000)   45.939
   0.941   (  45.939    0.000    0.000)   45.939
   1.403   (  67.739    0.000    0.000)   67.739
   4.551   (   0.987    0.000    0.000)    0.987
   4.551   (   0.987    0.000    0.000)    0.987
   6.633   (   3.793    0.000    0.000)    3.793
   6.633   (   3.793    0.000    0.000)    3.793
   6.653   (   5.698    0.000    0.000)    5.698
   6.817   (   9.849    0.000    0.000)    9.849
  11.081   (  -2.215    0.000    0.000)    2.215
  11.081   (  -2.215    0.000    0.000)    2.215
  14.252   (   5.385    0.000    0.000)    5.385
  17.325   (  -0.538    0.000    0.000)    0.538
  17.325   (  -0.538    0.000    0.000)    0.538
  19.823   (  -8.124    0.000    0.000)    8.124
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.822   (  41.518    0.000    0.000)   41.518
   1.822   (  41.518    0.000    0.000)   41.518
   2.674   (  58.152    0.000    0.000)   58.152
   4.579   (   1.777    0.000    0.000)    1.777
   4.579   (   1.777    0.000    0.000)    1.777
   6.743   (   7.094    0.000    0.000)    7.094
   6.743   (   7.094    0.000    0.000)    7.094
   6.812   (   9.840    0.000    0.000)    9.840
   7.105   (  18.515    0.000    0.000)   18.515
  11.016   (  -4.185    0.000    0.000)    4.185
  11.016   (  -4.185    0.000    0.000)    4.185
  14.419   (  11.450    0.000    0.000)   11.450
  17.310   (  -0.944    0.000    0.000)    0.944
  17.310   (  -0.944    0.000    0.000)    0.944
  19.580   ( -16.057    0.000    0.000)   16.057
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.588   (  34.482    0.000    0.000)   34.482
   2.588   (  34.482    0.000    0.000)   34.482
   3.705   (  44.338    0.000    0.000)   44.338
   4.620   (   2.218    0.000    0.000)    2.218
   4.620   (   2.218    0.000    0.000)    2.218
   6.910   (   9.338    0.000    0.000)    9.338
   6.910   (   9.338    0.000    0.000)    9.338
   7.030   (  11.477    0.000    0.000)   11.477
   7.537   (  23.823    0.000    0.000)   23.823
  10.917   (  -5.583    0.000    0.000)    5.583
  10.917   (  -5.583    0.000    0.000)    5.583
  14.718   (  18.577    0.000    0.000)   18.577
  17.288   (  -1.124    0.000    0.000)    1.124
  17.288   (  -1.124    0.000    0.000)    1.124
  19.182   ( -23.528    0.000    0.000)   23.528
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.191   (  25.451    0.000    0.000)   25.451
   3.191   (  25.451    0.000    0.000)   25.451
   4.447   (  29.817    0.000    0.000)   29.817
   4.665   (   2.221    0.000    0.000)    2.221
   4.665   (   2.221    0.000    0.000)    2.221
   7.105   (   9.792    0.000    0.000)    9.792
   7.105   (   9.792    0.000    0.000)    9.792
   7.254   (  10.452    0.000    0.000)   10.452
   8.019   (  23.123    0.000    0.000)   23.123
  10.799   (  -5.943    0.000    0.000)    5.943
  10.799   (  -5.943    0.000    0.000)    5.943
  15.168   (  26.329    0.000    0.000)   26.329
  17.266   (  -1.043    0.000    0.000)    1.043
  17.266   (  -1.043    0.000    0.000)    1.043
  18.642   ( -30.109    0.000    0.000)   30.109
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.601   (  15.353    0.000    0.000)   15.353
   3.601   (  15.353    0.000    0.000)   15.353
   4.705   (   1.719    0.000    0.000)    1.719
   4.705   (   1.719    0.000    0.000)    1.719
   4.912   (  16.893    0.000    0.000)   16.893
   7.285   (   7.694    0.000    0.000)    7.694
   7.285   (   7.694    0.000    0.000)    7.694
   7.435   (   7.219    0.000    0.000)    7.219
   8.422   (  16.238    0.000    0.000)   16.238
  10.689   (  -4.744    0.000    0.000)    4.744
  10.689   (  -4.744    0.000    0.000)    4.744
  15.767   (  33.144    0.000    0.000)   33.144
  17.248   (  -0.731    0.000    0.000)    0.731
  17.248   (  -0.731    0.000    0.000)    0.731
  17.985   ( -35.126    0.000    0.000)   35.126
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.806   (   5.081    0.000    0.000)    5.081
   3.806   (   5.081    0.000    0.000)    5.081
   4.730   (   0.670    0.000    0.000)    0.670
   4.730   (   0.670    0.000    0.000)    0.670
   5.134   (   5.462    0.000    0.000)    5.462
   7.396   (   2.961    0.000    0.000)    2.961
   7.396   (   2.961    0.000    0.000)    2.961
   7.534   (   2.568    0.000    0.000)    2.568
   8.646   (   5.747    0.000    0.000)    5.747
  10.621   (  -1.847    0.000    0.000)    1.847
  10.621   (  -1.847    0.000    0.000)    1.847
  16.475   (  36.548    0.000    0.000)   36.548
  17.238   (  -0.262    0.000    0.000)    0.262
  17.238   (  -0.262    0.000    0.000)    0.262
  17.253   ( -37.094    0.000    0.000)   37.094
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.124   (  26.973   26.973    0.000)   38.145
   1.322   (  31.834   31.834    0.000)   45.019
   2.061   (  47.640   47.640    0.000)   67.373
   4.563   (   1.206    1.206    0.000)    1.706
   4.594   (   2.433    2.433    0.000)    3.440
   6.532   (  -3.148   -3.148    0.000)    4.451
   6.670   (   3.691    3.691    0.000)    5.220
   6.728   (   2.079    2.079    0.000)    2.940
   6.890   (  13.282   13.282    0.000)   18.783
  11.057   (  -2.392   -2.392    0.000)    3.383
  11.133   (   1.629    1.629    0.000)    2.303
  14.351   (   7.274    7.274    0.000)   10.287
  17.186   (  -7.088   -7.088    0.000)   10.025
  17.319   (  -0.533   -0.533    0.000)    0.754
  19.823   (  -4.083   -4.083    0.000)    5.774
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.829   (  38.541    5.516    0.000)   38.934
   2.028   (  36.132   19.016    0.000)   40.830
   3.032   (  46.508   29.377    0.000)   55.010
   4.598   (   2.307    1.962    0.000)    3.029
   4.634   (   1.933    4.142    0.000)    4.571
   6.491   (  -1.180  -10.521    0.000)   10.587
   6.777   (   6.912    3.405    0.000)    7.705
   6.910   (  14.610   -7.248    0.000)   16.309
   7.229   (  20.367   19.974    0.000)   28.527
  10.987   (  -4.528   -2.913    0.000)    5.385
  11.121   (  -2.310    7.802    0.000)    8.137
  14.524   (  10.469    9.168    0.000)   13.916
  17.093   (  -2.981  -17.009    0.000)   17.269
  17.304   (  -0.936   -0.521    0.000)    1.071
  19.633   ( -14.303    2.145    0.000)   14.463
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.575   (  34.754   -0.188    0.000)   34.755
   2.714   (  31.398   11.827    0.000)   33.552
   3.862   (  35.565   13.822    0.000)   38.156
   4.654   (   3.212    3.522    0.000)    4.767
   4.683   (   3.209    5.021    0.000)    5.959
   6.479   (  -0.232  -20.916    0.000)   20.917
   6.940   (   9.130    3.002    0.000)    9.610
   7.246   (  17.538   -4.119    0.000)   18.015
   7.704   (  26.669   23.896    0.000)   35.808
  10.879   (  -6.060   -3.739    0.000)    7.120
  11.052   (  -4.430   11.474    0.000)   12.299
  14.787   (  16.145    6.549    0.000)   17.423
  17.047   (  -1.937  -20.986    0.000)   21.075
  17.283   (  -1.115   -0.504    0.000)    1.223
  19.267   ( -21.881    5.923    0.000)   22.668
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.191   (  26.167   -0.404    0.000)   26.170
   3.263   (  23.015    6.651    0.000)   23.957
   4.447   (  22.782    0.047    0.000)   22.782
   4.725   (   3.737    6.118    0.000)    7.169
   4.769   (   5.514    8.852    0.000)   10.429
   6.473   (  -0.507  -33.054    0.000)   33.058
   7.132   (   9.641    2.623    0.000)    9.991
   7.577   (  14.900    1.686    0.000)   14.995
   8.264   (  27.967   27.723    0.000)   39.379
  10.751   (  -6.482   -4.771    0.000)    8.048
  10.948   (  -5.876   13.321    0.000)   14.559
  15.178   (  22.853    1.036    0.000)   22.876
  17.010   (  -1.858  -22.477    0.000)   22.554
  17.261   (  -1.035   -0.487    0.000)    1.144
  18.769   ( -27.318    9.799    0.000)   29.023
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.608   (  15.279    0.222    0.000)   15.281
   3.628   (  13.383    2.369    0.000)   13.591
   4.774   (  10.001   -1.988    0.000)   10.197
   4.798   (   3.323    9.304    0.000)    9.879
   4.913   (   8.655    7.984    0.000)   11.775
   6.456   (  -1.147  -45.541    0.000)   45.555
   7.310   (   7.645    2.417    0.000)    8.018
   7.828   (   9.814    4.934    0.000)   10.984
   8.775   (  21.770   33.869    0.000)   40.263
  10.631   (  -5.204   -5.802    0.000)    7.793
  10.824   (  -6.147   12.625    0.000)   14.042
  15.692   (  27.740   -7.130    0.000)   28.642
  16.964   (  -3.087  -22.986    0.000)   23.193
  17.243   (  -0.726   -0.474    0.000)    0.867
  18.199   ( -28.389   15.904    0.000)   32.540
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.803   (   4.272   -0.275    0.000)    4.280
   3.807   (   4.787    0.002    0.000)    4.787
   4.847   (   1.385   11.567    0.000)   11.650
   4.895   (   2.909   11.344    0.000)   11.711
   5.054   (   4.017   -3.417    0.000)    5.274
   6.435   (  -0.715  -54.632    0.000)   54.637
   7.420   (   2.963    2.383    0.000)    3.802
   7.962   (   3.413    6.039    0.000)    6.937
   9.090   (   8.588   39.984    0.000)   40.896
  10.555   (  -2.035   -6.484    0.000)    6.796
  10.724   (  -3.052    9.961    0.000)   10.418
  16.225   (  21.646  -21.363    0.000)   30.413
  16.859   (  -7.547  -18.676    0.000)   20.143
  17.233   (  -0.260   -0.467    0.000)    0.534
  17.723   ( -15.262   24.426    0.000)   28.802
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.114   (  21.744   21.744    0.000)   30.751
   2.501   (  26.213   26.213    0.000)   37.071
   3.657   (  30.136   30.136    0.000)   42.618
   4.659   (   4.152    4.152    0.000)    5.872
   4.693   (   1.481    1.481    0.000)    2.095
   6.358   (  -4.738   -4.738    0.000)    6.700
   6.839   (   3.963    3.963    0.000)    5.604
   6.877   (   6.381    6.381    0.000)    9.025
   7.714   (  27.267   27.267    0.000)   38.562
  10.901   (  -5.540   -5.540    0.000)    7.834
  11.248   (   4.189    4.189    0.000)    5.923
  14.733   (  10.926   10.926    0.000)   15.451
  16.776   ( -13.102  -13.102    0.000)   18.529
  17.290   (  -0.915   -0.915    0.000)    1.294
  19.580   (  -7.864   -7.864    0.000)   11.121
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.626   (  27.665    6.036    0.000)   28.316
   3.015   (  23.838   16.868    0.000)   29.202
   4.169   (  19.849   14.931    0.000)   24.838
   4.726   (   3.110   -3.053    0.000)    4.358
   4.767   (   6.587    7.885    0.000)   10.274
   6.230   (  -8.137   -2.350    0.000)    8.470
   7.028   (   8.493    5.591    0.000)   10.168
   7.084   (  18.022   -7.548    0.000)   19.539
   8.306   (  31.357   33.225    0.000)   45.685
  10.769   (  -7.477   -7.154    0.000)   10.348
  11.290   (   0.185   10.865    0.000)   10.867
  14.956   (  11.613    9.561    0.000)   15.043
  16.578   (  -6.954  -23.430    0.000)   24.441
  17.269   (  -1.090   -0.886    0.000)    1.405
  19.329   ( -16.618   -0.811    0.000)   16.638
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.167   (  25.072   -1.771    0.000)   25.135
   3.425   (  16.630    8.576    0.000)   18.711
   4.444   (   8.824   -0.803    0.000)    8.861
   4.853   (   9.090   -6.429    0.000)   11.134
   4.919   (   8.269   13.363    0.000)   15.714
   6.040   ( -10.778   -4.312    0.000)   11.609
   7.207   (   9.145    4.782    0.000)   10.320
   7.463   (  18.073   -8.849    0.000)   20.123
   8.947   (  31.915   36.913    0.000)   48.797
  10.610   (  -8.104   -9.206    0.000)   12.265
  11.257   (  -3.550   15.935    0.000)   16.325
  15.210   (  13.922    2.150    0.000)   14.087
  16.482   (  -3.061  -28.310    0.000)   28.475
  17.248   (  -1.014   -0.857    0.000)    1.328
  18.939   ( -21.495    6.003    0.000)   22.318
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.589   (  16.174   -2.494    0.000)   16.366
   3.679   (   8.917    2.348    0.000)    9.221
   4.598   (   7.582  -12.436    0.000)   14.565
   5.044   (   9.428  -12.229    0.000)   15.442
   5.079   (   7.133   18.575    0.000)   19.897
   5.800   ( -13.009   -5.275    0.000)   14.037
   7.379   (   7.469    4.290    0.000)    8.613
   7.768   (  11.808   -7.532    0.000)   14.006
   9.548   (  26.988   39.883    0.000)   48.156
  10.458   (  -6.616  -11.314    0.000)   13.106
  11.146   (  -7.304   18.311    0.000)   19.714
  15.506   (  14.819  -10.781    0.000)   18.326
  16.443   (  -1.021  -27.839    0.000)   27.858
  17.230   (  -0.711   -0.835    0.000)    1.097
  18.505   ( -20.609   13.557    0.000)   24.668
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.793   (   2.716   -0.757    0.000)    2.820
   3.800   (   5.016   -0.881    0.000)    5.093
   4.733   (   4.482  -16.245    0.000)   16.852
   5.183   (   2.796   21.456    0.000)   21.638
   5.237   (   9.940  -30.152    0.000)   31.748
   5.533   ( -12.561    5.629    0.000)   13.764
   7.488   (   2.966    4.188    0.000)    5.132
   7.927   (   4.002   -6.619    0.000)    7.734
   9.970   (  12.834   44.240    0.000)   46.064
  10.362   (  -2.623  -12.756    0.000)   13.022
  10.996   (  -5.856   16.537    0.000)   17.543
  15.748   (   7.429  -24.526    0.000)   25.626
  16.433   (  -0.134  -24.088    0.000)   24.088
  17.220   (  -0.255   -0.823    0.000)    0.862
  18.182   (  -9.542   20.733    0.000)   22.823
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.824   (  13.423   13.423    0.000)   18.983
   3.339   (  14.548   14.548    0.000)   20.574
   4.414   (   8.188    8.188    0.000)   11.579
   4.576   (  -8.902   -8.902    0.000)   12.590
   4.974   (  12.634   12.634    0.000)   17.867
   6.285   (   2.730    2.730    0.000)    3.860
   7.049   (   5.986    5.986    0.000)    8.466
   7.160   (   7.442    7.442    0.000)   10.524
   8.996   (  35.061   35.061    0.000)   49.584
  10.598   (  -9.778   -9.778    0.000)   13.829
  11.466   (   6.568    6.568    0.000)    9.288
  15.141   (   8.231    8.231    0.000)   11.641
  16.164   ( -17.058  -17.058    0.000)   24.123
  17.249   (  -1.058   -1.058    0.000)    1.496
  19.227   (  -8.941   -8.941    0.000)   12.644
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.144   (  17.787    0.170    0.000)   17.788
   3.578   (   9.201    6.360    0.000)   11.185
   4.355   (  -7.694  -11.956    0.000)   14.218
   4.571   (   2.802  -12.512    0.000)   12.822
   5.255   (  14.657   19.689    0.000)   24.545
   6.250   (  -5.110   18.186    0.000)   18.890
   7.296   (  15.665   -6.722    0.000)   17.046
   7.319   (   8.200    6.183    0.000)   10.270
   9.690   (  33.622   36.537    0.000)   49.653
  10.388   ( -10.831  -12.811    0.000)   16.776
  11.550   (   1.673   12.922    0.000)   13.030
  15.260   (   3.211    2.657    0.000)    4.168
  15.922   (  -6.680  -26.778    0.000)   27.598
  17.228   (  -0.985   -1.025    0.000)    1.422
  18.995   ( -13.393   -0.250    0.000)   13.396
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.502   (  16.859   -6.251    0.000)   17.981
   3.715   (   4.662    1.189    0.000)    4.811
   4.244   (  -2.918  -25.018    0.000)   25.187
   4.565   (  -3.420  -19.909    0.000)   20.201
   5.525   (  11.482   25.354    0.000)   27.833
   6.136   (  -5.398   25.244    0.000)   25.814
   7.476   (   7.063    5.250    0.000)    8.800
   7.583   (  11.715  -10.273    0.000)   15.581
  10.182   (  -9.120  -16.121    0.000)   18.522
  10.321   (  28.760   36.764    0.000)   46.677
  11.521   (  -4.817   18.453    0.000)   19.071
  15.280   (  -0.327  -11.440    0.000)   11.444
  15.878   (   1.023  -27.858    0.000)   27.877
  17.211   (  -0.692   -1.000    0.000)    1.216
  18.727   ( -12.408    8.526    0.000)   15.056
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.761   (   7.328   -3.650    0.000)    8.186
   3.774   (   1.402   -1.116    0.000)    1.792
   4.252   (   3.527  -39.480    0.000)   39.637
   4.450   (  -6.792  -22.103    0.000)   23.123
   5.688   (   4.276   28.439    0.000)   28.759
   6.057   (  -2.121   28.723    0.000)   28.801
   7.582   (   2.956    5.001    0.000)    5.809
   7.743   (   4.057  -11.224    0.000)   11.935
  10.047   (  -3.717  -18.557    0.000)   18.926
  10.803   (  17.114   38.211    0.000)   41.868
  11.368   (  -8.682   20.295    0.000)   22.074
  15.273   (  -0.169  -22.550    0.000)   22.551
  15.911   (   1.304  -27.226    0.000)   27.258
  17.202   (  -0.248   -0.987    0.000)    1.017
  18.540   (  -5.213   14.922    0.000)   15.807
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.193   (   4.837    4.837    0.000)    6.841
   3.677   (   3.603    3.603    0.000)    5.095
   4.015   ( -16.640  -16.640    0.000)   23.533
   4.428   (  -5.787   -5.787    0.000)    8.184
   5.676   (  21.296   21.296    0.000)   30.117
   6.553   (   8.967    8.967    0.000)   12.681
   7.245   (   2.866    2.866    0.000)    4.053
   7.450   (   6.761    6.761    0.000)    9.561
  10.108   ( -14.679  -14.679    0.000)   20.759
  10.395   (  33.444   33.444    0.000)   47.296
  11.766   (   8.313    8.313    0.000)   11.757
  15.297   (  -1.246   -1.246    0.000)    1.762
  15.443   ( -18.469  -18.469    0.000)   26.119
  17.208   (  -0.957   -0.957    0.000)    1.353
  18.946   (  -4.180   -4.180    0.000)    5.911
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.321   (   6.892  -12.550    0.000)   14.318
   3.729   (   1.727    0.325    0.000)    1.757
   3.768   (  -7.066  -16.814    0.000)   18.238
   4.286   (  -8.468  -11.128    0.000)   13.984
   6.066   (  16.546   27.703    0.000)   32.268
   6.608   (  -1.254   20.049    0.000)   20.088
   7.378   (   7.948   -9.447    0.000)   12.346
   7.583   (   6.210    5.256    0.000)    8.136
   9.822   ( -13.065  -19.413    0.000)   23.400
  11.014   (  27.873   32.225    0.000)   42.606
  11.867   (   1.259   15.598    0.000)   15.649
  15.040   ( -14.178  -13.349    0.000)   19.473
  15.365   (   0.880  -21.816    0.000)   21.833
  17.191   (  -0.673   -0.935    0.000)    1.152
  18.852   (  -4.606    4.183    0.000)    6.222
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.407   (   1.247  -32.261    0.000)   32.285
   3.751   (   0.506   -1.150    0.000)    1.256
   3.752   (   3.756   -5.976    0.000)    7.059
   4.111   (  -7.367  -14.232    0.000)   16.026
   6.301   (   6.245   32.214    0.000)   32.814
   6.574   (  -1.224   21.948    0.000)   21.982
   7.497   (   3.197  -12.724    0.000)   13.120
   7.680   (   2.905    4.667    0.000)    5.497
   9.624   (  -5.646  -23.491    0.000)   24.160
  11.491   (  18.867   30.404    0.000)   35.782
  11.797   (  -8.267   21.925    0.000)   23.432
  14.844   (  -5.007  -20.345    0.000)   20.952
  15.377   (   0.279  -25.138    0.000)   25.139
  17.182   (  -0.242   -0.923    0.000)    0.955
  18.782   (  -1.955    9.319    0.000)    9.522
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.040   ( -12.888  -12.888    0.000)   18.227
   3.610   (  -0.945   -0.945    0.000)    1.337
   3.731   (  -0.027   -0.027    0.000)    0.038
   4.072   ( -10.533  -10.533    0.000)   14.895
   6.595   (  23.866   23.866    0.000)   33.751
   6.888   (   6.844    6.844    0.000)    9.678
   7.245   (  -2.573   -2.573    0.000)    3.638
   7.683   (   4.652    4.652    0.000)    6.579
   9.427   ( -19.095  -19.095    0.000)   27.004
  11.603   (  26.036   26.036    0.000)   36.820
  12.124   (   9.321    9.321    0.000)   13.182
  14.723   ( -16.820  -16.820    0.000)   23.787
  15.057   (  -9.557   -9.557    0.000)   13.516
  17.175   (  -0.659   -0.659    0.000)    0.932
  18.906   (   1.500    1.500    0.000)    2.122
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.862   (  -4.522  -20.701    0.000)   21.189
   3.669   (   5.510   -2.958    0.000)    6.254
   3.730   (  -0.046   -0.850    0.000)    0.851
   3.870   (  -8.906   -9.902    0.000)   13.318
   6.919   (  -0.840   12.390    0.000)   12.418
   6.956   (  10.205   32.335    0.000)   33.907
   7.259   (   1.669  -10.247    0.000)   10.382
   7.762   (   2.704    3.325    0.000)    4.286
   9.118   (  -9.553  -26.496    0.000)   28.165
  12.029   (  15.649   23.325    0.000)   28.088
  12.207   (  -1.798   17.728    0.000)   17.819
  14.460   (  -7.577  -17.643    0.000)   19.201
  14.936   (  -2.909  -18.111    0.000)   18.343
  17.166   (  -0.237   -0.651    0.000)    0.693
  18.922   (   0.305    4.862    0.000)    4.871
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 343
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.588   (  -6.739   -6.739    0.000)    9.531
   3.645   (   0.916    0.916    0.000)    1.296
   3.713   (  -5.638   -5.638    0.000)    7.973
   3.718   (  -0.312   -0.312    0.000)    0.441
   7.071   (   2.257    2.257    0.000)    3.192
   7.119   (  -2.392   -2.392    0.000)    3.383
   7.535   (  21.579   21.579    0.000)   30.518
   7.811   (   1.650    1.650    0.000)    2.334
   8.615   ( -19.985  -19.985    0.000)   28.263
  12.395   (  11.719   11.719    0.000)   16.573
  12.456   (   5.743    5.743    0.000)    8.122
  14.180   (  -8.517   -8.517    0.000)   12.045
  14.680   (  -6.826   -6.826    0.000)    9.653
  17.157   (  -0.234   -0.234    0.000)    0.331
  18.984   (   1.481    1.481    0.000)    2.094
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.435   (  22.142   22.142   22.142)   38.351
   1.435   (  22.142   22.142   22.142)   38.351
   2.541   (  38.020   38.020   38.020)   65.852
   4.594   (   1.400    1.400    1.400)    2.426
   4.594   (   1.400    1.400    1.400)    2.426
   6.413   (  -6.094   -6.094   -6.094)   10.555
   6.733   (   1.264    1.264    1.264)    2.189
   6.888   (   9.638    9.638    9.638)   16.694
   6.888   (   9.638    9.638    9.638)   16.694
  11.111   (   0.367    0.367    0.367)    0.636
  11.111   (   0.367    0.367    0.367)    0.636
  14.442   (   7.375    7.375    7.375)   12.773
  17.181   (  -4.891   -4.891   -4.891)    8.471
  17.181   (  -4.891   -4.891   -4.891)    8.471
  19.823   (  -2.691   -2.691   -2.691)    4.661
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.977   (  33.200   10.340   10.340)   36.278
   2.067   (  33.577   11.067   11.067)   37.045
   3.336   (  38.167   25.244   25.244)   52.261
   4.602   (   0.734    1.105    1.105)    1.726
   4.664   (   4.487    1.795    1.795)    5.155
   6.324   (  -2.970  -11.161  -11.161)   16.062
   6.773   (   2.043    2.052    2.052)    3.549
   7.123   (  20.016    6.800    6.800)   22.206
   7.245   (  19.527   11.515   11.515)   25.427
  11.046   (  -4.261    1.640    1.640)    4.851
  11.129   (  -0.994    3.019    3.019)    4.383
  14.611   (   9.646    7.780    7.780)   14.632
  17.037   (  -5.125   -9.001   -9.001)   13.722
  17.167   (  -0.907   -7.042   -7.042)   10.000
  19.668   ( -12.576    1.358    1.358)   12.722
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.652   (  32.372    3.395    3.395)   32.727
   2.722   (  30.870    6.253    6.253)   32.111
   4.000   (  27.311   11.557   11.557)   31.828
   4.627   (   1.870    0.221    0.221)    1.896
   4.786   (   7.806    6.083    6.083)   11.616
   6.284   (  -1.250  -16.404  -16.404)   23.232
   6.807   (   1.235   -2.025   -2.025)    3.119
   7.645   (  30.508    7.627    7.627)   32.359
   7.654   (  20.841   16.070   16.070)   30.836
  10.942   (  -6.001    1.457    1.457)    6.345
  11.081   (  -3.647    5.996    5.996)    9.230
  14.843   (  13.916    5.297    5.297)   15.804
  16.959   (  -3.004  -12.875  -12.875)   18.454
  17.146   (  -1.075   -6.989   -6.989)    9.942
  19.336   ( -20.153    4.957    4.957)   21.337
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.226   (  24.003    0.692    0.692)   24.023
   3.276   (  23.795    3.922    3.922)   24.433
   4.439   (  17.336   -0.704   -0.704)   17.364
   4.678   (   3.203    0.139    0.139)    3.209
   4.970   (   9.976   12.941   12.941)   20.844
   6.265   (  -0.829  -21.809  -21.809)   30.854
   6.818   (  -0.254   -9.047   -9.047)   12.797
   8.061   (  19.267   20.142   20.142)   34.389
   8.284   (  31.770   11.761   11.761)   35.861
  10.811   (  -6.942    0.863    0.863)    7.048
  10.987   (  -5.708    7.614    7.614)   12.187
  15.180   (  19.701    0.199    0.199)   19.703
  16.905   (  -2.737  -14.704  -14.704)   20.974
  17.125   (  -0.995   -6.930   -6.930)    9.850
  18.879   ( -24.761    8.835    8.835)   27.735
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.598   (  13.123   -0.887   -0.887)   13.182
   3.659   (  14.263    2.732    2.732)   14.777
   4.728   (  11.831   -6.395   -6.395)   14.891
   4.751   (   3.793    1.424    1.424)    4.295
   5.154   (   7.756   18.249   18.249)   26.948
   6.248   (  -0.867  -27.014  -27.014)   38.214
   6.797   (  -1.630  -17.055  -17.055)   24.174
   8.403   (  14.286   23.707   23.707)   36.444
   8.865   (  24.935   16.685   16.685)   34.330
  10.675   (  -6.169   -0.291   -0.291)    6.183
  10.857   (  -6.861    7.321    7.321)   12.421
  15.622   (  23.701   -6.664   -6.664)   25.505
  16.836   (  -4.506  -15.788  -15.788)   22.777
  17.108   (  -0.696   -6.879   -6.879)    9.753
  18.378   ( -23.923   14.290   14.290)   31.316
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.764   (   3.918   -1.960   -1.960)    4.799
   3.847   (   4.553    1.982    1.982)    5.347
   4.813   (   1.867    3.227    3.227)    4.931
   4.899   (   4.595   -8.178   -8.178)   12.444
   5.265   (   2.969   22.225   22.225)   31.571
   6.233   (  -0.439  -30.605  -30.605)   43.285
   6.766   (  -1.071  -23.164  -23.164)   32.776
   8.608   (   5.459   26.169   26.169)   37.409
   9.231   (  10.184   21.104   21.104)   31.535
  10.582   (  -2.631   -1.515   -1.515)    3.393
  10.739   (  -3.827    5.359    5.359)    8.490
  16.062   (  16.711  -15.174  -15.174)   27.198
  16.712   (  -6.979  -14.794  -14.794)   22.055
  17.098   (  -0.249   -6.850   -6.850)    9.691
  18.000   ( -11.235   20.466   20.466)   31.047
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.248   (  19.232   19.232   12.818)   30.067
   2.478   (  25.056   25.056    0.282)   35.435
   3.861   (  24.028   24.028   17.093)   38.038
   4.657   (   3.563    3.563   -0.482)    5.062
   4.683   (   1.515    1.515   -1.004)    2.367
   6.136   (  -7.407   -7.407  -15.204)   18.463
   6.882   (   6.181    6.181    6.222)   10.729
   7.183   (   8.151    8.151   23.474)   26.151
   7.692   (  26.825   26.825   -1.703)   37.975
  10.998   (  -4.641   -4.641    8.155)   10.468
  11.226   (   4.206    4.206   -2.135)    6.320
  14.791   (   9.382    9.382    5.290)   14.283
  16.772   ( -13.062  -13.062   -0.420)   18.477
  17.109   (  -1.369   -1.369  -15.794)   15.912
  19.617   (  -7.075   -7.075    2.170)   10.238
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.711   (  25.323    5.585    8.086)   27.163
   2.988   (  24.515   16.657   -1.415)   29.672
   4.232   (  12.835    9.989    5.509)   17.172
   4.663   (   0.233    1.904   -2.986)    3.549
   4.847   (  12.094    0.038    3.016)   12.465
   6.010   (  -5.608   -8.413  -18.577)   21.150
   6.912   (  -2.013   12.277    5.042)   13.424
   7.547   (  25.078   -4.634   28.617)   38.332
   8.274   (  30.930   33.164   -2.954)   45.445
  10.881   (  -6.950   -5.877    9.575)   13.211
  11.274   (   0.516   10.233   -1.436)   10.346
  14.978   (   9.534    7.518    2.045)   12.312
  16.564   (  -7.731  -22.079   -1.691)   23.455
  17.084   (  -1.251   -1.573  -16.442)   16.565
  19.390   ( -14.993   -0.515    4.780)   15.745
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.204   (  22.615   -1.434    3.483)   22.927
   3.423   (  18.411    8.969    0.275)   20.481
   4.418   (   7.081   -1.465   -2.122)    7.536
   4.704   (   4.127   -0.488   -6.101)    7.382
   5.107   (  12.513   -2.536    7.781)   14.952
   5.901   (  -5.205   -5.827  -16.798)   18.526
   6.840   (  -4.871    9.436   -7.206)   12.833
   8.069   (  25.099   -5.638   36.341)   44.525
   8.913   (  32.073   36.948   -3.426)   49.047
  10.726   (  -8.293   -7.578   10.028)   15.058
  11.246   (  -3.428   15.077   -0.715)   15.478
  15.187   (  11.645    0.574   -2.225)   11.869
  16.448   (  -4.246  -26.865   -3.805)   27.463
  17.058   (  -1.317   -1.844  -16.556)   16.710
  19.042   ( -18.999    6.204    8.580)   21.750
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.569   (  12.931   -2.098   -1.237)   13.158
   3.722   (  11.509    2.975    3.473)   12.384
   4.550   (   6.427  -10.948   -4.659)   13.524
   4.822   (   7.235    0.743   -8.979)   11.555
   5.300   (   5.280  -17.164    6.974)   19.264
   5.821   (  -1.865    8.200   -3.348)    9.051
   6.726   (  -6.127    4.857  -26.554)   27.681
   8.510   (  18.177   -3.662   42.780)   46.626
   9.523   (  27.668   39.030   -2.316)   47.898
  10.560   (  -7.756   -9.526    9.017)   15.238
  11.134   (  -7.612   17.520   -0.754)   19.117
  15.445   (  13.569  -10.401   -5.952)   18.103
  16.375   (  -3.461  -27.241   -7.244)   28.400
  17.031   (  -1.303   -2.330  -16.538)   16.752
  18.668   ( -17.139   13.346   13.605)   25.632
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.727   (   3.499   -1.826   -2.929)    4.915
   3.879   (   3.957    1.067    4.152)    5.833
   4.658   (   3.330  -14.752   -7.104)   16.709
   4.969   (   6.095    8.476   -9.081)   13.837
   5.300   (  -3.305  -44.774   -2.450)   44.963
   5.838   (   2.093   29.280   23.825)   37.807
   6.621   (  -3.456   -0.455  -47.829)   47.956
   8.769   (   6.968   -1.781   47.490)   48.032
   9.959   (  13.479   42.474   -0.575)   44.565
  10.440   (  -3.518  -11.149    6.981)   13.617
  10.974   (  -6.359   15.988   -1.724)   17.293
  15.681   (   7.890  -21.307   -6.778)   23.710
  16.307   (  -2.703  -24.334  -12.420)   27.454
  17.010   (  -0.682   -2.805  -16.265)   16.519
  18.409   (  -7.270   19.153   18.249)   27.436
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.896   (  12.621   12.621    7.063)   19.195
   3.325   (  15.878   15.878   -1.151)   22.484
   4.379   (   4.546    4.546   -2.150)    6.779
   4.580   (  -8.603   -8.603    0.365)   12.171
   4.939   (   9.567    9.567   -5.293)   14.528
   5.930   (  -0.990   -0.990  -26.615)   26.652
   7.199   (   8.761    8.761   16.313)   20.485
   7.538   (   9.400    9.400   28.918)   31.827
   8.961   (  34.901   34.901   -3.452)   49.478
  10.731   (  -8.853   -8.853   12.028)   17.362
  11.445   (   6.611    6.611   -2.026)    9.566
  15.118   (   5.878    5.878   -2.326)    8.632
  16.162   ( -16.996  -16.996   -0.211)   24.037
  17.055   (  -1.451   -1.451  -17.809)   17.927
  19.306   (  -7.529   -7.529    6.646)   12.552
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.193   (  16.206    0.727    4.822)   16.924
   3.594   (  10.662    7.359    1.354)   13.025
   4.346   (  -6.817   -9.439   -1.141)   11.699
   4.503   (   0.967  -13.554   -5.253)   14.568
   5.151   (  10.111   10.618  -10.496)   18.032
   5.953   (   3.815    5.354  -28.351)   29.103
   7.162   (  -7.531   20.877   12.782)   25.611
   7.943   (  24.404   -6.181   34.882)   43.018
   9.654   (  33.660   36.459   -3.605)   49.751
  10.533   ( -10.609  -11.359   13.367)   20.500
  11.532   (   1.736   12.872   -1.824)   13.116
  15.199   (   2.200    0.415   -6.064)    6.464
  15.907   (  -8.353  -26.134   -1.492)   27.477
  17.024   (  -1.623   -1.682  -18.569)   18.716
  19.112   ( -11.176    0.845    9.997)   15.019
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.506   (  13.709   -4.562    0.789)   14.469
   3.763   (   6.858    0.970    4.265)    8.134
   4.237   (  -3.199  -24.010   -0.798)   24.236
   4.495   (  -2.607  -18.604   -6.420)   19.852
   5.301   (   4.453   15.403  -17.152)   23.479
   6.108   (  11.638   11.005  -21.550)   26.851
   6.967   ( -11.627   19.094    1.803)   22.428
   8.395   (  19.458   -6.942   42.097)   46.893
  10.287   (  28.123   35.305   -2.975)   45.235
  10.320   (  -9.315  -12.890   12.593)   20.285
  11.502   (  -4.989   18.448   -1.875)   19.202
  15.226   (   1.056  -11.141   -5.804)   12.607
  15.806   (  -2.608  -28.702   -6.751)   29.600
  16.991   (  -1.558   -1.798  -19.393)   19.539
  18.891   (  -9.941    8.834   14.035)   19.335
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.687   (   3.979   -2.579   -2.810)    5.512
   3.882   (   4.098   -0.996    6.125)    7.436
   4.238   (   3.277  -38.384   -1.543)   38.555
   4.399   (  -6.062  -21.297   -4.783)   22.654
   5.346   (   0.726   19.306  -21.218)   28.696
   6.389   (  15.296   17.455   -7.336)   24.340
   6.703   ( -14.349   14.198  -14.508)   24.859
   8.681   (   7.924   -5.999   48.041)   49.058
  10.158   (  -4.836  -16.722   10.662)   20.413
  10.776   (  17.488   37.917   -2.676)   41.841
  11.343   (  -9.073   20.434   -2.472)   22.495
  15.249   (   0.933  -21.321   -2.670)   21.508
  15.775   (  -0.728  -27.885  -13.156)   30.841
  16.967   (  -0.715   -1.682  -19.984)   20.068
  18.745   (  -4.000   14.120   17.084)   22.522
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.243   (   4.161    4.161    4.839)    7.618
   3.708   (   4.116    4.116    2.867)    6.489
   4.016   ( -16.737  -16.737    0.073)   23.670
   4.368   (  -4.761   -4.761   -5.493)    8.689
   5.434   (  15.239   15.239  -21.878)   30.710
   6.055   (   5.562    5.562  -36.756)   37.588
   7.508   (   6.132    6.132   24.658)   26.138
   7.914   (   8.984    8.984   33.413)   35.747
  10.277   ( -12.536  -12.536   15.229)   23.371
  10.361   (  32.378   32.378   -3.017)   45.889
  11.748   (   8.431    8.431   -1.724)   12.048
  15.185   (  -3.089   -3.089  -10.979)   11.817
  15.443   ( -18.444  -18.444   -0.014)   26.083
  16.990   (  -1.698   -1.698  -19.788)   19.933
  19.089   (  -2.679   -2.679   12.373)   12.940
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.339   (   4.422  -13.776    1.679)   14.565
   3.747   (  -3.355   -6.753    1.043)    7.612
   3.795   (  -0.946   -8.478    3.234)    9.123
   4.250   (  -7.176   -9.403   -3.342)   12.291
   5.678   (   7.774   20.294  -29.921)   36.981
   6.226   (  12.101    2.327  -39.709)   41.577
   7.409   ( -10.452   22.040   24.879)   34.842
   8.255   (  18.832   -6.201   38.912)   43.672
  10.000   ( -13.549  -17.380   16.573)   27.574
  10.982   (  28.094   32.438   -3.235)   43.035
  11.851   (   1.214   15.888   -1.609)   16.016
  15.004   ( -11.803  -11.283   -4.990)   17.074
  15.257   (  -2.776  -24.921   -9.140)   26.689
  16.958   (  -1.361   -1.474  -20.817)   20.913
  19.024   (  -3.262    4.614   14.817)   15.858
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.382   (   0.237  -32.222   -2.276)   32.303
   3.712   (   0.483   -3.900   -0.533)    3.966
   3.851   (   4.046   -3.027    6.182)    7.984
   4.097   (  -6.796  -12.241   -1.197)   14.052
   5.759   (   1.407   20.838  -33.318)   39.323
   6.480   (   9.957   -1.415  -49.306)   50.322
   7.197   (  -8.340   27.526   31.267)   42.484
   8.552   (   8.880   -6.447   45.254)   46.565
   9.778   (  -6.973  -21.019   14.660)   26.558
  11.464   (  19.143   30.726   -2.718)   36.303
  11.778   (  -8.518   22.344   -1.953)   23.992
  14.835   (  -4.332  -19.917   -1.089)   20.412
  15.225   (  -0.815  -26.097  -14.715)   29.970
  16.938   (  -0.534   -1.173  -21.453)   21.492
  18.974   (  -1.402    8.861   16.390)   18.685
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.023   ( -14.352  -14.352   -1.820)   20.378
   3.646   (  -0.321   -0.321    3.523)    3.553
   3.765   (  -0.217   -0.217    3.151)    3.166
   4.070   (  -8.823   -8.823    0.015)   12.477
   6.042   (  13.028   13.028  -40.113)   44.142
   6.267   (   4.305    4.305  -44.109)   44.527
   7.723   (   6.152    6.152   33.909)   35.008
   8.228   (   6.338    6.338   36.667)   37.747
   9.639   ( -17.967  -17.967   19.545)   32.056
  11.577   (  26.457   26.457   -2.602)   37.506
  12.114   (   9.572    9.572   -1.008)   13.574
  14.721   ( -16.933  -16.933   -0.155)   23.948
  14.893   ( -10.382  -10.382  -15.993)   21.710
  16.933   (  -0.997   -0.997  -21.643)   21.688
  19.086   (   1.840    1.840   15.533)   15.749
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.825   (  -4.990  -21.477   -3.609)   22.342
   3.682   (   2.253   -0.535    1.846)    2.962
   3.782   (   2.814   -3.539    4.164)    6.146
   3.898   (  -7.917   -7.591    3.143)   11.410
   6.140   (   0.505   16.150  -43.175)   46.099
   6.427   (   6.789   -3.199  -49.871)   50.433
   7.717   (  -2.431   24.166   36.071)   43.486
   8.439   (   9.406   -4.389   40.872)   42.170
   9.328   ( -10.703  -23.406   19.545)   32.318
  12.013   (  16.153   24.007   -1.670)   28.984
  12.198   (  -2.015   18.343   -0.938)   18.477
  14.454   (  -7.606  -17.681   -0.568)   19.256
  14.763   (  -3.238  -19.176  -16.896)   25.762
  16.919   (  -0.354   -0.728  -22.041)   22.056
  19.107   (   0.432    4.597   15.944)   16.599
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.539   (  -7.119   -7.119   -4.869)   11.183
   3.684   (   0.849    0.849    3.804)    3.989
   3.723   (  -2.243   -2.243    2.367)    3.958
   3.791   (  -2.881   -2.881    5.753)    7.050
   6.364   (   3.169    3.169  -48.606)   48.812
   6.380   (   1.361    1.361  -48.153)   48.191
   8.140   (  15.011   15.011   35.200)   41.106
   8.402   (   2.229    2.229   38.290)   38.419
   8.876   ( -19.344  -19.344   24.016)   36.403
  12.393   (  12.312   12.312   -0.299)   17.414
  12.457   (   6.001    6.001    0.139)    8.488
  14.171   (  -8.714   -8.714   -0.875)   12.354
  14.489   (  -7.390   -7.390  -18.723)   21.443
  16.909   (  -0.249   -0.249  -22.244)   22.247
  19.166   (   1.389    1.389   15.703)   15.826
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.593   (  16.100   16.100   16.100)   27.886
   2.593   (  16.100   16.100   16.100)   27.886
   4.204   (  14.004   14.004   14.004)   24.255
   4.655   (   0.168    0.168    0.168)    0.292
   4.655   (   0.168    0.168    0.168)    0.292
   5.853   ( -12.532  -12.532  -12.532)   21.707
   7.122   (  14.566   14.566   14.566)   25.229
   7.700   (  16.587   16.587   16.587)   28.730
   7.700   (  16.587   16.587   16.587)   28.730
  11.162   (   1.387    1.387    1.387)    2.403
  11.162   (   1.387    1.387    1.387)    2.403
  14.914   (   6.237    6.237    6.237)   10.802
  16.760   (  -8.875   -8.875   -8.875)   15.371
  16.760   (  -8.875   -8.875   -8.875)   15.371
  19.608   (  -3.779   -3.779   -3.779)    6.546
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.910   (  18.601    9.596    9.596)   23.025
   3.036   (  22.714    8.694    8.694)   25.828
   4.346   (   1.745    4.415    4.415)    6.484
   4.621   (  -1.402   -0.698   -0.698)    1.715
   4.819   (  13.348   -5.713   -5.713)   15.603
   5.633   (  -9.648  -16.592  -16.592)   25.370
   7.232   (  -3.852   24.770   24.770)   35.241
   8.108   (  32.855   14.017   14.017)   38.372
   8.250   (  27.216   13.724   13.724)   33.427
  11.053   (  -6.538    4.219    4.219)    8.851
  11.265   (   3.060    2.598    2.598)    4.782
  15.024   (   4.797    2.072    2.072)    5.621
  16.468   ( -11.856  -10.474  -10.474)   18.973
  16.744   (  -0.845  -12.786  -12.786)   18.102
  19.462   ( -10.398    1.465    1.465)   10.603
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.279   (  16.934    3.114    3.114)   17.498
   3.470   (  19.825    5.313    5.313)   21.202
   4.369   (   1.672   -2.929   -2.929)    4.467
   4.612   (   0.911   -4.160   -4.160)    5.954
   5.115   (  15.408   -9.687   -9.687)   20.618
   5.453   (  -8.638  -18.711  -18.711)   27.835
   7.067   ( -11.246   25.034   25.034)   37.146
   8.776   (  24.850   15.504   15.504)   33.140
   8.842   (  37.624   14.821   14.821)   43.069
  10.901   (  -8.557    3.955    3.955)   10.222
  11.279   (  -1.913    6.087    6.087)    8.819
  15.123   (   5.755   -4.146   -4.146)    8.215
  16.274   (  -8.152  -15.526  -15.526)   23.422
  16.728   (  -0.770  -12.607  -12.607)   17.845
  19.216   ( -13.451    7.646    7.646)   17.259
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.553   (  10.006   -0.346   -0.346)   10.018
   3.799   (  12.316    3.929    3.929)   13.511
   4.420   (   3.074   -7.977   -7.977)   11.693
   4.676   (   5.796   -8.088   -8.088)   12.822
   5.276   (  -9.302  -19.578  -19.578)   29.208
   5.427   (  15.442  -14.324  -14.324)   25.472
   6.809   ( -13.877   26.141   26.141)   39.487
   9.229   (  19.765   17.240   17.240)   31.386
   9.551   (  31.969   16.545   16.545)   39.617
  10.717   (  -9.380    3.033    3.033)   10.314
  11.180   (  -7.907    7.953    7.953)   13.749
  15.275   (   9.461  -10.733  -10.733)   17.886
  16.114   (  -8.299  -19.530  -19.530)   28.839
  16.714   (  -0.533  -12.454  -12.454)   17.621
  18.957   ( -11.456   13.908   13.908)   22.762
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.681   (   3.029   -1.590   -1.590)    3.772
   3.956   (   3.610    3.315    3.315)    5.917
   4.471   (   1.529  -10.930  -10.930)   15.533
   4.837   (   9.187   -7.841   -7.841)   14.401
   5.081   (  -9.598  -22.345  -22.345)   33.026
   5.697   (   9.359  -22.692  -22.692)   33.428
   6.558   (  -8.898   32.811   32.811)   47.247
   9.530   (   8.860   19.482   19.482)   28.941
  10.063   (  16.461   19.323   19.323)   31.901
  10.557   (  -5.297    1.129    1.129)    5.532
  10.995   (  -8.179    6.886    6.886)   12.717
  15.472   (   8.145  -13.948  -13.948)   21.341
  15.947   (  -7.011  -23.014  -23.014)   33.293
  16.707   (  -0.190  -12.366  -12.366)   17.490
  18.790   (  -4.534   18.259   18.259)   26.217
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.095   (  10.460   10.460   12.288)   19.231
   3.314   (  16.170   16.170    0.668)   22.878
   4.362   (  -2.328   -2.328    0.337)    3.310
   4.585   (  -7.324   -7.324   -0.103)   10.359
   4.744   (   3.679    3.679  -12.497)   13.536
   5.385   (  -8.255   -8.255  -26.331)   28.803
   7.716   (  13.399   13.399   33.882)   38.821
   8.155   (  11.358   11.358   31.330)   35.208
   8.860   (  34.003   34.003   -6.427)   48.516
  11.013   (  -6.017   -6.017   14.649)   16.941
  11.384   (   6.746    6.746   -4.115)   10.389
  15.044   (  -0.314   -0.314   -5.201)    5.220
  16.156   ( -16.827  -16.827   -0.353)   23.799
  16.617   (  -2.855   -2.855  -23.975)   24.312
  19.437   (  -3.769   -3.769    5.425)    7.605
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.329   (  12.146    2.349    8.369)   14.936
   3.623   (  13.888    8.678    1.438)   16.439
   4.294   (  -4.136   -4.826   -3.974)    7.496
   4.421   (  -6.755  -17.255   -2.681)   18.723
   4.871   (   6.376   -4.574  -15.741)   17.589
   5.293   (  -0.826   -2.281  -33.958)   34.044
   7.668   ( -10.109   31.622   33.568)   47.212
   8.653   (  29.426   -8.127   35.528)   46.842
   9.545   (  33.816   36.165   -7.198)   50.032
  10.855   (  -9.398   -6.802   16.956)   20.545
  11.481   (   2.289   11.848   -2.961)   12.425
  15.020   (  -1.412   -6.297  -11.738)   13.395
  15.862   ( -12.457  -22.658   -3.166)   26.050
  16.574   (  -1.793   -4.539  -24.216)   24.702
  19.328   (  -6.571    3.506   10.402)   12.793
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.545   (   8.597   -1.253    3.142)    9.238
   3.860   (   9.827    1.347    4.805)   11.022
   4.204   (  -3.908  -19.823   -2.731)   20.389
   4.335   (  -2.758  -17.431   -9.288)   19.943
   4.948   (   1.138   -2.148  -16.829)   17.003
   5.359   (   6.777   -1.177  -42.005)   42.565
   7.440   ( -11.383   33.414   34.227)   49.169
   9.198   (  23.913   -8.387   37.591)   45.335
  10.191   (  29.936   35.660   -6.029)   46.949
  10.646   ( -11.137   -8.509   17.851)   22.695
  11.454   (  -5.326   17.136   -2.337)   18.096
  15.027   (   2.952  -13.747  -14.404)   20.129
  15.635   ( -10.809  -26.259   -9.801)   30.041
  16.542   (  -1.316   -5.714  -23.167)   23.897
  19.198   (  -5.852    9.870   15.090)   18.957
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.657   (   2.670   -1.427   -0.130)    3.030
   4.003   (   3.477    1.045    5.564)    6.644
   4.178   (   1.938  -35.461   -5.025)   35.867
   4.273   (  -3.687  -17.432   -7.687)   19.405
   4.947   (  -0.450    3.986  -18.134)   18.572
   5.493   (   4.561   -3.910  -48.384)   48.756
   7.259   (  -5.360   35.066   36.949)   51.220
   9.571   (  11.539   -6.501   39.849)   41.992
  10.443   (  -7.304  -10.815   16.520)   21.053
  10.703   (  18.894   35.691   -3.920)   40.574
  11.278   ( -10.352   19.496   -3.395)   22.334
  15.139   (   6.962  -18.808  -10.112)   22.460
  15.424   (  -9.233  -27.865  -19.640)   35.319
  16.524   (  -0.499   -6.386  -22.196)   23.101
  19.113   (  -2.303   13.644   17.993)   22.698
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.376   (   1.958    1.958    8.125)    8.584
   3.776   (   5.887    5.887    3.312)    8.960
   4.017   ( -16.754  -16.754   -0.033)   23.693
   4.235   (  -3.424   -3.424   -7.166)    8.649
   4.912   (   5.671    5.671  -27.474)   28.620
   5.282   (   1.014    1.014  -38.846)   38.872
   8.160   (   8.858    8.858   38.680)   40.658
   8.594   (  10.215   10.215   33.254)   36.256
  10.250   (  33.794   33.794   -7.310)   48.348
  10.671   ( -11.122  -11.122   21.635)   26.748
  11.696   (   8.772    8.772   -3.526)   12.897
  14.871   (  -8.110   -8.110  -19.803)   22.884
  15.442   ( -18.378  -18.378   -0.017)   25.991
  16.504   (  -2.816   -2.816  -26.944)   27.237
  19.363   (   0.246    0.246   13.598)   13.602
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.378   (  -2.207  -19.221    1.780)   19.429
   3.768   (  -5.339   -6.251    1.331)    8.328
   3.888   (   3.675   -3.557    5.364)    7.411
   4.169   (  -3.632   -4.797   -4.287)    7.388
   5.026   (   4.431    7.318  -33.472)   34.548
   5.340   (   4.718   -0.811  -45.354)   45.605
   8.072   (  -9.202   28.213   37.813)   48.067
   9.011   (  23.222   -9.003   36.182)   43.925
  10.422   ( -13.384  -13.319   23.924)   30.478
  10.882   (  28.737   32.872   -6.788)   44.187
  11.802   (   1.071   16.591   -3.239)   16.938
  14.732   (  -4.469  -15.365  -22.206)   27.371
  15.126   ( -13.673  -23.457   -3.186)   27.338
  16.457   (  -1.852   -3.176  -27.408)   27.654
  19.352   (  -1.002    5.609   16.312)   17.279
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.321   (  -2.147  -34.312   -3.930)   34.603
   3.724   (  -0.349   -1.014    1.562)    1.895
   3.977   (   4.995   -4.246    5.368)    8.473
   4.082   (  -4.975   -5.732    0.596)    7.614
   5.076   (   0.952    8.016  -34.197)   35.137
   5.436   (   3.353   -2.469  -52.137)   52.303
   7.922   (  -4.453   30.437   38.567)   49.331
   9.398   (  13.030   -9.728   38.965)   42.222
  10.170   (  -9.585  -15.818   23.386)   29.816
  11.379   (  19.929   31.316   -5.812)   37.572
  11.718   (  -9.305   23.314   -3.998)   25.418
  14.728   (   4.192  -22.202  -18.193)   29.008
  14.882   ( -10.455  -24.863  -10.085)   28.795
  16.431   (  -0.663   -3.218  -27.323)   27.520
  19.334   (  -0.549    8.527   17.852)   19.791
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.965   ( -17.785  -17.785   -3.977)   25.465
   3.748   (   0.831    0.831    6.450)    6.556
   3.841   (  -0.847   -0.847    3.915)    4.094
   4.085   (  -4.141   -4.141    1.737)    6.108
   5.178   (   6.244    6.244  -44.543)   45.410
   5.336   (   1.061    1.061  -47.489)   47.513
   8.472   (   7.724    7.724   40.245)   41.701
   8.945   (   6.973    6.973   33.933)   35.337
  10.127   ( -15.686  -15.686   27.351)   35.217
  11.493   (  27.538   27.538   -5.972)   39.399
  12.082   (  10.297   10.297   -2.168)   14.723
  14.439   ( -12.608  -12.608  -28.377)   33.513
  14.717   ( -17.243  -17.243   -0.267)   24.387
  16.408   (  -1.766   -1.766  -29.002)   29.110
  19.431   (   2.514    2.514   17.225)   17.588
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.721   (  -6.209  -23.750   -6.722)   25.452
   3.749   (  -0.365    2.216    4.701)    5.210
   3.869   (   2.611   -4.984    4.466)    7.183
   4.002   (  -3.578   -3.499    6.809)    8.450
   5.237   (   0.615    7.355  -46.263)   46.848
   5.388   (   2.536   -2.203  -52.716)   52.823
   8.469   (  -2.284   24.154   38.915)   45.859
   9.214   (  13.083   -7.821   36.002)   39.095
   9.825   ( -12.274  -18.195   28.437)   35.922
  11.955   (  17.672   25.900   -4.389)   31.660
  12.167   (  -2.867   20.187   -2.270)   20.516
  14.278   (  -3.050  -21.026  -29.166)   36.084
  14.446   (  -8.652  -18.945   -1.499)   20.881
  16.385   (  -0.580   -1.574  -29.552)   29.599
  19.461   (   0.666    4.343   17.692)   18.229
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.396   (  -8.246   -8.246   -9.151)   14.823
   3.792   (   0.702    0.702    6.731)    6.804
   3.804   (  -0.523   -0.523    5.464)    5.514
   3.936   (  -2.502   -2.502    8.304)    9.027
   5.342   (   1.848    1.848  -52.575)   52.640
   5.358   (   0.171    0.171  -52.731)   52.732
   8.864   (  12.465   12.465   37.745)   41.658
   9.136   (   2.416    2.416   34.023)   34.194
   9.463   ( -16.586  -16.586   32.264)   39.890
  12.377   (  14.411   14.411   -1.554)   20.440
  12.459   (   6.751    6.751   -0.001)    9.548
  13.955   (  -9.388   -9.388  -33.527)   36.060
  14.146   (  -9.270   -9.270   -1.537)   13.199
  16.365   (  -0.480   -0.480  -30.391)   30.398
  19.516   (   1.287    1.287   17.560)   17.654
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.370   (   9.865    9.865    9.865)   17.087
   3.370   (   9.865    9.865    9.865)   17.087
   4.337   (  -5.125   -5.125   -5.125)    8.877
   4.551   (  -4.213   -4.213   -4.213)    7.298
   4.551   (  -4.213   -4.213   -4.213)    7.298
   4.921   ( -18.162  -18.162  -18.162)   31.458
   8.418   (  25.208   25.208   25.208)   43.662
   8.768   (  17.997   17.997   17.997)   31.171
   8.768   (  17.997   17.997   17.997)   31.171
  11.280   (   2.547    2.547    2.547)    4.411
  11.280   (   2.547    2.547    2.547)    4.411
  14.908   (  -8.238   -8.238   -8.238)   14.269
  16.148   ( -11.192  -11.192  -11.192)   19.385
  16.148   ( -11.192  -11.192  -11.192)   19.385
  19.506   (   1.434    1.434    1.434)    2.485
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.500   (   7.129    7.165    7.165)   12.389
   3.664   (  14.154    3.846    3.846)   15.163
   4.203   (  -5.772   -5.058   -5.058)    9.191
   4.337   ( -14.394   -9.274   -9.274)   19.473
   4.571   (   3.268  -13.660  -13.660)   19.592
   4.719   (  -3.948  -18.436  -18.436)   26.370
   8.393   ( -10.980   36.632   36.632)   52.956
   9.327   (  26.055   12.345   12.345)   31.364
   9.404   (  39.410   12.707   12.707)   43.314
  11.130   (  -8.616    7.299    7.299)   13.446
  11.452   (   5.290    2.559    2.559)    6.410
  14.739   (  -7.530  -15.809  -15.809)   23.592
  15.748   ( -18.144   -9.982   -9.982)   22.989
  16.141   (  -0.338  -16.346  -16.346)   23.119
  19.503   (  -1.412    6.843    6.843)    9.779
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.622   (   4.519    3.583    3.583)    6.790
   3.912   (   9.257   -3.666   -3.666)   10.610
   4.083   (  -9.001  -12.816  -12.816)   20.236
   4.134   (  -1.391   -5.835   -5.835)    8.368
   4.648   (   4.097  -20.200  -20.200)   28.860
   4.667   (  -2.067  -14.761  -14.761)   20.978
   8.166   ( -10.499   36.439   36.439)   52.591
   9.805   (  21.416   12.081   12.081)   27.396
  10.149   (  34.015   13.103   13.103)   38.734
  10.936   ( -10.641    7.609    7.609)   15.134
  11.467   (  -4.367    6.316    6.316)    9.943
  14.656   (   0.219  -21.835  -21.835)   30.880
  15.405   ( -16.350  -13.768  -13.768)   25.425
  16.135   (  -0.223  -16.143  -16.143)   22.831
  19.465   (  -2.047   11.210   11.210)   15.984
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.677   (   1.213    1.501    1.501)    2.445
   3.908   (  -3.276  -25.166  -25.166)   35.741
   4.101   (   3.991    3.169    3.169)    6.001
   4.125   (   0.034   -5.539   -5.539)    7.833
   4.632   (  -1.205  -13.371  -13.371)   18.947
   4.717   (   2.156  -24.473  -24.473)   34.678
   8.009   (  -4.380   36.717   36.717)   52.110
  10.160   (  12.258   12.188   12.188)   21.151
  10.726   (  -8.619    7.450    7.450)   13.612
  10.734   (  22.433   13.826   13.826)   29.759
  11.281   ( -12.524    7.496    7.496)   16.408
  14.749   (   7.518  -27.021  -27.021)   38.945
  15.101   ( -12.611  -13.804  -13.804)   23.241
  16.132   (  -0.077  -16.029  -16.029)   22.668
  19.433   (  -0.901   13.617   13.617)   19.279
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.551   (  -2.482   -2.482    8.796)    9.471
   3.806   (   8.185    8.185   -1.430)   11.664
   4.001   ( -16.097  -16.097   -2.746)   22.929
   4.103   (  -4.210   -4.210   -5.475)    8.088
   4.425   (  -2.772   -2.772  -19.192)   19.589
   4.545   (  -3.562   -3.562  -32.840)   33.224
   8.996   (   9.343    9.343   43.265)   45.238
   9.226   (  10.343   10.343   29.567)   32.987
  10.065   (  33.972   33.972  -10.979)   49.281
  11.103   (  -7.160   -7.160   20.701)   23.045
  11.606   (   9.308    9.308   -5.452)   14.249
  14.419   ( -14.912  -14.912  -24.473)   32.306
  15.442   ( -18.296  -18.296   -0.006)   25.874
  15.957   (  -4.138   -4.138  -26.891)   27.521
  19.603   (   3.268    3.268   10.022)   11.036
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.389   ( -11.274  -28.129   -1.425)   30.338
   3.790   (  -3.132    3.912   -0.796)    5.074
   3.980   (   4.393   -3.997    2.717)    6.531
   4.080   (   2.626    0.911   -5.649)    6.295
   4.397   (  -0.292   -5.965  -26.848)   27.505
   4.496   (  -1.898   -5.840  -35.386)   35.915
   8.849   ( -10.214   30.686   38.216)   50.064
   9.703   (  24.202  -10.819   32.774)   42.153
  10.708   (  28.848   30.917   -9.526)   43.345
  10.915   ( -10.343   -6.504   23.615)   26.589
  11.725   (   1.074   16.494   -4.210)   17.057
  14.197   (  -6.092  -22.688  -29.519)   37.726
  15.079   ( -17.249  -18.159   -2.416)   25.162
  15.908   (  -1.367   -7.183  -26.435)   27.428
  19.643   (   1.043    6.693   12.395)   14.125
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.214   (  -4.881  -36.279   -7.468)   37.360
   3.757   (  -0.688    4.028    0.380)    4.104
   4.028   (   1.146   -3.903    0.747)    4.135
   4.148   (   2.221    0.980    4.082)    4.749
   4.410   (   1.175   -7.547  -32.783)   33.661
   4.461   (  -1.516   -6.952  -40.778)   41.394
   8.695   (  -4.291   31.098   37.486)   48.895
  10.105   (  14.093  -10.361   30.940)   35.542
  10.676   (  -9.856   -9.408   26.162)   29.498
  11.236   (  21.112   29.677   -7.739)   37.234
  11.623   ( -10.682   23.181   -5.088)   26.026
  14.173   (   1.946  -28.860  -31.248)   42.580
  14.790   (  -9.560  -17.439   -4.099)   20.305
  15.892   (  -0.374   -8.277  -25.608)   26.915
  19.654   (   0.166    8.461   13.631)   16.044
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.865   ( -21.741  -21.741   -5.863)   31.300
   3.880   (   1.859    1.859    4.662)    5.353
   3.907   (  -0.512   -0.512    2.068)    2.191
   4.119   (   2.090    2.090   -2.972)    4.192
   4.324   (  -1.959   -1.959  -34.742)   34.852
   4.408   (  -3.356   -3.356  -42.753)   43.016
   9.305   (   6.940    6.940   41.670)   42.811
   9.573   (   6.711    6.711   28.582)   30.116
  10.691   ( -12.664  -12.664   28.231)   33.433
  11.332   (  28.220   28.220   -9.989)   41.140
  12.026   (  11.435   11.435   -3.520)   16.551
  13.802   ( -15.175  -15.175  -33.778)   40.019
  14.711   ( -17.668  -17.668   -0.308)   24.988
  15.813   (  -2.892   -2.892  -29.721)   30.001
  19.740   (   3.161    3.161   13.184)   13.921
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.562   (  -7.900  -27.021   -8.952)   29.541
   3.849   (  -1.169    4.713    3.277)    5.858
   3.965   (   2.145   -2.413    4.586)    5.609
   4.143   (   0.259   -1.330    6.010)    6.161
   4.312   (   0.434   -2.685  -44.643)   44.726
   4.356   (  -1.727   -3.941  -47.480)   47.675
   9.246   (  -3.662   23.804   37.173)   44.293
   9.879   (  13.884  -10.547   30.368)   35.017
  10.431   ( -11.199  -14.122   31.778)   36.533
  11.819   (  19.226   27.565   -9.783)   35.002
  12.104   (  -4.448   23.151   -4.014)   23.914
  13.627   (  -2.823  -25.026  -33.436)   41.860
  14.422   (  -9.134  -18.433   -1.332)   20.615
  15.778   (  -0.806   -3.582  -30.329)   30.550
  19.779   (   0.873    4.215   13.547)   14.215
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.180   ( -10.093  -10.093  -12.167)   18.756
   3.927   (   1.195    1.195    5.161)    5.431
   3.935   (   0.702    0.702    5.570)    5.658
   4.106   (  -1.685   -1.685    8.191)    8.531
   4.286   (  -0.284   -0.284  -50.917)   50.919
   4.300   (  -1.613   -1.613  -50.713)   50.764
   9.622   (   9.340    9.340   35.774)   38.135
   9.753   (   2.236    2.236   27.460)   27.641
  10.141   ( -13.299  -13.299   36.128)   40.730
  12.303   (  18.997   18.997   -7.473)   27.886
  12.455   (   7.922    7.922   -0.518)   11.216
  13.213   ( -13.759  -13.759  -37.543)   42.286
  14.112   ( -10.081  -10.081   -1.829)   14.374
  15.734   (  -1.009   -1.009  -31.905)   31.936
  19.831   (   1.218    1.218   13.441)   13.551
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.707   (  -2.418   -2.418   -2.418)    4.188
   3.707   (  -2.418   -2.418   -2.418)    4.188
   3.724   ( -21.067  -21.067  -21.067)   36.489
   4.023   (  -3.014   -3.014   -3.014)    5.221
   4.177   (  -4.599   -4.599   -4.599)    7.966
   4.177   (  -4.599   -4.599   -4.599)    7.966
   9.774   (  14.826   14.826   14.826)   25.679
   9.774   (  14.826   14.826   14.826)   25.679
   9.891   (  22.632   22.632   22.632)   39.199
  11.477   (   4.216    4.216    4.216)    7.302
  11.477   (   4.216    4.216    4.216)    7.302
  13.938   ( -21.980  -21.980  -21.980)   38.071
  15.442   ( -12.144  -12.144  -12.144)   21.033
  15.442   ( -12.144  -12.144  -12.144)   21.033
  19.761   (   5.812    5.812    5.812)   10.067
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.233   ( -20.644  -19.037  -19.037)   33.927
   3.590   (  -4.106  -12.952  -12.952)   18.772
   3.806   (   3.122  -19.826  -19.826)   28.211
   3.979   (  -1.073   -0.468   -0.468)    1.261
   4.104   (  -2.296   -1.045   -1.045)    2.731
   4.196   (   2.514   -0.687   -0.687)    2.696
   9.553   ( -11.397   31.094   31.094)   45.427
  10.222   (  19.579    8.584    8.584)   23.037
  10.599   (  33.000    9.459    9.459)   35.609
  11.310   (  -9.014   11.152   11.152)   18.165
  11.690   (   3.740    4.845    4.845)    7.807
  13.608   (  -9.400  -27.834  -27.834)   40.469
  14.989   ( -20.503   -8.226   -8.226)   23.574
  15.442   (   0.015  -18.165  -18.165)   25.690
  19.843   (   2.613    7.549    7.549)   10.991
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.942   (  -7.625  -22.750  -22.750)   33.065
   3.540   (  -1.128  -19.431  -19.431)   27.503
   3.843   (   0.809  -24.481  -24.481)   34.631
   3.978   (   0.553   -0.744   -0.744)    1.189
   4.066   (  -1.361    1.144    1.144)    2.114
   4.240   (   1.318    2.367    2.367)    3.597
   9.384   (  -4.637   30.501   30.501)   43.384
  10.534   (  10.031    6.559    6.559)   13.662
  11.139   (  -6.506   13.087   13.087)   19.618
  11.169   (  23.195    7.837    7.837)   25.707
  11.595   ( -12.105    8.031    8.031)   16.599
  13.556   (   1.583  -29.297  -29.297)   41.463
  14.662   ( -10.157   -9.653   -9.653)   17.016
  15.443   (   0.007  -18.138  -18.138)   25.651
  19.874   (   0.686    8.396    8.396)   11.894
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.728   ( -24.698  -24.698   -8.178)   35.872
   3.461   (  -7.495   -7.495  -37.872)   39.327
   3.473   (  -7.974   -7.974  -34.853)   36.632
   3.997   (   1.576    1.576    3.243)    3.935
   4.115   (   0.242    0.242    5.944)    5.954
   4.223   (   2.247    2.247    1.959)    3.734
  10.079   (   4.717    4.717   33.520)   34.178
  10.087   (   5.713    5.713   22.631)   24.029
  11.044   (  22.275   22.275   -8.797)   32.707
  11.313   (  -3.652   -3.652   24.940)   25.469
  11.942   (  12.965   12.965   -4.814)   18.957
  13.146   ( -16.031  -16.031  -30.150)   37.722
  14.705   ( -18.102  -18.102   -0.275)   25.601
  15.232   (  -4.417   -4.417  -28.033)   28.721
  19.952   (   3.611    3.611    8.047)    9.531
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 1099
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.366   (  -9.861  -29.204  -10.613)   32.600
   3.298   (  -4.126   -6.473  -45.090)   45.739
   3.434   (  -0.770  -15.810  -40.644)   43.617
   4.031   (   1.495    3.746    5.461)    6.788
   4.103   (  -1.026    1.918    5.882)    6.272
   4.262   (   1.212   -0.090    4.309)    4.477
   9.905   (  -5.078   21.375   27.682)   35.340
  10.418   (  10.242  -11.441   22.485)   27.228
  11.101   (  -5.473   -8.106   34.551)   35.908
  11.548   (  16.027   21.569  -16.554)   31.561
  12.010   (  -5.413   26.489   -5.079)   27.509
  13.024   (   1.407  -23.422  -24.950)   34.250
  14.388   ( -10.415  -16.817   -2.267)   19.910
  15.189   (  -0.783   -7.598  -27.989)   29.013
  19.997   (   1.001    4.069    8.240)    9.244
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.923   ( -12.548  -12.548  -13.057)   22.031
   3.219   (  -3.581   -3.581  -50.911)   51.162
   3.220   (  -3.669   -3.669  -50.328)   50.594
   4.098   (   2.564    2.564    7.163)    8.029
   4.127   (   0.176    0.176    7.659)    7.663
   4.249   (  -0.726   -0.726    5.700)    5.792
  10.223   (   2.456    2.456   23.831)   24.083
  10.233   (   1.719    1.719   20.404)   20.548
  10.935   (  -6.286   -6.286   42.199)   43.126
  11.936   (  14.704   14.704  -30.897)   37.243
  12.437   (   9.383    9.383   -1.184)   13.322
  12.655   (  -9.159   -9.159  -14.620)   19.532
  14.076   ( -10.979  -10.979   -1.695)   15.619
  15.102   (  -1.815   -1.815  -30.903)   31.009
  20.047   (   1.139    1.139    8.089)    8.248
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.469   ( -19.870  -19.870  -19.870)   34.417
   2.869   ( -20.500  -20.500  -20.500)   35.507
   2.869   ( -20.500  -20.500  -20.500)   35.507
   4.072   (   3.916    3.916    3.916)    6.782
   4.240   (   2.551    2.551    2.551)    4.418
   4.240   (   2.551    2.551    2.551)    4.418
  10.475   (   7.930    7.930    7.930)   13.736
  10.475   (   7.930    7.930    7.930)   13.736
  10.961   (   9.841    9.841    9.841)   17.045
  11.839   (   8.232    8.232    8.232)   14.258
  11.839   (   8.232    8.232    8.232)   14.258
  12.661   ( -15.496  -15.496  -15.496)   26.841
  14.700   ( -12.361  -12.361  -12.361)   21.409
  14.700   ( -12.361  -12.361  -12.361)   21.409
  20.069   (   3.818    3.818    3.818)    6.613
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.096   ( -11.800  -19.077  -19.077)   29.447
   2.531   (  -8.809  -26.554  -26.554)   38.572
   2.708   (  -2.894  -30.563  -30.563)   43.319
   4.145   (   3.112    5.249    5.249)    8.049
   4.213   (  -2.004    4.584    4.584)    6.786
   4.317   (   2.012    1.293    1.293)    2.718
  10.348   (  -5.478   16.759   16.759)   24.326
  10.700   (   6.110    2.145    2.145)    6.822
  11.238   (  10.733   -6.526   -6.526)   14.155
  11.754   (  -2.611   16.917   16.917)   24.066
  11.974   (  -4.373   11.671   11.671)   17.075
  12.649   (   8.340  -13.014  -13.014)   20.206
  14.301   ( -12.708   -8.287   -8.287)   17.287
  14.698   (  -0.065  -18.638  -18.638)   26.359
  20.116   (   1.043    3.865    3.865)    5.564
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 595
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.677   ( -14.932  -14.932  -10.898)   23.763
   2.256   (  -7.699   -7.699  -43.611)   44.950
   2.272   (  -9.057   -9.057  -42.788)   44.664
   4.231   (   3.183    3.183    5.683)    7.249
   4.266   (   0.257    0.257    5.759)    5.771
   4.329   (   0.322    0.322    2.285)    2.330
  10.568   (   0.971    0.971   12.850)   12.923
  10.595   (  -1.312   -1.312   13.506)   13.633
  11.268   (   3.610    3.610  -30.064)   30.494
  11.768   (  -2.218   -2.218   38.707)   38.833
  12.410   (  10.879   10.879   -1.417)   15.450
  12.500   (   2.332    2.332   -3.948)    5.144
  14.046   ( -11.753  -11.753   -1.184)   16.663
  14.506   (  -2.899   -2.899  -28.110)   28.408
  20.162   (   1.040    1.040    3.698)    3.980
======================= 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 5.83e-04 3.36e-05 
Number of triplets: 231
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.483   ( -10.764  -10.764  -10.764)   18.644
   1.604   ( -17.315  -17.315  -17.315)   29.990
   1.604   ( -17.315  -17.315  -17.315)   29.990
   4.317   (   2.794    2.794    2.794)    4.839
   4.347   (   0.901    0.901    0.901)    1.560
   4.347   (   0.901    0.901    0.901)    1.560
  10.743   (   1.717    1.717    1.717)    2.974
  10.743   (   1.717    1.717    1.717)    2.974
  10.867   (  -6.032   -6.032   -6.032)   10.447
  12.387   (   7.721    7.721    7.721)   13.374
  12.387   (   7.721    7.721    7.721)   13.374
  12.463   (   3.047    3.047    3.047)    5.278
  14.030   (  -8.275   -8.275   -8.275)   14.333
  14.030   (  -8.275   -8.275   -8.275)   14.333
  20.207   (   0.964    0.964    0.964)    1.669
=================== 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    771.010    771.010    771.010     -0.000     -0.000      0.000 3/32955
   20.0    207.346    207.346    207.346     -0.000     -0.000      0.000 3/32955
   30.0     70.259     70.259     70.259     -0.000     -0.000      0.000 3/32955
   40.0     42.555     42.555     42.555     -0.000     -0.000      0.000 3/32955
   50.0     31.168     31.168     31.168     -0.000      0.000      0.000 3/32955
   60.0     24.678     24.678     24.678     -0.000      0.000      0.000 3/32955
   70.0     20.485     20.485     20.485     -0.000      0.000      0.000 3/32955
   80.0     17.588     17.588     17.588     -0.000      0.000      0.000 3/32955
   90.0     15.487     15.487     15.487      0.000      0.000      0.000 3/32955
  100.0     13.899     13.899     13.899      0.000      0.000      0.000 3/32955
  110.0     12.656     12.656     12.656      0.000      0.000      0.000 3/32955
  120.0     11.652     11.652     11.652      0.000      0.000      0.000 3/32955
  130.0     10.820     10.820     10.820      0.000      0.000      0.000 3/32955
  140.0     10.116     10.116     10.116      0.000      0.000      0.000 3/32955
  150.0      9.509      9.509      9.509      0.000      0.000      0.000 3/32955
  160.0      8.979      8.979      8.979      0.000      0.000      0.000 3/32955
  170.0      8.509      8.509      8.509      0.000      0.000      0.000 3/32955
  180.0      8.091      8.091      8.091      0.000      0.000      0.000 3/32955
  190.0      7.713      7.713      7.713      0.000      0.000      0.000 3/32955
  200.0      7.372      7.372      7.372      0.000      0.000      0.000 3/32955
  210.0      7.060      7.060      7.060      0.000      0.000      0.000 3/32955
  220.0      6.774      6.774      6.774      0.000      0.000      0.000 3/32955
  230.0      6.511      6.511      6.511      0.000      0.000      0.000 3/32955
  240.0      6.268      6.268      6.268      0.000      0.000      0.000 3/32955
  250.0      6.042      6.042      6.042      0.000      0.000      0.000 3/32955
  260.0      5.833      5.833      5.833      0.000      0.000      0.000 3/32955
  270.0      5.637      5.637      5.637      0.000      0.000      0.000 3/32955
  280.0      5.454      5.454      5.454      0.000      0.000      0.000 3/32955
  290.0      5.283      5.283      5.283      0.000      0.000      0.000 3/32955
  300.0      5.122      5.122      5.122      0.000      0.000      0.000 3/32955
  310.0      4.970      4.970      4.970      0.000      0.000      0.000 3/32955
  320.0      4.827      4.827      4.827      0.000      0.000      0.000 3/32955
  330.0      4.692      4.692      4.692      0.000      0.000      0.000 3/32955
  340.0      4.565      4.565      4.565      0.000      0.000      0.000 3/32955
  350.0      4.444      4.444      4.444      0.000      0.000      0.000 3/32955
  360.0      4.329      4.329      4.329      0.000      0.000      0.000 3/32955
  370.0      4.220      4.220      4.220      0.000      0.000      0.000 3/32955
  380.0      4.116      4.116      4.116      0.000      0.000      0.000 3/32955
  390.0      4.017      4.017      4.017      0.000      0.000      0.000 3/32955
  400.0      3.923      3.923      3.923      0.000      0.000      0.000 3/32955
  410.0      3.833      3.833      3.833      0.000      0.000      0.000 3/32955
  420.0      3.747      3.747      3.747      0.000      0.000      0.000 3/32955
  430.0      3.665      3.665      3.665      0.000      0.000      0.000 3/32955
  440.0      3.586      3.586      3.586      0.000      0.000      0.000 3/32955
  450.0      3.510      3.510      3.510      0.000      0.000      0.000 3/32955
  460.0      3.438      3.438      3.438      0.000      0.000      0.000 3/32955
  470.0      3.368      3.368      3.368      0.000      0.000      0.000 3/32955
  480.0      3.302      3.302      3.302      0.000      0.000      0.000 3/32955
  490.0      3.237      3.237      3.237      0.000      0.000      0.000 3/32955
  500.0      3.176      3.176      3.176      0.000      0.000      0.000 3/32955
  510.0      3.116      3.116      3.116      0.000      0.000      0.000 3/32955
  520.0      3.059      3.059      3.059      0.000      0.000      0.000 3/32955
  530.0      3.003      3.003      3.003      0.000      0.000      0.000 3/32955
  540.0      2.950      2.950      2.950      0.000      0.000      0.000 3/32955
  550.0      2.899      2.899      2.899      0.000      0.000      0.000 3/32955
  560.0      2.849      2.849      2.849      0.000      0.000      0.000 3/32955
  570.0      2.801      2.801      2.801      0.000      0.000      0.000 3/32955
  580.0      2.754      2.754      2.754      0.000      0.000      0.000 3/32955
  590.0      2.709      2.709      2.709      0.000      0.000      0.000 3/32955
  600.0      2.666      2.666      2.666      0.000      0.000      0.000 3/32955
  610.0      2.624      2.624      2.624      0.000      0.000      0.000 3/32955
  620.0      2.583      2.583      2.583      0.000      0.000      0.000 3/32955
  630.0      2.543      2.543      2.543      0.000      0.000      0.000 3/32955
  640.0      2.505      2.505      2.505      0.000      0.000      0.000 3/32955
  650.0      2.467      2.467      2.467      0.000      0.000      0.000 3/32955
  660.0      2.431      2.431      2.431      0.000      0.000      0.000 3/32955
  670.0      2.396      2.396      2.396      0.000      0.000      0.000 3/32955
  680.0      2.362      2.362      2.362      0.000      0.000      0.000 3/32955
  690.0      2.329      2.329      2.329      0.000      0.000      0.000 3/32955
  700.0      2.296      2.296      2.296      0.000      0.000      0.000 3/32955
  710.0      2.265      2.265      2.265      0.000      0.000      0.000 3/32955
  720.0      2.234      2.234      2.234      0.000      0.000      0.000 3/32955
  730.0      2.204      2.204      2.204      0.000      0.000      0.000 3/32955
  740.0      2.175      2.175      2.175      0.000      0.000      0.000 3/32955
  750.0      2.147      2.147      2.147      0.000      0.000      0.000 3/32955
  760.0      2.120      2.120      2.120      0.000      0.000      0.000 3/32955
  770.0      2.093      2.093      2.093      0.000      0.000      0.000 3/32955
  780.0      2.067      2.067      2.067      0.000      0.000      0.000 3/32955
  790.0      2.041      2.041      2.041      0.000      0.000      0.000 3/32955
  800.0      2.016      2.016      2.016      0.000      0.000      0.000 3/32955
  810.0      1.992      1.992      1.992      0.000      0.000      0.000 3/32955
  820.0      1.968      1.968      1.968      0.000      0.000      0.000 3/32955
  830.0      1.945      1.945      1.945      0.000      0.000      0.000 3/32955
  840.0      1.922      1.922      1.922      0.000      0.000      0.000 3/32955
  850.0      1.900      1.900      1.900      0.000      0.000      0.000 3/32955
  860.0      1.879      1.879      1.879      0.000      0.000      0.000 3/32955
  870.0      1.857      1.857      1.857      0.000      0.000      0.000 3/32955
  880.0      1.837      1.837      1.837      0.000      0.000      0.000 3/32955
  890.0      1.816      1.816      1.816      0.000      0.000      0.000 3/32955
  900.0      1.797      1.797      1.797      0.000      0.000      0.000 3/32955
  910.0      1.777      1.777      1.777      0.000      0.000      0.000 3/32955
  920.0      1.758      1.758      1.758      0.000      0.000      0.000 3/32955
  930.0      1.740      1.740      1.740      0.000      0.000      0.000 3/32955
  940.0      1.722      1.722      1.722      0.000      0.000      0.000 3/32955
  950.0      1.704      1.704      1.704      0.000      0.000      0.000 3/32955
  960.0      1.686      1.686      1.686      0.000      0.000      0.000 3/32955
  970.0      1.669      1.669      1.669      0.000      0.000      0.000 3/32955
  980.0      1.653      1.653      1.653      0.000      0.000      0.000 3/32955
  990.0      1.636      1.636      1.636      0.000      0.000      0.000 3/32955
 1000.0      1.620      1.620      1.620      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-09 04:50:08]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

