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

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

-------------------------[time 2026-01-07 22:08:34]-------------------------
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.356899010000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356899010000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356899010000000
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468
   *5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.356899010000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356899010000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356899010000000
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468 > 4
   *5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   13.070697030000002    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.070697030000002    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.070697030000002
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468 > 4
   83 Rb  0.33333333333333  0.00000000000000  0.00000000000000  85.468 > 4
   84 Rb  0.66666666666667  0.00000000000000  0.00000000000000  85.468 > 4
   85 Rb  0.00000000000000  0.33333333333333  0.00000000000000  85.468 > 4
   86 Rb  0.33333333333333  0.33333333333333  0.00000000000000  85.468 > 4
   87 Rb  0.66666666666667  0.33333333333333  0.00000000000000  85.468 > 4
   88 Rb  0.00000000000000  0.66666666666667  0.00000000000000  85.468 > 4
   89 Rb  0.33333333333333  0.66666666666667  0.00000000000000  85.468 > 4
   90 Rb  0.66666666666667  0.66666666666667  0.00000000000000  85.468 > 4
   91 Rb  0.00000000000000  0.00000000000000  0.33333333333333  85.468 > 4
   92 Rb  0.33333333333333  0.00000000000000  0.33333333333333  85.468 > 4
   93 Rb  0.66666666666667  0.00000000000000  0.33333333333333  85.468 > 4
   94 Rb  0.00000000000000  0.33333333333333  0.33333333333333  85.468 > 4
   95 Rb  0.33333333333333  0.33333333333333  0.33333333333333  85.468 > 4
   96 Rb  0.66666666666667  0.33333333333333  0.33333333333333  85.468 > 4
   97 Rb  0.00000000000000  0.66666666666667  0.33333333333333  85.468 > 4
   98 Rb  0.33333333333333  0.66666666666667  0.33333333333333  85.468 > 4
   99 Rb  0.66666666666667  0.66666666666667  0.33333333333333  85.468 > 4
  100 Rb  0.00000000000000  0.00000000000000  0.66666666666667  85.468 > 4
  101 Rb  0.33333333333333  0.00000000000000  0.66666666666667  85.468 > 4
  102 Rb  0.66666666666667  0.00000000000000  0.66666666666667  85.468 > 4
  103 Rb  0.00000000000000  0.33333333333333  0.66666666666667  85.468 > 4
  104 Rb  0.33333333333333  0.33333333333333  0.66666666666667  85.468 > 4
  105 Rb  0.66666666666667  0.33333333333333  0.66666666666667  85.468 > 4
  106 Rb  0.00000000000000  0.66666666666667  0.66666666666667  85.468 > 4
  107 Rb  0.33333333333333  0.66666666666667  0.66666666666667  85.468 > 4
  108 Rb  0.66666666666667  0.66666666666667  0.66666666666667  85.468 > 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.2563801    0.0000000    0.0000000
            0.0000000    2.2563801    0.0000000
            0.0000000    0.0000000    2.2563801
-------------------------- Born effective charges --------------------------
    1 F    -1.7792314    0.0000000    0.0000000
            0.0000000   -0.9484423    0.0000000
            0.0000000    0.0000000   -0.9484423
    2 F    -0.9484423    0.0000000    0.0000000
            0.0000000   -1.7792314    0.0000000
            0.0000000    0.0000000   -0.9484423
    3 F    -0.9484423    0.0000000    0.0000000
            0.0000000   -0.9484423    0.0000000
            0.0000000    0.0000000   -1.7792314
    4 Rb    1.2332581    0.0000000    0.0000000
            0.0000000    1.2332581    0.0000000
            0.0000000    0.0000000    1.2332581
    5 Yb    2.4428579    0.0000000    0.0000000
            0.0000000    2.4428579    0.0000000
            0.0000000    0.0000000    2.4428579
----------------------------------------------------------------------------
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.073
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.077
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) 
Permutation basis: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (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:08:39]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 22:08:39]-------------------------
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.356899010000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356899010000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356899010000000
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468
    5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.070697030000002    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.070697030000002    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.070697030000002
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468 > 82
   83 Rb  0.33333333333333  0.00000000000000  0.00000000000000  85.468 > 82
   84 Rb  0.66666666666667  0.00000000000000  0.00000000000000  85.468 > 82
   85 Rb  0.00000000000000  0.33333333333333  0.00000000000000  85.468 > 82
   86 Rb  0.33333333333333  0.33333333333333  0.00000000000000  85.468 > 82
   87 Rb  0.66666666666667  0.33333333333333  0.00000000000000  85.468 > 82
   88 Rb  0.00000000000000  0.66666666666667  0.00000000000000  85.468 > 82
   89 Rb  0.33333333333333  0.66666666666667  0.00000000000000  85.468 > 82
   90 Rb  0.66666666666667  0.66666666666667  0.00000000000000  85.468 > 82
   91 Rb  0.00000000000000  0.00000000000000  0.33333333333333  85.468 > 82
   92 Rb  0.33333333333333  0.00000000000000  0.33333333333333  85.468 > 82
   93 Rb  0.66666666666667  0.00000000000000  0.33333333333333  85.468 > 82
   94 Rb  0.00000000000000  0.33333333333333  0.33333333333333  85.468 > 82
   95 Rb  0.33333333333333  0.33333333333333  0.33333333333333  85.468 > 82
   96 Rb  0.66666666666667  0.33333333333333  0.33333333333333  85.468 > 82
   97 Rb  0.00000000000000  0.66666666666667  0.33333333333333  85.468 > 82
   98 Rb  0.33333333333333  0.66666666666667  0.33333333333333  85.468 > 82
   99 Rb  0.66666666666667  0.66666666666667  0.33333333333333  85.468 > 82
  100 Rb  0.00000000000000  0.00000000000000  0.66666666666667  85.468 > 82
  101 Rb  0.33333333333333  0.00000000000000  0.66666666666667  85.468 > 82
  102 Rb  0.66666666666667  0.00000000000000  0.66666666666667  85.468 > 82
  103 Rb  0.00000000000000  0.33333333333333  0.66666666666667  85.468 > 82
  104 Rb  0.33333333333333  0.33333333333333  0.66666666666667  85.468 > 82
  105 Rb  0.66666666666667  0.33333333333333  0.66666666666667  85.468 > 82
  106 Rb  0.00000000000000  0.66666666666667  0.66666666666667  85.468 > 82
  107 Rb  0.33333333333333  0.66666666666667  0.66666666666667  85.468 > 82
  108 Rb  0.66666666666667  0.66666666666667  0.66666666666667  85.468 > 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.2563801    0.0000000    0.0000000
            0.0000000    2.2563801    0.0000000
            0.0000000    0.0000000    2.2563801
-------------------------- Born effective charges --------------------------
    1 F    -1.7792314    0.0000000    0.0000000
            0.0000000   -0.9484423    0.0000000
            0.0000000    0.0000000   -0.9484423
    2 F    -0.9484423    0.0000000    0.0000000
            0.0000000   -1.7792314    0.0000000
            0.0000000    0.0000000   -0.9484423
    3 F    -0.9484423    0.0000000    0.0000000
            0.0000000   -0.9484423    0.0000000
            0.0000000    0.0000000   -1.7792314
    4 Rb    1.2332581    0.0000000    0.0000000
            0.0000000    1.2332581    0.0000000
            0.0000000    0.0000000    1.2332581
    5 Yb    2.4428579    0.0000000    0.0000000
            0.0000000    2.4428579    0.0000000
            0.0000000    0.0000000    2.4428579
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  0.0000]
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 82, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 109, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (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:08:43]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 22:08:43]-------------------------
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.356899010000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.356899010000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.356899010000000
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468
    5 Yb  0.50000000000000  0.50000000000000  0.50000000000000 173.054
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.070697030000002    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.070697030000002    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.070697030000002
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 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468 > 82
   83 Rb  0.33333333333333  0.00000000000000  0.00000000000000  85.468 > 82
   84 Rb  0.66666666666667  0.00000000000000  0.00000000000000  85.468 > 82
   85 Rb  0.00000000000000  0.33333333333333  0.00000000000000  85.468 > 82
   86 Rb  0.33333333333333  0.33333333333333  0.00000000000000  85.468 > 82
   87 Rb  0.66666666666667  0.33333333333333  0.00000000000000  85.468 > 82
   88 Rb  0.00000000000000  0.66666666666667  0.00000000000000  85.468 > 82
   89 Rb  0.33333333333333  0.66666666666667  0.00000000000000  85.468 > 82
   90 Rb  0.66666666666667  0.66666666666667  0.00000000000000  85.468 > 82
   91 Rb  0.00000000000000  0.00000000000000  0.33333333333333  85.468 > 82
   92 Rb  0.33333333333333  0.00000000000000  0.33333333333333  85.468 > 82
   93 Rb  0.66666666666667  0.00000000000000  0.33333333333333  85.468 > 82
   94 Rb  0.00000000000000  0.33333333333333  0.33333333333333  85.468 > 82
   95 Rb  0.33333333333333  0.33333333333333  0.33333333333333  85.468 > 82
   96 Rb  0.66666666666667  0.33333333333333  0.33333333333333  85.468 > 82
   97 Rb  0.00000000000000  0.66666666666667  0.33333333333333  85.468 > 82
   98 Rb  0.33333333333333  0.66666666666667  0.33333333333333  85.468 > 82
   99 Rb  0.66666666666667  0.66666666666667  0.33333333333333  85.468 > 82
  100 Rb  0.00000000000000  0.00000000000000  0.66666666666667  85.468 > 82
  101 Rb  0.33333333333333  0.00000000000000  0.66666666666667  85.468 > 82
  102 Rb  0.66666666666667  0.00000000000000  0.66666666666667  85.468 > 82
  103 Rb  0.00000000000000  0.33333333333333  0.66666666666667  85.468 > 82
  104 Rb  0.33333333333333  0.33333333333333  0.66666666666667  85.468 > 82
  105 Rb  0.66666666666667  0.33333333333333  0.66666666666667  85.468 > 82
  106 Rb  0.00000000000000  0.66666666666667  0.66666666666667  85.468 > 82
  107 Rb  0.33333333333333  0.66666666666667  0.66666666666667  85.468 > 82
  108 Rb  0.66666666666667  0.66666666666667  0.66666666666667  85.468 > 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.2563801    0.0000000    0.0000000
            0.0000000    2.2563801    0.0000000
            0.0000000    0.0000000    2.2563801
-------------------------- Born effective charges --------------------------
    1 F    -1.7792314    0.0000000    0.0000000
            0.0000000   -0.9484423    0.0000000
            0.0000000    0.0000000   -0.9484423
    2 F    -0.9484423    0.0000000    0.0000000
            0.0000000   -1.7792314    0.0000000
            0.0000000    0.0000000   -0.9484423
    3 F    -0.9484423    0.0000000    0.0000000
            0.0000000   -0.9484423    0.0000000
            0.0000000    0.0000000   -1.7792314
    4 Rb    1.2332581    0.0000000    0.0000000
            0.0000000    1.2332581    0.0000000
            0.0000000    0.0000000    1.2332581
    5 Yb    2.4428579    0.0000000    0.0000000
            0.0000000    2.4428579    0.0000000
            0.0000000    0.0000000    2.4428579
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 11 11 11 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.95, Number of G-points: 305, Lambda: 0.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 1.10e-04 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.158   (   0.000    0.000    0.000)    0.000
   2.158   (   0.000    0.000    0.000)    0.000
   2.158   (   0.000    0.000    0.000)    0.000
   3.203   (   0.000    0.000    0.000)    0.000
   3.203   (   0.000    0.000    0.000)    0.000
   3.203   (   0.000    0.000    0.000)    0.000
   4.253   (   0.000    0.000    0.000)    0.000
   4.253   (   0.000    0.000    0.000)    0.000
   4.253   (   0.000    0.000    0.000)    0.000
  11.660   (   0.000    0.000    0.000)    0.000
  11.660   (   0.000    0.000    0.000)    0.000
  11.660   (   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 1.10e-04 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.395   (  18.389    0.000    0.000)   18.389
   0.395   (  18.389    0.000    0.000)   18.389
   0.986   (  44.256    0.000    0.000)   44.256
   2.148   (  -0.928    0.000    0.000)    0.928
   2.148   (  -0.928    0.000    0.000)    0.928
   2.729   (   8.847    0.000    0.000)    8.847
   3.239   (   3.400    0.000    0.000)    3.400
   3.252   (   4.595    0.000    0.000)    4.595
   3.252   (   4.595    0.000    0.000)    4.595
   4.251   (  -0.133    0.000    0.000)    0.133
   4.251   (  -0.133    0.000    0.000)    0.133
   6.405   (   4.889    0.000    0.000)    4.889
  11.658   (  -0.171    0.000    0.000)    0.171
  11.658   (  -0.171    0.000    0.000)    0.171
  13.681   (   4.971    0.000    0.000)    4.971
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.759   (  16.339    0.000    0.000)   16.339
   0.759   (  16.339    0.000    0.000)   16.339
   1.798   (  32.903    0.000    0.000)   32.903
   2.116   (  -2.229    0.000    0.000)    2.229
   2.116   (  -2.229    0.000    0.000)    2.229
   2.987   (  15.168    0.000    0.000)   15.168
   3.336   (   5.542    0.000    0.000)    5.542
   3.387   (   8.134    0.000    0.000)    8.134
   3.387   (   8.134    0.000    0.000)    8.134
   4.245   (  -0.495    0.000    0.000)    0.495
   4.245   (  -0.495    0.000    0.000)    0.495
   6.548   (   8.570    0.000    0.000)    8.570
  11.653   (  -0.289    0.000    0.000)    0.289
  11.653   (  -0.289    0.000    0.000)    0.289
  13.827   (   8.712    0.000    0.000)    8.712
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.073   (  13.794    0.000    0.000)   13.794
   1.073   (  13.794    0.000    0.000)   13.794
   2.054   (  -3.801    0.000    0.000)    3.801
   2.054   (  -3.801    0.000    0.000)    3.801
   2.358   (  21.265    0.000    0.000)   21.265
   3.317   (  15.233    0.000    0.000)   15.233
   3.457   (   5.801    0.000    0.000)    5.801
   3.577   (   9.686    0.000    0.000)    9.686
   3.577   (   9.686    0.000    0.000)    9.686
   4.228   (  -1.295    0.000    0.000)    1.295
   4.228   (  -1.295    0.000    0.000)    1.295
   6.744   (   9.696    0.000    0.000)    9.696
  11.647   (  -0.316    0.000    0.000)    0.316
  11.647   (  -0.316    0.000    0.000)    0.316
  14.027   (   9.953    0.000    0.000)    9.953
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.333   (  10.995    0.000    0.000)   10.995
   1.333   (  10.995    0.000    0.000)   10.995
   1.959   (  -5.139    0.000    0.000)    5.139
   1.959   (  -5.139    0.000    0.000)    5.139
   2.707   (  12.476    0.000    0.000)   12.476
   3.565   (   4.284    0.000    0.000)    4.284
   3.582   (   9.706    0.000    0.000)    9.706
   3.776   (   8.994    0.000    0.000)    8.994
   3.776   (   8.994    0.000    0.000)    8.994
   4.188   (  -2.589    0.000    0.000)    2.589
   4.188   (  -2.589    0.000    0.000)    2.589
   6.930   (   7.545    0.000    0.000)    7.545
  11.641   (  -0.243    0.000    0.000)    0.243
  11.641   (  -0.243    0.000    0.000)    0.243
  14.219   (   7.909    0.000    0.000)    7.909
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.519   (   5.947    0.000    0.000)    5.947
   1.519   (   5.947    0.000    0.000)    5.947
   1.854   (  -4.027    0.000    0.000)    4.027
   1.854   (  -4.027    0.000    0.000)    4.027
   2.882   (   4.263    0.000    0.000)    4.263
   3.628   (   1.567    0.000    0.000)    1.567
   3.716   (   3.119    0.000    0.000)    3.119
   3.935   (   5.527    0.000    0.000)    5.527
   3.935   (   5.527    0.000    0.000)    5.527
   4.123   (  -3.145    0.000    0.000)    3.145
   4.123   (  -3.145    0.000    0.000)    3.145
   7.041   (   2.822    0.000    0.000)    2.822
  11.637   (  -0.091    0.000    0.000)    0.091
  11.637   (  -0.091    0.000    0.000)    0.091
  14.337   (   3.013    0.000    0.000)    3.013
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.559   (  13.053   13.053    0.000)   18.459
   0.895   (  19.249   19.249    0.000)   27.223
   1.189   (  26.571   26.571    0.000)   37.578
   2.140   (  -0.821   -0.821    0.000)    1.161
   2.257   (   4.649    4.649    0.000)    6.575
   2.652   (   0.830    0.830    0.000)    1.173
   3.157   (  -2.509   -2.509    0.000)    3.548
   3.298   (   4.372    4.372    0.000)    6.183
   3.404   (   9.078    9.078    0.000)   12.838
   4.251   (  -0.079   -0.079    0.000)    0.112
   4.289   (   2.053    2.053    0.000)    2.903
   6.450   (   4.228    4.228    0.000)    5.979
  11.656   (  -0.175   -0.175    0.000)    0.247
  11.726   (   3.227    3.227    0.000)    4.563
  13.687   (   3.038    3.038    0.000)    4.297
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.857   (  14.490    8.575    0.000)   16.837
   1.149   (  10.506   25.664    0.000)   27.732
   1.848   (  28.916    6.484    0.000)   29.634
   2.111   (  -2.044   -0.520    0.000)    2.110
   2.253   (  -2.986   11.785    0.000)   12.158
   2.855   (  15.394   -9.690    0.000)   18.190
   3.144   (   1.385   -9.353    0.000)    9.455
   3.427   (   7.780    3.793    0.000)    8.656
   3.611   (  10.371   12.874    0.000)   16.532
   4.246   (  -0.391    0.078    0.000)    0.399
   4.327   (   1.749    5.982    0.000)    6.232
   6.574   (   7.468    2.487    0.000)    7.871
  11.651   (  -0.295   -0.185    0.000)    0.349
  11.764   (   0.811    8.633    0.000)    8.671
  13.821   (   9.093    0.733    0.000)    9.123
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.145   (  12.906    6.362    0.000)   14.389
   1.370   (  10.530   20.735    0.000)   23.256
   2.053   (  -3.602   -0.109    0.000)    3.604
   2.166   (  -3.811   11.667    0.000)   12.274
   2.342   (  17.439   -1.272    0.000)   17.486
   3.123   (   6.059  -18.185    0.000)   19.167
   3.271   (  12.993   -5.942    0.000)   14.287
   3.610   (   9.315    3.068    0.000)    9.807
   3.822   (   9.502   14.600    0.000)   17.420
   4.231   (  -1.150    0.326    0.000)    1.195
   4.366   (   1.966   10.254    0.000)   10.440
   6.747   (   8.653    0.353    0.000)    8.660
  11.645   (  -0.323   -0.199    0.000)    0.380
  11.772   (   0.164   10.937    0.000)   10.939
  14.036   (  10.853    1.430    0.000)   10.947
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.391   (  10.599    5.173    0.000)   11.794
   1.582   (   9.685   16.677    0.000)   19.285
   1.962   (  -5.056    0.246    0.000)    5.062
   2.065   (  -5.581   12.127    0.000)   13.350
   2.629   (  10.224   -6.281    0.000)   12.000
   3.198   (   2.440  -24.773    0.000)   24.893
   3.536   (  10.297   -2.803    0.000)   10.671
   3.801   (   8.677    2.381    0.000)    8.998
   3.993   (   6.649   13.887    0.000)   15.397
   4.195   (  -2.414    0.648    0.000)    2.500
   4.405   (   1.621   15.213    0.000)   15.299
   6.915   (   6.914   -1.302    0.000)    7.036
  11.638   (  -0.249   -0.213    0.000)    0.327
  11.774   (   0.036   11.993    0.000)   11.993
  14.245   (   8.602    2.772    0.000)    9.038
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.575   (   6.250    5.034    0.000)    8.025
   1.771   (   8.339   12.660    0.000)   15.159
   1.854   (  -4.424    0.033    0.000)    4.424
   1.934   (  -6.936   14.677    0.000)   16.233
   2.773   (   3.493   -8.832    0.000)    9.497
   3.232   (   0.805  -28.002    0.000)   28.014
   3.680   (   3.443   -1.468    0.000)    3.743
   3.955   (   5.349    1.864    0.000)    5.664
   4.089   (   2.390   11.904    0.000)   12.142
   4.133   (  -3.015    0.978    0.000)    3.170
   4.428   (   0.593   19.164    0.000)   19.174
   7.018   (   2.629   -2.099    0.000)    3.364
  11.635   (  -0.093   -0.221    0.000)    0.240
  11.774   (   0.009   12.475    0.000)   12.475
  14.373   (   3.260    3.614    0.000)    4.867
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.074   (  11.449   11.449    0.000)   16.191
   1.528   (  10.653   10.653    0.000)   15.066
   2.091   (  -1.471   -1.471    0.000)    2.080
   2.092   (  14.004   14.004    0.000)   19.805
   2.499   (   4.855    4.855    0.000)    6.865
   2.750   (   5.742    5.742    0.000)    8.121
   2.991   (  -4.089   -4.089    0.000)    5.783
   3.540   (   6.828    6.828    0.000)    9.656
   3.888   (  13.238   13.238    0.000)   18.721
   4.247   (  -0.089   -0.089    0.000)    0.127
   4.451   (   5.790    5.790    0.000)    8.188
   6.650   (   4.701    4.701    0.000)    6.648
  11.646   (  -0.313   -0.313    0.000)    0.443
  11.934   (   6.718    6.718    0.000)    9.501
  13.905   (   7.594    7.594    0.000)   10.739
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.311   (  10.945    8.955    0.000)   14.141
   1.691   (   6.248    9.194    0.000)   11.116
   2.047   (  -2.880   -0.542    0.000)    2.931
   2.304   (   4.901   -1.220    0.000)    5.051
   2.503   (   1.833    8.790    0.000)    8.979
   2.798   (  -4.193   -3.038    0.000)    5.178
   3.175   (  17.315   -3.432    0.000)   17.652
   3.701   (   8.282    5.571    0.000)    9.981
   4.141   (  10.507   15.632    0.000)   18.835
   4.239   (  -0.724    0.389    0.000)    0.822
   4.575   (   6.220    9.222    0.000)   11.123
   6.764   (   5.945    1.358    0.000)    6.098
  11.639   (  -0.344   -0.338    0.000)    0.482
  12.037   (   3.400   12.900    0.000)   13.341
  14.115   (  11.618    6.650    0.000)   13.387
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.524   (   9.415    6.995    0.000)   11.729
   1.803   (   4.414    5.475    0.000)    7.033
   1.970   (  -4.482    0.572    0.000)    4.518
   2.339   (   0.645   -0.465    0.000)    0.795
   2.587   (   4.124   -2.080    0.000)    4.619
   2.704   (  -4.147   -6.129    0.000)    7.400
   3.485   (  11.362   -2.187    0.000)   11.571
   3.872   (   7.773    4.327    0.000)    8.897
   4.213   (  -1.905    1.007    0.000)    2.155
   4.312   (   5.741   16.320    0.000)   17.301
   4.704   (   5.716   12.692    0.000)   13.919
   6.884   (   5.193   -1.374    0.000)    5.372
  11.632   (  -0.265   -0.361    0.000)    0.448
  12.087   (   1.591   16.495    0.000)   16.571
  14.348   (   9.766    7.279    0.000)   12.181
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.701   (   7.266    6.051    0.000)    9.456
   1.862   (  -5.630    1.133    0.000)    5.743
   1.867   (   1.589    2.314    0.000)    2.808
   2.359   (   0.912    6.785    0.000)    6.846
   2.610   (  -1.526  -29.959    0.000)   29.997
   2.670   (   0.848    9.430    0.000)    9.468
   3.642   (   3.746   -2.258    0.000)    4.374
   4.010   (   4.852    3.371    0.000)    5.908
   4.161   (  -2.651    1.633    0.000)    3.114
   4.387   (   1.669   16.001    0.000)   16.088
   4.792   (   2.314   15.151    0.000)   15.327
   6.964   (   2.086   -2.774    0.000)    3.471
  11.628   (  -0.099   -0.376    0.000)    0.388
  12.108   (   0.486   18.042    0.000)   18.048
  14.494   (   3.739    8.026    0.000)    8.854
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.497   (   8.545    8.545    0.000)   12.085
   1.809   (   3.340    3.340    0.000)    4.723
   2.027   (  -1.412   -1.412    0.000)    1.997
   2.281   (  -1.592   -1.592    0.000)    2.251
   2.320   ( -12.091  -12.091    0.000)   17.100
   3.086   (   6.867    6.867    0.000)    9.711
   3.157   (  10.165   10.165    0.000)   14.376
   3.833   (   6.830    6.830    0.000)    9.659
   4.244   (  -0.049   -0.049    0.000)    0.069
   4.451   (  13.232   13.232    0.000)   18.712
   4.753   (   8.190    8.190    0.000)   11.583
   6.806   (   2.592    2.592    0.000)    3.666
  11.631   (  -0.372   -0.372    0.000)    0.525
  12.268   (   8.881    8.881    0.000)   12.560
  14.310   (  11.322   11.322    0.000)   16.011
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.660   (   6.956    5.618    0.000)    8.941
   1.857   (   1.515   -1.816    0.000)    2.365
   1.986   (  -2.515    0.972    0.000)    2.696
   2.097   (  -7.819  -17.173    0.000)   18.870
   2.269   (  -0.708   -8.648    0.000)    8.677
   3.083   (  -0.872   18.480    0.000)   18.500
   3.443   (  10.643   -1.702    0.000)   10.778
   3.975   (   6.433    5.285    0.000)    8.325
   4.233   (  -1.094    0.830    0.000)    1.373
   4.660   (   6.518   16.456    0.000)   17.700
   4.927   (   8.110    8.752    0.000)   11.932
   6.865   (   2.886   -0.327    0.000)    2.904
  11.624   (  -0.286   -0.398    0.000)    0.490
  12.407   (   4.630   13.330    0.000)   14.112
  14.546   (  10.296   11.283    0.000)   15.274
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.781   (   4.211    1.416    0.000)    4.442
   1.876   (   0.397  -10.173    0.000)   10.181
   1.927   (  -2.733    4.920    0.000)    5.628
   1.992   (  -2.383  -17.781    0.000)   17.940
   2.257   (  -0.363   -9.055    0.000)    9.062
   3.071   (  -0.286   18.569    0.000)   18.571
   3.587   (   3.313   -2.935    0.000)    4.426
   4.089   (   4.108    4.017    0.000)    5.746
   4.197   (  -2.090    1.732    0.000)    2.714
   4.740   (   1.680   17.330    0.000)   17.411
   5.055   (   3.430   10.090    0.000)   10.657
   6.912   (   1.304   -1.978    0.000)    2.369
  11.620   (  -0.107   -0.414    0.000)    0.427
  12.469   (   1.422   15.502    0.000)   15.567
  14.703   (   4.054   11.546    0.000)   12.237
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.600   ( -16.175  -16.175    0.000)   22.875
   1.751   (   3.051    3.051    0.000)    4.315
   1.967   (  -0.266   -0.266    0.000)    0.377
   2.009   (   1.098    1.098    0.000)    1.553
   2.144   (  -4.014   -4.014    0.000)    5.677
   3.344   (   4.792    4.792    0.000)    6.777
   3.439   (   3.315    3.315    0.000)    4.688
   4.083   (   4.876    4.876    0.000)    6.895
   4.243   (  -0.016   -0.016    0.000)    0.022
   4.957   (  10.528   10.528    0.000)   14.889
   5.084   (   7.053    7.053    0.000)    9.975
   6.870   (   0.638    0.638    0.000)    0.902
  11.617   (  -0.307   -0.307    0.000)    0.434
  12.626   (   7.625    7.625    0.000)   10.783
  14.784   (  10.428   10.428    0.000)   14.748
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.349   (  -6.601  -25.558    0.000)   26.397
   1.791   (   0.890   -0.012    0.000)    0.890
   1.964   (  -0.194    1.220    0.000)    1.235
   2.026   (   0.494    4.026    0.000)    4.056
   2.094   (  -1.004   -6.243    0.000)    6.323
   3.359   (  -0.324    9.290    0.000)    9.295
   3.526   (   2.526   -2.600    0.000)    3.626
   4.169   (   3.116    3.401    0.000)    4.612
   4.229   (  -1.267    1.184    0.000)    1.734
   5.077   (   2.158   13.925    0.000)   14.091
   5.217   (   4.102    5.588    0.000)    6.932
   6.884   (   0.456   -0.811    0.000)    0.930
  11.612   (  -0.115   -0.319    0.000)    0.340
  12.730   (   2.407    9.273    0.000)    9.580
  14.943   (   4.177   10.539    0.000)   11.337
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.931   ( -12.200  -12.200    0.000)   17.253
   1.789   (  -0.067   -0.067    0.000)    0.095
   1.981   (   0.272    0.272    0.000)    0.385
   2.006   (  -2.013   -2.013    0.000)    2.847
   2.084   (   1.408    1.408    0.000)    1.991
   3.475   (   1.517    1.517    0.000)    2.146
   3.491   (   0.064    0.064    0.000)    0.091
   4.224   (   1.753    1.753    0.000)    2.479
   4.242   (  -0.004   -0.004    0.000)    0.005
   5.280   (   4.312    4.312    0.000)    6.099
   5.297   (   2.757    2.757    0.000)    3.899
   6.875   (  -0.152   -0.152    0.000)    0.215
  11.607   (  -0.120   -0.120    0.000)    0.170
  12.857   (   2.981    2.981    0.000)    4.216
  15.105   (   4.233    4.233    0.000)    5.986
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.967   (  13.630   13.630   13.630)   23.608
   0.967   (  13.630   13.630   13.630)   23.608
   1.370   (  20.659   20.659   20.659)   35.783
   2.248   (   2.829    2.829    2.829)    4.900
   2.248   (   2.829    2.829    2.829)    4.900
   2.578   (  -1.643   -1.643   -1.643)    2.846
   3.086   (  -3.798   -3.798   -3.798)    6.579
   3.421   (   5.570    5.570    5.570)    9.648
   3.421   (   5.570    5.570    5.570)    9.648
   4.303   (   2.362    2.362    2.362)    4.092
   4.303   (   2.362    2.362    2.362)    4.092
   6.489   (   3.635    3.635    3.635)    6.296
  11.726   (   2.227    2.227    2.227)    3.858
  11.726   (   2.227    2.227    2.227)    3.858
  13.696   (   2.566    2.566    2.566)    4.444
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.142   (   9.850   13.389   13.389)   21.343
   1.251   (  10.451   14.013   14.013)   22.405
   1.922   (  25.546    8.157    8.157)   28.030
   2.220   (  -1.999    4.835    4.835)    7.124
   2.255   (  -3.803    5.109    5.109)    8.165
   2.760   (  16.083   -7.504   -7.504)   19.269
   3.035   (  -0.968   -7.588   -7.588)   10.775
   3.526   (   7.492    4.823    4.823)   10.132
   3.603   (   8.359    3.205    3.205)    9.509
   4.348   (   2.946    5.358    5.358)    8.130
   4.361   (   2.270    5.207    5.207)    7.706
   6.596   (   6.437    2.034    2.034)    7.050
  11.728   (   0.076    3.574    3.574)    5.055
  11.790   (   1.673    4.677    4.677)    6.823
  13.823   (   9.212    1.264    1.264)    9.384
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.359   (  10.558   10.304   10.304)   17.995
   1.458   (   9.517   11.887   11.887)   19.318
   2.162   (  -3.625    5.251    5.251)    8.263
   2.167   (  -4.236    5.175    5.175)    8.456
   2.350   (  15.627    1.699    1.699)   15.811
   3.041   (   1.336  -11.588  -11.588)   16.443
   3.097   (  14.116   -8.231   -8.231)   18.296
   3.730   (  11.784    5.250    5.250)   13.928
   3.753   (   5.818    0.613    0.613)    5.882
   4.400   (   1.315    7.702    7.702)   10.972
   4.424   (   4.175    9.130    9.130)   13.570
   6.746   (   7.641    0.007    0.007)    7.641
  11.729   (   0.092    3.914    3.914)    5.536
  11.808   (   0.322    6.630    6.630)    9.382
  14.048   (  11.538    1.663    1.663)   11.775
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.574   (   9.909    8.351    8.351)   15.417
   1.645   (   8.230    8.619    8.619)   14.707
   2.067   (  -5.456    5.782    5.782)    9.830
   2.075   (  -4.600    6.598    6.598)   10.403
   2.601   (   8.808   -1.859   -1.859)    9.192
   3.079   (   1.962  -14.685  -14.685)   20.860
   3.315   (   6.936  -10.464  -10.464)   16.343
   3.844   (   2.967   -2.042   -2.042)    4.140
   3.989   (  12.153    7.298    7.298)   15.944
   4.407   (  -0.875    9.290    9.290)   13.168
   4.511   (   3.780   12.429   12.429)   17.978
   6.897   (   6.284   -1.613   -1.613)    6.685
  11.731   (   0.076    4.234    4.234)    5.988
  11.810   (  -0.001    7.389    7.389)   10.450
  14.272   (   9.188    2.739    2.739)    9.971
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.769   (   8.694    7.149    7.149)   13.334
   1.794   (   5.477    3.238    3.238)    7.139
   1.935   (  -7.187    6.389    6.389)   11.545
   1.978   (  -4.207   11.316   11.316)   16.547
   2.723   (   2.952   -3.760   -3.760)    6.082
   3.111   (   0.889  -16.322  -16.322)   23.100
   3.404   (   1.965  -11.825  -11.825)   16.838
   3.883   (   0.859   -3.414   -3.414)    4.905
   4.211   (   8.699    8.281    8.281)   14.588
   4.353   (  -4.472   10.207   10.207)   15.112
   4.568   (   1.507   14.073   14.073)   19.959
   6.991   (   2.433   -2.412   -2.412)    4.190
  11.732   (   0.030    4.427    4.427)    6.261
  11.810   (  -0.018    7.660    7.660)   10.833
  14.408   (   3.475    3.457    3.457)    5.999
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.325   (   8.293    8.293   17.545)   21.104
   1.551   (  10.089   10.089    2.172)   14.433
   2.156   (   8.069    8.069    5.216)   12.547
   2.211   (   1.893    1.893   11.830)   12.129
   2.488   (   4.946    4.946   -1.009)    7.067
   2.686   (   7.682    7.682   -4.600)   11.798
   2.892   (  -5.084   -5.084   -7.843)   10.640
   3.605   (   5.402    5.402    4.731)    8.986
   3.747   (   7.492    7.492   -7.301)   12.868
   4.441   (   5.236    5.236   12.427)   14.465
   4.510   (   6.897    6.897    5.586)   11.241
   6.658   (   3.849    3.849    0.724)    5.491
  11.749   (   0.182    0.182    8.788)    8.792
  11.938   (   6.765    6.765    0.367)    9.574
  13.915   (   7.855    7.855    1.622)   11.226
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.504   (   8.676    6.555   13.438)   17.287
   1.710   (   5.971    8.291    1.826)   10.379
   2.138   (  -2.676   -2.001    7.113)    7.859
   2.326   (   1.377    2.474    2.492)    3.772
   2.526   (   4.320    5.278    3.005)    7.453
   2.754   (  -3.680   -3.388   -3.870)    6.324
   3.021   (  12.730   -0.872  -11.368)   17.090
   3.755   (   7.326    3.794    1.935)    8.474
   3.903   (   8.682    6.585   -7.750)   13.372
   4.548   (   4.330    8.671   15.536)   18.311
   4.655   (   7.134    9.994    7.533)   14.405
   6.752   (   5.042    0.683   -1.031)    5.191
  11.751   (   0.017    0.049    9.778)    9.778
  12.044   (   3.494   12.552    0.740)   13.050
  14.132   (  12.002    6.809    1.962)   13.938
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.678   (   7.707    4.268    9.882)   13.239
   1.810   (   3.679    5.090    0.738)    6.324
   2.078   (  -3.098   -0.015    9.789)   10.268
   2.321   (  -0.291    3.312   -1.133)    3.512
   2.607   (   0.976  -17.774   -0.266)   17.803
   2.703   (   0.247    7.675    2.728)    8.149
   3.215   (   6.090   -1.459  -19.308)   20.298
   3.853   (   2.427    3.590   -6.900)    8.148
   4.126   (  11.181    4.944    3.890)   12.830
   4.606   (   1.186   11.179   15.417)   19.081
   4.797   (   6.001   12.632    8.898)   16.576
   6.857   (   4.593   -1.952   -2.523)    5.592
  11.751   (  -0.010   -0.103   10.463)   10.464
  12.094   (   1.561   16.077    0.876)   16.176
  14.373   (  10.130    7.171    2.576)   12.676
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.807   (   3.735    0.808    3.533)    5.204
   1.870   (   2.435    2.742    1.130)    3.837
   2.007   (  -3.331    3.528   15.111)   15.871
   2.327   (   0.485    7.323   -2.727)    7.829
   2.584   (  -1.317  -29.048   -2.936)   29.226
   2.737   (   1.307   10.058    7.382)   12.545
   3.295   (   1.810   -2.382  -24.351)   24.534
   3.881   (   0.549    3.551   -9.511)   10.167
   4.310   (   5.343    2.452    9.305)   11.007
   4.605   (  -0.765   13.726   14.624)   20.071
   4.886   (   2.278   14.301    9.016)   17.059
   6.928   (   1.892   -3.323   -3.294)    5.047
  11.751   (  -0.006   -0.197   10.846)   10.848
  12.114   (   0.459   17.613    0.835)   17.638
  14.524   (   3.880    7.757    3.023)    9.185
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.642   (   6.359    6.359    9.866)   13.350
   1.822   (   3.286    3.286    1.248)    4.812
   2.094   (  -1.933   -1.933    4.069)    4.902
   2.310   ( -12.226  -12.226   -1.000)   17.318
   2.359   (   0.049    0.049    8.887)    8.887
   2.914   (   4.669    4.669  -14.280)   15.732
   3.101   (   7.850    7.850   -3.549)   11.656
   3.891   (   7.801    7.801    3.766)   11.657
   4.026   (   6.384    6.384  -10.283)   13.684
   4.739   (   8.435    8.435   16.864)   20.657
   4.839   (   8.324    8.324    7.898)   14.175
   6.778   (   1.840    1.840   -2.487)    3.600
  11.749   (  -0.152   -0.152   10.580)   10.582
  12.271   (   8.783    8.783    0.269)   12.424
  14.329   (  11.453   11.453    2.113)   16.334
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.758   (   4.473    3.214    6.176)    8.276
   1.865   (   1.077   -1.972    0.771)    2.375
   2.050   (  -3.165   -5.380    2.208)    6.621
   2.112   (  -5.167  -13.469    1.765)   14.534
   2.355   (  -0.676   -1.982   10.311)   10.521
   2.974   (   2.256    9.081  -10.231)   13.865
   3.221   (   3.188    1.570  -14.708)   15.131
   4.000   (   2.744    9.940   -2.269)   10.559
   4.204   (   9.230    2.894    2.191)    9.918
   4.861   (   3.405   12.599   14.389)   19.426
   5.016   (   7.954    8.386    8.370)   14.271
   6.825   (   2.364   -0.953   -3.772)    4.552
  11.746   (  -0.144   -0.308   11.005)   11.010
  12.408   (   4.536   13.106    0.104)   13.870
  14.569   (  10.431   11.162    2.297)   15.449
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.818   (   1.238    0.075    2.308)    2.620
   1.875   (   0.006  -12.798   -0.262)   12.801
   1.988   (  -1.286  -13.924    0.299)   13.986
   2.070   (  -0.639    0.417    7.360)    7.400
   2.339   (  -0.522   -1.634   10.279)   10.422
   3.011   (   1.028   10.431   -3.822)   11.157
   3.254   (   0.525   -1.683  -24.390)   24.453
   4.033   (   0.688   10.295   -5.300)   11.600
   4.347   (   3.755    1.356    9.700)   10.489
   4.900   (   0.713   14.214   12.694)   19.071
   5.137   (   3.144    9.630    7.748)   12.754
   6.864   (   1.131   -2.539   -4.453)    5.250
  11.744   (  -0.057   -0.402   11.222)   11.229
  12.469   (   1.378   15.227   -0.025)   15.289
  14.727   (   4.113   11.272    2.457)   12.248
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.591   ( -16.508  -16.508   -0.822)   23.361
   1.807   (   1.419    1.419    3.890)    4.377
   1.976   (  -0.019   -0.019    0.795)    0.795
   2.047   (  -0.619   -0.619    0.841)    1.214
   2.343   (   0.211    0.211   19.360)   19.362
   3.084   (   3.130    3.130  -20.953)   21.415
   3.243   (  -0.084   -0.084  -14.906)   14.907
   4.195   (   6.153    6.153    7.129)   11.249
   4.269   (   4.543    4.543    1.523)    6.602
   5.093   (   8.009    8.009   11.514)   16.152
   5.163   (   6.693    6.693    7.318)   11.964
   6.818   (   0.216    0.216   -4.845)    4.855
  11.740   (  -0.265   -0.265   11.192)   11.198
  12.623   (   7.460    7.460   -0.292)   10.554
  14.804   (  10.367   10.367    2.077)   14.808
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.337   (  -6.713  -25.691   -1.209)   26.581
   1.822   (   0.206    0.019    2.672)    2.680
   1.974   (  -0.088    1.251    0.972)    1.586
   2.046   (   0.177   -2.376    0.500)    2.434
   2.349   (   0.240    1.958   22.395)   22.482
   3.137   (   1.887    2.807  -18.511)   18.818
   3.222   (  -1.355   -1.433  -22.437)   22.523
   4.248   (   0.764    8.938    4.271)    9.935
   4.370   (   3.015    0.807    9.320)    9.829
   5.178   (   1.503   11.693    9.030)   14.850
   5.290   (   3.788    5.181    6.920)    9.438
   6.826   (   0.320   -1.165   -5.397)    5.531
  11.736   (  -0.103   -0.338   11.266)   11.271
  12.724   (   2.341    9.073   -0.523)    9.385
  14.963   (   4.163   10.363    1.986)   11.343
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.913   ( -12.413  -12.413   -1.617)   17.629
   1.821   (  -0.041   -0.041    2.912)    2.913
   1.992   (   0.267    0.267    0.963)    1.035
   2.006   (  -0.990   -0.990    0.637)    1.539
   2.387   (   1.105    1.105   25.479)   25.527
   3.170   (   1.009    1.009  -24.223)   24.265
   3.194   (  -1.201   -1.201  -23.298)   23.359
   4.373   (   2.200    2.200    9.490)    9.986
   4.384   (   1.180    1.180    8.967)    9.121
   5.352   (   3.616    3.616    6.716)    8.441
   5.363   (   2.579    2.579    6.265)    7.249
   6.812   (  -0.261   -0.261   -5.853)    5.865
  11.731   (  -0.130   -0.130   11.254)   11.256
  12.848   (   2.904    2.904   -0.805)    4.184
  15.122   (   4.176    4.176    1.745)    6.158
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.617   (   7.011    7.011    7.011)   12.144
   1.617   (   7.011    7.011    7.011)   12.144
   2.239   (   0.460    0.460    0.460)    0.797
   2.458   (   2.570    2.570    2.570)    4.452
   2.458   (   2.570    2.570    2.570)    4.452
   2.720   (  -7.872   -7.872   -7.872)   13.634
   2.738   (  13.356   13.356   13.356)   23.133
   3.691   (   2.971    2.971    2.971)    5.145
   3.691   (   2.971    2.971    2.971)    5.145
   4.670   (   9.143    9.143    9.143)   15.837
   4.670   (   9.143    9.143    9.143)   15.837
   6.682   (   1.623    1.623    1.623)    2.811
  11.949   (   4.804    4.804    4.804)    8.320
  11.949   (   4.804    4.804    4.804)    8.320
  14.005   (   7.508    7.508    7.508)   13.004
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.715   (   5.172    6.203    6.203)   10.184
   1.772   (   5.265    4.278    4.278)    8.020
   2.191   (  -2.677   -2.073   -2.073)    3.970
   2.413   (  -2.238    3.661    3.661)    5.640
   2.578   (  -5.769  -10.060  -10.060)   15.352
   2.602   (   7.192   -2.863   -2.863)    8.254
   2.987   (   6.764   12.366   12.366)   18.751
   3.732   (   1.072    0.843    0.843)    1.603
   3.890   (  15.067    2.098    2.098)   15.356
   4.822   (   5.453   11.574   11.574)   17.253
   4.866   (   9.146   10.898   10.898)   17.922
   6.727   (   2.700   -1.201   -1.201)    3.190
  11.963   (   0.706    7.054    7.054)   10.001
  12.095   (   4.999    6.213    6.213)   10.110
  14.219   (  12.130    6.780    6.780)   15.463
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.823   (   4.884    3.719    3.719)    7.178
   1.839   (   1.346    2.776    2.776)    4.151
   2.173   (   0.907   -2.664   -2.664)    3.875
   2.370   (  -1.747    3.287    3.287)    4.966
   2.478   (  -3.852  -11.922  -11.922)   17.294
   2.736   (   5.377   -8.160   -8.160)   12.732
   3.042   (  -0.072   12.189   12.189)   17.238
   3.743   (   0.171   -0.576   -0.576)    0.832
   4.220   (  14.491    3.541    3.541)   15.332
   4.897   (   1.867   13.022   13.022)   18.510
   5.032   (   6.431   11.820   11.820)   17.911
   6.789   (   3.000   -3.630   -3.630)    5.946
  11.976   (   0.552    7.205    7.205)   10.204
  12.166   (   2.104    8.550    8.550)   12.273
  14.467   (  10.621    6.732    6.732)   14.263
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.847   (  -0.040    1.349    1.349)    1.908
   1.902   (   2.233    1.732    1.732)    3.315
   2.202   (   1.092   -2.233   -2.233)    3.342
   2.342   (  -0.843    6.491    6.491)    9.218
   2.421   (  -1.409  -15.727  -15.727)   22.285
   2.817   (   2.163  -14.216  -14.216)   20.220
   3.025   (  -0.837   15.421   15.421)   21.825
   3.744   (  -0.023   -1.299   -1.299)    1.837
   4.442   (   5.950    3.777    3.777)    7.996
   4.909   (  -0.232   14.166   14.166)   20.035
   5.125   (   2.288   12.183   12.183)   17.381
   6.838   (   1.369   -4.924   -4.924)    7.097
  11.985   (   0.209    7.294    7.294)   10.318
  12.193   (   0.566    9.541    9.541)   13.505
  14.627   (   4.117    6.978    6.978)   10.693
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.807   (   3.581    3.581    6.034)    7.878
   1.859   (   3.187    3.187    2.206)    5.018
   2.125   (  -3.163   -3.163   -0.582)    4.511
   2.270   ( -12.490  -12.490   -3.452)   17.998
   2.494   (   1.074    1.074   -2.221)    2.691
   2.582   (   0.231    0.231  -15.569)   15.573
   3.215   (   4.990    4.990   18.139)   19.463
   3.890   (   8.356    8.356   -2.837)   12.152
   3.945   (   8.608    8.608    1.255)   12.238
   5.051   (   8.649    8.649   11.454)   16.757
   5.058   (   8.741    8.741   13.293)   18.152
   6.710   (  -0.227   -0.227   -3.744)    3.758
  12.022   (   0.950    0.950   14.163)   14.227
  12.279   (   8.516    8.516    0.462)   12.053
  14.414   (  11.332   11.332    6.313)   17.224
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.864   (   1.289    0.806    3.797)    4.090
   1.888   (  -0.385   -2.658    1.592)    3.122
   2.046   (  -5.283  -12.009   -1.538)   13.210
   2.113   (  -0.734   -7.269   -1.766)    7.516
   2.484   (  -1.304    3.388   -5.069)    6.235
   2.620   (   2.441   -3.207  -21.642)   22.014
   3.237   (  -1.082    2.078   19.804)   19.942
   3.928   (   0.414   17.208   -2.962)   17.466
   4.272   (  14.963    1.669    2.616)   15.282
   5.155   (   2.883   11.467   13.257)   17.764
   5.244   (   7.295    8.024   12.038)   16.202
   6.719   (   0.931   -2.796   -5.989)    6.675
  12.035   (   0.369    0.672   15.137)   15.157
  12.414   (   4.467   12.007    0.600)   12.825
  14.653   (  10.485   10.629    5.961)   16.076
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.846   (  -1.552  -18.363   -2.184)   18.557
   1.876   (  -0.374   -0.325    2.324)    2.377
   2.000   (  -0.106   -6.516    0.970)    6.589
   2.124   (   0.609   -4.070   -0.172)    4.119
   2.461   (  -0.657    4.090   -5.386)    6.795
   2.660   (   1.133   -3.566  -28.045)   28.293
   3.210   (  -0.904    0.328   22.800)   22.821
   3.932   (   0.081   17.815   -3.869)   18.230
   4.493   (   5.633    1.524    4.326)    7.264
   5.186   (   0.454   12.176   13.711)   18.342
   5.349   (   2.616    8.783   11.346)   14.584
   6.739   (   0.661   -4.247   -7.174)    8.363
  12.039   (   0.114    0.418   15.528)   15.534
  12.473   (   1.322   13.842    0.698)   13.922
  14.812   (   4.169   10.426    5.880)   12.675
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.564   ( -17.192  -17.192   -1.772)   24.378
   1.893   (   0.562    0.562    4.287)    4.360
   2.000   (   0.551    0.551    1.453)    1.648
   2.047   (  -0.672   -0.672   -0.538)    1.093
   2.537   (   1.106    1.106  -12.360)   12.459
   2.582   (   0.053    0.053  -25.769)   25.769
   3.206   (  -3.084   -3.084   21.379)   21.819
   4.288   (   8.475    8.475    0.361)   11.990
   4.302   (   7.471    7.471    2.907)   10.959
   5.364   (   6.106    6.106   11.100)   14.063
   5.364   (   5.829    5.829   13.216)   15.576
   6.679   (  -1.014   -1.014   -8.064)    8.191
  12.040   (   0.031    0.031   16.249)   16.249
  12.615   (   7.018    7.018   -0.490)    9.936
  14.878   (  10.032   10.032    5.246)   15.126
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.300   (  -7.004  -25.910   -2.277)   26.936
   1.893   (  -0.227    0.673    3.944)    4.008
   2.002   (   0.058    1.177    1.534)    1.934
   2.056   (   0.542   -2.371    0.419)    2.468
   2.534   (  -0.401    3.036  -17.146)   17.418
   2.609   (   1.183   -1.674  -29.543)   29.614
   3.150   (  -1.692   -4.158   25.415)   25.809
   4.304   (   0.204   15.193    1.157)   15.239
   4.514   (   5.397    0.581    4.281)    6.913
   5.409   (   0.873    8.880   12.022)   14.971
   5.485   (   3.065    4.390   10.859)   12.107
   6.670   (  -0.064   -2.239   -9.100)    9.372
  12.039   (  -0.024   -0.191   16.544)   16.546
  12.710   (   2.182    8.467   -0.831)    8.783
  15.033   (   4.064    9.782    4.882)   11.664
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.865   ( -13.030  -13.030   -2.891)   18.652
   1.901   (   0.171    0.171    4.534)    4.541
   2.018   (   0.200    0.200    1.385)    1.413
   2.026   (  -0.466   -0.466    1.016)    1.210
   2.581   (   0.745    0.745  -27.225)   27.245
   2.588   (   0.129    0.129  -30.193)   30.194
   3.078   (  -2.081   -2.081   31.361)   31.499
   4.520   (   2.728    2.728    4.331)    5.801
   4.521   (   2.661    2.661    3.971)    5.470
   5.540   (   2.361    2.361   10.533)   11.049
   5.542   (   2.234    2.234   10.140)   10.621
   6.642   (  -0.574   -0.574   -9.999)   10.032
  12.036   (  -0.109   -0.109   16.739)   16.740
  12.825   (   2.695    2.695   -1.357)    4.046
  15.184   (   3.983    3.983    4.338)    7.110
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.908   (   2.898    2.898    2.898)    5.019
   1.908   (   2.898    2.898    2.898)    5.019
   2.091   (  -2.298   -2.298   -2.298)    3.981
   2.114   ( -11.255  -11.255  -11.255)   19.494
   2.356   (  -6.475   -6.475   -6.475)   11.216
   2.356   (  -6.475   -6.475   -6.475)   11.216
   3.623   (  12.511   12.511   12.511)   21.670
   3.950   (   5.718    5.718    5.718)    9.904
   3.950   (   5.718    5.718    5.718)    9.904
   5.285   (   9.391    9.391    9.391)   16.266
   5.285   (   9.391    9.391    9.391)   16.266
   6.636   (  -3.031   -3.031   -3.031)    5.250
  12.290   (   5.613    5.613    5.613)    9.721
  12.290   (   5.613    5.613    5.613)    9.721
  14.593   (  10.497   10.497   10.497)   18.182
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.841   (  -8.932   -9.165   -9.165)   15.741
   1.951   (  -0.069    2.545    2.545)    3.599
   1.992   (   1.327   -1.660   -1.660)    2.697
   2.058   (   0.181   -2.549   -2.549)    3.609
   2.231   (  -4.433  -12.118  -12.118)   17.701
   2.287   (  -2.176  -10.048  -10.048)   14.376
   3.716   (   0.832   15.606   15.606)   22.085
   3.933   (  -0.909    8.861    8.861)   12.564
   4.305   (  16.515    0.807    0.807)   16.555
   5.409   (   3.638   10.832   10.832)   15.745
   5.469   (   6.940    8.930    8.930)   14.411
   6.593   (  -1.043   -5.586   -5.586)    7.969
  12.299   (   0.400    7.831    7.831)   11.083
  12.458   (   5.547    5.094    5.094)    9.092
  14.818   (   9.999    9.567    9.567)   16.824
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.701   (  -3.631  -14.607  -14.607)   20.974
   1.930   (  -0.903    2.444    2.444)    3.572
   2.025   (   1.071    0.793    0.793)    1.551
   2.093   (   2.390   -4.469   -4.469)    6.757
   2.157   (  -2.875  -13.709  -13.709)   19.600
   2.255   (  -0.885  -11.019  -11.019)   15.608
   3.722   (   0.024   15.541   15.541)   21.978
   3.918   (  -0.411    8.851    8.851)   12.524
   4.548   (   6.127    1.406    1.406)    6.442
   5.450   (   0.696   11.259   11.259)   15.938
   5.567   (   2.411    8.914    8.914)   12.834
   6.586   (   0.018   -6.965   -6.965)    9.851
  12.305   (   0.151    7.633    7.633)   10.796
  12.531   (   1.618    6.162    6.162)    8.864
  14.972   (   4.048    9.106    9.106)   13.499
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.511   ( -17.084  -17.084   -3.716)   24.444
   1.972   (   0.558    0.558    0.712)    1.063
   2.014   (  -1.960   -1.960   -3.588)    4.534
   2.031   (   2.065    2.065   -1.077)    3.113
   2.055   (  -5.823   -5.823  -21.049)   22.603
   2.154   (  -3.043   -3.043  -15.522)   16.107
   3.785   (  -0.751   -0.751   28.048)   28.068
   4.289   (   9.543    9.543   -0.161)   13.497
   4.329   (   8.099    8.099    0.148)   11.455
   5.592   (   5.704    5.704   10.124)   12.945
   5.611   (   5.574    5.574    9.998)   12.732
   6.502   (  -2.952   -2.952   -8.347)    9.333
  12.364   (   0.778    0.778   13.855)   13.898
  12.604   (   6.435    6.435   -0.532)    9.115
  15.024   (   9.252    9.252    8.432)   15.566
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.241   (  -7.366  -25.089   -3.598)   26.394
   1.938   (  -2.415   -3.092  -22.235)   22.579
   1.978   (  -0.016    1.101    2.305)    2.554
   2.015   (  -0.325   -7.488  -16.090)   17.750
   2.069   (   1.090    1.256   -2.920)    3.361
   2.118   (  -0.831   -1.637   -8.922)    9.109
   3.769   (  -0.468   -1.373   28.946)   28.983
   4.308   (  -0.027   17.982   -0.519)   17.989
   4.564   (   6.282    0.306    1.007)    6.370
   5.647   (   0.878    7.489   10.338)   12.796
   5.706   (   2.482    4.286    9.564)   10.770
   6.467   (  -0.629   -4.074   -9.757)   10.592
  12.371   (   0.098    0.757   14.129)   14.150
  12.693   (   2.089    7.187   -0.597)    7.509
  15.167   (   3.805    8.793    7.782)   12.343
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.798   ( -13.943  -13.943   -3.463)   20.020
   1.896   (  -1.851   -1.851  -27.991)   28.113
   1.899   (  -2.028   -2.028  -26.614)   26.768
   1.999   (   0.483    0.483    4.453)    4.505
   2.090   (   0.717    0.717   -2.182)    2.406
   2.102   (  -0.399   -0.399   -3.298)    3.346
   3.746   (  -0.636   -0.636   29.843)   29.857
   4.565   (   3.354    3.354    0.663)    4.789
   4.568   (   3.008    3.008    0.762)    4.322
   5.753   (   1.974    1.974    9.509)    9.910
   5.757   (   1.627    1.627    9.659)    9.929
   6.416   (  -1.048   -1.048  -11.103)   11.202
  12.377   (   0.026    0.026   14.866)   14.866
  12.794   (   2.421    2.421   -1.480)    3.730
  15.304   (   3.637    3.637    7.007)    8.692
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.357   ( -12.460  -12.460  -12.460)   21.581
   1.672   ( -13.422  -13.422  -13.422)   23.248
   1.672   ( -13.422  -13.422  -13.422)   23.248
   2.045   (   0.771    0.771    0.771)    1.336
   2.068   (   1.301    1.301    1.301)    2.253
   2.068   (   1.301    1.301    1.301)    2.253
   4.248   (   7.411    7.411    7.411)   12.837
   4.322   (   5.366    5.366    5.366)    9.293
   4.322   (   5.366    5.366    5.366)    9.293
   5.770   (   5.846    5.846    5.846)   10.126
   5.770   (   5.846    5.846    5.846)   10.126
   6.355   (  -5.247   -5.247   -5.247)    9.088
  12.594   (   3.795    3.795    3.795)    6.574
  12.594   (   3.795    3.795    3.795)    6.574
  15.207   (   8.308    8.308    8.308)   14.390
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.101   (  -8.134  -13.257  -13.257)   20.436
   1.417   (  -6.696  -17.487  -17.487)   25.621
   1.573   (  -1.569  -19.734  -19.734)   27.952
   2.045   (  -0.082    2.547    2.547)    3.603
   2.090   (   0.557    1.423    1.423)    2.088
   2.100   (   0.910    0.665    0.665)    1.309
   4.265   (  -0.118   10.215   10.215)   14.447
   4.310   (  -0.332    8.469    8.469)   11.981
   4.571   (   6.635   -0.110   -0.110)    6.637
   5.827   (   0.969    6.728    6.728)    9.563
   5.864   (   2.299    5.285    5.285)    7.819
   6.290   (  -1.277   -6.659   -6.659)    9.503
  12.588   (  -0.151    5.574    5.574)    7.884
  12.701   (   2.463    2.269    2.269)    4.045
  15.337   (   3.467    7.716    7.716)   11.449
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.728   ( -14.681  -14.681   -3.131)   20.997
   1.254   (  -4.454   -4.454  -30.210)   30.860
   1.273   (  -5.937   -5.937  -29.618)   30.785
   2.080   (   0.798    0.798    3.133)    3.330
   2.111   (   0.571    0.571    1.456)    1.665
   2.113   (   0.410    0.410    1.248)    1.376
   4.273   (  -0.357   -0.357   19.976)   19.982
   4.565   (   3.523    3.523   -0.389)    4.998
   4.569   (   3.121    3.121   -0.393)    4.432
   5.921   (   1.813    1.813    6.289)    6.791
   5.926   (   1.442    1.442    6.260)    6.584
   6.205   (  -1.658   -1.658   -8.578)    8.893
  12.628   (   0.239    0.239    8.958)    8.964
  12.766   (   2.168    2.168   -1.133)    3.269
  15.457   (   3.243    3.243    7.045)    8.407
======================= 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 1.10e-04 8.55e-05 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.641   (  -8.945   -8.945   -8.945)   15.493
   0.749   ( -13.807  -13.807  -13.807)   23.915
   0.749   ( -13.807  -13.807  -13.807)   23.915
   2.124   (   1.083    1.083    1.083)    1.875
   2.133   (   0.574    0.574    0.574)    0.994
   2.133   (   0.574    0.574    0.574)    0.994
   4.554   (   2.429    2.429    2.429)    4.207
   4.561   (   2.017    2.017    2.017)    3.493
   4.561   (   2.017    2.017    2.017)    3.493
   6.010   (   1.875    1.875    1.875)    3.248
   6.010   (   1.875    1.875    1.875)    3.248
   6.082   (  -2.677   -2.677   -2.677)    4.636
  12.750   (   1.204    1.204    1.204)    2.085
  12.750   (   1.204    1.204    1.204)    2.085
  15.568   (   2.989    2.989    2.989)    5.177
=================== 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     20.642     20.642     20.642     -0.000     -0.000      0.000 3/19965
   20.0      8.446      8.446      8.446     -0.000     -0.000      0.000 3/19965
   30.0      5.426      5.426      5.426     -0.000     -0.000      0.000 3/19965
   40.0      4.073      4.073      4.073     -0.000     -0.000      0.000 3/19965
   50.0      3.321      3.321      3.321      0.000     -0.000      0.000 3/19965
   60.0      2.838      2.838      2.838      0.000     -0.000      0.000 3/19965
   70.0      2.495      2.495      2.495      0.000     -0.000      0.000 3/19965
   80.0      2.235      2.235      2.235      0.000     -0.000      0.000 3/19965
   90.0      2.029      2.029      2.029      0.000     -0.000      0.000 3/19965
  100.0      1.860      1.860      1.860      0.000     -0.000      0.000 3/19965
  110.0      1.718      1.718      1.718      0.000     -0.000      0.000 3/19965
  120.0      1.596      1.596      1.596      0.000     -0.000      0.000 3/19965
  130.0      1.491      1.491      1.491      0.000     -0.000      0.000 3/19965
  140.0      1.399      1.399      1.399      0.000     -0.000      0.000 3/19965
  150.0      1.318      1.318      1.318      0.000     -0.000      0.000 3/19965
  160.0      1.246      1.246      1.246      0.000     -0.000      0.000 3/19965
  170.0      1.181      1.181      1.181      0.000     -0.000      0.000 3/19965
  180.0      1.122      1.122      1.122      0.000     -0.000      0.000 3/19965
  190.0      1.069      1.069      1.069      0.000     -0.000      0.000 3/19965
  200.0      1.021      1.021      1.021      0.000     -0.000      0.000 3/19965
  210.0      0.976      0.976      0.976      0.000     -0.000      0.000 3/19965
  220.0      0.936      0.936      0.936      0.000     -0.000      0.000 3/19965
  230.0      0.898      0.898      0.898      0.000     -0.000      0.000 3/19965
  240.0      0.863      0.863      0.863      0.000     -0.000      0.000 3/19965
  250.0      0.831      0.831      0.831      0.000     -0.000      0.000 3/19965
  260.0      0.801      0.801      0.801      0.000     -0.000      0.000 3/19965
  270.0      0.774      0.774      0.774      0.000     -0.000      0.000 3/19965
  280.0      0.748      0.748      0.748      0.000     -0.000      0.000 3/19965
  290.0      0.723      0.723      0.723      0.000     -0.000      0.000 3/19965
  300.0      0.700      0.700      0.700      0.000     -0.000      0.000 3/19965
  310.0      0.679      0.679      0.679      0.000     -0.000      0.000 3/19965
  320.0      0.659      0.659      0.659      0.000     -0.000      0.000 3/19965
  330.0      0.640      0.640      0.640      0.000     -0.000      0.000 3/19965
  340.0      0.622      0.622      0.622      0.000     -0.000      0.000 3/19965
  350.0      0.605      0.605      0.605      0.000     -0.000      0.000 3/19965
  360.0      0.588      0.588      0.588      0.000     -0.000      0.000 3/19965
  370.0      0.573      0.573      0.573      0.000     -0.000      0.000 3/19965
  380.0      0.559      0.559      0.559      0.000     -0.000      0.000 3/19965
  390.0      0.545      0.545      0.545      0.000     -0.000      0.000 3/19965
  400.0      0.532      0.532      0.532      0.000     -0.000      0.000 3/19965
  410.0      0.519      0.519      0.519      0.000     -0.000      0.000 3/19965
  420.0      0.507      0.507      0.507      0.000     -0.000      0.000 3/19965
  430.0      0.496      0.496      0.496      0.000     -0.000      0.000 3/19965
  440.0      0.485      0.485      0.485      0.000     -0.000      0.000 3/19965
  450.0      0.474      0.474      0.474      0.000     -0.000      0.000 3/19965
  460.0      0.464      0.464      0.464      0.000     -0.000      0.000 3/19965
  470.0      0.454      0.454      0.454      0.000     -0.000      0.000 3/19965
  480.0      0.445      0.445      0.445      0.000     -0.000      0.000 3/19965
  490.0      0.436      0.436      0.436      0.000     -0.000      0.000 3/19965
  500.0      0.428      0.428      0.428      0.000     -0.000      0.000 3/19965
  510.0      0.420      0.420      0.420      0.000     -0.000      0.000 3/19965
  520.0      0.412      0.412      0.412      0.000     -0.000      0.000 3/19965
  530.0      0.404      0.404      0.404      0.000     -0.000      0.000 3/19965
  540.0      0.397      0.397      0.397      0.000     -0.000      0.000 3/19965
  550.0      0.390      0.390      0.390      0.000     -0.000      0.000 3/19965
  560.0      0.383      0.383      0.383      0.000     -0.000      0.000 3/19965
  570.0      0.376      0.376      0.376      0.000     -0.000      0.000 3/19965
  580.0      0.370      0.370      0.370      0.000     -0.000      0.000 3/19965
  590.0      0.364      0.364      0.364      0.000     -0.000      0.000 3/19965
  600.0      0.358      0.358      0.358      0.000     -0.000      0.000 3/19965
  610.0      0.352      0.352      0.352      0.000     -0.000      0.000 3/19965
  620.0      0.346      0.346      0.346      0.000     -0.000      0.000 3/19965
  630.0      0.341      0.341      0.341      0.000     -0.000      0.000 3/19965
  640.0      0.336      0.336      0.336      0.000     -0.000      0.000 3/19965
  650.0      0.331      0.331      0.331      0.000     -0.000      0.000 3/19965
  660.0      0.326      0.326      0.326      0.000     -0.000      0.000 3/19965
  670.0      0.321      0.321      0.321      0.000     -0.000      0.000 3/19965
  680.0      0.316      0.316      0.316      0.000     -0.000      0.000 3/19965
  690.0      0.312      0.312      0.312      0.000     -0.000      0.000 3/19965
  700.0      0.307      0.307      0.307      0.000     -0.000      0.000 3/19965
  710.0      0.303      0.303      0.303      0.000     -0.000      0.000 3/19965
  720.0      0.299      0.299      0.299      0.000     -0.000      0.000 3/19965
  730.0      0.295      0.295      0.295      0.000     -0.000      0.000 3/19965
  740.0      0.291      0.291      0.291      0.000     -0.000      0.000 3/19965
  750.0      0.287      0.287      0.287      0.000     -0.000      0.000 3/19965
  760.0      0.283      0.283      0.283      0.000     -0.000      0.000 3/19965
  770.0      0.280      0.280      0.280      0.000     -0.000      0.000 3/19965
  780.0      0.276      0.276      0.276      0.000     -0.000      0.000 3/19965
  790.0      0.273      0.273      0.273      0.000     -0.000      0.000 3/19965
  800.0      0.269      0.269      0.269      0.000     -0.000      0.000 3/19965
  810.0      0.266      0.266      0.266      0.000     -0.000      0.000 3/19965
  820.0      0.263      0.263      0.263      0.000     -0.000      0.000 3/19965
  830.0      0.260      0.260      0.260      0.000     -0.000      0.000 3/19965
  840.0      0.256      0.256      0.256      0.000     -0.000      0.000 3/19965
  850.0      0.253      0.253      0.253      0.000     -0.000      0.000 3/19965
  860.0      0.251      0.251      0.251      0.000     -0.000      0.000 3/19965
  870.0      0.248      0.248      0.248      0.000     -0.000      0.000 3/19965
  880.0      0.245      0.245      0.245      0.000     -0.000      0.000 3/19965
  890.0      0.242      0.242      0.242      0.000     -0.000      0.000 3/19965
  900.0      0.239      0.239      0.239      0.000     -0.000      0.000 3/19965
  910.0      0.237      0.237      0.237      0.000     -0.000      0.000 3/19965
  920.0      0.234      0.234      0.234      0.000     -0.000      0.000 3/19965
  930.0      0.232      0.232      0.232      0.000     -0.000      0.000 3/19965
  940.0      0.229      0.229      0.229      0.000     -0.000      0.000 3/19965
  950.0      0.227      0.227      0.227      0.000     -0.000      0.000 3/19965
  960.0      0.225      0.225      0.225      0.000     -0.000      0.000 3/19965
  970.0      0.222      0.222      0.222      0.000     -0.000      0.000 3/19965
  980.0      0.220      0.220      0.220      0.000     -0.000      0.000 3/19965
  990.0      0.218      0.218      0.218      0.000     -0.000      0.000 3/19965
 1000.0      0.216      0.216      0.216      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:08:50]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

