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

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

-------------------------[time 2026-01-08 18:05: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    5.205557870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.205557870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    5.205557870000000
Atomic positions (fractional):
   *1 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
   *2 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411
   *3 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453
    4 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453
    5 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    5.205557870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.205557870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    5.205557870000000
Atomic positions (fractional):
   *1 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 1
   *2 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411 > 2
   *3 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453 > 3
    4 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 4
    5 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   15.616673609999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   15.616673609999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   15.616673609999999
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 Cd  0.16666666666667  0.16666666666667  0.16666666666667 112.411 > 2
   29 Cd  0.50000000000000  0.16666666666667  0.16666666666667 112.411 > 2
   30 Cd  0.83333333333333  0.16666666666667  0.16666666666667 112.411 > 2
   31 Cd  0.16666666666667  0.50000000000000  0.16666666666667 112.411 > 2
   32 Cd  0.50000000000000  0.50000000000000  0.16666666666667 112.411 > 2
   33 Cd  0.83333333333333  0.50000000000000  0.16666666666667 112.411 > 2
   34 Cd  0.16666666666667  0.83333333333333  0.16666666666667 112.411 > 2
   35 Cd  0.50000000000000  0.83333333333333  0.16666666666667 112.411 > 2
   36 Cd  0.83333333333333  0.83333333333333  0.16666666666667 112.411 > 2
   37 Cd  0.16666666666667  0.16666666666667  0.50000000000000 112.411 > 2
   38 Cd  0.50000000000000  0.16666666666667  0.50000000000000 112.411 > 2
   39 Cd  0.83333333333333  0.16666666666667  0.50000000000000 112.411 > 2
   40 Cd  0.16666666666667  0.50000000000000  0.50000000000000 112.411 > 2
   41 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411 > 2
   42 Cd  0.83333333333333  0.50000000000000  0.50000000000000 112.411 > 2
   43 Cd  0.16666666666667  0.83333333333333  0.50000000000000 112.411 > 2
   44 Cd  0.50000000000000  0.83333333333333  0.50000000000000 112.411 > 2
   45 Cd  0.83333333333333  0.83333333333333  0.50000000000000 112.411 > 2
   46 Cd  0.16666666666667  0.16666666666667  0.83333333333333 112.411 > 2
   47 Cd  0.50000000000000  0.16666666666667  0.83333333333333 112.411 > 2
   48 Cd  0.83333333333333  0.16666666666667  0.83333333333333 112.411 > 2
   49 Cd  0.16666666666667  0.50000000000000  0.83333333333333 112.411 > 2
   50 Cd  0.50000000000000  0.50000000000000  0.83333333333333 112.411 > 2
   51 Cd  0.83333333333333  0.50000000000000  0.83333333333333 112.411 > 2
   52 Cd  0.16666666666667  0.83333333333333  0.83333333333333 112.411 > 2
   53 Cd  0.50000000000000  0.83333333333333  0.83333333333333 112.411 > 2
   54 Cd  0.83333333333333  0.83333333333333  0.83333333333333 112.411 > 2
  *55 Cl  0.00000000000000  0.16666666666667  0.16666666666667  35.453 > 3
   56 Cl  0.33333333333333  0.16666666666667  0.16666666666667  35.453 > 3
   57 Cl  0.66666666666667  0.16666666666667  0.16666666666667  35.453 > 3
   58 Cl  0.00000000000000  0.50000000000000  0.16666666666667  35.453 > 3
   59 Cl  0.33333333333333  0.50000000000000  0.16666666666667  35.453 > 3
   60 Cl  0.66666666666667  0.50000000000000  0.16666666666667  35.453 > 3
   61 Cl  0.00000000000000  0.83333333333333  0.16666666666667  35.453 > 3
   62 Cl  0.33333333333333  0.83333333333333  0.16666666666667  35.453 > 3
   63 Cl  0.66666666666667  0.83333333333333  0.16666666666667  35.453 > 3
   64 Cl  0.00000000000000  0.16666666666667  0.50000000000000  35.453 > 3
   65 Cl  0.33333333333333  0.16666666666667  0.50000000000000  35.453 > 3
   66 Cl  0.66666666666667  0.16666666666667  0.50000000000000  35.453 > 3
   67 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453 > 3
   68 Cl  0.33333333333333  0.50000000000000  0.50000000000000  35.453 > 3
   69 Cl  0.66666666666667  0.50000000000000  0.50000000000000  35.453 > 3
   70 Cl  0.00000000000000  0.83333333333333  0.50000000000000  35.453 > 3
   71 Cl  0.33333333333333  0.83333333333333  0.50000000000000  35.453 > 3
   72 Cl  0.66666666666667  0.83333333333333  0.50000000000000  35.453 > 3
   73 Cl  0.00000000000000  0.16666666666667  0.83333333333333  35.453 > 3
   74 Cl  0.33333333333333  0.16666666666667  0.83333333333333  35.453 > 3
   75 Cl  0.66666666666667  0.16666666666667  0.83333333333333  35.453 > 3
   76 Cl  0.00000000000000  0.50000000000000  0.83333333333333  35.453 > 3
   77 Cl  0.33333333333333  0.50000000000000  0.83333333333333  35.453 > 3
   78 Cl  0.66666666666667  0.50000000000000  0.83333333333333  35.453 > 3
   79 Cl  0.00000000000000  0.83333333333333  0.83333333333333  35.453 > 3
   80 Cl  0.33333333333333  0.83333333333333  0.83333333333333  35.453 > 3
   81 Cl  0.66666666666667  0.83333333333333  0.83333333333333  35.453 > 3
   82 Cl  0.16666666666667  0.16666666666667  0.00000000000000  35.453 > 4
   83 Cl  0.50000000000000  0.16666666666667  0.00000000000000  35.453 > 4
   84 Cl  0.83333333333333  0.16666666666667  0.00000000000000  35.453 > 4
   85 Cl  0.16666666666667  0.50000000000000  0.00000000000000  35.453 > 4
   86 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 4
   87 Cl  0.83333333333333  0.50000000000000  0.00000000000000  35.453 > 4
   88 Cl  0.16666666666667  0.83333333333333  0.00000000000000  35.453 > 4
   89 Cl  0.50000000000000  0.83333333333333  0.00000000000000  35.453 > 4
   90 Cl  0.83333333333333  0.83333333333333  0.00000000000000  35.453 > 4
   91 Cl  0.16666666666667  0.16666666666667  0.33333333333333  35.453 > 4
   92 Cl  0.50000000000000  0.16666666666667  0.33333333333333  35.453 > 4
   93 Cl  0.83333333333333  0.16666666666667  0.33333333333333  35.453 > 4
   94 Cl  0.16666666666667  0.50000000000000  0.33333333333333  35.453 > 4
   95 Cl  0.50000000000000  0.50000000000000  0.33333333333333  35.453 > 4
   96 Cl  0.83333333333333  0.50000000000000  0.33333333333333  35.453 > 4
   97 Cl  0.16666666666667  0.83333333333333  0.33333333333333  35.453 > 4
   98 Cl  0.50000000000000  0.83333333333333  0.33333333333333  35.453 > 4
   99 Cl  0.83333333333333  0.83333333333333  0.33333333333333  35.453 > 4
  100 Cl  0.16666666666667  0.16666666666667  0.66666666666667  35.453 > 4
  101 Cl  0.50000000000000  0.16666666666667  0.66666666666667  35.453 > 4
  102 Cl  0.83333333333333  0.16666666666667  0.66666666666667  35.453 > 4
  103 Cl  0.16666666666667  0.50000000000000  0.66666666666667  35.453 > 4
  104 Cl  0.50000000000000  0.50000000000000  0.66666666666667  35.453 > 4
  105 Cl  0.83333333333333  0.50000000000000  0.66666666666667  35.453 > 4
  106 Cl  0.16666666666667  0.83333333333333  0.66666666666667  35.453 > 4
  107 Cl  0.50000000000000  0.83333333333333  0.66666666666667  35.453 > 4
  108 Cl  0.83333333333333  0.83333333333333  0.66666666666667  35.453 > 4
  109 Cl  0.16666666666667  0.00000000000000  0.16666666666667  35.453 > 5
  110 Cl  0.50000000000000  0.00000000000000  0.16666666666667  35.453 > 5
  111 Cl  0.83333333333333  0.00000000000000  0.16666666666667  35.453 > 5
  112 Cl  0.16666666666667  0.33333333333333  0.16666666666667  35.453 > 5
  113 Cl  0.50000000000000  0.33333333333333  0.16666666666667  35.453 > 5
  114 Cl  0.83333333333333  0.33333333333333  0.16666666666667  35.453 > 5
  115 Cl  0.16666666666667  0.66666666666667  0.16666666666667  35.453 > 5
  116 Cl  0.50000000000000  0.66666666666667  0.16666666666667  35.453 > 5
  117 Cl  0.83333333333333  0.66666666666667  0.16666666666667  35.453 > 5
  118 Cl  0.16666666666667  0.00000000000000  0.50000000000000  35.453 > 5
  119 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453 > 5
  120 Cl  0.83333333333333  0.00000000000000  0.50000000000000  35.453 > 5
  121 Cl  0.16666666666667  0.33333333333333  0.50000000000000  35.453 > 5
  122 Cl  0.50000000000000  0.33333333333333  0.50000000000000  35.453 > 5
  123 Cl  0.83333333333333  0.33333333333333  0.50000000000000  35.453 > 5
  124 Cl  0.16666666666667  0.66666666666667  0.50000000000000  35.453 > 5
  125 Cl  0.50000000000000  0.66666666666667  0.50000000000000  35.453 > 5
  126 Cl  0.83333333333333  0.66666666666667  0.50000000000000  35.453 > 5
  127 Cl  0.16666666666667  0.00000000000000  0.83333333333333  35.453 > 5
  128 Cl  0.50000000000000  0.00000000000000  0.83333333333333  35.453 > 5
  129 Cl  0.83333333333333  0.00000000000000  0.83333333333333  35.453 > 5
  130 Cl  0.16666666666667  0.33333333333333  0.83333333333333  35.453 > 5
  131 Cl  0.50000000000000  0.33333333333333  0.83333333333333  35.453 > 5
  132 Cl  0.83333333333333  0.33333333333333  0.83333333333333  35.453 > 5
  133 Cl  0.16666666666667  0.66666666666667  0.83333333333333  35.453 > 5
  134 Cl  0.50000000000000  0.66666666666667  0.83333333333333  35.453 > 5
  135 Cl  0.83333333333333  0.66666666666667  0.83333333333333  35.453 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.6012901    0.0000000    0.0000000
            0.0000000    3.6012901    0.0000000
            0.0000000    0.0000000    3.6012901
-------------------------- Born effective charges --------------------------
    1 Cs    1.3495474    0.0000000    0.0000000
            0.0000000    1.3495474    0.0000000
            0.0000000    0.0000000    1.3495474
    2 Cd    2.7294263    0.0000000    0.0000000
            0.0000000    2.7294263    0.0000000
            0.0000000    0.0000000    2.7294263
    3 Cl   -2.2701857    0.0000000    0.0000000
            0.0000000   -0.9043940    0.0000000
            0.0000000    0.0000000   -0.9043940
    4 Cl   -0.9043940    0.0000000    0.0000000
            0.0000000   -0.9043940    0.0000000
            0.0000000    0.0000000   -2.2701857
    5 Cl   -0.9043940    0.0000000    0.0000000
            0.0000000   -2.2701857    0.0000000
            0.0000000    0.0000000   -0.9043940
----------------------------------------------------------------------------
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.065
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.070
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 
Permutation basis: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) 
Force constants are written into "force_constants.hdf5".

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

Summary of calculation was written in "phonopy.yaml".
-------------------------[time 2026-01-08 18:05:38]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 18:05: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    5.205557870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.205557870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    5.205557870000000
Atomic positions (fractional):
    1 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
    2 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411
    3 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453
    4 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453
    5 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.616673609999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   15.616673609999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   15.616673609999999
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 Cd  0.16666666666667  0.16666666666667  0.16666666666667 112.411 > 28
   29 Cd  0.50000000000000  0.16666666666667  0.16666666666667 112.411 > 28
   30 Cd  0.83333333333333  0.16666666666667  0.16666666666667 112.411 > 28
   31 Cd  0.16666666666667  0.50000000000000  0.16666666666667 112.411 > 28
   32 Cd  0.50000000000000  0.50000000000000  0.16666666666667 112.411 > 28
   33 Cd  0.83333333333333  0.50000000000000  0.16666666666667 112.411 > 28
   34 Cd  0.16666666666667  0.83333333333333  0.16666666666667 112.411 > 28
   35 Cd  0.50000000000000  0.83333333333333  0.16666666666667 112.411 > 28
   36 Cd  0.83333333333333  0.83333333333333  0.16666666666667 112.411 > 28
   37 Cd  0.16666666666667  0.16666666666667  0.50000000000000 112.411 > 28
   38 Cd  0.50000000000000  0.16666666666667  0.50000000000000 112.411 > 28
   39 Cd  0.83333333333333  0.16666666666667  0.50000000000000 112.411 > 28
   40 Cd  0.16666666666667  0.50000000000000  0.50000000000000 112.411 > 28
   41 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411 > 28
   42 Cd  0.83333333333333  0.50000000000000  0.50000000000000 112.411 > 28
   43 Cd  0.16666666666667  0.83333333333333  0.50000000000000 112.411 > 28
   44 Cd  0.50000000000000  0.83333333333333  0.50000000000000 112.411 > 28
   45 Cd  0.83333333333333  0.83333333333333  0.50000000000000 112.411 > 28
   46 Cd  0.16666666666667  0.16666666666667  0.83333333333333 112.411 > 28
   47 Cd  0.50000000000000  0.16666666666667  0.83333333333333 112.411 > 28
   48 Cd  0.83333333333333  0.16666666666667  0.83333333333333 112.411 > 28
   49 Cd  0.16666666666667  0.50000000000000  0.83333333333333 112.411 > 28
   50 Cd  0.50000000000000  0.50000000000000  0.83333333333333 112.411 > 28
   51 Cd  0.83333333333333  0.50000000000000  0.83333333333333 112.411 > 28
   52 Cd  0.16666666666667  0.83333333333333  0.83333333333333 112.411 > 28
   53 Cd  0.50000000000000  0.83333333333333  0.83333333333333 112.411 > 28
   54 Cd  0.83333333333333  0.83333333333333  0.83333333333333 112.411 > 28
   55 Cl  0.00000000000000  0.16666666666667  0.16666666666667  35.453 > 55
   56 Cl  0.33333333333333  0.16666666666667  0.16666666666667  35.453 > 55
   57 Cl  0.66666666666667  0.16666666666667  0.16666666666667  35.453 > 55
   58 Cl  0.00000000000000  0.50000000000000  0.16666666666667  35.453 > 55
   59 Cl  0.33333333333333  0.50000000000000  0.16666666666667  35.453 > 55
   60 Cl  0.66666666666667  0.50000000000000  0.16666666666667  35.453 > 55
   61 Cl  0.00000000000000  0.83333333333333  0.16666666666667  35.453 > 55
   62 Cl  0.33333333333333  0.83333333333333  0.16666666666667  35.453 > 55
   63 Cl  0.66666666666667  0.83333333333333  0.16666666666667  35.453 > 55
   64 Cl  0.00000000000000  0.16666666666667  0.50000000000000  35.453 > 55
   65 Cl  0.33333333333333  0.16666666666667  0.50000000000000  35.453 > 55
   66 Cl  0.66666666666667  0.16666666666667  0.50000000000000  35.453 > 55
   67 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453 > 55
   68 Cl  0.33333333333333  0.50000000000000  0.50000000000000  35.453 > 55
   69 Cl  0.66666666666667  0.50000000000000  0.50000000000000  35.453 > 55
   70 Cl  0.00000000000000  0.83333333333333  0.50000000000000  35.453 > 55
   71 Cl  0.33333333333333  0.83333333333333  0.50000000000000  35.453 > 55
   72 Cl  0.66666666666667  0.83333333333333  0.50000000000000  35.453 > 55
   73 Cl  0.00000000000000  0.16666666666667  0.83333333333333  35.453 > 55
   74 Cl  0.33333333333333  0.16666666666667  0.83333333333333  35.453 > 55
   75 Cl  0.66666666666667  0.16666666666667  0.83333333333333  35.453 > 55
   76 Cl  0.00000000000000  0.50000000000000  0.83333333333333  35.453 > 55
   77 Cl  0.33333333333333  0.50000000000000  0.83333333333333  35.453 > 55
   78 Cl  0.66666666666667  0.50000000000000  0.83333333333333  35.453 > 55
   79 Cl  0.00000000000000  0.83333333333333  0.83333333333333  35.453 > 55
   80 Cl  0.33333333333333  0.83333333333333  0.83333333333333  35.453 > 55
   81 Cl  0.66666666666667  0.83333333333333  0.83333333333333  35.453 > 55
   82 Cl  0.16666666666667  0.16666666666667  0.00000000000000  35.453 > 82
   83 Cl  0.50000000000000  0.16666666666667  0.00000000000000  35.453 > 82
   84 Cl  0.83333333333333  0.16666666666667  0.00000000000000  35.453 > 82
   85 Cl  0.16666666666667  0.50000000000000  0.00000000000000  35.453 > 82
   86 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 82
   87 Cl  0.83333333333333  0.50000000000000  0.00000000000000  35.453 > 82
   88 Cl  0.16666666666667  0.83333333333333  0.00000000000000  35.453 > 82
   89 Cl  0.50000000000000  0.83333333333333  0.00000000000000  35.453 > 82
   90 Cl  0.83333333333333  0.83333333333333  0.00000000000000  35.453 > 82
   91 Cl  0.16666666666667  0.16666666666667  0.33333333333333  35.453 > 82
   92 Cl  0.50000000000000  0.16666666666667  0.33333333333333  35.453 > 82
   93 Cl  0.83333333333333  0.16666666666667  0.33333333333333  35.453 > 82
   94 Cl  0.16666666666667  0.50000000000000  0.33333333333333  35.453 > 82
   95 Cl  0.50000000000000  0.50000000000000  0.33333333333333  35.453 > 82
   96 Cl  0.83333333333333  0.50000000000000  0.33333333333333  35.453 > 82
   97 Cl  0.16666666666667  0.83333333333333  0.33333333333333  35.453 > 82
   98 Cl  0.50000000000000  0.83333333333333  0.33333333333333  35.453 > 82
   99 Cl  0.83333333333333  0.83333333333333  0.33333333333333  35.453 > 82
  100 Cl  0.16666666666667  0.16666666666667  0.66666666666667  35.453 > 82
  101 Cl  0.50000000000000  0.16666666666667  0.66666666666667  35.453 > 82
  102 Cl  0.83333333333333  0.16666666666667  0.66666666666667  35.453 > 82
  103 Cl  0.16666666666667  0.50000000000000  0.66666666666667  35.453 > 82
  104 Cl  0.50000000000000  0.50000000000000  0.66666666666667  35.453 > 82
  105 Cl  0.83333333333333  0.50000000000000  0.66666666666667  35.453 > 82
  106 Cl  0.16666666666667  0.83333333333333  0.66666666666667  35.453 > 82
  107 Cl  0.50000000000000  0.83333333333333  0.66666666666667  35.453 > 82
  108 Cl  0.83333333333333  0.83333333333333  0.66666666666667  35.453 > 82
  109 Cl  0.16666666666667  0.00000000000000  0.16666666666667  35.453 > 109
  110 Cl  0.50000000000000  0.00000000000000  0.16666666666667  35.453 > 109
  111 Cl  0.83333333333333  0.00000000000000  0.16666666666667  35.453 > 109
  112 Cl  0.16666666666667  0.33333333333333  0.16666666666667  35.453 > 109
  113 Cl  0.50000000000000  0.33333333333333  0.16666666666667  35.453 > 109
  114 Cl  0.83333333333333  0.33333333333333  0.16666666666667  35.453 > 109
  115 Cl  0.16666666666667  0.66666666666667  0.16666666666667  35.453 > 109
  116 Cl  0.50000000000000  0.66666666666667  0.16666666666667  35.453 > 109
  117 Cl  0.83333333333333  0.66666666666667  0.16666666666667  35.453 > 109
  118 Cl  0.16666666666667  0.00000000000000  0.50000000000000  35.453 > 109
  119 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453 > 109
  120 Cl  0.83333333333333  0.00000000000000  0.50000000000000  35.453 > 109
  121 Cl  0.16666666666667  0.33333333333333  0.50000000000000  35.453 > 109
  122 Cl  0.50000000000000  0.33333333333333  0.50000000000000  35.453 > 109
  123 Cl  0.83333333333333  0.33333333333333  0.50000000000000  35.453 > 109
  124 Cl  0.16666666666667  0.66666666666667  0.50000000000000  35.453 > 109
  125 Cl  0.50000000000000  0.66666666666667  0.50000000000000  35.453 > 109
  126 Cl  0.83333333333333  0.66666666666667  0.50000000000000  35.453 > 109
  127 Cl  0.16666666666667  0.00000000000000  0.83333333333333  35.453 > 109
  128 Cl  0.50000000000000  0.00000000000000  0.83333333333333  35.453 > 109
  129 Cl  0.83333333333333  0.00000000000000  0.83333333333333  35.453 > 109
  130 Cl  0.16666666666667  0.33333333333333  0.83333333333333  35.453 > 109
  131 Cl  0.50000000000000  0.33333333333333  0.83333333333333  35.453 > 109
  132 Cl  0.83333333333333  0.33333333333333  0.83333333333333  35.453 > 109
  133 Cl  0.16666666666667  0.66666666666667  0.83333333333333  35.453 > 109
  134 Cl  0.50000000000000  0.66666666666667  0.83333333333333  35.453 > 109
  135 Cl  0.83333333333333  0.66666666666667  0.83333333333333  35.453 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.6012901    0.0000000    0.0000000
            0.0000000    3.6012901    0.0000000
            0.0000000    0.0000000    3.6012901
-------------------------- Born effective charges --------------------------
    1 Cs    1.3495474    0.0000000    0.0000000
            0.0000000    1.3495474    0.0000000
            0.0000000    0.0000000    1.3495474
    2 Cd    2.7294263    0.0000000    0.0000000
            0.0000000    2.7294263    0.0000000
            0.0000000    0.0000000    2.7294263
    3 Cl   -2.2701857    0.0000000    0.0000000
            0.0000000   -0.9043940    0.0000000
            0.0000000    0.0000000   -0.9043940
    4 Cl   -0.9043940    0.0000000    0.0000000
            0.0000000   -0.9043940    0.0000000
            0.0000000    0.0000000   -2.2701857
    5 Cl   -0.9043940    0.0000000    0.0000000
            0.0000000   -2.2701857    0.0000000
            0.0000000    0.0000000   -0.9043940
----------------------------------------------------------------------------
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.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  0.0000]
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000000 (yyy) 0.00000000 (yyy) 0.00000000 (yyy)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (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 18:05:42]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 18:05: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    5.205557870000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    5.205557870000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    5.205557870000000
Atomic positions (fractional):
    1 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
    2 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411
    3 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453
    4 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453
    5 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   15.616673609999999    0.000000000000000    0.000000000000000
  b    0.000000000000000   15.616673609999999    0.000000000000000
  c    0.000000000000000    0.000000000000000   15.616673609999999
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 Cd  0.16666666666667  0.16666666666667  0.16666666666667 112.411 > 28
   29 Cd  0.50000000000000  0.16666666666667  0.16666666666667 112.411 > 28
   30 Cd  0.83333333333333  0.16666666666667  0.16666666666667 112.411 > 28
   31 Cd  0.16666666666667  0.50000000000000  0.16666666666667 112.411 > 28
   32 Cd  0.50000000000000  0.50000000000000  0.16666666666667 112.411 > 28
   33 Cd  0.83333333333333  0.50000000000000  0.16666666666667 112.411 > 28
   34 Cd  0.16666666666667  0.83333333333333  0.16666666666667 112.411 > 28
   35 Cd  0.50000000000000  0.83333333333333  0.16666666666667 112.411 > 28
   36 Cd  0.83333333333333  0.83333333333333  0.16666666666667 112.411 > 28
   37 Cd  0.16666666666667  0.16666666666667  0.50000000000000 112.411 > 28
   38 Cd  0.50000000000000  0.16666666666667  0.50000000000000 112.411 > 28
   39 Cd  0.83333333333333  0.16666666666667  0.50000000000000 112.411 > 28
   40 Cd  0.16666666666667  0.50000000000000  0.50000000000000 112.411 > 28
   41 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411 > 28
   42 Cd  0.83333333333333  0.50000000000000  0.50000000000000 112.411 > 28
   43 Cd  0.16666666666667  0.83333333333333  0.50000000000000 112.411 > 28
   44 Cd  0.50000000000000  0.83333333333333  0.50000000000000 112.411 > 28
   45 Cd  0.83333333333333  0.83333333333333  0.50000000000000 112.411 > 28
   46 Cd  0.16666666666667  0.16666666666667  0.83333333333333 112.411 > 28
   47 Cd  0.50000000000000  0.16666666666667  0.83333333333333 112.411 > 28
   48 Cd  0.83333333333333  0.16666666666667  0.83333333333333 112.411 > 28
   49 Cd  0.16666666666667  0.50000000000000  0.83333333333333 112.411 > 28
   50 Cd  0.50000000000000  0.50000000000000  0.83333333333333 112.411 > 28
   51 Cd  0.83333333333333  0.50000000000000  0.83333333333333 112.411 > 28
   52 Cd  0.16666666666667  0.83333333333333  0.83333333333333 112.411 > 28
   53 Cd  0.50000000000000  0.83333333333333  0.83333333333333 112.411 > 28
   54 Cd  0.83333333333333  0.83333333333333  0.83333333333333 112.411 > 28
   55 Cl  0.00000000000000  0.16666666666667  0.16666666666667  35.453 > 55
   56 Cl  0.33333333333333  0.16666666666667  0.16666666666667  35.453 > 55
   57 Cl  0.66666666666667  0.16666666666667  0.16666666666667  35.453 > 55
   58 Cl  0.00000000000000  0.50000000000000  0.16666666666667  35.453 > 55
   59 Cl  0.33333333333333  0.50000000000000  0.16666666666667  35.453 > 55
   60 Cl  0.66666666666667  0.50000000000000  0.16666666666667  35.453 > 55
   61 Cl  0.00000000000000  0.83333333333333  0.16666666666667  35.453 > 55
   62 Cl  0.33333333333333  0.83333333333333  0.16666666666667  35.453 > 55
   63 Cl  0.66666666666667  0.83333333333333  0.16666666666667  35.453 > 55
   64 Cl  0.00000000000000  0.16666666666667  0.50000000000000  35.453 > 55
   65 Cl  0.33333333333333  0.16666666666667  0.50000000000000  35.453 > 55
   66 Cl  0.66666666666667  0.16666666666667  0.50000000000000  35.453 > 55
   67 Cl  0.00000000000000  0.50000000000000  0.50000000000000  35.453 > 55
   68 Cl  0.33333333333333  0.50000000000000  0.50000000000000  35.453 > 55
   69 Cl  0.66666666666667  0.50000000000000  0.50000000000000  35.453 > 55
   70 Cl  0.00000000000000  0.83333333333333  0.50000000000000  35.453 > 55
   71 Cl  0.33333333333333  0.83333333333333  0.50000000000000  35.453 > 55
   72 Cl  0.66666666666667  0.83333333333333  0.50000000000000  35.453 > 55
   73 Cl  0.00000000000000  0.16666666666667  0.83333333333333  35.453 > 55
   74 Cl  0.33333333333333  0.16666666666667  0.83333333333333  35.453 > 55
   75 Cl  0.66666666666667  0.16666666666667  0.83333333333333  35.453 > 55
   76 Cl  0.00000000000000  0.50000000000000  0.83333333333333  35.453 > 55
   77 Cl  0.33333333333333  0.50000000000000  0.83333333333333  35.453 > 55
   78 Cl  0.66666666666667  0.50000000000000  0.83333333333333  35.453 > 55
   79 Cl  0.00000000000000  0.83333333333333  0.83333333333333  35.453 > 55
   80 Cl  0.33333333333333  0.83333333333333  0.83333333333333  35.453 > 55
   81 Cl  0.66666666666667  0.83333333333333  0.83333333333333  35.453 > 55
   82 Cl  0.16666666666667  0.16666666666667  0.00000000000000  35.453 > 82
   83 Cl  0.50000000000000  0.16666666666667  0.00000000000000  35.453 > 82
   84 Cl  0.83333333333333  0.16666666666667  0.00000000000000  35.453 > 82
   85 Cl  0.16666666666667  0.50000000000000  0.00000000000000  35.453 > 82
   86 Cl  0.50000000000000  0.50000000000000  0.00000000000000  35.453 > 82
   87 Cl  0.83333333333333  0.50000000000000  0.00000000000000  35.453 > 82
   88 Cl  0.16666666666667  0.83333333333333  0.00000000000000  35.453 > 82
   89 Cl  0.50000000000000  0.83333333333333  0.00000000000000  35.453 > 82
   90 Cl  0.83333333333333  0.83333333333333  0.00000000000000  35.453 > 82
   91 Cl  0.16666666666667  0.16666666666667  0.33333333333333  35.453 > 82
   92 Cl  0.50000000000000  0.16666666666667  0.33333333333333  35.453 > 82
   93 Cl  0.83333333333333  0.16666666666667  0.33333333333333  35.453 > 82
   94 Cl  0.16666666666667  0.50000000000000  0.33333333333333  35.453 > 82
   95 Cl  0.50000000000000  0.50000000000000  0.33333333333333  35.453 > 82
   96 Cl  0.83333333333333  0.50000000000000  0.33333333333333  35.453 > 82
   97 Cl  0.16666666666667  0.83333333333333  0.33333333333333  35.453 > 82
   98 Cl  0.50000000000000  0.83333333333333  0.33333333333333  35.453 > 82
   99 Cl  0.83333333333333  0.83333333333333  0.33333333333333  35.453 > 82
  100 Cl  0.16666666666667  0.16666666666667  0.66666666666667  35.453 > 82
  101 Cl  0.50000000000000  0.16666666666667  0.66666666666667  35.453 > 82
  102 Cl  0.83333333333333  0.16666666666667  0.66666666666667  35.453 > 82
  103 Cl  0.16666666666667  0.50000000000000  0.66666666666667  35.453 > 82
  104 Cl  0.50000000000000  0.50000000000000  0.66666666666667  35.453 > 82
  105 Cl  0.83333333333333  0.50000000000000  0.66666666666667  35.453 > 82
  106 Cl  0.16666666666667  0.83333333333333  0.66666666666667  35.453 > 82
  107 Cl  0.50000000000000  0.83333333333333  0.66666666666667  35.453 > 82
  108 Cl  0.83333333333333  0.83333333333333  0.66666666666667  35.453 > 82
  109 Cl  0.16666666666667  0.00000000000000  0.16666666666667  35.453 > 109
  110 Cl  0.50000000000000  0.00000000000000  0.16666666666667  35.453 > 109
  111 Cl  0.83333333333333  0.00000000000000  0.16666666666667  35.453 > 109
  112 Cl  0.16666666666667  0.33333333333333  0.16666666666667  35.453 > 109
  113 Cl  0.50000000000000  0.33333333333333  0.16666666666667  35.453 > 109
  114 Cl  0.83333333333333  0.33333333333333  0.16666666666667  35.453 > 109
  115 Cl  0.16666666666667  0.66666666666667  0.16666666666667  35.453 > 109
  116 Cl  0.50000000000000  0.66666666666667  0.16666666666667  35.453 > 109
  117 Cl  0.83333333333333  0.66666666666667  0.16666666666667  35.453 > 109
  118 Cl  0.16666666666667  0.00000000000000  0.50000000000000  35.453 > 109
  119 Cl  0.50000000000000  0.00000000000000  0.50000000000000  35.453 > 109
  120 Cl  0.83333333333333  0.00000000000000  0.50000000000000  35.453 > 109
  121 Cl  0.16666666666667  0.33333333333333  0.50000000000000  35.453 > 109
  122 Cl  0.50000000000000  0.33333333333333  0.50000000000000  35.453 > 109
  123 Cl  0.83333333333333  0.33333333333333  0.50000000000000  35.453 > 109
  124 Cl  0.16666666666667  0.66666666666667  0.50000000000000  35.453 > 109
  125 Cl  0.50000000000000  0.66666666666667  0.50000000000000  35.453 > 109
  126 Cl  0.83333333333333  0.66666666666667  0.50000000000000  35.453 > 109
  127 Cl  0.16666666666667  0.00000000000000  0.83333333333333  35.453 > 109
  128 Cl  0.50000000000000  0.00000000000000  0.83333333333333  35.453 > 109
  129 Cl  0.83333333333333  0.00000000000000  0.83333333333333  35.453 > 109
  130 Cl  0.16666666666667  0.33333333333333  0.83333333333333  35.453 > 109
  131 Cl  0.50000000000000  0.33333333333333  0.83333333333333  35.453 > 109
  132 Cl  0.83333333333333  0.33333333333333  0.83333333333333  35.453 > 109
  133 Cl  0.16666666666667  0.66666666666667  0.83333333333333  35.453 > 109
  134 Cl  0.50000000000000  0.66666666666667  0.83333333333333  35.453 > 109
  135 Cl  0.83333333333333  0.66666666666667  0.83333333333333  35.453 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.6012901    0.0000000    0.0000000
            0.0000000    3.6012901    0.0000000
            0.0000000    0.0000000    3.6012901
-------------------------- Born effective charges --------------------------
    1 Cs    1.3495474    0.0000000    0.0000000
            0.0000000    1.3495474    0.0000000
            0.0000000    0.0000000    1.3495474
    2 Cd    2.7294263    0.0000000    0.0000000
            0.0000000    2.7294263    0.0000000
            0.0000000    0.0000000    2.7294263
    3 Cl   -2.2701857    0.0000000    0.0000000
            0.0000000   -0.9043940    0.0000000
            0.0000000    0.0000000   -0.9043940
    4 Cl   -0.9043940    0.0000000    0.0000000
            0.0000000   -0.9043940    0.0000000
            0.0000000    0.0000000   -2.2701857
    5 Cl   -0.9043940    0.0000000    0.0000000
            0.0000000   -2.2701857    0.0000000
            0.0000000    0.0000000   -0.9043940
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (yyy) 0.00000000 (yyy) 0.00000000 (yyy)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 10 10 10 ]
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.80, Number of G-points: 305, Lambda: 0.16
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.72e-04 5.83e-04 5.83e-04 5.83e-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
   1.401   (   0.000    0.000    0.000)    0.000
   1.401   (   0.000    0.000    0.000)    0.000
   1.401   (   0.000    0.000    0.000)    0.000
   2.202   (   0.000    0.000    0.000)    0.000
   2.202   (   0.000    0.000    0.000)    0.000
   2.202   (   0.000    0.000    0.000)    0.000
   2.226   (   0.000    0.000    0.000)    0.000
   2.226   (   0.000    0.000    0.000)    0.000
   2.226   (   0.000    0.000    0.000)    0.000
   7.104   (   0.000    0.000    0.000)    0.000
   7.104   (   0.000    0.000    0.000)    0.000
   7.104   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/56) =======================
q-point: ( 0.10  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 105
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.324   (  16.282    0.000    0.000)   16.282
   0.324   (  16.282    0.000    0.000)   16.282
   0.689   (  29.390    0.000    0.000)   29.390
   1.404   (   0.212    0.000    0.000)    0.212
   1.404   (   0.212    0.000    0.000)    0.212
   1.623   (  23.169    0.000    0.000)   23.169
   2.214   (   0.852    0.000    0.000)    0.852
   2.214   (   0.852    0.000    0.000)    0.852
   2.222   (   1.938    0.000    0.000)    1.938
   2.242   (   1.912    0.000    0.000)    1.912
   2.242   (   1.912    0.000    0.000)    1.912
   3.575   (   0.317    0.000    0.000)    0.317
   7.104   (   0.005    0.000    0.000)    0.005
   7.104   (   0.005    0.000    0.000)    0.005
   8.088   (  -5.544    0.000    0.000)    5.544
======================= Grid point 2 (3/56) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.616   (  13.916    0.000    0.000)   13.916
   0.616   (  13.916    0.000    0.000)   13.916
   1.118   (  17.144    0.000    0.000)   17.144
   1.406   (  -0.105    0.000    0.000)    0.105
   1.406   (  -0.105    0.000    0.000)    0.105
   2.218   (  34.665    0.000    0.000)   34.665
   2.229   (   0.742    0.000    0.000)    0.742
   2.229   (   0.742    0.000    0.000)    0.742
   2.271   (   3.063    0.000    0.000)    3.063
   2.300   (   3.863    0.000    0.000)    3.863
   2.300   (   3.863    0.000    0.000)    3.863
   3.593   (   2.062    0.000    0.000)    2.062
   7.104   (   0.020    0.000    0.000)    0.020
   7.104   (   0.020    0.000    0.000)    0.020
   7.934   ( -10.161    0.000    0.000)   10.161
======================= Grid point 3 (4/56) =======================
q-point: ( 0.30  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 105
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.854   (  10.810    0.000    0.000)   10.810
   0.854   (  10.810    0.000    0.000)   10.810
   1.395   (  -1.066    0.000    0.000)    1.066
   1.395   (  -1.066    0.000    0.000)    1.066
   1.397   (  12.329    0.000    0.000)   12.329
   2.243   (   0.685    0.000    0.000)    0.685
   2.243   (   0.685    0.000    0.000)    0.685
   2.331   (   2.978    0.000    0.000)    2.978
   2.379   (   4.043    0.000    0.000)    4.043
   2.379   (   4.043    0.000    0.000)    4.043
   2.808   (  23.668    0.000    0.000)   23.668
   3.689   (   9.156    0.000    0.000)    9.156
   7.104   (   0.035    0.000    0.000)    0.035
   7.104   (   0.035    0.000    0.000)    0.035
   7.713   ( -12.279    0.000    0.000)   12.279
======================= Grid point 4 (5/56) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.027   (   6.905    0.000    0.000)    6.905
   1.027   (   6.905    0.000    0.000)    6.905
   1.365   (  -1.907    0.000    0.000)    1.907
   1.365   (  -1.907    0.000    0.000)    1.907
   1.587   (   7.074    0.000    0.000)    7.074
   2.254   (   0.422    0.000    0.000)    0.422
   2.254   (   0.422    0.000    0.000)    0.422
   2.379   (   1.801    0.000    0.000)    1.801
   2.445   (   2.540    0.000    0.000)    2.540
   2.445   (   2.540    0.000    0.000)    2.540
   3.080   (   6.270    0.000    0.000)    6.270
   3.934   (  13.471    0.000    0.000)   13.471
   7.105   (   0.029    0.000    0.000)    0.029
   7.105   (   0.029    0.000    0.000)    0.029
   7.495   (  -9.273    0.000    0.000)    9.273
======================= Grid point 5 (6/56) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 63
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.100   (  -0.000    0.000    0.000)    0.000
   1.100   (  -0.000    0.000    0.000)    0.000
   1.340   (   0.000    0.000    0.000)    0.000
   1.340   (   0.000    0.000    0.000)    0.000
   1.657   (   0.000    0.000    0.000)    0.000
   2.258   (   0.000    0.000    0.000)    0.000
   2.258   (   0.000    0.000    0.000)    0.000
   2.397   (   0.000    0.000    0.000)    0.000
   2.470   (   0.000    0.000    0.000)    0.000
   2.470   (   0.000    0.000    0.000)    0.000
   3.129   (   0.000    0.000    0.000)    0.000
   4.082   (   0.000    0.000    0.000)    0.000
   7.105   (   0.000    0.000    0.000)    0.000
   7.105   (   0.000    0.000    0.000)    0.000
   7.399   (   0.000    0.000    0.000)    0.000
======================= Grid point 12 (7/56) =======================
q-point: ( 0.10  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 180
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.456   (  11.279   11.279    0.000)   15.951
   0.534   (  11.012   11.012    0.000)   15.574
   0.897   (  17.566   17.566    0.000)   24.842
   1.408   (   0.371    0.371    0.000)    0.524
   1.604   (  11.836   11.836    0.000)   16.739
   1.641   (  11.345   11.345    0.000)   16.045
   2.150   (  -2.708   -2.708    0.000)    3.830
   2.233   (   1.521    1.521    0.000)    2.151
   2.250   (   1.207    1.207    0.000)    1.708
   2.262   (   2.485    2.485    0.000)    3.514
   2.304   (   4.963    4.963    0.000)    7.019
   3.577   (   0.286    0.286    0.000)    0.405
   6.931   (  -9.077   -9.077    0.000)   12.837
   7.104   (   0.014    0.014    0.000)    0.019
   8.180   (   1.710    1.710    0.000)    2.419
======================= Grid point 13 (8/56) =======================
q-point: ( 0.20  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.688   (  11.974    6.782    0.000)   13.761
   0.734   (  10.438    7.481    0.000)   12.842
   1.182   (  11.667    5.390    0.000)   12.852
   1.415   (   0.225    0.857    0.000)    0.886
   1.645   (   1.883   23.405    0.000)   23.481
   2.077   (   6.367   -8.708    0.000)   10.787
   2.195   (  14.674   -1.351    0.000)   14.736
   2.251   (   0.718    2.293    0.000)    2.403
   2.304   (   3.835    0.349    0.000)    3.851
   2.330   (   8.546    7.278    0.000)   11.225
   2.445   (  11.942    7.097    0.000)   13.891
   3.593   (   1.861    0.093    0.000)    1.863
   6.799   (  -5.689  -22.453    0.000)   23.163
   7.104   (   0.033    0.034    0.000)    0.048
   8.124   (  -6.228   11.228    0.000)   12.840
======================= Grid point 14 (9/56) =======================
q-point: ( 0.30  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.898   (   9.603    4.138    0.000)   10.457
   0.926   (   9.131    4.808    0.000)   10.320
   1.356   (   6.542    0.220    0.000)    6.545
   1.412   (  -0.531    1.709    0.000)    1.789
   1.694   (   3.244   23.336    0.000)   23.560
   2.102   (   0.641  -13.899    0.000)   13.913
   2.265   (   0.653    2.232    0.000)    2.326
   2.301   (   2.014    6.323    0.000)    6.636
   2.382   (   4.028    0.315    0.000)    4.041
   2.499   (   5.877    4.495    0.000)    7.399
   2.875   (  22.395    6.253    0.000)   23.252
   3.683   (   8.795   -0.535    0.000)    8.811
   6.691   (  -5.646  -31.188    0.000)   31.695
   7.105   (   0.047    0.059    0.000)    0.075
   7.977   (  -8.278   16.238    0.000)   18.227
======================= Grid point 15 (10/56) =======================
q-point: ( 0.40  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.049   (   5.904    2.111    0.000)    6.270
   1.068   (   5.345    3.108    0.000)    6.183
   1.395   (  -1.173    2.941    0.000)    3.167
   1.438   (   2.298    2.182    0.000)    3.169
   1.759   (   3.010   19.269    0.000)   19.502
   2.115   (   0.607  -17.222    0.000)   17.232
   2.275   (   0.404    2.185    0.000)    2.222
   2.331   (   1.174    9.079    0.000)    9.154
   2.448   (   2.532    0.293    0.000)    2.549
   2.584   (   2.967    5.055    0.000)    5.861
   3.137   (   6.059    5.545    0.000)    8.213
   3.926   (  13.516   -0.853    0.000)   13.543
   6.588   (  -4.613  -36.894    0.000)   37.181
   7.106   (   0.037    0.078    0.000)    0.086
   7.836   (  -5.763   18.955    0.000)   19.812
======================= Grid point 16 (11/56) =======================
q-point: (-0.50  0.10  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 180
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.110   (  -0.000    0.980    0.000)    0.980
   1.121   (   0.000    2.129    0.000)    2.129
   1.380   (   0.000    3.840    0.000)    3.840
   1.458   (   0.000    4.738    0.000)    4.738
   1.791   (   0.000   16.479    0.000)   16.479
   2.122   (   0.000  -18.434    0.000)   18.434
   2.279   (   0.000    2.168    0.000)    2.168
   2.343   (   0.000   10.124    0.000)   10.124
   2.473   (   0.000    0.286    0.000)    0.286
   2.612   (   0.000    5.098    0.000)    5.098
   3.184   (   0.000    5.411    0.000)    5.411
   4.075   (   0.000   -0.746    0.000)    0.746
   6.539   (   0.000  -38.957    0.000)   38.957
   7.106   (   0.000    0.085    0.000)    0.085
   7.777   (   0.000   19.775    0.000)   19.775
======================= Grid point 24 (12/56) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 186
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.842   (   8.482    8.482    0.000)   11.996
   0.863   (   6.832    6.832    0.000)    9.662
   1.287   (   5.138    5.138    0.000)    7.266
   1.436   (   1.211    1.211    0.000)    1.713
   1.878   (  -4.577   -4.577    0.000)    6.473
   2.143   (   4.189    4.189    0.000)    5.925
   2.236   (  15.721   15.721    0.000)   22.232
   2.311   (   2.395    2.395    0.000)    3.387
   2.314   (   2.034    2.034    0.000)    2.877
   2.550   (  14.510   14.510    0.000)   20.520
   2.569   (   9.895    9.895    0.000)   13.994
   3.603   (   1.254    1.254    0.000)    1.773
   6.400   ( -18.494  -18.494    0.000)   26.155
   7.105   (   0.065    0.065    0.000)    0.092
   8.239   (   0.900    0.900    0.000)    1.273
======================= Grid point 25 (13/56) =======================
q-point: ( 0.30  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.993   (   6.868    5.248    0.000)    8.643
   1.005   (   7.506    3.612    0.000)    8.330
   1.373   (   3.789    1.223    0.000)    3.981
   1.457   (   0.924    2.768    0.000)    2.918
   1.799   (  -3.005  -12.819    0.000)   13.166
   2.139   (  -0.982   13.229    0.000)   13.265
   2.324   (   0.583    3.773    0.000)    3.818
   2.390   (   3.992    0.526    0.000)    4.027
   2.480   (   8.515    6.466    0.000)   10.692
   2.653   (   3.907   20.025    0.000)   20.402
   3.015   (  22.135    7.814    0.000)   23.474
   3.673   (   7.488   -0.321    0.000)    7.495
   6.083   ( -14.640  -30.751    0.000)   34.058
   7.107   (   0.078    0.103    0.000)    0.129
   8.193   (  -4.627    5.968    0.000)    7.551
======================= Grid point 26 (14/56) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 312
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.097   (   3.816    2.638    0.000)    4.639
   1.128   (   4.727    2.819    0.000)    5.504
   1.436   (   2.851   -1.329    0.000)    3.145
   1.469   (   0.366    4.497    0.000)    4.511
   1.757   (  -1.560  -18.264    0.000)   18.331
   2.122   (  -0.635   13.282    0.000)   13.297
   2.334   (   0.366    3.722    0.000)    3.740
   2.455   (   2.513    0.492    0.000)    2.561
   2.604   (   4.673    4.222    0.000)    6.298
   2.707   (   1.661   23.584    0.000)   23.642
   3.272   (   5.666    7.976    0.000)    9.783
   3.904   (  13.639   -1.366    0.000)   13.708
   5.845   (  -9.601  -39.186    0.000)   40.345
   7.108   (   0.055    0.132    0.000)    0.143
   8.099   (  -4.349    8.339    0.000)    9.405
======================= Grid point 27 (15/56) =======================
q-point: (-0.50  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 180
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.135   (  -0.000    1.508    0.000)    1.508
   1.174   (   0.000    2.952    0.000)    2.952
   1.469   (   0.000   -1.705    0.000)    1.705
   1.472   (   0.000    5.352    0.000)    5.352
   1.740   (   0.000  -20.623    0.000)   20.623
   2.116   (   0.000   13.056    0.000)   13.056
   2.337   (   0.000    3.704    0.000)    3.704
   2.480   (   0.000    0.482    0.000)    0.482
   2.653   (   0.000   -0.176    0.000)    0.176
   2.720   (   0.000   28.567    0.000)   28.567
   3.316   (   0.000    7.617    0.000)    7.617
   4.056   (   0.000   -1.113    0.000)    1.113
   5.746   (   0.000  -42.437    0.000)   42.437
   7.108   (   0.000    0.143    0.000)    0.143
   8.053   (   0.000    8.967    0.000)    8.967
======================= Grid point 36 (16/56) =======================
q-point: ( 0.30  0.30  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 180
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.074   (   3.716    3.716    0.000)    5.255
   1.084   (   4.041    4.041    0.000)    5.715
   1.396   (   0.956    0.956    0.000)    1.353
   1.511   (   2.703    2.703    0.000)    3.823
   1.580   (  -9.162   -9.162    0.000)   12.958
   2.281   (   3.528    3.528    0.000)    4.989
   2.398   (   2.197    2.197    0.000)    3.108
   2.404   (   2.281    2.281    0.000)    3.226
   2.545   (   2.871    2.871    0.000)    4.061
   3.113   (  14.548   14.548    0.000)   20.574
   3.198   (  16.642   16.642    0.000)   23.536
   3.691   (   3.573    3.573    0.000)    5.053
   5.537   ( -25.576  -25.576    0.000)   36.171
   7.109   (   0.114    0.114    0.000)    0.161
   8.227   (  -1.506   -1.506    0.000)    2.130
======================= Grid point 37 (17/56) =======================
q-point: ( 0.40  0.30  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.144   (   2.081    2.093    0.000)    2.951
   1.145   (   3.415   -2.781    0.000)    4.404
   1.402   (  -0.444   -2.490    0.000)    2.530
   1.456   (  -3.740  -10.292    0.000)   10.950
   1.554   (   1.610    4.032    0.000)    4.341
   2.315   (   0.601    8.126    0.000)    8.148
   2.409   (   0.346    3.910    0.000)    3.926
   2.465   (   2.489    0.517    0.000)    2.542
   2.602   (   2.581   -1.198    0.000)    2.846
   3.254   (   4.766   25.677    0.000)   26.115
   3.437   (   3.974   10.277    0.000)   11.018
   3.875   (  13.523   -1.498    0.000)   13.606
   5.125   ( -16.524  -34.740    0.000)   38.470
   7.111   (   0.076    0.141    0.000)    0.160
   8.174   (  -3.080    0.290    0.000)    3.094
======================= Grid point 38 (18/56) =======================
q-point: (-0.50  0.30  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 180
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.164   (   0.000    1.418    0.000)    1.418
   1.183   (   0.000   -6.411    0.000)    6.411
   1.391   (  -0.000  -10.950    0.000)   10.950
   1.425   (   0.000   -2.450    0.000)    2.450
   1.570   (   0.000    4.503    0.000)    4.503
   2.319   (   0.000    9.374    0.000)    9.374
   2.413   (   0.000    3.908    0.000)    3.908
   2.490   (   0.000    0.505    0.000)    0.505
   2.629   (   0.000   -1.864    0.000)    1.864
   3.304   (   0.000   29.100    0.000)   29.100
   3.457   (   0.000    6.747    0.000)    6.747
   4.034   (   0.000   -1.190    0.000)    1.190
   4.952   (   0.000  -38.953    0.000)   38.953
   7.111   (   0.000    0.150    0.000)    0.150
   8.140   (   0.000    0.828    0.000)    0.828
======================= Grid point 48 (19/56) =======================
q-point: ( 0.40  0.40  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 186
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.007   (  -8.433   -8.433    0.000)   11.926
   1.175   (   1.127    1.127    0.000)    1.593
   1.356   (  -0.506   -0.506    0.000)    0.716
   1.390   (  -0.789   -0.789    0.000)    1.116
   1.616   (   2.255    2.255    0.000)    3.190
   2.445   (   4.966    4.966    0.000)    7.024
   2.472   (   1.479    1.479    0.000)    2.092
   2.475   (   1.339    1.339    0.000)    1.894
   2.586   (  -0.035   -0.035    0.000)    0.049
   3.508   (   5.668    5.668    0.000)    8.015
   3.727   (  10.132   10.132    0.000)   14.328
   3.879   (   5.455    5.455    0.000)    7.715
   4.557   ( -23.008  -23.008    0.000)   32.538
   7.113   (   0.091    0.091    0.000)    0.129
   8.146   (  -2.223   -2.223    0.000)    3.144
======================= Grid point 49 (20/56) =======================
q-point: (-0.50  0.40  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 180
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.916   (   0.000  -14.473    0.000)   14.473
   1.186   (   0.000    0.834    0.000)    0.834
   1.354   (   0.000    1.287    0.000)    1.287
   1.385   (   0.000   -1.584    0.000)    1.584
   1.638   (   0.000    2.443    0.000)    2.443
   2.477   (   0.000    2.493    0.000)    2.493
   2.490   (   0.000    7.981    0.000)    7.981
   2.498   (   0.000    0.324    0.000)    0.324
   2.595   (   0.000   -1.435    0.000)    1.435
   3.564   (   0.000    4.099    0.000)    4.099
   3.765   (   0.000   17.660    0.000)   17.660
   4.012   (   0.000   -0.938    0.000)    0.938
   4.306   (   0.000  -26.488    0.000)   26.488
   7.114   (   0.000    0.097    0.000)    0.097
   8.120   (   0.000   -1.938    0.000)    1.938
======================= Grid point 65 (21/56) =======================
q-point: (-0.50 -0.50  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 72
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.760   (   0.000    0.000    0.000)    0.000
   1.194   (   0.000    0.000    0.000)    0.000
   1.369   (   0.000    0.000    0.000)    0.000
   1.369   (   0.000    0.000    0.000)    0.000
   1.662   (   0.000    0.000    0.000)    0.000
   2.502   (   0.000    0.000    0.000)    0.000
   2.502   (   0.000    0.000    0.000)    0.000
   2.580   (   0.000    0.000    0.000)    0.000
   2.580   (   0.000    0.000    0.000)    0.000
   3.605   (   0.000    0.000    0.000)    0.000
   3.969   (   0.000    0.000    0.000)    0.000
   4.002   (  -0.000   -0.000    0.000)    0.000
   4.002   (   0.000    0.000    0.000)    0.000
   7.115   (   0.000    0.000    0.000)    0.000
   8.096   (   0.000    0.000    0.000)    0.000
======================= Grid point 133 (22/56) =======================
q-point: ( 0.10  0.10  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 110
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.608   (   8.298    8.298    8.298)   14.373
   0.608   (   8.298    8.298    8.298)   14.373
   1.040   (  11.432   11.432   11.432)   19.800
   1.607   (   9.900    9.900    9.900)   17.147
   1.641   (   7.370    7.370    7.370)   12.766
   1.641   (   7.370    7.370    7.370)   12.766
   2.107   (  -3.292   -3.292   -3.292)    5.701
   2.208   (   0.134    0.134    0.134)    0.232
   2.208   (   0.134    0.134    0.134)    0.232
   2.369   (   5.165    5.165    5.165)    8.947
   2.369   (   5.165    5.165    5.165)    8.947
   3.577   (   0.120    0.120    0.120)    0.208
   6.928   (  -6.210   -6.210   -6.210)   10.757
   6.928   (  -6.210   -6.210   -6.210)   10.757
   8.269   (   3.952    3.952    3.952)    6.845
======================= Grid point 134 (23/56) =======================
q-point: ( 0.20  0.10  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.773   (   9.266    5.571    5.571)   12.163
   0.787   (   9.603    5.096    5.096)   12.006
   1.221   (   7.121    3.248    3.248)    8.474
   1.639   (  -0.254   10.493   10.493)   14.842
   1.685   (   4.699   14.285   14.285)   20.741
   2.055   (  -1.885   -5.723   -5.723)    8.310
   2.066   (  12.981   -3.738   -3.738)   14.016
   2.218   (   0.829    0.146    0.146)    0.855
   2.372   (  15.833    3.812    3.812)   16.726
   2.429   (   9.633    6.844    6.844)   13.655
   2.513   (   7.450    5.880    5.880)   11.164
   3.590   (   1.641   -0.239   -0.239)    1.675
   6.723   (  -9.499  -11.487  -11.487)   18.818
   6.922   (  -0.399   -9.505   -9.505)   13.448
   8.266   (  -3.693    9.379    9.379)   13.768
======================= Grid point 135 (24/56) =======================
q-point: ( 0.30  0.10  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.943   (   8.098    3.164    3.164)    9.252
   0.967   (   8.568    3.640    3.640)    9.995
   1.327   (   4.245   -0.197   -0.197)    4.254
   1.629   (  -0.825   10.111   10.111)   14.323
   1.785   (   5.270   15.066   15.066)   21.948
   2.035   (  -0.342   -8.297   -8.297)   11.738
   2.195   (   4.297   -6.700   -6.700)   10.404
   2.235   (   0.863    0.563    0.563)    1.174
   2.432   (   1.126    8.383    8.383)   11.909
   2.645   (   6.034    7.029    7.029)   11.629
   2.911   (  22.636    4.073    4.073)   23.358
   3.675   (   8.528   -0.818   -0.818)    8.606
   6.549   (  -8.697  -18.176  -18.176)   27.136
   6.914   (  -0.402   -9.905   -9.905)   14.013
   8.161   (  -6.416   12.949   12.949)   19.404
======================= Grid point 136 (25/56) =======================
q-point: ( 0.40  0.10  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.071   (   4.891    1.617    1.617)    5.399
   1.097   (   4.572    2.591    2.591)    5.859
   1.392   (   2.542   -0.076   -0.076)    2.545
   1.610   (  -1.022   10.081   10.081)   14.293
   1.871   (   3.303   14.521   14.521)   20.800
   2.036   (   0.296  -10.093  -10.093)   14.277
   2.249   (   0.511    1.007    1.007)    1.514
   2.262   (   2.586   -7.826   -7.826)   11.366
   2.448   (   0.508    8.793    8.793)   12.446
   2.737   (   3.362    7.631    7.631)   11.304
   3.171   (   5.870    3.770    3.770)    7.930
   3.916   (  13.617   -0.961   -0.961)   13.685
   6.398   (  -6.488  -22.553  -22.553)   32.547
   6.908   (  -0.250  -10.241  -10.241)   14.485
   8.048   (  -4.698   14.791   14.791)   21.439
======================= Grid point 137 (26/56) =======================
q-point: (-0.50  0.10  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 165
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.120   (  -0.000    0.978    0.978)    1.383
   1.140   (   0.000    1.911    1.911)    2.702
   1.418   (   0.000    0.755    0.755)    1.068
   1.598   (   0.000   10.194   10.194)   14.417
   1.904   (   0.000   14.007   14.007)   19.809
   2.040   (   0.000  -10.720  -10.720)   15.161
   2.254   (   0.000    1.204    1.204)    1.702
   2.288   (   0.000   -8.357   -8.357)   11.818
   2.453   (   0.000    9.088    9.088)   12.852
   2.770   (   0.000    7.813    7.813)   11.049
   3.217   (   0.000    3.587    3.587)    5.072
   4.066   (   0.000   -0.814   -0.814)    1.151
   6.329   (   0.000  -24.165  -24.165)   34.175
   6.905   (   0.000  -10.373  -10.373)   14.669
   8.000   (   0.000   15.339   15.339)   21.692
======================= Grid point 145 (27/56) =======================
q-point: ( 0.20  0.20  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.892   (   6.174    6.174    2.882)    9.195
   0.895   (   6.912    6.912    3.516)   10.389
   1.284   (   3.061    3.061    0.902)    4.422
   1.698   (   2.339    2.339   20.561)   20.825
   1.867   (  -4.634   -4.634   -1.087)    6.643
   2.003   (   0.255    0.255   -9.578)    9.585
   2.131   (   4.941    4.941   -1.407)    7.128
   2.406   (  14.838   14.838    4.649)   21.493
   2.427   (   9.103    9.103    6.367)   14.361
   2.594   (  10.801   10.801    5.182)   16.130
   2.632   (   7.566    7.566    5.454)   12.010
   3.593   (   0.915    0.915   -0.907)    1.580
   6.393   ( -18.590  -18.590   -0.719)   26.300
   6.844   (  -1.665   -1.665  -22.520)   22.643
   8.374   (   1.552    1.552    9.850)   10.092
======================= Grid point 146 (28/56) =======================
q-point: ( 0.30  0.20  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 500
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.014   (   6.093    4.001    1.704)    7.486
   1.034   (   7.194    3.405    2.390)    8.310
   1.337   (   2.540    0.984   -0.642)    2.799
   1.724   (   0.048    2.775   15.597)   15.843
   1.792   (  -2.550  -11.931    0.840)   12.229
   2.020   (   2.087    1.291   -7.823)    8.199
   2.181   (   1.607    1.269   -5.703)    6.060
   2.444   (   0.826   18.325    7.025)   19.643
   2.645   (   5.429   15.216    7.227)   17.698
   2.739   (   5.817    5.877    8.801)   12.076
   3.026   (  22.282    7.271    1.251)   23.472
   3.658   (   7.394   -0.673   -1.462)    7.567
   6.059   ( -15.913  -28.589   -2.853)   32.844
   6.820   (  -1.055   -2.301  -24.679)   24.808
   8.346   (  -3.557    5.600   11.504)   13.280
======================= Grid point 147 (29/56) =======================
q-point: ( 0.40  0.20  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 500
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.109   (   3.556    2.265    1.004)    4.334
   1.149   (   4.157    2.484    1.938)    5.216
   1.387   (   2.625   -0.455   -1.191)    2.918
   1.716   (  -0.284   -4.083    7.430)    8.483
   1.754   (  -1.792   -8.735    4.318)    9.907
   2.073   (   2.730    1.706    1.426)    3.521
   2.203   (   0.710   -0.645  -12.278)   12.315
   2.451   (  -0.011   15.436    7.749)   17.272
   2.720   (   2.670   22.677    4.036)   23.187
   2.828   (   3.261    2.765   11.081)   11.878
   3.279   (   5.328    6.999    0.836)    8.836
   3.892   (  13.885   -1.460   -1.165)   14.010
   5.798   ( -10.597  -35.842   -5.543)   37.784
   6.803   (  -0.644   -2.691  -25.646)   25.795
   8.270   (  -3.581    7.590   12.809)   15.314
======================= Grid point 148 (30/56) =======================
q-point: (-0.50  0.20  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.144   (   0.000    1.438    0.916)    1.705
   1.188   (   0.000    2.709    1.432)    3.064
   1.418   (   0.000   -0.650   -1.113)    1.289
   1.722   (   0.000  -19.715   -2.157)   19.833
   1.726   (   0.000    5.254   12.199)   13.283
   2.103   (   0.000    2.630    8.534)    8.931
   2.209   (   0.000   -1.829  -18.171)   18.263
   2.450   (   0.000   14.370    7.972)   16.434
   2.746   (   0.000   24.807    3.040)   24.992
   2.860   (   0.000    1.966   11.795)   11.958
   3.319   (   0.000    6.576    0.510)    6.596
   4.046   (   0.000   -1.149   -0.923)    1.474
   5.688   (   0.000  -38.694   -6.730)   39.275
   6.797   (   0.000   -2.852  -25.954)   26.110
   8.232   (   0.000    8.134   13.295)   15.586
======================= Grid point 157 (31/56) =======================
q-point: ( 0.30  0.30  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.086   (   3.430    3.430    1.251)    5.010
   1.104   (   3.708    3.708    1.475)    5.447
   1.358   (   1.109    1.109   -1.085)    1.907
   1.566   (  -9.196   -9.196   -1.502)   13.091
   1.800   (   2.903    2.903   18.899)   19.340
   2.029   (   1.139    1.139  -15.722)   15.805
   2.196   (   0.437    0.437   -5.880)    5.913
   2.634   (   3.858    3.858   13.538)   14.596
   2.762   (   3.620    3.620   12.793)   13.779
   3.120   (  14.312   14.312    0.730)   20.254
   3.202   (  16.458   16.458    0.455)   23.279
   3.671   (   3.516    3.516   -1.990)    5.356
   5.530   ( -25.454  -25.454   -0.728)   36.004
   6.787   (  -1.351   -1.351  -28.663)   28.726
   8.384   (  -0.923   -0.923   11.932)   12.003
======================= Grid point 158 (32/56) =======================
q-point: ( 0.40  0.30  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 500
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.142   (   2.120   -0.958    0.322)    2.349
   1.169   (   2.906   -0.324    1.757)    3.411
   1.375   (   0.546   -0.915   -1.198)    1.603
   1.438   (  -4.230  -10.697   -1.468)   11.597
   1.832   (   0.547    5.257   15.865)   16.722
   2.077   (   3.433    0.181  -11.306)   11.817
   2.196   (  -0.481   -0.280  -10.432)   10.447
   2.658   (  -0.126    7.788   13.756)   15.808
   2.836   (   3.289   -0.133   13.361)   13.760
   3.257   (   4.609   25.268    0.313)   25.687
   3.431   (   3.407    9.992   -0.597)   10.574
   3.863   (  13.977   -1.453   -1.206)   14.104
   5.115   ( -16.775  -33.656   -1.063)   37.620
   6.766   (  -0.831   -1.501  -30.518)   30.567
   8.344   (  -2.463    0.659   12.772)   13.024
======================= Grid point 159 (33/56) =======================
q-point: (-0.50  0.30  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.162   (   0.000   -2.066   -0.028)    2.066
   1.201   (   0.000   -3.700    1.625)    4.041
   1.380   (  -0.000  -10.452   -1.061)   10.505
   1.397   (   0.000   -1.364   -1.298)    1.883
   1.835   (   0.000    5.999   13.258)   14.552
   2.123   (   0.000    0.719    0.201)    0.746
   2.184   (   0.000   -1.020  -20.194)   20.220
   2.654   (   0.000    8.297   13.763)   16.070
   2.870   (   0.000   -0.443   13.478)   13.486
   3.306   (   0.000   28.856    0.163)   28.856
   3.447   (   0.000    6.210   -1.053)    6.298
   4.025   (   0.000   -1.091   -0.867)    1.394
   4.939   (   0.000  -37.478   -1.426)   37.505
   6.757   (   0.000   -1.563  -31.139)   31.178
   8.316   (   0.000    1.151   13.168)   13.218
======================= Grid point 169 (34/56) =======================
q-point: ( 0.40  0.40  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.002   (  -8.876   -8.876   -0.512)   12.563
   1.189   (   1.083    1.083    1.346)    2.039
   1.349   (  -0.306   -0.306   -0.710)    0.831
   1.373   (  -0.233   -0.233   -0.965)    1.020
   1.915   (   2.847    2.847   21.894)   22.261
   2.101   (   2.399    2.399  -19.751)   20.040
   2.187   (  -0.764   -0.764  -12.209)   12.257
   2.768   (   3.212    3.212   14.420)   15.119
   2.837   (   0.843    0.843   13.781)   13.832
   3.495   (   5.100    5.100   -1.349)    7.337
   3.727   (  10.103   10.103    0.017)   14.287
   3.869   (   5.820    5.820   -0.966)    8.288
   4.560   ( -22.623  -22.623    0.332)   31.996
   6.742   (  -0.915   -0.915  -32.775)   32.801
   8.326   (  -1.672   -1.672   13.373)   13.581
======================= Grid point 170 (35/56) =======================
q-point: (-0.50  0.40  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.907   (   0.000  -14.710   -0.953)   14.741
   1.200   (   0.000    0.762    1.353)    1.553
   1.349   (   0.000    1.237   -0.512)    1.339
   1.373   (   0.000   -0.983   -0.842)    1.295
   1.940   (   0.000    4.293   22.300)   22.710
   2.144   (   0.000    1.377  -12.641)   12.715
   2.168   (  -0.000   -0.596  -20.896)   20.905
   2.789   (   0.000    5.640   14.284)   15.357
   2.857   (   0.000   -0.628   13.910)   13.924
   3.544   (   0.000    3.722   -1.953)    4.203
   3.764   (   0.000   17.742   -0.016)   17.742
   4.006   (   0.000   -0.755   -0.583)    0.954
   4.314   (   0.000  -25.840    0.730)   25.850
   6.733   (   0.000   -0.947  -33.528)   33.542
   8.306   (   0.000   -1.398   13.711)   13.782
======================= Grid point 186 (36/56) =======================
q-point: (-0.50 -0.50  0.10)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 105
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.747   (   0.000    0.000   -1.251)    1.251
   1.208   (   0.000    0.000    1.314)    1.314
   1.363   (   0.000    0.000   -0.606)    0.606
   1.363   (   0.000    0.000   -0.606)    0.606
   1.982   (   0.000    0.000   29.693)   29.693
   2.162   (   0.000    0.000  -21.085)   21.085
   2.162   (   0.000    0.000  -21.085)   21.085
   2.850   (   0.000    0.000   13.965)   13.965
   2.850   (   0.000    0.000   13.965)   13.965
   3.581   (   0.000    0.000   -2.341)    2.341
   3.988   (   0.000    0.000    1.910)    1.910
   3.998   (  -0.000   -0.000   -0.383)    0.383
   3.998   (   0.000    0.000   -0.383)    0.383
   6.723   (   0.000    0.000  -34.322)   34.322
   8.289   (   0.000    0.000   14.032)   14.032
======================= Grid point 266 (37/56) =======================
q-point: ( 0.20  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 116
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.966   (   4.602    4.602    4.602)    7.971
   0.966   (   4.602    4.602    4.602)    7.971
   1.313   (   1.612    1.612    1.612)    2.793
   1.827   (  -6.335   -6.335   -6.335)   10.972
   1.842   (  -2.474   -2.474   -2.474)    4.285
   1.842   (  -2.474   -2.474   -2.474)    4.285
   2.463   (  15.649   15.649   15.649)   27.105
   2.506   (  10.821   10.821   10.821)   18.742
   2.506   (  10.821   10.821   10.821)   18.742
   2.739   (   7.121    7.121    7.121)   12.335
   2.739   (   7.121    7.121    7.121)   12.335
   3.579   (  -0.084   -0.084   -0.084)    0.145
   6.374   ( -12.963  -12.963  -12.963)   22.452
   6.374   ( -12.963  -12.963  -12.963)   22.452
   8.511   (   3.681    3.681    3.681)    6.376
======================= Grid point 267 (38/56) =======================
q-point: ( 0.30  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.058   (   4.680    2.889    2.889)    6.212
   1.079   (   6.723    2.360    2.360)    7.505
   1.339   (   1.260    0.568    0.568)    1.494
   1.707   (  -4.994   -7.899   -7.899)   12.236
   1.797   (  -2.327   -3.911   -3.911)    6.001
   1.832   (   1.045   -7.961   -7.961)   11.307
   2.534   (   1.934   18.335   18.335)   26.002
   2.548   (   1.844   15.761   15.761)   22.366
   2.800   (   3.772    9.783    9.783)   14.340
   2.893   (   6.972    5.900    5.900)   10.874
   3.075   (  22.775    4.228    4.228)   23.547
   3.625   (   6.718   -1.604   -1.604)    7.090
   5.926   ( -21.401  -13.870  -13.870)   29.030
   6.350   (  -1.232  -19.180  -19.180)   27.152
   8.523   (  -1.794    5.909    5.909)    8.547
======================= Grid point 268 (39/56) =======================
q-point: ( 0.40  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 282
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.135   (   3.040    1.643    1.643)    3.826
   1.190   (   3.917    2.027    2.027)    4.854
   1.371   (   2.082   -0.711   -0.711)    2.312
   1.637   (  -2.284  -10.251  -10.251)   14.676
   1.756   (  -1.697   -4.033   -4.033)    5.950
   1.870   (   2.253  -10.742  -10.742)   15.358
   2.553   (   0.188   17.383   17.383)   24.584
   2.576   (   1.051   15.822   15.822)   22.400
   2.852   (   1.873   11.480   11.480)   16.342
   2.998   (   3.769    5.899    5.899)    9.155
   3.315   (   4.366    3.286    3.286)    6.376
   3.864   (  14.519   -1.603   -1.603)   14.695
   5.577   ( -14.104  -19.369  -19.369)   30.810
   6.330   (  -0.782  -19.460  -19.460)   27.532
   8.472   (  -2.685    7.275    7.275)   10.633
======================= Grid point 269 (40/56) =======================
q-point: (-0.50  0.20  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 165
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.166   (   0.000    1.194    1.194)    1.689
   1.225   (   0.000    2.153    2.153)    3.045
   1.398   (   0.000   -1.143   -1.143)    1.616
   1.616   (   0.000  -11.091  -11.091)   15.685
   1.739   (   0.000   -4.078   -4.078)    5.768
   1.896   (   0.000  -11.584  -11.584)   16.382
   2.552   (   0.000   17.032   17.032)   24.087
   2.587   (   0.000   15.809   15.809)   22.358
   2.871   (   0.000   12.038   12.038)   17.025
   3.034   (   0.000    5.905    5.905)    8.351
   3.346   (   0.000    2.686    2.686)    3.799
   4.025   (   0.000   -1.178   -1.178)    1.666
   5.430   (   0.000  -21.681  -21.681)   30.662
   6.322   (   0.000  -19.570  -19.570)   27.676
   8.443   (   0.000    7.710    7.710)   10.903
======================= Grid point 278 (41/56) =======================
q-point: ( 0.30  0.30  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.119   (   2.766    2.766    2.131)    4.455
   1.131   (   3.569    3.569    1.330)    5.220
   1.350   (   0.476    0.476    0.032)    0.675
   1.518   (  -9.110   -9.110   -3.870)   13.452
   1.740   (  -2.141   -2.141  -13.519)   13.854
   1.753   (  -1.499   -1.499  -12.490)   12.669
   2.590   (   1.840    1.840   33.466)   33.567
   2.872   (   4.585    4.585   11.036)   12.801
   2.977   (   4.258    4.258    9.300)   11.079
   3.142   (  13.663   13.663    1.753)   19.402
   3.214   (  15.912   15.912    0.754)   22.516
   3.623   (   3.234    3.234   -2.643)    5.282
   5.511   ( -25.125  -25.125   -1.208)   35.553
   6.151   (  -4.716   -4.716  -35.376)   35.999
   8.574   (  -0.072   -0.072    6.773)    6.774
======================= Grid point 279 (42/56) =======================
q-point: ( 0.40  0.30  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 500
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.154   (   0.977   -1.977    0.888)    2.378
   1.206   (   3.347   -0.615    1.892)    3.894
   1.357   (   0.094   -1.172   -0.808)    1.426
   1.394   (  -3.494   -9.286   -3.338)   10.468
   1.702   (  -1.605   -1.029  -14.023)   14.152
   1.760   (   1.414   -2.894  -18.987)   19.258
   2.615   (   0.729    1.657   34.153)   34.201
   2.901   (   0.014    8.578   11.548)   14.385
   3.058   (   3.376    1.024    9.476)   10.111
   3.266   (   4.239   23.819    0.678)   24.203
   3.416   (   1.966    8.883   -0.816)    9.134
   3.835   (  15.062   -1.264   -1.495)   15.188
   5.079   ( -18.049  -29.197   -3.067)   34.462
   6.094   (  -1.746   -6.902  -36.750)   37.433
   8.548   (  -1.838    1.071    7.534)    7.828
======================= Grid point 280 (43/56) =======================
q-point: (-0.50  0.30  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.162   (   0.000   -5.757   -0.224)    5.762
   1.240   (   0.000   -2.448    2.237)    3.316
   1.349   (   0.000   -8.041   -2.230)    8.345
   1.374   (   0.000   -1.210   -1.167)    1.681
   1.685   (   0.000   -0.821  -14.257)   14.281
   1.779   (   0.000   -2.890  -21.338)   21.533
   2.621   (   0.000    1.630   34.297)   34.335
   2.899   (   0.000    8.979   11.817)   14.841
   3.091   (   0.000    0.707    9.409)    9.436
   3.311   (  -0.000   27.909    0.518)   27.913
   3.420   (   0.000    4.381   -1.533)    4.641
   4.005   (   0.000   -0.859   -1.036)    1.346
   4.888   (   0.000  -31.611   -4.432)   31.920
   6.078   (   0.000   -7.471  -36.803)   37.553
   8.527   (   0.000    1.456    7.827)    7.961
======================= Grid point 290 (44/56) =======================
q-point: ( 0.40  0.40  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 282
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.988   ( -10.032  -10.032   -0.933)   14.217
   1.221   (   1.167    1.167    1.710)    2.377
   1.329   (   0.338    0.338   -1.314)    1.398
   1.355   (  -0.193   -0.193   -0.891)    0.932
   1.709   (   0.440    0.440  -20.701)   20.710
   1.726   (  -0.009   -0.009  -20.705)   20.705
   2.636   (   0.614    0.614   35.506)   35.517
   3.013   (   2.962    2.962   10.984)   11.756
   3.068   (   0.963    0.963   10.256)   10.346
   3.457   (   3.715    3.715   -2.479)    5.809
   3.727   (  10.023   10.023    0.060)   14.174
   3.847   (   6.621    6.621   -1.139)    9.432
   4.569   ( -21.597  -21.597    0.553)   30.547
   6.009   (  -2.619   -2.619  -41.279)   41.444
   8.539   (  -1.220   -1.220    7.952)    8.137
======================= Grid point 291 (45/56) =======================
q-point: (-0.50  0.40  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.882   (   0.000  -15.323   -1.628)   15.409
   1.233   (   0.000    0.630    1.888)    1.991
   1.336   (   0.000    1.078   -0.805)    1.345
   1.355   (   0.000   -0.754   -1.048)    1.291
   1.704   (   0.000    2.095  -21.138)   21.242
   1.742   (   0.000   -1.162  -23.133)   23.162
   2.641   (   0.000    0.575   35.875)   35.879
   3.030   (   0.000    5.009   10.815)   11.919
   3.089   (   0.000   -0.511   10.147)   10.160
   3.493   (   0.000    2.830   -3.268)    4.323
   3.764   (   0.000   17.863    0.001)   17.863
   3.993   (   0.000   -0.410   -0.731)    0.838
   4.332   (   0.000  -23.865    1.096)   23.890
   5.984   (   0.000   -2.889  -42.100)   42.199
   8.524   (   0.000   -0.997    8.149)    8.210
======================= Grid point 307 (46/56) =======================
q-point: (-0.50 -0.50  0.20)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.714   (   0.000    0.000   -2.119)    2.119
   1.240   (   0.000    0.000    1.881)    1.881
   1.347   (   0.000    0.000   -0.984)    0.984
   1.347   (   0.000    0.000   -0.984)    0.984
   1.731   (   0.000    0.000  -23.545)   23.545
   1.731   (   0.000    0.000  -23.545)   23.545
   2.646   (   0.000    0.000   36.273)   36.273
   3.082   (   0.000    0.000   10.167)   10.167
   3.082   (   0.000    0.000   10.167)   10.167
   3.521   (   0.000    0.000   -3.768)    3.768
   3.989   (   0.000    0.000   -0.558)    0.558
   3.989   (   0.000    0.000   -0.558)    0.558
   4.037   (  -0.000   -0.000    3.052)    3.052
   5.957   (   0.000    0.000  -43.170)   43.170
   8.511   (   0.000    0.000    8.316)    8.316
======================= Grid point 399 (47/56) =======================
q-point: ( 0.30  0.30  0.30)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 110
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.162   (   2.069    2.069    2.069)    3.584
   1.162   (   2.069    2.069    2.069)    3.584
   1.352   (  -0.039   -0.039   -0.039)    0.068
   1.395   (  -8.435   -8.435   -8.435)   14.610
   1.552   (  -6.671   -6.671   -6.671)   11.554
   1.552   (  -6.671   -6.671   -6.671)   11.554
   3.054   (   6.073    6.073    6.073)   10.519
   3.054   (   6.073    6.073    6.073)   10.519
   3.204   (   8.621    8.621    8.621)   14.932
   3.229   (  10.231   10.231   10.231)   17.721
   3.229   (  10.231   10.231   10.231)   17.721
   3.596   (   1.675    1.675    1.675)    2.901
   5.487   ( -16.885  -16.885  -16.885)   29.246
   5.487   ( -16.885  -16.885  -16.885)   29.246
   8.642   (   0.817    0.817    0.817)    1.415
======================= Grid point 400 (48/56) =======================
q-point: ( 0.40  0.30  0.30)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 280
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.160   (  -2.853   -2.636   -2.636)    4.694
   1.227   (   4.176   -0.853   -0.853)    4.347
   1.285   (  -1.286   -6.003   -6.003)    8.586
   1.347   (  -0.201   -0.533   -0.533)    0.781
   1.461   (  -2.851   -9.439   -9.439)   13.650
   1.472   (  -3.023   -9.299   -9.299)   13.493
   3.070   (  -0.377    8.927    8.927)   12.630
   3.170   (   4.373    6.223    6.223)    9.827
   3.238   (   0.102   13.611   13.611)   19.250
   3.309   (   4.592    9.049    9.049)   13.597
   3.424   (   0.133    4.876    4.876)    6.898
   3.811   (  15.893   -0.892   -0.892)   15.943
   4.952   ( -21.751  -12.867  -12.867)   28.360
   5.466   (  -0.795  -24.348  -24.348)   34.442
   8.631   (  -1.299    1.503    1.503)    2.491
======================= Grid point 401 (49/56) =======================
q-point: (-0.50  0.30  0.30)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 165
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.118   (   0.000   -7.348   -7.348)   10.392
   1.277   (   0.000   -0.423   -0.423)    0.599
   1.288   (   0.000   -2.136   -2.136)    3.020
   1.352   (  -0.000   -1.107   -1.107)    1.566
   1.427   (   0.000  -11.760  -11.760)   16.631
   1.439   (   0.000  -10.447  -10.447)   14.774
   3.064   (   0.000    9.167    9.167)   12.965
   3.221   (  -0.000    3.820    3.820)    5.402
   3.230   (   0.000   15.857   15.857)   22.426
   3.368   (  -0.000   13.118   13.118)   18.551
   3.402   (   0.000    0.199    0.199)    0.281
   3.988   (   0.000   -0.636   -0.636)    0.899
   4.722   (   0.000  -14.672  -14.672)   20.749
   5.458   (   0.000  -24.208  -24.208)   34.236
   8.615   (   0.000    1.761    1.761)    2.490
======================= Grid point 411 (50/56) =======================
q-point: ( 0.40  0.40  0.30)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.966   ( -11.260  -11.260   -1.403)   15.985
   1.229   (   1.431    1.431   -3.921)    4.413
   1.268   (   1.072    1.072  -10.827)   10.932
   1.337   (  -0.773   -0.773   -2.158)    2.419
   1.353   (  -2.817   -2.817   -9.848)   10.623
   1.362   (  -3.394   -3.394  -12.141)   13.056
   3.189   (   2.477    2.477    7.336)    8.130
   3.208   (   1.984    1.984   10.061)   10.445
   3.296   (   0.389    0.389   25.477)   25.483
   3.408   (   1.902    1.902   -2.074)    3.396
   3.729   (   9.928    9.928    0.106)   14.041
   3.830   (   7.173    7.173   -0.629)   10.164
   4.580   ( -20.289  -20.289    0.574)   28.699
   5.225   (  -4.591   -4.591  -38.676)   39.217
   8.630   (  -0.882   -0.882    1.903)    2.275
======================= Grid point 412 (51/56) =======================
q-point: (-0.50  0.40  0.30)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 300
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.847   (   0.000  -15.662   -1.986)   15.787
   1.250   (   0.000    0.439   -6.955)    6.969
   1.302   (   0.000   -3.912  -19.863)   20.245
   1.302   (  -0.000   -0.617  -10.487)   10.505
   1.328   (   0.000   -1.091   -4.871)    4.992
   1.334   (   0.000   -0.691   -1.116)    1.312
   3.193   (   0.000    4.723    8.440)    9.672
   3.247   (  -0.000   -0.204    6.361)    6.365
   3.295   (   0.000    0.221   28.448)   28.448
   3.427   (   0.000    1.722   -3.252)    3.680
   3.765   (   0.000   17.766    0.104)   17.766
   3.981   (   0.000   -0.221   -0.513)    0.559
   4.349   (   0.000  -20.028    0.431)   20.032
   5.188   (   0.000   -6.318  -38.853)   39.364
   8.618   (   0.000   -0.722    2.074)    2.196
======================= Grid point 428 (52/56) =======================
q-point: (-0.50 -0.50  0.30)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 105
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.671   (   0.000    0.000   -2.256)    2.256
   1.270   (   0.000    0.000  -24.060)   24.060
   1.270   (   0.000    0.000  -24.060)   24.060
   1.274   (  -0.000   -0.000    1.581)    1.581
   1.328   (   0.000    0.000   -1.006)    1.006
   1.328   (   0.000    0.000   -1.006)    1.006
   3.242   (   0.000    0.000    6.492)    6.492
   3.242   (   0.000    0.000    6.492)    6.492
   3.294   (   0.000    0.000   29.845)   29.845
   3.446   (   0.000    0.000   -3.735)    3.735
   3.979   (   0.000    0.000   -0.479)    0.479
   3.979   (   0.000    0.000   -0.479)    0.479
   4.097   (  -0.000   -0.000    3.008)    3.008
   5.136   (   0.000    0.000  -40.570)   40.570
   8.608   (   0.000    0.000    2.217)    2.217
======================= Grid point 532 (53/56) =======================
q-point: ( 0.40  0.40  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 116
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.901   (  -8.024   -8.024   -8.024)   13.897
   1.010   (  -8.778   -8.778   -8.778)   15.204
   1.010   (  -8.778   -8.778   -8.778)   15.204
   1.302   (  -0.148   -0.148   -0.148)    0.257
   1.302   (  -0.148   -0.148   -0.148)    0.257
   1.327   (  -0.613   -0.613   -0.613)    1.061
   3.294   (   2.559    2.559    2.559)    4.433
   3.294   (   2.559    2.559    2.559)    4.433
   3.376   (  -0.095   -0.095   -0.095)    0.165
   3.731   (   6.609    6.609    6.609)   11.447
   3.731   (   6.609    6.609    6.609)   11.447
   3.840   (   4.407    4.407    4.407)    7.632
   4.590   ( -12.678  -12.678  -12.678)   21.959
   4.590   ( -12.678  -12.678  -12.678)   21.959
   8.636   (  -0.603   -0.603   -0.603)    1.044
======================= Grid point 533 (54/56) =======================
q-point: (-0.50  0.40  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 165
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.783   (   0.000   -7.678   -7.678)   10.858
   0.903   (   0.000  -12.021  -12.021)   17.000
   0.972   (   0.000  -13.250  -13.250)   18.739
   1.294   (   0.000    0.569    0.569)    0.805
   1.312   (   0.000   -0.228   -0.228)    0.322
   1.318   (  -0.000   -0.631   -0.631)    0.893
   3.300   (   0.000    3.561    3.561)    5.036
   3.329   (   0.000    1.821    1.821)    2.575
   3.382   (   0.000   -0.813   -0.813)    1.150
   3.732   (   0.000    9.851    9.851)   13.932
   3.789   (   0.000    8.216    8.216)   11.619
   3.974   (   0.000   -0.225   -0.225)    0.318
   4.312   (   0.000   -6.991   -6.991)    9.887
   4.593   (   0.000  -18.777  -18.777)   26.555
   8.627   (   0.000   -0.487   -0.487)    0.688
======================= Grid point 549 (55/56) =======================
q-point: (-0.50 -0.50  0.40)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.633   (   0.000    0.000   -1.476)    1.476
   0.837   (   0.000    0.000  -19.548)   19.548
   0.837   (   0.000    0.000  -19.548)   19.548
   1.298   (   0.000    0.000    0.826)    0.826
   1.312   (   0.000    0.000   -0.618)    0.618
   1.312   (   0.000    0.000   -0.618)    0.618
   3.334   (   0.000    0.000    3.128)    3.128
   3.334   (   0.000    0.000    3.128)    3.128
   3.386   (   0.000    0.000   -2.287)    2.287
   3.758   (   0.000    0.000   18.054)   18.054
   3.971   (   0.000    0.000   -0.257)    0.257
   3.971   (   0.000    0.000   -0.257)    0.257
   4.145   (   0.000    0.000    1.837)    1.837
   4.469   (   0.000    0.000  -27.076)   27.076
   8.620   (   0.000    0.000   -0.385)    0.385
======================= Grid point 735 (56/56) =======================
q-point: (-0.50 -0.50 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 2.72e-04 5.83e-04 5.83e-04 5.83e-04 
Number of triplets: 28
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.618   (   0.000    0.000    0.000)    0.000
   0.618   (   0.000    0.000    0.000)    0.000
   0.618   (   0.000    0.000    0.000)    0.000
   1.305   (   0.000    0.000    0.000)    0.000
   1.305   (   0.000    0.000    0.000)    0.000
   1.305   (   0.000    0.000    0.000)    0.000
   3.364   (   0.000    0.000    0.000)    0.000
   3.364   (   0.000    0.000    0.000)    0.000
   3.364   (   0.000    0.000    0.000)    0.000
   3.969   (   0.000    0.000    0.000)    0.000
   3.969   (   0.000    0.000    0.000)    0.000
   3.969   (   0.000    0.000    0.000)    0.000
   4.164   (   0.000    0.000    0.000)    0.000
   4.164   (   0.000    0.000    0.000)    0.000
   8.613   (   0.000    0.000    0.000)    0.000
=================== 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/15000
   10.0     17.174     17.174     17.174     -0.000      0.000      0.000 3/15000
   20.0      4.690      4.690      4.690      0.000      0.000      0.000 3/15000
   30.0      3.115      3.115      3.115      0.000      0.000      0.000 3/15000
   40.0      2.547      2.547      2.547      0.000      0.000      0.000 3/15000
   50.0      2.214      2.214      2.214      0.000      0.000      0.000 3/15000
   60.0      1.976      1.976      1.976      0.000      0.000      0.000 3/15000
   70.0      1.789      1.789      1.789      0.000      0.000      0.000 3/15000
   80.0      1.635      1.635      1.635      0.000      0.000      0.000 3/15000
   90.0      1.504      1.504      1.504      0.000      0.000      0.000 3/15000
  100.0      1.392      1.392      1.392      0.000      0.000      0.000 3/15000
  110.0      1.294      1.294      1.294      0.000      0.000      0.000 3/15000
  120.0      1.208      1.208      1.208      0.000      0.000      0.000 3/15000
  130.0      1.132      1.132      1.132      0.000      0.000      0.000 3/15000
  140.0      1.065      1.065      1.065      0.000      0.000      0.000 3/15000
  150.0      1.004      1.004      1.004      0.000      0.000      0.000 3/15000
  160.0      0.950      0.950      0.950      0.000     -0.000      0.000 3/15000
  170.0      0.901      0.901      0.901      0.000     -0.000      0.000 3/15000
  180.0      0.857      0.857      0.857      0.000     -0.000      0.000 3/15000
  190.0      0.816      0.816      0.816      0.000     -0.000      0.000 3/15000
  200.0      0.779      0.779      0.779      0.000     -0.000      0.000 3/15000
  210.0      0.745      0.745      0.745      0.000     -0.000      0.000 3/15000
  220.0      0.714      0.714      0.714      0.000     -0.000      0.000 3/15000
  230.0      0.685      0.685      0.685      0.000     -0.000      0.000 3/15000
  240.0      0.659      0.659      0.659      0.000     -0.000      0.000 3/15000
  250.0      0.634      0.634      0.634      0.000     -0.000      0.000 3/15000
  260.0      0.611      0.611      0.611      0.000     -0.000      0.000 3/15000
  270.0      0.590      0.590      0.590      0.000     -0.000      0.000 3/15000
  280.0      0.570      0.570      0.570      0.000     -0.000      0.000 3/15000
  290.0      0.551      0.551      0.551      0.000     -0.000      0.000 3/15000
  300.0      0.534      0.534      0.534      0.000     -0.000      0.000 3/15000
  310.0      0.517      0.517      0.517      0.000     -0.000      0.000 3/15000
  320.0      0.502      0.502      0.502      0.000     -0.000      0.000 3/15000
  330.0      0.487      0.487      0.487      0.000     -0.000      0.000 3/15000
  340.0      0.473      0.473      0.473      0.000     -0.000      0.000 3/15000
  350.0      0.460      0.460      0.460      0.000     -0.000      0.000 3/15000
  360.0      0.448      0.448      0.448      0.000     -0.000      0.000 3/15000
  370.0      0.436      0.436      0.436      0.000     -0.000      0.000 3/15000
  380.0      0.425      0.425      0.425      0.000     -0.000      0.000 3/15000
  390.0      0.415      0.415      0.415      0.000     -0.000      0.000 3/15000
  400.0      0.405      0.405      0.405      0.000     -0.000      0.000 3/15000
  410.0      0.395      0.395      0.395      0.000     -0.000      0.000 3/15000
  420.0      0.386      0.386      0.386      0.000     -0.000      0.000 3/15000
  430.0      0.377      0.377      0.377      0.000     -0.000      0.000 3/15000
  440.0      0.369      0.369      0.369      0.000     -0.000      0.000 3/15000
  450.0      0.361      0.361      0.361      0.000     -0.000      0.000 3/15000
  460.0      0.353      0.353      0.353      0.000     -0.000      0.000 3/15000
  470.0      0.346      0.346      0.346      0.000     -0.000      0.000 3/15000
  480.0      0.339      0.339      0.339      0.000     -0.000      0.000 3/15000
  490.0      0.332      0.332      0.332      0.000     -0.000      0.000 3/15000
  500.0      0.325      0.325      0.325      0.000     -0.000      0.000 3/15000
  510.0      0.319      0.319      0.319      0.000     -0.000      0.000 3/15000
  520.0      0.313      0.313      0.313      0.000     -0.000      0.000 3/15000
  530.0      0.307      0.307      0.307      0.000     -0.000      0.000 3/15000
  540.0      0.302      0.302      0.302      0.000     -0.000      0.000 3/15000
  550.0      0.296      0.296      0.296      0.000     -0.000      0.000 3/15000
  560.0      0.291      0.291      0.291      0.000     -0.000      0.000 3/15000
  570.0      0.286      0.286      0.286      0.000     -0.000      0.000 3/15000
  580.0      0.281      0.281      0.281      0.000     -0.000      0.000 3/15000
  590.0      0.276      0.276      0.276      0.000     -0.000      0.000 3/15000
  600.0      0.272      0.272      0.272      0.000     -0.000      0.000 3/15000
  610.0      0.268      0.268      0.268      0.000     -0.000      0.000 3/15000
  620.0      0.263      0.263      0.263      0.000     -0.000      0.000 3/15000
  630.0      0.259      0.259      0.259      0.000     -0.000      0.000 3/15000
  640.0      0.255      0.255      0.255      0.000     -0.000      0.000 3/15000
  650.0      0.251      0.251      0.251      0.000     -0.000      0.000 3/15000
  660.0      0.247      0.247      0.247      0.000     -0.000      0.000 3/15000
  670.0      0.244      0.244      0.244      0.000     -0.000      0.000 3/15000
  680.0      0.240      0.240      0.240      0.000     -0.000      0.000 3/15000
  690.0      0.237      0.237      0.237      0.000     -0.000      0.000 3/15000
  700.0      0.234      0.234      0.234      0.000     -0.000      0.000 3/15000
  710.0      0.230      0.230      0.230      0.000     -0.000      0.000 3/15000
  720.0      0.227      0.227      0.227      0.000     -0.000      0.000 3/15000
  730.0      0.224      0.224      0.224      0.000     -0.000      0.000 3/15000
  740.0      0.221      0.221      0.221      0.000     -0.000      0.000 3/15000
  750.0      0.218      0.218      0.218      0.000     -0.000      0.000 3/15000
  760.0      0.215      0.215      0.215      0.000     -0.000      0.000 3/15000
  770.0      0.212      0.212      0.212      0.000     -0.000      0.000 3/15000
  780.0      0.210      0.210      0.210      0.000     -0.000      0.000 3/15000
  790.0      0.207      0.207      0.207      0.000     -0.000      0.000 3/15000
  800.0      0.205      0.205      0.205      0.000     -0.000      0.000 3/15000
  810.0      0.202      0.202      0.202      0.000     -0.000      0.000 3/15000
  820.0      0.200      0.200      0.200      0.000     -0.000      0.000 3/15000
  830.0      0.197      0.197      0.197      0.000     -0.000      0.000 3/15000
  840.0      0.195      0.195      0.195      0.000     -0.000      0.000 3/15000
  850.0      0.193      0.193      0.193      0.000     -0.000      0.000 3/15000
  860.0      0.190      0.190      0.190      0.000     -0.000      0.000 3/15000
  870.0      0.188      0.188      0.188      0.000     -0.000      0.000 3/15000
  880.0      0.186      0.186      0.186      0.000     -0.000      0.000 3/15000
  890.0      0.184      0.184      0.184      0.000     -0.000      0.000 3/15000
  900.0      0.182      0.182      0.182      0.000     -0.000      0.000 3/15000
  910.0      0.180      0.180      0.180      0.000     -0.000      0.000 3/15000
  920.0      0.178      0.178      0.178      0.000     -0.000      0.000 3/15000
  930.0      0.176      0.176      0.176      0.000     -0.000      0.000 3/15000
  940.0      0.174      0.174      0.174      0.000     -0.000      0.000 3/15000
  950.0      0.172      0.172      0.172      0.000     -0.000      0.000 3/15000
  960.0      0.171      0.171      0.171      0.000     -0.000      0.000 3/15000
  970.0      0.169      0.169      0.169      0.000     -0.000      0.000 3/15000
  980.0      0.167      0.167      0.167      0.000     -0.000      0.000 3/15000
  990.0      0.166      0.166      0.166      0.000     -0.000      0.000 3/15000
 1000.0      0.164      0.164      0.164      0.000     -0.000      0.000 3/15000

Thermal conductivity related properties were written into 
"kappa-m101010.hdf5".
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 18:05:49]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

