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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [3 3 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.406825950000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.406825950000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.406825950000000
Atomic positions (fractional):
   *1 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998
    2 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998
    3 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998
   *4 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
   *5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.406825950000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.406825950000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.406825950000000
Atomic positions (fractional):
   *1 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998 > 1
    2 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998 > 2
    3 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998 > 3
   *4 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 4
   *5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   13.220477850000000    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.220477850000000    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.220477850000000
Atomic positions (fractional):
   *1 F   0.00000000000000  0.16666666666667  0.16666666666667  18.998 > 1
    2 F   0.33333333333333  0.16666666666667  0.16666666666667  18.998 > 1
    3 F   0.66666666666667  0.16666666666667  0.16666666666667  18.998 > 1
    4 F   0.00000000000000  0.50000000000000  0.16666666666667  18.998 > 1
    5 F   0.33333333333333  0.50000000000000  0.16666666666667  18.998 > 1
    6 F   0.66666666666667  0.50000000000000  0.16666666666667  18.998 > 1
    7 F   0.00000000000000  0.83333333333333  0.16666666666667  18.998 > 1
    8 F   0.33333333333333  0.83333333333333  0.16666666666667  18.998 > 1
    9 F   0.66666666666667  0.83333333333333  0.16666666666667  18.998 > 1
   10 F   0.00000000000000  0.16666666666667  0.50000000000000  18.998 > 1
   11 F   0.33333333333333  0.16666666666667  0.50000000000000  18.998 > 1
   12 F   0.66666666666667  0.16666666666667  0.50000000000000  18.998 > 1
   13 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998 > 1
   14 F   0.33333333333333  0.50000000000000  0.50000000000000  18.998 > 1
   15 F   0.66666666666667  0.50000000000000  0.50000000000000  18.998 > 1
   16 F   0.00000000000000  0.83333333333333  0.50000000000000  18.998 > 1
   17 F   0.33333333333333  0.83333333333333  0.50000000000000  18.998 > 1
   18 F   0.66666666666667  0.83333333333333  0.50000000000000  18.998 > 1
   19 F   0.00000000000000  0.16666666666667  0.83333333333333  18.998 > 1
   20 F   0.33333333333333  0.16666666666667  0.83333333333333  18.998 > 1
   21 F   0.66666666666667  0.16666666666667  0.83333333333333  18.998 > 1
   22 F   0.00000000000000  0.50000000000000  0.83333333333333  18.998 > 1
   23 F   0.33333333333333  0.50000000000000  0.83333333333333  18.998 > 1
   24 F   0.66666666666667  0.50000000000000  0.83333333333333  18.998 > 1
   25 F   0.00000000000000  0.83333333333333  0.83333333333333  18.998 > 1
   26 F   0.33333333333333  0.83333333333333  0.83333333333333  18.998 > 1
   27 F   0.66666666666667  0.83333333333333  0.83333333333333  18.998 > 1
   28 F   0.16666666666667  0.00000000000000  0.16666666666667  18.998 > 2
   29 F   0.50000000000000  0.00000000000000  0.16666666666667  18.998 > 2
   30 F   0.83333333333333  0.00000000000000  0.16666666666667  18.998 > 2
   31 F   0.16666666666667  0.33333333333333  0.16666666666667  18.998 > 2
   32 F   0.50000000000000  0.33333333333333  0.16666666666667  18.998 > 2
   33 F   0.83333333333333  0.33333333333333  0.16666666666667  18.998 > 2
   34 F   0.16666666666667  0.66666666666667  0.16666666666667  18.998 > 2
   35 F   0.50000000000000  0.66666666666667  0.16666666666667  18.998 > 2
   36 F   0.83333333333333  0.66666666666667  0.16666666666667  18.998 > 2
   37 F   0.16666666666667  0.00000000000000  0.50000000000000  18.998 > 2
   38 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998 > 2
   39 F   0.83333333333333  0.00000000000000  0.50000000000000  18.998 > 2
   40 F   0.16666666666667  0.33333333333333  0.50000000000000  18.998 > 2
   41 F   0.50000000000000  0.33333333333333  0.50000000000000  18.998 > 2
   42 F   0.83333333333333  0.33333333333333  0.50000000000000  18.998 > 2
   43 F   0.16666666666667  0.66666666666667  0.50000000000000  18.998 > 2
   44 F   0.50000000000000  0.66666666666667  0.50000000000000  18.998 > 2
   45 F   0.83333333333333  0.66666666666667  0.50000000000000  18.998 > 2
   46 F   0.16666666666667  0.00000000000000  0.83333333333333  18.998 > 2
   47 F   0.50000000000000  0.00000000000000  0.83333333333333  18.998 > 2
   48 F   0.83333333333333  0.00000000000000  0.83333333333333  18.998 > 2
   49 F   0.16666666666667  0.33333333333333  0.83333333333333  18.998 > 2
   50 F   0.50000000000000  0.33333333333333  0.83333333333333  18.998 > 2
   51 F   0.83333333333333  0.33333333333333  0.83333333333333  18.998 > 2
   52 F   0.16666666666667  0.66666666666667  0.83333333333333  18.998 > 2
   53 F   0.50000000000000  0.66666666666667  0.83333333333333  18.998 > 2
   54 F   0.83333333333333  0.66666666666667  0.83333333333333  18.998 > 2
   55 F   0.16666666666667  0.16666666666667  0.00000000000000  18.998 > 3
   56 F   0.50000000000000  0.16666666666667  0.00000000000000  18.998 > 3
   57 F   0.83333333333333  0.16666666666667  0.00000000000000  18.998 > 3
   58 F   0.16666666666667  0.50000000000000  0.00000000000000  18.998 > 3
   59 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998 > 3
   60 F   0.83333333333333  0.50000000000000  0.00000000000000  18.998 > 3
   61 F   0.16666666666667  0.83333333333333  0.00000000000000  18.998 > 3
   62 F   0.50000000000000  0.83333333333333  0.00000000000000  18.998 > 3
   63 F   0.83333333333333  0.83333333333333  0.00000000000000  18.998 > 3
   64 F   0.16666666666667  0.16666666666667  0.33333333333333  18.998 > 3
   65 F   0.50000000000000  0.16666666666667  0.33333333333333  18.998 > 3
   66 F   0.83333333333333  0.16666666666667  0.33333333333333  18.998 > 3
   67 F   0.16666666666667  0.50000000000000  0.33333333333333  18.998 > 3
   68 F   0.50000000000000  0.50000000000000  0.33333333333333  18.998 > 3
   69 F   0.83333333333333  0.50000000000000  0.33333333333333  18.998 > 3
   70 F   0.16666666666667  0.83333333333333  0.33333333333333  18.998 > 3
   71 F   0.50000000000000  0.83333333333333  0.33333333333333  18.998 > 3
   72 F   0.83333333333333  0.83333333333333  0.33333333333333  18.998 > 3
   73 F   0.16666666666667  0.16666666666667  0.66666666666667  18.998 > 3
   74 F   0.50000000000000  0.16666666666667  0.66666666666667  18.998 > 3
   75 F   0.83333333333333  0.16666666666667  0.66666666666667  18.998 > 3
   76 F   0.16666666666667  0.50000000000000  0.66666666666667  18.998 > 3
   77 F   0.50000000000000  0.50000000000000  0.66666666666667  18.998 > 3
   78 F   0.83333333333333  0.50000000000000  0.66666666666667  18.998 > 3
   79 F   0.16666666666667  0.83333333333333  0.66666666666667  18.998 > 3
   80 F   0.50000000000000  0.83333333333333  0.66666666666667  18.998 > 3
   81 F   0.83333333333333  0.83333333333333  0.66666666666667  18.998 > 3
  *82 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 4
   83 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 4
   84 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 4
   85 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 4
   86 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 4
   87 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 4
   88 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 4
   89 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 4
   90 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 4
   91 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 4
   92 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 4
   93 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 4
   94 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 4
   95 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 4
   96 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 4
   97 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 4
   98 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 4
   99 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 4
  100 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 4
  101 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 4
  102 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 4
  103 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 4
  104 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 4
  105 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 4
  106 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 4
  107 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 4
  108 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 4
 *109 Yb  0.16666666666667  0.16666666666667  0.16666666666667 173.054 > 5
  110 Yb  0.50000000000000  0.16666666666667  0.16666666666667 173.054 > 5
  111 Yb  0.83333333333333  0.16666666666667  0.16666666666667 173.054 > 5
  112 Yb  0.16666666666667  0.50000000000000  0.16666666666667 173.054 > 5
  113 Yb  0.50000000000000  0.50000000000000  0.16666666666667 173.054 > 5
  114 Yb  0.83333333333333  0.50000000000000  0.16666666666667 173.054 > 5
  115 Yb  0.16666666666667  0.83333333333333  0.16666666666667 173.054 > 5
  116 Yb  0.50000000000000  0.83333333333333  0.16666666666667 173.054 > 5
  117 Yb  0.83333333333333  0.83333333333333  0.16666666666667 173.054 > 5
  118 Yb  0.16666666666667  0.16666666666667  0.50000000000000 173.054 > 5
  119 Yb  0.50000000000000  0.16666666666667  0.50000000000000 173.054 > 5
  120 Yb  0.83333333333333  0.16666666666667  0.50000000000000 173.054 > 5
  121 Yb  0.16666666666667  0.50000000000000  0.50000000000000 173.054 > 5
  122 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 5
  123 Yb  0.83333333333333  0.50000000000000  0.50000000000000 173.054 > 5
  124 Yb  0.16666666666667  0.83333333333333  0.50000000000000 173.054 > 5
  125 Yb  0.50000000000000  0.83333333333333  0.50000000000000 173.054 > 5
  126 Yb  0.83333333333333  0.83333333333333  0.50000000000000 173.054 > 5
  127 Yb  0.16666666666667  0.16666666666667  0.83333333333333 173.054 > 5
  128 Yb  0.50000000000000  0.16666666666667  0.83333333333333 173.054 > 5
  129 Yb  0.83333333333333  0.16666666666667  0.83333333333333 173.054 > 5
  130 Yb  0.16666666666667  0.50000000000000  0.83333333333333 173.054 > 5
  131 Yb  0.50000000000000  0.50000000000000  0.83333333333333 173.054 > 5
  132 Yb  0.83333333333333  0.50000000000000  0.83333333333333 173.054 > 5
  133 Yb  0.16666666666667  0.83333333333333  0.83333333333333 173.054 > 5
  134 Yb  0.50000000000000  0.83333333333333  0.83333333333333 173.054 > 5
  135 Yb  0.83333333333333  0.83333333333333  0.83333333333333 173.054 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.4871409    0.0000000    0.0000000
            0.0000000    2.4871409    0.0000000
            0.0000000    0.0000000    2.4871409
-------------------------- Born effective charges --------------------------
    1 F    -1.7151624    0.0000000    0.0000000
            0.0000000   -1.0824469    0.0000000
            0.0000000    0.0000000   -1.0824469
    2 F    -1.0824469    0.0000000    0.0000000
            0.0000000   -1.7151624    0.0000000
            0.0000000    0.0000000   -1.0824469
    3 F    -1.0824469    0.0000000    0.0000000
            0.0000000   -1.0824469    0.0000000
            0.0000000    0.0000000   -1.7151624
    4 Cs    1.3619498    0.0000000    0.0000000
            0.0000000    1.3619498    0.0000000
            0.0000000    0.0000000    1.3619498
    5 Yb    2.5181065    0.0000000    0.0000000
            0.0000000    2.5181065    0.0000000
            0.0000000    0.0000000    2.5181065
----------------------------------------------------------------------------
Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
-------------------------------- Symfc start -------------------------------
Symfc version 1.5.4 (https://github.com/symfc/symfc)
Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024)
Computing [2] order force constants.
Permutation basis: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), data: True
-----
Solver_atoms: 1 -- 135 / 135
Time (Solver_compr_matrix_reshape): 0.003
Solver_block: 80 / 80
 - Time: 0.094
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.098
--------------------------------- 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-07 22:52:53]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

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


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

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

-------------------------[time 2026-01-07 22:52:58]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [3 3 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: Pm-3m (221)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.406825950000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.406825950000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.406825950000000
Atomic positions (fractional):
    1 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998
    2 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998
    3 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998
    4 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
    5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.220477850000000    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.220477850000000    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.220477850000000
Atomic positions (fractional):
    1 F   0.00000000000000  0.16666666666667  0.16666666666667  18.998 > 1
    2 F   0.33333333333333  0.16666666666667  0.16666666666667  18.998 > 1
    3 F   0.66666666666667  0.16666666666667  0.16666666666667  18.998 > 1
    4 F   0.00000000000000  0.50000000000000  0.16666666666667  18.998 > 1
    5 F   0.33333333333333  0.50000000000000  0.16666666666667  18.998 > 1
    6 F   0.66666666666667  0.50000000000000  0.16666666666667  18.998 > 1
    7 F   0.00000000000000  0.83333333333333  0.16666666666667  18.998 > 1
    8 F   0.33333333333333  0.83333333333333  0.16666666666667  18.998 > 1
    9 F   0.66666666666667  0.83333333333333  0.16666666666667  18.998 > 1
   10 F   0.00000000000000  0.16666666666667  0.50000000000000  18.998 > 1
   11 F   0.33333333333333  0.16666666666667  0.50000000000000  18.998 > 1
   12 F   0.66666666666667  0.16666666666667  0.50000000000000  18.998 > 1
   13 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998 > 1
   14 F   0.33333333333333  0.50000000000000  0.50000000000000  18.998 > 1
   15 F   0.66666666666667  0.50000000000000  0.50000000000000  18.998 > 1
   16 F   0.00000000000000  0.83333333333333  0.50000000000000  18.998 > 1
   17 F   0.33333333333333  0.83333333333333  0.50000000000000  18.998 > 1
   18 F   0.66666666666667  0.83333333333333  0.50000000000000  18.998 > 1
   19 F   0.00000000000000  0.16666666666667  0.83333333333333  18.998 > 1
   20 F   0.33333333333333  0.16666666666667  0.83333333333333  18.998 > 1
   21 F   0.66666666666667  0.16666666666667  0.83333333333333  18.998 > 1
   22 F   0.00000000000000  0.50000000000000  0.83333333333333  18.998 > 1
   23 F   0.33333333333333  0.50000000000000  0.83333333333333  18.998 > 1
   24 F   0.66666666666667  0.50000000000000  0.83333333333333  18.998 > 1
   25 F   0.00000000000000  0.83333333333333  0.83333333333333  18.998 > 1
   26 F   0.33333333333333  0.83333333333333  0.83333333333333  18.998 > 1
   27 F   0.66666666666667  0.83333333333333  0.83333333333333  18.998 > 1
   28 F   0.16666666666667  0.00000000000000  0.16666666666667  18.998 > 28
   29 F   0.50000000000000  0.00000000000000  0.16666666666667  18.998 > 28
   30 F   0.83333333333333  0.00000000000000  0.16666666666667  18.998 > 28
   31 F   0.16666666666667  0.33333333333333  0.16666666666667  18.998 > 28
   32 F   0.50000000000000  0.33333333333333  0.16666666666667  18.998 > 28
   33 F   0.83333333333333  0.33333333333333  0.16666666666667  18.998 > 28
   34 F   0.16666666666667  0.66666666666667  0.16666666666667  18.998 > 28
   35 F   0.50000000000000  0.66666666666667  0.16666666666667  18.998 > 28
   36 F   0.83333333333333  0.66666666666667  0.16666666666667  18.998 > 28
   37 F   0.16666666666667  0.00000000000000  0.50000000000000  18.998 > 28
   38 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998 > 28
   39 F   0.83333333333333  0.00000000000000  0.50000000000000  18.998 > 28
   40 F   0.16666666666667  0.33333333333333  0.50000000000000  18.998 > 28
   41 F   0.50000000000000  0.33333333333333  0.50000000000000  18.998 > 28
   42 F   0.83333333333333  0.33333333333333  0.50000000000000  18.998 > 28
   43 F   0.16666666666667  0.66666666666667  0.50000000000000  18.998 > 28
   44 F   0.50000000000000  0.66666666666667  0.50000000000000  18.998 > 28
   45 F   0.83333333333333  0.66666666666667  0.50000000000000  18.998 > 28
   46 F   0.16666666666667  0.00000000000000  0.83333333333333  18.998 > 28
   47 F   0.50000000000000  0.00000000000000  0.83333333333333  18.998 > 28
   48 F   0.83333333333333  0.00000000000000  0.83333333333333  18.998 > 28
   49 F   0.16666666666667  0.33333333333333  0.83333333333333  18.998 > 28
   50 F   0.50000000000000  0.33333333333333  0.83333333333333  18.998 > 28
   51 F   0.83333333333333  0.33333333333333  0.83333333333333  18.998 > 28
   52 F   0.16666666666667  0.66666666666667  0.83333333333333  18.998 > 28
   53 F   0.50000000000000  0.66666666666667  0.83333333333333  18.998 > 28
   54 F   0.83333333333333  0.66666666666667  0.83333333333333  18.998 > 28
   55 F   0.16666666666667  0.16666666666667  0.00000000000000  18.998 > 55
   56 F   0.50000000000000  0.16666666666667  0.00000000000000  18.998 > 55
   57 F   0.83333333333333  0.16666666666667  0.00000000000000  18.998 > 55
   58 F   0.16666666666667  0.50000000000000  0.00000000000000  18.998 > 55
   59 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998 > 55
   60 F   0.83333333333333  0.50000000000000  0.00000000000000  18.998 > 55
   61 F   0.16666666666667  0.83333333333333  0.00000000000000  18.998 > 55
   62 F   0.50000000000000  0.83333333333333  0.00000000000000  18.998 > 55
   63 F   0.83333333333333  0.83333333333333  0.00000000000000  18.998 > 55
   64 F   0.16666666666667  0.16666666666667  0.33333333333333  18.998 > 55
   65 F   0.50000000000000  0.16666666666667  0.33333333333333  18.998 > 55
   66 F   0.83333333333333  0.16666666666667  0.33333333333333  18.998 > 55
   67 F   0.16666666666667  0.50000000000000  0.33333333333333  18.998 > 55
   68 F   0.50000000000000  0.50000000000000  0.33333333333333  18.998 > 55
   69 F   0.83333333333333  0.50000000000000  0.33333333333333  18.998 > 55
   70 F   0.16666666666667  0.83333333333333  0.33333333333333  18.998 > 55
   71 F   0.50000000000000  0.83333333333333  0.33333333333333  18.998 > 55
   72 F   0.83333333333333  0.83333333333333  0.33333333333333  18.998 > 55
   73 F   0.16666666666667  0.16666666666667  0.66666666666667  18.998 > 55
   74 F   0.50000000000000  0.16666666666667  0.66666666666667  18.998 > 55
   75 F   0.83333333333333  0.16666666666667  0.66666666666667  18.998 > 55
   76 F   0.16666666666667  0.50000000000000  0.66666666666667  18.998 > 55
   77 F   0.50000000000000  0.50000000000000  0.66666666666667  18.998 > 55
   78 F   0.83333333333333  0.50000000000000  0.66666666666667  18.998 > 55
   79 F   0.16666666666667  0.83333333333333  0.66666666666667  18.998 > 55
   80 F   0.50000000000000  0.83333333333333  0.66666666666667  18.998 > 55
   81 F   0.83333333333333  0.83333333333333  0.66666666666667  18.998 > 55
   82 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 82
   83 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 82
   84 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 82
   85 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 82
   86 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 82
   87 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 82
   88 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 82
   89 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 82
   90 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 82
   91 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 82
   92 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 82
   93 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 82
   94 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 82
   95 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 82
   96 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 82
   97 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 82
   98 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 82
   99 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 82
  100 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 82
  101 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 82
  102 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 82
  103 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 82
  104 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 82
  105 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 82
  106 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 82
  107 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 82
  108 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 82
  109 Yb  0.16666666666667  0.16666666666667  0.16666666666667 173.054 > 109
  110 Yb  0.50000000000000  0.16666666666667  0.16666666666667 173.054 > 109
  111 Yb  0.83333333333333  0.16666666666667  0.16666666666667 173.054 > 109
  112 Yb  0.16666666666667  0.50000000000000  0.16666666666667 173.054 > 109
  113 Yb  0.50000000000000  0.50000000000000  0.16666666666667 173.054 > 109
  114 Yb  0.83333333333333  0.50000000000000  0.16666666666667 173.054 > 109
  115 Yb  0.16666666666667  0.83333333333333  0.16666666666667 173.054 > 109
  116 Yb  0.50000000000000  0.83333333333333  0.16666666666667 173.054 > 109
  117 Yb  0.83333333333333  0.83333333333333  0.16666666666667 173.054 > 109
  118 Yb  0.16666666666667  0.16666666666667  0.50000000000000 173.054 > 109
  119 Yb  0.50000000000000  0.16666666666667  0.50000000000000 173.054 > 109
  120 Yb  0.83333333333333  0.16666666666667  0.50000000000000 173.054 > 109
  121 Yb  0.16666666666667  0.50000000000000  0.50000000000000 173.054 > 109
  122 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 109
  123 Yb  0.83333333333333  0.50000000000000  0.50000000000000 173.054 > 109
  124 Yb  0.16666666666667  0.83333333333333  0.50000000000000 173.054 > 109
  125 Yb  0.50000000000000  0.83333333333333  0.50000000000000 173.054 > 109
  126 Yb  0.83333333333333  0.83333333333333  0.50000000000000 173.054 > 109
  127 Yb  0.16666666666667  0.16666666666667  0.83333333333333 173.054 > 109
  128 Yb  0.50000000000000  0.16666666666667  0.83333333333333 173.054 > 109
  129 Yb  0.83333333333333  0.16666666666667  0.83333333333333 173.054 > 109
  130 Yb  0.16666666666667  0.50000000000000  0.83333333333333 173.054 > 109
  131 Yb  0.50000000000000  0.50000000000000  0.83333333333333 173.054 > 109
  132 Yb  0.83333333333333  0.50000000000000  0.83333333333333 173.054 > 109
  133 Yb  0.16666666666667  0.83333333333333  0.83333333333333 173.054 > 109
  134 Yb  0.50000000000000  0.83333333333333  0.83333333333333 173.054 > 109
  135 Yb  0.83333333333333  0.83333333333333  0.83333333333333 173.054 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.4871409    0.0000000    0.0000000
            0.0000000    2.4871409    0.0000000
            0.0000000    0.0000000    2.4871409
-------------------------- Born effective charges --------------------------
    1 F    -1.7151624    0.0000000    0.0000000
            0.0000000   -1.0824469    0.0000000
            0.0000000    0.0000000   -1.0824469
    2 F    -1.0824469    0.0000000    0.0000000
            0.0000000   -1.7151624    0.0000000
            0.0000000    0.0000000   -1.0824469
    3 F    -1.0824469    0.0000000    0.0000000
            0.0000000   -1.0824469    0.0000000
            0.0000000    0.0000000   -1.7151624
    4 Cs    1.3619498    0.0000000    0.0000000
            0.0000000    1.3619498    0.0000000
            0.0000000    0.0000000    1.3619498
    5 Yb    2.5181065    0.0000000    0.0000000
            0.0000000    2.5181065    0.0000000
            0.0000000    0.0000000    2.5181065
----------------------------------------------------------------------------
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.94, Number of G-points: 305, Lambda: 0.15
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 0.00e+00 0.00e+00 0.00e+00 8.55e-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
   2.408   (   0.000    0.000    0.000)    0.000
   2.408   (   0.000    0.000    0.000)    0.000
   2.408   (   0.000    0.000    0.000)    0.000
   4.126   (   0.000    0.000    0.000)    0.000
   4.126   (   0.000    0.000    0.000)    0.000
   4.126   (   0.000    0.000    0.000)    0.000
   4.660   (   0.000    0.000    0.000)    0.000
   4.660   (   0.000    0.000    0.000)    0.000
   4.660   (   0.000    0.000    0.000)    0.000
  10.569   (   0.000    0.000    0.000)    0.000
  10.569   (   0.000    0.000    0.000)    0.000
  10.569   (   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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.403   (  19.226    0.000    0.000)   19.226
   0.403   (  19.226    0.000    0.000)   19.226
   0.893   (  40.567    0.000    0.000)   40.567
   2.392   (  -1.613    0.000    0.000)    1.613
   2.392   (  -1.613    0.000    0.000)    1.613
   2.620   (  10.394    0.000    0.000)   10.394
   4.161   (   3.316    0.000    0.000)    3.316
   4.162   (   3.462    0.000    0.000)    3.462
   4.162   (   3.462    0.000    0.000)    3.462
   4.671   (   0.972    0.000    0.000)    0.972
   4.671   (   0.972    0.000    0.000)    0.972
   6.872   (   5.104    0.000    0.000)    5.104
  10.568   (  -0.150    0.000    0.000)    0.150
  10.568   (  -0.150    0.000    0.000)    0.150
  12.497   (   5.382    0.000    0.000)    5.382
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.787   (  17.882    0.000    0.000)   17.882
   0.787   (  17.882    0.000    0.000)   17.882
   1.634   (  30.810    0.000    0.000)   30.810
   2.340   (  -3.432    0.000    0.000)    3.432
   2.340   (  -3.432    0.000    0.000)    3.432
   2.912   (  16.965    0.000    0.000)   16.965
   4.254   (   5.446    0.000    0.000)    5.446
   4.262   (   5.992    0.000    0.000)    5.992
   4.262   (   5.992    0.000    0.000)    5.992
   4.697   (   1.479    0.000    0.000)    1.479
   4.697   (   1.479    0.000    0.000)    1.479
   7.015   (   8.280    0.000    0.000)    8.280
  10.563   (  -0.255    0.000    0.000)    0.255
  10.563   (  -0.255    0.000    0.000)    0.255
  12.654   (   9.563    0.000    0.000)    9.563
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.137   (  15.976    0.000    0.000)   15.976
   1.137   (  15.976    0.000    0.000)   15.976
   2.165   (  20.991    0.000    0.000)   20.991
   2.249   (  -5.426    0.000    0.000)    5.426
   2.249   (  -5.426    0.000    0.000)    5.426
   3.270   (  16.545    0.000    0.000)   16.545
   4.373   (   5.749    0.000    0.000)    5.749
   4.399   (   6.955    0.000    0.000)    6.955
   4.399   (   6.955    0.000    0.000)    6.955
   4.726   (   1.166    0.000    0.000)    1.166
   4.726   (   1.166    0.000    0.000)    1.166
   7.193   (   8.504    0.000    0.000)    8.504
  10.558   (  -0.282    0.000    0.000)    0.282
  10.558   (  -0.282    0.000    0.000)    0.282
  12.873   (  11.077    0.000    0.000)   11.077
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.445   (  13.775    0.000    0.000)   13.775
   1.445   (  13.775    0.000    0.000)   13.775
   2.116   (  -7.424    0.000    0.000)    7.424
   2.116   (  -7.424    0.000    0.000)    7.424
   2.511   (  12.670    0.000    0.000)   12.670
   3.558   (  10.864    0.000    0.000)   10.864
   4.479   (   4.275    0.000    0.000)    4.275
   4.538   (   6.275    0.000    0.000)    6.275
   4.538   (   6.275    0.000    0.000)    6.275
   4.739   (  -0.093    0.000    0.000)    0.093
   4.739   (  -0.093    0.000    0.000)    0.093
   7.348   (   6.107    0.000    0.000)    6.107
  10.552   (  -0.218    0.000    0.000)    0.218
  10.552   (  -0.218    0.000    0.000)    0.218
  13.085   (   8.878    0.000    0.000)    8.878
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.701   (  10.702    0.000    0.000)   10.702
   1.701   (  10.702    0.000    0.000)   10.702
   1.946   (  -8.589    0.000    0.000)    8.589
   1.946   (  -8.589    0.000    0.000)    8.589
   2.687   (   4.329    0.000    0.000)    4.329
   3.709   (   3.649    0.000    0.000)    3.649
   4.541   (   1.570    0.000    0.000)    1.570
   4.649   (   4.417    0.000    0.000)    4.417
   4.649   (   4.417    0.000    0.000)    4.417
   4.717   (  -2.119    0.000    0.000)    2.119
   4.717   (  -2.119    0.000    0.000)    2.119
   7.435   (   2.189    0.000    0.000)    2.189
  10.549   (  -0.082    0.000    0.000)    0.082
  10.549   (  -0.082    0.000    0.000)    0.082
  13.216   (   3.393    0.000    0.000)    3.393
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.567   (  13.321   13.321    0.000)   18.839
   0.790   (  17.372   17.372    0.000)   24.567
   1.130   (  25.930   25.930    0.000)   36.671
   2.377   (  -1.490   -1.490    0.000)    2.107
   2.540   (   6.294    6.294    0.000)    8.901
   2.556   (   2.095    2.095    0.000)    2.963
   4.102   (  -1.478   -1.478    0.000)    2.091
   4.198   (   3.361    3.361    0.000)    4.754
   4.278   (   6.904    6.904    0.000)    9.764
   4.674   (   0.850    0.850    0.000)    1.203
   4.682   (   1.014    1.014    0.000)    1.435
   6.925   (   4.586    4.586    0.000)    6.486
  10.566   (  -0.160   -0.160    0.000)    0.227
  10.639   (   3.424    3.424    0.000)    4.843
  12.504   (   3.336    3.336    0.000)    4.717
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.875   (  15.639    7.936    0.000)   17.538
   1.071   (  12.761   19.351    0.000)   23.180
   1.735   (  27.600   10.535    0.000)   29.542
   2.328   (  -3.209   -1.142    0.000)    3.407
   2.462   (  -4.309   11.075    0.000)   11.884
   2.875   (  18.239   -2.609    0.000)   18.424
   4.102   (   1.315   -7.120    0.000)    7.241
   4.295   (   5.826    3.095    0.000)    6.598
   4.432   (   7.586    9.197    0.000)   11.922
   4.709   (   1.557    1.135    0.000)    1.926
   4.710   (   2.650    1.086    0.000)    2.864
   7.047   (   6.957    2.874    0.000)    7.527
  10.561   (  -0.273   -0.189    0.000)    0.332
  10.678   (   0.832    9.067    0.000)    9.105
  12.649   (  10.004    0.916    0.000)   10.046
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.192   (  14.804    5.043    0.000)   15.640
   1.337   (  12.897   14.640    0.000)   19.511
   2.204   (  16.535    3.272    0.000)   16.856
   2.243   (  -5.141   -0.630    0.000)    5.179
   2.372   (  -3.125   12.738    0.000)   13.116
   3.254   (  17.332   -1.112    0.000)   17.368
   4.148   (   2.937  -11.871    0.000)   12.229
   4.428   (   6.770    2.746    0.000)    7.305
   4.575   (   5.959    8.293    0.000)   10.212
   4.740   (   1.269    1.316    0.000)    1.828
   4.780   (   3.972    4.064    0.000)    5.683
   7.196   (   7.063    0.274    0.000)    7.068
  10.555   (  -0.301   -0.226    0.000)    0.377
  10.687   (   0.166   11.370    0.000)   11.371
  12.885   (  12.064    1.725    0.000)   12.187
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.480   (  13.068    3.297    0.000)   13.477
   1.596   (  11.994   10.704    0.000)   16.076
   2.116   (  -7.118   -0.024    0.000)    7.118
   2.212   (  -6.869    8.103    0.000)   10.623
   2.544   (  10.999    5.279    0.000)   12.200
   3.554   (  11.186   -0.374    0.000)   11.192
   4.212   (   2.922  -14.835    0.000)   15.120
   4.563   (   6.102    2.390    0.000)    6.553
   4.670   (   3.243    4.857    0.000)    5.840
   4.755   (   0.026    1.537    0.000)    1.537
   4.864   (   3.823    8.919    0.000)    9.703
   7.324   (   5.090   -2.048    0.000)    5.486
  10.550   (  -0.233   -0.261    0.000)    0.350
  10.689   (   0.050   12.397    0.000)   12.397
  13.116   (   9.628    3.193    0.000)   10.144
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.727   (  10.580    2.313    0.000)   10.830
   1.825   (   9.763    6.279    0.000)   11.608
   1.951   (  -8.585    0.486    0.000)    8.599
   2.052   (  -8.108   11.775    0.000)   14.297
   2.699   (   3.830    2.779    0.000)    4.733
   3.707   (   3.700   -0.384    0.000)    3.720
   4.257   (   1.253  -15.896    0.000)   15.945
   4.671   (   4.266    2.069    0.000)    4.742
   4.713   (   0.978    2.722    0.000)    2.893
   4.736   (  -1.986    1.787    0.000)    2.672
   4.922   (   1.542   11.705    0.000)   11.806
   7.397   (   1.836   -3.332    0.000)    3.804
  10.546   (  -0.088   -0.282    0.000)    0.295
  10.689   (   0.017   12.873    0.000)   12.873
  13.257   (   3.659    4.101    0.000)    5.496
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.083   (  11.554   11.554    0.000)   16.339
   1.386   (  11.345   11.345    0.000)   16.045
   2.035   (  15.393   15.393    0.000)   21.769
   2.291   (  -2.566   -2.566    0.000)    3.629
   2.710   (   7.082    7.082    0.000)   10.016
   2.898   (  10.347   10.347    0.000)   14.633
   3.980   (  -4.462   -4.462    0.000)    6.311
   4.384   (   5.378    5.378    0.000)    7.606
   4.628   (   9.017    9.017    0.000)   12.752
   4.741   (   1.779    1.779    0.000)    2.516
   4.742   (   2.352    2.352    0.000)    3.326
   7.121   (   4.043    4.043    0.000)    5.717
  10.556   (  -0.321   -0.321    0.000)    0.454
  10.854   (   6.950    6.950    0.000)    9.828
  12.743   (   8.458    8.458    0.000)   11.962
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.330   (  12.041    7.996    0.000)   14.455
   1.590   (   9.131    9.215    0.000)   12.973
   2.221   (  -4.285   -1.575    0.000)    4.565
   2.266   (   6.086    2.149    0.000)    6.454
   2.769   (   2.910   23.570    0.000)   23.749
   3.247   (  17.840    0.771    0.000)   17.857
   3.922   (  -1.158   -9.660    0.000)    9.729
   4.507   (   6.258    4.771    0.000)    7.869
   4.739   (   2.188    7.347    0.000)    7.666
   4.777   (   1.565    2.117    0.000)    2.633
   4.858   (   8.145    3.369    0.000)    8.814
   7.205   (   3.941    0.664    0.000)    3.997
  10.549   (  -0.354   -0.384    0.000)    0.523
  10.958   (   3.433   13.329    0.000)   13.764
  12.974   (  12.936    7.489    0.000)   14.947
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.571   (  11.128    5.224    0.000)   12.294
   1.768   (   7.996    5.931    0.000)    9.956
   2.113   (  -6.132   -0.316    0.000)    6.140
   2.300   (  -1.674   -0.399    0.000)    1.721
   2.855   (   4.607   23.097    0.000)   23.552
   3.542   (   9.978   -1.040    0.000)   10.032
   3.932   (   2.114  -11.421    0.000)   11.616
   4.632   (   5.612    4.132    0.000)    6.969
   4.769   (   0.912    4.846    0.000)    4.931
   4.799   (   0.371    2.534    0.000)    2.561
   5.012   (   6.200    4.756    0.000)    7.814
   7.277   (   2.891   -2.119    0.000)    3.584
  10.542   (  -0.274   -0.443    0.000)    0.521
  11.008   (   1.602   16.978    0.000)   17.053
  13.230   (  10.881    8.170    0.000)   13.607
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.786   (   9.655    3.064    0.000)   10.130
   1.902   (   4.163    2.657    0.000)    4.939
   1.968   (  -7.950    1.153    0.000)    8.034
   2.246   (  -2.399    2.495    0.000)    3.462
   2.927   (   1.954   20.064    0.000)   20.159
   3.662   (   2.367   -5.338    0.000)    5.839
   3.985   (   1.956   -8.863    0.000)    9.077
   4.731   (   3.827    3.536    0.000)    5.211
   4.779   (   0.182    3.568    0.000)    3.572
   4.787   (  -1.596    3.010    0.000)    3.407
   5.101   (   2.273    5.251    0.000)    5.722
   7.319   (   1.074   -3.600    0.000)    3.757
  10.538   (  -0.103   -0.478    0.000)    0.489
  11.029   (   0.494   18.535    0.000)   18.542
  13.391   (   4.167    8.951    0.000)    9.873
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.507   (   8.802    8.802    0.000)   12.448
   1.745   (   6.453    6.453    0.000)    9.125
   2.175   (  -2.879   -2.879    0.000)    4.071
   2.284   (  -0.444   -0.444    0.000)    0.628
   3.215   (  14.269   14.269    0.000)   20.180
   3.308   (   8.123    8.123    0.000)   11.487
   3.764   (  -5.056   -5.056    0.000)    7.150
   4.616   (   5.534    5.534    0.000)    7.826
   4.821   (   2.019    2.019    0.000)    2.856
   4.850   (   2.745    2.745    0.000)    3.882
   4.941   (   5.364    5.364    0.000)    7.586
   7.224   (   1.088    1.088    0.000)    1.538
  10.540   (  -0.423   -0.423    0.000)    0.598
  11.193   (   9.059    9.059    0.000)   12.812
  13.189   (  12.603   12.603    0.000)   17.823
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.684   (   8.281    5.570    0.000)    9.980
   1.867   (   5.382    4.125    0.000)    6.781
   2.101   (  -4.329   -0.935    0.000)    4.429
   2.255   (  -2.125   -3.155    0.000)    3.803
   3.330   (   3.134   19.213    0.000)   19.467
   3.498   (   2.579   -2.083    0.000)    3.315
   3.770   (   6.300   -4.112    0.000)    7.523
   4.726   (   4.890    4.735    0.000)    6.807
   4.853   (   0.906    2.587    0.000)    2.741
   4.870   (  -0.199    4.983    0.000)    4.987
   5.059   (   5.355    0.163    0.000)    5.357
   7.248   (   1.176   -0.621    0.000)    1.330
  10.532   (  -0.327   -0.487    0.000)    0.586
  11.333   (   4.688   13.613    0.000)   14.397
  13.449   (  11.398   12.484    0.000)   16.904
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.841   (   6.436    2.068    0.000)    6.760
   1.956   (   2.687    2.616    0.000)    3.750
   2.000   (  -5.069    1.921    0.000)    5.421
   2.210   (  -1.608   -3.811    0.000)    4.137
   3.379   (   1.441   19.833    0.000)   19.885
   3.491   (  -1.013   -8.779    0.000)    8.837
   3.897   (   3.785   -1.322    0.000)    4.009
   4.810   (   3.158    3.956    0.000)    5.062
   4.853   (  -0.993    3.238    0.000)    3.386
   4.859   (  -0.501    4.305    0.000)    4.334
   5.141   (   2.207   -0.934    0.000)    2.396
   7.267   (   0.530   -1.317    0.000)    1.419
  10.528   (  -0.123   -0.525    0.000)    0.539
  11.395   (   1.440   15.821    0.000)   15.886
  13.620   (   4.476   12.698    0.000)   13.464
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.793   (   4.741    4.741    0.000)    6.705
   1.946   (   3.556    3.556    0.000)    5.030
   2.074   (  -1.612   -1.612    0.000)    2.280
   2.187   (  -3.215   -3.215    0.000)    4.547
   3.374   (  -4.204   -4.204    0.000)    5.945
   3.693   (   8.020    8.020    0.000)   11.342
   3.747   (   3.777    3.777    0.000)    5.341
   4.818   (   4.050    4.050    0.000)    5.727
   4.898   (   1.580    1.580    0.000)    2.234
   4.959   (   2.399    2.399    0.000)    3.392
   5.050   (   0.416    0.416    0.000)    0.588
   7.249   (   0.529    0.529    0.000)    0.748
  10.523   (  -0.375   -0.375    0.000)    0.531
  11.553   (   7.742    7.742    0.000)   10.948
  13.707   (  11.440   11.440    0.000)   16.179
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.870   (   2.359    0.880    0.000)    2.517
   2.008   (   2.068    2.332    0.000)    3.117
   2.041   (  -1.272    1.805    0.000)    2.208
   2.130   (  -1.917   -3.507    0.000)    3.997
   3.310   (  -1.669   -6.819    0.000)    7.020
   3.721   (   0.069   11.543    0.000)   11.543
   3.882   (   3.798   -0.639    0.000)    3.851
   4.886   (   2.350    3.186    0.000)    3.959
   4.913   (  -0.245    2.347    0.000)    2.360
   4.957   (  -0.884    4.864    0.000)    4.944
   5.092   (   1.898   -3.035    0.000)    3.580
   7.260   (   0.387    0.413    0.000)    0.567
  10.518   (  -0.141   -0.404    0.000)    0.428
  11.657   (   2.441    9.420    0.000)    9.731
  13.880   (   4.563   11.484    0.000)   12.357
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.881   (   0.228    0.228    0.000)    0.323
   2.046   (   1.255    1.255    0.000)    1.776
   2.068   (   0.701    0.701    0.000)    0.991
   2.076   (  -1.618   -1.618    0.000)    2.289
   3.214   (  -2.394   -2.394    0.000)    3.385
   3.867   (   1.516    1.516    0.000)    2.143
   3.871   (   1.339    1.339    0.000)    1.893
   4.935   (   1.477    1.477    0.000)    2.089
   4.944   (   0.593    0.593    0.000)    0.839
   5.034   (   1.039    1.039    0.000)    1.469
   5.046   (  -0.091   -0.091    0.000)    0.129
   7.271   (   0.390    0.390    0.000)    0.552
  10.512   (  -0.152   -0.152    0.000)    0.214
  11.784   (   3.024    3.024    0.000)    4.277
  14.053   (   4.591    4.591    0.000)    6.493
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.884   (  13.018   13.018   13.018)   22.548
   0.884   (  13.018   13.018   13.018)   22.548
   1.320   (  20.045   20.045   20.045)   34.719
   2.495   (  -0.470   -0.470   -0.470)    0.815
   2.525   (   3.736    3.736    3.736)    6.471
   2.525   (   3.736    3.736    3.736)    6.471
   4.056   (  -2.301   -2.301   -2.301)    3.986
   4.283   (   3.571    3.571    3.571)    6.185
   4.283   (   3.571    3.571    3.571)    6.185
   4.702   (   2.213    2.213    2.213)    3.834
   4.702   (   2.213    2.213    2.213)    3.834
   6.971   (   4.017    4.017    4.017)    6.958
  10.639   (   2.357    2.357    2.357)    4.082
  10.639   (   2.357    2.357    2.357)    4.082
  12.515   (   2.845    2.845    2.845)    4.927
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.105   (  12.304   11.745   11.745)   20.671
   1.158   (  11.395   11.845   11.845)   20.260
   1.831   (  24.537    9.734    9.734)   28.134
   2.415   (  -4.038    2.949    2.949)    5.805
   2.476   (  -3.259    6.553    6.553)    9.823
   2.853   (  18.823   -1.550   -1.550)   18.950
   4.033   (   0.042   -5.057   -5.057)    7.151
   4.343   (   4.360    1.708    1.708)    4.984
   4.403   (   5.465    0.473    0.473)    5.506
   4.744   (   2.337    4.074    4.074)    6.218
   4.786   (   5.211    4.863    4.863)    8.629
   7.074   (   5.737    2.376    2.376)    6.649
  10.641   (   0.084    3.748    3.748)    5.302
  10.705   (   1.709    4.897    4.897)    7.133
  12.653   (  10.163    1.493    1.493)   10.380
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.369   (  12.915    8.896    8.896)   18.030
   1.392   (  11.294    9.012    9.012)   17.029
   2.218   (   8.024    0.811    0.811)    8.105
   2.367   (   4.084    7.628    7.628)   11.534
   2.389   (  -5.218    6.881    6.881)   11.042
   3.245   (  17.811   -0.682   -0.682)   17.837
   4.053   (   1.742   -7.899   -7.899)   11.306
   4.450   (   5.704   -0.075   -0.075)    5.705
   4.501   (   3.937   -1.914   -1.914)    4.778
   4.795   (   2.480    5.196    5.196)    7.756
   4.906   (   6.013    7.927    7.927)   12.721
   7.195   (   5.662   -0.094   -0.094)    5.663
  10.643   (   0.096    4.065    4.065)    5.750
  10.723   (   0.313    6.872    6.872)    9.724
  12.900   (  12.820    1.976    1.976)   13.121
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.619   (  10.616    5.998    5.998)   13.589
   1.628   (  11.976    6.750    6.750)   15.315
   2.169   (  -6.047    1.545    1.545)    6.430
   2.261   (  -7.081    7.497    7.497)   12.749
   2.580   (   9.878    5.133    5.133)   12.258
   3.550   (  11.337   -0.334   -0.334)   11.347
   4.094   (   1.969  -10.047  -10.047)   14.344
   4.565   (   5.092   -0.906   -0.906)    5.250
   4.566   (   2.371   -3.426   -3.426)    5.394
   4.839   (   1.600    6.273    6.273)    9.014
   5.018   (   4.591    9.750    9.750)   14.533
   7.297   (   4.062   -2.389   -2.389)    5.283
  10.645   (   0.077    4.364    4.364)    6.172
  10.725   (   0.013    7.621    7.621)   10.778
  13.146   (  10.261    3.153    3.153)   11.187
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.822   (   8.521    2.680    2.680)    9.326
   1.852   (   9.085    3.921    3.921)   10.644
   2.028   (  -7.033    4.291    4.291)    9.289
   2.104   (  -7.428    9.602    9.602)   15.479
   2.719   (   3.443    3.334    3.334)    5.839
   3.705   (   3.733   -0.425   -0.425)    3.781
   4.125   (   0.838  -11.173  -11.173)   15.823
   4.598   (   0.799   -4.068   -4.068)    5.808
   4.644   (   2.203   -1.397   -1.397)    2.959
   4.858   (   0.374    7.234    7.234)   10.237
   5.085   (   1.691   10.493   10.493)   14.936
   7.356   (   1.471   -3.697   -3.697)    5.431
  10.646   (   0.029    4.544    4.544)    6.426
  10.725   (  -0.007    7.893    7.893)   11.163
  13.297   (   3.889    3.923    3.923)    6.775
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.289   (   9.412    9.412   15.015)   20.065
   1.430   (  10.643   10.643    4.187)   15.623
   2.090   (  12.135   12.135    4.891)   17.845
   2.401   (  -1.116   -1.116   10.879)   10.993
   2.719   (   8.039    8.039    1.583)   11.479
   2.886   (  10.409   10.409   -1.146)   14.765
   3.933   (  -4.091   -4.091   -4.076)    7.077
   4.361   (   1.935    1.935   -1.345)    3.050
   4.463   (   3.375    3.375   -8.356)    9.623
   4.843   (   4.506    4.506    8.473)   10.601
   4.902   (   6.824    6.824    8.377)   12.779
   7.132   (   3.004    3.004    0.904)    4.343
  10.663   (   0.166    0.166    9.183)    9.186
  10.858   (   6.994    6.994    0.354)    9.897
  12.755   (   8.758    8.758    1.910)   12.532
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.494   (  10.139    6.748   12.131)   17.190
   1.622   (   8.543    8.807    2.990)   12.628
   2.217   (   0.025   -1.018   -0.821)    1.308
   2.407   (   0.977    1.462   14.414)   14.521
   2.786   (   3.185   23.758    1.714)   24.032
   3.242   (  17.991    0.775   -0.321)   18.010
   3.874   (  -1.586   -8.564   -4.465)    9.788
   4.429   (   4.218    0.243   -4.821)    6.410
   4.526   (   2.799    2.067  -10.891)   11.433
   4.922   (   3.272    6.297   10.068)   12.318
   5.055   (   7.392    6.495   10.583)   14.451
   7.192   (   2.652   -0.198   -1.323)    2.970
  10.664   (  -0.003    0.007   10.167)   10.167
  10.965   (   3.529   12.973    0.734)   13.465
  12.994   (  13.357    7.662    2.305)   15.570
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.695   (   9.081    4.220    8.566)   13.178
   1.790   (   7.646    5.803    2.174)    9.841
   2.163   (  -4.105   -1.093    2.141)    4.757
   2.385   (  -2.970    1.950   11.395)   11.936
   2.875   (   4.596   22.367    2.241)   22.944
   3.538   (  10.009   -1.111   -0.395)   10.078
   3.869   (   1.010  -10.297   -6.019)   11.970
   4.509   (   2.860   -0.673   -7.633)    8.179
   4.588   (   3.418    1.493   -9.253)    9.976
   4.978   (   2.065    6.669   10.663)   12.746
   5.189   (   5.294    6.287   11.263)   13.943
   7.239   (   1.866   -2.898   -3.541)    4.941
  10.664   (  -0.023   -0.167   10.842)   10.843
  11.015   (   1.575   16.555    0.877)   16.653
  13.259   (  11.266    8.044    2.963)   14.157
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.854   (   5.512    1.415    2.738)    6.315
   1.921   (   4.306    2.748    2.109)    5.526
   2.070   (  -4.180    1.482    7.324)    8.562
   2.316   (  -2.653    4.737   10.197)   11.552
   2.946   (   1.897   19.479    2.290)   19.705
   3.659   (   2.405   -5.592   -0.295)    6.094
   3.900   (   1.320   -7.871   -7.968)   11.278
   4.542   (   0.639   -1.221   -9.885)    9.980
   4.651   (   1.924    1.398   -6.364)    6.794
   5.005   (   0.619    6.717   11.015)   12.916
   5.265   (   1.917    6.195   11.274)   13.006
   7.266   (   0.694   -4.372   -4.886)    6.593
  10.663   (  -0.010   -0.274   11.219)   11.223
  11.035   (   0.470   18.100    0.844)   18.126
  13.425   (   4.314    8.645    3.430)   10.252
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.643   (   7.425    7.425   10.036)   14.525
   1.770   (   6.146    6.146    2.350)    9.004
   2.183   (  -2.093   -2.093   -0.601)    3.020
   2.419   (  -0.516   -0.516   14.855)   14.873
   3.236   (  13.761   13.761    2.089)   19.573
   3.295   (   8.050    8.050   -1.242)   11.452
   3.726   (  -5.164   -5.164   -3.618)    8.150
   4.469   (   3.259    3.259   -8.935)   10.054
   4.567   (   1.963    1.963  -13.791)   14.068
   5.034   (   4.409    4.409   12.418)   13.895
   5.178   (   5.426    5.426   12.229)   14.437
   7.189   (  -0.051   -0.051   -3.311)    3.312
  10.661   (  -0.200   -0.200   10.984)   10.988
  11.196   (   8.963    8.963    0.258)   12.678
  13.212   (  12.736   12.736    2.463)   18.180
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.787   (   6.147    4.431    6.931)   10.270
   1.887   (   5.236    3.998    1.940)    6.868
   2.136   (  -2.299   -1.378    1.069)    2.886
   2.387   (  -2.410   -1.029   15.544)   15.764
   3.331   (   2.975   17.650    0.205)   17.901
   3.493   (   2.184   -1.896   -0.427)    2.923
   3.718   (   5.333   -3.832   -4.913)    8.202
   4.534   (   2.102    2.961  -11.525)   12.084
   4.619   (   3.423    1.375  -12.317)   12.857
   5.102   (   2.238    5.281   12.988)   14.198
   5.275   (   3.923    2.333   13.025)   13.802
   7.191   (   0.240   -1.599   -5.371)    5.609
  10.657   (  -0.178   -0.381   11.421)   11.428
  11.334   (   4.597   13.389    0.097)   14.156
  13.475   (  11.528   12.337    2.636)   17.089
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.882   (   2.685    1.279    2.326)    3.776
   1.975   (   2.702    2.510    1.857)    4.129
   2.095   (  -1.356    0.739    4.220)    4.493
   2.337   (  -1.763   -0.641   16.510)   16.616
   3.380   (   1.496   18.346    0.277)   18.409
   3.480   (  -1.152   -8.835   -1.056)    8.972
   3.825   (   3.082   -1.186   -6.826)    7.583
   4.553   (   0.232    2.398  -13.520)   13.733
   4.685   (   2.087    1.736   -8.938)    9.341
   5.130   (   0.569    5.335   13.138)   14.191
   5.333   (   1.519    0.744   12.930)   13.040
   7.196   (   0.177   -2.179   -6.646)    6.997
  10.655   (  -0.069   -0.492   11.646)   11.657
  11.394   (   1.396   15.544   -0.029)   15.606
  13.647   (   4.531   12.381    2.783)   13.475
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.866   (   2.941    2.941    4.860)    6.397
   1.963   (   3.427    3.427    1.558)    5.091
   2.112   (  -0.926   -0.926    1.101)    1.710
   2.357   (  -1.612   -1.612   18.928)   19.065
   3.362   (  -4.258   -4.258   -1.205)    6.142
   3.674   (   6.648    6.648   -1.812)    9.575
   3.690   (   3.243    3.243   -5.366)    7.059
   4.616   (   3.537    3.537  -10.957)   12.045
   4.646   (   2.078    2.078  -13.993)   14.299
   5.190   (   3.026    3.026   13.643)   14.299
   5.302   (   0.734    0.734   14.089)   14.128
   7.173   (  -0.210   -0.210   -7.135)    7.142
  10.650   (  -0.320   -0.320   11.636)   11.645
  11.550   (   7.576    7.576   -0.301)   10.718
  13.730   (  11.355   11.355    2.360)   16.231
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.903   (   0.793    0.733    2.639)    2.851
   2.023   (   2.065    2.146    1.470)    3.321
   2.096   (  -0.626   -0.561    1.398)    1.631
   2.329   (  -0.855    0.080   21.530)   21.547
   3.297   (  -1.670   -6.801   -1.242)    7.113
   3.684   (   0.102    9.290   -3.487)    9.923
   3.809   (   2.789   -0.693   -6.822)    7.403
   4.638   (   0.125    5.210  -12.880)   13.894
   4.715   (   2.416    1.059  -10.230)   10.565
   5.224   (   0.507    3.696   13.788)   14.283
   5.318   (   0.693   -1.544   13.869)   13.972
   7.174   (   0.124   -0.157   -8.141)    8.143
  10.645   (  -0.123   -0.407   11.727)   11.735
  11.651   (   2.374    9.216   -0.534)    9.532
  13.902   (   4.540   11.273    2.233)   12.356
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.913   (   0.191    0.191    2.908)    2.921
   2.059   (   1.214    1.214    1.238)    2.117
   2.081   (  -0.759   -0.759    1.105)    1.541
   2.335   (   0.267    0.267   23.114)   23.117
   3.202   (  -2.368   -2.368   -1.168)    3.546
   3.791   (   1.221    1.221   -7.201)    7.405
   3.799   (   0.555    0.555   -6.737)    6.782
   4.728   (   1.564    1.564  -10.294)   10.529
   4.729   (   1.432    1.432  -10.750)   10.940
   5.276   (   1.082    1.082   13.605)   13.691
   5.293   (  -0.461   -0.461   13.779)   13.794
   7.177   (   0.202    0.202   -8.897)    8.902
  10.639   (  -0.156   -0.156   11.737)   11.739
  11.776   (   2.945    2.945   -0.822)    4.245
  14.072   (   4.522    4.522    1.944)    6.684
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.549   (   8.285    8.285    8.285)   14.350
   1.549   (   8.285    8.285    8.285)   14.350
   2.178   (   2.406    2.406    2.406)    4.167
   2.691   (  11.967   11.967   11.967)   20.728
   2.850   (   6.355    6.355    6.355)   11.008
   2.850   (   6.355    6.355    6.355)   11.008
   3.837   (  -4.736   -4.736   -4.736)    8.203
   4.347   (  -1.389   -1.389   -1.389)    2.405
   4.347   (  -1.389   -1.389   -1.389)    2.405
   5.043   (   7.789    7.789    7.789)   13.491
   5.043   (   7.789    7.789    7.789)   13.491
   7.149   (   0.511    0.511    0.511)    0.885
  10.868   (   4.945    4.945    4.945)    8.564
  10.868   (   4.945    4.945    4.945)    8.564
  12.857   (   8.420    8.420    8.420)   14.584
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.707   (   6.522    6.288    6.288)   11.028
   1.707   (   7.878    6.749    6.749)   12.376
   2.189   (  -0.861   -1.743   -1.743)    2.611
   2.789   (  -3.528   11.458   11.458)   16.584
   2.869   (   7.392   14.553   14.553)   21.868
   3.243   (  18.050    0.728    0.728)   18.079
   3.754   (  -3.145   -6.628   -6.628)    9.888
   4.351   (   0.434   -4.432   -4.432)    6.282
   4.355   (   3.338   -3.477   -3.477)    5.943
   5.134   (   3.496    9.633    9.633)   14.065
   5.241   (   8.733    7.617    7.617)   13.867
   7.151   (  -0.217   -2.532   -2.532)    3.588
  10.882   (   0.677    7.271    7.271)   10.305
  11.017   (   5.078    6.396    6.396)   10.373
  13.092   (  13.497    7.634    7.634)   17.284
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.822   (   4.625    4.066    4.066)    7.379
   1.865   (   6.920    4.582    4.582)    9.480
   2.161   (  -1.597   -1.393   -1.393)    2.536
   2.710   (  -3.922   12.739   12.739)   18.437
   3.002   (   5.281   13.382   13.382)   19.648
   3.525   (   9.242   -0.902   -0.902)    9.330
   3.708   (  -1.360   -8.514   -8.514)   12.117
   4.360   (   0.417   -5.941   -5.941)    8.412
   4.463   (   6.119   -3.315   -3.315)    7.708
   5.190   (   1.953    9.445    9.445)   13.499
   5.393   (   5.774    8.012    8.012)   12.718
   7.144   (  -0.389   -5.262   -5.262)    7.451
  10.894   (   0.528    7.417    7.417)   10.503
  11.088   (   2.139    8.788    8.788)   12.610
  13.365   (  11.764    7.541    7.541)   15.879
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.891   (   1.887    1.687    1.687)    3.042
   1.973   (   3.013    2.673    2.673)    4.834
   2.133   (  -0.861    0.377    0.377)    1.013
   2.644   (  -1.985   14.190   14.190)   20.165
   3.078   (   1.887   12.752   12.752)   18.133
   3.642   (   2.576   -2.918   -2.918)    4.864
   3.692   (  -0.311   -9.586   -9.586)   13.560
   4.367   (   0.157   -6.721   -6.721)    9.506
   4.565   (   2.908   -2.157   -2.157)    4.214
   5.215   (   0.579    9.147    9.147)   12.948
   5.474   (   2.009    8.130    8.130)   11.672
   7.138   (  -0.148   -6.855   -6.855)    9.696
  10.902   (   0.199    7.505    7.505)   10.615
  11.114   (   0.584    9.792    9.792)   13.860
  13.539   (   4.551    7.762    7.762)   11.883
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.821   (   5.347    5.347    6.753)   10.139
   1.838   (   5.346    5.346    4.119)    8.610
   2.153   (  -1.559   -1.559   -1.916)    2.920
   2.846   (   0.370    0.370   23.017)   23.023
   3.259   (   7.796    7.796   -2.150)   11.233
   3.311   (  12.598   12.598    5.813)   18.741
   3.624   (  -5.374   -5.374   -5.881)    9.610
   4.312   (   0.915    0.915   -6.133)    6.268
   4.322   (  -0.005   -0.005   -9.347)    9.347
   5.290   (   5.123    5.123   11.825)   13.868
   5.388   (   6.885    6.885    7.987)   12.593
   7.093   (  -2.909   -2.909   -5.714)    7.041
  10.941   (   0.931    0.931   14.638)   14.697
  11.203   (   8.701    8.701    0.442)   12.313
  13.307   (  12.574   12.574    7.100)   19.147
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.902   (   2.488    3.767    4.370)    6.283
   1.948   (   4.999    3.274    3.758)    7.060
   2.131   (  -0.451   -1.105   -0.915)    1.504
   2.833   (  -1.355    2.380   23.218)   23.379
   3.344   (   3.643   11.296    1.176)   11.927
   3.489   (  -0.452   -1.263    0.167)    1.352
   3.586   (   3.534   -1.377   -6.567)    7.583
   4.302   (  -0.584    0.986  -10.318)   10.381
   4.422   (   7.065   -1.021   -6.813)    9.867
   5.357   (   1.848    6.635   11.315)   13.247
   5.513   (   4.920    3.782    9.622)   11.449
   7.040   (  -2.092   -4.458   -8.793)   10.078
  10.953   (   0.348    0.653   15.639)   15.656
  11.339   (   4.536   12.276    0.598)   13.100
  13.568   (  11.538   11.711    6.656)   17.737
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.931   (   0.562    2.075    2.589)    3.365
   2.026   (   2.164    2.304    2.860)    4.263
   2.129   (   0.069   -0.383    0.578)    0.697
   2.806   (  -0.869    4.018   23.511)   23.868
   3.408   (   2.380    8.969    2.661)    9.654
   3.451   (  -1.810   -4.273   -1.065)    4.761
   3.641   (   1.128   -0.776   -9.841)    9.936
   4.293   (  -0.218    0.574  -11.341)   11.357
   4.544   (   3.520   -0.219   -4.781)    5.942
   5.378   (   0.391    6.594   10.713)   12.586
   5.583   (   1.722    2.393   10.493)   10.900
   7.011   (  -0.705   -4.983  -10.674)   11.801
  10.957   (   0.106    0.398   16.037)   16.042
  11.398   (   1.344   14.135    0.707)   14.216
  13.742   (   4.570   11.413    6.515)   13.913
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.961   (   1.437    1.437    4.462)    4.903
   2.009   (   3.085    3.085    2.816)    5.193
   2.122   (   0.062    0.062    0.245)    0.260
   2.862   (   0.755    0.755   25.849)   25.872
   3.326   (  -4.407   -4.407   -2.249)    6.625
   3.541   (   1.885    1.885   -8.440)    8.851
   3.636   (   2.557    2.557   -0.721)    3.688
   4.393   (   3.650    3.650  -10.185)   11.419
   4.422   (   3.752    3.752   -7.834)    9.462
   5.458   (   3.029    3.029   11.980)   12.723
   5.562   (   1.477    1.477   10.851)   11.050
   6.969   (  -2.315   -2.315  -12.115)   12.550
  10.957   (   0.006    0.006   16.804)   16.804
  11.542   (   7.131    7.131   -0.506)   10.097
  13.813   (  10.941   10.941    5.815)   16.530
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.971   (   0.016    1.605    3.751)    4.080
   2.067   (   1.996    1.697    2.632)    3.713
   2.122   (  -0.235   -0.347    1.360)    1.423
   2.870   (   0.158    2.585   27.174)   27.297
   3.261   (  -1.657   -6.653   -2.266)    7.221
   3.583   (   2.257    1.566   -7.067)    7.582
   3.636   (  -1.716    0.060   -7.119)    7.323
   4.401   (  -0.058    8.386   -9.863)   12.946
   4.545   (   3.931    0.214   -6.339)    7.462
   5.487   (   0.306    3.992   11.506)   12.182
   5.587   (   0.788   -1.327   11.733)   11.834
   6.941   (  -0.602   -1.911  -13.931)   14.074
  10.957   (  -0.033   -0.222   17.121)   17.122
  11.636   (   2.214    8.599   -0.849)    8.920
  13.979   (   4.410   10.595    5.373)   12.672
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.992   (   0.483    0.483    4.555)    4.606
   2.096   (   1.107    1.107    2.287)    2.771
   2.114   (  -0.490   -0.490    2.061)    2.175
   2.911   (   1.161    1.161   29.940)   29.985
   3.168   (  -2.267   -2.267   -2.013)    3.785
   3.592   (   0.498    0.498  -11.301)   11.323
   3.615   (  -1.581   -1.581  -10.052)   10.297
   4.545   (   2.467    2.467   -7.180)    7.982
   4.552   (   1.892    1.892   -6.971)    7.467
   5.541   (   0.990    0.990   11.697)   11.780
   5.559   (  -0.670   -0.670   11.692)   11.731
   6.921   (  -0.332   -0.332  -15.421)   15.428
  10.953   (  -0.122   -0.122   17.343)   17.344
  11.752   (   2.732    2.732   -1.386)    4.104
  14.140   (   4.292    4.292    4.740)    7.701
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.932   (   4.467    4.467    4.467)    7.737
   1.932   (   4.467    4.467    4.467)    7.737
   2.117   (  -1.270   -1.270   -1.270)    2.199
   3.211   (   4.182    4.182    4.182)    7.243
   3.211   (   4.182    4.182    4.182)    7.243
   3.501   (  -5.667   -5.667   -5.667)    9.816
   3.505   (  11.612   11.612   11.612)   20.113
   4.211   (  -2.113   -2.113   -2.113)    3.661
   4.211   (  -2.113   -2.113   -2.113)    3.661
   5.509   (   6.601    6.601    6.601)   11.433
   5.509   (   6.601    6.601    6.601)   11.433
   6.966   (  -6.218   -6.218   -6.218)   10.771
  11.213   (   5.733    5.733    5.733)    9.930
  11.213   (   5.733    5.733    5.733)    9.930
  13.505   (  11.595   11.595   11.595)   20.083
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.986   (   1.021    3.966    3.966)    5.702
   2.025   (   3.951    3.288    3.288)    6.102
   2.114   (   1.042   -0.498   -0.498)    1.258
   3.165   (  -1.916    7.543    7.543)   10.838
   3.362   (   5.430   -0.627   -0.627)    5.502
   3.402   (  -3.729   -6.199   -6.199)    9.526
   3.648   (   1.643   10.653   10.653)   15.156
   4.148   (  -2.320   -2.968   -2.968)    4.795
   4.319   (   8.987   -3.530   -3.530)   10.280
   5.559   (   1.116    8.028    8.028)   11.408
   5.666   (   5.823    5.187    5.187)    9.366
   6.852   (  -4.501   -8.747   -8.747)   13.163
  11.223   (   0.382    8.030    8.030)   11.363
  11.383   (   5.651    5.197    5.197)    9.271
  13.749   (  10.929   10.487   10.487)   18.423
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.990   (  -0.107    3.103    3.103)    4.390
   2.084   (   1.562    2.599    2.599)    3.994
   2.140   (   0.916    0.527    0.527)    1.181
   3.137   (  -0.778    7.989    7.989)   11.324
   3.351   (  -1.232   -6.767   -6.767)    9.649
   3.449   (   2.669   -4.477   -4.477)    6.871
   3.635   (  -1.309   11.560   11.560)   16.401
   4.115   (  -0.820   -3.950   -3.950)    5.645
   4.475   (   4.459   -2.238   -2.238)    5.468
   5.567   (   0.037    7.622    7.622)   10.779
   5.747   (   1.980    5.089    5.089)    7.464
   6.789   (  -1.561   -9.940   -9.940)   14.143
  11.228   (   0.144    7.836    7.836)   11.082
  11.456   (   1.648    6.283    6.283)    9.037
  13.914   (   4.402    9.915    9.915)   14.697
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.051   (   1.052    1.052    4.230)    4.484
   2.074   (   2.636    2.636    3.334)    5.002
   2.127   (   1.636    1.636    0.259)    2.329
   3.235   (   0.991    0.991    6.944)    7.084
   3.271   (  -4.546   -4.546   -3.111)    7.142
   3.370   (   0.252    0.252   -7.484)    7.492
   3.755   (  -1.123   -1.123   14.536)   14.622
   4.229   (   4.808    4.808   -5.746)    8.902
   4.286   (   3.986    3.986   -5.479)    7.861
   5.675   (   2.850    2.850    8.965)    9.829
   5.742   (   2.911    2.911    6.430)    7.635
   6.702   (  -5.443   -5.443  -12.968)   15.080
  11.288   (   0.784    0.784   14.213)   14.256
  11.530   (   6.544    6.544   -0.550)    9.271
  13.970   (  10.025   10.025    9.176)   16.888
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.052   (  -0.166    2.648    3.995)    4.796
   2.126   (   1.671    1.422    2.915)    3.648
   2.153   (   0.563    0.713    1.573)    1.816
   3.206   (  -1.520   -5.797   -2.891)    6.654
   3.241   (  -0.124    3.776    4.128)    5.596
   3.386   (   0.809   -2.188   -9.623)    9.901
   3.706   (  -1.962   -4.387   15.968)   16.675
   4.230   (  -0.510   11.899   -6.833)   13.731
   4.442   (   4.999   -1.075   -3.922)    6.445
   5.692   (  -0.163    4.486    8.305)    9.441
   5.791   (   1.398   -0.392    7.833)    7.966
   6.631   (  -1.607   -5.061  -15.200)   16.101
  11.295   (   0.096    0.771   14.498)   14.519
  11.620   (   2.120    7.298   -0.613)    7.624
  14.124   (   4.100    9.466    8.417)   13.314
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.090   (   0.867    0.867    4.636)    4.796
   2.149   (   0.966    0.966    2.734)    3.056
   2.162   (  -0.049   -0.049    2.431)    2.432
   3.123   (  -2.053   -2.053   -2.268)    3.684
   3.316   (   2.637    2.637   -3.764)    5.299
   3.357   (  -0.397   -0.397  -10.842)   10.857
   3.595   (  -4.311   -4.311   20.600)   21.483
   4.423   (   2.900    2.900   -4.868)    6.365
   4.431   (   2.145    2.145   -4.901)    5.765
   5.751   (   0.874    0.874    8.579)    8.668
   5.767   (  -0.598   -0.598    8.463)    8.505
   6.572   (  -1.122   -1.122  -17.506)   17.578
  11.300   (   0.024    0.024   15.258)   15.258
  11.722   (   2.451    2.451   -1.512)    3.782
  14.269   (   3.893    3.893    7.537)    9.334
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.123   (   1.394    1.394    1.394)    2.414
   2.138   (   2.362    2.362    2.362)    4.090
   2.138   (   2.362    2.362    2.362)    4.090
   3.196   (  -3.783   -3.783   -3.783)    6.552
   3.256   (  -2.480   -2.480   -2.480)    4.296
   3.256   (  -2.480   -2.480   -2.480)    4.296
   4.057   (   6.201    6.201    6.201)   10.741
   4.192   (   1.535    1.535    1.535)    2.659
   4.192   (   1.535    1.535    1.535)    2.659
   5.825   (   3.572    3.572    3.572)    6.187
   5.825   (   3.572    3.572    3.572)    6.187
   6.471   (  -8.672   -8.672   -8.672)   15.021
  11.520   (   3.861    3.861    3.861)    6.687
  11.520   (   3.861    3.861    3.861)    6.687
  14.166   (   8.934    8.934    8.934)   15.475
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.129   (   0.018    3.209    3.209)    4.539
   2.178   (   1.037    1.873    1.873)    2.844
   2.183   (   1.278    1.376    1.376)    2.328
   3.139   (  -1.462   -3.493   -3.493)    5.152
   3.212   (  -1.033   -3.098   -3.098)    4.501
   3.229   (  -0.791   -5.161   -5.161)    7.342
   4.063   (  -0.735    8.247    8.247)   11.686
   4.140   (  -1.319    4.155    4.155)    6.023
   4.377   (   5.657   -2.391   -2.391)    6.591
   5.829   (  -0.506    4.912    4.912)    6.965
   5.900   (   1.829    2.478    2.478)    3.953
   6.358   (  -2.452  -10.086  -10.086)   14.473
  11.514   (  -0.156    5.677    5.677)    8.030
  11.628   (   2.498    2.298    2.298)    4.099
  14.304   (   3.706    8.250    8.250)   12.242
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.175   (   1.065    1.065    3.413)    3.731
   2.201   (   0.813    0.813    2.144)    2.433
   2.207   (   0.486    0.486    1.794)    1.921
   3.079   (  -1.713   -1.713   -1.825)    3.032
   3.165   (  -1.172   -1.172   -6.790)    6.990
   3.167   (  -1.310   -1.310   -7.208)    7.443
   4.040   (  -1.512   -1.512   18.100)   18.226
   4.340   (   2.903    2.903   -3.164)    5.184
   4.346   (   2.301    2.301   -3.275)    4.617
   5.892   (   0.763    0.763    5.140)    5.253
   5.906   (  -0.213   -0.213    4.890)    4.899
   6.238   (  -2.123   -2.123  -13.848)   14.170
  11.554   (   0.241    0.241    9.110)    9.117
  11.693   (   2.193    2.193   -1.158)    3.310
  14.431   (   3.447    3.447    7.496)    8.942
======================= 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 0.00e+00 0.00e+00 0.00e+00 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.224   (   1.200    1.200    1.200)    2.078
   2.232   (   0.709    0.709    0.709)    1.228
   2.232   (   0.709    0.709    0.709)    1.228
   3.050   (  -1.059   -1.059   -1.059)    1.834
   3.069   (  -2.208   -2.208   -2.208)    3.824
   3.069   (  -2.208   -2.208   -2.208)    3.824
   4.292   (   1.695    1.695    1.695)    2.936
   4.299   (   1.185    1.185    1.185)    2.053
   4.299   (   1.185    1.185    1.185)    2.053
   5.963   (   1.031    1.031    1.031)    1.786
   5.963   (   1.031    1.031    1.031)    1.786
   6.043   (  -4.071   -4.071   -4.071)    7.051
  11.676   (   1.215    1.215    1.215)    2.105
  11.676   (   1.215    1.215    1.215)    2.105
  14.547   (   3.162    3.162    3.162)    5.477
=================== 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    807.963    807.963    807.963      0.000      0.000      0.000 3/19965
   20.0    117.270    117.270    117.270     -0.000     -0.000      0.000 3/19965
   30.0     42.509     42.509     42.509     -0.000     -0.000      0.000 3/19965
   40.0     24.064     24.064     24.064     -0.000     -0.000      0.000 3/19965
   50.0     16.690     16.690     16.690     -0.000     -0.000      0.000 3/19965
   60.0     12.848     12.848     12.848     -0.000     -0.000      0.000 3/19965
   70.0     10.497     10.497     10.497     -0.000     -0.000      0.000 3/19965
   80.0      8.898      8.898      8.898     -0.000     -0.000      0.000 3/19965
   90.0      7.732      7.732      7.732     -0.000     -0.000      0.000 3/19965
  100.0      6.838      6.838      6.838     -0.000     -0.000      0.000 3/19965
  110.0      6.129      6.129      6.129     -0.000     -0.000      0.000 3/19965
  120.0      5.552      5.552      5.552     -0.000     -0.000      0.000 3/19965
  130.0      5.074      5.074      5.074     -0.000     -0.000      0.000 3/19965
  140.0      4.671      4.671      4.671     -0.000     -0.000      0.000 3/19965
  150.0      4.327      4.327      4.327     -0.000     -0.000      0.000 3/19965
  160.0      4.030      4.030      4.030     -0.000     -0.000      0.000 3/19965
  170.0      3.771      3.771      3.771     -0.000     -0.000      0.000 3/19965
  180.0      3.543      3.543      3.543     -0.000     -0.000      0.000 3/19965
  190.0      3.342      3.342      3.342     -0.000     -0.000      0.000 3/19965
  200.0      3.162      3.162      3.162     -0.000     -0.000      0.000 3/19965
  210.0      3.001      3.001      3.001     -0.000     -0.000      0.000 3/19965
  220.0      2.856      2.856      2.856     -0.000     -0.000      0.000 3/19965
  230.0      2.724      2.724      2.724     -0.000     -0.000      0.000 3/19965
  240.0      2.604      2.604      2.604     -0.000     -0.000      0.000 3/19965
  250.0      2.494      2.494      2.494     -0.000     -0.000      0.000 3/19965
  260.0      2.393      2.393      2.393     -0.000     -0.000      0.000 3/19965
  270.0      2.300      2.300      2.300     -0.000     -0.000      0.000 3/19965
  280.0      2.215      2.215      2.215     -0.000     -0.000      0.000 3/19965
  290.0      2.135      2.135      2.135     -0.000     -0.000      0.000 3/19965
  300.0      2.061      2.061      2.061     -0.000     -0.000      0.000 3/19965
  310.0      1.992      1.992      1.992     -0.000     -0.000      0.000 3/19965
  320.0      1.927      1.927      1.927     -0.000     -0.000      0.000 3/19965
  330.0      1.867      1.867      1.867     -0.000     -0.000      0.000 3/19965
  340.0      1.810      1.810      1.810     -0.000     -0.000      0.000 3/19965
  350.0      1.757      1.757      1.757     -0.000     -0.000      0.000 3/19965
  360.0      1.707      1.707      1.707     -0.000     -0.000      0.000 3/19965
  370.0      1.660      1.660      1.660     -0.000     -0.000      0.000 3/19965
  380.0      1.615      1.615      1.615     -0.000     -0.000      0.000 3/19965
  390.0      1.572      1.572      1.572     -0.000     -0.000      0.000 3/19965
  400.0      1.532      1.532      1.532     -0.000     -0.000      0.000 3/19965
  410.0      1.494      1.494      1.494     -0.000     -0.000      0.000 3/19965
  420.0      1.458      1.458      1.458     -0.000     -0.000      0.000 3/19965
  430.0      1.423      1.423      1.423     -0.000     -0.000      0.000 3/19965
  440.0      1.390      1.390      1.390     -0.000     -0.000      0.000 3/19965
  450.0      1.358      1.358      1.358     -0.000     -0.000      0.000 3/19965
  460.0      1.328      1.328      1.328     -0.000     -0.000      0.000 3/19965
  470.0      1.300      1.300      1.300     -0.000     -0.000      0.000 3/19965
  480.0      1.272      1.272      1.272     -0.000     -0.000      0.000 3/19965
  490.0      1.246      1.246      1.246     -0.000     -0.000      0.000 3/19965
  500.0      1.220      1.220      1.220     -0.000     -0.000      0.000 3/19965
  510.0      1.196      1.196      1.196     -0.000     -0.000      0.000 3/19965
  520.0      1.173      1.173      1.173     -0.000     -0.000      0.000 3/19965
  530.0      1.150      1.150      1.150     -0.000     -0.000      0.000 3/19965
  540.0      1.129      1.129      1.129     -0.000     -0.000      0.000 3/19965
  550.0      1.108      1.108      1.108     -0.000     -0.000      0.000 3/19965
  560.0      1.088      1.088      1.088     -0.000     -0.000      0.000 3/19965
  570.0      1.069      1.069      1.069     -0.000     -0.000      0.000 3/19965
  580.0      1.050      1.050      1.050     -0.000     -0.000      0.000 3/19965
  590.0      1.032      1.032      1.032     -0.000     -0.000      0.000 3/19965
  600.0      1.015      1.015      1.015     -0.000     -0.000      0.000 3/19965
  610.0      0.998      0.998      0.998     -0.000     -0.000      0.000 3/19965
  620.0      0.982      0.982      0.982     -0.000     -0.000      0.000 3/19965
  630.0      0.966      0.966      0.966     -0.000     -0.000      0.000 3/19965
  640.0      0.951      0.951      0.951     -0.000     -0.000      0.000 3/19965
  650.0      0.936      0.936      0.936     -0.000     -0.000      0.000 3/19965
  660.0      0.922      0.922      0.922     -0.000     -0.000      0.000 3/19965
  670.0      0.908      0.908      0.908     -0.000     -0.000      0.000 3/19965
  680.0      0.894      0.894      0.894     -0.000     -0.000      0.000 3/19965
  690.0      0.881      0.881      0.881     -0.000     -0.000      0.000 3/19965
  700.0      0.869      0.869      0.869     -0.000     -0.000      0.000 3/19965
  710.0      0.856      0.856      0.856     -0.000     -0.000      0.000 3/19965
  720.0      0.844      0.844      0.844     -0.000     -0.000      0.000 3/19965
  730.0      0.833      0.833      0.833     -0.000     -0.000      0.000 3/19965
  740.0      0.821      0.821      0.821     -0.000     -0.000      0.000 3/19965
  750.0      0.810      0.810      0.810     -0.000     -0.000      0.000 3/19965
  760.0      0.799      0.799      0.799     -0.000     -0.000      0.000 3/19965
  770.0      0.789      0.789      0.789     -0.000     -0.000      0.000 3/19965
  780.0      0.779      0.779      0.779     -0.000     -0.000      0.000 3/19965
  790.0      0.769      0.769      0.769     -0.000     -0.000      0.000 3/19965
  800.0      0.759      0.759      0.759     -0.000     -0.000      0.000 3/19965
  810.0      0.750      0.750      0.750     -0.000     -0.000      0.000 3/19965
  820.0      0.741      0.741      0.741     -0.000     -0.000      0.000 3/19965
  830.0      0.732      0.732      0.732     -0.000     -0.000      0.000 3/19965
  840.0      0.723      0.723      0.723     -0.000     -0.000      0.000 3/19965
  850.0      0.714      0.714      0.714     -0.000     -0.000      0.000 3/19965
  860.0      0.706      0.706      0.706     -0.000     -0.000      0.000 3/19965
  870.0      0.698      0.698      0.698     -0.000     -0.000      0.000 3/19965
  880.0      0.690      0.690      0.690     -0.000     -0.000      0.000 3/19965
  890.0      0.682      0.682      0.682     -0.000     -0.000      0.000 3/19965
  900.0      0.674      0.674      0.674     -0.000     -0.000      0.000 3/19965
  910.0      0.667      0.667      0.667     -0.000     -0.000      0.000 3/19965
  920.0      0.660      0.660      0.660     -0.000     -0.000      0.000 3/19965
  930.0      0.653      0.653      0.653     -0.000     -0.000      0.000 3/19965
  940.0      0.646      0.646      0.646     -0.000     -0.000      0.000 3/19965
  950.0      0.639      0.639      0.639     -0.000     -0.000      0.000 3/19965
  960.0      0.632      0.632      0.632     -0.000     -0.000      0.000 3/19965
  970.0      0.626      0.626      0.626     -0.000     -0.000      0.000 3/19965
  980.0      0.619      0.619      0.619     -0.000     -0.000      0.000 3/19965
  990.0      0.613      0.613      0.613     -0.000     -0.000      0.000 3/19965
 1000.0      0.607      0.607      0.607     -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-07 22:53:17]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

