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

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

-------------------------[time 2026-01-07 21:04:04]-------------------------
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.066743169999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.066743169999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.066743169999999
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098
   *5 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.066743169999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.066743169999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.066743169999999
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098 > 4
   *5 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 5
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   12.200229509999996    0.000000000000000    0.000000000000000
  b    0.000000000000000   12.200229509999996    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.200229509999996
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098 > 4
   83 K   0.33333333333333  0.00000000000000  0.00000000000000  39.098 > 4
   84 K   0.66666666666667  0.00000000000000  0.00000000000000  39.098 > 4
   85 K   0.00000000000000  0.33333333333333  0.00000000000000  39.098 > 4
   86 K   0.33333333333333  0.33333333333333  0.00000000000000  39.098 > 4
   87 K   0.66666666666667  0.33333333333333  0.00000000000000  39.098 > 4
   88 K   0.00000000000000  0.66666666666667  0.00000000000000  39.098 > 4
   89 K   0.33333333333333  0.66666666666667  0.00000000000000  39.098 > 4
   90 K   0.66666666666667  0.66666666666667  0.00000000000000  39.098 > 4
   91 K   0.00000000000000  0.00000000000000  0.33333333333333  39.098 > 4
   92 K   0.33333333333333  0.00000000000000  0.33333333333333  39.098 > 4
   93 K   0.66666666666667  0.00000000000000  0.33333333333333  39.098 > 4
   94 K   0.00000000000000  0.33333333333333  0.33333333333333  39.098 > 4
   95 K   0.33333333333333  0.33333333333333  0.33333333333333  39.098 > 4
   96 K   0.66666666666667  0.33333333333333  0.33333333333333  39.098 > 4
   97 K   0.00000000000000  0.66666666666667  0.33333333333333  39.098 > 4
   98 K   0.33333333333333  0.66666666666667  0.33333333333333  39.098 > 4
   99 K   0.66666666666667  0.66666666666667  0.33333333333333  39.098 > 4
  100 K   0.00000000000000  0.00000000000000  0.66666666666667  39.098 > 4
  101 K   0.33333333333333  0.00000000000000  0.66666666666667  39.098 > 4
  102 K   0.66666666666667  0.00000000000000  0.66666666666667  39.098 > 4
  103 K   0.00000000000000  0.33333333333333  0.66666666666667  39.098 > 4
  104 K   0.33333333333333  0.33333333333333  0.66666666666667  39.098 > 4
  105 K   0.66666666666667  0.33333333333333  0.66666666666667  39.098 > 4
  106 K   0.00000000000000  0.66666666666667  0.66666666666667  39.098 > 4
  107 K   0.33333333333333  0.66666666666667  0.66666666666667  39.098 > 4
  108 K   0.66666666666667  0.66666666666667  0.66666666666667  39.098 > 4
 *109 Zn  0.16666666666667  0.16666666666667  0.16666666666667  65.380 > 5
  110 Zn  0.50000000000000  0.16666666666667  0.16666666666667  65.380 > 5
  111 Zn  0.83333333333333  0.16666666666667  0.16666666666667  65.380 > 5
  112 Zn  0.16666666666667  0.50000000000000  0.16666666666667  65.380 > 5
  113 Zn  0.50000000000000  0.50000000000000  0.16666666666667  65.380 > 5
  114 Zn  0.83333333333333  0.50000000000000  0.16666666666667  65.380 > 5
  115 Zn  0.16666666666667  0.83333333333333  0.16666666666667  65.380 > 5
  116 Zn  0.50000000000000  0.83333333333333  0.16666666666667  65.380 > 5
  117 Zn  0.83333333333333  0.83333333333333  0.16666666666667  65.380 > 5
  118 Zn  0.16666666666667  0.16666666666667  0.50000000000000  65.380 > 5
  119 Zn  0.50000000000000  0.16666666666667  0.50000000000000  65.380 > 5
  120 Zn  0.83333333333333  0.16666666666667  0.50000000000000  65.380 > 5
  121 Zn  0.16666666666667  0.50000000000000  0.50000000000000  65.380 > 5
  122 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 5
  123 Zn  0.83333333333333  0.50000000000000  0.50000000000000  65.380 > 5
  124 Zn  0.16666666666667  0.83333333333333  0.50000000000000  65.380 > 5
  125 Zn  0.50000000000000  0.83333333333333  0.50000000000000  65.380 > 5
  126 Zn  0.83333333333333  0.83333333333333  0.50000000000000  65.380 > 5
  127 Zn  0.16666666666667  0.16666666666667  0.83333333333333  65.380 > 5
  128 Zn  0.50000000000000  0.16666666666667  0.83333333333333  65.380 > 5
  129 Zn  0.83333333333333  0.16666666666667  0.83333333333333  65.380 > 5
  130 Zn  0.16666666666667  0.50000000000000  0.83333333333333  65.380 > 5
  131 Zn  0.50000000000000  0.50000000000000  0.83333333333333  65.380 > 5
  132 Zn  0.83333333333333  0.50000000000000  0.83333333333333  65.380 > 5
  133 Zn  0.16666666666667  0.83333333333333  0.83333333333333  65.380 > 5
  134 Zn  0.50000000000000  0.83333333333333  0.83333333333333  65.380 > 5
  135 Zn  0.83333333333333  0.83333333333333  0.83333333333333  65.380 > 5
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.4022284    0.0000000    0.0000000
            0.0000000    2.4022284    0.0000000
            0.0000000    0.0000000    2.4022284
-------------------------- Born effective charges --------------------------
    1 F    -1.7002887    0.0000000    0.0000000
            0.0000000   -0.8973290    0.0000000
            0.0000000    0.0000000   -0.8973290
    2 F    -0.8973290    0.0000000    0.0000000
            0.0000000   -1.7002887    0.0000000
            0.0000000    0.0000000   -0.8973290
    3 F    -0.8973290    0.0000000    0.0000000
            0.0000000   -0.8973290    0.0000000
            0.0000000    0.0000000   -1.7002887
    4 K     1.1962251    0.0000000    0.0000000
            0.0000000    1.1962251    0.0000000
            0.0000000    0.0000000    1.1962251
    5 Zn    2.2987216    0.0000000    0.0000000
            0.0000000    2.2987216    0.0000000
            0.0000000    0.0000000    2.2987216
----------------------------------------------------------------------------
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.059
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.063
--------------------------------- Symfc end --------------------------------
Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) 
Permutation basis: 405/405
Permutation basis: 3615/3615
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 104
Number of blocks in projector: 104
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 2
--- Eigsh_solver_block: 1 / 2 ---
Block_size: 74
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 2 ---
Block_size: 30
Use standard eigh solver.
Tree of FC basis block matrices:
- (104, 100), data: False
|-- (30, 30), data: True
|-- (74, 70), data: True
Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) 
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 21:04:09]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 21:04: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.066743169999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.066743169999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.066743169999999
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098
    5 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   12.200229509999996    0.000000000000000    0.000000000000000
  b    0.000000000000000   12.200229509999996    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.200229509999996
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098 > 82
   83 K   0.33333333333333  0.00000000000000  0.00000000000000  39.098 > 82
   84 K   0.66666666666667  0.00000000000000  0.00000000000000  39.098 > 82
   85 K   0.00000000000000  0.33333333333333  0.00000000000000  39.098 > 82
   86 K   0.33333333333333  0.33333333333333  0.00000000000000  39.098 > 82
   87 K   0.66666666666667  0.33333333333333  0.00000000000000  39.098 > 82
   88 K   0.00000000000000  0.66666666666667  0.00000000000000  39.098 > 82
   89 K   0.33333333333333  0.66666666666667  0.00000000000000  39.098 > 82
   90 K   0.66666666666667  0.66666666666667  0.00000000000000  39.098 > 82
   91 K   0.00000000000000  0.00000000000000  0.33333333333333  39.098 > 82
   92 K   0.33333333333333  0.00000000000000  0.33333333333333  39.098 > 82
   93 K   0.66666666666667  0.00000000000000  0.33333333333333  39.098 > 82
   94 K   0.00000000000000  0.33333333333333  0.33333333333333  39.098 > 82
   95 K   0.33333333333333  0.33333333333333  0.33333333333333  39.098 > 82
   96 K   0.66666666666667  0.33333333333333  0.33333333333333  39.098 > 82
   97 K   0.00000000000000  0.66666666666667  0.33333333333333  39.098 > 82
   98 K   0.33333333333333  0.66666666666667  0.33333333333333  39.098 > 82
   99 K   0.66666666666667  0.66666666666667  0.33333333333333  39.098 > 82
  100 K   0.00000000000000  0.00000000000000  0.66666666666667  39.098 > 82
  101 K   0.33333333333333  0.00000000000000  0.66666666666667  39.098 > 82
  102 K   0.66666666666667  0.00000000000000  0.66666666666667  39.098 > 82
  103 K   0.00000000000000  0.33333333333333  0.66666666666667  39.098 > 82
  104 K   0.33333333333333  0.33333333333333  0.66666666666667  39.098 > 82
  105 K   0.66666666666667  0.33333333333333  0.66666666666667  39.098 > 82
  106 K   0.00000000000000  0.66666666666667  0.66666666666667  39.098 > 82
  107 K   0.33333333333333  0.66666666666667  0.66666666666667  39.098 > 82
  108 K   0.66666666666667  0.66666666666667  0.66666666666667  39.098 > 82
  109 Zn  0.16666666666667  0.16666666666667  0.16666666666667  65.380 > 109
  110 Zn  0.50000000000000  0.16666666666667  0.16666666666667  65.380 > 109
  111 Zn  0.83333333333333  0.16666666666667  0.16666666666667  65.380 > 109
  112 Zn  0.16666666666667  0.50000000000000  0.16666666666667  65.380 > 109
  113 Zn  0.50000000000000  0.50000000000000  0.16666666666667  65.380 > 109
  114 Zn  0.83333333333333  0.50000000000000  0.16666666666667  65.380 > 109
  115 Zn  0.16666666666667  0.83333333333333  0.16666666666667  65.380 > 109
  116 Zn  0.50000000000000  0.83333333333333  0.16666666666667  65.380 > 109
  117 Zn  0.83333333333333  0.83333333333333  0.16666666666667  65.380 > 109
  118 Zn  0.16666666666667  0.16666666666667  0.50000000000000  65.380 > 109
  119 Zn  0.50000000000000  0.16666666666667  0.50000000000000  65.380 > 109
  120 Zn  0.83333333333333  0.16666666666667  0.50000000000000  65.380 > 109
  121 Zn  0.16666666666667  0.50000000000000  0.50000000000000  65.380 > 109
  122 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 109
  123 Zn  0.83333333333333  0.50000000000000  0.50000000000000  65.380 > 109
  124 Zn  0.16666666666667  0.83333333333333  0.50000000000000  65.380 > 109
  125 Zn  0.50000000000000  0.83333333333333  0.50000000000000  65.380 > 109
  126 Zn  0.83333333333333  0.83333333333333  0.50000000000000  65.380 > 109
  127 Zn  0.16666666666667  0.16666666666667  0.83333333333333  65.380 > 109
  128 Zn  0.50000000000000  0.16666666666667  0.83333333333333  65.380 > 109
  129 Zn  0.83333333333333  0.16666666666667  0.83333333333333  65.380 > 109
  130 Zn  0.16666666666667  0.50000000000000  0.83333333333333  65.380 > 109
  131 Zn  0.50000000000000  0.50000000000000  0.83333333333333  65.380 > 109
  132 Zn  0.83333333333333  0.50000000000000  0.83333333333333  65.380 > 109
  133 Zn  0.16666666666667  0.83333333333333  0.83333333333333  65.380 > 109
  134 Zn  0.50000000000000  0.83333333333333  0.83333333333333  65.380 > 109
  135 Zn  0.83333333333333  0.83333333333333  0.83333333333333  65.380 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            2.4022284    0.0000000    0.0000000
            0.0000000    2.4022284    0.0000000
            0.0000000    0.0000000    2.4022284
-------------------------- Born effective charges --------------------------
    1 F    -1.7002887    0.0000000    0.0000000
            0.0000000   -0.8973290    0.0000000
            0.0000000    0.0000000   -0.8973290
    2 F    -0.8973290    0.0000000    0.0000000
            0.0000000   -1.7002887    0.0000000
            0.0000000    0.0000000   -0.8973290
    3 F    -0.8973290    0.0000000    0.0000000
            0.0000000   -0.8973290    0.0000000
            0.0000000    0.0000000   -1.7002887
    4 K     1.1962251    0.0000000    0.0000000
            0.0000000    1.1962251    0.0000000
            0.0000000    0.0000000    1.1962251
    5 Zn    2.2987216    0.0000000    0.0000000
            0.0000000    2.2987216    0.0000000
            0.0000000    0.0000000    2.2987216
----------------------------------------------------------------------------
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 21:04:13]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-07 21:04: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.066743169999999    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.066743169999999    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.066743169999999
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098
    5 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   12.200229509999996    0.000000000000000    0.000000000000000
  b    0.000000000000000   12.200229509999996    0.000000000000000
  c    0.000000000000000    0.000000000000000   12.200229509999996
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 K   0.00000000000000  0.00000000000000  0.00000000000000  39.098 > 82
   83 K   0.33333333333333  0.00000000000000  0.00000000000000  39.098 > 82
   84 K   0.66666666666667  0.00000000000000  0.00000000000000  39.098 > 82
   85 K   0.00000000000000  0.33333333333333  0.00000000000000  39.098 > 82
   86 K   0.33333333333333  0.33333333333333  0.00000000000000  39.098 > 82
   87 K   0.66666666666667  0.33333333333333  0.00000000000000  39.098 > 82
   88 K   0.00000000000000  0.66666666666667  0.00000000000000  39.098 > 82
   89 K   0.33333333333333  0.66666666666667  0.00000000000000  39.098 > 82
   90 K   0.66666666666667  0.66666666666667  0.00000000000000  39.098 > 82
   91 K   0.00000000000000  0.00000000000000  0.33333333333333  39.098 > 82
   92 K   0.33333333333333  0.00000000000000  0.33333333333333  39.098 > 82
   93 K   0.66666666666667  0.00000000000000  0.33333333333333  39.098 > 82
   94 K   0.00000000000000  0.33333333333333  0.33333333333333  39.098 > 82
   95 K   0.33333333333333  0.33333333333333  0.33333333333333  39.098 > 82
   96 K   0.66666666666667  0.33333333333333  0.33333333333333  39.098 > 82
   97 K   0.00000000000000  0.66666666666667  0.33333333333333  39.098 > 82
   98 K   0.33333333333333  0.66666666666667  0.33333333333333  39.098 > 82
   99 K   0.66666666666667  0.66666666666667  0.33333333333333  39.098 > 82
  100 K   0.00000000000000  0.00000000000000  0.66666666666667  39.098 > 82
  101 K   0.33333333333333  0.00000000000000  0.66666666666667  39.098 > 82
  102 K   0.66666666666667  0.00000000000000  0.66666666666667  39.098 > 82
  103 K   0.00000000000000  0.33333333333333  0.66666666666667  39.098 > 82
  104 K   0.33333333333333  0.33333333333333  0.66666666666667  39.098 > 82
  105 K   0.66666666666667  0.33333333333333  0.66666666666667  39.098 > 82
  106 K   0.00000000000000  0.66666666666667  0.66666666666667  39.098 > 82
  107 K   0.33333333333333  0.66666666666667  0.66666666666667  39.098 > 82
  108 K   0.66666666666667  0.66666666666667  0.66666666666667  39.098 > 82
  109 Zn  0.16666666666667  0.16666666666667  0.16666666666667  65.380 > 109
  110 Zn  0.50000000000000  0.16666666666667  0.16666666666667  65.380 > 109
  111 Zn  0.83333333333333  0.16666666666667  0.16666666666667  65.380 > 109
  112 Zn  0.16666666666667  0.50000000000000  0.16666666666667  65.380 > 109
  113 Zn  0.50000000000000  0.50000000000000  0.16666666666667  65.380 > 109
  114 Zn  0.83333333333333  0.50000000000000  0.16666666666667  65.380 > 109
  115 Zn  0.16666666666667  0.83333333333333  0.16666666666667  65.380 > 109
  116 Zn  0.50000000000000  0.83333333333333  0.16666666666667  65.380 > 109
  117 Zn  0.83333333333333  0.83333333333333  0.16666666666667  65.380 > 109
  118 Zn  0.16666666666667  0.16666666666667  0.50000000000000  65.380 > 109
  119 Zn  0.50000000000000  0.16666666666667  0.50000000000000  65.380 > 109
  120 Zn  0.83333333333333  0.16666666666667  0.50000000000000  65.380 > 109
  121 Zn  0.16666666666667  0.50000000000000  0.50000000000000  65.380 > 109
  122 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 109
  123 Zn  0.83333333333333  0.50000000000000  0.50000000000000  65.380 > 109
  124 Zn  0.16666666666667  0.83333333333333  0.50000000000000  65.380 > 109
  125 Zn  0.50000000000000  0.83333333333333  0.50000000000000  65.380 > 109
  126 Zn  0.83333333333333  0.83333333333333  0.50000000000000  65.380 > 109
  127 Zn  0.16666666666667  0.16666666666667  0.83333333333333  65.380 > 109
  128 Zn  0.50000000000000  0.16666666666667  0.83333333333333  65.380 > 109
  129 Zn  0.83333333333333  0.16666666666667  0.83333333333333  65.380 > 109
  130 Zn  0.16666666666667  0.50000000000000  0.83333333333333  65.380 > 109
  131 Zn  0.50000000000000  0.50000000000000  0.83333333333333  65.380 > 109
  132 Zn  0.83333333333333  0.50000000000000  0.83333333333333  65.380 > 109
  133 Zn  0.16666666666667  0.83333333333333  0.83333333333333  65.380 > 109
  134 Zn  0.50000000000000  0.83333333333333  0.83333333333333  65.380 > 109
  135 Zn  0.83333333333333  0.83333333333333  0.83333333333333  65.380 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            2.4022284    0.0000000    0.0000000
            0.0000000    2.4022284    0.0000000
            0.0000000    0.0000000    2.4022284
-------------------------- Born effective charges --------------------------
    1 F    -1.7002887    0.0000000    0.0000000
            0.0000000   -0.8973290    0.0000000
            0.0000000    0.0000000   -0.8973290
    2 F    -0.8973290    0.0000000    0.0000000
            0.0000000   -1.7002887    0.0000000
            0.0000000    0.0000000   -0.8973290
    3 F    -0.8973290    0.0000000    0.0000000
            0.0000000   -0.8973290    0.0000000
            0.0000000    0.0000000   -1.7002887
    4 K     1.1962251    0.0000000    0.0000000
            0.0000000    1.1962251    0.0000000
            0.0000000    0.0000000    1.1962251
    5 Zn    2.2987216    0.0000000    0.0000000
            0.0000000    2.2987216    0.0000000
            0.0000000    0.0000000    2.2987216
----------------------------------------------------------------------------
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 (zz) -0.00000000 (zz) 
--------------------------- 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 ... [ 12 12 12 ]
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: 1.02, Number of G-points: 305, Lambda: 0.16
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/84) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 84
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
  -0.000   (   0.000    0.000    0.000)    0.000
   3.859   (   0.000    0.000    0.000)    0.000
   3.859   (   0.000    0.000    0.000)    0.000
   3.859   (   0.000    0.000    0.000)    0.000
   4.489   (   0.000    0.000    0.000)    0.000
   4.489   (   0.000    0.000    0.000)    0.000
   4.489   (   0.000    0.000    0.000)    0.000
   5.437   (   0.000    0.000    0.000)    0.000
   5.437   (   0.000    0.000    0.000)    0.000
   5.437   (   0.000    0.000    0.000)    0.000
  12.501   (   0.000    0.000    0.000)    0.000
  12.501   (   0.000    0.000    0.000)    0.000
  12.501   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/84) =======================
q-point: ( 0.08  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 168
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.597   (  28.702    0.000    0.000)   28.702
   0.597   (  28.702    0.000    0.000)   28.702
   1.189   (  56.423    0.000    0.000)   56.423
   3.844   (  -1.565    0.000    0.000)    1.565
   3.844   (  -1.565    0.000    0.000)    1.565
   4.312   (  11.763    0.000    0.000)   11.763
   4.511   (   2.187    0.000    0.000)    2.187
   4.511   (   2.187    0.000    0.000)    2.187
   4.511   (   2.146    0.000    0.000)    2.146
   5.430   (  -0.668    0.000    0.000)    0.668
   5.430   (  -0.668    0.000    0.000)    0.668
   8.206   (   0.315    0.000    0.000)    0.315
  12.494   (  -0.649    0.000    0.000)    0.649
  12.494   (  -0.649    0.000    0.000)    0.649
  14.520   (  -6.836    0.000    0.000)    6.836
======================= Grid point 2 (3/84) =======================
q-point: ( 0.17  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.168   (  26.914    0.000    0.000)   26.914
   1.168   (  26.914    0.000    0.000)   26.914
   2.282   (  49.661    0.000    0.000)   49.661
   3.791   (  -3.700    0.000    0.000)    3.700
   3.791   (  -3.700    0.000    0.000)    3.700
   4.572   (   3.660    0.000    0.000)    3.660
   4.579   (   4.448    0.000    0.000)    4.448
   4.579   (   4.448    0.000    0.000)    4.448
   4.645   (  19.722    0.000    0.000)   19.722
   5.409   (  -1.392    0.000    0.000)    1.392
   5.409   (  -1.392    0.000    0.000)    1.392
   8.225   (   1.933    0.000    0.000)    1.933
  12.475   (  -1.136    0.000    0.000)    1.136
  12.475   (  -1.136    0.000    0.000)    1.136
  14.315   ( -12.877    0.000    0.000)   12.877
======================= Grid point 3 (4/84) =======================
q-point: ( 0.25  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 168
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.693   (  24.174    0.000    0.000)   24.174
   1.693   (  24.174    0.000    0.000)   24.174
   3.206   (  40.307    0.000    0.000)   40.307
   3.690   (  -6.205    0.000    0.000)    6.205
   3.690   (  -6.205    0.000    0.000)    6.205
   4.654   (   4.141    0.000    0.000)    4.141
   4.690   (   6.233    0.000    0.000)    6.233
   4.690   (   6.233    0.000    0.000)    6.233
   5.068   (  20.016    0.000    0.000)   20.016
   5.372   (  -2.237    0.000    0.000)    2.237
   5.372   (  -2.237    0.000    0.000)    2.237
   8.302   (   5.989    0.000    0.000)    5.989
  12.449   (  -1.332    0.000    0.000)    1.332
  12.449   (  -1.332    0.000    0.000)    1.332
  14.005   ( -17.025    0.000    0.000)   17.025
======================= Grid point 4 (5/84) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 196
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.153   (  20.589    0.000    0.000)   20.589
   2.153   (  20.589    0.000    0.000)   20.589
   3.538   (  -8.482    0.000    0.000)    8.482
   3.538   (  -8.482    0.000    0.000)    8.482
   3.930   (  30.268    0.000    0.000)   30.268
   4.734   (   3.516    0.000    0.000)    3.516
   4.827   (   6.958    0.000    0.000)    6.958
   4.827   (   6.958    0.000    0.000)    6.958
   5.316   (  -3.216    0.000    0.000)    3.216
   5.316   (  -3.216    0.000    0.000)    3.216
   5.404   (  11.772    0.000    0.000)   11.772
   8.478   (  10.904    0.000    0.000)   10.904
  12.423   (  -1.172    0.000    0.000)    1.172
  12.423   (  -1.172    0.000    0.000)    1.172
  13.643   ( -17.506    0.000    0.000)   17.506
======================= Grid point 5 (6/84) =======================
q-point: ( 0.42  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 168
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.525   (  15.095    0.000    0.000)   15.095
   2.525   (  15.095    0.000    0.000)   15.095
   3.353   (  -9.033    0.000    0.000)    9.033
   3.353   (  -9.033    0.000    0.000)    9.033
   4.435   (  18.274    0.000    0.000)   18.274
   4.792   (   2.001    0.000    0.000)    2.001
   4.965   (   6.240    0.000    0.000)    6.240
   4.965   (   6.240    0.000    0.000)    6.240
   5.241   (  -3.988    0.000    0.000)    3.988
   5.241   (  -3.988    0.000    0.000)    3.988
   5.541   (   2.196    0.000    0.000)    2.196
   8.708   (   9.922    0.000    0.000)    9.922
  12.404   (  -0.684    0.000    0.000)    0.684
  12.404   (  -0.684    0.000    0.000)    0.684
  13.329   ( -11.910    0.000    0.000)   11.910
======================= Grid point 6 (7/84) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 112
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.707   (   0.000    0.000    0.000)    0.000
   2.707   (   0.000    0.000    0.000)    0.000
   3.233   (   0.000    0.000    0.000)    0.000
   3.233   (   0.000    0.000    0.000)    0.000
   4.634   (  -0.000    0.000    0.000)    0.000
   4.813   (   0.000    0.000    0.000)    0.000
   5.050   (  -0.000    0.000    0.000)    0.000
   5.050   (  -0.000    0.000    0.000)    0.000
   5.181   (   0.000    0.000    0.000)    0.000
   5.181   (   0.000    0.000    0.000)    0.000
   5.549   (   0.000    0.000    0.000)    0.000
   8.820   (   0.000    0.000    0.000)    0.000
  12.396   (   0.000    0.000    0.000)    0.000
  12.396   (   0.000    0.000    0.000)    0.000
  13.199   (   0.000    0.000    0.000)    0.000
======================= Grid point 14 (8/84) =======================
q-point: ( 0.08  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 294
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.843   (  20.249   20.249    0.000)   28.636
   0.922   (  21.567   21.567    0.000)   30.500
   1.625   (  38.125   38.125    0.000)   53.917
   3.830   (  -1.386   -1.386    0.000)    1.960
   3.968   (   4.667    4.667    0.000)    6.600
   4.207   (   0.802    0.802    0.000)    1.134
   4.435   (  -2.734   -2.734    0.000)    3.867
   4.531   (   2.011    2.011    0.000)    2.845
   4.625   (   6.428    6.428    0.000)    9.091
   5.423   (  -0.651   -0.651    0.000)    0.921
   5.513   (   4.181    4.181    0.000)    5.913
   8.203   (  -0.001   -0.001    0.000)    0.001
  12.307   (  -9.449   -9.449    0.000)   13.362
  12.487   (  -0.652   -0.652    0.000)    0.923
  14.600   (   0.427    0.427    0.000)    0.603
======================= Grid point 15 (9/84) =======================
q-point: ( 0.17  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.310   (  23.909   12.897    0.000)   27.166
   1.377   (  22.802   16.666    0.000)   28.243
   2.500   (  43.360   20.625    0.000)   48.015
   3.783   (  -3.403   -0.896    0.000)    3.519
   3.910   (  -5.821   10.793    0.000)   12.262
   4.387   (   0.544   -9.998    0.000)   10.013
   4.477   (  16.825  -10.934    0.000)   20.066
   4.594   (   4.170    1.547    0.000)    4.448
   4.771   (   7.762   10.189    0.000)   12.809
   5.403   (  -1.364   -0.606    0.000)    1.492
   5.580   (   3.037   12.368    0.000)   12.735
   8.216   (   1.545   -0.766    0.000)    1.724
  12.165   (  -5.442  -22.760    0.000)   23.402
  12.468   (  -1.142   -0.662    0.000)    1.320
  14.497   (  -9.371   11.056    0.000)   14.493
======================= Grid point 16 (10/84) =======================
q-point: ( 0.25  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.792   (  22.741    9.247    0.000)   24.549
   1.845   (  22.487   12.673    0.000)   25.813
   3.322   (  35.949   11.069    0.000)   37.615
   3.688   (  -5.893   -0.255    0.000)    5.898
   3.776   (  -7.338    8.205    0.000)   11.008
   4.391   (   0.801  -14.735    0.000)   14.756
   4.700   (   5.967    0.981    0.000)    6.047
   4.842   (  16.202  -13.060    0.000)   20.810
   4.937   (   8.320   13.437    0.000)   15.804
   5.366   (  -2.203   -0.543    0.000)    2.269
   5.655   (   4.488   18.777    0.000)   19.306
   8.283   (   5.485   -1.742    0.000)    5.755
  12.059   (  -5.188  -30.736    0.000)   31.171
  12.442   (  -1.339   -0.675    0.000)    1.500
  14.252   ( -13.756   17.300    0.000)   22.102
======================= Grid point 17 (11/84) =======================
q-point: ( 0.33  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.231   (  19.845    7.275    0.000)   21.136
   2.287   (  20.400   11.215    0.000)   23.280
   3.542   (  -8.258    0.332    0.000)    8.264
   3.604   (  -9.376    6.413    0.000)   11.360
   3.954   (  25.323    2.338    0.000)   25.431
   4.419   (   1.737  -18.996    0.000)   19.075
   4.832   (   6.766    0.503    0.000)    6.784
   5.085   (   6.978   -2.728    0.000)    7.493
   5.118   (   9.740   10.011    0.000)   13.967
   5.311   (  -3.182   -0.474    0.000)    3.218
   5.749   (   3.685   20.625    0.000)   20.951
   8.449   (  10.442   -2.790    0.000)   10.808
  11.949   (  -5.398  -36.763    0.000)   37.157
  12.416   (  -1.178   -0.688    0.000)    1.365
  13.963   ( -13.730   22.108    0.000)   26.025
======================= Grid point 18 (12/84) =======================
q-point: ( 0.42  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.593   (  14.974    6.403    0.000)   16.285
   2.669   (  16.408   12.152    0.000)   20.418
   3.358   (  -9.127    0.564    0.000)    9.145
   3.395   ( -10.689    4.458    0.000)   11.582
   4.341   (  11.618  -11.368    0.000)   16.255
   4.461   (   2.806  -16.887    0.000)   17.118
   4.967   (   6.148    0.199    0.000)    6.151
   5.174   (   2.292    2.032    0.000)    3.063
   5.237   (  -3.973   -0.417    0.000)    3.995
   5.309   (   8.639   16.174    0.000)   18.336
   5.772   (  -1.588   18.442    0.000)   18.511
   8.671   (   9.695   -3.507    0.000)   10.310
  11.847   (  -4.138  -41.199    0.000)   41.406
  12.396   (  -0.688   -0.699    0.000)    0.981
  13.724   (  -8.739   25.720    0.000)   27.164
======================= Grid point 19 (13/84) =======================
q-point: (-0.50  0.08  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 294
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.778   (   0.000    6.698    0.000)    6.698
   2.885   (   0.000   15.977    0.000)   15.977
   3.234   (   0.000    0.079    0.000)    0.079
   3.238   (   0.000    0.404    0.000)    0.404
   4.431   (  -0.000  -25.442    0.000)   25.442
   4.520   (   0.000   -7.942    0.000)    7.942
   5.051   (  -0.000    0.145    0.000)    0.145
   5.176   (   0.000   -0.443    0.000)    0.443
   5.196   (   0.000    4.198    0.000)    4.198
   5.421   (  -0.000   20.285    0.000)   20.285
   5.732   (   0.000   16.310    0.000)   16.310
   8.782   (   0.000   -3.720    0.000)    3.720
  11.801   (   0.000  -42.907    0.000)   42.907
  12.389   (   0.000   -0.703    0.000)    0.703
  13.631   (   0.000   27.101    0.000)   27.101
======================= Grid point 28 (14/84) =======================
q-point: ( 0.17  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 301
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.646   (  18.767   18.767    0.000)   26.540
   1.741   (  18.187   18.187    0.000)   25.720
   3.048   (  29.868   29.868    0.000)   42.239
   3.749   (  -2.505   -2.505    0.000)    3.542
   4.111   (  -1.734   -1.734    0.000)    2.453
   4.249   (   1.509    1.509    0.000)    2.134
   4.307   (  -0.297   -0.297    0.000)    0.420
   4.644   (   3.368    3.368    0.000)    4.763
   4.993   (  11.368   11.368    0.000)   16.077
   5.384   (  -1.284   -1.284    0.000)    1.817
   5.828   (  11.463   11.463    0.000)   16.211
   8.205   (   0.104    0.104    0.000)    0.148
  11.732   ( -18.445  -18.445    0.000)   26.086
  12.449   (  -1.159   -1.159    0.000)    1.639
  14.616   (   0.240    0.240    0.000)    0.340
======================= Grid point 29 (15/84) =======================
q-point: ( 0.25  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.045   (  19.578   14.669    0.000)   24.464
   2.123   (  19.138   13.408    0.000)   23.367
   3.609   (  22.172   14.419    0.000)   26.448
   3.675   (  -4.840   -1.142    0.000)    4.973
   3.996   (  -4.253    5.537    0.000)    6.982
   4.115   (  -7.663   -3.383    0.000)    8.376
   4.613   (  17.837   -9.412    0.000)   20.168
   4.732   (   5.121    2.252    0.000)    5.594
   5.234   (  12.142   15.503    0.000)   19.692
   5.349   (  -2.103   -1.168    0.000)    2.406
   6.047   (   9.699   19.058    0.000)   21.384
   8.239   (   3.654   -2.316    0.000)    4.326
  11.402   ( -14.225  -31.127    0.000)   34.223
  12.423   (  -1.359   -1.182    0.000)    1.801
  14.533   (  -7.390    8.769    0.000)   11.468
======================= Grid point 30 (16/84) =======================
q-point: ( 0.33  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 518
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.432   (  17.934   11.790    0.000)   21.463
   2.520   (  19.349    9.839    0.000)   21.707
   3.549   (  -7.377    0.327    0.000)    7.384
   3.755   (  -6.448    5.106    0.000)    8.225
   3.972   (  -4.575  -11.148    0.000)   12.051
   4.114   (  12.325    3.246    0.000)   12.745
   4.849   (   6.113    1.191    0.000)    6.228
   4.945   (  13.708   -5.501    0.000)   14.771
   5.296   (  -3.074   -1.025    0.000)    3.241
   5.491   (  12.940   20.294    0.000)   24.068
   6.198   (   4.234   22.987    0.000)   23.374
   8.368   (   8.882   -4.941    0.000)   10.164
  11.138   ( -11.528  -40.170    0.000)   41.792
  12.396   (  -1.196   -1.207    0.000)    1.699
  14.352   (  -9.477   14.433    0.000)   17.266
======================= Grid point 31 (17/84) =======================
q-point: ( 0.42  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.769   (  14.582   10.159    0.000)   17.771
   2.906   (  18.073    8.158    0.000)   19.829
   3.377   (  -9.170    1.298    0.000)    9.261
   3.547   ( -12.595    9.865    0.000)   15.998
   3.908   (  -1.749  -23.549    0.000)   23.614
   4.319   (   6.698   -1.630    0.000)    6.893
   4.974   (   5.805    0.448    0.000)    5.822
   5.155   (   6.745   -1.748    0.000)    6.967
   5.224   (  -3.902   -0.877    0.000)    3.999
   5.757   (  12.506   24.654    0.000)   27.645
   6.204   (  -3.468   23.140    0.000)   23.399
   8.567   (   8.998   -6.578    0.000)   11.146
  10.939   (  -7.460  -45.633    0.000)   46.239
  12.376   (  -0.699   -1.225    0.000)    1.411
  14.180   (  -6.502   17.660    0.000)   18.818
======================= Grid point 32 (18/84) =======================
q-point: (-0.50  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 301
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.966   (   0.000   11.137    0.000)   11.137
   3.236   (   0.000    0.123    0.000)    0.123
   3.236   (  -0.000   -1.573    0.000)    1.573
   3.273   (   0.000   19.428    0.000)   19.428
   3.891   (   0.000  -25.574    0.000)   25.574
   4.386   (   0.000   -3.949    0.000)    3.949
   5.055   (  -0.000    0.261    0.000)    0.261
   5.162   (   0.000   -0.889    0.000)    0.889
   5.224   (   0.000   -0.350    0.000)    0.350
   5.935   (  -0.000   28.893    0.000)   28.893
   6.119   (   0.000   20.050    0.000)   20.050
   8.671   (   0.000   -7.007    0.000)    7.007
  10.857   (   0.000  -47.546    0.000)   47.546
  12.369   (   0.000   -1.232    0.000)    1.232
  14.110   (   0.000   18.714    0.000)   18.714
======================= Grid point 42 (19/84) =======================
q-point: ( 0.25  0.25  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 294
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.366   (  16.155   16.155    0.000)   22.846
   2.390   (  12.904   12.904    0.000)   18.249
   3.635   (  -2.911   -2.911    0.000)    4.116
   3.795   ( -10.443  -10.443    0.000)   14.768
   3.828   (   5.576    5.576    0.000)    7.885
   4.438   (   4.285    4.285    0.000)    6.061
   4.475   (  11.189   11.189    0.000)   15.824
   4.793   (   3.728    3.728    0.000)    5.272
   5.318   (  -1.940   -1.940    0.000)    2.744
   5.569   (  17.045   17.045    0.000)   24.106
   6.419   (  16.440   16.440    0.000)   23.250
   8.207   (  -0.168   -0.168    0.000)    0.238
  10.818   ( -25.712  -25.712    0.000)   36.362
  12.396   (  -1.387   -1.387    0.000)    1.962
  14.611   (  -0.564   -0.564    0.000)    0.798
======================= Grid point 43 (20/84) =======================
q-point: ( 0.33  0.25  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.659   (  13.245    3.470    0.000)   13.692
   2.690   (  15.207   12.903    0.000)   19.943
   3.551   (  -5.308   -0.266    0.000)    5.315
   3.615   (  -6.642  -12.572    0.000)   14.219
   3.848   (  -3.819   -1.015    0.000)    3.951
   4.384   (  -0.534   16.347    0.000)   16.356
   4.861   (  16.416   -2.469    0.000)   16.601
   4.883   (   4.872    2.124    0.000)    5.315
   5.269   (  -2.863   -1.704    0.000)    3.332
   5.924   (  17.301   21.398    0.000)   27.517
   6.677   (   7.660   23.040    0.000)   24.281
   8.257   (   5.427   -5.475    0.000)    7.709
  10.351   ( -19.670  -35.957    0.000)   40.986
  12.369   (  -1.222   -1.417    0.000)    1.871
  14.544   (  -5.072    4.732    0.000)    6.937
======================= Grid point 44 (21/84) =======================
q-point: ( 0.42  0.25  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.916   (  10.998   -9.001    0.000)   14.212
   2.982   (  13.094    9.844    0.000)   16.382
   3.416   (  -7.927    2.653    0.000)    8.359
   3.524   (  -2.075  -10.774    0.000)   10.972
   3.711   (  -8.318    3.833    0.000)    9.159
   4.400   (   1.215   10.201    0.000)   10.273
   4.986   (   5.038    0.791    0.000)    5.100
   5.116   (   8.376   -1.978    0.000)    8.606
   5.201   (  -3.692   -1.388    0.000)    3.944
   6.264   (  15.577   23.889    0.000)   28.519
   6.713   (  -3.813   26.246    0.000)   26.522
   8.407   (   7.653   -8.945    0.000)   11.772
  10.024   ( -11.810  -42.853    0.000)   44.451
  12.348   (  -0.714   -1.439    0.000)    1.607
  14.439   (  -4.358    7.922    0.000)    9.042
======================= Grid point 45 (22/84) =======================
q-point: (-0.50  0.25  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 294
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.047   (   0.000  -19.846    0.000)   19.846
   3.212   (  -0.000   12.364    0.000)   12.364
   3.238   (   0.000    0.121    0.000)    0.121
   3.529   (   0.000   -6.365    0.000)    6.365
   3.580   (   0.000    8.631    0.000)    8.631
   4.414   (   0.000    7.865    0.000)    7.865
   5.061   (  -0.000    0.316    0.000)    0.316
   5.140   (   0.000   -1.279    0.000)    1.279
   5.202   (   0.000   -1.523    0.000)    1.523
   6.530   (  -0.000   19.840    0.000)   19.840
   6.565   (   0.000   31.550    0.000)   31.550
   8.499   (   0.000   -9.713    0.000)    9.713
   9.896   (   0.000  -45.403    0.000)   45.403
  12.341   (   0.000   -1.447    0.000)    1.447
  14.391   (   0.000    8.980    0.000)    8.980
======================= Grid point 56 (23/84) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 301
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.663   (  -2.558   -2.558    0.000)    3.618
   2.943   (  11.417   11.417    0.000)   16.146
   3.449   (  -3.146   -3.146    0.000)    4.450
   3.534   (  -1.469   -1.469    0.000)    2.078
   3.800   (  -3.611   -3.611    0.000)    5.107
   4.645   (   6.177    6.177    0.000)    8.735
   4.890   (   7.217    7.217    0.000)   10.206
   4.936   (   3.107    3.107    0.000)    4.394
   5.226   (  -2.490   -2.490    0.000)    3.522
   6.348   (  19.684   19.684    0.000)   27.838
   7.089   (  15.292   15.292    0.000)   21.626
   8.172   (  -1.827   -1.827    0.000)    2.584
   9.687   ( -28.338  -28.338    0.000)   40.076
  12.341   (  -1.248   -1.248    0.000)    1.765
  14.572   (  -1.270   -1.270    0.000)    1.797
======================= Grid point 57 (24/84) =======================
q-point: ( 0.42  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.594   (  -3.160  -18.458    0.000)   18.727
   3.146   (   7.998    5.792    0.000)    9.875
   3.459   (   2.978    1.400    0.000)    3.290
   3.491   (  -2.639    4.573    0.000)    5.280
   3.710   (  -4.885   -2.552    0.000)    5.511
   4.683   (   0.029   15.252    0.000)   15.252
   5.007   (   3.670    1.247    0.000)    3.876
   5.080   (   8.156   -1.383    0.000)    8.272
   5.167   (  -3.185   -1.881    0.000)    3.699
   6.717   (  15.818   20.235    0.000)   25.684
   7.241   (  -0.409   23.736    0.000)   23.740
   8.209   (   4.639  -10.015    0.000)   11.037
   9.215   ( -16.939  -35.172    0.000)   39.039
  12.320   (  -0.730   -1.268    0.000)    1.463
  14.526   (  -2.416    1.209    0.000)    2.702
======================= Grid point 58 (25/84) =======================
q-point: (-0.50  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 301
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.559   (   0.000  -24.168    0.000)   24.168
   3.240   (  -0.000    0.087    0.000)    0.087
   3.451   (   0.000   10.501    0.000)   10.501
   3.514   (   0.000    1.819    0.000)    1.819
   3.638   (   0.000   -1.140    0.000)    1.140
   4.682   (   0.000   16.070    0.000)   16.070
   5.068   (  -0.000    0.286    0.000)    0.286
   5.112   (   0.000   -1.423    0.000)    1.423
   5.168   (   0.000   -1.640    0.000)    1.640
   6.927   (  -0.000   18.753    0.000)   18.753
   7.179   (   0.000   27.204    0.000)   27.204
   8.276   (   0.000  -11.867    0.000)   11.867
   9.032   (   0.000  -37.981    0.000)   37.981
  12.312   (   0.000   -1.276    0.000)    1.276
  14.497   (   0.000    2.062    0.000)    2.062
======================= Grid point 70 (26/84) =======================
q-point: ( 0.42  0.42  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 294
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.266   ( -12.083  -12.083    0.000)   17.087
   3.222   (   2.048    2.048    0.000)    2.897
   3.509   (   2.854    2.854    0.000)    4.036
   3.583   (   3.795    3.795    0.000)    5.368
   3.645   (  -3.265   -3.265    0.000)    4.617
   4.941   (   7.478    7.478    0.000)   10.575
   5.038   (   1.749    1.749    0.000)    2.473
   5.080   (   2.561    2.561    0.000)    3.622
   5.125   (  -2.178   -2.178    0.000)    3.079
   7.090   (  15.752   15.752    0.000)   22.277
   7.606   (   9.363    9.363    0.000)   13.242
   8.033   (  -5.132   -5.132    0.000)    7.257
   8.622   ( -21.632  -21.632    0.000)   30.593
  12.299   (  -0.742   -0.742    0.000)    1.049
  14.520   (  -1.086   -1.086    0.000)    1.536
======================= Grid point 71 (27/84) =======================
q-point: (-0.50  0.42  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 294
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.137   (   0.000  -15.580    0.000)   15.580
   3.242   (   0.000    0.043    0.000)    0.043
   3.554   (   0.000    1.693    0.000)    1.693
   3.596   (  -0.000   -2.061    0.000)    2.061
   3.625   (   0.000    6.070    0.000)    6.070
   4.991   (   0.000   12.415    0.000)   12.415
   5.073   (  -0.000    0.171    0.000)    0.171
   5.086   (   0.000   -1.015    0.000)    1.015
   5.140   (   0.000   -1.029    0.000)    1.029
   7.291   (  -0.000   16.444    0.000)   16.444
   7.636   (   0.000   16.507    0.000)   16.507
   8.019   (  -0.000  -12.941    0.000)   12.941
   8.385   (   0.000  -23.954    0.000)   23.954
  12.291   (   0.000   -0.746    0.000)    0.746
  14.505   (   0.000   -0.617    0.000)    0.617
======================= Grid point 90 (28/84) =======================
q-point: (-0.50 -0.50  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 112
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.970   (   0.000    0.000    0.000)    0.000
   3.242   (   0.000    0.000    0.000)    0.000
   3.573   (   0.000    0.000    0.000)    0.000
   3.573   (   0.000    0.000    0.000)    0.000
   3.689   (   0.000    0.000    0.000)    0.000
   5.074   (   0.000    0.000    0.000)    0.000
   5.074   (   0.000    0.000    0.000)    0.000
   5.129   (   0.000    0.000    0.000)    0.000
   5.129   (   0.000    0.000    0.000)    0.000
   7.524   (  -0.000   -0.000    0.000)    0.000
   7.822   (   0.000    0.000    0.000)    0.000
   7.822   (   0.000    0.000    0.000)    0.000
   8.115   (   0.000    0.000    0.000)    0.000
  12.283   (   0.000    0.000    0.000)    0.000
  14.495   (   0.000    0.000    0.000)    0.000
======================= Grid point 183 (29/84) =======================
q-point: ( 0.08  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 182
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.089   (  16.825   16.825   16.825)   29.143
   1.089   (  16.825   16.825   16.825)   29.143
   1.965   (  30.579   30.579   30.579)   52.964
   3.949   (   2.390    2.390    2.390)    4.140
   3.949   (   2.390    2.390    2.390)    4.140
   4.108   (  -2.525   -2.525   -2.525)    4.374
   4.365   (  -3.991   -3.991   -3.991)    6.912
   4.636   (   4.358    4.358    4.358)    7.548
   4.636   (   4.358    4.358    4.358)    7.548
   5.511   (   2.881    2.881    2.881)    4.991
   5.511   (   2.881    2.881    2.881)    4.991
   8.198   (  -0.260   -0.260   -0.260)    0.451
  12.299   (  -6.563   -6.563   -6.563)   11.368
  12.299   (  -6.563   -6.563   -6.563)   11.368
  14.680   (   2.817    2.817    2.817)    4.879
======================= Grid point 184 (30/84) =======================
q-point: ( 0.17  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 474
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.466   (  20.645   12.411   12.411)   27.098
   1.506   (  21.370   12.320   12.320)   27.572
   2.709   (  38.385   19.596   19.596)   47.344
   3.880   (  -6.607    3.004    3.004)    7.855
   3.892   (  -3.886    4.247    4.247)    7.153
   4.300   (  -2.307   -6.940   -6.940)   10.083
   4.332   (  16.025   -9.928   -9.928)   21.306
   4.694   (   4.805    5.411    5.411)    9.036
   4.797   (   8.109    4.339    4.339)   10.169
   5.507   (  -0.190    5.163    5.163)    7.304
   5.624   (   5.383    6.440    6.440)   10.579
   8.204   (   1.153   -0.983   -0.983)    1.806
  12.082   (  -8.998  -11.929  -11.929)   19.119
  12.278   (  -1.269   -9.583   -9.583)   13.612
  14.637   (  -6.536    9.277    9.277)   14.657
======================= Grid point 185 (31/84) =======================
q-point: ( 0.25  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.901   (  21.280    8.978    8.978)   24.780
   1.950   (  21.582    9.892    9.892)   25.719
   3.441   (  31.434   10.903   10.903)   35.012
   3.746   (  -6.208    2.849    2.849)    7.401
   3.786   (  -6.520    4.127    4.127)    8.751
   4.272   (  -0.428  -10.079  -10.079)   14.260
   4.599   (   8.531  -10.829  -10.829)   17.530
   4.871   (  12.598    4.497    4.497)   14.112
   4.957   (   7.297    4.670    4.670)    9.842
   5.505   (   0.003    6.747    6.747)    9.541
   5.743   (   6.203   10.566   10.566)   16.179
   8.262   (   5.004   -1.933   -1.933)    5.702
  11.911   (  -8.064  -18.530  -18.530)   27.417
  12.249   (  -1.495   -9.726   -9.726)   13.835
  14.445   ( -11.413   14.264   14.264)   23.177
======================= Grid point 186 (32/84) =======================
q-point: ( 0.33  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 474
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.323   (  19.596    7.268    7.268)   22.128
   2.377   (  19.815    8.729    8.729)   23.346
   3.575   (  -8.313    1.146    1.146)    8.469
   3.625   (  -9.126    3.835    3.835)   10.616
   4.014   (  21.926    5.872    5.872)   23.446
   4.277   (   0.700  -12.650  -12.650)   17.903
   4.714   (   3.457  -10.657  -10.657)   15.463
   5.087   (   5.150    3.505    3.505)    7.148
   5.146   (  13.233    5.966    5.966)   15.694
   5.509   (   0.389    9.316    9.316)   13.180
   5.855   (   3.675   11.902   11.902)   17.229
   8.418   (  10.008   -2.928   -2.928)   10.831
  11.747   (  -7.806  -23.338  -23.338)   33.915
  12.220   (  -1.321   -9.877   -9.877)   14.030
  14.201   ( -11.696   17.885   17.885)   27.867
======================= Grid point 187 (33/84) =======================
q-point: ( 0.42  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.691   (  15.848    6.984    6.984)   18.674
   2.751   (  16.330    8.660    8.660)   20.412
   3.382   ( -10.309    1.771    1.771)   10.609
   3.420   ( -10.479    3.101    3.101)   11.360
   4.293   (   0.721  -14.310  -14.310)   20.251
   4.362   (  12.043   -0.058   -0.058)   12.043
   4.748   (  -0.035  -11.579  -11.579)   16.375
   5.164   (   2.445    1.635    1.635)    3.365
   5.406   (  11.891   10.790   10.790)   19.346
   5.519   (   0.523   12.047   12.047)   17.046
   5.863   (  -2.892   10.936   10.936)   15.734
   8.634   (   9.483   -3.603   -3.603)   10.765
  11.604   (  -5.651  -26.534  -26.534)   37.948
  12.198   (  -0.774   -9.993   -9.993)   14.153
  13.997   (  -7.451   20.294   20.294)   29.651
======================= Grid point 188 (34/84) =======================
q-point: (-0.50  0.08  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 276
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.898   (   0.000    8.210    8.210)   11.611
   2.980   (   0.000   11.293   11.293)   15.971
   3.214   (   0.000   -0.594   -0.594)    0.840
   3.268   (   0.000    1.566    1.566)    2.215
   4.301   (   0.000  -14.888  -14.888)   21.055
   4.494   (  -0.000   -1.967   -1.967)    2.782
   4.736   (   0.000  -12.019  -12.019)   16.998
   5.189   (   0.000    0.762    0.762)    1.078
   5.525   (   0.000   13.214   13.214)   18.687
   5.569   (  -0.000   14.397   14.397)   20.360
   5.794   (   0.000    8.666    8.666)   12.255
   8.742   (   0.000   -3.802   -3.802)    5.377
  11.541   (   0.000  -27.701  -27.701)   39.175
  12.190   (   0.000  -10.036  -10.036)   14.193
  13.918   (   0.000   21.153   21.153)   29.915
======================= Grid point 197 (35/84) =======================
q-point: ( 0.17  0.17  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.769   (  16.896   16.896   10.535)   26.113
   1.818   (  16.888   16.888    7.320)   24.980
   3.219   (  27.150   27.150   15.831)   41.532
   3.824   (  -3.352   -3.352    6.988)    8.444
   4.091   (  -1.779   -1.779   -2.044)    3.241
   4.154   (   2.298    2.298   -7.728)    8.383
   4.186   (  -1.947   -1.947   -8.740)    9.163
   4.779   (   4.603    4.603   10.007)   11.938
   4.955   (   9.319    9.319   -1.263)   13.239
   5.554   (   1.248    1.248   13.250)   13.367
   5.830   (  11.358   11.358    0.334)   16.066
   8.187   (  -0.252   -0.252   -1.584)    1.624
  11.723   ( -18.488  -18.488   -0.892)   26.162
  12.189   (  -2.303   -2.303  -21.771)   22.014
  14.750   (   0.970    0.970    9.854)    9.949
======================= Grid point 198 (36/84) =======================
q-point: ( 0.25  0.17  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.135   (  18.386   13.651    7.535)   24.107
   2.185   (  18.750   11.926    5.995)   23.016
   3.633   (   4.629    2.787    2.039)    5.775
   3.816   (   6.182    7.820   11.973)   15.580
   3.991   (  -0.239    4.544    0.389)    4.567
   4.039   (  -6.168   -1.972   -6.017)    8.840
   4.400   (   9.655   -7.410  -12.904)   17.738
   4.925   (   9.737    2.360   11.756)   15.446
   5.146   (   9.332   11.057   -2.212)   14.637
   5.588   (   2.433    3.496   15.534)   16.107
   6.055   (  10.091   17.539    1.046)   20.262
   8.213   (   3.235   -2.603   -2.414)    4.803
  11.374   ( -15.555  -28.973   -3.016)   33.022
  12.146   (  -2.003   -2.840  -23.298)   23.556
  14.690   (  -6.080    8.110   12.214)   15.872
======================= Grid point 199 (37/84) =======================
q-point: ( 0.33  0.17  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.509   (  17.806   10.781    6.001)   21.664
   2.573   (  18.687    8.692    5.296)   21.279
   3.565   (  -6.292   -1.357    0.983)    6.511
   3.757   (  -8.685    7.218   -0.007)   11.293
   3.933   (  -3.645  -17.144   -4.339)   18.056
   4.210   (  12.612   10.086   10.392)   19.204
   4.537   (   4.261   -5.349  -18.330)   19.564
   5.113   (   6.543    0.853    7.694)   10.136
   5.362   (  12.041   10.104    5.255)   16.574
   5.662   (   5.001    9.319   14.315)   17.798
   6.213   (   4.439   21.010    1.856)   21.554
   8.334   (   8.519   -5.142   -3.268)   10.473
  11.083   ( -12.821  -36.855   -6.043)   39.486
  12.107   (  -1.703   -3.195  -23.912)   24.185
  14.535   (  -8.251   13.011   14.530)   21.177
======================= Grid point 200 (38/84) =======================
q-point: ( 0.42  0.17  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.853   (  15.363    9.050    6.373)   18.935
   2.939   (  16.767    6.767    3.476)   18.412
   3.410   (  -8.688    1.040    3.028)    9.260
   3.542   ( -11.978    8.019   -0.459)   14.422
   3.881   (  -1.579  -21.959   -3.065)   22.228
   4.403   (   6.202    2.791    8.338)   10.759
   4.594   (   1.600   -3.975  -22.753)   23.153
   5.190   (   1.772    1.213    4.100)    4.628
   5.574   (   6.689    2.631   18.862)   20.185
   5.813   (  10.010   20.387    6.047)   23.502
   6.219   (  -3.723   21.581    1.889)   21.981
   8.527   (   8.825   -6.709   -3.832)   11.729
  10.861   (  -8.273  -41.638   -8.357)   43.267
  12.079   (  -1.005   -3.458  -24.231)   24.497
  14.384   (  -5.727   15.795   16.195)   23.336
======================= Grid point 201 (39/84) =======================
q-point: (-0.50  0.17  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.076   (   0.000   10.043    8.779)   13.338
   3.209   (   0.000   -0.336   -0.468)    0.576
   3.286   (  -0.000   -0.654    2.502)    2.586
   3.292   (   0.000   15.773    2.424)   15.958
   3.865   (   0.000  -23.497   -2.897)   23.675
   4.466   (   0.000   -0.236    7.846)    7.849
   4.610   (   0.000   -2.869  -24.964)   25.129
   5.206   (   0.000    1.170    3.846)    4.020
   5.636   (   0.000    1.123   22.218)   22.247
   5.986   (  -0.000   25.368    5.426)   25.942
   6.119   (   0.000   19.455    0.080)   19.455
   8.629   (   0.000   -7.115   -3.994)    8.160
  10.771   (   0.000  -43.305   -9.249)   44.281
  12.068   (   0.000   -3.563  -24.331)   24.591
  14.323   (   0.000   16.696   16.813)   23.694
======================= Grid point 211 (40/84) =======================
q-point: ( 0.25  0.25  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.430   (  12.289   12.289    3.809)   17.792
   2.443   (  15.637   15.637    6.186)   22.963
   3.613   (  -2.848   -2.848   -0.855)    4.117
   3.766   ( -10.203  -10.203   -3.049)   14.748
   3.963   (   2.665    2.665    8.733)    9.511
   4.231   (   1.376    1.376  -13.647)   13.785
   4.413   (  11.044   11.044   -1.213)   15.666
   5.006   (   6.038    6.038   14.540)   16.862
   5.347   (   7.776    7.776    2.309)   11.236
   5.709   (   9.237    9.237   12.431)   18.032
   6.410   (  16.063   16.063   -0.809)   22.731
   8.175   (  -0.478   -0.478   -3.024)    3.099
  10.807   ( -25.698  -25.698   -0.977)   36.355
  12.095   (  -2.274   -2.274  -26.343)   26.539
  14.767   (  -0.132   -0.132   12.528)   12.529
======================= Grid point 212 (41/84) =======================
q-point: ( 0.33  0.25  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.686   (  12.664    3.216    2.639)   13.330
   2.763   (  15.181   12.303    5.686)   20.350
   3.529   (  -5.132   -2.344   -1.420)    5.819
   3.609   (  -5.255  -11.851   -1.570)   13.058
   3.892   (  -7.709    4.660    1.249)    9.094
   4.309   (   5.866   -1.221   -7.614)    9.689
   4.573   (   3.588    7.696   -6.547)   10.722
   5.150   (   6.338    3.270   13.469)   15.241
   5.478   (   5.911    2.544   14.758)   16.100
   5.963   (  14.604   18.397    4.244)   23.869
   6.665   (   7.700   22.193   -1.054)   23.514
   8.219   (   5.163   -5.652   -3.644)    8.478
  10.335   ( -20.181  -35.018   -1.679)   40.452
  12.052   (  -1.904   -2.425  -27.774)   27.945
  14.713   (  -4.362    4.574   13.676)   15.065
======================= Grid point 213 (42/84) =======================
q-point: ( 0.42  0.25  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.930   (  10.263   -9.539    1.321)   14.074
   3.051   (  12.692    9.126    5.308)   16.509
   3.424   (  -5.259    0.416    0.257)    5.281
   3.529   (  -2.400   -7.664    0.051)    8.031
   3.708   (  -9.402    6.391   -0.803)   11.397
   4.426   (   5.026   -1.031   -1.728)    5.414
   4.592   (  -0.909    3.241  -14.853)   15.230
   5.233   (   2.079    3.482   10.140)   10.921
   5.591   (   4.456   -0.056   21.673)   22.126
   6.271   (  14.716   22.585    0.714)   26.966
   6.704   (  -3.621   25.261   -0.773)   25.531
   8.364   (   7.536   -9.029   -4.049)   12.438
   9.997   ( -12.225  -41.381   -2.747)   43.236
  12.020   (  -1.131   -2.526  -28.617)   28.751
  14.621   (  -3.839    7.426   14.753)   16.957
======================= Grid point 214 (43/84) =======================
q-point: (-0.50  0.25  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.050   (   0.000  -20.277    0.257)   20.278
   3.253   (  -0.000    0.225    1.247)    1.267
   3.306   (   0.000   11.881    6.913)   13.746
   3.543   (   0.000   -4.880    1.333)    5.059
   3.553   (  -0.000    9.251   -2.751)    9.652
   4.491   (   0.000    2.315    9.804)   10.073
   4.568   (  -0.000   -1.505  -28.002)   28.042
   5.252   (   0.000    3.717    9.595)   10.290
   5.639   (   0.000   -0.385   22.773)   22.777
   6.522   (  -0.000   19.446   -0.792)   19.462
   6.564   (   0.000   29.869    0.116)   29.870
   8.455   (   0.000   -9.765   -4.159)   10.614
   9.865   (   0.000  -43.747   -3.222)   43.866
  12.008   (   0.000   -2.564  -28.914)   29.027
  14.579   (   0.000    8.372   15.192)   17.346
======================= Grid point 225 (44/84) =======================
q-point: ( 0.33  0.33  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.678   (  -3.325   -3.325    1.391)    4.904
   3.014   (  11.413   11.413    5.369)   17.011
   3.442   (  -2.281   -2.281   -0.936)    3.359
   3.515   (  -0.983   -0.983   -1.148)    1.803
   3.918   (  -2.690   -2.690    7.208)    8.150
   4.305   (   2.550    2.550  -21.323)   21.626
   4.680   (   1.489    1.489   -7.048)    7.356
   5.251   (   5.833    5.833   19.133)   20.835
   5.521   (   2.378    2.378   19.042)   19.336
   6.351   (  18.688   18.688    0.389)   26.432
   7.066   (  15.001   15.001   -2.209)   21.329
   8.131   (  -1.949   -1.949   -3.946)    4.814
   9.680   ( -28.188  -28.188   -0.674)   39.870
  12.006   (  -1.993   -1.993  -29.608)   29.742
  14.743   (  -0.947   -0.947   13.978)   14.042
======================= Grid point 226 (45/84) =======================
q-point: ( 0.42  0.33  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.596   (  -3.610  -18.867    0.132)   19.210
   3.205   (   6.555    5.305    4.611)    9.611
   3.464   (   3.619    3.015    0.409)    4.728
   3.512   (   0.495    3.167    1.392)    3.495
   3.809   (  -6.805    3.496    5.689)    9.534
   4.405   (   6.292   -0.144  -19.908)   20.880
   4.648   (  -3.954    0.940  -11.828)   12.507
   5.343   (   2.629    6.961   18.438)   19.883
   5.586   (   3.445   -0.270   22.101)   22.370
   6.706   (  15.333   19.640   -1.065)   24.939
   7.216   (  -0.375   23.188   -2.408)   23.316
   8.166   (   4.609   -9.996   -4.113)   11.751
   9.209   ( -17.007  -34.553   -0.678)   38.518
  11.973   (  -1.180   -2.039  -30.684)   30.774
  14.705   (  -2.091    1.292   14.571)   14.776
======================= Grid point 227 (46/84) =======================
q-point: (-0.50  0.33  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.556   (   0.000  -24.321   -0.304)   24.323
   3.270   (   0.000    0.962    2.713)    2.879
   3.519   (   0.000    7.153    3.960)    8.176
   3.537   (   0.000    1.724    2.203)    2.797
   3.706   (   0.000    7.572    3.567)    8.370
   4.534   (   0.000    1.277   -2.814)    3.090
   4.543   (  -0.000   -1.050  -29.172)   29.191
   5.367   (   0.000    7.348   17.695)   19.160
   5.627   (   0.000   -0.668   22.854)   22.864
   6.910   (  -0.000   18.281   -1.643)   18.355
   7.155   (   0.000   26.518   -2.256)   26.613
   8.233   (   0.000  -11.788   -4.139)   12.494
   9.024   (   0.000  -37.189   -0.782)   37.197
  11.960   (   0.000   -2.055  -31.065)   31.133
  14.680   (   0.000    2.063   14.853)   14.995
======================= Grid point 239 (47/84) =======================
q-point: ( 0.42  0.42  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.262   ( -12.250  -12.250   -0.377)   17.328
   3.272   (   1.606    1.606    4.309)    4.871
   3.525   (   3.049    3.049    1.397)    4.533
   3.581   (   2.084    2.084    1.063)    3.133
   3.862   (   1.806    1.806   15.190)   15.403
   4.438   (   3.394    3.394  -27.758)   28.169
   4.621   (  -3.263   -3.263  -18.082)   18.661
   5.481   (   5.177    5.177   22.101)   23.282
   5.590   (   1.038    1.038   22.346)   22.394
   7.069   (  15.332   15.332   -2.065)   21.781
   7.574   (   9.248    9.248   -3.047)   13.428
   7.992   (  -5.021   -5.021   -3.984)    8.143
   8.623   ( -21.419  -21.419    0.118)   30.291
  11.939   (  -1.200   -1.200  -31.872)   31.917
  14.702   (  -0.887   -0.887   14.855)   14.908
======================= Grid point 240 (48/84) =======================
q-point: (-0.50  0.42  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.132   (   0.000  -15.688   -0.542)   15.697
   3.286   (   0.000    0.565    4.133)    4.171
   3.574   (   0.000    1.524    1.877)    2.418
   3.586   (  -0.000    0.794    1.194)    1.434
   3.881   (   0.000    7.938   18.090)   19.755
   4.526   (   0.000   -0.570  -29.709)   29.714
   4.535   (  -0.000   -0.916  -20.153)   20.174
   5.527   (   0.000    7.180   22.093)   23.230
   5.615   (   0.000   -0.463   22.821)   22.826
   7.264   (  -0.000   15.962   -2.585)   16.170
   7.604   (   0.000   16.330   -3.094)   16.620
   7.979   (  -0.000  -12.688   -3.818)   13.250
   8.388   (   0.000  -23.635    0.314)   23.637
  11.926   (   0.000   -1.206  -32.290)   32.312
  14.690   (   0.000   -0.462   15.032)   15.039
======================= Grid point 259 (49/84) =======================
q-point: (-0.50 -0.50  0.08)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 168
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.963   (   0.000    0.000   -0.661)    0.661
   3.292   (   0.000    0.000    4.669)    4.669
   3.590   (   0.000    0.000    1.688)    1.688
   3.590   (   0.000    0.000    1.688)    1.688
   3.972   (   0.000    0.000   26.141)   26.141
   4.520   (   0.000    0.000  -29.863)   29.863
   4.520   (   0.000    0.000  -29.863)   29.863
   5.610   (   0.000    0.000   22.791)   22.791
   5.610   (   0.000    0.000   22.791)   22.791
   7.489   (  -0.000   -0.000   -3.349)    3.349
   7.788   (   0.000    0.000   -3.283)    3.283
   7.788   (   0.000    0.000   -3.283)    3.283
   8.122   (   0.000    0.000    0.660)    0.660
  11.914   (   0.000    0.000  -32.718)   32.718
  14.682   (   0.000    0.000   15.173)   15.173
======================= Grid point 366 (50/84) =======================
q-point: ( 0.17  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 189
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.027   (  13.421   13.421   13.421)   23.246
   2.027   (  13.421   13.421   13.421)   23.246
   3.586   (  14.620   14.620   14.620)   25.323
   4.005   (  -7.638   -7.638   -7.638)   13.230
   4.011   (   6.923    6.923    6.923)   11.990
   4.033   (  -0.647   -0.647   -0.647)    1.120
   4.033   (  -0.647   -0.647   -0.647)    1.120
   4.987   (   7.182    7.182    7.182)   12.440
   4.987   (   7.182    7.182    7.182)   12.440
   5.845   (   7.837    7.837    7.837)   13.574
   5.845   (   7.837    7.837    7.837)   13.574
   8.152   (  -1.433   -1.433   -1.433)    2.481
  11.697   ( -12.931  -12.931  -12.931)   22.397
  11.697   ( -12.931  -12.931  -12.931)   22.397
  14.907   (   4.222    4.222    4.222)    7.313
======================= Grid point 367 (51/84) =======================
q-point: ( 0.25  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.318   (  15.415   10.371   10.371)   21.278
   2.348   (  16.951    9.088    9.088)   21.272
   3.658   (  -3.629   -0.516   -0.516)    3.702
   3.860   (  -6.390   -9.899   -9.899)   15.388
   3.929   (  -5.779    1.233    1.233)    6.037
   4.127   (   6.897   -5.654   -5.654)   10.560
   4.299   (  13.891   17.612   17.612)   28.518
   5.153   (   6.497    6.554    6.554)   11.319
   5.175   (  13.847    8.512    8.512)   18.349
   5.885   (   2.637   11.052   11.052)   15.851
   6.127   (  12.466    7.933    7.933)   16.772
   8.149   (   1.706   -3.428   -3.428)    5.139
  11.228   ( -21.283  -14.180  -14.180)   29.243
  11.661   (  -1.866  -18.781  -18.781)   26.626
  14.918   (  -2.737    8.429    8.429)   12.231
======================= Grid point 368 (52/84) =======================
q-point: ( 0.33  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 476
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.645   (  16.110    7.349    7.349)   19.172
   2.707   (  17.566    6.605    6.605)   19.895
   3.566   (  -4.622   -1.204   -1.204)    4.925
   3.723   (  -7.794   -7.106   -7.106)   12.718
   3.815   (  -4.558   -3.012   -3.012)    6.238
   4.211   (   1.178  -12.077  -12.077)   17.120
   4.526   (   8.225   18.072   18.072)   26.849
   5.251   (   3.218    6.404    6.404)    9.611
   5.520   (  18.356   10.049   10.049)   23.215
   5.948   (   3.337   11.171   11.171)   16.147
   6.321   (   5.588   10.506   10.506)   15.874
   8.240   (   7.226   -5.654   -5.654)   10.777
  10.831   ( -17.239  -20.340  -20.340)   33.535
  11.624   (  -1.660  -18.962  -18.962)   26.868
  14.823   (  -5.775   11.851   11.851)   17.728
======================= Grid point 369 (53/84) =======================
q-point: ( 0.42  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.968   (  15.225    4.996    4.996)   16.785
   3.039   (  13.954    4.811    4.811)   15.524
   3.473   (  -4.557    1.701    1.701)    5.153
   3.528   ( -11.028   -1.383   -1.383)   11.200
   3.752   (  -1.788   -9.131   -9.131)   13.036
   4.196   (  -1.819  -13.693  -13.693)   19.449
   4.648   (   3.852   13.681   13.681)   19.727
   5.294   (   1.140    6.053    6.053)    8.636
   5.882   (  16.265   11.123   11.123)   22.627
   6.009   (   2.275   11.334   11.334)   16.190
   6.336   (  -4.049   11.681   11.681)   17.008
   8.414   (   8.271   -7.029   -7.029)   12.931
  10.537   ( -10.760  -24.468  -24.468)   36.237
  11.596   (  -0.978  -19.101  -19.101)   27.030
  14.711   (  -4.403   13.978   13.978)   20.252
======================= Grid point 370 (54/84) =======================
q-point: (-0.50  0.17  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 277
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.217   (   0.000    0.980    0.980)    1.387
   3.228   (  -0.000    5.680    5.680)    8.033
   3.312   (   0.000   -1.477   -1.477)    2.089
   3.399   (   0.000    7.155    7.155)   10.119
   3.735   (   0.000  -10.637  -10.637)   15.043
   4.171   (   0.000  -15.077  -15.077)   21.322
   4.687   (   0.000   12.076   12.076)   17.078
   5.305   (   0.000    5.895    5.895)    8.337
   6.033   (   0.000   11.399   11.399)   16.121
   6.156   (  -0.000    8.578    8.578)   12.132
   6.186   (   0.000   14.736   14.736)   20.840
   8.511   (   0.000   -7.372   -7.372)   10.426
  10.421   (   0.000  -25.972  -25.972)   36.731
  11.586   (   0.000  -19.153  -19.153)   27.086
  14.663   (   0.000   14.697   14.697)   20.785
======================= Grid point 380 (55/84) =======================
q-point: ( 0.25  0.25  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 474
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.540   (  10.508   10.508    6.791)   16.338
   2.570   (  13.274   13.274    5.836)   19.659
   3.608   (  -3.721   -3.721   -0.028)    5.262
   3.651   (  -9.489   -9.489   -8.503)   15.886
   3.968   (  -1.467   -1.467   -6.101)    6.444
   3.991   (  -2.367   -2.367   -8.820)    9.434
   4.672   (  13.412   13.412   24.042)   30.623
   5.295   (   8.309    8.309   12.951)   17.488
   5.404   (  13.112   13.112    3.337)   18.842
   6.021   (   4.659    4.659   17.076)   18.304
   6.390   (  14.892   14.892   -0.918)   21.080
   8.092   (  -1.517   -1.517   -4.715)    5.180
  10.779   ( -25.658  -25.658   -1.725)   36.327
  11.447   (  -5.237   -5.237  -34.594)   35.378
  15.016   (   1.337    1.337   10.098)   10.274
======================= Grid point 381 (56/84) =======================
q-point: ( 0.33  0.25  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.762   (  10.910    1.929    4.625)   12.006
   2.861   (  14.579    9.851    3.022)   17.853
   3.480   (  -6.386   -8.746   -6.126)   12.442
   3.545   (  -2.667   -3.671   -2.540)    5.200
   3.836   (  -7.322    2.623   -4.510)    8.991
   4.027   (   1.148   -6.974  -14.624)   16.242
   4.828   (   3.349    9.065   24.695)   26.519
   5.420   (   3.711   10.462   12.231)   16.518
   5.733   (  16.844    9.449    9.189)   21.388
   6.130   (   6.430    9.120   12.294)   16.603
   6.643   (   8.156   18.773   -0.728)   20.481
   8.115   (   4.199   -6.143   -6.273)    9.732
  10.277   ( -22.291  -30.860   -4.325)   38.314
  11.370   (  -2.787   -7.547  -35.948)   36.837
  14.993   (  -2.860    4.803   11.993)   13.232
======================= Grid point 382 (57/84) =======================
q-point: ( 0.42  0.25  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.964   (   7.906  -11.121    1.657)   13.745
   3.137   (  11.407    6.768    2.164)   13.439
   3.387   (  -2.701   -3.699   -5.145)    6.888
   3.515   (  -0.119    0.336   -1.015)    1.076
   3.681   (  -7.409    1.997   -0.960)    7.733
   4.021   (  -1.520   -5.402  -21.705)   22.419
   4.860   (   0.500    5.209   24.865)   25.410
   5.465   (   1.107   11.045   11.460)   15.955
   6.001   (   8.086    1.601   17.791)   19.608
   6.305   (  10.842   17.682    3.226)   20.991
   6.694   (  -2.878   21.342    0.430)   21.539
   8.246   (   7.135   -9.227   -7.331)   13.776
   9.899   ( -13.788  -35.247   -7.423)   38.569
  11.327   (  -1.456   -8.549  -35.936)   36.968
  14.927   (  -2.914    7.012   13.385)   15.389
======================= Grid point 383 (58/84) =======================
q-point: (-0.50  0.25  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.050   (   0.000  -21.802   -0.802)   21.817
   3.277   (   0.000    2.413    0.607)    2.488
   3.359   (   0.000    7.862   -4.018)    8.829
   3.524   (   0.000    4.609    2.217)    5.114
   3.581   (   0.000   -1.373    2.343)    2.716
   3.996   (   0.000   -4.439  -26.696)   27.063
   4.864   (   0.000    3.554   24.942)   25.194
   5.475   (   0.000   11.104   11.211)   15.779
   6.077   (   0.000   -0.551   19.953)   19.961
   6.498   (  -0.000   18.373   -1.478)   18.433
   6.588   (   0.000   23.615    2.886)   23.791
   8.333   (   0.000   -9.875   -7.581)   12.450
   9.750   (   0.000  -36.996   -8.681)   38.001
  11.311   (   0.000   -8.860  -35.876)   36.954
  14.894   (   0.000    7.763   13.895)   15.916
======================= Grid point 394 (59/84) =======================
q-point: ( 0.33  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 476
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.716   (  -5.503   -5.503    2.235)    8.097
   3.087   (  11.343   11.343    0.140)   16.043
   3.384   (   0.200    0.200   -6.197)    6.204
   3.502   (  -0.850   -0.850    0.102)    1.206
   3.881   (  -2.661   -2.661   -8.736)    9.512
   3.902   (  -1.604   -1.604  -16.019)   16.179
   4.904   (  -0.132   -0.132   26.609)   26.609
   5.622   (   7.449    7.449   16.696)   19.742
   5.889   (   6.419    6.419   15.695)   18.131
   6.378   (  14.498   14.498    2.938)   20.712
   7.002   (  14.105   14.105   -3.876)   20.320
   8.017   (  -2.367   -2.367   -6.883)    7.654
   9.661   ( -27.776  -27.776   -1.189)   39.299
  11.259   (  -3.935   -3.935  -40.909)   41.286
  15.034   (  -0.232   -0.232   12.757)   12.761
======================= Grid point 395 (60/84) =======================
q-point: ( 0.42  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.598   (  -4.822  -19.830    0.013)   20.408
   3.285   (   6.531    6.468    1.443)    9.304
   3.438   (   4.466    6.000   -5.169)    9.092
   3.532   (   3.751    1.766    0.032)    4.147
   3.744   (  -7.505    3.367   -7.683)   11.255
   3.927   (   0.760   -4.041  -22.753)   23.122
   4.872   (  -2.108   -2.795   28.571)   28.784
   5.710   (   1.876   11.571   16.815)   20.498
   6.008   (   4.734   -0.199   19.003)   19.585
   6.678   (  13.707   17.710   -1.480)   22.444
   7.147   (  -0.135   21.321   -4.069)   21.706
   8.046   (   4.451   -9.923   -7.382)   13.144
   9.187   ( -17.355  -32.498   -1.479)   36.871
  11.199   (  -1.969   -4.531  -42.326)   42.613
  15.010   (  -1.529    1.532   13.516)   13.688
======================= Grid point 396 (61/84) =======================
q-point: (-0.50  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 506
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.546   (   0.000  -24.549   -0.736)   24.560
   3.339   (   0.000    3.209    3.347)    4.637
   3.554   (   0.000   10.872  -14.289)   17.955
   3.593   (   0.000    2.296    1.009)    2.508
   3.600   (  -0.000    1.696    3.781)    4.144
   3.932   (   0.000   -2.240  -29.592)   29.676
   4.846   (   0.000   -4.106   29.937)   30.218
   5.726   (   0.000   12.243   16.659)   20.674
   6.060   (   0.000   -0.939   19.374)   19.397
   6.862   (  -0.000   17.005   -2.908)   17.252
   7.094   (   0.000   24.068   -3.505)   24.322
   8.112   (   0.000  -11.552   -7.489)   13.767
   8.998   (   0.000  -34.553   -1.865)   34.603
  11.178   (   0.000   -4.722  -42.655)   42.916
  14.991   (   0.000    2.148   13.828)   13.994
======================= Grid point 408 (62/84) =======================
q-point: ( 0.42  0.42  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 474
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.251   ( -12.677  -12.677   -0.699)   17.942
   3.376   (   2.403    2.403    4.853)    5.925
   3.547   (   4.306    4.306   -1.360)    6.240
   3.609   (   4.347    4.347   -0.267)    6.153
   3.792   (   0.033    0.033  -18.717)   18.717
   3.872   (   0.168    0.168  -25.050)   25.051
   4.790   (  -4.287   -4.287   31.916)   32.487
   5.898   (   5.767    5.767   18.492)   20.211
   6.013   (   1.251    1.251   18.961)   19.044
   7.009   (  14.122   14.122   -3.564)   20.287
   7.485   (   8.906    8.906   -5.561)   13.768
   7.876   (  -4.735   -4.735   -7.131)    9.782
   8.626   ( -20.831  -20.831    0.192)   29.460
  11.130   (  -2.242   -2.242  -44.555)   44.667
  15.014   (  -0.567   -0.567   13.876)   13.900
======================= Grid point 409 (63/84) =======================
q-point: (-0.50  0.42  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.116   (   0.000  -15.961   -0.964)   15.990
   3.397   (   0.000    2.247    6.105)    6.506
   3.625   (   0.000    0.984    1.478)    1.775
   3.627   (  -0.000    0.853    3.018)    3.136
   3.776   (  -0.000    9.632  -24.239)   26.082
   3.899   (   0.000   -1.024  -30.684)   30.701
   4.737   (   0.000   -5.528   34.791)   35.227
   5.944   (   0.000    8.338   18.565)   20.352
   6.043   (   0.000   -0.644   18.930)   18.941
   7.190   (  -0.000   14.639   -4.539)   15.326
   7.513   (   0.000   15.729   -5.594)   16.695
   7.867   (  -0.000  -11.983   -6.890)   13.823
   8.397   (   0.000  -22.693    0.508)   22.698
  11.107   (   0.000   -2.328  -45.135)   45.195
  15.005   (   0.000   -0.224   14.056)   14.058
======================= Grid point 428 (64/84) =======================
q-point: (-0.50 -0.50  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 172
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.944   (   0.000    0.000   -1.156)    1.156
   3.422   (   0.000    0.000    7.685)    7.685
   3.635   (   0.000    0.000    1.831)    1.831
   3.635   (   0.000    0.000    1.831)    1.831
   3.888   (   0.000    0.000  -30.734)   30.734
   3.888   (   0.000    0.000  -30.734)   30.734
   4.669   (   0.000    0.000   39.481)   39.481
   6.036   (   0.000    0.000   18.767)   18.767
   6.036   (   0.000    0.000   18.767)   18.767
   7.393   (  -0.000   -0.000   -5.823)    5.823
   7.691   (   0.000    0.000   -5.995)    5.995
   7.691   (   0.000    0.000   -5.995)    5.995
   8.141   (   0.000    0.000    1.141)    1.141
  11.083   (   0.000    0.000  -45.794)   45.794
  15.000   (   0.000    0.000   14.189)   14.189
======================= Grid point 549 (65/84) =======================
q-point: ( 0.25  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 182
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.697   (   7.769    7.769    7.769)   13.457
   2.697   (   7.769    7.769    7.769)   13.457
   3.446   ( -10.334  -10.334  -10.334)   17.898
   3.589   (  -2.192   -2.192   -2.192)    3.797
   3.859   (  -4.180   -4.180   -4.180)    7.240
   3.859   (  -4.180   -4.180   -4.180)    7.240
   5.170   (  18.958   18.958   18.958)   32.837
   5.520   (   9.842    9.842    9.842)   17.047
   5.520   (   9.842    9.842    9.842)   17.047
   6.380   (   8.666    8.666    8.666)   15.010
   6.380   (   8.666    8.666    8.666)   15.010
   8.000   (  -3.598   -3.598   -3.598)    6.232
  10.740   ( -17.751  -17.751  -17.751)   30.746
  10.740   ( -17.751  -17.751  -17.751)   30.746
  15.155   (   3.515    3.515    3.515)    6.089
======================= Grid point 550 (66/84) =======================
q-point: ( 0.33  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 474
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.858   (   6.802    2.331    2.331)    7.559
   2.888   (  11.150    1.597    1.597)   11.377
   3.267   (  -6.303  -11.328  -11.328)   17.215
   3.532   (  -2.638    0.030    0.030)    2.638
   3.766   (  -4.415   -2.394   -2.394)    5.564
   3.810   (  -2.407   -6.504   -6.504)    9.508
   5.347   (   2.349   20.583   20.583)   29.203
   5.650   (   4.085   12.169   12.169)   17.687
   5.893   (  20.879    7.725    7.725)   23.565
   6.398   (   1.537   10.744   10.744)   15.272
   6.679   (  10.480    6.486    6.486)   13.927
   7.976   (   1.885   -6.860   -6.860)    9.882
  10.108   ( -27.495  -14.900  -14.900)   34.641
  10.699   (  -1.826  -25.545  -25.545)   36.172
  15.175   (  -1.025    5.630    5.630)    8.028
======================= Grid point 551 (67/84) =======================
q-point: ( 0.42  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.942   (   1.116   -8.178   -8.178)   11.620
   3.125   (  10.599   -2.864   -2.864)   11.347
   3.215   (   1.027   -7.014   -7.014)    9.972
   3.505   (   0.689    0.317    0.317)    0.822
   3.679   (  -3.765    0.443    0.443)    3.817
   3.731   (  -5.621   -5.826   -5.826)    9.974
   5.364   (  -0.009   19.176   19.176)   27.120
   5.700   (   1.152   13.105   13.105)   18.569
   6.275   (  15.613    7.935    7.935)   19.227
   6.433   (   1.517    9.294    9.294)   13.232
   6.773   (  -0.977    9.028    9.028)   12.805
   8.073   (   6.228   -9.373   -9.373)   14.646
   9.645   ( -16.782  -18.780  -18.780)   31.417
  10.668   (  -1.078  -25.501  -25.501)   36.080
  15.139   (  -1.884    7.105    7.105)   10.223
======================= Grid point 552 (68/84) =======================
q-point: (-0.50  0.25  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 276
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.944   (   0.000  -13.208  -13.208)   18.679
   3.242   (   0.000   -5.236   -5.236)    7.406
   3.246   (   0.000   -3.364   -3.364)    4.758
   3.587   (  -0.000   -9.432   -9.432)   13.338
   3.592   (   0.000    3.440    3.440)    4.866
   3.634   (   0.000    2.490    2.490)    3.521
   5.362   (   0.000   18.693   18.693)   26.435
   5.710   (   0.000   13.326   13.326)   18.846
   6.450   (   0.000    8.849    8.849)   12.514
   6.473   (   0.000    6.694    6.694)    9.467
   6.720   (   0.000   10.827   10.827)   15.312
   8.152   (   0.000   -9.896   -9.896)   13.995
   9.464   (   0.000  -20.332  -20.332)   28.753
  10.657   (   0.000  -25.485  -25.485)   36.041
  15.117   (   0.000    7.631    7.631)   10.792
======================= Grid point 563 (69/84) =======================
q-point: ( 0.33  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.764   (  -8.708   -8.708    2.220)   12.513
   2.990   (   5.985    5.985   -8.530)   12.017
   3.137   (  -0.615   -0.615  -16.669)   16.692
   3.529   (  -0.322   -0.322    2.343)    2.387
   3.717   (  -2.191   -2.191   -2.841)    4.203
   3.740   (  -1.630   -1.630   -3.682)    4.344
   5.544   (   2.236    2.236   32.238)   32.393
   5.926   (   8.946    8.946   12.856)   18.037
   6.106   (  14.368   14.368    4.865)   20.894
   6.527   (   4.597    4.597   12.405)   14.006
   6.914   (  12.519   12.519   -4.487)   18.265
   7.862   (  -3.291   -3.291   -7.909)    9.177
   9.633   ( -27.202  -27.202   -1.412)   38.495
  10.401   (  -6.712   -6.712  -41.498)   42.570
  15.235   (   0.713    0.713    6.672)    6.748
======================= Grid point 564 (70/84) =======================
q-point: ( 0.42  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.592   (  -6.596  -20.253   -0.841)   21.317
   3.117   (   5.649    1.569  -16.637)   17.640
   3.173   (   2.548    0.400  -20.428)   20.590
   3.527   (   0.354    1.535    2.461)    2.922
   3.683   (  -1.412    0.431    0.568)    1.581
   3.704   (  -1.909    0.638    0.230)    2.026
   5.550   (  -0.652    1.941   33.688)   33.751
   6.019   (   1.916   14.856   13.202)   19.967
   6.354   (   7.476    1.044   13.967)   15.876
   6.664   (   8.804   12.992    1.141)   15.736
   7.061   (   0.924   16.875   -3.796)   17.321
   7.873   (   3.909   -9.714   -9.238)   13.963
   9.142   ( -18.617  -27.613   -3.212)   33.457
  10.319   (  -2.218  -10.071  -41.786)   43.039
  15.227   (  -0.918    1.934    7.511)    7.810
======================= Grid point 565 (71/84) =======================
q-point: (-0.50  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.522   (   0.000  -23.645   -1.915)   23.722
   3.177   (   0.000   -0.976  -18.628)   18.654
   3.209   (   0.000   -0.445  -23.860)   23.864
   3.538   (  -0.000    1.208   -0.485)    1.302
   3.658   (   0.000    2.904    3.506)    4.553
   3.683   (   0.000    2.072    4.115)    4.607
   5.540   (   0.000    1.595   34.122)   34.160
   6.036   (   0.000   15.558   13.495)   20.595
   6.423   (   0.000   -1.800   16.099)   16.199
   6.796   (  -0.000   15.301   -3.403)   15.675
   7.030   (   0.000   18.110   -2.089)   18.230
   7.935   (   0.000  -11.122   -9.584)   14.682
   8.939   (   0.000  -28.467   -4.435)   28.810
  10.297   (   0.000  -10.874  -41.453)   42.856
  15.215   (   0.000    2.373    7.832)    8.183
======================= Grid point 577 (72/84) =======================
q-point: ( 0.42  0.42  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.234   ( -13.164  -13.164   -0.918)   18.640
   3.164   (   2.428    2.428  -25.652)   25.880
   3.190   (   1.021    1.021  -27.435)   27.473
   3.558   (   1.268    1.268    5.286)    5.582
   3.699   (   0.838    0.838    2.446)    2.718
   3.713   (   0.205    0.205    3.211)    3.224
   5.544   (  -1.303   -1.303   37.833)   37.878
   6.239   (   6.107    6.107   14.777)   17.116
   6.367   (   1.135    1.135   15.495)   15.578
   6.931   (  12.044   12.044   -3.757)   17.442
   7.353   (   8.350    8.350   -7.129)   13.794
   7.709   (  -4.400   -4.400   -8.848)   10.818
   8.631   ( -20.017  -20.017    0.203)   28.308
  10.189   (  -3.579   -3.579  -45.840)   46.119
  15.239   (  -0.230   -0.230    7.943)    7.949
======================= Grid point 578 (73/84) =======================
q-point: (-0.50  0.42  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.094   (   0.000  -16.236   -1.177)   16.279
   3.189   (  -0.000    1.418  -28.456)   28.492
   3.207   (   0.000    0.054  -30.407)   30.407
   3.572   (  -0.000    1.580    5.961)    6.167
   3.707   (   0.000    1.811    3.723)    4.140
   3.716   (  -0.000    1.072    4.142)    4.279
   5.526   (   0.000   -1.652   39.046)   39.080
   6.287   (   0.000    8.714   14.995)   17.343
   6.392   (   0.000   -1.115   15.200)   15.240
   7.086   (  -0.000   12.807   -5.318)   13.868
   7.382   (   0.000   14.423   -6.968)   16.018
   7.705   (  -0.000  -10.955   -8.728)   14.007
   8.407   (   0.000  -21.110    0.450)   21.115
  10.154   (   0.000   -4.011  -46.383)   46.556
  15.234   (   0.000    0.019    8.127)    8.127
======================= Grid point 597 (74/84) =======================
q-point: (-0.50 -0.50  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 168
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.918   (   0.000    0.000   -1.354)    1.354
   3.208   (   0.000    0.000  -32.050)   32.050
   3.208   (   0.000    0.000  -32.050)   32.050
   3.590   (  -0.000   -0.000    8.315)    8.315
   3.727   (   0.000    0.000    3.733)    3.733
   3.727   (   0.000    0.000    3.733)    3.733
   5.505   (   0.000    0.000   40.508)   40.508
   6.380   (   0.000    0.000   14.884)   14.884
   6.380   (   0.000    0.000   14.884)   14.884
   7.261   (  -0.000   -0.000   -6.762)    6.762
   7.549   (   0.000    0.000   -7.704)    7.704
   7.549   (   0.000    0.000   -7.704)    7.704
   8.167   (   0.000    0.000    1.313)    1.313
  10.114   (   0.000    0.000  -47.236)   47.236
  15.232   (   0.000    0.000    8.265)    8.265
======================= Grid point 732 (75/84) =======================
q-point: ( 0.33  0.33  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 189
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.772   ( -11.166  -11.166  -11.166)   19.339
   2.819   (  -4.936   -4.936   -4.936)    8.549
   2.819   (  -4.936   -4.936   -4.936)    8.549
   3.577   (   1.939    1.939    1.939)    3.359
   3.713   (   0.745    0.745    0.745)    1.290
   3.713   (   0.745    0.745    0.745)    1.290
   6.129   (  12.198   12.198   12.198)   21.128
   6.144   (  10.062   10.062   10.062)   17.427
   6.144   (  10.062   10.062   10.062)   17.427
   6.828   (   5.516    5.516    5.516)    9.554
   6.828   (   5.516    5.516    5.516)    9.554
   7.716   (  -5.440   -5.440   -5.440)    9.422
   9.605   ( -18.163  -18.163  -18.163)   31.460
   9.605   ( -18.163  -18.163  -18.163)   31.460
  15.316   (   1.689    1.689    1.689)    2.925
======================= Grid point 733 (76/84) =======================
q-point: ( 0.42  0.33  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 468
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.515   (  -9.407  -10.357  -10.357)   17.408
   2.738   (  -1.611  -12.602  -12.602)   17.894
   2.833   (   0.785  -12.650  -12.650)   17.907
   3.605   (   0.902    4.095    4.095)    5.861
   3.724   (   0.512    2.591    2.591)    3.700
   3.743   (   1.357    2.369    2.369)    3.615
   6.116   (  -1.079   16.841   16.841)   23.842
   6.267   (   3.739   13.327   13.327)   19.214
   6.535   (  14.408    4.935    4.935)   16.010
   6.771   (  -1.796    7.303    7.303)   10.483
   7.037   (   4.114    3.506    3.506)    6.443
   7.680   (   2.174   -9.237   -9.237)   13.242
   9.010   ( -21.773  -12.194  -12.194)   27.774
   9.584   (  -0.743  -26.178  -26.178)   37.029
  15.325   (  -0.327    2.429    2.429)    3.451
======================= Grid point 734 (77/84) =======================
q-point: (-0.50  0.33  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 277
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.416   (   0.000  -11.825  -11.825)   16.723
   2.727   (   0.000  -16.020  -16.020)   22.656
   2.841   (   0.000  -13.518  -13.518)   19.117
   3.615   (   0.000    4.883    4.883)    6.906
   3.729   (   0.000    3.287    3.287)    4.648
   3.758   (   0.000    3.010    3.010)    4.257
   6.104   (   0.000   16.719   16.719)   23.645
   6.304   (   0.000   14.535   14.535)   20.555
   6.699   (   0.000    4.442    4.442)    6.282
   6.753   (   0.000    5.928    5.928)    8.383
   7.049   (   0.000    5.134    5.134)    7.261
   7.728   (   0.000  -10.383  -10.383)   14.684
   8.773   (   0.000  -13.323  -13.323)   18.841
   9.577   (   0.000  -26.016  -26.016)   36.792
  15.319   (   0.000    2.706    2.706)    3.827
======================= Grid point 746 (78/84) =======================
q-point: ( 0.42  0.42  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 474
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.213   ( -13.283  -13.283   -1.239)   18.826
   2.643   (  -3.296   -3.296  -24.558)   24.997
   2.648   (  -3.451   -3.451  -22.929)   23.442
   3.677   (   2.644    2.644    5.494)    6.645
   3.770   (   1.601    1.601    3.671)    4.313
   3.784   (   1.479    1.479    3.223)    3.843
   6.250   (   1.416    1.416   29.521)   29.589
   6.506   (   6.519    6.519   11.365)   14.634
   6.637   (   3.908    3.908    9.557)   11.040
   6.886   (   6.029    6.029    1.231)    8.615
   7.200   (   7.557    7.557   -7.609)   13.119
   7.525   (  -4.295   -4.295   -8.801)   10.694
   8.634   ( -19.189  -19.189    0.158)   27.138
   9.306   (  -4.944   -4.944  -39.121)   39.741
  15.346   (   0.111    0.111    2.852)    2.857
======================= Grid point 747 (79/84) =======================
q-point: (-0.50  0.42  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 504
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.069   (   0.000  -16.029   -1.278)   16.080
   2.590   (   0.000   -1.945  -28.205)   28.272
   2.637   (   0.000   -6.390  -24.464)   25.285
   3.706   (  -0.000    3.637    6.102)    7.104
   3.785   (   0.000    2.043    3.645)    4.178
   3.800   (   0.000    1.032    3.778)    3.917
   6.258   (   0.000    2.337   30.668)   30.757
   6.566   (   0.000    8.447   12.632)   15.196
   6.669   (   0.000   -1.776   11.829)   11.962
   6.982   (  -0.000   10.783   -4.502)   11.685
   7.239   (   0.000   11.295   -6.573)   13.069
   7.519   (  -0.000   -9.722   -9.238)   13.412
   8.411   (   0.000  -18.246   -0.266)   18.248
   9.264   (   0.000   -6.891  -39.027)   39.631
  15.345   (   0.000    0.274    3.022)    3.034
======================= Grid point 766 (80/84) =======================
q-point: (-0.50 -0.50  0.33)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 172
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.891   (   0.000    0.000   -1.189)    1.189
   2.573   (   0.000    0.000  -28.959)   28.959
   2.573   (   0.000    0.000  -28.959)   28.959
   3.748   (  -0.000   -0.000    6.776)    6.776
   3.807   (   0.000    0.000    3.780)    3.780
   3.807   (   0.000    0.000    3.780)    3.780
   6.275   (   0.000    0.000   33.599)   33.599
   6.650   (   0.000    0.000   11.526)   11.526
   6.650   (   0.000    0.000   11.526)   11.526
   7.128   (  -0.000   -0.000   -5.889)    5.889
   7.384   (   0.000    0.000   -8.272)    8.272
   7.384   (   0.000    0.000   -8.272)    8.272
   8.192   (   0.000    0.000    1.134)    1.134
   9.203   (   0.000    0.000  -40.421)   40.421
  15.345   (   0.000    0.000    3.153)    3.153
======================= Grid point 915 (81/84) =======================
q-point: ( 0.42  0.42  0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 182
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.146   (  -8.368   -8.368   -8.368)   14.494
   2.269   ( -10.509  -10.509  -10.509)   18.202
   2.269   ( -10.509  -10.509  -10.509)   18.202
   3.771   (   3.492    3.492    3.492)    6.049
   3.834   (   1.953    1.953    1.953)    3.382
   3.834   (   1.953    1.953    1.953)    3.382
   6.697   (   6.724    6.724    6.724)   11.646
   6.709   (   7.827    7.827    7.827)   13.557
   6.709   (   7.827    7.827    7.827)   13.557
   7.047   (   1.795    1.795    1.795)    3.110
   7.047   (   1.795    1.795    1.795)    3.110
   7.368   (  -5.561   -5.561   -5.561)    9.632
   8.637   ( -12.356  -12.356  -12.356)   21.401
   8.637   ( -12.356  -12.356  -12.356)   21.401
  15.374   (   0.382    0.382    0.382)    0.662
======================= Grid point 916 (82/84) =======================
q-point: (-0.50  0.42  0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 276
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.017   (   0.000   -7.175   -7.175)   10.147
   2.129   (   0.000  -11.874  -11.874)   16.792
   2.225   (  -0.000  -14.741  -14.741)   20.847
   3.810   (  -0.000    3.938    3.938)    5.569
   3.848   (   0.000    2.262    2.262)    3.199
   3.859   (   0.000    1.775    1.775)    2.511
   6.691   (   0.000   11.535   11.535)   16.313
   6.834   (  -0.000   10.692   10.692)   15.120
   6.857   (   0.000    3.602    3.602)    5.094
   6.935   (   0.000    2.932    2.932)    4.147
   7.150   (   0.000    0.055    0.055)    0.077
   7.335   (  -0.000   -8.433   -8.433)   11.926
   8.363   (   0.000   -6.840   -6.840)    9.673
   8.638   (   0.000  -18.349  -18.349)   25.950
  15.376   (   0.000    0.488    0.488)    0.690
======================= Grid point 935 (83/84) =======================
q-point: (-0.50 -0.50  0.42)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 168
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.871   (   0.000    0.000   -0.694)    0.694
   2.068   (  -0.000   -0.000  -18.840)   18.840
   2.068   (   0.000    0.000  -18.840)   18.840
   3.858   (  -0.000   -0.000    3.749)    3.749
   3.872   (   0.000    0.000    2.343)    2.343
   3.872   (   0.000    0.000    2.343)    2.343
   6.848   (  -0.000   -0.000   21.439)   21.439
   6.858   (   0.000    0.000    8.971)    8.971
   6.858   (   0.000    0.000    8.971)    8.971
   7.031   (  -0.000   -0.000   -3.414)    3.414
   7.216   (   0.000    0.000   -7.917)    7.917
   7.216   (   0.000    0.000   -7.917)    7.917
   8.211   (   0.000    0.000    0.653)    0.653
   8.509   (   0.000    0.000  -25.999)   25.999
  15.379   (   0.000    0.000    0.573)    0.573
======================= Grid point 1194 (84/84) =======================
q-point: (-0.50 -0.50 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.64e-04 5.98e-04 
Number of triplets: 44
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.864   (   0.000    0.000    0.000)    0.000
   1.864   (   0.000    0.000    0.000)    0.000
   1.864   (   0.000    0.000    0.000)    0.000
   3.897   (   0.000    0.000    0.000)    0.000
   3.897   (   0.000    0.000    0.000)    0.000
   3.897   (   0.000    0.000    0.000)    0.000
   6.995   (  -0.000   -0.000   -0.000)    0.000
   6.995   (  -0.000   -0.000   -0.000)    0.000
   6.995   (  -0.000   -0.000   -0.000)    0.000
   7.090   (   0.000    0.000    0.000)    0.000
   7.090   (   0.000    0.000    0.000)    0.000
   7.090   (   0.000    0.000    0.000)    0.000
   8.218   (   0.000    0.000    0.000)    0.000
   8.218   (   0.000    0.000    0.000)    0.000
  15.383   (   0.000    0.000    0.000)    0.000
=================== End of collection of collisions ===================
----------- Thermal conductivity (W/m-k) with tetrahedron method -----------
#  T(K)        xx         yy         zz         yz         xz         xy        #ipm
    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/25920
   10.0   1821.043   1821.043   1821.043     -0.000     -0.000      0.000 3/25920
   20.0    197.534    197.534    197.534     -0.000      0.000      0.000 3/25920
   30.0     67.887     67.887     67.887     -0.000      0.000      0.000 3/25920
   40.0     37.779     37.779     37.779     -0.000      0.000      0.000 3/25920
   50.0     25.992     25.992     25.992     -0.000      0.000      0.000 3/25920
   60.0     19.938     19.938     19.938     -0.000      0.000      0.000 3/25920
   70.0     16.275     16.275     16.275     -0.000      0.000      0.000 3/25920
   80.0     13.812     13.812     13.812     -0.000      0.000      0.000 3/25920
   90.0     12.035     12.035     12.035     -0.000      0.000      0.000 3/25920
  100.0     10.686     10.686     10.686     -0.000      0.000      0.000 3/25920
  110.0      9.623      9.623      9.623     -0.000      0.000      0.000 3/25920
  120.0      8.761      8.761      8.761     -0.000      0.000      0.000 3/25920
  130.0      8.047      8.047      8.047     -0.000      0.000      0.000 3/25920
  140.0      7.444      7.444      7.444     -0.000      0.000      0.000 3/25920
  150.0      6.928      6.928      6.928     -0.000      0.000      0.000 3/25920
  160.0      6.480      6.480      6.480     -0.000      0.000      0.000 3/25920
  170.0      6.088      6.088      6.088     -0.000      0.000      0.000 3/25920
  180.0      5.742      5.742      5.742     -0.000      0.000      0.000 3/25920
  190.0      5.434      5.434      5.434     -0.000      0.000      0.000 3/25920
  200.0      5.158      5.158      5.158     -0.000      0.000      0.000 3/25920
  210.0      4.909      4.909      4.909     -0.000      0.000      0.000 3/25920
  220.0      4.683      4.683      4.683     -0.000      0.000      0.000 3/25920
  230.0      4.477      4.477      4.477     -0.000      0.000      0.000 3/25920
  240.0      4.289      4.289      4.289     -0.000      0.000      0.000 3/25920
  250.0      4.116      4.116      4.116     -0.000      0.000      0.000 3/25920
  260.0      3.957      3.957      3.957     -0.000      0.000      0.000 3/25920
  270.0      3.810      3.810      3.810     -0.000      0.000      0.000 3/25920
  280.0      3.673      3.673      3.673     -0.000      0.000      0.000 3/25920
  290.0      3.546      3.546      3.546     -0.000      0.000      0.000 3/25920
  300.0      3.427      3.427      3.427     -0.000      0.000      0.000 3/25920
  310.0      3.316      3.316      3.316     -0.000      0.000      0.000 3/25920
  320.0      3.213      3.213      3.213     -0.000      0.000      0.000 3/25920
  330.0      3.115      3.115      3.115     -0.000      0.000      0.000 3/25920
  340.0      3.023      3.023      3.023     -0.000      0.000      0.000 3/25920
  350.0      2.937      2.937      2.937     -0.000      0.000      0.000 3/25920
  360.0      2.855      2.855      2.855     -0.000      0.000      0.000 3/25920
  370.0      2.778      2.778      2.778     -0.000      0.000      0.000 3/25920
  380.0      2.705      2.705      2.705     -0.000      0.000      0.000 3/25920
  390.0      2.635      2.635      2.635     -0.000      0.000      0.000 3/25920
  400.0      2.569      2.569      2.569     -0.000      0.000      0.000 3/25920
  410.0      2.507      2.507      2.507     -0.000      0.000      0.000 3/25920
  420.0      2.447      2.447      2.447     -0.000      0.000      0.000 3/25920
  430.0      2.390      2.390      2.390     -0.000      0.000      0.000 3/25920
  440.0      2.336      2.336      2.336     -0.000      0.000      0.000 3/25920
  450.0      2.284      2.284      2.284     -0.000      0.000      0.000 3/25920
  460.0      2.234      2.234      2.234     -0.000      0.000      0.000 3/25920
  470.0      2.187      2.187      2.187     -0.000      0.000      0.000 3/25920
  480.0      2.141      2.141      2.141     -0.000      0.000      0.000 3/25920
  490.0      2.098      2.098      2.098     -0.000      0.000      0.000 3/25920
  500.0      2.056      2.056      2.056     -0.000      0.000      0.000 3/25920
  510.0      2.015      2.015      2.015     -0.000      0.000      0.000 3/25920
  520.0      1.977      1.977      1.977     -0.000      0.000      0.000 3/25920
  530.0      1.940      1.940      1.940     -0.000      0.000      0.000 3/25920
  540.0      1.904      1.904      1.904     -0.000      0.000      0.000 3/25920
  550.0      1.869      1.869      1.869     -0.000      0.000      0.000 3/25920
  560.0      1.836      1.836      1.836     -0.000      0.000      0.000 3/25920
  570.0      1.804      1.804      1.804     -0.000      0.000      0.000 3/25920
  580.0      1.772      1.772      1.772     -0.000      0.000      0.000 3/25920
  590.0      1.742      1.742      1.742     -0.000      0.000      0.000 3/25920
  600.0      1.713      1.713      1.713     -0.000      0.000      0.000 3/25920
  610.0      1.685      1.685      1.685     -0.000      0.000      0.000 3/25920
  620.0      1.658      1.658      1.658     -0.000      0.000      0.000 3/25920
  630.0      1.632      1.632      1.632     -0.000      0.000      0.000 3/25920
  640.0      1.606      1.606      1.606     -0.000      0.000      0.000 3/25920
  650.0      1.582      1.582      1.582     -0.000      0.000      0.000 3/25920
  660.0      1.558      1.558      1.558     -0.000      0.000      0.000 3/25920
  670.0      1.535      1.535      1.535     -0.000      0.000      0.000 3/25920
  680.0      1.512      1.512      1.512     -0.000      0.000      0.000 3/25920
  690.0      1.490      1.490      1.490     -0.000      0.000      0.000 3/25920
  700.0      1.469      1.469      1.469     -0.000      0.000      0.000 3/25920
  710.0      1.448      1.448      1.448     -0.000      0.000      0.000 3/25920
  720.0      1.428      1.428      1.428     -0.000      0.000      0.000 3/25920
  730.0      1.409      1.409      1.409     -0.000      0.000      0.000 3/25920
  740.0      1.390      1.390      1.390     -0.000      0.000      0.000 3/25920
  750.0      1.371      1.371      1.371     -0.000      0.000      0.000 3/25920
  760.0      1.353      1.353      1.353     -0.000      0.000      0.000 3/25920
  770.0      1.336      1.336      1.336     -0.000      0.000      0.000 3/25920
  780.0      1.318      1.318      1.318     -0.000      0.000      0.000 3/25920
  790.0      1.302      1.302      1.302     -0.000      0.000      0.000 3/25920
  800.0      1.286      1.286      1.286     -0.000      0.000      0.000 3/25920
  810.0      1.270      1.270      1.270     -0.000      0.000      0.000 3/25920
  820.0      1.254      1.254      1.254     -0.000      0.000      0.000 3/25920
  830.0      1.239      1.239      1.239     -0.000      0.000      0.000 3/25920
  840.0      1.224      1.224      1.224     -0.000      0.000      0.000 3/25920
  850.0      1.210      1.210      1.210     -0.000      0.000      0.000 3/25920
  860.0      1.196      1.196      1.196     -0.000      0.000      0.000 3/25920
  870.0      1.182      1.182      1.182     -0.000      0.000      0.000 3/25920
  880.0      1.169      1.169      1.169     -0.000      0.000      0.000 3/25920
  890.0      1.156      1.156      1.156     -0.000      0.000      0.000 3/25920
  900.0      1.143      1.143      1.143     -0.000      0.000      0.000 3/25920
  910.0      1.130      1.130      1.130     -0.000      0.000      0.000 3/25920
  920.0      1.118      1.118      1.118     -0.000      0.000      0.000 3/25920
  930.0      1.106      1.106      1.106     -0.000      0.000      0.000 3/25920
  940.0      1.094      1.094      1.094     -0.000      0.000      0.000 3/25920
  950.0      1.083      1.083      1.083     -0.000      0.000      0.000 3/25920
  960.0      1.072      1.072      1.072     -0.000      0.000      0.000 3/25920
  970.0      1.061      1.061      1.061     -0.000      0.000      0.000 3/25920
  980.0      1.050      1.050      1.050     -0.000      0.000      0.000 3/25920
  990.0      1.039      1.039      1.039     -0.000      0.000      0.000 3/25920
 1000.0      1.029      1.029      1.029     -0.000      0.000      0.000 3/25920

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

