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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [3 3 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    4.356066310000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356066310000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356066310000000
Atomic positions (fractional):
   *1 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036
   *2 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383
   *3 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999
    4 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999
    5 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.356066310000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356066310000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356066310000000
Atomic positions (fractional):
   *1 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036 > 1
   *2 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383 > 2
   *3 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999 > 3
    4 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999 > 4
    5 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   13.068198930000001    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.068198930000001    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.068198930000001
Atomic positions (fractional):
   *1 Pa  0.16666666666667  0.16666666666667  0.16666666666667 231.036 > 1
    2 Pa  0.50000000000000  0.16666666666667  0.16666666666667 231.036 > 1
    3 Pa  0.83333333333333  0.16666666666667  0.16666666666667 231.036 > 1
    4 Pa  0.16666666666667  0.50000000000000  0.16666666666667 231.036 > 1
    5 Pa  0.50000000000000  0.50000000000000  0.16666666666667 231.036 > 1
    6 Pa  0.83333333333333  0.50000000000000  0.16666666666667 231.036 > 1
    7 Pa  0.16666666666667  0.83333333333333  0.16666666666667 231.036 > 1
    8 Pa  0.50000000000000  0.83333333333333  0.16666666666667 231.036 > 1
    9 Pa  0.83333333333333  0.83333333333333  0.16666666666667 231.036 > 1
   10 Pa  0.16666666666667  0.16666666666667  0.50000000000000 231.036 > 1
   11 Pa  0.50000000000000  0.16666666666667  0.50000000000000 231.036 > 1
   12 Pa  0.83333333333333  0.16666666666667  0.50000000000000 231.036 > 1
   13 Pa  0.16666666666667  0.50000000000000  0.50000000000000 231.036 > 1
   14 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036 > 1
   15 Pa  0.83333333333333  0.50000000000000  0.50000000000000 231.036 > 1
   16 Pa  0.16666666666667  0.83333333333333  0.50000000000000 231.036 > 1
   17 Pa  0.50000000000000  0.83333333333333  0.50000000000000 231.036 > 1
   18 Pa  0.83333333333333  0.83333333333333  0.50000000000000 231.036 > 1
   19 Pa  0.16666666666667  0.16666666666667  0.83333333333333 231.036 > 1
   20 Pa  0.50000000000000  0.16666666666667  0.83333333333333 231.036 > 1
   21 Pa  0.83333333333333  0.16666666666667  0.83333333333333 231.036 > 1
   22 Pa  0.16666666666667  0.50000000000000  0.83333333333333 231.036 > 1
   23 Pa  0.50000000000000  0.50000000000000  0.83333333333333 231.036 > 1
   24 Pa  0.83333333333333  0.50000000000000  0.83333333333333 231.036 > 1
   25 Pa  0.16666666666667  0.83333333333333  0.83333333333333 231.036 > 1
   26 Pa  0.50000000000000  0.83333333333333  0.83333333333333 231.036 > 1
   27 Pa  0.83333333333333  0.83333333333333  0.83333333333333 231.036 > 1
  *28 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383 > 2
   29 Tl  0.33333333333333  0.00000000000000  0.00000000000000 204.383 > 2
   30 Tl  0.66666666666667  0.00000000000000  0.00000000000000 204.383 > 2
   31 Tl  0.00000000000000  0.33333333333333  0.00000000000000 204.383 > 2
   32 Tl  0.33333333333333  0.33333333333333  0.00000000000000 204.383 > 2
   33 Tl  0.66666666666667  0.33333333333333  0.00000000000000 204.383 > 2
   34 Tl  0.00000000000000  0.66666666666667  0.00000000000000 204.383 > 2
   35 Tl  0.33333333333333  0.66666666666667  0.00000000000000 204.383 > 2
   36 Tl  0.66666666666667  0.66666666666667  0.00000000000000 204.383 > 2
   37 Tl  0.00000000000000  0.00000000000000  0.33333333333333 204.383 > 2
   38 Tl  0.33333333333333  0.00000000000000  0.33333333333333 204.383 > 2
   39 Tl  0.66666666666667  0.00000000000000  0.33333333333333 204.383 > 2
   40 Tl  0.00000000000000  0.33333333333333  0.33333333333333 204.383 > 2
   41 Tl  0.33333333333333  0.33333333333333  0.33333333333333 204.383 > 2
   42 Tl  0.66666666666667  0.33333333333333  0.33333333333333 204.383 > 2
   43 Tl  0.00000000000000  0.66666666666667  0.33333333333333 204.383 > 2
   44 Tl  0.33333333333333  0.66666666666667  0.33333333333333 204.383 > 2
   45 Tl  0.66666666666667  0.66666666666667  0.33333333333333 204.383 > 2
   46 Tl  0.00000000000000  0.00000000000000  0.66666666666667 204.383 > 2
   47 Tl  0.33333333333333  0.00000000000000  0.66666666666667 204.383 > 2
   48 Tl  0.66666666666667  0.00000000000000  0.66666666666667 204.383 > 2
   49 Tl  0.00000000000000  0.33333333333333  0.66666666666667 204.383 > 2
   50 Tl  0.33333333333333  0.33333333333333  0.66666666666667 204.383 > 2
   51 Tl  0.66666666666667  0.33333333333333  0.66666666666667 204.383 > 2
   52 Tl  0.00000000000000  0.66666666666667  0.66666666666667 204.383 > 2
   53 Tl  0.33333333333333  0.66666666666667  0.66666666666667 204.383 > 2
   54 Tl  0.66666666666667  0.66666666666667  0.66666666666667 204.383 > 2
  *55 O   0.16666666666667  0.16666666666667  0.00000000000000  15.999 > 3
   56 O   0.50000000000000  0.16666666666667  0.00000000000000  15.999 > 3
   57 O   0.83333333333333  0.16666666666667  0.00000000000000  15.999 > 3
   58 O   0.16666666666667  0.50000000000000  0.00000000000000  15.999 > 3
   59 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999 > 3
   60 O   0.83333333333333  0.50000000000000  0.00000000000000  15.999 > 3
   61 O   0.16666666666667  0.83333333333333  0.00000000000000  15.999 > 3
   62 O   0.50000000000000  0.83333333333333  0.00000000000000  15.999 > 3
   63 O   0.83333333333333  0.83333333333333  0.00000000000000  15.999 > 3
   64 O   0.16666666666667  0.16666666666667  0.33333333333333  15.999 > 3
   65 O   0.50000000000000  0.16666666666667  0.33333333333333  15.999 > 3
   66 O   0.83333333333333  0.16666666666667  0.33333333333333  15.999 > 3
   67 O   0.16666666666667  0.50000000000000  0.33333333333333  15.999 > 3
   68 O   0.50000000000000  0.50000000000000  0.33333333333333  15.999 > 3
   69 O   0.83333333333333  0.50000000000000  0.33333333333333  15.999 > 3
   70 O   0.16666666666667  0.83333333333333  0.33333333333333  15.999 > 3
   71 O   0.50000000000000  0.83333333333333  0.33333333333333  15.999 > 3
   72 O   0.83333333333333  0.83333333333333  0.33333333333333  15.999 > 3
   73 O   0.16666666666667  0.16666666666667  0.66666666666667  15.999 > 3
   74 O   0.50000000000000  0.16666666666667  0.66666666666667  15.999 > 3
   75 O   0.83333333333333  0.16666666666667  0.66666666666667  15.999 > 3
   76 O   0.16666666666667  0.50000000000000  0.66666666666667  15.999 > 3
   77 O   0.50000000000000  0.50000000000000  0.66666666666667  15.999 > 3
   78 O   0.83333333333333  0.50000000000000  0.66666666666667  15.999 > 3
   79 O   0.16666666666667  0.83333333333333  0.66666666666667  15.999 > 3
   80 O   0.50000000000000  0.83333333333333  0.66666666666667  15.999 > 3
   81 O   0.83333333333333  0.83333333333333  0.66666666666667  15.999 > 3
   82 O   0.16666666666667  0.00000000000000  0.16666666666667  15.999 > 4
   83 O   0.50000000000000  0.00000000000000  0.16666666666667  15.999 > 4
   84 O   0.83333333333333  0.00000000000000  0.16666666666667  15.999 > 4
   85 O   0.16666666666667  0.33333333333333  0.16666666666667  15.999 > 4
   86 O   0.50000000000000  0.33333333333333  0.16666666666667  15.999 > 4
   87 O   0.83333333333333  0.33333333333333  0.16666666666667  15.999 > 4
   88 O   0.16666666666667  0.66666666666667  0.16666666666667  15.999 > 4
   89 O   0.50000000000000  0.66666666666667  0.16666666666667  15.999 > 4
   90 O   0.83333333333333  0.66666666666667  0.16666666666667  15.999 > 4
   91 O   0.16666666666667  0.00000000000000  0.50000000000000  15.999 > 4
   92 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999 > 4
   93 O   0.83333333333333  0.00000000000000  0.50000000000000  15.999 > 4
   94 O   0.16666666666667  0.33333333333333  0.50000000000000  15.999 > 4
   95 O   0.50000000000000  0.33333333333333  0.50000000000000  15.999 > 4
   96 O   0.83333333333333  0.33333333333333  0.50000000000000  15.999 > 4
   97 O   0.16666666666667  0.66666666666667  0.50000000000000  15.999 > 4
   98 O   0.50000000000000  0.66666666666667  0.50000000000000  15.999 > 4
   99 O   0.83333333333333  0.66666666666667  0.50000000000000  15.999 > 4
  100 O   0.16666666666667  0.00000000000000  0.83333333333333  15.999 > 4
  101 O   0.50000000000000  0.00000000000000  0.83333333333333  15.999 > 4
  102 O   0.83333333333333  0.00000000000000  0.83333333333333  15.999 > 4
  103 O   0.16666666666667  0.33333333333333  0.83333333333333  15.999 > 4
  104 O   0.50000000000000  0.33333333333333  0.83333333333333  15.999 > 4
  105 O   0.83333333333333  0.33333333333333  0.83333333333333  15.999 > 4
  106 O   0.16666666666667  0.66666666666667  0.83333333333333  15.999 > 4
  107 O   0.50000000000000  0.66666666666667  0.83333333333333  15.999 > 4
  108 O   0.83333333333333  0.66666666666667  0.83333333333333  15.999 > 4
  109 O   0.00000000000000  0.16666666666667  0.16666666666667  15.999 > 5
  110 O   0.33333333333333  0.16666666666667  0.16666666666667  15.999 > 5
  111 O   0.66666666666667  0.16666666666667  0.16666666666667  15.999 > 5
  112 O   0.00000000000000  0.50000000000000  0.16666666666667  15.999 > 5
  113 O   0.33333333333333  0.50000000000000  0.16666666666667  15.999 > 5
  114 O   0.66666666666667  0.50000000000000  0.16666666666667  15.999 > 5
  115 O   0.00000000000000  0.83333333333333  0.16666666666667  15.999 > 5
  116 O   0.33333333333333  0.83333333333333  0.16666666666667  15.999 > 5
  117 O   0.66666666666667  0.83333333333333  0.16666666666667  15.999 > 5
  118 O   0.00000000000000  0.16666666666667  0.50000000000000  15.999 > 5
  119 O   0.33333333333333  0.16666666666667  0.50000000000000  15.999 > 5
  120 O   0.66666666666667  0.16666666666667  0.50000000000000  15.999 > 5
  121 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999 > 5
  122 O   0.33333333333333  0.50000000000000  0.50000000000000  15.999 > 5
  123 O   0.66666666666667  0.50000000000000  0.50000000000000  15.999 > 5
  124 O   0.00000000000000  0.83333333333333  0.50000000000000  15.999 > 5
  125 O   0.33333333333333  0.83333333333333  0.50000000000000  15.999 > 5
  126 O   0.66666666666667  0.83333333333333  0.50000000000000  15.999 > 5
  127 O   0.00000000000000  0.16666666666667  0.83333333333333  15.999 > 5
  128 O   0.33333333333333  0.16666666666667  0.83333333333333  15.999 > 5
  129 O   0.66666666666667  0.16666666666667  0.83333333333333  15.999 > 5
  130 O   0.00000000000000  0.50000000000000  0.83333333333333  15.999 > 5
  131 O   0.33333333333333  0.50000000000000  0.83333333333333  15.999 > 5
  132 O   0.66666666666667  0.50000000000000  0.83333333333333  15.999 > 5
  133 O   0.00000000000000  0.83333333333333  0.83333333333333  15.999 > 5
  134 O   0.33333333333333  0.83333333333333  0.83333333333333  15.999 > 5
  135 O   0.66666666666667  0.83333333333333  0.83333333333333  15.999 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            6.0396110    0.0000000    0.0000000
            0.0000000    6.0396110    0.0000000
            0.0000000    0.0000000    6.0396110
-------------------------- Born effective charges --------------------------
    1 Pa    6.9426443    0.0000000    0.0000000
            0.0000000    6.9426443    0.0000000
            0.0000000    0.0000000    6.9426443
    2 Tl    1.9946292    0.0000000    0.0000000
            0.0000000    1.9946292    0.0000000
            0.0000000    0.0000000    1.9946292
    3 O    -1.6778173    0.0000000    0.0000000
            0.0000000   -1.6778173    0.0000000
            0.0000000    0.0000000   -5.5816389
    4 O    -1.6778173    0.0000000    0.0000000
            0.0000000   -5.5816389    0.0000000
            0.0000000    0.0000000   -1.6778173
    5 O    -5.5816389    0.0000000    0.0000000
            0.0000000   -1.6778173    0.0000000
            0.0000000    0.0000000   -1.6778173
----------------------------------------------------------------------------
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.004
Solver_block: 80 / 80
 - Time: 0.079
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.084
--------------------------------- 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:07:24]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-09 04:07:24]-------------------------
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    4.356066310000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356066310000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356066310000000
Atomic positions (fractional):
    1 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036
    2 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383
    3 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999
    4 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999
    5 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.068198930000001    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.068198930000001    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.068198930000001
Atomic positions (fractional):
    1 Pa  0.16666666666667  0.16666666666667  0.16666666666667 231.036 > 1
    2 Pa  0.50000000000000  0.16666666666667  0.16666666666667 231.036 > 1
    3 Pa  0.83333333333333  0.16666666666667  0.16666666666667 231.036 > 1
    4 Pa  0.16666666666667  0.50000000000000  0.16666666666667 231.036 > 1
    5 Pa  0.50000000000000  0.50000000000000  0.16666666666667 231.036 > 1
    6 Pa  0.83333333333333  0.50000000000000  0.16666666666667 231.036 > 1
    7 Pa  0.16666666666667  0.83333333333333  0.16666666666667 231.036 > 1
    8 Pa  0.50000000000000  0.83333333333333  0.16666666666667 231.036 > 1
    9 Pa  0.83333333333333  0.83333333333333  0.16666666666667 231.036 > 1
   10 Pa  0.16666666666667  0.16666666666667  0.50000000000000 231.036 > 1
   11 Pa  0.50000000000000  0.16666666666667  0.50000000000000 231.036 > 1
   12 Pa  0.83333333333333  0.16666666666667  0.50000000000000 231.036 > 1
   13 Pa  0.16666666666667  0.50000000000000  0.50000000000000 231.036 > 1
   14 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036 > 1
   15 Pa  0.83333333333333  0.50000000000000  0.50000000000000 231.036 > 1
   16 Pa  0.16666666666667  0.83333333333333  0.50000000000000 231.036 > 1
   17 Pa  0.50000000000000  0.83333333333333  0.50000000000000 231.036 > 1
   18 Pa  0.83333333333333  0.83333333333333  0.50000000000000 231.036 > 1
   19 Pa  0.16666666666667  0.16666666666667  0.83333333333333 231.036 > 1
   20 Pa  0.50000000000000  0.16666666666667  0.83333333333333 231.036 > 1
   21 Pa  0.83333333333333  0.16666666666667  0.83333333333333 231.036 > 1
   22 Pa  0.16666666666667  0.50000000000000  0.83333333333333 231.036 > 1
   23 Pa  0.50000000000000  0.50000000000000  0.83333333333333 231.036 > 1
   24 Pa  0.83333333333333  0.50000000000000  0.83333333333333 231.036 > 1
   25 Pa  0.16666666666667  0.83333333333333  0.83333333333333 231.036 > 1
   26 Pa  0.50000000000000  0.83333333333333  0.83333333333333 231.036 > 1
   27 Pa  0.83333333333333  0.83333333333333  0.83333333333333 231.036 > 1
   28 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383 > 28
   29 Tl  0.33333333333333  0.00000000000000  0.00000000000000 204.383 > 28
   30 Tl  0.66666666666667  0.00000000000000  0.00000000000000 204.383 > 28
   31 Tl  0.00000000000000  0.33333333333333  0.00000000000000 204.383 > 28
   32 Tl  0.33333333333333  0.33333333333333  0.00000000000000 204.383 > 28
   33 Tl  0.66666666666667  0.33333333333333  0.00000000000000 204.383 > 28
   34 Tl  0.00000000000000  0.66666666666667  0.00000000000000 204.383 > 28
   35 Tl  0.33333333333333  0.66666666666667  0.00000000000000 204.383 > 28
   36 Tl  0.66666666666667  0.66666666666667  0.00000000000000 204.383 > 28
   37 Tl  0.00000000000000  0.00000000000000  0.33333333333333 204.383 > 28
   38 Tl  0.33333333333333  0.00000000000000  0.33333333333333 204.383 > 28
   39 Tl  0.66666666666667  0.00000000000000  0.33333333333333 204.383 > 28
   40 Tl  0.00000000000000  0.33333333333333  0.33333333333333 204.383 > 28
   41 Tl  0.33333333333333  0.33333333333333  0.33333333333333 204.383 > 28
   42 Tl  0.66666666666667  0.33333333333333  0.33333333333333 204.383 > 28
   43 Tl  0.00000000000000  0.66666666666667  0.33333333333333 204.383 > 28
   44 Tl  0.33333333333333  0.66666666666667  0.33333333333333 204.383 > 28
   45 Tl  0.66666666666667  0.66666666666667  0.33333333333333 204.383 > 28
   46 Tl  0.00000000000000  0.00000000000000  0.66666666666667 204.383 > 28
   47 Tl  0.33333333333333  0.00000000000000  0.66666666666667 204.383 > 28
   48 Tl  0.66666666666667  0.00000000000000  0.66666666666667 204.383 > 28
   49 Tl  0.00000000000000  0.33333333333333  0.66666666666667 204.383 > 28
   50 Tl  0.33333333333333  0.33333333333333  0.66666666666667 204.383 > 28
   51 Tl  0.66666666666667  0.33333333333333  0.66666666666667 204.383 > 28
   52 Tl  0.00000000000000  0.66666666666667  0.66666666666667 204.383 > 28
   53 Tl  0.33333333333333  0.66666666666667  0.66666666666667 204.383 > 28
   54 Tl  0.66666666666667  0.66666666666667  0.66666666666667 204.383 > 28
   55 O   0.16666666666667  0.16666666666667  0.00000000000000  15.999 > 55
   56 O   0.50000000000000  0.16666666666667  0.00000000000000  15.999 > 55
   57 O   0.83333333333333  0.16666666666667  0.00000000000000  15.999 > 55
   58 O   0.16666666666667  0.50000000000000  0.00000000000000  15.999 > 55
   59 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999 > 55
   60 O   0.83333333333333  0.50000000000000  0.00000000000000  15.999 > 55
   61 O   0.16666666666667  0.83333333333333  0.00000000000000  15.999 > 55
   62 O   0.50000000000000  0.83333333333333  0.00000000000000  15.999 > 55
   63 O   0.83333333333333  0.83333333333333  0.00000000000000  15.999 > 55
   64 O   0.16666666666667  0.16666666666667  0.33333333333333  15.999 > 55
   65 O   0.50000000000000  0.16666666666667  0.33333333333333  15.999 > 55
   66 O   0.83333333333333  0.16666666666667  0.33333333333333  15.999 > 55
   67 O   0.16666666666667  0.50000000000000  0.33333333333333  15.999 > 55
   68 O   0.50000000000000  0.50000000000000  0.33333333333333  15.999 > 55
   69 O   0.83333333333333  0.50000000000000  0.33333333333333  15.999 > 55
   70 O   0.16666666666667  0.83333333333333  0.33333333333333  15.999 > 55
   71 O   0.50000000000000  0.83333333333333  0.33333333333333  15.999 > 55
   72 O   0.83333333333333  0.83333333333333  0.33333333333333  15.999 > 55
   73 O   0.16666666666667  0.16666666666667  0.66666666666667  15.999 > 55
   74 O   0.50000000000000  0.16666666666667  0.66666666666667  15.999 > 55
   75 O   0.83333333333333  0.16666666666667  0.66666666666667  15.999 > 55
   76 O   0.16666666666667  0.50000000000000  0.66666666666667  15.999 > 55
   77 O   0.50000000000000  0.50000000000000  0.66666666666667  15.999 > 55
   78 O   0.83333333333333  0.50000000000000  0.66666666666667  15.999 > 55
   79 O   0.16666666666667  0.83333333333333  0.66666666666667  15.999 > 55
   80 O   0.50000000000000  0.83333333333333  0.66666666666667  15.999 > 55
   81 O   0.83333333333333  0.83333333333333  0.66666666666667  15.999 > 55
   82 O   0.16666666666667  0.00000000000000  0.16666666666667  15.999 > 82
   83 O   0.50000000000000  0.00000000000000  0.16666666666667  15.999 > 82
   84 O   0.83333333333333  0.00000000000000  0.16666666666667  15.999 > 82
   85 O   0.16666666666667  0.33333333333333  0.16666666666667  15.999 > 82
   86 O   0.50000000000000  0.33333333333333  0.16666666666667  15.999 > 82
   87 O   0.83333333333333  0.33333333333333  0.16666666666667  15.999 > 82
   88 O   0.16666666666667  0.66666666666667  0.16666666666667  15.999 > 82
   89 O   0.50000000000000  0.66666666666667  0.16666666666667  15.999 > 82
   90 O   0.83333333333333  0.66666666666667  0.16666666666667  15.999 > 82
   91 O   0.16666666666667  0.00000000000000  0.50000000000000  15.999 > 82
   92 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999 > 82
   93 O   0.83333333333333  0.00000000000000  0.50000000000000  15.999 > 82
   94 O   0.16666666666667  0.33333333333333  0.50000000000000  15.999 > 82
   95 O   0.50000000000000  0.33333333333333  0.50000000000000  15.999 > 82
   96 O   0.83333333333333  0.33333333333333  0.50000000000000  15.999 > 82
   97 O   0.16666666666667  0.66666666666667  0.50000000000000  15.999 > 82
   98 O   0.50000000000000  0.66666666666667  0.50000000000000  15.999 > 82
   99 O   0.83333333333333  0.66666666666667  0.50000000000000  15.999 > 82
  100 O   0.16666666666667  0.00000000000000  0.83333333333333  15.999 > 82
  101 O   0.50000000000000  0.00000000000000  0.83333333333333  15.999 > 82
  102 O   0.83333333333333  0.00000000000000  0.83333333333333  15.999 > 82
  103 O   0.16666666666667  0.33333333333333  0.83333333333333  15.999 > 82
  104 O   0.50000000000000  0.33333333333333  0.83333333333333  15.999 > 82
  105 O   0.83333333333333  0.33333333333333  0.83333333333333  15.999 > 82
  106 O   0.16666666666667  0.66666666666667  0.83333333333333  15.999 > 82
  107 O   0.50000000000000  0.66666666666667  0.83333333333333  15.999 > 82
  108 O   0.83333333333333  0.66666666666667  0.83333333333333  15.999 > 82
  109 O   0.00000000000000  0.16666666666667  0.16666666666667  15.999 > 109
  110 O   0.33333333333333  0.16666666666667  0.16666666666667  15.999 > 109
  111 O   0.66666666666667  0.16666666666667  0.16666666666667  15.999 > 109
  112 O   0.00000000000000  0.50000000000000  0.16666666666667  15.999 > 109
  113 O   0.33333333333333  0.50000000000000  0.16666666666667  15.999 > 109
  114 O   0.66666666666667  0.50000000000000  0.16666666666667  15.999 > 109
  115 O   0.00000000000000  0.83333333333333  0.16666666666667  15.999 > 109
  116 O   0.33333333333333  0.83333333333333  0.16666666666667  15.999 > 109
  117 O   0.66666666666667  0.83333333333333  0.16666666666667  15.999 > 109
  118 O   0.00000000000000  0.16666666666667  0.50000000000000  15.999 > 109
  119 O   0.33333333333333  0.16666666666667  0.50000000000000  15.999 > 109
  120 O   0.66666666666667  0.16666666666667  0.50000000000000  15.999 > 109
  121 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999 > 109
  122 O   0.33333333333333  0.50000000000000  0.50000000000000  15.999 > 109
  123 O   0.66666666666667  0.50000000000000  0.50000000000000  15.999 > 109
  124 O   0.00000000000000  0.83333333333333  0.50000000000000  15.999 > 109
  125 O   0.33333333333333  0.83333333333333  0.50000000000000  15.999 > 109
  126 O   0.66666666666667  0.83333333333333  0.50000000000000  15.999 > 109
  127 O   0.00000000000000  0.16666666666667  0.83333333333333  15.999 > 109
  128 O   0.33333333333333  0.16666666666667  0.83333333333333  15.999 > 109
  129 O   0.66666666666667  0.16666666666667  0.83333333333333  15.999 > 109
  130 O   0.00000000000000  0.50000000000000  0.83333333333333  15.999 > 109
  131 O   0.33333333333333  0.50000000000000  0.83333333333333  15.999 > 109
  132 O   0.66666666666667  0.50000000000000  0.83333333333333  15.999 > 109
  133 O   0.00000000000000  0.83333333333333  0.83333333333333  15.999 > 109
  134 O   0.33333333333333  0.83333333333333  0.83333333333333  15.999 > 109
  135 O   0.66666666666667  0.83333333333333  0.83333333333333  15.999 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            6.0396110    0.0000000    0.0000000
            0.0000000    6.0396110    0.0000000
            0.0000000    0.0000000    6.0396110
-------------------------- Born effective charges --------------------------
    1 Pa    6.9426443    0.0000000    0.0000000
            0.0000000    6.9426443    0.0000000
            0.0000000    0.0000000    6.9426443
    2 Tl    1.9946292    0.0000000    0.0000000
            0.0000000    1.9946292    0.0000000
            0.0000000    0.0000000    1.9946292
    3 O    -1.6778173    0.0000000    0.0000000
            0.0000000   -1.6778173    0.0000000
            0.0000000    0.0000000   -5.5816389
    4 O    -1.6778173    0.0000000    0.0000000
            0.0000000   -5.5816389    0.0000000
            0.0000000    0.0000000   -1.6778173
    5 O    -5.5816389    0.0000000    0.0000000
            0.0000000   -1.6778173    0.0000000
            0.0000000    0.0000000   -1.6778173
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 28, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 55, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000000 (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:07:28]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-09 04:07:28]-------------------------
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    4.356066310000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356066310000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356066310000000
Atomic positions (fractional):
    1 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036
    2 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383
    3 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999
    4 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999
    5 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.068198930000001    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.068198930000001    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.068198930000001
Atomic positions (fractional):
    1 Pa  0.16666666666667  0.16666666666667  0.16666666666667 231.036 > 1
    2 Pa  0.50000000000000  0.16666666666667  0.16666666666667 231.036 > 1
    3 Pa  0.83333333333333  0.16666666666667  0.16666666666667 231.036 > 1
    4 Pa  0.16666666666667  0.50000000000000  0.16666666666667 231.036 > 1
    5 Pa  0.50000000000000  0.50000000000000  0.16666666666667 231.036 > 1
    6 Pa  0.83333333333333  0.50000000000000  0.16666666666667 231.036 > 1
    7 Pa  0.16666666666667  0.83333333333333  0.16666666666667 231.036 > 1
    8 Pa  0.50000000000000  0.83333333333333  0.16666666666667 231.036 > 1
    9 Pa  0.83333333333333  0.83333333333333  0.16666666666667 231.036 > 1
   10 Pa  0.16666666666667  0.16666666666667  0.50000000000000 231.036 > 1
   11 Pa  0.50000000000000  0.16666666666667  0.50000000000000 231.036 > 1
   12 Pa  0.83333333333333  0.16666666666667  0.50000000000000 231.036 > 1
   13 Pa  0.16666666666667  0.50000000000000  0.50000000000000 231.036 > 1
   14 Pa  0.50000000000000  0.50000000000000  0.50000000000000 231.036 > 1
   15 Pa  0.83333333333333  0.50000000000000  0.50000000000000 231.036 > 1
   16 Pa  0.16666666666667  0.83333333333333  0.50000000000000 231.036 > 1
   17 Pa  0.50000000000000  0.83333333333333  0.50000000000000 231.036 > 1
   18 Pa  0.83333333333333  0.83333333333333  0.50000000000000 231.036 > 1
   19 Pa  0.16666666666667  0.16666666666667  0.83333333333333 231.036 > 1
   20 Pa  0.50000000000000  0.16666666666667  0.83333333333333 231.036 > 1
   21 Pa  0.83333333333333  0.16666666666667  0.83333333333333 231.036 > 1
   22 Pa  0.16666666666667  0.50000000000000  0.83333333333333 231.036 > 1
   23 Pa  0.50000000000000  0.50000000000000  0.83333333333333 231.036 > 1
   24 Pa  0.83333333333333  0.50000000000000  0.83333333333333 231.036 > 1
   25 Pa  0.16666666666667  0.83333333333333  0.83333333333333 231.036 > 1
   26 Pa  0.50000000000000  0.83333333333333  0.83333333333333 231.036 > 1
   27 Pa  0.83333333333333  0.83333333333333  0.83333333333333 231.036 > 1
   28 Tl  0.00000000000000  0.00000000000000  0.00000000000000 204.383 > 28
   29 Tl  0.33333333333333  0.00000000000000  0.00000000000000 204.383 > 28
   30 Tl  0.66666666666667  0.00000000000000  0.00000000000000 204.383 > 28
   31 Tl  0.00000000000000  0.33333333333333  0.00000000000000 204.383 > 28
   32 Tl  0.33333333333333  0.33333333333333  0.00000000000000 204.383 > 28
   33 Tl  0.66666666666667  0.33333333333333  0.00000000000000 204.383 > 28
   34 Tl  0.00000000000000  0.66666666666667  0.00000000000000 204.383 > 28
   35 Tl  0.33333333333333  0.66666666666667  0.00000000000000 204.383 > 28
   36 Tl  0.66666666666667  0.66666666666667  0.00000000000000 204.383 > 28
   37 Tl  0.00000000000000  0.00000000000000  0.33333333333333 204.383 > 28
   38 Tl  0.33333333333333  0.00000000000000  0.33333333333333 204.383 > 28
   39 Tl  0.66666666666667  0.00000000000000  0.33333333333333 204.383 > 28
   40 Tl  0.00000000000000  0.33333333333333  0.33333333333333 204.383 > 28
   41 Tl  0.33333333333333  0.33333333333333  0.33333333333333 204.383 > 28
   42 Tl  0.66666666666667  0.33333333333333  0.33333333333333 204.383 > 28
   43 Tl  0.00000000000000  0.66666666666667  0.33333333333333 204.383 > 28
   44 Tl  0.33333333333333  0.66666666666667  0.33333333333333 204.383 > 28
   45 Tl  0.66666666666667  0.66666666666667  0.33333333333333 204.383 > 28
   46 Tl  0.00000000000000  0.00000000000000  0.66666666666667 204.383 > 28
   47 Tl  0.33333333333333  0.00000000000000  0.66666666666667 204.383 > 28
   48 Tl  0.66666666666667  0.00000000000000  0.66666666666667 204.383 > 28
   49 Tl  0.00000000000000  0.33333333333333  0.66666666666667 204.383 > 28
   50 Tl  0.33333333333333  0.33333333333333  0.66666666666667 204.383 > 28
   51 Tl  0.66666666666667  0.33333333333333  0.66666666666667 204.383 > 28
   52 Tl  0.00000000000000  0.66666666666667  0.66666666666667 204.383 > 28
   53 Tl  0.33333333333333  0.66666666666667  0.66666666666667 204.383 > 28
   54 Tl  0.66666666666667  0.66666666666667  0.66666666666667 204.383 > 28
   55 O   0.16666666666667  0.16666666666667  0.00000000000000  15.999 > 55
   56 O   0.50000000000000  0.16666666666667  0.00000000000000  15.999 > 55
   57 O   0.83333333333333  0.16666666666667  0.00000000000000  15.999 > 55
   58 O   0.16666666666667  0.50000000000000  0.00000000000000  15.999 > 55
   59 O   0.50000000000000  0.50000000000000  0.00000000000000  15.999 > 55
   60 O   0.83333333333333  0.50000000000000  0.00000000000000  15.999 > 55
   61 O   0.16666666666667  0.83333333333333  0.00000000000000  15.999 > 55
   62 O   0.50000000000000  0.83333333333333  0.00000000000000  15.999 > 55
   63 O   0.83333333333333  0.83333333333333  0.00000000000000  15.999 > 55
   64 O   0.16666666666667  0.16666666666667  0.33333333333333  15.999 > 55
   65 O   0.50000000000000  0.16666666666667  0.33333333333333  15.999 > 55
   66 O   0.83333333333333  0.16666666666667  0.33333333333333  15.999 > 55
   67 O   0.16666666666667  0.50000000000000  0.33333333333333  15.999 > 55
   68 O   0.50000000000000  0.50000000000000  0.33333333333333  15.999 > 55
   69 O   0.83333333333333  0.50000000000000  0.33333333333333  15.999 > 55
   70 O   0.16666666666667  0.83333333333333  0.33333333333333  15.999 > 55
   71 O   0.50000000000000  0.83333333333333  0.33333333333333  15.999 > 55
   72 O   0.83333333333333  0.83333333333333  0.33333333333333  15.999 > 55
   73 O   0.16666666666667  0.16666666666667  0.66666666666667  15.999 > 55
   74 O   0.50000000000000  0.16666666666667  0.66666666666667  15.999 > 55
   75 O   0.83333333333333  0.16666666666667  0.66666666666667  15.999 > 55
   76 O   0.16666666666667  0.50000000000000  0.66666666666667  15.999 > 55
   77 O   0.50000000000000  0.50000000000000  0.66666666666667  15.999 > 55
   78 O   0.83333333333333  0.50000000000000  0.66666666666667  15.999 > 55
   79 O   0.16666666666667  0.83333333333333  0.66666666666667  15.999 > 55
   80 O   0.50000000000000  0.83333333333333  0.66666666666667  15.999 > 55
   81 O   0.83333333333333  0.83333333333333  0.66666666666667  15.999 > 55
   82 O   0.16666666666667  0.00000000000000  0.16666666666667  15.999 > 82
   83 O   0.50000000000000  0.00000000000000  0.16666666666667  15.999 > 82
   84 O   0.83333333333333  0.00000000000000  0.16666666666667  15.999 > 82
   85 O   0.16666666666667  0.33333333333333  0.16666666666667  15.999 > 82
   86 O   0.50000000000000  0.33333333333333  0.16666666666667  15.999 > 82
   87 O   0.83333333333333  0.33333333333333  0.16666666666667  15.999 > 82
   88 O   0.16666666666667  0.66666666666667  0.16666666666667  15.999 > 82
   89 O   0.50000000000000  0.66666666666667  0.16666666666667  15.999 > 82
   90 O   0.83333333333333  0.66666666666667  0.16666666666667  15.999 > 82
   91 O   0.16666666666667  0.00000000000000  0.50000000000000  15.999 > 82
   92 O   0.50000000000000  0.00000000000000  0.50000000000000  15.999 > 82
   93 O   0.83333333333333  0.00000000000000  0.50000000000000  15.999 > 82
   94 O   0.16666666666667  0.33333333333333  0.50000000000000  15.999 > 82
   95 O   0.50000000000000  0.33333333333333  0.50000000000000  15.999 > 82
   96 O   0.83333333333333  0.33333333333333  0.50000000000000  15.999 > 82
   97 O   0.16666666666667  0.66666666666667  0.50000000000000  15.999 > 82
   98 O   0.50000000000000  0.66666666666667  0.50000000000000  15.999 > 82
   99 O   0.83333333333333  0.66666666666667  0.50000000000000  15.999 > 82
  100 O   0.16666666666667  0.00000000000000  0.83333333333333  15.999 > 82
  101 O   0.50000000000000  0.00000000000000  0.83333333333333  15.999 > 82
  102 O   0.83333333333333  0.00000000000000  0.83333333333333  15.999 > 82
  103 O   0.16666666666667  0.33333333333333  0.83333333333333  15.999 > 82
  104 O   0.50000000000000  0.33333333333333  0.83333333333333  15.999 > 82
  105 O   0.83333333333333  0.33333333333333  0.83333333333333  15.999 > 82
  106 O   0.16666666666667  0.66666666666667  0.83333333333333  15.999 > 82
  107 O   0.50000000000000  0.66666666666667  0.83333333333333  15.999 > 82
  108 O   0.83333333333333  0.66666666666667  0.83333333333333  15.999 > 82
  109 O   0.00000000000000  0.16666666666667  0.16666666666667  15.999 > 109
  110 O   0.33333333333333  0.16666666666667  0.16666666666667  15.999 > 109
  111 O   0.66666666666667  0.16666666666667  0.16666666666667  15.999 > 109
  112 O   0.00000000000000  0.50000000000000  0.16666666666667  15.999 > 109
  113 O   0.33333333333333  0.50000000000000  0.16666666666667  15.999 > 109
  114 O   0.66666666666667  0.50000000000000  0.16666666666667  15.999 > 109
  115 O   0.00000000000000  0.83333333333333  0.16666666666667  15.999 > 109
  116 O   0.33333333333333  0.83333333333333  0.16666666666667  15.999 > 109
  117 O   0.66666666666667  0.83333333333333  0.16666666666667  15.999 > 109
  118 O   0.00000000000000  0.16666666666667  0.50000000000000  15.999 > 109
  119 O   0.33333333333333  0.16666666666667  0.50000000000000  15.999 > 109
  120 O   0.66666666666667  0.16666666666667  0.50000000000000  15.999 > 109
  121 O   0.00000000000000  0.50000000000000  0.50000000000000  15.999 > 109
  122 O   0.33333333333333  0.50000000000000  0.50000000000000  15.999 > 109
  123 O   0.66666666666667  0.50000000000000  0.50000000000000  15.999 > 109
  124 O   0.00000000000000  0.83333333333333  0.50000000000000  15.999 > 109
  125 O   0.33333333333333  0.83333333333333  0.50000000000000  15.999 > 109
  126 O   0.66666666666667  0.83333333333333  0.50000000000000  15.999 > 109
  127 O   0.00000000000000  0.16666666666667  0.83333333333333  15.999 > 109
  128 O   0.33333333333333  0.16666666666667  0.83333333333333  15.999 > 109
  129 O   0.66666666666667  0.16666666666667  0.83333333333333  15.999 > 109
  130 O   0.00000000000000  0.50000000000000  0.83333333333333  15.999 > 109
  131 O   0.33333333333333  0.50000000000000  0.83333333333333  15.999 > 109
  132 O   0.66666666666667  0.50000000000000  0.83333333333333  15.999 > 109
  133 O   0.00000000000000  0.83333333333333  0.83333333333333  15.999 > 109
  134 O   0.33333333333333  0.83333333333333  0.83333333333333  15.999 > 109
  135 O   0.66666666666667  0.83333333333333  0.83333333333333  15.999 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            6.0396110    0.0000000    0.0000000
            0.0000000    6.0396110    0.0000000
            0.0000000    0.0000000    6.0396110
-------------------------- Born effective charges --------------------------
    1 Pa    6.9426443    0.0000000    0.0000000
            0.0000000    6.9426443    0.0000000
            0.0000000    0.0000000    6.9426443
    2 Tl    1.9946292    0.0000000    0.0000000
            0.0000000    1.9946292    0.0000000
            0.0000000    0.0000000    1.9946292
    3 O    -1.6778173    0.0000000    0.0000000
            0.0000000   -1.6778173    0.0000000
            0.0000000    0.0000000   -5.5816389
    4 O    -1.6778173    0.0000000    0.0000000
            0.0000000   -5.5816389    0.0000000
            0.0000000    0.0000000   -1.6778173
    5 O    -5.5816389    0.0000000    0.0000000
            0.0000000   -1.6778173    0.0000000
            0.0000000    0.0000000   -1.6778173
----------------------------------------------------------------------------
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 ... [ 11 11 11 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.95, Number of G-points: 305, Lambda: 0.24
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/56) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 56
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   1.454   (   0.000    0.000    0.000)    0.000
   1.454   (   0.000    0.000    0.000)    0.000
   1.454   (   0.000    0.000    0.000)    0.000
   5.358   (   0.000    0.000    0.000)    0.000
   5.358   (   0.000    0.000    0.000)    0.000
   5.358   (   0.000    0.000    0.000)    0.000
   5.528   (   0.000    0.000    0.000)    0.000
   5.528   (   0.000    0.000    0.000)    0.000
   5.528   (   0.000    0.000    0.000)    0.000
  13.225   (   0.000    0.000    0.000)    0.000
  13.225   (   0.000    0.000    0.000)    0.000
  13.225   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/56) =======================
q-point: ( 0.09  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.413   (  18.737    0.000    0.000)   18.737
   0.413   (  18.737    0.000    0.000)   18.737
   0.989   (  35.968    0.000    0.000)   35.968
   1.487   (   3.118    0.000    0.000)    3.118
   1.487   (   3.118    0.000    0.000)    3.118
   2.206   (  24.743    0.000    0.000)   24.743
   5.383   (   2.308    0.000    0.000)    2.308
   5.383   (   2.308    0.000    0.000)    2.308
   5.557   (   2.766    0.000    0.000)    2.766
   5.557   (   2.766    0.000    0.000)    2.766
   5.577   (   4.553    0.000    0.000)    4.553
   7.243   (  19.697    0.000    0.000)   19.697
  13.229   (   0.385    0.000    0.000)    0.385
  13.229   (   0.385    0.000    0.000)    0.385
  19.460   (   5.704    0.000    0.000)    5.704
======================= Grid point 2 (3/56) =======================
q-point: ( 0.18  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.762   (  14.348    0.000    0.000)   14.348
   0.762   (  14.348    0.000    0.000)   14.348
   1.455   (  12.637    0.000    0.000)   12.637
   1.582   (   5.867    0.000    0.000)    5.867
   1.582   (   5.867    0.000    0.000)    5.867
   2.960   (  42.023    0.000    0.000)   42.023
   5.445   (   3.394    0.000    0.000)    3.394
   5.445   (   3.394    0.000    0.000)    3.394
   5.641   (   5.104    0.000    0.000)    5.104
   5.641   (   5.104    0.000    0.000)    5.104
   5.707   (   7.476    0.000    0.000)    7.476
   7.790   (  30.924    0.000    0.000)   30.924
  13.241   (   0.671    0.000    0.000)    0.671
  13.241   (   0.671    0.000    0.000)    0.671
  19.627   (   9.908    0.000    0.000)    9.908
======================= Grid point 3 (4/56) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.005   (   8.922    0.000    0.000)    8.922
   1.005   (   8.922    0.000    0.000)    8.922
   1.644   (   6.792    0.000    0.000)    6.792
   1.721   (   7.080    0.000    0.000)    7.080
   1.721   (   7.080    0.000    0.000)    7.080
   3.792   (  35.749    0.000    0.000)   35.749
   5.518   (   3.378    0.000    0.000)    3.378
   5.518   (   3.378    0.000    0.000)    3.378
   5.758   (   5.816    0.000    0.000)    5.816
   5.758   (   5.816    0.000    0.000)    5.816
   5.872   (   7.893    0.000    0.000)    7.893
   8.452   (  30.687    0.000    0.000)   30.687
  13.256   (   0.764    0.000    0.000)    0.764
  13.256   (   0.764    0.000    0.000)    0.764
  19.853   (  11.207    0.000    0.000)   11.207
======================= Grid point 4 (5/56) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.143   (   4.575    0.000    0.000)    4.575
   1.143   (   4.575    0.000    0.000)    4.575
   1.757   (   4.099    0.000    0.000)    4.099
   1.858   (   5.561    0.000    0.000)    5.561
   1.858   (   5.561    0.000    0.000)    5.561
   4.409   (  22.899    0.000    0.000)   22.899
   5.580   (   2.518    0.000    0.000)    2.518
   5.580   (   2.518    0.000    0.000)    2.518
   5.869   (   4.468    0.000    0.000)    4.468
   5.869   (   4.468    0.000    0.000)    4.468
   6.019   (   5.869    0.000    0.000)    5.869
   9.009   (  21.541    0.000    0.000)   21.541
  13.271   (   0.607    0.000    0.000)    0.607
  13.271   (   0.607    0.000    0.000)    0.607
  20.069   (   8.837    0.000    0.000)    8.837
======================= Grid point 5 (6/56) =======================
q-point: ( 0.45  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.204   (   1.400    0.000    0.000)    1.400
   1.204   (   1.400    0.000    0.000)    1.400
   1.815   (   1.413    0.000    0.000)    1.413
   1.940   (   2.056    0.000    0.000)    2.056
   1.940   (   2.056    0.000    0.000)    2.056
   4.732   (   7.854    0.000    0.000)    7.854
   5.618   (   0.938    0.000    0.000)    0.938
   5.618   (   0.938    0.000    0.000)    0.938
   5.935   (   1.659    0.000    0.000)    1.659
   5.935   (   1.659    0.000    0.000)    1.659
   6.105   (   2.155    0.000    0.000)    2.155
   9.318   (   7.645    0.000    0.000)    7.645
  13.280   (   0.232    0.000    0.000)    0.232
  13.280   (   0.232    0.000    0.000)    0.232
  20.200   (   3.353    0.000    0.000)    3.353
======================= Grid point 12 (7/56) =======================
q-point: ( 0.09  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.572   (  12.333   12.333    0.000)   17.442
   0.764   (  10.901   10.901    0.000)   15.416
   1.224   (  20.286   20.286    0.000)   28.688
   1.522   (   3.328    3.328    0.000)    4.706
   1.810   (  17.017   17.017    0.000)   24.066
   2.175   (  13.323   13.323    0.000)   18.841
   5.301   (  -5.029   -5.029    0.000)    7.112
   5.412   (   2.579    2.579    0.000)    3.647
   5.541   (   2.946    2.946    0.000)    4.166
   5.584   (   2.596    2.596    0.000)    3.672
   5.732   (   8.982    8.982    0.000)   12.703
   7.481   (  19.269   19.269    0.000)   27.250
  13.061   (  -7.867   -7.867    0.000)   11.125
  13.233   (   0.367    0.367    0.000)    0.519
  19.635   (  11.264   11.264    0.000)   15.929
======================= Grid point 13 (8/56) =======================
q-point: ( 0.18  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.834   (  11.658    6.155    0.000)   13.183
   0.925   (   6.753    7.095    0.000)    9.795
   1.511   (   6.704    3.906    0.000)    7.759
   1.623   (   6.175    3.851    0.000)    7.278
   1.971   (   7.878   34.152    0.000)   35.048
   2.913   (  42.719   -1.519    0.000)   42.746
   5.227   (  -1.986  -14.396    0.000)   14.532
   5.481   (   3.757    3.258    0.000)    4.972
   5.640   (   6.012    4.151    0.000)    7.306
   5.664   (   4.913    2.190    0.000)    5.379
   5.935   (  10.138   12.054    0.000)   15.750
   7.976   (  27.008   15.823    0.000)   31.301
  12.966   (  -2.243  -20.696    0.000)   20.817
  13.244   (   0.639    0.317    0.000)    0.714
  19.877   (  12.114   19.918    0.000)   23.312
======================= Grid point 14 (9/56) =======================
q-point: ( 0.27  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.036   (   7.596    2.769    0.000)    8.085
   1.057   (   5.633    2.593    0.000)    6.201
   1.586   (   2.129   -1.338    0.000)    2.514
   1.767   (   7.201    4.239    0.000)    8.357
   2.175   (  10.275   34.331    0.000)   35.836
   3.725   (  33.373   -5.135    0.000)   33.765
   5.219   (   1.214  -19.293    0.000)   19.331
   5.560   (   3.634    3.903    0.000)    5.333
   5.774   (   6.371    4.199    0.000)    7.631
   5.779   (   5.750    1.904    0.000)    6.056
   6.140   (   9.231   14.198    0.000)   16.935
   8.552   (  26.791    8.972    0.000)   28.253
  12.940   (  -0.601  -26.691    0.000)   26.698
  13.259   (   0.728    0.249    0.000)    0.769
  20.136   (  12.343   24.231    0.000)   27.193
======================= Grid point 15 (10/56) =======================
q-point: ( 0.36  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.154   (   3.550    0.537    0.000)    3.590
   1.156   (   4.009    1.138    0.000)    4.167
   1.624   (   1.553   -4.677    0.000)    4.928
   1.904   (   5.504    4.223    0.000)    6.938
   2.364   (   7.351   35.862    0.000)   36.608
   4.271   (  18.663  -13.285    0.000)   22.909
   5.274   (   3.765  -18.121    0.000)   18.508
   5.627   (   2.658    4.294    0.000)    5.050
   5.889   (   4.473    1.846    0.000)    4.839
   5.889   (   4.416    3.208    0.000)    5.458
   6.309   (   6.602   16.027    0.000)   17.334
   9.041   (  19.076    3.186    0.000)   19.341
  12.934   (  -0.128  -29.614    0.000)   29.614
  13.273   (   0.578    0.183    0.000)    0.607
  20.369   (   9.419   26.379    0.000)   28.010
======================= Grid point 16 (11/56) =======================
q-point: ( 0.45  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.203   (   1.182   -0.094    0.000)    1.185
   1.210   (   1.246    0.514    0.000)    1.348
   1.648   (   0.637   -6.484    0.000)    6.516
   1.984   (   2.012    4.088    0.000)    4.556
   2.470   (   2.620   37.275    0.000)   37.367
   4.515   (   5.392  -21.577    0.000)   22.240
   5.348   (   2.419  -13.641    0.000)   13.854
   5.666   (   0.983    4.471    0.000)    4.578
   5.951   (   1.457    1.802    0.000)    2.318
   5.955   (   1.668    1.867    0.000)    2.504
   6.405   (   2.439   17.372    0.000)   17.542
   9.317   (   6.838    0.191    0.000)    6.840
  12.932   (  -0.014  -30.917    0.000)   30.917
  13.281   (   0.221    0.143    0.000)    0.263
  20.508   (   3.530   27.286    0.000)   27.514
======================= Grid point 24 (12/56) =======================
q-point: ( 0.18  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.979   (   6.915    6.915    0.000)    9.779
   1.021   (   3.454    3.454    0.000)    4.885
   1.566   (   0.973    0.973    0.000)    1.376
   1.736   (   6.729    6.729    0.000)    9.517
   2.600   (  17.153   17.153    0.000)   24.258
   3.126   (  26.657   26.657    0.000)   37.699
   4.968   (  -9.593   -9.593    0.000)   13.566
   5.570   (   4.868    4.868    0.000)    6.885
   5.732   (   4.260    4.260    0.000)    6.025
   5.760   (   6.800    6.800    0.000)    9.617
   6.169   (  10.430   10.430    0.000)   14.750
   8.368   (  20.472   20.472    0.000)   28.952
  12.568   ( -15.716  -15.716    0.000)   22.226
  13.253   (   0.552    0.552    0.000)    0.781
  20.318   (  20.870   20.870    0.000)   29.515
======================= Grid point 25 (13/56) =======================
q-point: ( 0.27  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.092   (   3.446    1.255    0.000)    3.667
   1.105   (   4.950    3.423    0.000)    6.018
   1.560   (  -1.034   -1.585    0.000)    1.892
   1.887   (   7.229    6.842    0.000)    9.954
   2.852   (  11.333   22.018    0.000)   24.764
   3.716   (  22.352   10.103    0.000)   24.529
   4.857   (  -0.070  -14.645    0.000)   14.646
   5.669   (   4.454    6.248    0.000)    7.674
   5.836   (   5.428    3.477    0.000)    6.446
   5.902   (   6.218    7.758    0.000)    9.942
   6.373   (   9.158    7.986    0.000)   12.151
   8.791   (  19.326   13.134    0.000)   23.367
  12.318   (  -8.730  -29.751    0.000)   31.005
  13.266   (   0.628    0.433    0.000)    0.763
  20.722   (  17.543   29.317    0.000)   34.164
======================= Grid point 26 (14/56) =======================
q-point: ( 0.36  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.159   (   2.772    0.084    0.000)    2.773
   1.185   (   2.804    1.542    0.000)    3.200
   1.539   (  -0.875   -3.952    0.000)    4.048
   2.020   (   5.248    6.526    0.000)    8.374
   3.100   (  11.879   25.161    0.000)   27.824
   3.955   (   1.955   -8.776    0.000)    8.991
   4.976   (   9.719  -10.654    0.000)   14.420
   5.749   (   3.053    7.065    0.000)    7.696
   5.943   (   4.457    3.203    0.000)    5.488
   5.996   (   2.539    7.083    0.000)    7.524
   6.552   (   7.708    6.769    0.000)   10.258
   9.145   (  13.914    6.749    0.000)   15.464
  12.185   (  -4.358  -38.825    0.000)   39.069
  13.278   (   0.498    0.318    0.000)    0.591
  21.035   (  11.999   34.455    0.000)   36.484
======================= Grid point 27 (15/56) =======================
q-point: ( 0.45  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.201   (   1.077   -0.087    0.000)    1.080
   1.223   (   0.904    0.756    0.000)    1.178
   1.526   (  -0.293   -5.664    0.000)    5.671
   2.096   (   1.881    6.235    0.000)    6.512
   3.307   (   6.390   37.495    0.000)   38.036
   3.916   (  -2.908  -30.528    0.000)   30.666
   5.152   (   5.179   -6.826    0.000)    8.568
   5.793   (   1.099    7.374    0.000)    7.456
   6.009   (   1.693    3.212    0.000)    3.631
   6.018   (   0.073    4.951    0.000)    4.952
   6.672   (   3.215    7.502    0.000)    8.162
   9.348   (   5.085    3.240    0.000)    6.029
  12.126   (  -1.339  -43.079    0.000)   43.100
  13.285   (   0.190    0.247    0.000)    0.311
  21.208   (   4.294   36.721    0.000)   36.971
======================= Grid point 36 (16/56) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.121   (   1.651    1.651    0.000)    2.335
   1.170   (   2.721    2.721    0.000)    3.848
   1.515   (  -2.732   -2.732    0.000)    3.863
   2.034   (   6.864    6.864    0.000)    9.707
   3.024   (   1.245    1.245    0.000)    1.760
   4.136   (  20.274   20.274    0.000)   28.672
   4.634   (  -4.611   -4.611    0.000)    6.521
   5.802   (   5.972    5.972    0.000)    8.446
   5.920   (   4.546    4.546    0.000)    6.429
   6.078   (   8.190    8.190    0.000)   11.583
   6.491   (   4.286    4.286    0.000)    6.062
   9.070   (  12.963   12.963    0.000)   18.333
  11.773   ( -21.789  -21.789    0.000)   30.815
  13.276   (   0.491    0.491    0.000)    0.694
  21.293   (  24.464   24.464    0.000)   34.598
======================= Grid point 37 (17/56) =======================
q-point: ( 0.36  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.162   (   2.113    0.259    0.000)    2.129
   1.216   (   1.663    1.342    0.000)    2.137
   1.455   (  -2.773   -4.062    0.000)    4.918
   2.159   (   4.811    6.351    0.000)    7.968
   3.065   (   3.369  -17.425    0.000)   17.748
   4.240   (  -5.520   24.997    0.000)   25.600
   4.794   (  14.125   -7.324    0.000)   15.911
   5.904   (   3.715    7.388    0.000)    8.269
   6.016   (   4.271    3.698    0.000)    5.650
   6.187   (   1.467   11.076    0.000)   11.173
   6.603   (   6.926   -0.992    0.000)    6.996
   9.310   (   9.595    8.574    0.000)   12.868
  11.415   ( -12.538  -33.283    0.000)   35.566
  13.285   (   0.389    0.359    0.000)    0.529
  21.719   (  15.909   29.676    0.000)   33.672
======================= Grid point 38 (18/56) =======================
q-point: ( 0.45  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.200   (   1.237    0.042    0.000)    1.238
   1.239   (   0.555    0.700    0.000)    0.893
   1.409   (  -1.355   -5.487    0.000)    5.652
   2.228   (   1.702    6.017    0.000)    6.253
   3.141   (   2.721  -31.640    0.000)   31.756
   4.124   (  -3.772   28.450    0.000)   28.699
   5.012   (   5.594   -7.130    0.000)    9.062
   5.956   (   1.248    7.825    0.000)    7.924
   6.082   (   1.719    3.591    0.000)    3.981
   6.175   (  -1.166   10.378    0.000)   10.443
   6.733   (   3.921   -0.681    0.000)    3.979
   9.452   (   3.580    6.267    0.000)    7.218
  11.243   (  -4.056  -39.564    0.000)   39.772
  13.291   (   0.148    0.278    0.000)    0.315
  21.944   (   5.509   32.232    0.000)   32.699
======================= Grid point 48 (19/56) =======================
q-point: ( 0.36  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.174   (   1.028    1.028    0.000)    1.453
   1.239   (   0.866    0.866    0.000)    1.225
   1.372   (  -3.731   -3.731    0.000)    5.276
   2.273   (   4.389    4.389    0.000)    6.206
   2.728   ( -12.714  -12.714    0.000)   17.980
   4.611   (   2.113    2.113    0.000)    2.988
   4.710   (   6.323    6.323    0.000)    8.942
   6.038   (   4.953    4.953    0.000)    7.005
   6.092   (   3.452    3.452    0.000)    4.882
   6.414   (   7.375    7.375    0.000)   10.430
   6.573   (   1.214    1.214    0.000)    1.717
   9.479   (   7.016    7.016    0.000)    9.922
  10.839   ( -21.382  -21.382    0.000)   30.239
  13.292   (   0.284    0.284    0.000)    0.401
  22.235   (  19.261   19.261    0.000)   27.240
======================= Grid point 49 (20/56) =======================
q-point: ( 0.45  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.203   (   1.383    0.218    0.000)    1.400
   1.251   (   0.296    0.466    0.000)    0.552
   1.306   (  -2.295   -4.176    0.000)    4.765
   2.336   (   1.544    4.142    0.000)    4.420
   2.557   (  -3.838  -22.790    0.000)   23.111
   4.556   (  -2.129   12.977    0.000)   13.151
   4.856   (   3.511   -7.363    0.000)    8.157
   6.103   (   1.405    5.835    0.000)    6.002
   6.151   (   1.715    2.811    0.000)    3.292
   6.433   (  -1.175   13.043    0.000)   13.096
   6.689   (   4.856   -2.537    0.000)    5.479
   9.584   (   2.686    5.757    0.000)    6.353
  10.538   (  -7.289  -26.974    0.000)   27.941
  13.296   (   0.108    0.218    0.000)    0.244
  22.507   (   6.617   21.048    0.000)   22.064
======================= Grid point 60 (21/56) =======================
q-point: ( 0.45  0.45  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.211   (   0.645    0.645    0.000)    0.912
   1.240   (  -2.022   -2.022    0.000)    2.859
   1.258   (   0.161    0.161    0.000)    0.227
   2.229   (  -8.119   -8.119    0.000)   11.481
   2.395   (   1.455    1.455    0.000)    2.058
   4.703   (   1.511    1.511    0.000)    2.137
   4.742   (  -1.858   -1.858    0.000)    2.627
   6.187   (   1.919    1.919    0.000)    2.713
   6.193   (   1.289    1.289    0.000)    1.823
   6.643   (   3.044    3.044    0.000)    4.305
   6.654   (   1.864    1.864    0.000)    2.636
   9.671   (   2.254    2.254    0.000)    3.187
  10.151   (  -9.521   -9.521    0.000)   13.465
  13.300   (   0.083    0.083    0.000)    0.117
  22.805   (   7.233    7.233    0.000)   10.229
======================= Grid point 133 (22/56) =======================
q-point: ( 0.09  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.827   (   7.626    7.626    7.626)   13.209
   0.827   (   7.626    7.626    7.626)   13.209
   1.384   (  12.154   12.154   12.154)   21.052
   1.854   (  12.793   12.793   12.793)   22.158
   1.854   (  12.793   12.793   12.793)   22.158
   2.167   (  11.501   11.501   11.501)   19.921
   5.334   (  -1.379   -1.379   -1.379)    2.389
   5.334   (  -1.379   -1.379   -1.379)    2.389
   5.352   (  -5.695   -5.695   -5.695)    9.865
   5.817   (   9.031    9.031    9.031)   15.642
   5.817   (   9.031    9.031    9.031)   15.642
   7.692   (  17.599   17.599   17.599)   30.483
  13.066   (  -5.098   -5.098   -5.098)    8.829
  13.066   (  -5.098   -5.098   -5.098)    8.829
  19.811   (  13.069   13.069   13.069)   22.637
======================= Grid point 134 (23/56) =======================
q-point: ( 0.18  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.961   (   6.712    4.923    4.923)    9.671
   0.962   (   5.258    4.213    4.213)    7.947
   1.531   (   1.777    1.371    1.371)    2.630
   1.973   (   7.033   17.284   17.284)   25.435
   2.063   (   9.610   20.706   20.706)   30.819
   2.899   (  42.091    0.854    0.854)   42.109
   5.189   (  -6.228   -7.925   -7.925)   12.822
   5.345   (   3.750   -3.714   -3.714)    6.454
   5.439   (   5.884   -6.227   -6.227)   10.592
   5.928   (   5.245   12.121   12.121)   17.926
   6.082   (  12.334   10.391   10.391)   19.185
   8.135   (  23.673   13.693   13.693)   30.584
  12.899   (  -5.455  -10.568  -10.568)   15.910
  13.079   (   0.744   -7.862   -7.862)   11.143
  20.098   (  14.106   18.322   18.322)   29.503
======================= Grid point 135 (24/56) =======================
q-point: ( 0.27  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.067   (   4.721    1.702    1.702)    5.299
   1.088   (   5.269    2.374    2.374)    6.248
   1.533   (  -0.384   -2.037   -2.037)    2.906
   2.131   (   7.615   17.292   17.292)   25.613
   2.295   (  11.298   20.833   20.833)   31.555
   3.681   (  31.281   -3.355   -3.355)   31.638
   5.085   (  -3.797  -13.335  -13.335)   19.237
   5.473   (   8.044   -2.078   -2.078)    8.564
   5.572   (   6.488   -6.777   -6.777)   11.573
   6.033   (   4.649   11.748   11.748)   17.252
   6.320   (  10.233   12.329   12.329)   20.216
   8.635   (  23.231    7.511    7.511)   25.544
  12.819   (  -2.585  -16.568  -16.568)   23.573
  13.096   (   0.825   -7.870   -7.870)   11.161
  20.390   (  13.469   21.977   21.977)   33.873
======================= Grid point 136 (25/56) =======================
q-point: ( 0.36  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.152   (   3.225    0.287    0.287)    3.251
   1.177   (   3.174    1.163    1.163)    3.574
   1.531   (   0.150   -4.270   -4.270)    6.041
   2.274   (   5.674   17.342   17.342)   25.173
   2.501   (   7.955   21.809   21.809)   31.852
   4.185   (  17.182   -8.952   -8.952)   21.343
   5.026   (  -1.923  -17.536  -17.536)   24.874
   5.654   (   8.460    2.436    2.436)    9.135
   5.695   (   4.992   -7.262   -7.262)   11.419
   6.115   (   3.075   10.690   10.690)   15.427
   6.499   (   6.720   13.339   13.339)   20.026
   9.061   (  16.707    2.062    2.062)   16.960
  12.780   (  -1.255  -20.039  -20.039)   28.367
  13.111   (   0.640   -7.891   -7.891)   11.178
  20.640   (   9.969   23.906   23.906)   35.248
======================= Grid point 137 (26/56) =======================
q-point: ( 0.45  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.198   (   1.114   -0.192   -0.192)    1.147
   1.221   (   1.051    0.714    0.714)    1.458
   1.535   (   0.173   -5.534   -5.534)    7.828
   2.356   (   2.072   17.453   17.453)   24.769
   2.615   (   2.823   22.728   22.728)   32.267
   4.414   (   5.263  -12.708  -12.708)   18.726
   5.001   (  -0.570  -19.864  -19.864)   28.098
   5.769   (   1.867   -7.525   -7.525)   10.804
   5.789   (   3.652    5.739    5.739)    8.899
   6.158   (   1.025    9.808    9.808)   13.908
   6.594   (   2.336   13.747   13.747)   19.581
   9.304   (   6.044   -0.854   -0.854)    6.163
  12.763   (  -0.389  -21.671  -21.671)   30.650
  13.121   (   0.240   -7.910   -7.910)   11.189
  20.786   (   3.691   24.711   24.711)   35.141
======================= Grid point 145 (27/56) =======================
q-point: ( 0.18  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.042   (   3.066    3.066    2.018)    4.782
   1.046   (   4.312    4.312    3.546)    7.054
   1.534   (  -1.214   -1.214   -1.261)    2.131
   2.143   (   7.907    7.907   31.783)   33.693
   2.635   (  16.764   16.764    3.257)   23.931
   3.144   (  26.386   26.386    3.344)   37.465
   4.953   ( -10.997  -10.997   -2.359)   15.730
   5.385   (   2.849    2.849  -11.873)   12.538
   5.393   (   3.457    3.457  -14.071)   14.896
   6.133   (   7.521    7.521   18.150)   21.036
   6.318   (  11.791   11.791    9.696)   19.289
   8.474   (  17.582   17.582    9.226)   26.522
  12.569   ( -15.915  -15.915    0.123)   22.507
  13.027   (  -0.233   -0.233  -19.239)   19.242
  20.524   (  20.854   20.854   17.827)   34.462
======================= Grid point 146 (28/56) =======================
q-point: ( 0.27  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.103   (   2.900    1.564    1.187)    3.502
   1.134   (   3.852    2.308    1.783)    4.831
   1.495   (  -2.117   -1.973   -3.108)    4.247
   2.312   (   7.823    7.223   32.117)   33.836
   2.890   (  11.602   20.280    3.798)   23.671
   3.707   (  21.095   11.201    0.035)   23.885
   4.774   (  -5.652  -14.160   -8.209)   17.316
   5.487   (   6.299    1.057  -12.790)   14.296
   5.518   (   8.303    2.933  -12.860)   15.586
   6.257   (   4.625    9.546   17.816)   20.735
   6.550   (  10.119    9.146   12.239)   18.326
   8.834   (  16.291   10.754    3.765)   19.880
  12.296   ( -10.241  -28.079   -2.377)   29.983
  13.030   (   0.339   -0.923  -20.437)   20.460
  20.933   (  17.855   27.696   18.658)   37.868
======================= Grid point 147 (29/56) =======================
q-point: ( 0.36  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.160   (   2.417    0.435    0.336)    2.479
   1.201   (   2.478    1.221    1.118)    2.980
   1.457   (  -1.391   -3.136   -4.220)    5.438
   2.455   (   5.606    6.554   32.920)   34.032
   3.143   (  12.073   23.030    4.564)   26.400
   3.945   (   2.665   -6.577   -0.947)    7.159
   4.729   (   0.864  -10.016  -20.665)   22.980
   5.597   (   4.276   -1.231  -15.315)   15.948
   5.732   (  10.480    4.500   -1.981)   11.576
   6.329   (   2.240    9.643   16.367)   19.128
   6.731   (   7.021    7.720   13.288)   16.896
   9.132   (  11.738    4.723   -1.152)   12.705
  12.133   (  -5.566  -36.647   -5.538)   37.479
  13.038   (   0.338   -1.408  -20.832)   20.882
  21.251   (  12.181   32.005   19.343)   39.330
======================= Grid point 148 (30/56) =======================
q-point: ( 0.45  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.197   (   0.975    0.143   -0.077)    0.989
   1.236   (   0.851    0.742    0.933)    1.465
   1.438   (  -0.483   -4.203   -4.794)    6.394
   2.537   (   2.039    6.183   33.708)   34.331
   3.353   (   6.538   35.363    5.266)   36.346
   3.916   (  -2.759  -29.097   -0.267)   29.228
   4.762   (   1.249   -5.268  -30.101)   30.584
   5.660   (   1.561   -2.192  -15.983)   16.208
   5.897   (   4.386    5.054    5.338)    8.560
   6.356   (   0.521    9.017   15.476)   17.919
   6.834   (   2.619    7.673   12.855)   15.198
   9.304   (   4.323    1.334   -3.930)    5.993
  12.058   (  -1.780  -40.847   -7.387)   41.548
  13.043   (   0.135   -1.675  -21.002)   21.069
  21.426   (   4.344   33.922   19.638)   39.437
======================= Grid point 157 (31/56) =======================
q-point: ( 0.27  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.134   (   1.554    1.554    1.217)    2.512
   1.185   (   2.342    2.342    1.055)    3.476
   1.447   (  -2.549   -2.549   -3.595)    5.091
   2.464   (   7.012    7.012   33.175)   34.625
   3.045   (   1.062    1.062    1.950)    2.462
   4.107   (  18.116   18.116   -1.809)   25.683
   4.537   (  -7.121   -7.121   -9.282)   13.696
   5.565   (   5.844    5.844  -16.029)   18.034
   5.626   (   7.367    7.367  -14.125)   17.552
   6.426   (   6.403    6.403   19.384)   21.395
   6.703   (   5.881    5.881   14.309)   16.551
   9.059   (  10.394   10.394   -1.057)   14.737
  11.759   ( -22.287  -22.287   -1.247)   31.544
  13.020   (  -0.210   -0.210  -23.131)   23.133
  21.479   (  23.621   23.621   16.662)   37.330
======================= Grid point 158 (32/56) =======================
q-point: ( 0.36  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.170   (   1.809    0.536    0.830)    2.062
   1.227   (   1.624    1.187    0.882)    2.196
   1.394   (  -2.362   -2.926   -3.377)    5.054
   2.592   (   4.992    6.361   34.205)   35.147
   3.082   (   3.048  -17.750    1.634)   18.083
   4.204   (  -5.004   19.742   -2.951)   20.579
   4.576   (   6.951   -5.220  -19.273)   21.143
   5.665   (   3.544    7.943  -17.405)   19.457
   5.834   (  10.802    5.017   -5.776)   13.237
   6.516   (   2.227    8.296   18.443)   20.346
   6.816   (   5.224    1.035   15.472)   16.363
   9.252   (   7.695    6.415   -5.443)   11.401
  11.385   ( -13.446  -33.170   -2.896)   35.909
  13.018   (  -0.049   -0.667  -24.409)   24.418
  21.892   (  15.481   28.006   15.611)   35.605
======================= Grid point 159 (33/56) =======================
q-point: ( 0.45  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.203   (   1.056    0.363    0.370)    1.176
   1.251   (   0.621    0.627    0.954)    1.300
   1.355   (  -1.198   -3.711   -3.020)    4.932
   2.665   (   1.845    5.981   35.095)   35.649
   3.149   (   2.322  -31.856    0.751)   31.949
   4.102   (  -3.141   23.410   -1.398)   23.661
   4.673   (   2.153   -4.258  -28.095)   28.498
   5.714   (   1.152    8.015  -18.298)   20.010
   6.003   (   4.463    5.037    2.458)    7.164
   6.534   (   0.083    8.154   18.141)   19.890
   6.906   (   2.629   -0.117   14.067)   14.312
   9.366   (   2.898    4.338   -7.972)    9.528
  11.199   (  -4.454  -39.251   -4.337)   39.740
  13.017   (  -0.006   -0.932  -25.029)   25.046
  22.112   (   5.372   30.169   15.127)   34.174
======================= Grid point 169 (34/56) =======================
q-point: ( 0.36  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.185   (   1.016    1.016    1.054)    1.782
   1.249   (   0.849    0.849    0.803)    1.445
   1.336   (  -2.631   -2.631   -2.027)    4.237
   2.710   (   4.653    4.653   35.673)   36.275
   2.737   ( -13.136  -13.136    0.886)   18.598
   4.417   (   0.234    0.234  -17.656)   17.660
   4.549   (   2.912    2.912  -14.097)   14.687
   5.874   (   8.322    8.322   -9.710)   15.258
   5.946   (   7.053    7.053   -7.968)   12.767
   6.663   (   4.910    4.910   17.053)   18.412
   6.817   (   0.598    0.598   17.075)   17.096
   9.381   (   5.470    5.470   -9.191)   12.013
  10.801   ( -21.970  -21.970   -3.481)   31.265
  13.007   (  -0.388   -0.388  -26.334)   26.340
  22.381   (  18.318   18.318   13.227)   29.087
======================= Grid point 170 (35/56) =======================
q-point: ( 0.45  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.211   (   1.188    0.428    0.824)    1.508
   1.261   (   0.405    0.343    0.845)    0.998
   1.287   (  -1.809   -2.590   -0.827)    3.266
   2.558   (  -4.099  -23.163    0.115)   23.524
   2.778   (   1.692    4.533   36.805)   37.122
   4.399   (  -0.523    6.323  -14.297)   15.642
   4.577   (  -0.107   -4.693  -23.729)   24.189
   5.956   (   1.464   12.647  -10.363)   16.416
   6.098   (   4.504    3.802   -0.663)    5.932
   6.698   (  -0.208    7.386   17.067)   18.598
   6.869   (   2.610   -2.459   14.548)   14.984
   9.463   (   2.136    4.450  -11.284)   12.317
  10.490   (  -7.608  -27.411   -4.497)   28.800
  13.002   (  -0.132   -0.588  -27.311)   27.317
  22.641   (   6.311   19.835   12.077)   24.065
======================= Grid point 181 (36/56) =======================
q-point: ( 0.45  0.45  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.222   (   0.628    0.628    0.996)    1.335
   1.243   (  -1.303   -1.303    0.488)    1.906
   1.269   (   0.110    0.110    1.030)    1.041
   2.223   (  -8.367   -8.367   -0.593)   11.848
   2.845   (   1.724    1.724   38.433)   38.511
   4.455   (   0.854    0.854  -21.899)   21.932
   4.495   (  -2.623   -2.623  -21.648)   21.963
   6.147   (   3.659    3.659   -2.559)    5.773
   6.157   (   2.690    2.690   -2.395)    4.496
   6.817   (   2.084    2.084   13.890)   14.199
   6.833   (   0.427    0.427   14.233)   14.245
   9.532   (   1.797    1.797  -13.053)   13.298
  10.095   (  -9.797   -9.797   -5.302)   14.835
  12.993   (  -0.207   -0.207  -28.490)   28.492
  22.921   (   6.836    6.836   10.569)   14.323
======================= Grid point 266 (37/56) =======================
q-point: ( 0.18  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.099   (   2.548    2.548    2.548)    4.413
   1.099   (   2.548    2.548    2.548)    4.413
   1.499   (  -2.413   -2.413   -2.413)    4.180
   2.732   (  12.458   12.458   12.458)   21.578
   2.732   (  12.458   12.458   12.458)   21.578
   3.389   (  21.583   21.583   21.583)   37.383
   4.806   ( -11.748  -11.748  -11.748)   20.348
   5.231   (  -1.197   -1.197   -1.197)    2.073
   5.231   (  -1.197   -1.197   -1.197)    2.073
   6.490   (  11.144   11.144   11.144)   19.302
   6.490   (  11.144   11.144   11.144)   19.302
   8.702   (  11.541   11.541   11.541)   19.989
  12.573   ( -10.895  -10.895  -10.895)   18.870
  12.573   ( -10.895  -10.895  -10.895)   18.870
  20.965   (  22.488   22.488   22.488)   38.950
======================= Grid point 267 (38/56) =======================
q-point: ( 0.27  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.137   (   1.822    1.835    1.835)    3.170
   1.167   (   3.166    1.501    1.501)    3.812
   1.440   (  -2.871   -2.466   -2.466)    4.517
   2.871   (   7.305   14.793   14.793)   22.160
   3.028   (  12.900    8.464    8.464)   17.598
   3.826   (  14.863   14.911   14.911)   25.799
   4.535   ( -12.954  -13.353  -13.353)   22.900
   5.264   (   3.128   -5.693   -5.693)    8.638
   5.374   (  14.155   -2.392   -2.392)   14.554
   6.593   (   3.487   13.765   13.765)   19.777
   6.781   (  12.006    9.510    9.510)   18.028
   8.929   (   9.882    4.871    4.871)   12.046
  12.162   ( -16.791  -13.682  -13.682)   25.619
  12.577   (   0.225  -17.144  -17.144)   24.246
  21.411   (  19.558   25.107   25.107)   40.538
======================= Grid point 268 (39/56) =======================
q-point: ( 0.36  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.174   (   1.700    0.955    0.955)    2.171
   1.224   (   2.213    1.083    1.083)    2.691
   1.389   (  -1.940   -2.604   -2.604)    4.163
   3.024   (   6.924   14.953   14.953)   22.251
   3.316   (  14.540    7.275    7.275)   17.812
   3.967   (  -0.126    7.934    7.934)   11.221
   4.291   (  -9.989  -17.079  -17.079)   26.137
   5.332   (   3.041   -8.174   -8.174)   11.952
   5.726   (  16.680    0.833    0.833)   16.722
   6.640   (   1.151   13.213   13.213)   18.721
   6.985   (   7.455    9.745    9.745)   15.669
   9.106   (   6.934   -1.027   -1.027)    7.084
  11.881   ( -10.157  -21.132  -21.132)   31.564
  12.582   (   0.167  -17.785  -17.785)   25.153
  21.757   (  13.124   27.032   27.032)   40.419
======================= Grid point 269 (40/56) =======================
q-point: ( 0.45  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.202   (   0.800    0.572    0.572)    1.138
   1.256   (   0.780    0.907    0.907)    1.501
   1.360   (  -0.741   -2.857   -2.857)    4.108
   3.136   (   3.185   16.522   16.522)   23.582
   3.609   (  12.416   13.485   13.485)   22.757
   3.873   (  -7.649   -2.421   -2.421)    8.380
   4.138   (  -4.098  -20.856  -20.856)   29.778
   5.378   (   1.207   -9.373   -9.373)   13.311
   5.982   (   6.667    3.194    3.194)    8.053
   6.650   (   0.078   12.604   12.604)   17.825
   7.090   (   2.521    9.786    9.786)   14.068
   9.208   (   2.596   -4.461   -4.461)    6.822
  11.739   (  -3.400  -25.295  -25.295)   35.934
  12.584   (   0.060  -18.167  -18.167)   25.692
  21.945   (   4.620   27.915   27.915)   39.747
======================= Grid point 278 (41/56) =======================
q-point: ( 0.27  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.170   (   1.317    1.317    2.154)    2.848
   1.205   (   2.111    2.111    0.832)    3.099
   1.388   (  -2.434   -2.434   -2.277)    4.127
   3.075   (   8.023    8.023   17.348)   20.729
   3.100   (   0.437    0.437    3.071)    3.133
   4.192   (  12.761   12.761   14.949)   23.434
   4.261   ( -12.226  -12.226  -16.244)   23.723
   5.298   (   5.884    5.884   -7.830)   11.426
   5.407   (   8.629    8.629   -7.041)   14.089
   6.805   (   6.263    6.263   16.588)   18.804
   6.963   (   8.376    8.376    9.834)   15.396
   9.024   (   4.086    4.086   -2.332)    6.231
  11.723   ( -23.629  -23.629   -2.100)   33.483
  12.397   (  -3.533   -3.533  -34.264)   34.626
  21.916   (  22.201   22.201   23.435)   39.178
======================= Grid point 279 (42/56) =======================
q-point: ( 0.36  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.196   (   1.228    1.039    1.629)    2.289
   1.247   (   1.730    1.087    0.962)    2.258
   1.342   (  -1.878   -1.879   -1.729)    3.170
   3.112   (   1.847  -15.292    3.238)   15.740
   3.249   (   7.540    4.813   16.104)   18.422
   4.046   (  -7.596   -6.788  -22.213)   24.438
   4.289   (  -1.521   10.994   11.139)   15.724
   5.359   (   1.618   12.780  -10.935)   16.897
   5.754   (  18.531    1.692   -2.231)   18.741
   6.875   (   1.140    9.214   15.470)   18.042
   7.120   (   6.168    3.379   12.377)   14.235
   9.098   (   2.930    0.257   -8.838)    9.315
  11.292   ( -16.577  -31.281   -6.262)   35.951
  12.352   (  -1.320   -6.230  -36.747)   37.295
  22.307   (  14.721   24.364   22.504)   36.287
======================= Grid point 280 (43/56) =======================
q-point: ( 0.45  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.219   (   0.759    0.929    1.162)    1.670
   1.273   (   0.694    0.680    1.072)    1.447
   1.312   (  -0.948   -1.845   -1.225)    2.409
   3.162   (   2.069  -30.922    1.093)   31.011
   3.355   (   2.198    6.347   19.766)   20.876
   3.968   (  -0.208   -2.329  -33.121)   33.203
   4.211   (  -4.012   11.241   14.671)   18.913
   5.383   (   0.593   12.420  -12.709)   17.780
   6.040   (   7.430    2.425    1.162)    7.902
   6.879   (  -0.183    9.380   14.706)   17.444
   7.209   (   2.243    1.787   13.228)   13.535
   9.142   (   1.168   -1.651  -12.655)   12.815
  11.056   (  -5.807  -36.295   -9.836)   38.050
  12.335   (  -0.401   -7.398  -37.515)   38.240
  22.516   (   5.120   25.454   22.048)   34.062
======================= Grid point 290 (44/56) =======================
q-point: ( 0.36  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.217   (   0.971    0.971    1.881)    2.329
   1.268   (   0.905    0.905    1.040)    1.649
   1.309   (  -1.320   -1.320   -0.601)    1.961
   2.762   ( -14.250  -14.250    1.390)   20.200
   3.397   (   7.305    7.305   16.024)   19.066
   3.938   (  -3.409   -3.409  -26.590)   27.024
   4.331   (  -4.191   -4.191    5.657)    8.193
   5.725   (  11.920   11.920   -4.649)   17.487
   5.819   (   9.532    9.532   -4.431)   14.190
   7.016   (   3.880    3.880   15.849)   16.772
   7.156   (   1.261    1.261   14.371)   14.481
   9.114   (   1.130    1.130  -15.695)   15.777
  10.701   ( -23.572  -23.572   -5.918)   33.857
  12.263   (  -2.834   -2.834  -42.832)   43.019
  22.736   (  16.173   16.173   19.386)   29.982
======================= Grid point 291 (45/56) =======================
q-point: ( 0.45  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.237   (   0.854    0.767    1.635)    1.998
   1.277   (  -0.297   -0.679    0.645)    0.982
   1.290   (  -0.257    0.189    0.959)    1.010
   2.561   (  -4.831  -24.182    0.084)   24.660
   3.516   (   3.209    8.443   16.651)   18.943
   3.916   (   0.964   -2.590  -32.456)   32.574
   4.213   (  -5.054   -7.117    6.351)   10.795
   5.788   (   0.622   20.752   -5.887)   21.579
   6.083   (   7.759    1.655   -0.780)    7.972
   7.042   (  -0.181    6.185   15.255)   16.462
   7.193   (   1.487   -2.374   14.969)   15.228
   9.135   (   0.620    0.578  -19.425)   19.443
  10.358   (  -8.564  -28.327   -7.874)   30.622
  12.225   (  -0.901   -3.626  -44.990)   45.144
  22.966   (   5.604   16.963   17.820)   25.233
======================= Grid point 302 (46/56) =======================
q-point: ( 0.45  0.45  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.251   (   0.578    0.578    1.742)    1.924
   1.263   (  -0.534   -0.534    1.381)    1.574
   1.297   (   0.204    0.204    1.454)    1.483
   2.205   (  -9.012   -9.012   -1.039)   12.787
   3.692   (   7.152    7.152   16.707)   19.530
   3.878   (  -0.332   -0.332  -31.667)   31.671
   4.029   (  -8.557   -8.557    0.926)   12.137
   6.100   (   4.915    4.915   -1.951)    7.219
   6.112   (   3.825    3.825   -1.984)    5.761
   7.129   (   1.458    1.458   14.785)   14.928
   7.151   (  -0.534   -0.534   14.874)   14.893
   9.149   (   0.508    0.508  -22.780)   22.792
   9.941   ( -10.552  -10.552   -9.062)   17.459
  12.175   (  -1.179   -1.179  -47.848)   47.877
  23.207   (   5.879    5.879   15.659)   17.730
======================= Grid point 399 (47/56) =======================
q-point: ( 0.27  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.220   (   1.532    1.532    1.532)    2.653
   1.220   (   1.532    1.532    1.532)    2.653
   1.346   (  -1.726   -1.726   -1.726)    2.990
   3.165   (   0.180    0.180    0.180)    0.312
   3.165   (   0.180    0.180    0.180)    0.312
   3.909   ( -16.353  -16.353  -16.353)   28.324
   4.636   (  17.687   17.687   17.687)   30.635
   5.309   (   4.390    4.390    4.390)    7.603
   5.309   (   4.390    4.390    4.390)    7.603
   7.107   (   8.190    8.190    8.190)   14.185
   7.107   (   8.190    8.190    8.190)   14.185
   8.966   (  -3.022   -3.022   -3.022)    5.234
  11.676   ( -17.697  -17.697  -17.697)   30.652
  11.676   ( -17.697  -17.697  -17.697)   30.652
  22.395   (  21.261   21.261   21.261)   36.826
======================= Grid point 400 (48/56) =======================
q-point: ( 0.36  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.233   (   0.721    1.731    1.731)    2.552
   1.266   (   1.871    0.934    0.934)    2.290
   1.316   (  -1.178   -0.837   -0.837)    1.670
   3.123   (  -0.391   -8.104   -8.104)   11.467
   3.217   (   1.802   -6.242   -6.242)    9.009
   3.615   ( -11.159  -16.313  -16.313)   25.627
   4.784   (   0.438   22.030   22.030)   31.158
   5.268   (  -1.186    6.758    6.758)    9.630
   5.725   (  21.770   -0.751   -0.751)   21.796
   7.150   (   0.342   10.733   10.733)   15.182
   7.312   (   7.363    6.043    6.043)   11.280
   8.909   (  -2.204   -8.577   -8.577)   12.328
  11.075   ( -23.333  -17.354  -17.354)   33.864
  11.633   (  -1.757  -27.040  -27.040)   38.281
  22.770   (  14.200   20.748   20.748)   32.597
======================= Grid point 401 (49/56) =======================
q-point: ( 0.45  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.248   (   0.577    1.525    1.525)    2.232
   1.293   (   0.681    0.881    0.881)    1.420
   1.297   (  -0.629   -0.295   -0.295)    0.755
   3.138   (   1.009  -14.519  -14.519)   20.558
   3.237   (   0.288  -17.152  -17.152)   24.258
   3.460   (  -3.597   -8.423   -8.423)   12.444
   4.766   (  -0.984   23.100   23.100)   32.683
   5.253   (  -0.321    5.813    5.813)    8.227
   6.057   (   8.549    0.585    0.585)    8.589
   7.143   (  -0.489   10.540   10.540)   14.914
   7.414   (   2.448    5.941    5.941)    8.751
   8.879   (  -0.674  -11.486  -11.486)   16.257
  10.738   (  -8.371  -22.218  -22.218)   32.518
  11.607   (  -0.657  -28.020  -28.020)   39.632
  22.972   (   4.939   20.494   20.494)   29.400
======================= Grid point 411 (50/56) =======================
q-point: ( 0.36  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.260   (   0.892    0.892    2.201)    2.537
   1.289   (   1.184    1.184    0.882)    1.892
   1.305   (  -0.293   -0.293    0.121)    0.432
   2.791   ( -15.596  -15.596    1.318)   22.095
   3.211   (   0.280    0.280  -20.805)   20.808
   3.366   (  -7.609   -7.609  -26.986)   29.053
   4.920   (  -1.361   -1.361   37.418)   37.467
   5.661   (  12.906   12.906   -1.840)   18.344
   5.747   (  10.582   10.582   -2.678)   15.204
   7.306   (   3.478    3.478   11.774)   12.760
   7.396   (   3.325    3.325    8.390)    9.618
   8.763   (  -5.132   -5.132  -16.918)   18.409
  10.567   ( -25.730  -25.730   -6.536)   36.970
  11.313   (  -7.288   -7.288  -45.903)   47.046
  23.138   (  13.928   13.928   18.096)   26.748
======================= Grid point 412 (51/56) =======================
q-point: ( 0.45  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.276   (   0.608    1.019    1.980)    2.308
   1.295   (  -0.345    0.123    1.044)    1.106
   1.312   (   0.494    0.839    0.967)    1.372
   2.561   (  -5.809  -25.397   -0.089)   26.053
   3.192   (  -1.232   -0.235  -32.743)   32.767
   3.292   (  -0.721   -7.097  -23.620)   24.674
   4.885   (  -1.198   -2.344   40.041)   40.128
   5.696   (  -0.080   24.858   -3.229)   25.067
   6.066   (   9.050    0.251   -0.837)    9.093
   7.319   (  -0.523    6.365   11.171)   12.867
   7.460   (   1.873   -0.928   10.136)   10.350
   8.696   (  -1.465   -5.820  -21.390)   22.216
  10.171   ( -10.309  -27.826   -9.884)   31.277
  11.225   (  -1.894  -10.537  -48.081)   49.258
  23.336   (   4.845   13.802   16.718)   22.214
======================= Grid point 423 (52/56) =======================
q-point: ( 0.45  0.45  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.291   (   0.496    0.496    2.012)    2.131
   1.297   (  -0.039   -0.039    1.750)    1.751
   1.326   (   0.397    0.397    1.238)    1.360
   2.181   (  -9.798   -9.798   -1.198)   13.908
   3.185   (  -0.980   -0.980  -32.284)   32.314
   3.191   (  -1.768   -1.768  -32.965)   33.059
   4.838   (  -1.359   -1.359   42.926)   42.969
   6.063   (   5.212    5.212   -1.597)    7.542
   6.073   (   4.269    4.269   -1.696)    6.271
   7.405   (   1.232    1.232   11.286)   11.420
   7.426   (  -0.624   -0.624   11.177)   11.212
   8.626   (  -1.371   -1.371  -26.125)   26.197
   9.736   ( -11.570  -11.570  -10.079)   19.218
  11.093   (  -2.869   -2.869  -53.269)   53.423
  23.532   (   4.804    4.804   14.729)   16.220
======================= Grid point 532 (53/56) =======================
q-point: ( 0.36  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.303   (   1.139    1.139    1.139)    1.973
   1.303   (   1.139    1.139    1.139)    1.973
   1.313   (   0.444    0.444    0.444)    0.769
   2.801   ( -10.223  -10.223  -10.223)   17.707
   2.801   ( -10.223  -10.223  -10.223)   17.707
   2.872   ( -15.557  -15.557  -15.557)   26.946
   5.549   (  11.148   11.148   11.148)   19.309
   5.703   (   6.789    6.789    6.789)   11.760
   5.703   (   6.789    6.789    6.789)   11.760
   7.503   (   4.537    4.537    4.537)    7.859
   7.503   (   4.537    4.537    4.537)    7.859
   8.456   ( -11.452  -11.452  -11.452)   19.836
  10.443   ( -20.208  -20.208  -20.208)   35.001
  10.443   ( -20.208  -20.208  -20.208)   35.001
  23.460   (  12.218   12.218   12.218)   21.163
======================= Grid point 533 (54/56) =======================
q-point: ( 0.45  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.313   (   0.414    1.358    1.358)    1.965
   1.318   (   0.011    1.033    1.033)    1.460
   1.329   (   0.658    0.806    0.806)    1.317
   2.545   (  -7.792  -13.979  -13.979)   21.250
   2.569   (  -5.626  -12.515  -12.515)   18.572
   2.835   (   0.247  -18.859  -18.859)   26.672
   5.581   (  -0.324   15.339   15.339)   21.695
   5.680   (  -0.563   11.140   11.140)   15.764
   6.049   (   9.610   -0.703   -0.703)    9.662
   7.506   (  -0.713    6.664    6.664)    9.452
   7.601   (   2.345    3.176    3.176)    5.067
   8.304   (  -3.344  -14.639  -14.639)   20.971
   9.919   ( -13.358  -16.478  -16.478)   26.861
  10.368   (  -1.907  -29.030  -29.030)   41.099
  23.634   (   4.258   11.312   11.312)   16.554
======================= Grid point 544 (55/56) =======================
q-point: ( 0.45  0.45  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.330   (   0.406    0.406    1.634)    1.731
   1.332   (   0.292    0.292    1.456)    1.514
   1.346   (   0.545    0.545    0.685)    1.032
   2.158   ( -10.438  -10.438   -0.980)   14.795
   2.553   (  -3.564   -3.564  -26.813)   27.283
   2.556   (  -3.725   -3.725  -26.388)   26.909
   5.603   (  -0.570   -0.570   29.131)   29.142
   6.034   (   5.212    5.212   -1.153)    7.461
   6.042   (   4.441    4.441   -1.219)    6.398
   7.594   (   1.107    1.107    6.799)    6.977
   7.611   (  -0.128   -0.128    6.472)    6.474
   8.119   (  -3.567   -3.567  -20.905)   21.505
   9.543   ( -12.531  -12.531   -7.863)   19.388
  10.065   (  -4.886   -4.886  -42.548)   43.105
  23.795   (   3.944    3.944    9.973)   11.426
======================= Grid point 665 (56/56) =======================
q-point: ( 0.45  0.45  0.45)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.00e-05 3.36e-05 3.36e-05 3.36e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.354   (   0.508    0.508    0.508)    0.881
   1.355   (   0.444    0.444    0.444)    0.768
   1.355   (   0.444    0.444    0.444)    0.768
   2.132   (  -6.761   -6.761   -6.761)   11.711
   2.152   (  -7.680   -7.680   -7.680)   13.302
   2.152   (  -7.680   -7.680   -7.680)   13.302
   6.012   (   3.653    3.653    3.653)    6.327
   6.024   (   2.853    2.853    2.853)    4.942
   6.024   (   2.853    2.853    2.853)    4.942
   7.688   (   1.438    1.438    1.438)    2.491
   7.688   (   1.438    1.438    1.438)    2.491
   7.817   (  -6.590   -6.590   -6.590)   11.415
   9.426   (  -9.731   -9.731   -9.731)   16.855
   9.426   (  -9.731   -9.731   -9.731)   16.855
  23.937   (   3.477    3.477    3.477)    6.021
=================== 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/19965
   10.0    130.367    130.367    130.367      0.000     -0.000      0.000 3/19965
   20.0     32.020     32.020     32.020      0.000     -0.000      0.000 3/19965
   30.0     22.136     22.136     22.136      0.000     -0.000      0.000 3/19965
   40.0     18.979     18.979     18.979      0.000     -0.000      0.000 3/19965
   50.0     16.875     16.875     16.875      0.000     -0.000      0.000 3/19965
   60.0     15.149     15.149     15.149      0.000     -0.000      0.000 3/19965
   70.0     13.693     13.693     13.693      0.000     -0.000      0.000 3/19965
   80.0     12.461     12.461     12.461      0.000     -0.000      0.000 3/19965
   90.0     11.419     11.419     11.419      0.000     -0.000      0.000 3/19965
  100.0     10.533     10.533     10.533      0.000     -0.000      0.000 3/19965
  110.0      9.775      9.775      9.775      0.000     -0.000      0.000 3/19965
  120.0      9.121      9.121      9.121      0.000     -0.000      0.000 3/19965
  130.0      8.553      8.553      8.553      0.000     -0.000      0.000 3/19965
  140.0      8.055      8.055      8.055      0.000     -0.000      0.000 3/19965
  150.0      7.615      7.615      7.615      0.000     -0.000      0.000 3/19965
  160.0      7.224      7.224      7.224      0.000     -0.000      0.000 3/19965
  170.0      6.873      6.873      6.873      0.000     -0.000      0.000 3/19965
  180.0      6.557      6.557      6.557      0.000     -0.000      0.000 3/19965
  190.0      6.271      6.271      6.271      0.000     -0.000      0.000 3/19965
  200.0      6.010      6.010      6.010      0.000     -0.000      0.000 3/19965
  210.0      5.771      5.771      5.771      0.000     -0.000      0.000 3/19965
  220.0      5.551      5.551      5.551      0.000     -0.000      0.000 3/19965
  230.0      5.349      5.349      5.349      0.000     -0.000      0.000 3/19965
  240.0      5.161      5.161      5.161      0.000     -0.000      0.000 3/19965
  250.0      4.986      4.986      4.986      0.000     -0.000      0.000 3/19965
  260.0      4.823      4.823      4.823      0.000     -0.000      0.000 3/19965
  270.0      4.671      4.671      4.671      0.000     -0.000      0.000 3/19965
  280.0      4.528      4.528      4.528      0.000     -0.000      0.000 3/19965
  290.0      4.394      4.394      4.394      0.000     -0.000      0.000 3/19965
  300.0      4.268      4.268      4.268      0.000     -0.000      0.000 3/19965
  310.0      4.149      4.149      4.149      0.000     -0.000      0.000 3/19965
  320.0      4.036      4.036      4.036      0.000     -0.000      0.000 3/19965
  330.0      3.929      3.929      3.929      0.000     -0.000      0.000 3/19965
  340.0      3.828      3.828      3.828      0.000     -0.000      0.000 3/19965
  350.0      3.731      3.731      3.731      0.000     -0.000      0.000 3/19965
  360.0      3.640      3.640      3.640      0.000     -0.000      0.000 3/19965
  370.0      3.553      3.553      3.553      0.000     -0.000      0.000 3/19965
  380.0      3.469      3.469      3.469      0.000     -0.000      0.000 3/19965
  390.0      3.390      3.390      3.390      0.000     -0.000      0.000 3/19965
  400.0      3.314      3.314      3.314      0.000     -0.000      0.000 3/19965
  410.0      3.241      3.241      3.241      0.000     -0.000      0.000 3/19965
  420.0      3.172      3.172      3.172      0.000     -0.000      0.000 3/19965
  430.0      3.105      3.105      3.105      0.000     -0.000      0.000 3/19965
  440.0      3.041      3.041      3.041      0.000     -0.000      0.000 3/19965
  450.0      2.979      2.979      2.979      0.000     -0.000      0.000 3/19965
  460.0      2.920      2.920      2.920      0.000     -0.000      0.000 3/19965
  470.0      2.863      2.863      2.863      0.000     -0.000      0.000 3/19965
  480.0      2.809      2.809      2.809      0.000     -0.000      0.000 3/19965
  490.0      2.756      2.756      2.756      0.000     -0.000      0.000 3/19965
  500.0      2.705      2.705      2.705      0.000     -0.000      0.000 3/19965
  510.0      2.656      2.656      2.656      0.000     -0.000      0.000 3/19965
  520.0      2.609      2.609      2.609      0.000     -0.000      0.000 3/19965
  530.0      2.563      2.563      2.563      0.000     -0.000      0.000 3/19965
  540.0      2.519      2.519      2.519      0.000     -0.000      0.000 3/19965
  550.0      2.476      2.476      2.476      0.000     -0.000      0.000 3/19965
  560.0      2.435      2.435      2.435      0.000     -0.000      0.000 3/19965
  570.0      2.395      2.395      2.395      0.000     -0.000      0.000 3/19965
  580.0      2.356      2.356      2.356      0.000     -0.000      0.000 3/19965
  590.0      2.319      2.319      2.319      0.000     -0.000      0.000 3/19965
  600.0      2.282      2.282      2.282      0.000     -0.000      0.000 3/19965
  610.0      2.247      2.247      2.247      0.000     -0.000      0.000 3/19965
  620.0      2.213      2.213      2.213      0.000     -0.000      0.000 3/19965
  630.0      2.180      2.180      2.180      0.000     -0.000      0.000 3/19965
  640.0      2.147      2.147      2.147      0.000     -0.000      0.000 3/19965
  650.0      2.116      2.116      2.116      0.000     -0.000      0.000 3/19965
  660.0      2.086      2.086      2.086      0.000     -0.000      0.000 3/19965
  670.0      2.056      2.056      2.056      0.000     -0.000      0.000 3/19965
  680.0      2.027      2.027      2.027      0.000     -0.000      0.000 3/19965
  690.0      1.999      1.999      1.999      0.000     -0.000      0.000 3/19965
  700.0      1.972      1.972      1.972      0.000     -0.000      0.000 3/19965
  710.0      1.945      1.945      1.945      0.000     -0.000      0.000 3/19965
  720.0      1.920      1.920      1.920      0.000     -0.000      0.000 3/19965
  730.0      1.894      1.894      1.894      0.000     -0.000      0.000 3/19965
  740.0      1.870      1.870      1.870      0.000     -0.000      0.000 3/19965
  750.0      1.846      1.846      1.846      0.000     -0.000      0.000 3/19965
  760.0      1.823      1.823      1.823      0.000     -0.000      0.000 3/19965
  770.0      1.800      1.800      1.800      0.000     -0.000      0.000 3/19965
  780.0      1.778      1.778      1.778      0.000     -0.000      0.000 3/19965
  790.0      1.756      1.756      1.756      0.000     -0.000      0.000 3/19965
  800.0      1.735      1.735      1.735      0.000     -0.000      0.000 3/19965
  810.0      1.714      1.714      1.714      0.000     -0.000      0.000 3/19965
  820.0      1.694      1.694      1.694      0.000     -0.000      0.000 3/19965
  830.0      1.674      1.674      1.674      0.000     -0.000      0.000 3/19965
  840.0      1.655      1.655      1.655      0.000     -0.000      0.000 3/19965
  850.0      1.636      1.636      1.636      0.000     -0.000      0.000 3/19965
  860.0      1.618      1.618      1.618      0.000     -0.000      0.000 3/19965
  870.0      1.600      1.600      1.600      0.000     -0.000      0.000 3/19965
  880.0      1.582      1.582      1.582      0.000     -0.000      0.000 3/19965
  890.0      1.565      1.565      1.565      0.000     -0.000      0.000 3/19965
  900.0      1.548      1.548      1.548      0.000     -0.000      0.000 3/19965
  910.0      1.532      1.532      1.532      0.000     -0.000      0.000 3/19965
  920.0      1.515      1.515      1.515      0.000     -0.000      0.000 3/19965
  930.0      1.500      1.500      1.500      0.000     -0.000      0.000 3/19965
  940.0      1.484      1.484      1.484      0.000     -0.000      0.000 3/19965
  950.0      1.469      1.469      1.469      0.000     -0.000      0.000 3/19965
  960.0      1.454      1.454      1.454      0.000     -0.000      0.000 3/19965
  970.0      1.439      1.439      1.439      0.000     -0.000      0.000 3/19965
  980.0      1.425      1.425      1.425      0.000     -0.000      0.000 3/19965
  990.0      1.411      1.411      1.411      0.000     -0.000      0.000 3/19965
 1000.0      1.397      1.397      1.397      0.000     -0.000      0.000 3/19965

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

