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

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

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


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

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

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


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

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

-------------------------[time 2026-01-08 21:14:43]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [3 3 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: Pm-3m (221)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.512828610000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.512828610000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.512828610000000
Atomic positions (fractional):
    1 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
    2 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078
    3 H   0.50000000000000  0.00000000000000  0.50000000000000   1.008
    4 H   0.00000000000000  0.50000000000000  0.50000000000000   1.008
    5 H   0.50000000000000  0.50000000000000  0.00000000000000   1.008
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.538485829999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.538485829999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.538485829999999
Atomic positions (fractional):
    1 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 1
    2 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 1
    3 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 1
    4 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 1
    5 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 1
    6 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 1
    7 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 1
    8 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 1
    9 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 1
   10 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 1
   11 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 1
   12 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 1
   13 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 1
   14 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 1
   15 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 1
   16 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 1
   17 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 1
   18 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 1
   19 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 1
   20 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 1
   21 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 1
   22 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 1
   23 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 1
   24 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 1
   25 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 1
   26 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 1
   27 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 1
   28 Ca  0.16666666666667  0.16666666666667  0.16666666666667  40.078 > 28
   29 Ca  0.50000000000000  0.16666666666667  0.16666666666667  40.078 > 28
   30 Ca  0.83333333333333  0.16666666666667  0.16666666666667  40.078 > 28
   31 Ca  0.16666666666667  0.50000000000000  0.16666666666667  40.078 > 28
   32 Ca  0.50000000000000  0.50000000000000  0.16666666666667  40.078 > 28
   33 Ca  0.83333333333333  0.50000000000000  0.16666666666667  40.078 > 28
   34 Ca  0.16666666666667  0.83333333333333  0.16666666666667  40.078 > 28
   35 Ca  0.50000000000000  0.83333333333333  0.16666666666667  40.078 > 28
   36 Ca  0.83333333333333  0.83333333333333  0.16666666666667  40.078 > 28
   37 Ca  0.16666666666667  0.16666666666667  0.50000000000000  40.078 > 28
   38 Ca  0.50000000000000  0.16666666666667  0.50000000000000  40.078 > 28
   39 Ca  0.83333333333333  0.16666666666667  0.50000000000000  40.078 > 28
   40 Ca  0.16666666666667  0.50000000000000  0.50000000000000  40.078 > 28
   41 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 28
   42 Ca  0.83333333333333  0.50000000000000  0.50000000000000  40.078 > 28
   43 Ca  0.16666666666667  0.83333333333333  0.50000000000000  40.078 > 28
   44 Ca  0.50000000000000  0.83333333333333  0.50000000000000  40.078 > 28
   45 Ca  0.83333333333333  0.83333333333333  0.50000000000000  40.078 > 28
   46 Ca  0.16666666666667  0.16666666666667  0.83333333333333  40.078 > 28
   47 Ca  0.50000000000000  0.16666666666667  0.83333333333333  40.078 > 28
   48 Ca  0.83333333333333  0.16666666666667  0.83333333333333  40.078 > 28
   49 Ca  0.16666666666667  0.50000000000000  0.83333333333333  40.078 > 28
   50 Ca  0.50000000000000  0.50000000000000  0.83333333333333  40.078 > 28
   51 Ca  0.83333333333333  0.50000000000000  0.83333333333333  40.078 > 28
   52 Ca  0.16666666666667  0.83333333333333  0.83333333333333  40.078 > 28
   53 Ca  0.50000000000000  0.83333333333333  0.83333333333333  40.078 > 28
   54 Ca  0.83333333333333  0.83333333333333  0.83333333333333  40.078 > 28
   55 H   0.16666666666667  0.00000000000000  0.16666666666667   1.008 > 55
   56 H   0.50000000000000  0.00000000000000  0.16666666666667   1.008 > 55
   57 H   0.83333333333333  0.00000000000000  0.16666666666667   1.008 > 55
   58 H   0.16666666666667  0.33333333333333  0.16666666666667   1.008 > 55
   59 H   0.50000000000000  0.33333333333333  0.16666666666667   1.008 > 55
   60 H   0.83333333333333  0.33333333333333  0.16666666666667   1.008 > 55
   61 H   0.16666666666667  0.66666666666667  0.16666666666667   1.008 > 55
   62 H   0.50000000000000  0.66666666666667  0.16666666666667   1.008 > 55
   63 H   0.83333333333333  0.66666666666667  0.16666666666667   1.008 > 55
   64 H   0.16666666666667  0.00000000000000  0.50000000000000   1.008 > 55
   65 H   0.50000000000000  0.00000000000000  0.50000000000000   1.008 > 55
   66 H   0.83333333333333  0.00000000000000  0.50000000000000   1.008 > 55
   67 H   0.16666666666667  0.33333333333333  0.50000000000000   1.008 > 55
   68 H   0.50000000000000  0.33333333333333  0.50000000000000   1.008 > 55
   69 H   0.83333333333333  0.33333333333333  0.50000000000000   1.008 > 55
   70 H   0.16666666666667  0.66666666666667  0.50000000000000   1.008 > 55
   71 H   0.50000000000000  0.66666666666667  0.50000000000000   1.008 > 55
   72 H   0.83333333333333  0.66666666666667  0.50000000000000   1.008 > 55
   73 H   0.16666666666667  0.00000000000000  0.83333333333333   1.008 > 55
   74 H   0.50000000000000  0.00000000000000  0.83333333333333   1.008 > 55
   75 H   0.83333333333333  0.00000000000000  0.83333333333333   1.008 > 55
   76 H   0.16666666666667  0.33333333333333  0.83333333333333   1.008 > 55
   77 H   0.50000000000000  0.33333333333333  0.83333333333333   1.008 > 55
   78 H   0.83333333333333  0.33333333333333  0.83333333333333   1.008 > 55
   79 H   0.16666666666667  0.66666666666667  0.83333333333333   1.008 > 55
   80 H   0.50000000000000  0.66666666666667  0.83333333333333   1.008 > 55
   81 H   0.83333333333333  0.66666666666667  0.83333333333333   1.008 > 55
   82 H   0.00000000000000  0.16666666666667  0.16666666666667   1.008 > 82
   83 H   0.33333333333333  0.16666666666667  0.16666666666667   1.008 > 82
   84 H   0.66666666666667  0.16666666666667  0.16666666666667   1.008 > 82
   85 H   0.00000000000000  0.50000000000000  0.16666666666667   1.008 > 82
   86 H   0.33333333333333  0.50000000000000  0.16666666666667   1.008 > 82
   87 H   0.66666666666667  0.50000000000000  0.16666666666667   1.008 > 82
   88 H   0.00000000000000  0.83333333333333  0.16666666666667   1.008 > 82
   89 H   0.33333333333333  0.83333333333333  0.16666666666667   1.008 > 82
   90 H   0.66666666666667  0.83333333333333  0.16666666666667   1.008 > 82
   91 H   0.00000000000000  0.16666666666667  0.50000000000000   1.008 > 82
   92 H   0.33333333333333  0.16666666666667  0.50000000000000   1.008 > 82
   93 H   0.66666666666667  0.16666666666667  0.50000000000000   1.008 > 82
   94 H   0.00000000000000  0.50000000000000  0.50000000000000   1.008 > 82
   95 H   0.33333333333333  0.50000000000000  0.50000000000000   1.008 > 82
   96 H   0.66666666666667  0.50000000000000  0.50000000000000   1.008 > 82
   97 H   0.00000000000000  0.83333333333333  0.50000000000000   1.008 > 82
   98 H   0.33333333333333  0.83333333333333  0.50000000000000   1.008 > 82
   99 H   0.66666666666667  0.83333333333333  0.50000000000000   1.008 > 82
  100 H   0.00000000000000  0.16666666666667  0.83333333333333   1.008 > 82
  101 H   0.33333333333333  0.16666666666667  0.83333333333333   1.008 > 82
  102 H   0.66666666666667  0.16666666666667  0.83333333333333   1.008 > 82
  103 H   0.00000000000000  0.50000000000000  0.83333333333333   1.008 > 82
  104 H   0.33333333333333  0.50000000000000  0.83333333333333   1.008 > 82
  105 H   0.66666666666667  0.50000000000000  0.83333333333333   1.008 > 82
  106 H   0.00000000000000  0.83333333333333  0.83333333333333   1.008 > 82
  107 H   0.33333333333333  0.83333333333333  0.83333333333333   1.008 > 82
  108 H   0.66666666666667  0.83333333333333  0.83333333333333   1.008 > 82
  109 H   0.16666666666667  0.16666666666667  0.00000000000000   1.008 > 109
  110 H   0.50000000000000  0.16666666666667  0.00000000000000   1.008 > 109
  111 H   0.83333333333333  0.16666666666667  0.00000000000000   1.008 > 109
  112 H   0.16666666666667  0.50000000000000  0.00000000000000   1.008 > 109
  113 H   0.50000000000000  0.50000000000000  0.00000000000000   1.008 > 109
  114 H   0.83333333333333  0.50000000000000  0.00000000000000   1.008 > 109
  115 H   0.16666666666667  0.83333333333333  0.00000000000000   1.008 > 109
  116 H   0.50000000000000  0.83333333333333  0.00000000000000   1.008 > 109
  117 H   0.83333333333333  0.83333333333333  0.00000000000000   1.008 > 109
  118 H   0.16666666666667  0.16666666666667  0.33333333333333   1.008 > 109
  119 H   0.50000000000000  0.16666666666667  0.33333333333333   1.008 > 109
  120 H   0.83333333333333  0.16666666666667  0.33333333333333   1.008 > 109
  121 H   0.16666666666667  0.50000000000000  0.33333333333333   1.008 > 109
  122 H   0.50000000000000  0.50000000000000  0.33333333333333   1.008 > 109
  123 H   0.83333333333333  0.50000000000000  0.33333333333333   1.008 > 109
  124 H   0.16666666666667  0.83333333333333  0.33333333333333   1.008 > 109
  125 H   0.50000000000000  0.83333333333333  0.33333333333333   1.008 > 109
  126 H   0.83333333333333  0.83333333333333  0.33333333333333   1.008 > 109
  127 H   0.16666666666667  0.16666666666667  0.66666666666667   1.008 > 109
  128 H   0.50000000000000  0.16666666666667  0.66666666666667   1.008 > 109
  129 H   0.83333333333333  0.16666666666667  0.66666666666667   1.008 > 109
  130 H   0.16666666666667  0.50000000000000  0.66666666666667   1.008 > 109
  131 H   0.50000000000000  0.50000000000000  0.66666666666667   1.008 > 109
  132 H   0.83333333333333  0.50000000000000  0.66666666666667   1.008 > 109
  133 H   0.16666666666667  0.83333333333333  0.66666666666667   1.008 > 109
  134 H   0.50000000000000  0.83333333333333  0.66666666666667   1.008 > 109
  135 H   0.83333333333333  0.83333333333333  0.66666666666667   1.008 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.7365871    0.0000000    0.0000000
            0.0000000    3.7365871    0.0000000
            0.0000000    0.0000000    3.7365871
-------------------------- Born effective charges --------------------------
    1 Cs    1.1016009    0.0000000    0.0000000
            0.0000000    1.1016009    0.0000000
            0.0000000    0.0000000    1.1016009
    2 Ca    1.9271004    0.0000000    0.0000000
            0.0000000    1.9271004    0.0000000
            0.0000000    0.0000000    1.9271004
    3 H    -0.9573206    0.0000000    0.0000000
            0.0000000   -1.1140600    0.0000000
            0.0000000    0.0000000   -0.9573206
    4 H    -1.1140600    0.0000000    0.0000000
            0.0000000   -0.9573206    0.0000000
            0.0000000    0.0000000   -0.9573206
    5 H    -0.9573206    0.0000000    0.0000000
            0.0000000   -0.9573206    0.0000000
            0.0000000    0.0000000   -1.1140600
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: -0.00000000 (xxx) -0.00000000 (xxx) -0.00000000 (xxx)
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
--------------------------- 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.92, Number of G-points: 305, Lambda: 0.19
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/56) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
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
   3.500   (   0.000    0.000    0.000)    0.000
   3.500   (   0.000    0.000    0.000)    0.000
   3.500   (   0.000    0.000    0.000)    0.000
  13.322   (   0.000    0.000    0.000)    0.000
  13.322   (   0.000    0.000    0.000)    0.000
  13.322   (   0.000    0.000    0.000)    0.000
  14.075   (   0.000    0.000    0.000)    0.000
  14.075   (   0.000    0.000    0.000)    0.000
  14.075   (   0.000    0.000    0.000)    0.000
  32.391   (   0.000    0.000    0.000)    0.000
  32.391   (   0.000    0.000    0.000)    0.000
  32.391   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/56) =======================
q-point: ( 0.09  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.507   (  24.345    0.000    0.000)   24.345
   0.507   (  24.345    0.000    0.000)   24.345
   0.825   (  37.195    0.000    0.000)   37.195
   3.497   (  -0.304    0.000    0.000)    0.304
   3.497   (  -0.304    0.000    0.000)    0.304
   3.717   (  17.620    0.000    0.000)   17.620
  13.438   (  11.049    0.000    0.000)   11.049
  13.438   (  11.049    0.000    0.000)   11.049
  13.482   (  15.360    0.000    0.000)   15.360
  14.151   (   7.292    0.000    0.000)    7.292
  14.151   (   7.292    0.000    0.000)    7.292
  20.606   (   7.342    0.000    0.000)    7.342
  32.391   (   0.001    0.000    0.000)    0.001
  32.391   (   0.001    0.000    0.000)    0.001
  35.170   (   7.626    0.000    0.000)    7.626
======================= Grid point 2 (3/56) =======================
q-point: ( 0.18  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.967   (  20.999    0.000    0.000)   20.999
   0.967   (  20.999    0.000    0.000)   20.999
   1.459   (  25.654    0.000    0.000)   25.654
   3.487   (  -0.766    0.000    0.000)    0.766
   3.487   (  -0.766    0.000    0.000)    0.766
   4.199   (  28.497    0.000    0.000)   28.497
  13.738   (  17.769    0.000    0.000)   17.769
  13.738   (  17.769    0.000    0.000)   17.769
  13.902   (  25.030    0.000    0.000)   25.030
  14.356   (  12.592    0.000    0.000)   12.592
  14.356   (  12.592    0.000    0.000)   12.592
  20.804   (  11.586    0.000    0.000)   11.586
  32.391   (   0.011    0.000    0.000)    0.011
  32.391   (   0.011    0.000    0.000)    0.011
  35.386   (  13.260    0.000    0.000)   13.260
======================= Grid point 3 (4/56) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.343   (  16.195    0.000    0.000)   16.195
   1.343   (  16.195    0.000    0.000)   16.195
   1.874   (  16.121    0.000    0.000)   16.121
   3.466   (  -1.244    0.000    0.000)    1.244
   3.466   (  -1.244    0.000    0.000)    1.244
   4.788   (  28.161    0.000    0.000)   28.161
  14.109   (  17.987    0.000    0.000)   17.987
  14.109   (  17.987    0.000    0.000)   17.987
  14.433   (  26.191    0.000    0.000)   26.191
  14.634   (  14.378    0.000    0.000)   14.378
  14.634   (  14.378    0.000    0.000)   14.378
  21.044   (  11.554    0.000    0.000)   11.554
  32.391   (   0.025    0.000    0.000)    0.025
  32.391   (   0.025    0.000    0.000)    0.025
  35.678   (  14.999    0.000    0.000)   14.999
======================= Grid point 4 (5/56) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.612   (  10.343    0.000    0.000)   10.343
   1.612   (  10.343    0.000    0.000)   10.343
   2.126   (   9.094    0.000    0.000)    9.094
   3.439   (  -1.295    0.000    0.000)    1.295
   3.439   (  -1.295    0.000    0.000)    1.295
   5.278   (  19.534    0.000    0.000)   19.534
  14.425   (  12.670    0.000    0.000)   12.670
  14.425   (  12.670    0.000    0.000)   12.670
  14.903   (  19.336    0.000    0.000)   19.336
  14.904   (  11.597    0.000    0.000)   11.597
  14.904   (  11.597    0.000    0.000)   11.597
  21.246   (   8.117    0.000    0.000)    8.117
  32.392   (   0.029    0.000    0.000)    0.029
  32.392   (   0.029    0.000    0.000)    0.029
  35.956   (  11.820    0.000    0.000)   11.820
======================= Grid point 5 (6/56) =======================
q-point: ( 0.45  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.754   (   3.590    0.000    0.000)    3.590
   1.754   (   3.590    0.000    0.000)    3.590
   2.247   (   2.975    0.000    0.000)    2.975
   3.419   (  -0.575    0.000    0.000)    0.575
   3.419   (  -0.575    0.000    0.000)    0.575
   5.548   (   6.894    0.000    0.000)    6.894
  14.600   (   4.487    0.000    0.000)    4.487
  14.600   (   4.487    0.000    0.000)    4.487
  15.072   (   4.465    0.000    0.000)    4.465
  15.072   (   4.465    0.000    0.000)    4.465
  15.175   (   7.072    0.000    0.000)    7.072
  21.359   (   2.876    0.000    0.000)    2.876
  32.392   (   0.013    0.000    0.000)    0.013
  32.392   (   0.013    0.000    0.000)    0.013
  36.126   (   4.482    0.000    0.000)    4.482
======================= Grid point 12 (7/56) =======================
q-point: ( 0.09  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.705   (  16.365   16.365    0.000)   23.144
   0.736   (  16.239   16.239    0.000)   22.965
   1.115   (  23.879   23.879    0.000)   33.770
   3.496   (  -0.089   -0.089    0.000)    0.126
   3.684   (   7.735    7.735    0.000)   10.938
   3.718   (  10.422   10.422    0.000)   14.739
  13.333   (  -0.395   -0.395    0.000)    0.559
  13.552   (  10.940   10.940    0.000)   15.472
  13.819   (  22.797   22.797    0.000)   32.239
  14.154   (   4.175    4.175    0.000)    5.905
  14.225   (   7.193    7.193    0.000)   10.173
  20.723   (   8.171    8.171    0.000)   11.556
  32.390   (  -0.056   -0.056    0.000)    0.079
  32.428   (   1.900    1.900    0.000)    2.687
  35.221   (   6.640    6.640    0.000)    9.390
======================= Grid point 13 (8/56) =======================
q-point: ( 0.18  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.062   (  17.721    8.735    0.000)   19.757
   1.074   (  16.661    8.842    0.000)   18.862
   1.583   (  20.346   10.790    0.000)   23.030
   3.492   (  -0.372    0.518    0.000)    0.638
   3.692   (   0.071   19.586    0.000)   19.586
   4.206   (  29.177    0.824    0.000)   29.189
  13.427   (   8.960  -14.878    0.000)   17.368
  13.849   (  17.615   10.683    0.000)   20.601
  14.294   (  17.443   12.599    0.000)   21.517
  14.381   (  23.131   14.697    0.000)   27.405
  14.427   (  12.413    6.911    0.000)   14.207
  20.894   (   8.688    6.996    0.000)   11.154
  32.389   (  -0.084   -0.209    0.000)    0.225
  32.452   (   0.790    5.002    0.000)    5.064
  35.437   (  14.047    6.004    0.000)   15.276
======================= Grid point 14 (9/56) =======================
q-point: ( 0.27  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.387   (  14.052    3.854    0.000)   14.571
   1.389   (  14.332    4.305    0.000)   14.965
   1.913   (  12.562    3.211    0.000)   12.966
   3.481   (  -0.752    1.394    0.000)    1.584
   3.695   (   0.117   21.613    0.000)   21.614
   4.804   (  28.495    1.591    0.000)   28.539
  13.654   (  12.459  -26.493    0.000)   29.277
  14.217   (  17.837   10.402    0.000)   20.649
  14.604   (  15.450   -0.187    0.000)   15.452
  14.702   (  14.187    6.522    0.000)   15.614
  14.917   (  25.923   27.107    0.000)   37.507
  21.062   (   7.639    1.317    0.000)    7.752
  32.387   (  -0.077   -0.406    0.000)    0.414
  32.465   (   0.529    6.650    0.000)    6.671
  35.751   (  16.160    7.558    0.000)   17.840
======================= Grid point 15 (10/56) =======================
q-point: ( 0.36  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.627   (   9.545    1.400    0.000)    9.648
   1.629   (   9.231    1.547    0.000)    9.360
   2.101   (   6.451   -2.038    0.000)    6.765
   3.464   (  -0.852    2.314    0.000)    2.466
   3.696   (   0.010   23.696    0.000)   23.696
   5.299   (  19.674    1.987    0.000)   19.774
  13.889   (  10.029  -33.875    0.000)   35.329
  14.530   (  12.537   10.137    0.000)   16.122
  14.896   (  12.607   -0.456    0.000)   12.615
  14.968   (  11.482    6.204    0.000)   13.050
  15.368   (  17.976   26.061    0.000)   31.659
  21.191   (   4.965   -4.899    0.000)    6.975
  32.386   (  -0.049   -0.586    0.000)    0.588
  32.474   (   0.383    7.663    0.000)    7.672
  36.050   (  12.700    9.403    0.000)   15.802
======================= Grid point 16 (11/56) =======================
q-point: ( 0.45  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.755   (   3.198    0.115    0.000)    3.200
   1.761   (   3.463    0.681    0.000)    3.529
   2.183   (   1.888   -5.143    0.000)    5.478
   3.450   (  -0.388    2.935    0.000)    2.961
   3.696   (  -0.015   25.029    0.000)   25.029
   5.571   (   6.931    2.140    0.000)    7.254
  14.032   (   3.766  -37.276    0.000)   37.466
  14.704   (   4.431    9.946    0.000)   10.888
  15.079   (   4.893    0.369    0.000)    4.907
  15.134   (   4.431    6.061    0.000)    7.508
  15.617   (   6.359   24.107    0.000)   24.932
  21.258   (   1.680   -8.759    0.000)    8.918
  32.385   (  -0.016   -0.691    0.000)    0.692
  32.480   (   0.145    8.204    0.000)    8.205
  36.232   (   4.801   10.514    0.000)   11.558
======================= Grid point 24 (12/56) =======================
q-point: ( 0.18  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.273   (  10.445   10.445    0.000)   14.772
   1.277   (  11.559   11.559    0.000)   16.346
   1.812   (  10.387   10.387    0.000)   14.689
   3.505   (   0.701    0.701    0.000)    0.991
   4.174   (  16.038   16.038    0.000)   22.682
   4.276   (  16.211   16.211    0.000)   22.926
  13.222   (  -5.377   -5.377    0.000)    7.604
  14.140   (  17.278   17.278    0.000)   24.434
  14.398   (   7.058    7.058    0.000)    9.982
  14.621   (  11.882   11.882    0.000)   16.803
  14.913   (  27.809   27.809    0.000)   39.327
  21.018   (   4.583    4.583    0.000)    6.482
  32.383   (  -0.338   -0.338    0.000)    0.478
  32.550   (   4.249    4.249    0.000)    6.010
  35.648   (  14.675   14.675    0.000)   20.754
======================= Grid point 25 (13/56) =======================
q-point: ( 0.27  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.482   (  10.036    5.387    0.000)   11.391
   1.499   (  10.039    6.097    0.000)   11.745
   1.977   (   5.976    2.730    0.000)    6.570
   3.518   (   0.514    2.150    0.000)    2.211
   4.258   (   2.129   31.966    0.000)   32.037
   4.850   (  29.047    2.821    0.000)   29.184
  13.192   (   1.670  -18.516    0.000)   18.591
  14.502   (  17.538   16.979    0.000)   24.410
  14.599   (  12.355    0.428    0.000)   12.363
  14.884   (  13.591   11.094    0.000)   17.544
  15.423   (  21.775   21.703    0.000)   30.744
  21.076   (   1.466   -0.111    0.000)    1.470
  32.376   (  -0.350   -0.666    0.000)    0.752
  32.620   (   2.789    7.923    0.000)    8.399
  35.992   (  18.202   16.169    0.000)   24.346
======================= Grid point 26 (14/56) =======================
q-point: ( 0.36  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.665   (   7.796    2.285    0.000)    8.124
   1.668   (   6.554    2.195    0.000)    6.912
   2.057   (   2.177   -2.198    0.000)    3.093
   3.525   (   0.213    3.540    0.000)    3.547
   4.294   (   1.411   33.228    0.000)   33.258
   5.353   (  19.967    3.215    0.000)   20.224
  13.255   (   3.668  -27.888    0.000)   28.129
  14.809   (  12.234   16.702    0.000)   20.703
  14.863   (  12.779   -2.991    0.000)   13.125
  15.140   (  11.132   10.419    0.000)   15.247
  15.763   (  11.729   11.828    0.000)   16.658
  21.089   (   0.050   -4.131    0.000)    4.131
  32.369   (  -0.254   -0.963    0.000)    0.996
  32.666   (   1.707   10.380    0.000)   10.520
  36.333   (  14.535   18.253    0.000)   23.333
======================= Grid point 27 (15/56) =======================
q-point: ( 0.45  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.758   (   2.263    0.137    0.000)    2.267
   1.780   (   3.230    1.259    0.000)    3.467
   2.077   (   0.156   -5.007    0.000)    5.009
   3.527   (   0.038    4.402    0.000)    4.402
   4.314   (   0.491   33.837    0.000)   33.841
   5.629   (   7.017    3.416    0.000)    7.805
  13.313   (   1.700  -32.645    0.000)   32.689
  14.978   (   4.285   16.432    0.000)   16.981
  15.063   (   5.839   -2.397    0.000)    6.312
  15.302   (   4.340   10.141    0.000)   11.030
  15.905   (   3.018    2.699    0.000)    4.049
  21.087   (  -0.107   -6.095    0.000)    6.096
  32.366   (  -0.092   -1.137    0.000)    1.141
  32.689   (   0.596   11.595    0.000)   11.611
  36.541   (   5.499   19.639    0.000)   20.394
======================= Grid point 36 (16/56) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.599   (   6.048    6.048    0.000)    8.554
   1.619   (   5.470    5.470    0.000)    7.736
   2.014   (   0.903    0.903    0.000)    1.277
   3.561   (   2.061    2.061    0.000)    2.914
   4.871   (  17.093   17.093    0.000)   24.174
   4.937   (  15.699   15.699    0.000)   22.201
  12.913   (  -9.550   -9.550    0.000)   13.505
  14.645   (   4.573    4.573    0.000)    6.467
  14.860   (  17.447   17.447    0.000)   24.674
  15.127   (  12.552   12.552    0.000)   17.751
  15.751   (  10.156   10.156    0.000)   14.362
  21.064   (  -0.808   -0.808    0.000)    1.143
  32.362   (  -0.700   -0.700    0.000)    0.990
  32.764   (   6.099    6.099    0.000)    8.625
  36.371   (  20.228   20.228    0.000)   28.606
======================= Grid point 37 (17/56) =======================
q-point: ( 0.36  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.711   (   3.589    1.956    0.000)    4.087
   1.718   (   5.583    2.958    0.000)    6.318
   2.014   (  -0.864   -2.121    0.000)    2.290
   3.597   (   1.379    3.377    0.000)    3.648
   4.952   (   2.040   30.499    0.000)   30.567
   5.420   (  20.260    3.250    0.000)   20.519
  12.787   (  -3.594  -18.465    0.000)   18.812
  14.788   (   9.312   -3.649    0.000)   10.001
  15.165   (  12.052   17.597    0.000)   21.328
  15.366   (  10.485   11.403    0.000)   15.490
  15.826   (  -2.122   -5.215    0.000)    5.630
  21.055   (   0.051    0.981    0.000)    0.982
  32.349   (  -0.517   -1.015    0.000)    1.139
  32.863   (   3.755    8.790    0.000)    9.559
  36.752   (  16.359   22.166    0.000)   27.549
======================= Grid point 38 (18/56) =======================
q-point: ( 0.45  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.760   (   1.234    0.088    0.000)    1.237
   1.810   (   2.945    1.568    0.000)    3.336
   1.988   (  -1.271   -3.768    0.000)    3.977
   3.616   (   0.472    4.117    0.000)    4.144
   4.980   (   0.693   30.785    0.000)   30.793
   5.700   (   7.120    3.448    0.000)    7.911
  12.747   (  -0.782  -23.294    0.000)   23.307
  14.975   (   6.999   -6.186    0.000)    9.341
  15.329   (   4.114   17.474    0.000)   17.952
  15.521   (   4.209   10.947    0.000)   11.728
  15.727   (  -5.086  -18.569    0.000)   19.253
  21.061   (   0.292    3.373    0.000)    3.386
  32.342   (  -0.189   -1.200    0.000)    1.215
  32.914   (   1.281   10.193    0.000)   10.274
  36.987   (   6.213   23.396    0.000)   24.207
======================= Grid point 48 (19/56) =======================
q-point: ( 0.36  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.744   (   1.270    1.270    0.000)    1.796
   1.781   (   3.197    3.197    0.000)    4.522
   1.970   (  -2.238   -2.238    0.000)    3.165
   3.655   (   2.270    2.270    0.000)    3.210
   5.467   (  11.855   11.855    0.000)   16.765
   5.486   (  11.047   11.047    0.000)   15.623
  12.505   (  -9.830   -9.830    0.000)   13.902
  14.753   (   1.061    1.061    0.000)    1.501
  15.478   (  12.527   12.527    0.000)   17.716
  15.579   (   9.269    9.269    0.000)   13.109
  15.596   ( -15.890  -15.890    0.000)   22.471
  21.113   (   3.923    3.923    0.000)    5.548
  32.331   (  -0.753   -0.753    0.000)    1.064
  33.009   (   5.580    5.580    0.000)    7.891
  37.170   (  17.875   17.875    0.000)   25.279
======================= Grid point 49 (20/56) =======================
q-point: ( 0.45  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.761   (   0.434    0.035    0.000)    0.435
   1.841   (   2.362    1.421    0.000)    2.757
   1.925   (  -1.952   -2.418    0.000)    3.107
   3.686   (   0.782    2.732    0.000)    2.842
   5.508   (   0.808   20.806    0.000)   20.822
   5.761   (   7.221    2.443    0.000)    7.623
  12.378   (  -2.998  -13.386    0.000)   13.717
  14.839   (   6.224   -6.311    0.000)    8.864
  15.273   ( -13.094  -23.264    0.000)   26.696
  15.645   (   4.012   13.062    0.000)   13.665
  15.720   (   4.040    8.304    0.000)    9.234
  21.179   (   2.021    6.925    0.000)    7.214
  32.320   (  -0.275   -0.890    0.000)    0.932
  33.085   (   1.905    6.576    0.000)    6.846
  37.426   (   6.794   18.780    0.000)   19.971
======================= Grid point 60 (21/56) =======================
q-point: ( 0.45  0.45  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.762   (   0.007    0.007    0.000)    0.011
   1.865   (   1.025    1.025    0.000)    1.450
   1.889   (  -1.291   -1.291    0.000)    1.825
   3.723   (   0.936    0.936    0.000)    1.324
   5.794   (   4.160    4.160    0.000)    5.883
   5.795   (   4.003    4.003    0.000)    5.660
  12.201   (  -4.297   -4.297    0.000)    6.077
  14.764   (  -0.123   -0.123    0.000)    0.173
  14.895   ( -12.607  -12.607    0.000)   17.829
  15.828   (   4.516    4.516    0.000)    6.386
  15.839   (   3.377    3.377    0.000)    4.776
  21.290   (   3.166    3.166    0.000)    4.478
  32.308   (  -0.326   -0.326    0.000)    0.461
  33.175   (   2.252    2.252    0.000)    3.184
  37.696   (   7.135    7.135    0.000)   10.091
======================= Grid point 133 (22/56) =======================
q-point: ( 0.09  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.875   (  12.614   12.614   12.614)   21.848
   0.875   (  12.614   12.614   12.614)   21.848
   1.316   (  17.369   17.369   17.369)   30.084
   3.659   (   5.041    5.041    5.041)    8.731
   3.718   (   7.067    7.067    7.067)   12.241
   3.718   (   7.067    7.067    7.067)   12.241
  13.209   (  -4.150   -4.150   -4.150)    7.188
  13.841   (  11.153   11.153   11.153)   19.318
  13.841   (  11.153   11.153   11.153)   19.318
  14.283   (  10.477   10.477   10.477)   18.147
  14.283   (  10.477   10.477   10.477)   18.147
  20.827   (   7.664    7.664    7.664)   13.274
  32.430   (   1.419    1.419    1.419)    2.457
  32.430   (   1.419    1.419    1.419)    2.457
  35.274   (   6.449    6.449    6.449)   11.170
======================= Grid point 134 (23/56) =======================
q-point: ( 0.18  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.150   (  13.841    7.568    7.568)   17.496
   1.162   (  15.067    7.878    7.878)   18.739
   1.675   (  16.061    8.039    8.039)   19.677
   3.686   (   1.261    9.908    9.908)   14.068
   3.717   (  -0.242   10.962   10.962)   15.505
   4.215   (  29.626    1.047    1.047)   29.663
  13.227   (   5.645  -12.580  -12.580)   18.665
  13.979   (   9.850    4.316    4.316)   11.587
  14.183   (  15.018   -1.593   -1.593)   15.186
  14.528   (  16.064   15.907   15.907)   27.642
  14.704   (  26.388   19.397   19.397)   38.063
  20.971   (   6.298    5.885    5.885)   10.436
  32.437   (   0.423    2.296    2.296)    3.275
  32.467   (   1.065    2.889    2.889)    4.223
  35.492   (  14.535    6.157    6.157)   16.943
======================= Grid point 135 (24/56) =======================
q-point: ( 0.27  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.415   (  12.084    3.473    3.473)   13.044
   1.449   (  12.861    4.484    4.484)   14.340
   1.934   (   9.664    1.774    1.774)    9.984
   3.707   (  -0.742   11.389   11.389)   16.124
   3.717   (   1.662   12.043   12.043)   17.112
   4.821   (  28.751    1.643    1.643)   28.845
  13.405   (  10.976  -19.705  -19.705)   29.951
  14.207   (  11.946   -1.230   -1.230)   12.072
  14.466   (  13.086   -6.443   -6.443)   15.945
  14.887   (  18.487   16.830   16.830)   30.137
  15.252   (  26.236   23.122   23.122)   41.923
  21.076   (   4.010    0.773    0.773)    4.157
  32.446   (   0.482    2.682    2.682)    3.823
  32.480   (   0.374    3.954    3.954)    5.605
  35.821   (  17.085    7.288    7.288)   19.953
======================= Grid point 136 (25/56) =======================
q-point: ( 0.36  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.624   (   8.408    0.850    0.850)    8.494
   1.667   (   8.512    2.417    2.417)    9.172
   2.074   (   4.581   -2.237   -2.237)    5.568
   3.689   (  -0.934   11.724   11.724)   16.607
   3.748   (   1.333   13.984   13.984)   19.821
   5.319   (  19.782    1.954    1.954)   19.974
  13.626   (   9.986  -23.980  -23.980)   35.352
  14.428   (   9.404   -3.225   -3.225)   10.451
  14.698   (   9.462   -9.721   -9.721)   16.689
  15.229   (  14.436   16.763   16.763)   27.756
  15.710   (  18.335   23.696   23.696)   38.199
  21.133   (   1.797   -5.252   -5.252)    7.642
  32.455   (   0.382    3.042    3.042)    4.319
  32.485   (   0.183    4.358    4.358)    6.166
  36.138   (  13.448    8.801    8.801)   18.324
======================= Grid point 137 (26/56) =======================
q-point: ( 0.45  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.743   (   3.101   -0.348   -0.348)    3.140
   1.784   (   2.987    1.418    1.418)    3.598
   2.130   (   1.212   -4.301   -4.301)    6.202
   3.674   (  -0.438   11.909   11.909)   16.847
   3.767   (   0.504   15.115   15.115)   21.381
   5.593   (   6.960    2.073    2.073)    7.553
  13.773   (   3.957  -25.702  -25.702)   36.563
  14.563   (   3.543   -3.678   -3.678)    6.293
  14.831   (   3.468  -11.741  -11.741)   16.963
  15.435   (   5.439   16.667   16.667)   24.190
  15.965   (   6.524   23.533   23.533)   33.914
  21.154   (   0.422   -9.232   -9.232)   13.062
  32.461   (   0.146    3.258    3.258)    4.610
  32.487   (   0.065    4.503    4.503)    6.369
  36.330   (   5.076    9.749    9.749)   14.692
======================= Grid point 145 (27/56) =======================
q-point: ( 0.18  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.332   (   9.356    9.356    5.558)   14.351
   1.346   (  10.140   10.140    5.882)   15.499
   1.846   (   7.778    7.778    2.983)   11.397
   3.726   (   1.536    1.536   20.958)   21.070
   4.198   (  16.323   16.323    2.389)   23.207
   4.283   (  16.331   16.331    0.620)   23.104
  13.042   (  -5.139   -5.139  -14.582)   16.293
  14.008   (   2.959    2.959   -2.757)    5.011
  14.303   (   9.487    9.487  -18.500)   22.852
  14.864   (  17.106   17.106   32.200)   40.275
  15.202   (  27.773   27.773   17.332)   42.931
  21.064   (   2.504    2.504    3.280)    4.826
  32.450   (   0.213    0.213    5.848)    5.856
  32.557   (   4.304    4.304    0.619)    6.118
  35.705   (  14.838   14.838    6.169)   21.872
======================= Grid point 146 (28/56) =======================
q-point: ( 0.27  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.512   (   8.356    5.680    3.189)   10.595
   1.550   (   9.550    5.278    4.166)   11.679
   1.969   (   4.379    1.483   -0.707)    4.678
   3.755   (   1.223    2.507   22.208)   22.382
   4.278   (   2.224   32.048    1.932)   32.183
   4.867   (  29.221    2.798    1.675)   29.402
  13.017   (   2.145  -16.774  -16.128)   23.368
  14.137   (   9.286   -3.494  -12.243)   15.758
  14.496   (   8.862    6.607  -23.861)   26.297
  15.218   (  17.667   15.674   33.014)   40.592
  15.738   (  24.026   23.498   20.688)   39.464
  21.068   (  -1.798   -1.881   -1.407)    2.958
  32.453   (   0.105   -0.243    6.967)    6.972
  32.626   (   2.651    7.965    0.599)    8.416
  36.055   (  18.630   15.796    6.560)   25.291
======================= Grid point 147 (29/56) =======================
q-point: ( 0.36  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.659   (   6.014    2.463    0.868)    6.556
   1.716   (   6.744    2.367    3.111)    7.795
   2.027   (   1.576   -2.210   -2.491)    3.684
   3.773   (   0.536    3.505   22.837)   23.111
   4.317   (   1.579   32.941    2.302)   33.059
   5.373   (  20.018    3.144    1.859)   20.348
  13.094   (   4.583  -25.628  -15.784)   30.446
  14.350   (  11.092   -5.949  -18.937)   22.738
  14.635   (   4.474    2.961  -26.567)   27.103
  15.547   (  14.122   14.379   29.632)   35.836
  16.134   (  14.823   16.676   25.069)   33.560
  21.012   (  -3.183   -5.832   -7.349)    9.907
  32.455   (   0.048   -0.678    7.813)    7.843
  32.667   (   1.515   10.422    0.260)   10.535
  36.405   (  14.955   17.386    7.302)   24.067
======================= Grid point 148 (30/56) =======================
q-point: ( 0.45  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.744   (   2.216    0.429   -0.430)    2.298
   1.813   (   2.588    1.418    2.231)    3.699
   2.042   (   0.137   -4.076   -2.816)    4.956
   3.779   (   0.103    4.182   22.952)   23.330
   4.339   (   0.573   33.269    2.652)   33.380
   5.649   (   7.029    3.310    1.912)    8.001
  13.169   (   2.177  -30.276  -14.668)   33.713
  14.551   (   7.965   -3.051  -23.882)   25.359
  14.662   (  -1.859   -2.632  -26.654)   26.848
  15.753   (   5.567   14.418   26.379)   30.573
  16.329   (   4.708   10.657   28.846)   31.110
  20.961   (  -1.489   -7.906  -11.694)   14.194
  32.456   (   0.014   -0.940    8.268)    8.321
  32.688   (   0.515   11.615   -0.027)   11.626
  36.619   (   5.665   18.507    7.835)   20.880
======================= Grid point 157 (31/56) =======================
q-point: ( 0.27  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.629   (   5.617    5.617    2.945)    8.472
   1.654   (   5.021    5.021    3.206)    7.791
   1.989   (   0.488    0.488   -2.154)    2.262
   3.805   (   2.300    2.300   22.664)   22.897
   4.898   (  17.002   17.002    2.659)   24.191
   4.947   (  15.748   15.748    0.959)   22.292
  12.762   (  -8.539   -8.539  -14.762)   19.072
  14.124   (   2.530    2.530  -23.408)   23.680
  14.637   (   6.082    6.082  -33.235)   34.330
  15.526   (  14.910   14.910   37.724)   43.217
  16.143   (  15.801   15.801   24.271)   32.991
  21.005   (  -3.950   -3.950   -6.060)    8.242
  32.447   (  -0.389   -0.389    7.796)    7.816
  32.769   (   5.971    5.971    0.487)    8.458
  36.428   (  20.055   20.055    5.922)   28.973
======================= Grid point 158 (32/56) =======================
q-point: ( 0.36  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.718   (   3.076    3.181    1.379)    4.635
   1.757   (   5.028    1.702    2.978)    6.087
   1.989   (  -0.518   -1.624   -2.087)    2.694
   3.843   (   1.428    3.308   22.857)   23.140
   4.971   (   1.985   30.351    1.851)   30.472
   5.438   (  20.271    3.148    1.716)   20.586
  12.657   (  -2.550  -17.277  -12.971)   21.754
  14.235   (   8.045   -4.574  -30.550)   31.921
  14.706   (   0.357    2.040  -38.920)   38.975
  15.809   (  12.681   11.795   35.688)   39.669
  16.358   (   5.654    5.764   30.782)   31.823
  20.929   (  -3.121   -1.996  -12.202)   12.752
  32.439   (  -0.333   -0.837    8.339)    8.388
  32.865   (   3.542    8.655    0.124)    9.353
  36.808   (  16.358   21.466    5.699)   27.584
======================= Grid point 159 (33/56) =======================
q-point: ( 0.45  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.758   (   0.971    0.888    0.238)    1.337
   1.839   (   2.590    1.173    2.425)    3.737
   1.972   (  -0.876   -2.824   -1.203)    3.192
   3.862   (   0.457    3.933   22.848)   23.188
   4.998   (   0.683   30.540    1.831)   30.603
   5.717   (   7.114    3.318    1.708)    8.034
  12.634   (  -0.297  -22.068  -11.295)   24.792
  14.415   (   8.331   -7.788  -35.022)   36.832
  14.652   (  -4.866   -2.676  -40.696)   41.073
  16.007   (   5.917   10.719   31.611)   33.899
  16.401   (  -0.093   -1.639   36.958)   36.994
  20.887   (  -1.030    0.704  -16.647)   16.693
  32.434   (  -0.129   -1.105    8.606)    8.677
  32.912   (   1.184   10.040   -0.201)   10.111
  37.043   (   6.233   22.405    5.696)   23.943
======================= Grid point 169 (34/56) =======================
q-point: ( 0.36  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.766   (   1.049    1.049    2.056)    2.535
   1.800   (   3.015    3.015    1.850)    4.648
   1.959   (  -1.443   -1.443   -0.824)    2.201
   3.900   (   2.221    2.221   22.835)   23.050
   5.487   (  11.652   11.652    2.042)   16.604
   5.498   (  11.074   11.074    1.152)   15.704
  12.398   (  -8.717   -8.717  -10.646)   16.289
  14.192   (   0.834    0.834  -36.569)   36.588
  14.673   (  -5.616   -5.616  -45.912)   46.594
  16.026   (   9.630    9.630   35.695)   38.205
  16.396   (  -0.947   -0.947   35.737)   35.762
  20.927   (   1.270    1.270  -17.991)   18.081
  32.423   (  -0.674   -0.674    8.622)    8.674
  33.006   (   5.350    5.350   -0.271)    7.572
  37.214   (  17.493   17.493    4.606)   25.163
======================= Grid point 170 (35/56) =======================
q-point: ( 0.45  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.775   (   0.132    0.709    1.395)    1.570
   1.860   (   2.456    0.961    1.852)    3.223
   1.927   (  -1.594   -1.650    0.343)    2.320
   3.930   (   0.762    2.706   22.823)   22.996
   5.523   (   0.765   20.681    1.491)   20.749
   5.776   (   7.198    2.320    1.491)    7.708
  12.289   (  -2.497  -12.237   -8.804)   15.281
  14.271   (   6.194   -5.949  -40.002)   40.913
  14.508   (  -8.951  -10.707  -45.919)   47.993
  16.185   (   5.528    6.979   33.021)   34.200
  16.352   (  -2.657   -1.858   37.521)   37.661
  20.955   (   1.005    4.808  -21.691)   22.240
  32.414   (  -0.257   -0.875    8.731)    8.778
  33.079   (   1.798    6.333   -0.629)    6.613
  37.466   (   6.679   18.177    4.108)   19.797
======================= Grid point 181 (36/56) =======================
q-point: ( 0.45  0.45  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.784   (   0.202    0.202    2.176)    2.195
   1.880   (   0.972    0.972    1.416)    1.973
   1.901   (  -1.074   -1.074    1.222)    1.949
   3.968   (   0.951    0.951   22.915)   22.955
   5.807   (   4.053    4.053    1.396)    5.899
   5.808   (   4.014    4.014    1.273)    5.818
  12.131   (  -3.782   -3.782   -6.914)    8.741
  14.204   (  -0.003   -0.003  -41.737)   41.737
  14.295   (  -8.458   -8.458  -43.639)   45.249
  16.288   (   3.316    3.316   32.196)   32.535
  16.330   (  -0.885   -0.885   33.115)   33.139
  21.036   (   2.384    2.384  -24.636)   24.865
  32.401   (  -0.332   -0.332    8.746)    8.758
  33.164   (   2.138    2.138   -1.006)    3.187
  37.728   (   6.940    6.940    3.321)   10.361
======================= Grid point 266 (37/56) =======================
q-point: ( 0.18  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.476   (   7.400    7.400    7.400)   12.817
   1.476   (   7.400    7.400    7.400)   12.817
   1.909   (   2.801    2.801    2.801)    4.852
   4.247   (  13.687   13.687   13.687)   23.707
   4.297   (  11.308   11.308   11.308)   19.586
   4.297   (  11.308   11.308   11.308)   19.586
  12.770   ( -10.452  -10.452  -10.452)   18.103
  14.026   (  -2.454   -2.454   -2.454)    4.250
  14.026   (  -2.454   -2.454   -2.454)    4.250
  15.535   (  24.340   24.340   24.340)   42.158
  15.535   (  24.340   24.340   24.340)   42.158
  21.091   (  -1.591   -1.591   -1.591)    2.756
  32.574   (   3.322    3.322    3.322)    5.753
  32.574   (   3.322    3.322    3.322)    5.753
  35.916   (  14.671   14.671   14.671)   25.411
======================= Grid point 267 (38/56) =======================
q-point: ( 0.27  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.601   (   5.656    5.344    5.344)    9.440
   1.643   (   7.887    4.673    4.673)   10.290
   1.953   (   1.588   -0.804   -0.804)    1.953
   4.304   (   0.027   17.056   17.056)   24.120
   4.353   (   4.511   18.072   18.072)   25.953
   4.914   (  29.603    2.796    2.796)   29.865
  12.650   (  -1.517  -17.526  -17.526)   24.833
  13.974   (   2.142   -8.534   -8.534)   12.257
  14.097   (   3.266  -11.466  -11.466)   16.541
  15.839   (  14.186   26.514   26.514)   40.090
  16.144   (  26.820   20.160   20.160)   39.143
  20.993   (  -7.663   -6.457   -6.457)   11.920
  32.598   (   1.190    4.644    4.644)    6.675
  32.660   (   2.992    4.417    4.417)    6.926
  36.268   (  19.115   14.595   14.595)   28.132
======================= Grid point 268 (39/56) =======================
q-point: ( 0.36  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.698   (   3.751    2.795    2.795)    5.449
   1.781   (   5.556    2.999    2.999)    6.990
   1.976   (   0.799   -2.265   -2.265)    3.301
   4.299   (  -0.448   17.380   17.380)   24.583
   4.433   (   3.237   18.561   18.561)   26.448
   5.423   (  20.139    3.034    3.034)   20.591
  12.683   (   3.726  -22.275  -22.275)   31.721
  14.069   (   6.084  -14.083  -14.083)   20.825
  14.144   (   1.308  -17.235  -17.235)   24.409
  16.090   (  10.241   23.244   23.244)   34.430
  16.592   (  17.252   19.223   19.223)   32.198
  20.815   (  -8.919  -11.397  -11.397)   18.422
  32.620   (   0.930    4.823    4.823)    6.884
  32.701   (   1.257    5.710    5.710)    8.172
  36.631   (  15.666   15.065   15.065)   26.445
======================= Grid point 269 (40/56) =======================
q-point: ( 0.45  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.749   (   1.278    0.894    0.894)    1.798
   1.858   (   2.001    2.136    2.136)    3.623
   1.987   (   0.273   -2.335   -2.335)    3.313
   4.290   (  -0.273   17.548   17.548)   24.818
   4.478   (   1.157   18.820   18.820)   26.641
   5.701   (   7.059    3.117    3.117)    8.322
  12.755   (   2.407  -24.003  -24.003)   34.030
  14.155   (   0.070  -21.224  -21.224)   30.015
  14.172   (   3.079  -16.230  -16.230)   23.159
  16.235   (   3.764   20.837   20.837)   29.708
  16.827   (   5.925   18.452   18.452)   26.759
  20.676   (  -3.958  -14.640  -14.640)   21.079
  32.633   (   0.351    4.930    4.930)    6.981
  32.717   (   0.341    6.193    6.193)    8.765
  36.857   (   5.976   15.544   15.544)   22.780
======================= Grid point 278 (41/56) =======================
q-point: ( 0.27  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.709   (   4.583    4.583    4.695)    8.003
   1.732   (   4.386    4.386    4.131)    7.452
   1.942   (  -0.188   -0.188   -2.311)    2.326
   4.383   (   2.571    2.571   32.262)   32.467
   4.972   (  15.999   15.999    1.388)   22.668
   4.974   (  16.829   16.829    4.760)   24.272
  12.379   (  -8.753   -8.753  -20.541)   23.982
  13.783   (  -6.627   -6.627  -12.026)   15.247
  14.008   (  -2.159   -2.159  -27.379)   27.549
  16.244   (  14.613   14.613   32.736)   38.713
  16.561   (  19.886   19.886   16.742)   32.730
  20.811   ( -10.957  -10.957  -13.115)   20.300
  32.634   (   0.730    0.730    9.626)    9.682
  32.783   (   5.623    5.623    0.838)    7.997
  36.620   (  19.233   19.233   13.167)   30.219
======================= Grid point 279 (42/56) =======================
q-point: ( 0.36  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.764   (   1.347    3.584    3.024)    4.879
   1.826   (   4.397    1.309    3.631)    5.851
   1.946   (   0.501   -0.540   -1.895)    2.033
   4.423   (   1.395    3.635   32.223)   32.458
   5.025   (   2.138   29.829    3.535)   30.114
   5.487   (  20.270    3.083    3.062)   20.731
  12.291   (  -0.734  -15.230  -21.390)   26.269
  13.726   (   0.241  -14.655  -20.323)   25.057
  13.955   (  -3.167   -7.735  -34.193)   35.200
  16.471   (   8.333   14.753   29.015)   33.601
  16.888   (  12.010   10.060   21.325)   26.461
  20.588   ( -10.199  -10.335  -20.768)   25.340
  32.644   (   0.330   -0.072   10.754)   10.759
  32.869   (   3.019    8.106    0.386)    8.658
  36.989   (  16.079   19.435   12.301)   28.064
======================= Grid point 280 (43/56) =======================
q-point: ( 0.45  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.777   (   0.212    1.852    1.686)    2.514
   1.894   (   1.997    1.235    2.763)    3.626
   1.953   (   0.049   -0.987   -0.524)    1.119
   4.441   (   0.386    4.442   32.010)   32.319
   5.055   (   0.772   29.673    3.983)   29.949
   5.766   (   7.101    3.160    3.033)    8.343
  12.306   (   1.071  -18.741  -19.944)   27.389
  13.755   (   1.626  -17.255  -30.141)   34.768
  13.893   (  -2.285  -13.711  -33.473)   36.244
  16.584   (   2.814   13.869   24.867)   28.612
  17.043   (   3.654    3.489   25.828)   26.317
  20.437   (  -4.068   -8.128  -26.294)   27.821
  32.648   (   0.094   -0.599   11.342)   11.358
  32.908   (   0.932    9.409   -0.072)    9.456
  37.222   (   6.191   19.680   11.994)   23.864
======================= Grid point 290 (44/56) =======================
q-point: ( 0.36  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.819   (   0.643    0.643    2.968)    3.104
   1.851   (   2.561    2.561    3.075)    4.751
   1.948   (   0.296    0.296   -0.119)    0.435
   4.484   (   2.381    2.381   32.796)   32.969
   5.532   (  11.300   11.300    2.252)   16.138
   5.549   (  11.192   11.192    4.033)   16.334
  12.085   (  -5.644   -5.644  -19.461)   21.034
  13.516   (  -6.201   -6.201  -29.550)   30.824
  13.783   (  -9.251   -9.251  -41.411)   43.428
  16.697   (   7.931    7.931   29.807)   31.847
  17.007   (   2.403    2.403   24.626)   24.860
  20.422   (  -5.891   -5.891  -30.866)   31.970
  32.639   (  -0.332   -0.332   11.578)   11.587
  32.999   (   4.732    4.732   -0.452)    6.707
  37.363   (  16.307   16.307   10.213)   25.222
======================= Grid point 291 (45/56) =======================
q-point: ( 0.45  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.817   (  -0.203    1.854    2.670)    3.257
   1.908   (   2.381    0.389    2.641)    3.577
   1.944   (  -0.848   -0.157    1.317)    1.575
   4.517   (   0.819    3.099   33.045)   33.200
   5.570   (   0.751   20.388    3.231)   20.657
   5.821   (   7.130    2.181    3.051)    8.056
  12.031   (  -0.697   -8.668  -16.251)   18.431
  13.468   (   0.859  -10.976  -38.109)   39.667
  13.609   (  -6.775  -13.877  -42.097)   44.840
  16.799   (   2.417    7.711   26.830)   28.020
  17.021   (  -0.091   -3.884   28.320)   28.586
  20.348   (  -1.684   -1.495  -37.060)   37.128
  32.633   (  -0.159   -0.730   11.897)   11.920
  33.061   (   1.511    5.664   -1.043)    5.954
  37.599   (   6.310   16.462    9.191)   19.882
======================= Grid point 302 (46/56) =======================
q-point: ( 0.45  0.45  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.844   (   0.696    0.696    3.552)    3.686
   1.919   (   0.836    0.836    2.409)    2.683
   1.936   (  -0.685   -0.685    2.190)    2.394
   4.561   (   1.148    1.148   33.632)   33.671
   5.849   (   4.114    4.114    2.923)    6.512
   5.853   (   3.793    3.793    3.174)    6.232
  11.933   (  -1.875   -1.875  -12.281)   12.564
  13.327   (  -2.641   -2.641  -43.185)   43.346
  13.380   (  -7.540   -7.540  -45.246)   46.486
  16.900   (   2.398    2.398   27.114)   27.325
  16.949   (  -2.255   -2.255   27.051)   27.239
  20.342   (   0.220    0.220  -42.577)   42.578
  32.622   (  -0.308   -0.308   12.044)   12.052
  33.136   (   1.833    1.833   -1.694)    3.096
  37.838   (   6.376    6.376    7.747)   11.888
======================= Grid point 399 (47/56) =======================
q-point: ( 0.27  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.808   (   3.894    3.894    3.894)    6.745
   1.808   (   3.894    3.894    3.894)    6.745
   1.906   (  -0.919   -0.919   -0.919)    1.592
   4.998   (  11.398   11.398   11.398)   19.742
   4.998   (  11.398   11.398   11.398)   19.742
   5.105   (  13.332   13.332   13.332)   23.092
  12.032   ( -12.491  -12.491  -12.491)   21.635
  13.572   ( -12.862  -12.862  -12.862)   22.278
  13.572   ( -12.862  -12.862  -12.862)   22.278
  16.822   (  16.898   16.898   16.898)   29.268
  16.822   (  16.898   16.898   16.898)   29.268
  20.492   ( -17.805  -17.805  -17.805)   30.840
  32.801   (   3.753    3.753    3.753)    6.500
  32.801   (   3.753    3.753    3.753)    6.500
  36.941   (  17.850   17.850   17.850)   30.917
======================= Grid point 400 (48/56) =======================
q-point: ( 0.36  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.831   (  -0.266    3.398    3.398)    4.813
   1.894   (   3.546    2.497    2.497)    5.005
   1.920   (   1.887   -0.227   -0.227)    1.914
   5.012   (   0.377   17.285   17.285)   24.448
   5.154   (   4.038   16.340   16.340)   23.458
   5.556   (  20.070    3.593    3.593)   20.703
  11.894   (  -1.661  -15.970  -15.970)   22.646
  13.361   ( -10.497  -20.867  -20.867)   31.322
  13.367   (  -5.354  -19.870  -19.870)   28.605
  16.955   (   3.984   18.655   18.655)   26.681
  17.216   (  14.721   11.589   11.589)   22.031
  20.141   ( -15.685  -22.198  -22.198)   35.093
  32.818   (   0.731    4.729    4.729)    6.727
  32.896   (   3.059    3.561    3.561)    5.893
  37.289   (  15.334   16.691   16.691)   28.148
======================= Grid point 401 (49/56) =======================
q-point: ( 0.45  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.823   (  -0.300    2.686    2.686)    3.810
   1.944   (   1.301    1.979    1.979)    3.087
   1.952   (   0.916    0.441    0.441)    1.108
   5.015   (   0.029   17.763   17.763)   25.120
   5.209   (   1.412   16.831   16.831)   23.845
   5.833   (   7.079    3.423    3.423)    8.576
  11.913   (   1.795  -16.858  -16.858)   23.908
  13.186   (  -5.448  -25.829  -25.829)   36.931
  13.317   (  -0.654  -24.818  -24.818)   35.104
  16.998   (   0.787   16.356   16.356)   23.144
  17.414   (   4.910   10.579   10.579)   15.746
  19.910   (  -6.224  -23.851  -23.851)   34.300
  32.829   (   0.276    4.469    4.469)    6.326
  32.933   (   0.818    4.084    4.084)    5.834
  37.513   (   5.986   16.137   16.137)   23.593
======================= Grid point 411 (50/56) =======================
q-point: ( 0.36  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.879   (   0.550    0.550    2.792)    2.898
   1.917   (   2.023    2.023    3.274)    4.348
   1.955   (   1.922    1.922    0.750)    2.820
   5.138   (   3.000    3.000   30.345)   30.640
   5.585   (  12.020   12.020    2.849)   17.236
   5.646   (  10.778   10.778    5.512)   16.208
  11.691   (  -4.170   -4.170  -17.374)   18.347
  12.987   ( -14.115  -14.115  -23.872)   31.118
  13.041   ( -14.546  -14.546  -31.549)   37.663
  17.208   (   6.263    6.263   20.852)   22.655
  17.386   (   6.352    6.352   13.060)   15.852
  19.750   ( -15.436  -15.436  -33.827)   40.259
  32.850   (   0.272    0.272    8.732)    8.741
  32.989   (   3.916    3.916   -0.485)    5.559
  37.616   (  14.538   14.538   14.348)   25.071
======================= Grid point 412 (51/56) =======================
q-point: ( 0.45  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.879   (  -0.096    2.643    3.316)    4.242
   1.959   (   1.532    0.548    2.321)    2.834
   1.977   (   0.264    0.976    1.838)    2.098
   5.177   (   0.977    4.439   30.750)   31.084
   5.653   (   1.119   19.891    4.989)   20.538
   5.896   (   7.005    2.581    4.179)    8.555
  11.684   (   1.499   -5.534  -16.435)   17.407
  12.750   (  -7.878  -16.922  -33.050)   37.957
  12.856   (  -3.932  -20.883  -32.284)   38.650
  17.250   (  -0.134    8.869   17.931)   20.005
  17.481   (   2.606   -2.573   17.109)   17.497
  19.543   (  -5.061  -12.043  -40.432)   42.490
  32.851   (   0.001   -0.232    9.020)    9.023
  33.040   (   1.178    4.577   -0.987)    4.828
  37.829   (   5.728   14.092   13.103)   20.077
======================= Grid point 423 (52/56) =======================
q-point: ( 0.45  0.45  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.920   (   1.185    1.185    3.790)    4.144
   1.971   (   0.672    0.672    2.609)    2.777
   1.984   (  -0.273   -0.273    2.355)    2.387
   5.239   (   1.622    1.622   32.008)   32.090
   5.926   (   4.474    4.474    4.593)    7.819
   5.936   (   3.544    3.544    4.972)    7.060
  11.660   (   1.296    1.296  -13.986)   14.106
  12.516   (  -7.648   -7.648  -36.514)   38.083
  12.523   (  -8.473   -8.473  -38.842)   40.648
  17.360   (   1.665    1.665   18.456)   18.605
  17.403   (  -2.287   -2.287   18.066)   18.353
  19.404   (  -2.801   -2.801  -48.151)   48.313
  32.845   (  -0.218   -0.218    9.388)    9.393
  33.100   (   1.431    1.431   -1.845)    2.739
  38.036   (   5.572    5.572   11.456)   13.904
======================= Grid point 532 (53/56) =======================
q-point: ( 0.36  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.924   (   1.592    1.592    1.592)    2.758
   1.976   (   1.852    1.852    1.852)    3.208
   1.976   (   1.852    1.852    1.852)    3.208
   5.640   (   9.607    9.607    9.607)   16.640
   5.640   (   9.607    9.607    9.607)   16.640
   5.775   (   8.596    8.596    8.596)   14.889
  11.451   (  -5.746   -5.746   -5.746)    9.952
  12.531   ( -20.015  -20.015  -20.015)   34.668
  12.531   ( -20.015  -20.015  -20.015)   34.668
  17.538   (   7.234    7.234    7.234)   12.529
  17.538   (   7.234    7.234    7.234)   12.529
  19.148   ( -23.936  -23.936  -23.936)   41.459
  32.980   (   1.984    1.984    1.984)    3.437
  32.980   (   1.984    1.984    1.984)    3.437
  37.902   (  12.871   12.871   12.871)   22.293
======================= Grid point 533 (54/56) =======================
q-point: ( 0.45  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.943   (   0.428    2.932    2.932)    4.169
   1.998   (   0.660    1.495    1.495)    2.215
   2.011   (   0.893    1.286    1.286)    2.026
   5.673   (   0.884   14.257   14.257)   20.182
   5.790   (   2.445   11.519   11.519)   16.472
   5.980   (   6.712    3.855    3.855)    8.647
  11.437   (   2.229   -6.916   -6.916)   10.031
  12.177   ( -11.515  -22.769  -22.769)   34.197
  12.283   (  -6.014  -24.972  -24.972)   35.824
  17.528   (  -1.655    9.833    9.833)   14.004
  17.698   (   3.906    4.007    4.007)    6.882
  18.832   (  -7.419  -27.177  -27.177)   39.143
  32.975   (  -0.136    2.707    2.707)    3.831
  33.033   (   1.174    1.060    1.060)    1.904
  38.092   (   5.151   11.922   11.922)   17.630
======================= Grid point 544 (55/56) =======================
q-point: ( 0.45  0.45  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.989   (   1.390    1.390    2.936)    3.534
   2.018   (   0.529    0.529    1.937)    2.076
   2.025   (   0.192    0.192    1.626)    1.649
   5.795   (   2.302    2.302   22.149)   22.387
   6.025   (   5.207    5.207    4.804)    8.792
   6.043   (   3.464    3.464    5.250)    7.180
  11.410   (   2.743    2.743   -9.736)   10.480
  11.854   ( -10.372  -10.372  -26.729)   30.490
  11.880   ( -11.534  -11.534  -26.151)   30.821
  17.648   (   1.093    1.093   10.323)   10.438
  17.681   (  -1.414   -1.414    9.623)    9.828
  18.509   (  -6.193   -6.193  -37.983)   38.980
  32.987   (  -0.119   -0.119    4.767)    4.769
  33.066   (   1.061    1.061   -1.410)    2.059
  38.269   (   4.811    4.811   10.687)   12.669
======================= Grid point 665 (56/56) =======================
q-point: ( 0.45  0.45  0.45)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.98e-04 1.15e-04 1.15e-04 1.15e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.033   (   1.309    1.309    1.309)    2.267
   2.046   (   0.527    0.527    0.527)    0.914
   2.046   (   0.527    0.527    0.527)    0.914
   6.100   (   4.688    4.688    4.688)    8.119
   6.100   (   4.688    4.688    4.688)    8.119
   6.128   (   2.983    2.983    2.983)    5.166
  11.306   (  -0.267   -0.267   -0.267)    0.462
  11.480   ( -11.250  -11.250  -11.250)   19.485
  11.480   ( -11.250  -11.250  -11.250)   19.485
  17.784   (   1.608    1.608    1.608)    2.785
  17.784   (   1.608    1.608    1.608)    2.785
  17.984   ( -11.480  -11.480  -11.480)   19.883
  33.046   (   0.385    0.385    0.385)    0.667
  33.046   (   0.385    0.385    0.385)    0.667
  38.428   (   4.366    4.366    4.366)    7.563
=================== 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   1025.490   1025.490   1025.490      0.000     -0.000      0.000 3/19965
   20.0    757.442    757.442    757.442      0.000     -0.000      0.000 3/19965
   30.0    438.720    438.720    438.720     -0.000     -0.000      0.000 3/19965
   40.0    236.988    236.988    236.988     -0.000     -0.000      0.000 3/19965
   50.0    144.730    144.730    144.730     -0.000     -0.000      0.000 3/19965
   60.0     98.955     98.955     98.955     -0.000     -0.000      0.000 3/19965
   70.0     72.149     72.149     72.149     -0.000     -0.000      0.000 3/19965
   80.0     54.448     54.448     54.448     -0.000     -0.000      0.000 3/19965
   90.0     42.176     42.176     42.176      0.000     -0.000      0.000 3/19965
  100.0     33.571     33.571     33.571      0.000     -0.000      0.000 3/19965
  110.0     27.489     27.489     27.489      0.000     -0.000      0.000 3/19965
  120.0     23.122     23.122     23.122      0.000     -0.000      0.000 3/19965
  130.0     19.911     19.911     19.911      0.000     -0.000      0.000 3/19965
  140.0     17.482     17.482     17.482      0.000     -0.000      0.000 3/19965
  150.0     15.590     15.590     15.590      0.000     -0.000      0.000 3/19965
  160.0     14.073     14.073     14.073      0.000     -0.000      0.000 3/19965
  170.0     12.822     12.822     12.822      0.000     -0.000      0.000 3/19965
  180.0     11.767     11.767     11.767      0.000     -0.000      0.000 3/19965
  190.0     10.858     10.858     10.858      0.000     -0.000      0.000 3/19965
  200.0     10.061     10.061     10.061      0.000     -0.000      0.000 3/19965
  210.0      9.354      9.354      9.354      0.000     -0.000      0.000 3/19965
  220.0      8.720      8.720      8.720      0.000     -0.000      0.000 3/19965
  230.0      8.147      8.147      8.147      0.000     -0.000      0.000 3/19965
  240.0      7.627      7.627      7.627      0.000     -0.000      0.000 3/19965
  250.0      7.153      7.153      7.153      0.000     -0.000      0.000 3/19965
  260.0      6.721      6.721      6.721      0.000     -0.000      0.000 3/19965
  270.0      6.325      6.325      6.325      0.000     -0.000      0.000 3/19965
  280.0      5.964      5.964      5.964      0.000     -0.000      0.000 3/19965
  290.0      5.632      5.632      5.632      0.000     -0.000      0.000 3/19965
  300.0      5.329      5.329      5.329      0.000     -0.000      0.000 3/19965
  310.0      5.050      5.050      5.050      0.000     -0.000      0.000 3/19965
  320.0      4.794      4.794      4.794      0.000     -0.000      0.000 3/19965
  330.0      4.559      4.559      4.559      0.000     -0.000      0.000 3/19965
  340.0      4.342      4.342      4.342      0.000     -0.000      0.000 3/19965
  350.0      4.143      4.143      4.143      0.000     -0.000      0.000 3/19965
  360.0      3.959      3.959      3.959      0.000     -0.000      0.000 3/19965
  370.0      3.788      3.788      3.788      0.000     -0.000      0.000 3/19965
  380.0      3.631      3.631      3.631      0.000     -0.000      0.000 3/19965
  390.0      3.485      3.485      3.485      0.000     -0.000      0.000 3/19965
  400.0      3.349      3.349      3.349      0.000     -0.000      0.000 3/19965
  410.0      3.223      3.223      3.223      0.000     -0.000      0.000 3/19965
  420.0      3.105      3.105      3.105      0.000     -0.000      0.000 3/19965
  430.0      2.996      2.996      2.996      0.000     -0.000      0.000 3/19965
  440.0      2.893      2.893      2.893      0.000     -0.000      0.000 3/19965
  450.0      2.797      2.797      2.797      0.000     -0.000      0.000 3/19965
  460.0      2.707      2.707      2.707      0.000     -0.000      0.000 3/19965
  470.0      2.623      2.623      2.623      0.000     -0.000      0.000 3/19965
  480.0      2.543      2.543      2.543      0.000     -0.000      0.000 3/19965
  490.0      2.469      2.469      2.469      0.000     -0.000      0.000 3/19965
  500.0      2.398      2.398      2.398      0.000     -0.000      0.000 3/19965
  510.0      2.331      2.331      2.331      0.000     -0.000      0.000 3/19965
  520.0      2.268      2.268      2.268      0.000     -0.000      0.000 3/19965
  530.0      2.209      2.209      2.209      0.000     -0.000      0.000 3/19965
  540.0      2.152      2.152      2.152      0.000     -0.000      0.000 3/19965
  550.0      2.098      2.098      2.098      0.000     -0.000      0.000 3/19965
  560.0      2.047      2.047      2.047      0.000     -0.000      0.000 3/19965
  570.0      1.999      1.999      1.999      0.000     -0.000      0.000 3/19965
  580.0      1.953      1.953      1.953      0.000     -0.000      0.000 3/19965
  590.0      1.908      1.908      1.908      0.000     -0.000      0.000 3/19965
  600.0      1.866      1.866      1.866      0.000     -0.000      0.000 3/19965
  610.0      1.826      1.826      1.826      0.000     -0.000      0.000 3/19965
  620.0      1.788      1.788      1.788      0.000     -0.000      0.000 3/19965
  630.0      1.751      1.751      1.751      0.000     -0.000      0.000 3/19965
  640.0      1.715      1.715      1.715      0.000     -0.000      0.000 3/19965
  650.0      1.682      1.682      1.682      0.000     -0.000      0.000 3/19965
  660.0      1.649      1.649      1.649      0.000     -0.000      0.000 3/19965
  670.0      1.618      1.618      1.618      0.000     -0.000      0.000 3/19965
  680.0      1.588      1.588      1.588      0.000     -0.000      0.000 3/19965
  690.0      1.559      1.559      1.559      0.000     -0.000      0.000 3/19965
  700.0      1.531      1.531      1.531      0.000     -0.000      0.000 3/19965
  710.0      1.504      1.504      1.504      0.000     -0.000      0.000 3/19965
  720.0      1.479      1.479      1.479      0.000     -0.000      0.000 3/19965
  730.0      1.454      1.454      1.454      0.000     -0.000      0.000 3/19965
  740.0      1.430      1.430      1.430      0.000     -0.000      0.000 3/19965
  750.0      1.406      1.406      1.406      0.000     -0.000      0.000 3/19965
  760.0      1.384      1.384      1.384      0.000     -0.000      0.000 3/19965
  770.0      1.362      1.362      1.362      0.000     -0.000      0.000 3/19965
  780.0      1.341      1.341      1.341      0.000     -0.000      0.000 3/19965
  790.0      1.321      1.321      1.321      0.000     -0.000      0.000 3/19965
  800.0      1.301      1.301      1.301      0.000     -0.000      0.000 3/19965
  810.0      1.282      1.282      1.282      0.000     -0.000      0.000 3/19965
  820.0      1.263      1.263      1.263      0.000     -0.000      0.000 3/19965
  830.0      1.245      1.245      1.245      0.000     -0.000      0.000 3/19965
  840.0      1.228      1.228      1.228      0.000     -0.000      0.000 3/19965
  850.0      1.211      1.211      1.211      0.000     -0.000      0.000 3/19965
  860.0      1.195      1.195      1.195      0.000     -0.000      0.000 3/19965
  870.0      1.179      1.179      1.179      0.000     -0.000      0.000 3/19965
  880.0      1.163      1.163      1.163      0.000     -0.000      0.000 3/19965
  890.0      1.148      1.148      1.148      0.000     -0.000      0.000 3/19965
  900.0      1.133      1.133      1.133      0.000     -0.000      0.000 3/19965
  910.0      1.119      1.119      1.119      0.000     -0.000      0.000 3/19965
  920.0      1.105      1.105      1.105      0.000     -0.000      0.000 3/19965
  930.0      1.091      1.091      1.091      0.000     -0.000      0.000 3/19965
  940.0      1.078      1.078      1.078      0.000     -0.000      0.000 3/19965
  950.0      1.065      1.065      1.065      0.000     -0.000      0.000 3/19965
  960.0      1.052      1.052      1.052      0.000     -0.000      0.000 3/19965
  970.0      1.040      1.040      1.040      0.000     -0.000      0.000 3/19965
  980.0      1.028      1.028      1.028      0.000     -0.000      0.000 3/19965
  990.0      1.016      1.016      1.016      0.000     -0.000      0.000 3/19965
 1000.0      1.005      1.005      1.005      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-08 21:14:49]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

