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

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

-------------------------[time 2026-01-07 22:15: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.487526879999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.487526879999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.487526879999999
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 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078
   *5 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.487526879999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.487526879999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.487526879999999
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 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 4
   *5 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   13.462580639999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.462580639999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.462580639999997
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 Ca  0.16666666666667  0.16666666666667  0.16666666666667  40.078 > 4
   83 Ca  0.50000000000000  0.16666666666667  0.16666666666667  40.078 > 4
   84 Ca  0.83333333333333  0.16666666666667  0.16666666666667  40.078 > 4
   85 Ca  0.16666666666667  0.50000000000000  0.16666666666667  40.078 > 4
   86 Ca  0.50000000000000  0.50000000000000  0.16666666666667  40.078 > 4
   87 Ca  0.83333333333333  0.50000000000000  0.16666666666667  40.078 > 4
   88 Ca  0.16666666666667  0.83333333333333  0.16666666666667  40.078 > 4
   89 Ca  0.50000000000000  0.83333333333333  0.16666666666667  40.078 > 4
   90 Ca  0.83333333333333  0.83333333333333  0.16666666666667  40.078 > 4
   91 Ca  0.16666666666667  0.16666666666667  0.50000000000000  40.078 > 4
   92 Ca  0.50000000000000  0.16666666666667  0.50000000000000  40.078 > 4
   93 Ca  0.83333333333333  0.16666666666667  0.50000000000000  40.078 > 4
   94 Ca  0.16666666666667  0.50000000000000  0.50000000000000  40.078 > 4
   95 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 4
   96 Ca  0.83333333333333  0.50000000000000  0.50000000000000  40.078 > 4
   97 Ca  0.16666666666667  0.83333333333333  0.50000000000000  40.078 > 4
   98 Ca  0.50000000000000  0.83333333333333  0.50000000000000  40.078 > 4
   99 Ca  0.83333333333333  0.83333333333333  0.50000000000000  40.078 > 4
  100 Ca  0.16666666666667  0.16666666666667  0.83333333333333  40.078 > 4
  101 Ca  0.50000000000000  0.16666666666667  0.83333333333333  40.078 > 4
  102 Ca  0.83333333333333  0.16666666666667  0.83333333333333  40.078 > 4
  103 Ca  0.16666666666667  0.50000000000000  0.83333333333333  40.078 > 4
  104 Ca  0.50000000000000  0.50000000000000  0.83333333333333  40.078 > 4
  105 Ca  0.83333333333333  0.50000000000000  0.83333333333333  40.078 > 4
  106 Ca  0.16666666666667  0.83333333333333  0.83333333333333  40.078 > 4
  107 Ca  0.50000000000000  0.83333333333333  0.83333333333333  40.078 > 4
  108 Ca  0.83333333333333  0.83333333333333  0.83333333333333  40.078 > 4
 *109 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 5
  110 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 5
  111 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 5
  112 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 5
  113 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 5
  114 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 5
  115 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 5
  116 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 5
  117 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 5
  118 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 5
  119 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 5
  120 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 5
  121 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 5
  122 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 5
  123 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 5
  124 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 5
  125 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 5
  126 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 5
  127 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 5
  128 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 5
  129 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 5
  130 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 5
  131 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 5
  132 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 5
  133 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 5
  134 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 5
  135 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.4101213    0.0000000    0.0000000
            0.0000000    2.4101213    0.0000000
            0.0000000    0.0000000    2.4101213
-------------------------- Born effective charges --------------------------
    1 F    -1.6095942    0.0000000    0.0000000
            0.0000000   -1.0844960    0.0000000
            0.0000000    0.0000000   -1.0844960
    2 F    -1.0844960    0.0000000    0.0000000
            0.0000000   -1.6095942    0.0000000
            0.0000000    0.0000000   -1.0844960
    3 F    -1.0844960    0.0000000    0.0000000
            0.0000000   -1.0844960    0.0000000
            0.0000000    0.0000000   -1.6095942
    4 Ca    2.4388573    0.0000000    0.0000000
            0.0000000    2.4388573    0.0000000
            0.0000000    0.0000000    2.4388573
    5 Cs    1.3397288    0.0000000    0.0000000
            0.0000000    1.3397288    0.0000000
            0.0000000    0.0000000    1.3397288
----------------------------------------------------------------------------
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.139
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.149
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) 
Force constants are written into "force_constants.hdf5".

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

Summary of calculation was written in "phonopy.yaml".
-------------------------[time 2026-01-07 22:15:39]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 22:15:40]-------------------------
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.487526879999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.487526879999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.487526879999999
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 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078
    5 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.462580639999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.462580639999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.462580639999997
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 Ca  0.16666666666667  0.16666666666667  0.16666666666667  40.078 > 82
   83 Ca  0.50000000000000  0.16666666666667  0.16666666666667  40.078 > 82
   84 Ca  0.83333333333333  0.16666666666667  0.16666666666667  40.078 > 82
   85 Ca  0.16666666666667  0.50000000000000  0.16666666666667  40.078 > 82
   86 Ca  0.50000000000000  0.50000000000000  0.16666666666667  40.078 > 82
   87 Ca  0.83333333333333  0.50000000000000  0.16666666666667  40.078 > 82
   88 Ca  0.16666666666667  0.83333333333333  0.16666666666667  40.078 > 82
   89 Ca  0.50000000000000  0.83333333333333  0.16666666666667  40.078 > 82
   90 Ca  0.83333333333333  0.83333333333333  0.16666666666667  40.078 > 82
   91 Ca  0.16666666666667  0.16666666666667  0.50000000000000  40.078 > 82
   92 Ca  0.50000000000000  0.16666666666667  0.50000000000000  40.078 > 82
   93 Ca  0.83333333333333  0.16666666666667  0.50000000000000  40.078 > 82
   94 Ca  0.16666666666667  0.50000000000000  0.50000000000000  40.078 > 82
   95 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 82
   96 Ca  0.83333333333333  0.50000000000000  0.50000000000000  40.078 > 82
   97 Ca  0.16666666666667  0.83333333333333  0.50000000000000  40.078 > 82
   98 Ca  0.50000000000000  0.83333333333333  0.50000000000000  40.078 > 82
   99 Ca  0.83333333333333  0.83333333333333  0.50000000000000  40.078 > 82
  100 Ca  0.16666666666667  0.16666666666667  0.83333333333333  40.078 > 82
  101 Ca  0.50000000000000  0.16666666666667  0.83333333333333  40.078 > 82
  102 Ca  0.83333333333333  0.16666666666667  0.83333333333333  40.078 > 82
  103 Ca  0.16666666666667  0.50000000000000  0.83333333333333  40.078 > 82
  104 Ca  0.50000000000000  0.50000000000000  0.83333333333333  40.078 > 82
  105 Ca  0.83333333333333  0.50000000000000  0.83333333333333  40.078 > 82
  106 Ca  0.16666666666667  0.83333333333333  0.83333333333333  40.078 > 82
  107 Ca  0.50000000000000  0.83333333333333  0.83333333333333  40.078 > 82
  108 Ca  0.83333333333333  0.83333333333333  0.83333333333333  40.078 > 82
  109 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 109
  110 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 109
  111 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 109
  112 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 109
  113 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 109
  114 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 109
  115 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 109
  116 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 109
  117 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 109
  118 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 109
  119 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 109
  120 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 109
  121 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 109
  122 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 109
  123 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 109
  124 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 109
  125 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 109
  126 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 109
  127 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 109
  128 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 109
  129 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 109
  130 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 109
  131 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 109
  132 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 109
  133 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 109
  134 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 109
  135 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.4101213    0.0000000    0.0000000
            0.0000000    2.4101213    0.0000000
            0.0000000    0.0000000    2.4101213
-------------------------- Born effective charges --------------------------
    1 F    -1.6095942    0.0000000    0.0000000
            0.0000000   -1.0844960    0.0000000
            0.0000000    0.0000000   -1.0844960
    2 F    -1.0844960    0.0000000    0.0000000
            0.0000000   -1.6095942    0.0000000
            0.0000000    0.0000000   -1.0844960
    3 F    -1.0844960    0.0000000    0.0000000
            0.0000000   -1.0844960    0.0000000
            0.0000000    0.0000000   -1.6095942
    4 Ca    2.4388573    0.0000000    0.0000000
            0.0000000    2.4388573    0.0000000
            0.0000000    0.0000000    2.4388573
    5 Cs    1.3397288    0.0000000    0.0000000
            0.0000000    1.3397288    0.0000000
            0.0000000    0.0000000    1.3397288
----------------------------------------------------------------------------
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:15:44]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 22:15:44]-------------------------
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.487526879999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.487526879999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.487526879999999
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 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078
    5 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.462580639999997    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.462580639999997    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.462580639999997
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 Ca  0.16666666666667  0.16666666666667  0.16666666666667  40.078 > 82
   83 Ca  0.50000000000000  0.16666666666667  0.16666666666667  40.078 > 82
   84 Ca  0.83333333333333  0.16666666666667  0.16666666666667  40.078 > 82
   85 Ca  0.16666666666667  0.50000000000000  0.16666666666667  40.078 > 82
   86 Ca  0.50000000000000  0.50000000000000  0.16666666666667  40.078 > 82
   87 Ca  0.83333333333333  0.50000000000000  0.16666666666667  40.078 > 82
   88 Ca  0.16666666666667  0.83333333333333  0.16666666666667  40.078 > 82
   89 Ca  0.50000000000000  0.83333333333333  0.16666666666667  40.078 > 82
   90 Ca  0.83333333333333  0.83333333333333  0.16666666666667  40.078 > 82
   91 Ca  0.16666666666667  0.16666666666667  0.50000000000000  40.078 > 82
   92 Ca  0.50000000000000  0.16666666666667  0.50000000000000  40.078 > 82
   93 Ca  0.83333333333333  0.16666666666667  0.50000000000000  40.078 > 82
   94 Ca  0.16666666666667  0.50000000000000  0.50000000000000  40.078 > 82
   95 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 82
   96 Ca  0.83333333333333  0.50000000000000  0.50000000000000  40.078 > 82
   97 Ca  0.16666666666667  0.83333333333333  0.50000000000000  40.078 > 82
   98 Ca  0.50000000000000  0.83333333333333  0.50000000000000  40.078 > 82
   99 Ca  0.83333333333333  0.83333333333333  0.50000000000000  40.078 > 82
  100 Ca  0.16666666666667  0.16666666666667  0.83333333333333  40.078 > 82
  101 Ca  0.50000000000000  0.16666666666667  0.83333333333333  40.078 > 82
  102 Ca  0.83333333333333  0.16666666666667  0.83333333333333  40.078 > 82
  103 Ca  0.16666666666667  0.50000000000000  0.83333333333333  40.078 > 82
  104 Ca  0.50000000000000  0.50000000000000  0.83333333333333  40.078 > 82
  105 Ca  0.83333333333333  0.50000000000000  0.83333333333333  40.078 > 82
  106 Ca  0.16666666666667  0.83333333333333  0.83333333333333  40.078 > 82
  107 Ca  0.50000000000000  0.83333333333333  0.83333333333333  40.078 > 82
  108 Ca  0.83333333333333  0.83333333333333  0.83333333333333  40.078 > 82
  109 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 109
  110 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 109
  111 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 109
  112 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 109
  113 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 109
  114 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 109
  115 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 109
  116 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 109
  117 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 109
  118 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 109
  119 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 109
  120 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 109
  121 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 109
  122 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 109
  123 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 109
  124 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 109
  125 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 109
  126 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 109
  127 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 109
  128 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 109
  129 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 109
  130 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 109
  131 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 109
  132 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 109
  133 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 109
  134 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 109
  135 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.4101213    0.0000000    0.0000000
            0.0000000    2.4101213    0.0000000
            0.0000000    0.0000000    2.4101213
-------------------------- Born effective charges --------------------------
    1 F    -1.6095942    0.0000000    0.0000000
            0.0000000   -1.0844960    0.0000000
            0.0000000    0.0000000   -1.0844960
    2 F    -1.0844960    0.0000000    0.0000000
            0.0000000   -1.6095942    0.0000000
            0.0000000    0.0000000   -1.0844960
    3 F    -1.0844960    0.0000000    0.0000000
            0.0000000   -1.0844960    0.0000000
            0.0000000    0.0000000   -1.6095942
    4 Ca    2.4388573    0.0000000    0.0000000
            0.0000000    2.4388573    0.0000000
            0.0000000    0.0000000    2.4388573
    5 Cs    1.3397288    0.0000000    0.0000000
            0.0000000    1.3397288    0.0000000
            0.0000000    0.0000000    1.3397288
----------------------------------------------------------------------------
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.93, 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 2.98e-04 0.00e+00 
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.622   (   0.000    0.000    0.000)    0.000
   2.622   (   0.000    0.000    0.000)    0.000
   2.622   (   0.000    0.000    0.000)    0.000
   3.660   (   0.000    0.000    0.000)    0.000
   3.660   (   0.000    0.000    0.000)    0.000
   3.660   (   0.000    0.000    0.000)    0.000
   5.365   (   0.000    0.000    0.000)    0.000
   5.365   (   0.000    0.000    0.000)    0.000
   5.365   (   0.000    0.000    0.000)    0.000
  10.845   (   0.000    0.000    0.000)    0.000
  10.845   (   0.000    0.000    0.000)    0.000
  10.845   (   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 2.98e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.474   (  22.806    0.000    0.000)   22.806
   0.474   (  22.806    0.000    0.000)   22.806
   0.959   (  41.587    0.000    0.000)   41.587
   2.630   (   0.740    0.000    0.000)    0.740
   2.630   (   0.740    0.000    0.000)    0.740
   3.477   (  21.397    0.000    0.000)   21.397
   3.691   (   3.065    0.000    0.000)    3.065
   3.691   (   3.065    0.000    0.000)    3.065
   3.693   (   3.231    0.000    0.000)    3.231
   5.349   (  -1.649    0.000    0.000)    1.649
   5.349   (  -1.649    0.000    0.000)    1.649
   7.059   (  12.716    0.000    0.000)   12.716
  10.844   (  -0.092    0.000    0.000)    0.092
  10.844   (  -0.092    0.000    0.000)    0.092
  13.371   (  -5.799    0.000    0.000)    5.799
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.912   (  20.147    0.000    0.000)   20.147
   0.912   (  20.147    0.000    0.000)   20.147
   1.642   (  26.407    0.000    0.000)   26.407
   2.650   (   1.098    0.000    0.000)    1.098
   2.650   (   1.098    0.000    0.000)    1.098
   3.780   (   5.567    0.000    0.000)    5.567
   3.780   (   5.567    0.000    0.000)    5.567
   3.783   (   5.295    0.000    0.000)    5.295
   4.029   (  29.724    0.000    0.000)   29.724
   5.300   (  -3.133    0.000    0.000)    3.133
   5.300   (  -3.133    0.000    0.000)    3.133
   7.471   (  28.480    0.000    0.000)   28.480
  10.842   (  -0.154    0.000    0.000)    0.154
  10.842   (  -0.154    0.000    0.000)    0.154
  13.189   ( -12.374    0.000    0.000)   12.374
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.280   (  15.897    0.000    0.000)   15.897
   1.280   (  15.897    0.000    0.000)   15.897
   2.070   (  16.735    0.000    0.000)   16.735
   2.671   (   0.953    0.000    0.000)    0.953
   2.671   (   0.953    0.000    0.000)    0.953
   3.896   (   5.577    0.000    0.000)    5.577
   3.909   (   6.865    0.000    0.000)    6.865
   3.909   (   6.865    0.000    0.000)    6.865
   4.563   (  21.438    0.000    0.000)   21.438
   5.225   (  -4.140    0.000    0.000)    4.140
   5.225   (  -4.140    0.000    0.000)    4.140
   8.206   (  42.966    0.000    0.000)   42.966
  10.838   (  -0.165    0.000    0.000)    0.165
  10.838   (  -0.165    0.000    0.000)    0.165
  12.864   ( -19.727    0.000    0.000)   19.727
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.547   (  10.278    0.000    0.000)   10.278
   1.547   (  10.278    0.000    0.000)   10.278
   2.339   (   9.994    0.000    0.000)    9.994
   2.686   (   0.510    0.000    0.000)    0.510
   2.686   (   0.510    0.000    0.000)    0.510
   3.997   (   4.138    0.000    0.000)    4.138
   4.045   (   6.146    0.000    0.000)    6.146
   4.045   (   6.146    0.000    0.000)    6.146
   4.884   (  10.695    0.000    0.000)   10.695
   5.139   (  -4.011    0.000    0.000)    4.011
   5.139   (  -4.011    0.000    0.000)    4.011
   9.132   (  45.792    0.000    0.000)   45.792
  10.835   (  -0.124    0.000    0.000)    0.124
  10.835   (  -0.124    0.000    0.000)    0.124
  12.398   ( -25.593    0.000    0.000)   25.593
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.689   (   3.580    0.000    0.000)    3.580
   1.689   (   3.580    0.000    0.000)    3.580
   2.475   (   3.400    0.000    0.000)    3.400
   2.692   (   0.128    0.000    0.000)    0.128
   2.692   (   0.128    0.000    0.000)    0.128
   4.055   (   1.518    0.000    0.000)    1.518
   4.138   (   2.620    0.000    0.000)    2.620
   4.138   (   2.620    0.000    0.000)    2.620
   5.019   (   3.100    0.000    0.000)    3.100
   5.077   (  -1.812    0.000    0.000)    1.812
   5.077   (  -1.812    0.000    0.000)    1.812
   9.913   (  25.689    0.000    0.000)   25.689
  10.834   (  -0.046    0.000    0.000)    0.046
  10.834   (  -0.046    0.000    0.000)    0.046
  11.911   ( -17.873    0.000    0.000)   17.873
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.665   (  15.669   15.669    0.000)   22.160
   0.826   (  17.358   17.358    0.000)   24.548
   1.225   (  25.348   25.348    0.000)   35.847
   2.641   (   1.072    1.072    0.000)    1.515
   2.791   (   7.680    7.680    0.000)   10.862
   3.435   (   9.138    9.138    0.000)   12.923
   3.633   (  -1.377   -1.377    0.000)    1.947
   3.721   (   2.951    2.951    0.000)    4.174
   3.815   (   7.634    7.634    0.000)   10.796
   5.331   (  -1.706   -1.706    0.000)    2.412
   5.536   (   8.575    8.575    0.000)   12.126
   7.108   (   8.265    8.265    0.000)   11.689
  10.797   (  -2.262   -2.262    0.000)    3.199
  10.843   (  -0.098   -0.098    0.000)    0.139
  13.347   (  -3.857   -3.857    0.000)    5.455
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.011   (  17.223    8.951    0.000)   19.410
   1.125   (  13.868   14.490    0.000)   20.057
   1.730   (  21.261    8.290    0.000)   22.820
   2.671   (   1.791    2.104    0.000)    2.764
   2.855   (   1.358   17.009    0.000)   17.063
   3.615   (   0.712   -6.365    0.000)    6.405
   3.807   (   5.372    2.637    0.000)    5.985
   3.851   (  19.825   -4.918    0.000)   20.426
   4.036   (  17.831   10.515    0.000)   20.700
   5.281   (  -3.242   -1.872    0.000)    3.743
   5.617   (   0.560   23.998    0.000)   24.004
   7.455   (  26.254   -0.334    0.000)   26.256
  10.771   (  -0.599   -5.380    0.000)    5.413
  10.841   (  -0.163   -0.113    0.000)    0.199
  13.199   ( -10.822   -0.117    0.000)   10.823
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.332   (  14.039    4.817    0.000)   14.842
   1.391   (  11.993    8.268    0.000)   14.566
   2.073   (  13.130    0.299    0.000)   13.134
   2.710   (   1.956    3.826    0.000)    4.297
   2.893   (   2.472   18.838    0.000)   19.000
   3.637   (   1.408  -12.371    0.000)   12.451
   3.932   (   6.667    2.225    0.000)    7.028
   4.099   (   8.757    4.398    0.000)    9.800
   4.502   (  22.063    3.218    0.000)   22.297
   5.203   (  -4.295   -2.138    0.000)    4.797
   5.586   (  -3.366   28.662    0.000)   28.859
   8.157   (  41.710   -4.628    0.000)   41.965
  10.764   (  -0.292   -6.404    0.000)    6.410
  10.837   (  -0.175   -0.134    0.000)    0.221
  12.908   ( -17.854    3.471    0.000)   18.188
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.567   (   8.953    1.834    0.000)    9.139
   1.597   (   8.004    3.899    0.000)    8.903
   2.278   (   7.369   -4.573    0.000)    8.673
   2.746   (   1.517    5.859    0.000)    6.052
   2.947   (   2.525   21.256    0.000)   21.406
   3.666   (   1.297  -17.982    0.000)   18.029
   4.064   (   6.042    1.904    0.000)    6.334
   4.247   (   5.708    5.115    0.000)    7.665
   4.871   (  14.809    7.030    0.000)   16.393
   5.114   (  -4.187   -2.479    0.000)    4.866
   5.487   (  -6.244   26.955    0.000)   27.669
   9.057   (  44.372   -7.246    0.000)   44.960
  10.754   (  -0.874   -7.209    0.000)    7.261
  10.834   (  -0.132   -0.153    0.000)    0.202
  12.488   ( -22.679    7.928    0.000)   24.025
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.690   (   3.054    0.069    0.000)    3.055
   1.707   (   2.779    1.643    0.000)    3.229
   2.377   (   2.428   -6.718    0.000)    7.144
   2.768   (   0.574    7.310    0.000)    7.333
   2.985   (   1.044   22.934    0.000)   22.958
   3.685   (   0.500  -21.585    0.000)   21.591
   4.157   (   2.610    1.800    0.000)    3.170
   4.325   (   1.986    4.277    0.000)    4.715
   5.048   (  -1.908   -2.783    0.000)    3.374
   5.126   (  10.593   13.855    0.000)   17.441
   5.338   (  -8.257   22.058    0.000)   23.553
   9.807   (  24.403  -10.151    0.000)   26.430
  10.721   (  -2.126   -9.316    0.000)    9.556
  10.832   (  -0.049   -0.164    0.000)    0.171
  12.075   ( -14.195   13.965    0.000)   19.913
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.230   (  11.641   11.641    0.000)   16.463
   1.368   (   9.743    9.743    0.000)   13.779
   1.923   (   9.293    9.293    0.000)   13.142
   2.733   (   3.904    3.904    0.000)    5.521
   3.105   (   4.810    4.810    0.000)    6.803
   3.577   (   0.042    0.042    0.000)    0.060
   3.872   (   9.305    9.305    0.000)   13.159
   3.884   (   4.804    4.804    0.000)    6.793
   4.348   (  19.739   19.739    0.000)   27.915
   5.225   (  -3.568   -3.568    0.000)    5.046
   6.066   (  17.356   17.356    0.000)   24.545
   7.548   (  12.567   12.567    0.000)   17.773
  10.677   (  -3.513   -3.513    0.000)    4.969
  10.837   (  -0.189   -0.189    0.000)    0.268
  13.145   (  -5.647   -5.647    0.000)    7.987
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.451   (   9.688    6.322    0.000)   11.569
   1.543   (   7.736    6.479    0.000)   10.091
   2.080   (   6.170    0.319    0.000)    6.179
   2.823   (   4.805    7.189    0.000)    8.646
   3.150   (   1.676   -0.845    0.000)    1.877
   3.549   (  -2.125    7.635    0.000)    7.926
   3.996   (   6.045    3.983    0.000)    7.239
   4.095   (   9.775   -1.261    0.000)    9.856
   4.793   (  22.335   21.832    0.000)   31.233
   5.139   (  -4.767   -4.105    0.000)    6.291
   6.247   (   0.679   35.095    0.000)   35.102
   8.043   (  35.274   -5.428    0.000)   35.689
  10.626   (  -1.812   -6.370    0.000)    6.623
  10.833   (  -0.204   -0.224    0.000)    0.303
  12.963   ( -12.283    1.070    0.000)   12.329
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.611   (   5.931    2.212    0.000)    6.330
   1.677   (   5.264    3.842    0.000)    6.517
   2.179   (   3.675   -4.998    0.000)    6.203
   2.916   (   4.054   10.613    0.000)   11.361
   3.194   (   2.502   -9.773    0.000)   10.088
   3.503   (  -2.277   12.228    0.000)   12.438
   4.118   (   5.685    3.260    0.000)    6.554
   4.245   (   5.161   -3.300    0.000)    6.126
   5.040   (  -4.744   -4.824    0.000)    6.766
   5.209   (  18.537   23.405    0.000)   29.857
   6.153   (  -8.492   38.528    0.000)   39.453
   8.849   (  40.557  -12.796    0.000)   42.528
  10.588   (  -2.380   -8.291    0.000)    8.625
  10.830   (  -0.155   -0.256    0.000)    0.299
  12.666   ( -16.089    8.291    0.000)   18.100
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.691   (   1.948   -0.013    0.000)    1.948
   1.750   (   1.829    2.469    0.000)    3.072
   2.230   (   1.324   -7.502    0.000)    7.618
   2.976   (   1.589   12.763    0.000)   12.861
   3.238   (   1.414  -18.989    0.000)   19.042
   3.465   (  -1.236   18.553    0.000)   18.594
   4.207   (   2.576    3.000    0.000)    3.954
   4.312   (   1.579   -4.270    0.000)    4.553
   4.964   (  -2.229   -5.511    0.000)    5.944
   5.533   (  12.228   23.655    0.000)   26.628
   5.953   (  -9.436   38.888    0.000)   40.016
   9.524   (  21.242  -17.173    0.000)   27.316
  10.522   (  -3.435   -9.446    0.000)   10.051
  10.827   (  -0.057   -0.275    0.000)    0.281
  12.385   (  -9.005   14.908    0.000)   17.417
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.569   (   5.037    5.037    0.000)    7.124
   1.660   (   5.159    5.159    0.000)    7.295
   2.084   (   0.051    0.051    0.000)    0.072
   2.990   (   8.800    8.800    0.000)   12.446
   3.029   (  -6.764   -6.764    0.000)    9.565
   3.714   (   4.831    4.831    0.000)    6.832
   4.088   (   4.989    4.989    0.000)    7.055
   4.110   (   2.503    2.503    0.000)    3.540
   5.040   (  -5.579   -5.579    0.000)    7.890
   5.271   (  23.546   23.546    0.000)   33.298
   6.886   (  22.149   22.149    0.000)   31.323
   8.050   (  11.303   11.303    0.000)   15.985
  10.518   (  -4.345   -4.345    0.000)    6.145
  10.828   (  -0.242   -0.242    0.000)    0.342
  12.929   (  -4.601   -4.601    0.000)    6.507
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.650   (   2.928    1.600    0.000)    3.336
   1.751   (   3.702    3.456    0.000)    5.065
   2.084   (   0.032   -4.186    0.000)    4.187
   2.919   (  -4.009  -12.997    0.000)   13.601
   3.159   (   7.339   12.792    0.000)   14.748
   3.757   (   0.539    8.969    0.000)    8.985
   4.164   (   2.176   -4.157    0.000)    4.692
   4.191   (   4.936    3.793    0.000)    6.225
   4.922   (  -5.777   -6.744    0.000)    8.880
   5.707   (  19.026   24.015    0.000)   30.638
   7.017   (  -5.189   44.645    0.000)   44.946
   8.566   (  33.566  -14.048    0.000)   36.387
  10.437   (  -4.085   -6.530    0.000)    7.703
  10.824   (  -0.184   -0.277    0.000)    0.333
  12.787   (  -8.615    3.212    0.000)    9.194
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.689   (   0.935   -0.125    0.000)    0.943
   1.803   (   1.332    2.708    0.000)    3.018
   2.086   (   0.104   -6.443    0.000)    6.444
   2.866   (  -1.320  -16.406    0.000)   16.459
   3.266   (   2.865   15.422    0.000)   15.686
   3.760   (  -0.023    9.878    0.000)    9.878
   4.191   (   0.594   -7.018    0.000)    7.044
   4.272   (   2.496    3.231    0.000)    4.083
   4.826   (  -2.908   -8.010    0.000)    8.521
   6.011   (   9.345   22.708    0.000)   24.556
   6.865   (  -6.313   48.640    0.000)   49.048
   9.141   (  17.936  -19.763    0.000)   26.689
  10.352   (  -3.492   -7.431    0.000)    8.211
  10.821   (  -0.068   -0.299    0.000)    0.306
  12.633   (  -4.950    8.992    0.000)   10.264
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.675   (   0.859    0.859    0.000)    1.215
   1.815   (   2.746    2.746    0.000)    3.883
   2.013   (  -2.719   -2.719    0.000)    3.845
   2.702   (  -8.267   -8.267    0.000)   11.691
   3.409   (  11.206   11.206    0.000)   15.848
   3.866   (   2.528    2.528    0.000)    3.575
   4.090   (  -2.874   -2.874    0.000)    4.065
   4.267   (   3.637    3.637    0.000)    5.143
   4.775   (  -7.428   -7.428    0.000)   10.505
   6.146   (  18.738   18.738    0.000)   26.499
   7.772   (  20.386   20.386    0.000)   28.830
   8.404   (   5.689    5.689    0.000)    8.045
  10.315   (  -5.707   -5.707    0.000)    8.071
  10.819   (  -0.211   -0.211    0.000)    0.298
  12.795   (  -1.953   -1.953    0.000)    2.762
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.686   (   0.262   -0.144    0.000)    0.299
   1.856   (   1.093    2.359    0.000)    2.600
   1.977   (  -0.890   -4.233    0.000)    4.326
   2.588   (  -2.878  -10.746    0.000)   11.125
   3.580   (   4.768   14.923    0.000)   15.666
   3.889   (   0.248    3.553    0.000)    3.562
   4.051   (  -0.979   -6.227    0.000)    6.303
   4.331   (   2.269    2.484    0.000)    3.364
   4.645   (  -4.302   -9.599    0.000)   10.519
   6.432   (   8.152   18.212    0.000)   19.953
   7.810   (  -4.313   41.907    0.000)   42.128
   8.758   (  15.297  -17.134    0.000)   22.969
  10.211   (  -3.729   -6.691    0.000)    7.660
  10.816   (  -0.078   -0.227    0.000)    0.240
  12.744   (  -1.974    2.340    0.000)    3.061
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.684   (  -0.062   -0.062    0.000)    0.087
   1.892   (   1.055    1.055    0.000)    1.491
   1.918   (  -1.457   -1.457    0.000)    2.060
   2.438   (  -3.858   -3.858    0.000)    5.455
   3.832   (   8.323    8.323    0.000)   11.770
   3.927   (   0.658    0.658    0.000)    0.931
   3.959   (  -2.390   -2.390    0.000)    3.380
   4.369   (   1.316    1.316    0.000)    1.862
   4.465   (  -6.792   -6.792    0.000)    9.606
   6.709   (   7.851    7.851    0.000)   11.104
   8.423   (   9.844    9.844    0.000)   13.922
   8.520   (   0.762    0.762    0.000)    1.077
  10.090   (  -4.212   -4.212    0.000)    5.956
  10.813   (  -0.085   -0.085    0.000)    0.120
  12.757   (  -0.254   -0.254    0.000)    0.359
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.932   (  12.825   12.825   12.825)   22.213
   0.932   (  12.825   12.825   12.825)   22.213
   1.422   (  18.275   18.275   18.275)   31.653
   2.799   (   5.279    5.279    5.279)    9.144
   2.799   (   5.279    5.279    5.279)    9.144
   3.416   (   6.820    6.820    6.820)   11.813
   3.575   (  -2.828   -2.828   -2.828)    4.898
   3.839   (   5.660    5.660    5.660)    9.803
   3.839   (   5.660    5.660    5.660)    9.803
   5.516   (   4.988    4.988    4.988)    8.640
   5.516   (   4.988    4.988    4.988)    8.640
   7.151   (   6.296    6.296    6.296)   10.905
  10.798   (  -1.449   -1.449   -1.449)    2.510
  10.798   (  -1.449   -1.449   -1.449)    2.510
  13.328   (  -2.899   -2.899   -2.899)    5.021
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.169   (  12.843    9.435    9.435)   18.520
   1.207   (  12.156    9.031    9.031)   17.632
   1.804   (  16.533    6.631    6.631)   19.007
   2.821   (   1.468    7.382    7.382)   10.542
   2.908   (   3.376    9.507    9.507)   13.863
   3.541   (  -0.627   -5.652   -5.652)    8.018
   3.728   (  12.623   -0.050   -0.050)   12.623
   4.006   (  18.827    5.188    5.188)   20.206
   4.035   (  10.412    6.406    6.406)   13.802
   5.459   (  -3.667    8.103    8.103)   12.032
   5.684   (   4.287   11.317   11.317)   16.569
   7.443   (  23.721   -0.163   -0.163)   23.722
  10.754   (  -1.369   -2.644   -2.644)    3.982
  10.799   (   0.087   -2.032   -2.032)    2.875
  13.206   (  -9.333   -0.069   -0.069)    9.333
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.421   (  11.420    5.450    5.450)   13.778
   1.435   (  10.159    5.321    5.321)   12.642
   2.068   (   9.924   -0.412   -0.412)    9.941
   2.856   (   1.835    7.558    7.558)   10.845
   2.995   (   5.224   13.203   13.203)   19.389
   3.544   (   0.781   -8.539   -8.539)   12.101
   3.848   (   2.849   -1.788   -1.788)    3.809
   4.249   (  10.515    8.293    8.293)   15.751
   4.485   (  22.833    3.114    3.114)   23.254
   5.370   (  -4.932    7.579    7.579)   11.799
   5.699   (  -2.402   16.072   16.072)   22.856
   8.107   (  40.367   -4.674   -4.674)   40.905
  10.734   (  -0.838   -3.909   -3.909)    5.591
  10.801   (   0.095   -1.813   -1.813)    2.566
  12.948   ( -16.054    3.283    3.283)   16.712
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.607   (   6.549    2.166    2.166)    7.230
   1.616   (   7.590    2.810    2.810)    8.567
   2.223   (   5.643   -4.267   -4.267)    8.261
   2.891   (   1.510    8.095    8.095)   11.547
   3.103   (   4.904   16.926   16.926)   24.434
   3.565   (   1.119  -11.017  -11.017)   15.621
   3.885   (   1.033   -6.110   -6.110)    8.703
   4.448   (   8.746    9.862    9.862)   16.463
   4.880   (  16.640    5.310    5.310)   18.256
   5.266   (  -5.143    6.866    6.866)   10.987
   5.600   (  -7.129   15.472   15.472)   23.013
   8.983   (  43.133   -7.160   -7.160)   44.306
  10.711   (  -1.708   -4.968   -4.968)    7.230
  10.803   (   0.073   -1.613   -1.613)    2.283
  12.571   ( -20.177    7.351    7.351)   22.697
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.696   (   2.193    0.346    0.346)    2.247
   1.721   (   2.638    1.472    1.472)    3.361
   2.301   (   1.957   -5.775   -5.775)    8.398
   2.913   (   0.577    8.572    8.572)   12.136
   3.176   (   1.977   19.634   19.634)   27.836
   3.583   (   0.507  -12.481  -12.481)   17.659
   3.896   (   0.153   -9.650   -9.650)   13.648
   4.583   (   3.894   11.125   11.125)   16.208
   5.179   (  13.159    8.875    8.875)   18.184
   5.180   (  -2.660    6.018    6.018)    8.917
   5.419   ( -10.485   12.586   12.586)   20.657
   9.707   (  23.276   -9.643   -9.643)   26.976
  10.657   (  -3.157   -6.449   -6.449)    9.652
  10.804   (   0.027   -1.495   -1.495)    2.115
  12.213   ( -11.831   12.103   12.103)   20.807
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.344   (   8.946    8.946    8.655)   15.330
   1.407   (   8.882    8.882    3.751)   13.109
   1.949   (   6.594    6.594    2.123)    9.564
   2.902   (   3.747    3.747   12.919)   13.963
   3.103   (   4.856    4.856   -0.222)    6.871
   3.457   (  -2.254   -2.254   -8.109)    8.713
   3.890   (  10.672   10.672    6.352)   16.375
   4.049   (   6.215    6.215   12.568)   15.336
   4.331   (  18.875   18.875   -1.176)   26.718
   5.477   (  -2.663   -2.663   21.613)   21.939
   6.042   (  17.307   17.307   -2.347)   24.588
   7.526   (  11.249   11.249   -1.626)   15.991
  10.679   (  -3.488   -3.488    0.258)    4.940
  10.785   (  -0.127   -0.127   -4.500)    4.504
  13.162   (  -4.757   -4.757    1.137)    6.823
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.519   (   7.780    5.081    5.434)   10.764
   1.568   (   7.056    6.122    2.346)    9.632
   2.058   (   4.318   -0.583   -1.775)    4.704
   2.981   (   3.675    6.586   10.095)   12.601
   3.169   (   2.963   -3.578    2.317)    5.191
   3.441   (   0.759    2.040   -9.545)    9.790
   3.952   (  -0.942   11.387    4.256)   12.193
   4.298   (  13.242   -0.998   14.678)   19.794
   4.761   (  21.629   20.822   -1.628)   30.067
   5.401   (  -4.468   -1.725   21.994)   22.509
   6.239   (   1.415   32.956   -0.424)   32.989
   7.990   (  34.029   -5.697   -4.966)   34.858
  10.623   (  -2.253   -5.784   -0.336)    6.217
  10.784   (  -0.007   -0.351   -4.293)    4.307
  13.003   ( -10.894    1.353    3.494)   11.521
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.644   (   4.479    1.948    2.795)    5.627
   1.692   (   5.005    3.825    1.483)    6.471
   2.131   (   2.976   -4.535   -3.704)    6.568
   3.046   (   2.761    7.598    6.575)   10.420
   3.217   (   1.254  -13.318    1.680)   13.482
   3.497   (   4.754    9.311   -2.829)   10.830
   3.905   (  -3.531    7.576   -2.593)    8.752
   4.534   (  10.036   -0.223   16.961)   19.709
   5.134   (  11.762   12.867    3.705)   17.823
   5.332   (   1.124    6.592   16.226)   17.550
   6.149   (  -8.697   36.261    0.082)   37.290
   8.778   (  39.772  -12.585   -6.861)   42.276
  10.574   (  -2.981   -7.424   -1.444)    8.129
  10.784   (  -0.015   -0.552   -3.914)    3.952
  12.738   ( -14.357    7.928    6.524)   17.650
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.702   (   1.346    0.318    1.034)    1.727
   1.762   (   1.815    2.444    1.237)    3.287
   2.175   (   1.197   -6.286   -4.248)    7.681
   3.090   (   1.331    9.227    5.587)   10.869
   3.223   (  -0.231  -18.838   -2.086)   18.954
   3.606   (   4.675   16.788   12.162)   21.251
   3.822   (  -3.774    0.976  -15.503)   15.986
   4.692   (   4.740    0.782   19.070)   19.665
   5.152   (  -2.604   -4.111   18.303)   18.939
   5.543   (  12.033   21.656    1.460)   24.817
   5.940   ( -10.143   37.388   -0.973)   38.751
   9.437   (  20.569  -16.348   -8.339)   27.566
  10.498   (  -3.728   -8.215   -2.419)    9.340
  10.784   (  -0.017   -0.693   -3.652)    3.717
  12.490   (  -7.850   13.626    9.457)   18.350
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.616   (   4.227    4.227    4.017)    7.202
   1.679   (   4.890    4.890    1.816)    7.150
   2.048   (  -0.303   -0.303   -2.845)    2.877
   3.022   (  -7.053   -7.053   -0.707)   10.000
   3.131   (   7.270    7.270    7.787)   12.897
   3.500   (   2.918    2.918  -14.715)   15.283
   4.156   (   3.159    3.159   10.150)   11.089
   4.312   (   6.192    6.192   15.870)   18.126
   5.156   (  11.322   11.322    3.876)   16.475
   5.394   (   6.427    6.427   17.514)   19.732
   6.864   (  22.247   22.247   -2.201)   31.539
   7.989   (  10.759   10.759   -5.937)   16.333
  10.520   (  -4.424   -4.424    0.156)    6.259
  10.779   (  -0.228   -0.228   -4.561)    4.573
  12.980   (  -3.854   -3.854    4.596)    7.130
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.682   (   2.187    1.582    2.815)    3.900
   1.765   (   3.519    3.395    1.361)    5.076
   2.048   (   0.326   -3.492   -2.814)    4.497
   2.907   (  -4.184  -12.917   -1.203)   13.631
   3.255   (   4.341   10.010    2.468)   11.187
   3.584   (   5.759    0.852  -15.116)   16.198
   4.103   (  -5.394    8.920   11.437)   15.475
   4.504   (   9.821   -1.563   17.901)   20.478
   5.149   (  -5.327   -4.945   19.731)   21.027
   5.704   (  17.223   22.085    0.290)   28.009
   6.993   (  -5.277   44.553   -2.256)   44.921
   8.501   (  33.560  -13.742   -6.294)   36.807
  10.434   (  -4.415   -6.325   -0.315)    7.720
  10.774   (  -0.202   -0.435   -4.567)    4.592
  12.856   (  -7.642    3.258    6.308)   10.431
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.708   (   0.577    0.269    1.806)    1.915
   1.815   (   1.289    2.664    1.148)    3.174
   2.057   (   0.324   -5.213   -2.268)    5.695
   2.852   (  -1.358  -16.224   -1.424)   16.343
   3.306   (   0.986   10.763   -2.269)   11.043
   3.717   (   6.005   -1.253  -18.906)   19.876
   3.982   (  -5.464   10.295   16.915)   20.541
   4.670   (   5.521   -2.036   20.483)   21.312
   5.038   (  -4.288   -6.884   20.170)   21.739
   5.985   (   8.731   21.106   -2.403)   22.967
   6.843   (  -6.153   48.717   -2.144)   49.150
   9.074   (  17.783  -18.685   -6.458)   26.591
  10.343   (  -3.648   -7.065   -0.876)    8.000
  10.771   (  -0.095   -0.566   -4.570)    4.606
  12.720   (  -4.394    8.401    7.877)   12.327
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.705   (   0.704    0.704    2.771)    2.944
   1.827   (   2.648    2.648    1.182)    3.927
   1.993   (  -1.922   -1.922   -1.551)    3.129
   2.690   (  -8.272   -8.272   -1.179)   11.758
   3.445   (   7.231    7.231   -4.921)   11.349
   3.594   (   1.583    1.583  -18.918)   19.050
   4.197   (  -0.017   -0.017   19.424)   19.424
   4.521   (   3.946    3.946   17.228)   18.109
   5.022   (  -7.087   -7.087   21.316)   23.554
   6.111   (  17.397   17.397   -3.292)   24.822
   7.753   (  20.482   20.482   -1.769)   29.020
   8.344   (   6.344    6.344   -5.669)   10.613
  10.313   (  -5.834   -5.834   -0.276)    8.255
  10.766   (  -0.362   -0.362   -4.921)    4.948
  12.868   (  -1.634   -1.634    6.721)    7.107
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.712   (   0.127    0.133    2.487)    2.494
   1.866   (   1.048    2.295    1.034)    2.727
   1.969   (  -0.513   -3.322   -0.472)    3.394
   2.577   (  -2.869  -10.668   -1.160)   11.108
   3.507   (   0.433    8.494  -12.554)   15.164
   3.672   (   3.846   -2.574  -21.181)   21.680
   4.155   (  -2.312    7.114   26.477)   27.514
   4.632   (   5.320   -1.358   18.885)   19.667
   4.881   (  -5.792   -8.293   20.762)   23.095
   6.375   (   7.430   16.734   -5.541)   19.129
   7.791   (  -4.089   42.065   -1.844)   42.304
   8.714   (  15.598  -15.920   -4.173)   22.675
  10.205   (  -3.838   -6.703   -0.568)    7.745
  10.761   (  -0.155   -0.453   -5.149)    5.171
  12.824   (  -1.764    2.211    7.295)    7.824
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.714   (   0.025    0.025    2.842)    2.843
   1.902   (   1.018    1.018    0.927)    1.713
   1.924   (  -1.112   -1.112    0.605)    1.685
   2.428   (  -3.820   -3.820   -1.029)    5.499
   3.625   (   1.443    1.443  -19.814)   19.918
   3.633   (   0.427    0.427  -20.927)   20.935
   4.273   (   3.951    3.951   33.822)   34.281
   4.628   (   1.321    1.321   17.480)   17.580
   4.722   (  -6.612   -6.612   19.666)   21.776
   6.627   (   7.030    7.030   -7.971)   12.743
   8.408   (   9.883    9.883   -1.494)   14.056
   8.497   (   1.522    1.522   -2.227)    3.097
  10.083   (  -4.278   -4.278   -0.686)    6.089
  10.754   (  -0.186   -0.186   -5.507)    5.513
  12.836   (  -0.231   -0.231    7.259)    7.266
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.509   (   6.542    6.542    6.542)   11.331
   1.509   (   6.542    6.542    6.542)   11.331
   1.990   (   1.692    1.692    1.692)    2.931
   3.096   (   2.761    2.761    2.761)    4.782
   3.096   (   2.761    2.761    2.761)    4.782
   3.312   (  -5.783   -5.783   -5.783)   10.017
   4.223   (  18.220   18.220   18.220)   31.558
   4.324   (   9.949    9.949    9.949)   17.232
   4.324   (   9.949    9.949    9.949)   17.232
   5.973   (  10.014   10.014   10.014)   17.344
   5.973   (  10.014   10.014   10.014)   17.344
   7.536   (   4.998    4.998    4.998)    8.656
  10.686   (  -2.136   -2.136   -2.136)    3.699
  10.686   (  -2.136   -2.136   -2.136)    3.699
  13.162   (  -1.824   -1.824   -1.824)    3.159
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.632   (   4.871    4.469    4.469)    7.980
   1.633   (   5.968    4.682    4.682)    8.914
   2.019   (   1.403   -1.924   -1.924)    3.061
   3.084   (  -0.488    1.524    1.524)    2.209
   3.209   (   5.865   -2.943   -2.943)    7.192
   3.224   (  -3.044   -6.905   -6.905)   10.228
   4.242   (  -2.914   19.807   19.807)   28.162
   4.599   (  12.948    8.986    8.986)   18.143
   4.792   (  24.610   10.295   10.295)   28.594
   5.870   (  -5.687   17.141   17.141)   24.900
   6.308   (   5.871   13.377   13.377)   19.808
   7.862   (  28.739   -6.695   -6.695)   30.259
  10.601   (  -3.807   -2.366   -2.366)    5.068
  10.696   (   0.472   -3.332   -3.332)    4.736
  13.073   (  -7.022    2.590    2.590)    7.919
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.703   (   2.192    2.694    2.694)    4.396
   1.739   (   4.257    3.064    3.064)    6.075
   2.054   (   1.955   -3.659   -3.659)    5.531
   3.081   (   0.207   -3.977   -3.977)    5.628
   3.179   (  -1.664   -4.598   -4.598)    6.712
   3.332   (   5.895   -7.340   -7.340)   11.937
   4.171   (  -3.766   20.274   20.274)   28.918
   4.837   (  10.192    9.268    9.268)   16.604
   5.241   (  20.036   10.794   10.794)   25.189
   5.751   (  -5.765   17.338   17.338)   25.189
   6.247   (  -9.435   16.312   16.312)   24.923
   8.584   (  37.672  -11.860  -11.860)   41.237
  10.521   (  -4.449   -3.920   -3.920)    7.108
  10.705   (   0.360   -3.251   -3.251)    4.611
  12.891   ( -10.190    7.492    7.492)   14.700
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.728   (   0.490    1.407    1.407)    2.049
   1.798   (   1.524    2.239    2.239)    3.514
   2.087   (   1.006   -4.154   -4.154)    5.960
   3.089   (   0.435  -10.004  -10.004)   14.154
   3.153   (  -0.776   -0.336   -0.336)    0.910
   3.426   (   2.751  -11.089  -11.089)   15.921
   4.107   (  -2.009   22.048   22.048)   31.246
   4.992   (   4.405    8.898    8.898)   13.332
   5.609   (  15.536   10.449   10.449)   21.441
   5.656   (  -2.927   17.890   17.890)   25.468
   6.002   ( -12.989   16.395   16.395)   26.576
   9.204   (  19.013  -14.044  -14.044)   27.495
  10.424   (  -4.096   -4.727   -4.727)    7.840
  10.710   (   0.134   -3.202   -3.202)    4.531
  12.715   (  -5.503   11.417   11.417)   17.058
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.709   (   3.341    3.341    4.748)    6.699
   1.729   (   4.200    4.200    3.052)    6.678
   1.988   (  -0.846   -0.846   -2.853)    3.093
   2.996   (  -7.889   -7.889   -2.184)   11.369
   3.146   (   1.532    1.532   -4.394)    4.899
   3.216   (   0.290    0.290  -12.320)   12.327
   4.580   (   5.878    5.878   26.032)   27.327
   4.638   (   6.923    6.923   15.653)   18.463
   5.173   (  20.620   20.620    0.410)   29.164
   5.955   (  -2.779   -2.779   35.390)   35.607
   6.800   (  22.556   22.556   -3.951)   32.143
   7.826   (   8.761    8.761   -9.329)   15.509
  10.524   (  -4.634   -4.634    0.264)    6.559
  10.663   (  -0.522   -0.522   -6.413)    6.455
  13.085   (  -1.731   -1.731    4.891)    5.469
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.751   (   0.865    2.007    3.804)    4.387
   1.804   (   3.073    3.127    2.380)    4.989
   1.991   (   1.006   -2.322   -2.595)    3.625
   2.867   (  -4.729  -12.325   -3.035)   13.545
   3.166   (   0.157    1.391   -8.237)    8.355
   3.248   (   2.577   -2.554  -16.088)   16.492
   4.538   (  -3.300   15.275   23.531)   28.248
   4.885   (  11.144   -1.772   18.647)   21.795
   5.530   (  13.123   12.402    8.625)   20.011
   5.886   (  -2.176    2.580   27.828)   28.032
   6.931   (  -5.562   43.642   -3.613)   44.143
   8.325   (  33.680  -12.722  -10.571)   37.522
  10.423   (  -5.434   -5.458   -0.818)    7.746
  10.658   (  -0.080   -1.519   -6.390)    6.569
  13.008   (  -5.257    3.603    7.708)   10.001
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.756   (  -0.045    1.356    2.779)    3.093
   1.848   (   1.142    2.509    2.015)    3.415
   2.013   (   0.752   -3.055   -1.836)    3.643
   2.804   (  -1.521  -14.925   -3.510)   15.407
   3.164   (  -0.096    0.862   -9.171)    9.211
   3.294   (   1.409   -3.868  -19.739)   20.163
   4.484   (  -1.551   15.382   24.333)   28.829
   5.055   (   4.905   -1.077   16.503)   17.250
   5.642   (  -0.159   -4.500   38.275)   38.539
   5.934   (   4.722   14.243   -0.410)   15.011
   6.782   (  -5.786   48.711   -3.638)   49.188
   8.896   (  17.490  -15.649  -10.610)   25.755
  10.315   (  -4.083   -5.836   -1.884)    7.368
  10.658   (  -0.013   -2.028   -6.139)    6.465
  12.911   (  -3.162    7.329    9.934)   12.744
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.778   (   0.528    0.528    4.152)    4.219
   1.861   (   2.393    2.393    2.038)    3.951
   1.964   (  -0.411   -0.411   -1.175)    1.311
   2.655   (  -8.176   -8.176   -2.298)   11.788
   3.198   (   0.685    0.685  -15.395)   15.426
   3.211   (  -0.393   -0.393  -18.126)   18.135
   4.769   (   3.964    3.964   28.299)   28.849
   4.861   (   4.038    4.038   15.742)   16.746
   5.603   (  -3.400   -3.400   34.947)   35.276
   6.037   (  10.853   10.853   -1.202)   15.396
   7.703   (  20.768   20.768   -3.085)   29.532
   8.188   (   8.288    8.288   -9.295)   14.959
  10.305   (  -6.173   -6.173   -0.463)    8.743
  10.636   (  -0.818   -0.818   -7.595)    7.683
  13.033   (  -0.805   -0.805    8.559)    8.634
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.780   (  -0.082    0.876    3.924)    4.021
   1.896   (   0.933    2.120    1.803)    2.935
   1.965   (   0.199   -1.696    0.113)    1.712
   2.543   (  -2.825  -10.318   -2.144)   10.910
   3.184   (  -0.539    0.917  -17.065)   17.098
   3.232   (   0.816   -2.389  -21.147)   21.297
   4.752   (  -1.297   11.433   24.691)   27.241
   5.005   (   5.323   -2.930   17.249)   18.288
   5.510   (  -3.771   -7.655   43.398)   44.229
   6.216   (   5.200   12.348   -9.716)   16.550
   7.739   (  -3.466   42.544   -3.220)   42.806
   8.601   (  16.377  -12.663   -6.596)   21.727
  10.189   (  -4.164   -6.655   -1.023)    7.917
  10.624   (  -0.289   -1.289   -7.917)    8.026
  13.003   (  -1.276    2.015    9.426)    9.723
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.791   (   0.260    0.260    4.512)    4.527
   1.928   (   0.922    0.922    1.624)    2.083
   1.943   (  -0.490   -0.490    1.261)    1.439
   2.399   (  -3.691   -3.691   -1.775)    5.514
   3.197   (  -0.210   -0.210  -21.295)   21.297
   3.198   (  -0.330   -0.330  -21.094)   21.099
   4.934   (   3.837    3.837   21.349)   22.028
   4.968   (   1.293    1.293   15.482)   15.589
   5.363   (  -5.506   -5.506   49.485)   50.093
   6.399   (   4.974    4.974  -14.049)   15.712
   8.366   (   9.997    9.997   -2.552)   14.367
   8.436   (   3.375    3.375   -3.579)    5.966
  10.064   (  -4.457   -4.457   -1.155)    6.408
  10.606   (  -0.477   -0.477   -8.707)    8.733
  13.015   (  -0.161   -0.161    9.467)    9.470
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.796   (   3.347    3.347    3.347)    5.798
   1.796   (   3.347    3.347    3.347)    5.798
   1.941   (  -1.412   -1.412   -1.412)    2.445
   2.902   (  -7.247   -7.247   -7.247)   12.553
   3.049   (  -4.022   -4.022   -4.022)    6.966
   3.049   (  -4.022   -4.022   -4.022)    6.966
   4.932   (   8.983    8.983    8.983)   15.560
   4.932   (   8.983    8.983    8.983)   15.560
   5.300   (  14.771   14.771   14.771)   25.584
   6.710   (  13.805   13.805   13.805)   23.910
   6.710   (  13.805   13.805   13.805)   23.910
   7.691   (   0.915    0.915    0.915)    1.586
  10.530   (  -3.177   -3.177   -3.177)    5.502
  10.530   (  -3.177   -3.177   -3.177)    5.502
  13.151   (   1.294    1.294    1.294)    2.242
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.825   (  -0.251    3.186    3.186)    4.513
   1.856   (   2.524    2.665    2.665)    4.536
   1.949   (   1.901   -1.397   -1.397)    2.742
   2.763   (  -5.562   -7.639   -7.639)   12.150
   2.991   (  -1.635   -8.309   -8.309)   11.864
   3.004   (  -1.281   -6.866   -6.866)    9.793
   4.923   (  -2.564   14.897   14.897)   21.223
   5.160   (   9.099    6.355    6.355)   12.789
   5.635   (  15.666    7.166    7.166)   18.659
   6.615   (  -4.221   23.893   23.893)   34.053
   6.943   (  -5.507   18.184   18.184)   26.299
   8.098   (  33.319  -10.932  -10.932)   36.731
  10.392   (  -7.017   -2.704   -2.704)    7.991
  10.535   (   0.221   -5.160   -5.160)    7.301
  13.137   (  -2.361    4.448    4.448)    6.719
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.814   (  -0.430    2.756    2.756)    3.921
   1.892   (   0.937    2.285    2.285)    3.365
   1.985   (   1.098   -0.966   -0.966)    1.752
   2.687   (  -1.864   -8.575   -8.575)   12.270
   2.974   (  -0.296  -10.747  -10.747)   15.201
   2.987   (  -0.417   -7.754   -7.754)   10.974
   4.878   (  -1.364   15.087   15.087)   21.380
   5.289   (   3.358    5.712    5.712)    8.749
   5.886   (   7.683    2.277    2.277)    8.330
   6.552   (  -1.764   24.614   24.614)   34.854
   6.783   (  -6.341   21.698   21.698)   31.333
   8.668   (  17.339  -11.259  -11.259)   23.541
  10.262   (  -4.604   -3.559   -3.559)    6.821
  10.539   (   0.082   -5.309   -5.309)    7.508
  13.085   (  -1.834    6.649    6.649)    9.580
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.864   (   0.439    0.439    4.148)    4.194
   1.906   (   2.080    2.080    2.293)    3.729
   1.947   (   1.196    1.196   -0.567)    1.784
   2.596   (  -7.428   -7.428   -3.640)   11.117
   2.885   (  -3.369   -3.369  -13.462)   14.280
   2.902   (  -3.259   -3.259  -13.017)   13.809
   5.147   (   3.741    3.741   12.270)   13.362
   5.195   (   3.184    3.184   14.082)   14.784
   5.826   (   9.906    9.906   -4.615)   14.749
   6.645   (  -4.973   -4.973   49.943)   50.435
   7.634   (  21.214   21.214   -3.552)   30.211
   7.991   (  11.019   11.019   -9.481)   18.241
  10.295   (  -6.619   -6.619   -0.501)    9.373
  10.469   (  -1.599   -1.599   -8.716)    9.004
  13.182   (   0.308    0.308    5.619)    5.635
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.862   (  -0.200    1.832    3.972)    4.379
   1.936   (   0.791    1.892    2.056)    2.904
   1.973   (   0.881   -0.220    0.636)    1.109
   2.491   (  -2.690   -9.107   -2.981)    9.953
   2.841   (  -1.002   -3.729  -16.722)   17.162
   2.861   (  -0.980   -4.666  -14.429)   15.196
   5.120   (  -1.485    9.122   13.486)   16.349
   5.304   (   3.268   -2.549   11.805)   12.511
   5.983   (   4.568    6.088  -10.681)   13.116
   6.555   (  -2.924   -4.341   53.511)   53.766
   7.667   (  -2.615   43.188   -3.603)   43.417
   8.463   (  17.316   -8.467   -6.644)   20.389
  10.164   (  -4.742   -6.144   -1.463)    7.898
  10.452   (  -0.326   -3.098   -8.858)    9.390
  13.171   (  -0.714    2.008    6.550)    6.888
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.886   (   0.547    0.547    4.588)    4.652
   1.965   (   0.798    0.798    1.871)    2.184
   1.973   (   0.083    0.083    1.555)    1.559
   2.359   (  -3.446   -3.446   -2.007)    5.271
   2.794   (  -1.104   -1.104  -17.991)   18.059
   2.799   (  -1.518   -1.518  -17.883)   18.011
   5.244   (   1.153    1.153   11.621)   11.735
   5.253   (   0.306    0.306   12.257)   12.265
   6.082   (   2.855    2.855  -16.099)   16.598
   6.485   (  -2.075   -2.075   56.264)   56.340
   8.309   (  10.172   10.172   -2.854)   14.665
   8.360   (   5.448    5.448   -3.744)    8.566
  10.039   (  -4.694   -4.694   -1.259)    6.757
  10.412   (  -0.920   -0.920  -10.248)   10.330
  13.185   (  -0.041   -0.041    6.723)    6.723
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.930   (   1.075    1.075    1.075)    1.863
   1.948   (   1.811    1.811    1.811)    3.137
   1.948   (   1.811    1.811    1.811)    3.137
   2.503   (  -5.329   -5.329   -5.329)    9.230
   2.687   (  -7.016   -7.016   -7.016)   12.151
   2.687   (  -7.016   -7.016   -7.016)   12.151
   5.352   (   4.793    4.793    4.793)    8.302
   5.352   (   4.793    4.793    4.793)    8.302
   5.769   (   0.242    0.242    0.242)    0.419
   7.567   (  13.517   13.517   13.517)   23.413
   7.567   (  13.517   13.517   13.517)   23.413
   7.874   (   6.163    6.163    6.163)   10.674
  10.286   (  -4.812   -4.812   -4.812)    8.335
  10.286   (  -4.812   -4.812   -4.812)    8.335
  13.252   (   1.444    1.444    1.444)    2.501
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.932   (  -0.051    2.757    2.757)    3.899
   1.974   (   0.686    1.667    1.667)    2.455
   1.988   (   1.178    0.789    0.789)    1.623
   2.415   (  -2.510   -4.803   -4.803)    7.242
   2.578   (  -2.671   -8.468   -8.468)   12.270
   2.618   (  -1.440   -9.325   -9.325)   13.266
   5.334   (  -1.463    7.820    7.820)   11.156
   5.460   (   2.736    2.876    2.876)    4.901
   5.803   (   1.619   -5.208   -5.208)    7.541
   7.525   (  -1.127   22.084   22.084)   31.252
   7.652   (  -2.055   19.572   19.572)   27.755
   8.345   (  18.028   -4.753   -4.753)   19.241
  10.120   (  -5.613   -3.447   -3.447)    7.434
  10.280   (  -0.142   -7.271   -7.271)   10.284
  13.259   (  -0.209    2.215    2.215)    3.140
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.968   (   0.775    0.775    3.319)    3.495
   2.000   (   0.685    0.685    1.503)    1.788
   2.002   (   0.520    0.520    1.273)    1.470
   2.321   (  -3.029   -3.029   -1.658)    4.593
   2.485   (  -2.181   -2.181  -12.053)   12.441
   2.492   (  -2.692   -2.692  -12.025)   12.613
   5.435   (   0.966    0.966    7.198)    7.326
   5.450   (  -0.237   -0.237    7.190)    7.198
   5.776   (   0.567    0.567  -12.999)   13.024
   7.540   (  -1.069   -1.069   45.634)   45.659
   8.256   (  10.360   10.360   -2.251)   14.824
   8.292   (   7.050    7.050   -2.796)   10.355
  10.016   (  -4.915   -4.915   -0.963)    7.018
  10.196   (  -1.442   -1.442  -10.891)   11.081
  13.277   (   0.107    0.107    2.443)    2.448
======================= 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 2.98e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.014   (   1.060    1.060    1.060)    1.837
   2.022   (   0.596    0.596    0.596)    1.033
   2.022   (   0.596    0.596    0.596)    1.033
   2.291   (  -1.652   -1.652   -1.652)    2.861
   2.324   (  -3.748   -3.748   -3.748)    6.492
   2.324   (  -3.748   -3.748   -3.748)    6.492
   5.533   (   1.370    1.370    1.370)    2.374
   5.533   (   1.370    1.370    1.370)    2.374
   5.601   (  -2.907   -2.907   -2.907)    5.035
   8.224   (   6.705    6.705    6.705)   11.614
   8.224   (   6.705    6.705    6.705)   11.614
   8.261   (   4.362    4.362    4.362)    7.555
  10.002   (  -3.485   -3.485   -3.485)    6.036
  10.002   (  -3.485   -3.485   -3.485)    6.036
  13.299   (   0.218    0.218    0.218)    0.378
=================== 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    398.312    398.312    398.312      0.000     -0.000      0.000 3/19965
   20.0     48.039     48.039     48.039      0.000     -0.000      0.000 3/19965
   30.0     16.633     16.633     16.633      0.000     -0.000      0.000 3/19965
   40.0      9.790      9.790      9.790      0.000     -0.000      0.000 3/19965
   50.0      7.440      7.440      7.440      0.000     -0.000      0.000 3/19965
   60.0      6.440      6.440      6.440      0.000     -0.000      0.000 3/19965
   70.0      5.931      5.931      5.931      0.000     -0.000      0.000 3/19965
   80.0      5.611      5.611      5.611      0.000     -0.000      0.000 3/19965
   90.0      5.363      5.363      5.363      0.000     -0.000      0.000 3/19965
  100.0      5.143      5.143      5.143      0.000     -0.000      0.000 3/19965
  110.0      4.935      4.935      4.935      0.000     -0.000      0.000 3/19965
  120.0      4.736      4.736      4.736      0.000     -0.000      0.000 3/19965
  130.0      4.544      4.544      4.544     -0.000     -0.000      0.000 3/19965
  140.0      4.359      4.359      4.359     -0.000     -0.000      0.000 3/19965
  150.0      4.183      4.183      4.183     -0.000     -0.000      0.000 3/19965
  160.0      4.015      4.015      4.015     -0.000     -0.000      0.000 3/19965
  170.0      3.856      3.856      3.856     -0.000     -0.000      0.000 3/19965
  180.0      3.706      3.706      3.706     -0.000     -0.000      0.000 3/19965
  190.0      3.565      3.565      3.565     -0.000     -0.000      0.000 3/19965
  200.0      3.432      3.432      3.432     -0.000     -0.000      0.000 3/19965
  210.0      3.307      3.307      3.307     -0.000     -0.000      0.000 3/19965
  220.0      3.189      3.189      3.189     -0.000     -0.000      0.000 3/19965
  230.0      3.078      3.078      3.078     -0.000     -0.000      0.000 3/19965
  240.0      2.974      2.974      2.974     -0.000     -0.000      0.000 3/19965
  250.0      2.876      2.876      2.876     -0.000     -0.000      0.000 3/19965
  260.0      2.783      2.783      2.783     -0.000     -0.000      0.000 3/19965
  270.0      2.696      2.696      2.696     -0.000     -0.000      0.000 3/19965
  280.0      2.613      2.613      2.613     -0.000     -0.000      0.000 3/19965
  290.0      2.535      2.535      2.535     -0.000     -0.000      0.000 3/19965
  300.0      2.461      2.461      2.461     -0.000     -0.000      0.000 3/19965
  310.0      2.391      2.391      2.391     -0.000     -0.000      0.000 3/19965
  320.0      2.325      2.325      2.325     -0.000     -0.000      0.000 3/19965
  330.0      2.262      2.262      2.262     -0.000     -0.000      0.000 3/19965
  340.0      2.202      2.202      2.202     -0.000     -0.000      0.000 3/19965
  350.0      2.145      2.145      2.145     -0.000     -0.000      0.000 3/19965
  360.0      2.091      2.091      2.091     -0.000     -0.000      0.000 3/19965
  370.0      2.039      2.039      2.039     -0.000     -0.000      0.000 3/19965
  380.0      1.990      1.990      1.990     -0.000     -0.000      0.000 3/19965
  390.0      1.943      1.943      1.943     -0.000     -0.000      0.000 3/19965
  400.0      1.898      1.898      1.898     -0.000     -0.000      0.000 3/19965
  410.0      1.855      1.855      1.855     -0.000     -0.000      0.000 3/19965
  420.0      1.814      1.814      1.814     -0.000     -0.000      0.000 3/19965
  430.0      1.774      1.774      1.774     -0.000     -0.000      0.000 3/19965
  440.0      1.736      1.736      1.736     -0.000     -0.000      0.000 3/19965
  450.0      1.700      1.700      1.700     -0.000     -0.000      0.000 3/19965
  460.0      1.665      1.665      1.665     -0.000     -0.000      0.000 3/19965
  470.0      1.632      1.632      1.632     -0.000     -0.000      0.000 3/19965
  480.0      1.600      1.600      1.600     -0.000     -0.000      0.000 3/19965
  490.0      1.569      1.569      1.569     -0.000     -0.000      0.000 3/19965
  500.0      1.539      1.539      1.539     -0.000     -0.000      0.000 3/19965
  510.0      1.510      1.510      1.510     -0.000     -0.000      0.000 3/19965
  520.0      1.482      1.482      1.482     -0.000     -0.000      0.000 3/19965
  530.0      1.456      1.456      1.456     -0.000     -0.000      0.000 3/19965
  540.0      1.430      1.430      1.430     -0.000     -0.000      0.000 3/19965
  550.0      1.405      1.405      1.405     -0.000     -0.000      0.000 3/19965
  560.0      1.381      1.381      1.381     -0.000     -0.000      0.000 3/19965
  570.0      1.357      1.357      1.357     -0.000     -0.000      0.000 3/19965
  580.0      1.335      1.335      1.335     -0.000     -0.000      0.000 3/19965
  590.0      1.313      1.313      1.313     -0.000     -0.000      0.000 3/19965
  600.0      1.292      1.292      1.292     -0.000     -0.000      0.000 3/19965
  610.0      1.271      1.271      1.271     -0.000     -0.000      0.000 3/19965
  620.0      1.252      1.252      1.252     -0.000     -0.000      0.000 3/19965
  630.0      1.232      1.232      1.232     -0.000     -0.000      0.000 3/19965
  640.0      1.214      1.214      1.214     -0.000     -0.000      0.000 3/19965
  650.0      1.196      1.196      1.196     -0.000     -0.000      0.000 3/19965
  660.0      1.178      1.178      1.178     -0.000     -0.000      0.000 3/19965
  670.0      1.161      1.161      1.161     -0.000     -0.000      0.000 3/19965
  680.0      1.144      1.144      1.144     -0.000     -0.000      0.000 3/19965
  690.0      1.128      1.128      1.128     -0.000     -0.000      0.000 3/19965
  700.0      1.113      1.113      1.113     -0.000     -0.000      0.000 3/19965
  710.0      1.097      1.097      1.097     -0.000     -0.000      0.000 3/19965
  720.0      1.082      1.082      1.082     -0.000     -0.000      0.000 3/19965
  730.0      1.068      1.068      1.068     -0.000     -0.000      0.000 3/19965
  740.0      1.054      1.054      1.054     -0.000     -0.000      0.000 3/19965
  750.0      1.040      1.040      1.040     -0.000     -0.000      0.000 3/19965
  760.0      1.027      1.027      1.027     -0.000     -0.000      0.000 3/19965
  770.0      1.014      1.014      1.014     -0.000     -0.000      0.000 3/19965
  780.0      1.001      1.001      1.001     -0.000     -0.000      0.000 3/19965
  790.0      0.989      0.989      0.989     -0.000     -0.000      0.000 3/19965
  800.0      0.976      0.976      0.976     -0.000     -0.000      0.000 3/19965
  810.0      0.965      0.965      0.965     -0.000     -0.000      0.000 3/19965
  820.0      0.953      0.953      0.953     -0.000     -0.000      0.000 3/19965
  830.0      0.942      0.942      0.942     -0.000     -0.000      0.000 3/19965
  840.0      0.931      0.931      0.931     -0.000     -0.000      0.000 3/19965
  850.0      0.920      0.920      0.920     -0.000     -0.000      0.000 3/19965
  860.0      0.910      0.910      0.910     -0.000     -0.000      0.000 3/19965
  870.0      0.899      0.899      0.899     -0.000     -0.000      0.000 3/19965
  880.0      0.889      0.889      0.889     -0.000     -0.000      0.000 3/19965
  890.0      0.879      0.879      0.879     -0.000     -0.000      0.000 3/19965
  900.0      0.870      0.870      0.870     -0.000     -0.000      0.000 3/19965
  910.0      0.860      0.860      0.860     -0.000     -0.000      0.000 3/19965
  920.0      0.851      0.851      0.851     -0.000     -0.000      0.000 3/19965
  930.0      0.842      0.842      0.842     -0.000     -0.000      0.000 3/19965
  940.0      0.833      0.833      0.833     -0.000     -0.000      0.000 3/19965
  950.0      0.825      0.825      0.825     -0.000     -0.000      0.000 3/19965
  960.0      0.816      0.816      0.816     -0.000     -0.000      0.000 3/19965
  970.0      0.808      0.808      0.808     -0.000     -0.000      0.000 3/19965
  980.0      0.800      0.800      0.800     -0.000     -0.000      0.000 3/19965
  990.0      0.792      0.792      0.792     -0.000     -0.000      0.000 3/19965
 1000.0      0.784      0.784      0.784     -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:16:01]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

