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

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

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

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

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

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


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

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

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


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

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

-------------------------[time 2026-01-07 22:53:13]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: conductivity-RTA
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py.yaml".
Supercell (dim): [3 3 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: Pm-3m (221)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    4.470904309999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.470904309999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.470904309999999
Atomic positions (fractional):
    1 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998
    2 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998
    3 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998
    4 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411
    5 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.412712929999998    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.412712929999998    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.412712929999998
Atomic positions (fractional):
    1 F   0.00000000000000  0.16666666666667  0.16666666666667  18.998 > 1
    2 F   0.33333333333333  0.16666666666667  0.16666666666667  18.998 > 1
    3 F   0.66666666666667  0.16666666666667  0.16666666666667  18.998 > 1
    4 F   0.00000000000000  0.50000000000000  0.16666666666667  18.998 > 1
    5 F   0.33333333333333  0.50000000000000  0.16666666666667  18.998 > 1
    6 F   0.66666666666667  0.50000000000000  0.16666666666667  18.998 > 1
    7 F   0.00000000000000  0.83333333333333  0.16666666666667  18.998 > 1
    8 F   0.33333333333333  0.83333333333333  0.16666666666667  18.998 > 1
    9 F   0.66666666666667  0.83333333333333  0.16666666666667  18.998 > 1
   10 F   0.00000000000000  0.16666666666667  0.50000000000000  18.998 > 1
   11 F   0.33333333333333  0.16666666666667  0.50000000000000  18.998 > 1
   12 F   0.66666666666667  0.16666666666667  0.50000000000000  18.998 > 1
   13 F   0.00000000000000  0.50000000000000  0.50000000000000  18.998 > 1
   14 F   0.33333333333333  0.50000000000000  0.50000000000000  18.998 > 1
   15 F   0.66666666666667  0.50000000000000  0.50000000000000  18.998 > 1
   16 F   0.00000000000000  0.83333333333333  0.50000000000000  18.998 > 1
   17 F   0.33333333333333  0.83333333333333  0.50000000000000  18.998 > 1
   18 F   0.66666666666667  0.83333333333333  0.50000000000000  18.998 > 1
   19 F   0.00000000000000  0.16666666666667  0.83333333333333  18.998 > 1
   20 F   0.33333333333333  0.16666666666667  0.83333333333333  18.998 > 1
   21 F   0.66666666666667  0.16666666666667  0.83333333333333  18.998 > 1
   22 F   0.00000000000000  0.50000000000000  0.83333333333333  18.998 > 1
   23 F   0.33333333333333  0.50000000000000  0.83333333333333  18.998 > 1
   24 F   0.66666666666667  0.50000000000000  0.83333333333333  18.998 > 1
   25 F   0.00000000000000  0.83333333333333  0.83333333333333  18.998 > 1
   26 F   0.33333333333333  0.83333333333333  0.83333333333333  18.998 > 1
   27 F   0.66666666666667  0.83333333333333  0.83333333333333  18.998 > 1
   28 F   0.16666666666667  0.00000000000000  0.16666666666667  18.998 > 28
   29 F   0.50000000000000  0.00000000000000  0.16666666666667  18.998 > 28
   30 F   0.83333333333333  0.00000000000000  0.16666666666667  18.998 > 28
   31 F   0.16666666666667  0.33333333333333  0.16666666666667  18.998 > 28
   32 F   0.50000000000000  0.33333333333333  0.16666666666667  18.998 > 28
   33 F   0.83333333333333  0.33333333333333  0.16666666666667  18.998 > 28
   34 F   0.16666666666667  0.66666666666667  0.16666666666667  18.998 > 28
   35 F   0.50000000000000  0.66666666666667  0.16666666666667  18.998 > 28
   36 F   0.83333333333333  0.66666666666667  0.16666666666667  18.998 > 28
   37 F   0.16666666666667  0.00000000000000  0.50000000000000  18.998 > 28
   38 F   0.50000000000000  0.00000000000000  0.50000000000000  18.998 > 28
   39 F   0.83333333333333  0.00000000000000  0.50000000000000  18.998 > 28
   40 F   0.16666666666667  0.33333333333333  0.50000000000000  18.998 > 28
   41 F   0.50000000000000  0.33333333333333  0.50000000000000  18.998 > 28
   42 F   0.83333333333333  0.33333333333333  0.50000000000000  18.998 > 28
   43 F   0.16666666666667  0.66666666666667  0.50000000000000  18.998 > 28
   44 F   0.50000000000000  0.66666666666667  0.50000000000000  18.998 > 28
   45 F   0.83333333333333  0.66666666666667  0.50000000000000  18.998 > 28
   46 F   0.16666666666667  0.00000000000000  0.83333333333333  18.998 > 28
   47 F   0.50000000000000  0.00000000000000  0.83333333333333  18.998 > 28
   48 F   0.83333333333333  0.00000000000000  0.83333333333333  18.998 > 28
   49 F   0.16666666666667  0.33333333333333  0.83333333333333  18.998 > 28
   50 F   0.50000000000000  0.33333333333333  0.83333333333333  18.998 > 28
   51 F   0.83333333333333  0.33333333333333  0.83333333333333  18.998 > 28
   52 F   0.16666666666667  0.66666666666667  0.83333333333333  18.998 > 28
   53 F   0.50000000000000  0.66666666666667  0.83333333333333  18.998 > 28
   54 F   0.83333333333333  0.66666666666667  0.83333333333333  18.998 > 28
   55 F   0.16666666666667  0.16666666666667  0.00000000000000  18.998 > 55
   56 F   0.50000000000000  0.16666666666667  0.00000000000000  18.998 > 55
   57 F   0.83333333333333  0.16666666666667  0.00000000000000  18.998 > 55
   58 F   0.16666666666667  0.50000000000000  0.00000000000000  18.998 > 55
   59 F   0.50000000000000  0.50000000000000  0.00000000000000  18.998 > 55
   60 F   0.83333333333333  0.50000000000000  0.00000000000000  18.998 > 55
   61 F   0.16666666666667  0.83333333333333  0.00000000000000  18.998 > 55
   62 F   0.50000000000000  0.83333333333333  0.00000000000000  18.998 > 55
   63 F   0.83333333333333  0.83333333333333  0.00000000000000  18.998 > 55
   64 F   0.16666666666667  0.16666666666667  0.33333333333333  18.998 > 55
   65 F   0.50000000000000  0.16666666666667  0.33333333333333  18.998 > 55
   66 F   0.83333333333333  0.16666666666667  0.33333333333333  18.998 > 55
   67 F   0.16666666666667  0.50000000000000  0.33333333333333  18.998 > 55
   68 F   0.50000000000000  0.50000000000000  0.33333333333333  18.998 > 55
   69 F   0.83333333333333  0.50000000000000  0.33333333333333  18.998 > 55
   70 F   0.16666666666667  0.83333333333333  0.33333333333333  18.998 > 55
   71 F   0.50000000000000  0.83333333333333  0.33333333333333  18.998 > 55
   72 F   0.83333333333333  0.83333333333333  0.33333333333333  18.998 > 55
   73 F   0.16666666666667  0.16666666666667  0.66666666666667  18.998 > 55
   74 F   0.50000000000000  0.16666666666667  0.66666666666667  18.998 > 55
   75 F   0.83333333333333  0.16666666666667  0.66666666666667  18.998 > 55
   76 F   0.16666666666667  0.50000000000000  0.66666666666667  18.998 > 55
   77 F   0.50000000000000  0.50000000000000  0.66666666666667  18.998 > 55
   78 F   0.83333333333333  0.50000000000000  0.66666666666667  18.998 > 55
   79 F   0.16666666666667  0.83333333333333  0.66666666666667  18.998 > 55
   80 F   0.50000000000000  0.83333333333333  0.66666666666667  18.998 > 55
   81 F   0.83333333333333  0.83333333333333  0.66666666666667  18.998 > 55
   82 Cd  0.16666666666667  0.16666666666667  0.16666666666667 112.411 > 82
   83 Cd  0.50000000000000  0.16666666666667  0.16666666666667 112.411 > 82
   84 Cd  0.83333333333333  0.16666666666667  0.16666666666667 112.411 > 82
   85 Cd  0.16666666666667  0.50000000000000  0.16666666666667 112.411 > 82
   86 Cd  0.50000000000000  0.50000000000000  0.16666666666667 112.411 > 82
   87 Cd  0.83333333333333  0.50000000000000  0.16666666666667 112.411 > 82
   88 Cd  0.16666666666667  0.83333333333333  0.16666666666667 112.411 > 82
   89 Cd  0.50000000000000  0.83333333333333  0.16666666666667 112.411 > 82
   90 Cd  0.83333333333333  0.83333333333333  0.16666666666667 112.411 > 82
   91 Cd  0.16666666666667  0.16666666666667  0.50000000000000 112.411 > 82
   92 Cd  0.50000000000000  0.16666666666667  0.50000000000000 112.411 > 82
   93 Cd  0.83333333333333  0.16666666666667  0.50000000000000 112.411 > 82
   94 Cd  0.16666666666667  0.50000000000000  0.50000000000000 112.411 > 82
   95 Cd  0.50000000000000  0.50000000000000  0.50000000000000 112.411 > 82
   96 Cd  0.83333333333333  0.50000000000000  0.50000000000000 112.411 > 82
   97 Cd  0.16666666666667  0.83333333333333  0.50000000000000 112.411 > 82
   98 Cd  0.50000000000000  0.83333333333333  0.50000000000000 112.411 > 82
   99 Cd  0.83333333333333  0.83333333333333  0.50000000000000 112.411 > 82
  100 Cd  0.16666666666667  0.16666666666667  0.83333333333333 112.411 > 82
  101 Cd  0.50000000000000  0.16666666666667  0.83333333333333 112.411 > 82
  102 Cd  0.83333333333333  0.16666666666667  0.83333333333333 112.411 > 82
  103 Cd  0.16666666666667  0.50000000000000  0.83333333333333 112.411 > 82
  104 Cd  0.50000000000000  0.50000000000000  0.83333333333333 112.411 > 82
  105 Cd  0.83333333333333  0.50000000000000  0.83333333333333 112.411 > 82
  106 Cd  0.16666666666667  0.83333333333333  0.83333333333333 112.411 > 82
  107 Cd  0.50000000000000  0.83333333333333  0.83333333333333 112.411 > 82
  108 Cd  0.83333333333333  0.83333333333333  0.83333333333333 112.411 > 82
  109 Cs  0.00000000000000  0.00000000000000  0.00000000000000 132.905 > 109
  110 Cs  0.33333333333333  0.00000000000000  0.00000000000000 132.905 > 109
  111 Cs  0.66666666666667  0.00000000000000  0.00000000000000 132.905 > 109
  112 Cs  0.00000000000000  0.33333333333333  0.00000000000000 132.905 > 109
  113 Cs  0.33333333333333  0.33333333333333  0.00000000000000 132.905 > 109
  114 Cs  0.66666666666667  0.33333333333333  0.00000000000000 132.905 > 109
  115 Cs  0.00000000000000  0.66666666666667  0.00000000000000 132.905 > 109
  116 Cs  0.33333333333333  0.66666666666667  0.00000000000000 132.905 > 109
  117 Cs  0.66666666666667  0.66666666666667  0.00000000000000 132.905 > 109
  118 Cs  0.00000000000000  0.00000000000000  0.33333333333333 132.905 > 109
  119 Cs  0.33333333333333  0.00000000000000  0.33333333333333 132.905 > 109
  120 Cs  0.66666666666667  0.00000000000000  0.33333333333333 132.905 > 109
  121 Cs  0.00000000000000  0.33333333333333  0.33333333333333 132.905 > 109
  122 Cs  0.33333333333333  0.33333333333333  0.33333333333333 132.905 > 109
  123 Cs  0.66666666666667  0.33333333333333  0.33333333333333 132.905 > 109
  124 Cs  0.00000000000000  0.66666666666667  0.33333333333333 132.905 > 109
  125 Cs  0.33333333333333  0.66666666666667  0.33333333333333 132.905 > 109
  126 Cs  0.66666666666667  0.66666666666667  0.33333333333333 132.905 > 109
  127 Cs  0.00000000000000  0.00000000000000  0.66666666666667 132.905 > 109
  128 Cs  0.33333333333333  0.00000000000000  0.66666666666667 132.905 > 109
  129 Cs  0.66666666666667  0.00000000000000  0.66666666666667 132.905 > 109
  130 Cs  0.00000000000000  0.33333333333333  0.66666666666667 132.905 > 109
  131 Cs  0.33333333333333  0.33333333333333  0.66666666666667 132.905 > 109
  132 Cs  0.66666666666667  0.33333333333333  0.66666666666667 132.905 > 109
  133 Cs  0.00000000000000  0.66666666666667  0.66666666666667 132.905 > 109
  134 Cs  0.33333333333333  0.66666666666667  0.66666666666667 132.905 > 109
  135 Cs  0.66666666666667  0.66666666666667  0.66666666666667 132.905 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.6841910    0.0000000    0.0000000
            0.0000000    2.6841910    0.0000000
            0.0000000    0.0000000    2.6841910
-------------------------- Born effective charges --------------------------
    1 F    -1.7199014    0.0000000    0.0000000
            0.0000000   -1.0332342    0.0000000
            0.0000000    0.0000000   -1.0332342
    2 F    -1.0332342    0.0000000    0.0000000
            0.0000000   -1.7199014    0.0000000
            0.0000000    0.0000000   -1.0332342
    3 F    -1.0332342    0.0000000    0.0000000
            0.0000000   -1.0332342    0.0000000
            0.0000000    0.0000000   -1.7199014
    4 Cd    2.4345063    0.0000000    0.0000000
            0.0000000    2.4345063    0.0000000
            0.0000000    0.0000000    2.4345063
    5 Cs    1.3518636    0.0000000    0.0000000
            0.0000000    1.3518636    0.0000000
            0.0000000    0.0000000    1.3518636
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (yyy) 0.00000000 (yyy) 0.00000000 (yyy)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 11 11 11 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.93, 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 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 56
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   2.318   (   0.000    0.000    0.000)    0.000
   2.318   (   0.000    0.000    0.000)    0.000
   2.318   (   0.000    0.000    0.000)    0.000
   3.449   (   0.000    0.000    0.000)    0.000
   3.449   (   0.000    0.000    0.000)    0.000
   3.449   (   0.000    0.000    0.000)    0.000
   4.053   (   0.000    0.000    0.000)    0.000
   4.053   (   0.000    0.000    0.000)    0.000
   4.053   (   0.000    0.000    0.000)    0.000
  11.430   (   0.000    0.000    0.000)    0.000
  11.430   (   0.000    0.000    0.000)    0.000
  11.430   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/56) =======================
q-point: ( 0.09  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.396   (  19.155    0.000    0.000)   19.155
   0.396   (  19.155    0.000    0.000)   19.155
   0.870   (  38.994    0.000    0.000)   38.994
   2.304   (  -1.381    0.000    0.000)    1.381
   2.304   (  -1.381    0.000    0.000)    1.381
   2.615   (  16.226    0.000    0.000)   16.226
   3.474   (   2.445    0.000    0.000)    2.445
   3.474   (   2.445    0.000    0.000)    2.445
   3.481   (   3.013    0.000    0.000)    3.013
   4.055   (   0.172    0.000    0.000)    0.172
   4.055   (   0.172    0.000    0.000)    0.172
   6.419   (   1.823    0.000    0.000)    1.823
  11.428   (  -0.144    0.000    0.000)    0.144
  11.428   (  -0.144    0.000    0.000)    0.144
  12.869   (  -5.769    0.000    0.000)    5.769
======================= Grid point 2 (3/56) =======================
q-point: ( 0.18  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.773   (  17.776    0.000    0.000)   17.776
   0.773   (  17.776    0.000    0.000)   17.776
   1.537   (  26.497    0.000    0.000)   26.497
   2.260   (  -2.949    0.000    0.000)    2.949
   2.260   (  -2.949    0.000    0.000)    2.949
   3.086   (  28.552    0.000    0.000)   28.552
   3.544   (   4.279    0.000    0.000)    4.279
   3.544   (   4.279    0.000    0.000)    4.279
   3.564   (   4.939    0.000    0.000)    4.939
   4.059   (   0.164    0.000    0.000)    0.164
   4.059   (   0.164    0.000    0.000)    0.164
   6.480   (   4.291    0.000    0.000)    4.291
  11.424   (  -0.240    0.000    0.000)    0.240
  11.424   (  -0.240    0.000    0.000)    0.240
  12.701   ( -10.512    0.000    0.000)   10.512
======================= Grid point 3 (4/56) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.114   (  15.674    0.000    0.000)   15.674
   1.114   (  15.674    0.000    0.000)   15.674
   1.972   (  17.068    0.000    0.000)   17.068
   2.183   (  -4.648    0.000    0.000)    4.648
   2.183   (  -4.648    0.000    0.000)    4.648
   3.641   (   4.984    0.000    0.000)    4.984
   3.641   (   4.984    0.000    0.000)    4.984
   3.670   (   5.204    0.000    0.000)    5.204
   3.684   (  28.009    0.000    0.000)   28.009
   4.060   (  -0.120    0.000    0.000)    0.120
   4.060   (  -0.120    0.000    0.000)    0.120
   6.599   (   7.407    0.000    0.000)    7.407
  11.419   (  -0.258    0.000    0.000)    0.258
  11.419   (  -0.258    0.000    0.000)    0.258
  12.457   ( -12.897    0.000    0.000)   12.897
======================= Grid point 4 (5/56) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.406   (  12.916    0.000    0.000)   12.916
   1.406   (  12.916    0.000    0.000)   12.916
   2.072   (  -6.172    0.000    0.000)    6.172
   2.072   (  -6.172    0.000    0.000)    6.172
   2.249   (  10.326    0.000    0.000)   10.326
   3.737   (   4.252    0.000    0.000)    4.252
   3.737   (   4.252    0.000    0.000)    4.252
   3.765   (   3.863    0.000    0.000)    3.863
   4.053   (  -0.534    0.000    0.000)    0.534
   4.053   (  -0.534    0.000    0.000)    0.534
   4.159   (  17.832    0.000    0.000)   17.832
   6.769   (   8.691    0.000    0.000)    8.691
  11.414   (  -0.195    0.000    0.000)    0.195
  11.414   (  -0.195    0.000    0.000)    0.195
  12.204   ( -11.240    0.000    0.000)   11.240
======================= Grid point 5 (6/56) =======================
q-point: ( 0.45  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.625   (   7.636    0.000    0.000)    7.636
   1.625   (   7.636    0.000    0.000)    7.636
   1.945   (  -5.372    0.000    0.000)    5.372
   1.945   (  -5.372    0.000    0.000)    5.372
   2.391   (   3.545    0.000    0.000)    3.545
   3.802   (   1.811    0.000    0.000)    1.811
   3.802   (   1.811    0.000    0.000)    1.811
   3.820   (   1.417    0.000    0.000)    1.417
   4.042   (  -0.426    0.000    0.000)    0.426
   4.042   (  -0.426    0.000    0.000)    0.426
   4.398   (   5.748    0.000    0.000)    5.748
   6.911   (   4.147    0.000    0.000)    4.147
  11.412   (  -0.072    0.000    0.000)    0.072
  11.412   (  -0.072    0.000    0.000)    0.072
  12.035   (  -4.598    0.000    0.000)    4.598
======================= Grid point 12 (7/56) =======================
q-point: ( 0.09  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.557   (  13.328   13.328    0.000)   18.849
   0.707   (  15.632   15.632    0.000)   22.106
   1.134   (  24.896   24.896    0.000)   35.208
   2.292   (  -1.221   -1.221    0.000)    1.726
   2.484   (   7.787    7.787    0.000)   11.013
   2.557   (   5.691    5.691    0.000)    8.048
   3.422   (  -1.552   -1.552    0.000)    2.194
   3.499   (   2.380    2.380    0.000)    3.366
   3.578   (   5.982    5.982    0.000)    8.460
   4.057   (   0.159    0.159    0.000)    0.225
   4.095   (   2.303    2.303    0.000)    3.257
   6.433   (   1.318    1.318    0.000)    1.863
  11.254   (  -8.591   -8.591    0.000)   12.149
  11.427   (  -0.153   -0.153    0.000)    0.216
  12.958   (   1.374    1.374    0.000)    1.943
======================= Grid point 13 (8/56) =======================
q-point: ( 0.18  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.861   (  15.547    8.028    0.000)   17.498
   0.984   (  13.176   14.712    0.000)   19.750
   1.644   (  21.760    9.902    0.000)   23.907
   2.253   (  -2.632   -0.741    0.000)    2.734
   2.438   (  -3.148   17.420    0.000)   17.702
   3.017   (  28.541   -5.163    0.000)   29.005
   3.415   (   0.895   -6.071    0.000)    6.136
   3.567   (   4.172    2.206    0.000)    4.719
   3.713   (   7.041    8.327    0.000)   10.905
   4.060   (   0.145    0.127    0.000)    0.192
   4.137   (   2.149    6.112    0.000)    6.479
   6.477   (   3.261   -0.183    0.000)    3.266
  11.130   (  -4.675  -21.068    0.000)   21.580
  11.422   (  -0.255   -0.177    0.000)    0.311
  12.891   (  -6.738   11.597    0.000)   13.412
======================= Grid point 14 (9/56) =======================
q-point: ( 0.27  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.168   (  14.403    5.059    0.000)   15.266
   1.251   (  12.777   10.296    0.000)   16.409
   1.985   (  11.583    1.411    0.000)   11.668
   2.184   (  -4.173    0.053    0.000)    4.173
   2.390   (  -1.102   19.545    0.000)   19.576
   3.414   (   3.116  -12.307    0.000)   12.695
   3.625   (  24.174   -3.889    0.000)   24.485
   3.661   (   4.874    1.983    0.000)    5.262
   3.856   (   6.902    9.563    0.000)   11.793
   4.061   (  -0.141    0.090    0.000)    0.168
   4.193   (   3.753   10.091    0.000)   10.766
   6.572   (   6.141   -2.471    0.000)    6.620
  11.040   (  -4.333  -28.561    0.000)   28.888
  11.417   (  -0.275   -0.208    0.000)    0.345
  12.721   (  -9.265   17.700    0.000)   19.978
======================= Grid point 15 (10/56) =======================
q-point: ( 0.36  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.437   (  11.876    2.918    0.000)   12.229
   1.492   (  10.626    6.487    0.000)   12.449
   2.085   (  -5.472    1.216    0.000)    5.606
   2.097   (  -0.051    0.496    0.000)    0.499
   2.423   (   4.086   18.333    0.000)   18.783
   3.454   (   1.533  -15.472    0.000)   15.547
   3.756   (   4.187    1.797    0.000)    4.556
   3.957   (   6.633    4.569    0.000)    8.055
   4.011   (   8.048    2.219    0.000)    8.348
   4.054   (  -0.562    0.049    0.000)    0.564
   4.321   (   8.959   11.654    0.000)   14.700
   6.719   (   7.702   -4.762    0.000)    9.055
  10.954   (  -4.004  -33.706    0.000)   33.943
  11.412   (  -0.208   -0.238    0.000)    0.316
  12.542   (  -7.710   21.971    0.000)   23.284
======================= Grid point 16 (11/56) =======================
q-point: ( 0.45  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.634   (   6.487    0.789    0.000)    6.535
   1.662   (   5.163    2.866    0.000)    5.905
   1.976   (  -4.319    2.960    0.000)    5.236
   2.054   (  -2.388    6.378    0.000)    6.811
   2.498   (   2.281   13.386    0.000)   13.579
   3.476   (   0.586  -17.360    0.000)   17.370
   3.820   (   1.806    1.729    0.000)    2.500
   4.039   (   2.068    6.823    0.000)    7.129
   4.042   (  -0.455   -0.012    0.000)    0.455
   4.094   (   1.355    5.144    0.000)    5.320
   4.476   (   4.431    7.945    0.000)    9.097
   6.847   (   3.800   -6.096    0.000)    7.184
  10.891   (  -1.797  -36.507    0.000)   36.551
  11.409   (  -0.077   -0.256    0.000)    0.268
  12.429   (  -2.971   24.325    0.000)   24.506
======================= Grid point 24 (12/56) =======================
q-point: ( 0.18  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.066   (  11.451   11.451    0.000)   16.194
   1.241   (  10.721   10.721    0.000)   15.162
   1.866   (  10.019   10.019    0.000)   14.169
   2.229   (  -1.669   -1.669    0.000)    2.360
   2.872   (   9.449    9.449    0.000)   13.363
   2.985   (  15.470   15.470    0.000)   21.878
   3.316   (  -3.020   -3.020    0.000)    4.271
   3.630   (   3.870    3.870    0.000)    5.473
   3.888   (   8.708    8.708    0.000)   12.314
   4.063   (   0.102    0.102    0.000)    0.144
   4.265   (   6.253    6.253    0.000)    8.843
   6.478   (   0.513    0.513    0.000)    0.726
  10.742   ( -16.368  -16.368    0.000)   23.148
  11.418   (  -0.296   -0.296    0.000)    0.419
  13.032   (   2.173    2.173    0.000)    3.074
======================= Grid point 25 (13/56) =======================
q-point: ( 0.27  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.302   (  11.324    7.646    0.000)   13.663
   1.443   (   9.373    8.287    0.000)   12.511
   2.021   (   5.220    1.796    0.000)    5.520
   2.184   (  -2.703   -0.031    0.000)    2.703
   2.880   (   0.198   22.877    0.000)   22.878
   3.207   (  -1.915   -4.816    0.000)    5.183
   3.561   (  23.765   -1.996    0.000)   23.849
   3.718   (   4.545    3.452    0.000)    5.707
   4.058   (   7.802   10.384    0.000)   12.988
   4.063   (  -0.176    0.069    0.000)    0.189
   4.401   (   7.483    9.934    0.000)   12.437
   6.509   (   2.723   -3.403    0.000)    4.358
  10.456   ( -12.114  -27.236    0.000)   29.809
  11.411   (  -0.320   -0.350    0.000)    0.475
  13.009   (  -3.550    9.749    0.000)   10.375
======================= Grid point 26 (14/56) =======================
q-point: ( 0.36  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.513   (   9.037    4.128    0.000)    9.935
   1.617   (   7.432    5.653    0.000)    9.337
   2.085   (   1.390   -1.322    0.000)    1.918
   2.121   (  -3.374    2.331    0.000)    4.101
   2.903   (   1.718   23.181    0.000)   23.244
   3.160   (  -2.100  -10.575    0.000)   10.781
   3.807   (   3.969    3.063    0.000)    5.013
   3.932   (  12.002   -2.204    0.000)   12.202
   4.055   (  -0.613    0.048    0.000)    0.615
   4.198   (   5.654   11.871    0.000)   13.148
   4.568   (   8.222   12.041    0.000)   14.580
   6.590   (   4.960   -7.412    0.000)    8.918
  10.244   (  -8.591  -34.456    0.000)   35.511
  11.405   (  -0.243   -0.401    0.000)    0.469
  12.921   (  -4.399   14.467    0.000)   15.121
======================= Grid point 27 (15/56) =======================
q-point: ( 0.45  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.652   (   4.033    0.868    0.000)    4.125
   1.729   (   3.056    3.666    0.000)    4.773
   2.060   (  -2.079    5.108    0.000)    5.515
   2.094   (  -0.076   -1.341    0.000)    1.343
   2.935   (   0.995   25.028    0.000)   25.048
   3.128   (  -0.912  -15.928    0.000)   15.954
   3.868   (   1.784    2.893    0.000)    3.399
   4.041   (  -0.530   -0.020    0.000)    0.530
   4.077   (   3.127   -1.502    0.000)    3.469
   4.278   (   2.069   12.500    0.000)   12.670
   4.693   (   3.351   12.189    0.000)   12.641
   6.680   (   2.837   -9.811    0.000)   10.213
  10.119   (  -3.328  -38.015    0.000)   38.161
  11.402   (  -0.090   -0.432    0.000)    0.441
  12.853   (  -1.885   16.682    0.000)   16.788
======================= Grid point 36 (16/56) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.460   (   7.464    7.464    0.000)   10.556
   1.596   (   6.960    6.960    0.000)    9.843
   2.048   (   0.742    0.742    0.000)    1.050
   2.182   (  -0.211   -0.211    0.000)    0.299
   2.954   (  -5.477   -5.477    0.000)    7.746
   3.401   (   6.953    6.953    0.000)    9.833
   3.622   (  12.635   12.635    0.000)   17.868
   3.795   (   4.029    4.029    0.000)    5.698
   4.062   (  -0.162   -0.162    0.000)    0.229
   4.279   (  11.150   11.150    0.000)   15.768
   4.599   (   9.815    9.815    0.000)   13.881
   6.449   (  -2.036   -2.036    0.000)    2.879
   9.959   ( -21.422  -21.422    0.000)   30.296
  11.403   (  -0.381   -0.381    0.000)    0.538
  13.124   (   2.232    2.232    0.000)    3.156
======================= Grid point 37 (17/56) =======================
q-point: ( 0.36  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.592   (   5.292    3.482    0.000)    6.334
   1.725   (   5.484    4.995    0.000)    7.417
   2.052   (  -0.233   -1.804    0.000)    1.819
   2.176   (  -0.306    2.962    0.000)    2.978
   2.858   (  -3.597  -10.867    0.000)   11.447
   3.426   (  -0.542   18.187    0.000)   18.195
   3.874   (   3.558    3.458    0.000)    4.962
   3.892   (   9.847   -1.470    0.000)    9.956
   4.055   (  -0.585   -0.067    0.000)    0.589
   4.486   (   8.349   16.296    0.000)   18.310
   4.801   (   9.372   10.800    0.000)   14.299
   6.438   (   0.919   -7.087    0.000)    7.146
   9.591   ( -14.513  -29.123    0.000)   32.539
  11.396   (  -0.290   -0.437    0.000)    0.525
  13.129   (  -1.047    6.410    0.000)    6.495
======================= Grid point 38 (18/56) =======================
q-point: ( 0.45  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.667   (   1.911    0.590    0.000)    2.000
   1.807   (   2.221    3.889    0.000)    4.479
   2.045   (  -0.273   -2.880    0.000)    2.893
   2.171   (  -0.155    5.510    0.000)    5.513
   2.809   (  -1.218  -13.830    0.000)   13.883
   3.416   (  -0.288   19.275    0.000)   19.277
   3.931   (   1.722    3.129    0.000)    3.571
   4.015   (   2.713   -4.156    0.000)    4.963
   4.041   (  -0.610   -0.025    0.000)    0.610
   4.601   (   2.857   18.874    0.000)   19.089
   4.940   (   3.640   11.687    0.000)   12.241
   6.469   (   1.370  -10.383    0.000)   10.473
   9.384   (  -5.366  -33.402    0.000)   33.830
  11.392   (  -0.108   -0.471    0.000)    0.484
  13.106   (  -0.804    8.494    0.000)    8.532
======================= Grid point 48 (19/56) =======================
q-point: ( 0.36  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.649   (   2.048    2.048    0.000)    2.896
   1.819   (   4.066    4.066    0.000)    5.750
   2.016   (  -1.637   -1.637    0.000)    2.315
   2.234   (   2.532    2.532    0.000)    3.581
   2.672   (  -7.244   -7.244    0.000)   10.245
   3.664   (   5.492    5.492    0.000)    7.767
   3.889   (   1.398    1.398    0.000)    1.977
   3.941   (   2.954    2.954    0.000)    4.178
   4.051   (  -0.362   -0.362    0.000)    0.512
   4.810   (  13.687   13.687    0.000)   19.356
   5.007   (   9.442    9.442    0.000)   13.353
   6.321   (  -4.059   -4.059    0.000)    5.740
   9.083   ( -20.296  -20.296    0.000)   28.702
  11.388   (  -0.334   -0.334    0.000)    0.473
  13.204   (   1.616    1.616    0.000)    2.285
======================= Grid point 49 (20/56) =======================
q-point: ( 0.45  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.676   (   0.647    0.309    0.000)    0.717
   1.881   (   1.780    3.265    0.000)    3.719
   1.990   (  -0.750   -2.318    0.000)    2.437
   2.272   (   0.994    4.107    0.000)    4.226
   2.570   (  -2.595   -9.457    0.000)    9.806
   3.707   (   0.334    9.750    0.000)    9.756
   3.925   (   1.304   -4.153    0.000)    4.353
   3.989   (   1.555    2.403    0.000)    2.862
   4.040   (  -0.588   -0.050    0.000)    0.590
   4.996   (   4.488   18.157    0.000)   18.703
   5.157   (   4.185    9.343    0.000)   10.237
   6.275   (  -0.740   -8.258    0.000)    8.291
   8.793   (  -7.551  -23.962    0.000)   25.124
  11.384   (  -0.125   -0.361    0.000)    0.382
  13.217   (   0.089    3.029    0.000)    3.030
======================= Grid point 60 (21/56) =======================
q-point: ( 0.45  0.45  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.680   (   0.093    0.093    0.000)    0.131
   1.932   (   1.476    1.476    0.000)    2.088
   1.957   (  -0.920   -0.920    0.000)    1.301
   2.330   (   1.505    1.505    0.000)    2.129
   2.435   (  -3.516   -3.516    0.000)    4.973
   3.829   (   2.298    2.298    0.000)    3.249
   3.866   (  -1.101   -1.101    0.000)    1.557
   4.025   (   1.075    1.075    0.000)    1.520
   4.038   (  -0.201   -0.201    0.000)    0.285
   5.267   (   7.218    7.218    0.000)   10.208
   5.301   (   4.188    4.188    0.000)    5.923
   6.153   (  -3.335   -3.335    0.000)    4.716
   8.448   (  -9.073   -9.073    0.000)   12.831
  11.378   (  -0.135   -0.135    0.000)    0.191
  13.250   (   0.577    0.577    0.000)    0.816
======================= Grid point 133 (22/56) =======================
q-point: ( 0.09  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.804   (  11.878   11.878   11.878)   20.572
   0.804   (  11.878   11.878   11.878)   20.572
   1.334   (  18.676   18.676   18.676)   32.348
   2.471   (   4.785    4.785    4.785)    8.287
   2.471   (   4.785    4.785    4.785)    8.287
   2.504   (   2.737    2.737    2.737)    4.740
   3.374   (  -2.499   -2.499   -2.499)    4.328
   3.586   (   3.708    3.708    3.708)    6.423
   3.586   (   3.708    3.708    3.708)    6.423
   4.102   (   1.995    1.995    1.995)    3.455
   4.102   (   1.995    1.995    1.995)    3.455
   6.442   (   0.878    0.878    0.878)    1.520
  11.252   (  -5.797   -5.797   -5.797)   10.041
  11.252   (  -5.797   -5.797   -5.797)   10.041
  13.045   (   3.676    3.676    3.676)    6.367
======================= Grid point 134 (23/56) =======================
q-point: ( 0.18  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.029   (  12.476    9.694    9.694)   18.536
   1.064   (  11.766    9.507    9.507)   17.867
   1.733   (  17.124    7.830    7.830)   20.392
   2.428   (  -1.803    8.712    8.712)   12.453
   2.431   (  -2.617    8.010    8.010)   11.626
   2.972   (  28.195   -3.506   -3.506)   28.628
   3.342   (  -0.664   -4.997   -4.997)    7.098
   3.635   (   3.705    3.829    3.829)    6.561
   3.721   (   6.720    2.424    2.424)    7.543
   4.124   (   1.407    3.673    3.673)    5.382
   4.168   (   3.359    4.028    4.028)    6.613
   6.472   (   2.310   -0.474   -0.474)    2.405
  11.056   (  -8.126  -10.658  -10.658)   17.124
  11.247   (  -0.299   -8.651   -8.651)   12.239
  13.036   (  -3.970    9.741    9.741)   14.337
======================= Grid point 135 (24/56) =======================
q-point: ( 0.27  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.287   (  12.443    6.725    6.725)   15.661
   1.299   (  11.266    6.768    6.768)   14.783
   1.979   (   6.830   -0.039   -0.039)    6.831
   2.365   (  -3.911    8.543    8.543)   12.699
   2.432   (   2.549   12.067   12.067)   17.254
   3.343   (   0.678   -7.558   -7.558)   10.710
   3.478   (  16.454   -5.831   -5.831)   18.405
   3.777   (  12.633    4.740    4.740)   14.301
   3.848   (   5.534    1.177    1.177)    5.779
   4.158   (   1.891    5.579    5.579)    8.113
   4.249   (   4.910    6.464    6.464)   10.377
   6.544   (   4.932   -2.619   -2.619)    6.168
  10.904   (  -7.050  -16.664  -16.664)   24.598
  11.241   (  -0.326   -8.712   -8.712)   12.325
  12.917   (  -7.060   14.335   14.335)   21.466
======================= Grid point 136 (25/56) =======================
q-point: ( 0.36  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.512   (   9.294    3.913    3.913)   10.817
   1.520   (  10.117    4.288    4.288)   11.795
   2.030   (  -0.680   -1.661   -1.661)    2.446
   2.276   (  -4.575    9.618    9.618)   14.351
   2.517   (   4.770   12.550   12.550)   18.378
   3.362   (   1.027   -9.539   -9.539)   13.530
   3.640   (   3.417   -5.452   -5.452)    8.433
   3.939   (   3.326   -0.758   -0.758)    3.495
   4.063   (  10.653    4.405    4.405)   12.341
   4.197   (   1.828    7.827    7.827)   11.218
   4.382   (   8.018    6.885    6.885)   12.614
   6.668   (   6.751   -4.811   -4.811)    9.585
  10.770   (  -5.896  -20.658  -20.658)   29.804
  11.235   (  -0.249   -8.771   -8.771)   12.406
  12.776   (  -6.141   17.351   17.351)   25.295
======================= Grid point 137 (26/56) =======================
q-point: ( 0.45  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.658   (   4.348    1.238    1.238)    4.687
   1.677   (   4.563    2.020    2.020)    5.384
   1.995   (  -1.719    0.561    0.561)    1.894
   2.195   (  -2.662   11.284   11.284)   16.178
   2.591   (   2.042   11.256   11.256)   16.049
   3.379   (   0.464  -10.594  -10.594)   14.989
   3.680   (   0.879   -5.680   -5.680)    8.080
   3.983   (   1.055   -2.024   -2.024)    3.050
   4.186   (   2.525    6.626    6.626)    9.705
   4.226   (   0.784    9.200    9.200)   13.034
   4.518   (   3.863    5.008    5.008)    8.068
   6.783   (   3.470   -6.086   -6.086)    9.281
  10.680   (  -2.486  -22.701  -22.701)   32.200
  11.231   (  -0.093   -8.807   -8.807)   12.455
  12.686   (  -2.396   18.880   18.880)   26.808
======================= Grid point 145 (27/56) =======================
q-point: ( 0.18  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.206   (   9.391    9.391   10.456)   16.903
   1.285   (   9.931    9.931    4.194)   14.657
   1.894   (   6.598    6.598    2.326)    9.616
   2.420   (  -0.226   -0.226   18.140)   18.143
   2.857   (   9.351    9.351   -1.434)   13.301
   2.969   (  15.957   15.957   -0.671)   22.576
   3.250   (  -3.565   -3.565   -5.203)    7.245
   3.698   (   3.643    3.643    5.234)    7.344
   3.817   (   5.650    5.650   -3.473)    8.712
   4.184   (   3.097    3.097    8.235)    9.327
   4.285   (   6.351    6.351    2.028)    9.208
   6.463   (  -0.209   -0.209   -1.479)    1.508
  10.739   ( -16.385  -16.385   -0.251)   23.173
  11.175   (  -1.355   -1.355  -20.103)   20.194
  13.170   (   2.802    2.802   10.205)   10.948
======================= Grid point 146 (28/56) =======================
q-point: ( 0.27  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.402   (   9.472    6.352    7.664)   13.741
   1.473   (   8.699    7.800    2.836)   12.023
   1.982   (   2.268    0.179   -1.802)    2.902
   2.415   (  -0.499    1.665   19.223)   19.302
   2.885   (   1.315   21.167    0.820)   21.224
   3.169   (  -2.349   -5.709   -3.741)    7.219
   3.427   (  15.169    0.599   -7.821)   17.077
   3.827   (   9.137    1.919    5.748)   10.964
   3.935   (   7.047    5.796   -2.703)    9.516
   4.256   (   3.847    5.043    9.782)   11.658
   4.422   (   7.324    9.158    2.152)   11.923
   6.475   (   1.709   -3.843   -3.211)    5.291
  10.436   ( -13.314  -25.157   -2.185)   28.546
  11.155   (  -0.804   -1.928  -21.356)   21.458
  13.166   (  -2.470    8.974   12.199)   15.345
======================= Grid point 147 (29/56) =======================
q-point: ( 0.36  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.574   (   7.101    3.335    4.776)    9.184
   1.635   (   6.896    5.420    1.665)    8.928
   2.002   (   0.081   -1.178   -2.680)    2.929
   2.393   (  -1.582    3.712   19.134)   19.555
   2.931   (   2.790   20.211    3.441)   20.691
   3.125   (  -1.815   -9.510   -3.636)   10.342
   3.589   (   3.451   -0.115  -12.343)   12.817
   3.944   (   2.717    1.459   -2.995)    4.299
   4.156   (  11.078    4.333    8.753)   14.769
   4.329   (   2.995    6.758    9.741)   12.228
   4.581   (   7.683   11.059    1.429)   13.541
   6.537   (   4.134   -7.632   -5.087)   10.061
  10.201   (  -9.603  -31.335   -4.747)   33.115
  11.141   (  -0.537   -2.263  -21.684)   21.808
  13.098   (  -3.558   12.919   13.946)   19.341
======================= Grid point 148 (30/56) =======================
q-point: ( 0.45  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.677   (   2.692    0.829    2.053)    3.485
   1.738   (   2.818    3.562    0.929)    4.636
   1.998   (  -0.219   -0.206   -1.657)    1.684
   2.363   (  -1.009    5.545   19.067)   19.883
   2.979   (   1.553   22.077    5.399)   22.781
   3.097   (  -0.892  -14.294   -3.378)   14.715
   3.627   (   0.731   -0.943  -14.560)   14.609
   3.976   (   0.732    1.750   -4.036)    4.460
   4.309   (   3.959    2.430   12.696)   13.519
   4.368   (   0.734    8.999    7.866)   11.975
   4.697   (   3.067   11.416    0.522)   11.833
   6.615   (   2.559   -9.907   -6.180)   11.954
  10.061   (  -3.713  -34.430   -6.286)   35.196
  11.134   (  -0.195   -2.442  -21.785)   21.922
  13.042   (  -1.569   14.766   14.911)   21.044
======================= Grid point 157 (31/56) =======================
q-point: ( 0.27  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.533   (   6.198    6.198    5.760)   10.488
   1.618   (   6.564    6.564    2.028)    9.501
   1.983   (  -0.099   -0.099   -2.901)    2.904
   2.447   (   1.260    1.260   20.766)   20.842
   2.940   (  -5.699   -5.699   -1.406)    8.181
   3.271   (   4.854    4.854  -10.012)   12.139
   3.543   (   8.353    8.353   -5.162)   12.891
   3.905   (   6.068    6.068    7.983)   11.720
   4.067   (   7.765    7.765    0.492)   10.992
   4.375   (   6.732    6.732    7.816)   12.318
   4.608   (   9.309    9.309    0.956)   13.200
   6.404   (  -2.746   -2.746   -4.376)    5.851
   9.956   ( -21.406  -21.406   -0.266)   30.274
  11.126   (  -1.084   -1.084  -24.133)   24.182
  13.277   (   2.410    2.410   12.198)   12.665
======================= Grid point 158 (32/56) =======================
q-point: ( 0.36  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.640   (   3.993    2.867    3.852)    6.245
   1.738   (   5.121    4.814    1.281)    7.144
   1.982   (   0.145   -0.736   -2.815)    2.913
   2.463   (   0.351    3.256   21.819)   22.063
   2.842   (  -3.687  -10.827   -1.652)   11.556
   3.328   (   2.155    9.877   -8.102)   12.955
   3.624   (   0.547    2.607  -11.915)   12.209
   4.016   (   3.912    5.742    5.316)    8.748
   4.232   (   6.711    3.156   11.158)   13.398
   4.522   (   6.972   12.779    3.618)   15.000
   4.801   (   9.001   10.352    0.022)   13.718
   6.378   (   0.301   -7.514   -5.755)    9.470
   9.583   ( -14.902  -28.346   -0.837)   32.035
  11.108   (  -0.704   -1.224  -25.190)   25.230
  13.288   (  -0.707    5.990   12.797)   14.147
======================= Grid point 159 (33/56) =======================
q-point: ( 0.45  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.692   (   1.229    0.671    2.278)    2.674
   1.814   (   2.038    3.816    0.717)    4.385
   1.987   (   0.226   -0.824   -1.989)    2.165
   2.465   (  -0.016    4.868   22.257)   22.783
   2.792   (  -1.259  -13.529   -1.795)   13.705
   3.366   (   1.233   12.090   -3.330)   12.601
   3.614   (  -0.700   -0.573  -17.241)   17.264
   4.066   (   1.179    6.696    2.989)    7.427
   4.317   (   1.889   -0.564   15.011)   15.140
   4.626   (   2.746   16.211    2.431)   16.620
   4.934   (   3.472   11.367   -0.552)   11.898
   6.401   (   1.178  -10.644   -6.554)   12.555
   9.370   (  -5.557  -32.270   -1.525)   32.781
  11.098   (  -0.254   -1.292  -25.635)   25.669
  13.270   (  -0.653    7.780   13.248)   15.377
======================= Grid point 169 (34/56) =======================
q-point: ( 0.36  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.685   (   1.513    1.513    3.083)    3.753
   1.828   (   3.867    3.867    0.908)    5.543
   1.972   (  -0.285   -0.285   -1.973)    2.014
   2.522   (   2.488    2.488   23.513)   23.774
   2.656   (  -7.164   -7.164   -1.470)   10.238
   3.439   (   3.274    3.274  -15.448)   16.127
   3.640   (  -1.182   -1.182  -13.731)   13.832
   4.133   (   4.722    4.722   12.324)   14.017
   4.283   (   1.723    1.723   13.035)   13.260
   4.809   (  13.063   13.063   -0.151)   18.474
   5.001   (   9.220    9.220   -0.644)   13.055
   6.252   (  -4.477   -4.477   -6.649)    9.181
   9.082   ( -20.222  -20.222   -0.087)   28.598
  11.088   (  -0.774   -0.774  -26.578)   26.601
  13.360   (   1.588    1.588   12.582)   12.781
======================= Grid point 170 (35/56) =======================
q-point: ( 0.45  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.703   (   0.398    0.393    2.513)    2.575
   1.887   (   1.647    3.180    0.553)    3.624
   1.969   (  -0.091   -0.891   -0.902)    1.271
   2.556   (  -1.799   -6.589    3.643)    7.741
   2.559   (   0.180    1.352   19.322)   19.370
   3.499   (   2.537    2.862  -14.607)   15.099
   3.596   (  -2.451   -1.307  -17.457)   17.676
   4.198   (   1.662    5.282   11.728)   12.969
   4.303   (   0.360   -0.663   14.740)   14.759
   4.990   (   4.466   17.547   -0.595)   18.117
   5.147   (   4.051    9.219   -1.015)   10.121
   6.201   (  -0.856   -8.606   -7.181)   11.241
   8.792   (  -7.588  -23.633   -0.122)   24.821
  11.077   (  -0.276   -0.802  -27.165)   27.179
  13.373   (   0.116    2.809   12.627)   12.936
======================= Grid point 181 (36/56) =======================
q-point: ( 0.45  0.45  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.708   (   0.122    0.122    2.656)    2.662
   1.936   (   1.412    1.412    0.400)    2.036
   1.955   (  -0.405   -0.405    0.053)    0.575
   2.424   (  -3.453   -3.453   -1.103)    5.006
   2.619   (   1.620    1.620   25.873)   25.975
   3.535   (   1.277    1.277  -18.807)   18.894
   3.565   (  -1.560   -1.560  -18.448)   18.580
   4.270   (   1.824    1.824   14.046)   14.281
   4.293   (  -0.327   -0.327   14.221)   14.228
   5.256   (   7.204    7.204   -1.113)   10.249
   5.290   (   4.163    4.163   -1.056)    5.982
   6.072   (  -3.533   -3.533   -7.784)    9.249
   8.450   (  -9.014   -9.014    0.217)   12.750
  11.066   (  -0.285   -0.285  -27.806)   27.809
  13.403   (   0.538    0.538   12.457)   12.480
======================= Grid point 266 (37/56) =======================
q-point: ( 0.18  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.401   (   7.630    7.630    7.630)   13.216
   1.401   (   7.630    7.630    7.630)   13.216
   1.938   (   1.595    1.595    1.595)    2.763
   2.818   (   5.249    5.249    5.249)    9.092
   2.818   (   5.249    5.249    5.249)    9.092
   3.062   (  14.994   14.994   14.994)   25.970
   3.142   (  -5.039   -5.039   -5.039)    8.727
   3.801   (   3.610    3.610    3.610)    6.253
   3.801   (   3.610    3.610    3.610)    6.253
   4.346   (   5.698    5.698    5.698)    9.869
   4.346   (   5.698    5.698    5.698)    9.869
   6.423   (  -2.202   -2.202   -2.202)    3.815
  10.732   ( -11.097  -11.097  -11.097)   19.220
  10.732   ( -11.097  -11.097  -11.097)   19.220
  13.341   (   5.579    5.579    5.579)    9.663
======================= Grid point 267 (38/56) =======================
q-point: ( 0.27  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.550   (   7.397    5.786    5.786)   11.030
   1.554   (   6.679    5.585    5.585)   10.344
   1.956   (   0.337   -1.029   -1.029)    1.493
   2.768   (  -2.459    8.882    8.882)   12.799
   2.963   (   5.958    7.033    7.033)   11.594
   3.052   (  -3.707   -6.719   -6.719)   10.199
   3.384   (   9.227    6.814    6.814)   13.342
   3.874   (   2.624    2.322    2.322)    4.204
   3.987   (  15.821    4.763    4.763)   17.196
   4.432   (   4.390    6.927    6.927)   10.735
   4.496   (   7.402    5.562    5.562)   10.801
   6.387   (  -0.983   -5.177   -5.177)    7.387
  10.322   ( -18.283  -11.670  -11.670)   24.630
  10.723   (  -0.470  -16.495  -16.495)   23.332
  13.395   (   0.188    8.994    8.994)   12.721
======================= Grid point 268 (39/56) =======================
q-point: ( 0.36  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.667   (   4.192    3.656    3.656)    6.656
   1.687   (   5.783    3.808    3.808)    7.902
   1.962   (   0.530   -1.452   -1.452)    2.121
   2.721   (  -2.059    9.450    9.450)   13.521
   2.993   (  -2.140   -8.429   -8.429)   12.111
   3.094   (   6.705    8.729    8.729)   14.048
   3.446   (  -0.925    2.212    2.212)    3.262
   3.912   (   1.243    1.479    1.479)    2.434
   4.316   (  13.473    6.697    6.697)   16.469
   4.513   (   3.327    7.053    7.053)   10.514
   4.645   (   6.735    5.721    5.721)   10.527
   6.396   (   1.883   -8.390   -8.390)   12.014
  10.000   ( -12.976  -16.567  -16.567)   26.783
  10.714   (  -0.364  -16.554  -16.554)   23.414
  13.372   (  -1.867   11.524   11.524)   16.404
======================= Grid point 269 (40/56) =======================
q-point: ( 0.45  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.720   (   1.192    1.908    1.908)    2.950
   1.772   (   2.233    2.547    2.547)    4.238
   1.976   (   0.545   -0.756   -0.756)    1.200
   2.690   (  -0.876   10.212   10.212)   14.469
   2.965   (  -0.664   -9.567   -9.567)   13.546
   3.218   (   4.714   14.693   14.693)   21.308
   3.397   (  -2.938   -4.916   -4.916)    7.548
   3.928   (   0.358    0.972    0.972)    1.420
   4.498   (   4.534    7.823    7.823)   11.956
   4.561   (   1.218    7.067    7.067)   10.068
   4.741   (   2.285    5.388    5.388)    7.955
   6.442   (   1.804  -10.310  -10.310)   14.692
   9.813   (  -4.907  -19.263  -19.263)   27.680
  10.709   (  -0.137  -16.591  -16.591)   23.463
  13.339   (  -1.000   12.790   12.790)   18.115
======================= Grid point 278 (41/56) =======================
q-point: ( 0.27  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.651   (   4.714    4.714    5.505)    8.645
   1.675   (   5.527    5.527    3.497)    8.564
   1.939   (  -0.514   -0.514   -1.667)    1.819
   2.843   (   1.183    1.183   12.451)   12.563
   2.892   (  -6.325   -6.325   -3.566)    9.630
   3.068   (   2.427    2.427   -8.918)    9.556
   3.577   (   4.812    4.812   13.041)   14.710
   4.054   (   8.420    8.420    6.127)   13.392
   4.080   (  12.182   12.182    0.778)   17.246
   4.543   (   4.651    4.651    8.054)   10.398
   4.641   (   7.881    7.881    2.241)   11.369
   6.281   (  -4.706   -4.706   -7.335)    9.905
   9.949   ( -21.362  -21.362   -0.452)   30.213
  10.548   (  -3.578   -3.578  -30.634)   31.049
  13.517   (   3.283    3.283   10.048)   11.069
======================= Grid point 279 (42/56) =======================
q-point: ( 0.36  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.726   (   2.365    2.405    4.417)    5.558
   1.775   (   4.165    4.314    2.251)    6.406
   1.942   (   0.925   -0.449   -1.399)    1.736
   2.783   (  -3.814   -8.848   -3.309)   10.187
   2.851   (  -0.654    3.456    8.711)    9.395
   3.138   (   4.643   -1.398  -10.271)   11.358
   3.580   (  -2.592    2.749   13.386)   13.909
   4.104   (   1.294   16.835    3.090)   17.165
   4.426   (  13.650    4.388    7.671)   16.261
   4.629   (   3.825    6.060    6.501)    9.676
   4.804   (   7.737    8.803    0.482)   11.730
   6.216   (  -1.394   -8.764   -9.762)   13.192
   9.552   ( -16.516  -24.888   -2.581)   29.981
  10.504   (  -1.259   -5.645  -31.645)   32.169
  13.546   (   0.069    5.690   11.153)   12.521
======================= Grid point 280 (43/56) =======================
q-point: ( 0.45  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.751   (   0.407    1.204    3.280)    3.518
   1.836   (   1.609    3.473    1.430)    4.087
   1.965   (   0.853   -0.335   -0.498)    1.043
   2.732   (  -1.226  -11.770   -4.074)   12.515
   2.836   (  -0.540    5.810    7.533)    9.529
   3.233   (   3.624   -3.700  -14.299)   15.208
   3.514   (  -2.816    2.647   17.766)   18.181
   4.120   (   0.345   17.117    2.013)   17.238
   4.592   (   3.333    0.021   11.848)   12.308
   4.696   (   2.286   10.374    4.322)   11.469
   4.918   (   2.849   10.135   -0.909)   10.567
   6.216   (   0.677  -11.426  -11.139)   15.971
   9.313   (  -6.290  -27.479   -4.622)   28.566
  10.488   (  -0.372   -6.364  -31.577)   32.214
  13.539   (  -0.360    6.931   11.788)   13.680
======================= Grid point 290 (44/56) =======================
q-point: ( 0.36  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.761   (   0.983    0.983    4.212)    4.435
   1.854   (   3.327    3.327    1.637)    4.982
   1.946   (   0.695    0.695   -0.856)    1.303
   2.612   (  -6.808   -6.808   -2.834)   10.037
   2.931   (   3.200    3.200    6.407)    7.844
   3.122   (   0.737    0.737  -15.086)   15.122
   3.554   (  -3.632   -3.632   13.818)   14.742
   4.365   (   5.999    5.999   10.106)   13.195
   4.515   (   4.279    4.279    9.492)   11.257
   4.806   (  10.642   10.642    0.060)   15.051
   4.983   (   8.473    8.473   -1.094)   12.033
   6.063   (  -5.641   -5.641  -11.405)   13.918
   9.079   ( -20.022  -20.022   -0.149)   28.316
  10.431   (  -2.177   -2.177  -35.762)   35.894
  13.616   (   1.659    1.659   11.235)   11.477
======================= Grid point 291 (45/56) =======================
q-point: ( 0.45  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.771   (   0.146    0.727    3.918)    3.988
   1.903   (   1.283    2.939    1.010)    3.362
   1.961   (   0.553   -0.085    0.018)    0.560
   2.517   (  -2.416   -8.731   -2.451)    9.385
   2.964   (   0.381    6.786    1.947)    7.070
   3.171   (   2.784   -2.391  -17.592)   17.970
   3.474   (  -3.168   -4.678   17.808)   18.682
   4.417   (   0.769    9.621    9.444)   13.504
   4.572   (   1.311   -1.336   11.745)   11.894
   4.973   (   4.537   15.478   -1.035)   16.162
   5.115   (   3.510    8.898   -2.067)    9.786
   5.998   (  -1.115   -9.584  -12.268)   15.608
   8.787   (  -7.759  -22.505   -0.382)   23.808
  10.404   (  -0.633   -2.569  -36.736)   36.831
  13.631   (   0.180    2.516   11.367)   11.643
======================= Grid point 302 (46/56) =======================
q-point: ( 0.45  0.45  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.781   (   0.240    0.240    4.212)    4.225
   1.948   (   1.229    1.229    0.741)    1.890
   1.961   (   0.018    0.018    0.464)    0.464
   2.393   (  -3.254   -3.254   -1.901)    4.979
   3.091   (   4.127    4.127   -6.640)    8.841
   3.143   (   0.293    0.293  -19.174)   19.178
   3.367   (  -4.634   -4.634   24.233)   25.104
   4.527   (   2.013    2.013   11.128)   11.487
   4.552   (  -0.330   -0.330   11.198)   11.208
   5.222   (   7.238    7.238   -2.237)   10.478
   5.259   (   4.099    4.099   -1.993)    6.130
   5.851   (  -4.130   -4.130  -13.386)   14.605
   8.456   (  -8.856   -8.856    0.364)   12.530
  10.371   (  -0.749   -0.749  -38.203)   38.218
  13.658   (   0.478    0.478   11.233)   11.254
======================= Grid point 399 (47/56) =======================
q-point: ( 0.27  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.753   (   4.132    4.132    4.132)    7.157
   1.753   (   4.132    4.132    4.132)    7.157
   1.913   (  -0.744   -0.744   -0.744)    1.288
   2.788   (  -6.143   -6.143   -6.143)   10.640
   2.935   (  -1.826   -1.826   -1.826)    3.163
   2.935   (  -1.826   -1.826   -1.826)    3.163
   3.978   (  13.699   13.699   13.699)   23.727
   4.145   (   7.764    7.764    7.764)   13.447
   4.145   (   7.764    7.764    7.764)   13.447
   4.697   (   5.140    5.140    5.140)    8.903
   4.697   (   5.140    5.140    5.140)    8.903
   6.124   (  -7.529   -7.529   -7.529)   13.041
   9.939   ( -14.368  -14.368  -14.368)   24.886
   9.939   ( -14.368  -14.368  -14.368)   24.886
  13.666   (   4.622    4.622    4.622)    8.005
======================= Grid point 400 (48/56) =======================
q-point: ( 0.36  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.810   (   1.057    3.326    3.326)    4.821
   1.828   (   3.171    3.043    3.043)    5.346
   1.922   (   1.834   -0.454   -0.454)    1.944
   2.676   (  -4.400   -6.603   -6.603)   10.323
   2.871   (  -2.364   -2.416   -2.416)    4.155
   2.950   (   1.149   -7.336   -7.336)   10.438
   4.030   (  -0.889   17.365   17.365)   24.573
   4.181   (   0.873   11.112   11.112)   15.739
   4.545   (  17.246    4.205    4.205)   18.242
   4.757   (   2.561    4.958    4.958)    7.465
   4.837   (   6.391    3.710    3.710)    8.269
   6.002   (  -3.880  -10.663  -10.663)   15.571
   9.435   ( -20.499  -11.297  -11.297)   25.989
   9.930   (  -0.386  -21.247  -21.247)   30.050
  13.718   (   1.034    5.805    5.805)    8.275
======================= Grid point 401 (49/56) =======================
q-point: ( 0.45  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.814   (  -0.112    2.625    2.625)    3.715
   1.874   (   1.224    2.348    2.348)    3.539
   1.961   (   1.283    0.074    0.074)    1.288
   2.616   (  -1.471   -7.284   -7.284)   10.406
   2.838   (  -0.865   -2.982   -2.982)    4.305
   2.974   (   0.873   -9.934   -9.934)   14.076
   4.008   (  -0.762   17.648   17.648)   24.969
   4.190   (   0.136   11.414   11.414)   16.143
   4.777   (   5.899    5.584    5.584)    9.857
   4.794   (   0.961    4.408    4.408)    6.308
   4.925   (   1.662    3.508    3.508)    5.233
   5.970   (  -0.010  -12.561  -12.561)   17.764
   9.140   (  -7.732  -13.746  -13.746)   20.921
   9.924   (  -0.145  -21.213  -21.213)   30.001
  13.725   (  -0.021    6.457    6.457)    9.132
======================= Grid point 411 (50/56) =======================
q-point: ( 0.36  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.848   (   0.735    0.735    4.225)    4.351
   1.892   (   2.605    2.605    1.956)    4.172
   1.933   (   1.540    1.540   -0.491)    2.233
   2.542   (  -5.818   -5.818   -4.041)    9.166
   2.827   (  -2.096   -2.096   -9.216)    9.681
   2.847   (  -2.557   -2.557  -11.543)   12.096
   4.099   (  -0.615   -0.615   31.104)   31.116
   4.538   (   7.486    7.486    6.921)   12.648
   4.650   (  10.388   10.388    3.154)   15.026
   4.832   (   4.106    4.106    3.277)    6.668
   4.959   (   6.977    6.977   -1.170)    9.937
   5.811   (  -7.305   -7.305  -12.749)   16.409
   9.076   ( -19.759  -19.759   -0.165)   27.943
   9.702   (  -4.323   -4.323  -34.644)   35.179
  13.792   (   1.871    1.871    6.067)    6.619
======================= Grid point 412 (51/56) =======================
q-point: ( 0.45  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.853   (  -0.008    1.272    3.936)    4.136
   1.926   (   0.853    2.528    1.254)    2.948
   1.966   (   1.116    0.314    0.400)    1.227
   2.459   (  -2.155   -7.276   -3.272)    8.264
   2.787   (  -1.135   -1.582  -11.470)   11.634
   2.831   (  -0.016   -4.597  -15.023)   15.711
   4.083   (  -0.595   -0.676   32.718)   32.730
   4.577   (   0.212   14.452    6.211)   15.731
   4.782   (   1.905   -1.958    8.923)    9.332
   4.951   (   5.074   11.094   -1.014)   12.241
   5.064   (   2.138    8.197   -2.825)    8.929
   5.727   (  -1.242  -10.797  -13.653)   17.451
   8.772   (  -8.364  -19.555   -1.378)   21.313
   9.658   (  -0.811   -6.665  -34.967)   35.606
  13.811   (   0.272    2.337    6.248)    6.676
======================= Grid point 423 (52/56) =======================
q-point: ( 0.45  0.45  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.869   (   0.381    0.381    4.278)    4.312
   1.965   (   0.964    0.964    0.917)    1.643
   1.972   (   0.311    0.311    0.641)    0.778
   2.350   (  -2.915   -2.915   -2.147)    4.648
   2.770   (  -0.681   -0.681  -16.417)   16.445
   2.773   (  -0.931   -0.931  -16.794)   16.846
   4.066   (  -0.606   -0.606   34.875)   34.885
   4.723   (   1.999    1.999    8.153)    8.629
   4.749   (  -0.526   -0.526    8.317)    8.351
   5.164   (   7.461    7.461   -3.558)   11.135
   5.211   (   4.035    4.035   -2.715)    6.320
   5.556   (  -5.044   -5.044  -14.735)   16.371
   8.464   (  -8.650   -8.650    0.396)   12.239
   9.584   (  -1.446   -1.446  -37.712)   37.767
  13.836   (   0.437    0.437    6.179)    6.210
======================= Grid point 532 (53/56) =======================
q-point: ( 0.36  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.922   (   0.983    0.983    0.983)    1.702
   1.930   (   1.828    1.828    1.828)    3.167
   1.930   (   1.828    1.828    1.828)    3.167
   2.453   (  -4.390   -4.390   -4.390)    7.604
   2.666   (  -6.143   -6.143   -6.143)   10.640
   2.666   (  -6.143   -6.143   -6.143)   10.640
   4.647   (   7.711    7.711    7.711)   13.356
   4.647   (   7.711    7.711    7.711)   13.356
   4.666   (   8.584    8.584    8.584)   14.868
   4.938   (   2.810    2.810    2.810)    4.867
   4.938   (   2.810    2.810    2.810)    4.867
   5.576   (  -9.386   -9.386   -9.386)   16.257
   9.073   ( -13.052  -13.052  -13.052)   22.607
   9.073   ( -13.052  -13.052  -13.052)   22.607
  13.872   (   2.110    2.110    2.110)    3.654
======================= Grid point 533 (54/56) =======================
q-point: ( 0.45  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.920   (  -0.090    2.366    2.366)    3.347
   1.954   (   0.647    1.509    1.509)    2.231
   1.976   (   1.307    0.507    0.507)    1.491
   2.384   (  -1.936   -4.066   -4.066)    6.068
   2.585   (  -1.917   -8.680   -8.680)   12.424
   2.586   (  -1.828   -7.747   -7.747)   11.108
   4.631   (  -0.887   12.331   12.331)   17.461
   4.696   (   1.155   11.729   11.729)   16.627
   4.927   (  -0.250    2.409    2.409)    3.416
   4.940   (   7.985    2.350    2.350)    8.649
   5.032   (   0.372    1.984    1.984)    2.830
   5.469   (  -1.170  -11.209  -11.209)   15.894
   8.696   (  -9.980   -8.136   -8.136)   15.231
   9.071   (  -0.048  -19.368  -19.368)   27.391
  13.894   (   0.364    2.258    2.258)    3.214
======================= Grid point 544 (55/56) =======================
q-point: ( 0.45  0.45  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.946   (   0.461    0.461    3.093)    3.161
   1.983   (   0.686    0.686    0.783)    1.247
   1.985   (   0.552    0.552    0.545)    0.952
   2.310   (  -2.413   -2.413   -1.765)    3.842
   2.479   (  -2.322   -2.322  -11.645)   12.099
   2.481   (  -2.476   -2.476  -11.488)   12.010
   4.691   (   0.601    0.601   24.967)   24.981
   4.863   (   1.813    1.813    5.785)    6.328
   4.899   (  -1.470   -1.470    6.929)    7.234
   5.073   (   6.914    6.914   -5.375)   11.157
   5.151   (   4.090    4.090   -3.132)    6.577
   5.292   (  -4.949   -4.949  -10.322)   12.472
   8.471   (  -8.460   -8.460    0.302)   11.967
   8.895   (  -2.187   -2.187  -28.712)   28.878
  13.918   (   0.420    0.420    2.235)    2.312
======================= Grid point 665 (56/56) =======================
q-point: ( 0.45  0.45  0.45)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 2.72e-04 0.00e+00 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.988   (   0.780    0.780    0.780)    1.350
   1.994   (   0.429    0.429    0.429)    0.743
   1.994   (   0.429    0.429    0.429)    0.743
   2.279   (  -1.346   -1.346   -1.346)    2.331
   2.320   (  -3.834   -3.834   -3.834)    6.640
   2.320   (  -3.834   -3.834   -3.834)    6.640
   4.958   (   1.801    1.801    1.801)    3.120
   4.958   (   1.801    1.801    1.801)    3.120
   4.985   (   1.077    1.077    1.077)    1.865
   5.089   (   2.329    2.329    2.329)    4.033
   5.089   (   2.329    2.329    2.329)    4.033
   5.151   (  -2.626   -2.626   -2.626)    4.548
   8.475   (  -5.527   -5.527   -5.527)    9.573
   8.475   (  -5.527   -5.527   -5.527)    9.573
  13.942   (   0.416    0.416    0.416)    0.720
=================== End of collection of collisions ===================
----------- Thermal conductivity (W/m-k) with tetrahedron method -----------
#  T(K)        xx         yy         zz         yz         xz         xy        #ipm
    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/19965
   10.0    351.535    351.535    351.535      0.000     -0.000      0.000 3/19965
   20.0     49.130     49.130     49.130     -0.000     -0.000      0.000 3/19965
   30.0     18.443     18.443     18.443     -0.000     -0.000      0.000 3/19965
   40.0     11.174     11.174     11.174     -0.000     -0.000      0.000 3/19965
   50.0      8.294      8.294      8.294     -0.000     -0.000      0.000 3/19965
   60.0      6.757      6.757      6.757     -0.000     -0.000      0.000 3/19965
   70.0      5.777      5.777      5.777     -0.000     -0.000      0.000 3/19965
   80.0      5.078      5.078      5.078     -0.000     -0.000      0.000 3/19965
   90.0      4.545      4.545      4.545     -0.000     -0.000      0.000 3/19965
  100.0      4.120      4.120      4.120     -0.000     -0.000      0.000 3/19965
  110.0      3.771      3.771      3.771     -0.000     -0.000      0.000 3/19965
  120.0      3.478      3.478      3.478     -0.000     -0.000      0.000 3/19965
  130.0      3.228      3.228      3.228      0.000     -0.000      0.000 3/19965
  140.0      3.011      3.011      3.011      0.000     -0.000      0.000 3/19965
  150.0      2.822      2.822      2.822      0.000     -0.000      0.000 3/19965
  160.0      2.655      2.655      2.655      0.000     -0.000      0.000 3/19965
  170.0      2.507      2.507      2.507      0.000     -0.000      0.000 3/19965
  180.0      2.374      2.374      2.374      0.000     -0.000      0.000 3/19965
  190.0      2.255      2.255      2.255      0.000     -0.000      0.000 3/19965
  200.0      2.146      2.146      2.146      0.000     -0.000      0.000 3/19965
  210.0      2.048      2.048      2.048      0.000     -0.000      0.000 3/19965
  220.0      1.959      1.959      1.959      0.000     -0.000      0.000 3/19965
  230.0      1.876      1.876      1.876      0.000     -0.000      0.000 3/19965
  240.0      1.801      1.801      1.801      0.000     -0.000      0.000 3/19965
  250.0      1.731      1.731      1.731      0.000     -0.000      0.000 3/19965
  260.0      1.666      1.666      1.666      0.000     -0.000      0.000 3/19965
  270.0      1.606      1.606      1.606      0.000     -0.000      0.000 3/19965
  280.0      1.550      1.550      1.550      0.000     -0.000      0.000 3/19965
  290.0      1.498      1.498      1.498      0.000     -0.000      0.000 3/19965
  300.0      1.449      1.449      1.449      0.000     -0.000      0.000 3/19965
  310.0      1.404      1.404      1.404      0.000     -0.000      0.000 3/19965
  320.0      1.361      1.361      1.361      0.000     -0.000      0.000 3/19965
  330.0      1.320      1.320      1.320      0.000     -0.000      0.000 3/19965
  340.0      1.282      1.282      1.282      0.000     -0.000      0.000 3/19965
  350.0      1.246      1.246      1.246      0.000     -0.000      0.000 3/19965
  360.0      1.212      1.212      1.212      0.000     -0.000      0.000 3/19965
  370.0      1.180      1.180      1.180      0.000     -0.000      0.000 3/19965
  380.0      1.150      1.150      1.150      0.000     -0.000      0.000 3/19965
  390.0      1.120      1.120      1.120      0.000     -0.000      0.000 3/19965
  400.0      1.093      1.093      1.093      0.000     -0.000      0.000 3/19965
  410.0      1.067      1.067      1.067      0.000     -0.000      0.000 3/19965
  420.0      1.042      1.042      1.042      0.000     -0.000      0.000 3/19965
  430.0      1.018      1.018      1.018      0.000     -0.000      0.000 3/19965
  440.0      0.995      0.995      0.995      0.000     -0.000      0.000 3/19965
  450.0      0.973      0.973      0.973      0.000     -0.000      0.000 3/19965
  460.0      0.952      0.952      0.952      0.000     -0.000      0.000 3/19965
  470.0      0.932      0.932      0.932      0.000     -0.000      0.000 3/19965
  480.0      0.913      0.913      0.913      0.000     -0.000      0.000 3/19965
  490.0      0.894      0.894      0.894      0.000     -0.000      0.000 3/19965
  500.0      0.877      0.877      0.877      0.000     -0.000      0.000 3/19965
  510.0      0.860      0.860      0.860      0.000     -0.000      0.000 3/19965
  520.0      0.843      0.843      0.843      0.000     -0.000      0.000 3/19965
  530.0      0.828      0.828      0.828      0.000     -0.000      0.000 3/19965
  540.0      0.812      0.812      0.812      0.000     -0.000      0.000 3/19965
  550.0      0.798      0.798      0.798      0.000     -0.000      0.000 3/19965
  560.0      0.784      0.784      0.784      0.000     -0.000      0.000 3/19965
  570.0      0.770      0.770      0.770      0.000     -0.000      0.000 3/19965
  580.0      0.757      0.757      0.757      0.000     -0.000      0.000 3/19965
  590.0      0.744      0.744      0.744      0.000     -0.000      0.000 3/19965
  600.0      0.732      0.732      0.732      0.000     -0.000      0.000 3/19965
  610.0      0.720      0.720      0.720      0.000     -0.000      0.000 3/19965
  620.0      0.708      0.708      0.708      0.000     -0.000      0.000 3/19965
  630.0      0.697      0.697      0.697      0.000     -0.000      0.000 3/19965
  640.0      0.686      0.686      0.686      0.000     -0.000      0.000 3/19965
  650.0      0.676      0.676      0.676      0.000     -0.000      0.000 3/19965
  660.0      0.666      0.666      0.666      0.000     -0.000      0.000 3/19965
  670.0      0.656      0.656      0.656      0.000     -0.000      0.000 3/19965
  680.0      0.646      0.646      0.646      0.000     -0.000      0.000 3/19965
  690.0      0.637      0.637      0.637      0.000     -0.000      0.000 3/19965
  700.0      0.628      0.628      0.628      0.000     -0.000      0.000 3/19965
  710.0      0.619      0.619      0.619      0.000     -0.000      0.000 3/19965
  720.0      0.610      0.610      0.610      0.000     -0.000      0.000 3/19965
  730.0      0.602      0.602      0.602      0.000     -0.000      0.000 3/19965
  740.0      0.594      0.594      0.594      0.000     -0.000      0.000 3/19965
  750.0      0.586      0.586      0.586      0.000     -0.000      0.000 3/19965
  760.0      0.579      0.579      0.579      0.000     -0.000      0.000 3/19965
  770.0      0.571      0.571      0.571      0.000     -0.000      0.000 3/19965
  780.0      0.564      0.564      0.564      0.000     -0.000      0.000 3/19965
  790.0      0.557      0.557      0.557      0.000     -0.000      0.000 3/19965
  800.0      0.550      0.550      0.550      0.000     -0.000      0.000 3/19965
  810.0      0.543      0.543      0.543      0.000     -0.000      0.000 3/19965
  820.0      0.536      0.536      0.536      0.000     -0.000      0.000 3/19965
  830.0      0.530      0.530      0.530      0.000     -0.000      0.000 3/19965
  840.0      0.524      0.524      0.524      0.000     -0.000      0.000 3/19965
  850.0      0.518      0.518      0.518      0.000     -0.000      0.000 3/19965
  860.0      0.512      0.512      0.512      0.000     -0.000      0.000 3/19965
  870.0      0.506      0.506      0.506      0.000     -0.000      0.000 3/19965
  880.0      0.500      0.500      0.500      0.000     -0.000      0.000 3/19965
  890.0      0.494      0.494      0.494      0.000     -0.000      0.000 3/19965
  900.0      0.489      0.489      0.489      0.000     -0.000      0.000 3/19965
  910.0      0.484      0.484      0.484      0.000     -0.000      0.000 3/19965
  920.0      0.478      0.478      0.478      0.000     -0.000      0.000 3/19965
  930.0      0.473      0.473      0.473      0.000     -0.000      0.000 3/19965
  940.0      0.468      0.468      0.468      0.000     -0.000      0.000 3/19965
  950.0      0.463      0.463      0.463      0.000     -0.000      0.000 3/19965
  960.0      0.458      0.458      0.458      0.000     -0.000      0.000 3/19965
  970.0      0.454      0.454      0.454      0.000     -0.000      0.000 3/19965
  980.0      0.449      0.449      0.449      0.000     -0.000      0.000 3/19965
  990.0      0.445      0.445      0.445      0.000     -0.000      0.000 3/19965
 1000.0      0.440      0.440      0.440      0.000     -0.000      0.000 3/19965

Thermal conductivity related properties were written into 
"kappa-m111111.hdf5".
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-07 22:53:19]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|

