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

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

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

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

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


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

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

-------------------------[time 2026-01-08 05:42:24]-------------------------
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.449197749999998    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.449197749999998    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.449197749999998
Atomic positions (fractional):
    1 H   0.50000000000000  0.00000000000000  0.50000000000000   1.008
    2 H   0.50000000000000  0.50000000000000  0.00000000000000   1.008
    3 H   0.00000000000000  0.50000000000000  0.50000000000000   1.008
    4 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078
    5 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.347593249999994    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.347593249999994    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.347593249999994
Atomic positions (fractional):
    1 H   0.16666666666667  0.00000000000000  0.16666666666667   1.008 > 1
    2 H   0.50000000000000  0.00000000000000  0.16666666666667   1.008 > 1
    3 H   0.83333333333333  0.00000000000000  0.16666666666667   1.008 > 1
    4 H   0.16666666666667  0.33333333333333  0.16666666666667   1.008 > 1
    5 H   0.50000000000000  0.33333333333333  0.16666666666667   1.008 > 1
    6 H   0.83333333333333  0.33333333333333  0.16666666666667   1.008 > 1
    7 H   0.16666666666667  0.66666666666667  0.16666666666667   1.008 > 1
    8 H   0.50000000000000  0.66666666666667  0.16666666666667   1.008 > 1
    9 H   0.83333333333333  0.66666666666667  0.16666666666667   1.008 > 1
   10 H   0.16666666666667  0.00000000000000  0.50000000000000   1.008 > 1
   11 H   0.50000000000000  0.00000000000000  0.50000000000000   1.008 > 1
   12 H   0.83333333333333  0.00000000000000  0.50000000000000   1.008 > 1
   13 H   0.16666666666667  0.33333333333333  0.50000000000000   1.008 > 1
   14 H   0.50000000000000  0.33333333333333  0.50000000000000   1.008 > 1
   15 H   0.83333333333333  0.33333333333333  0.50000000000000   1.008 > 1
   16 H   0.16666666666667  0.66666666666667  0.50000000000000   1.008 > 1
   17 H   0.50000000000000  0.66666666666667  0.50000000000000   1.008 > 1
   18 H   0.83333333333333  0.66666666666667  0.50000000000000   1.008 > 1
   19 H   0.16666666666667  0.00000000000000  0.83333333333333   1.008 > 1
   20 H   0.50000000000000  0.00000000000000  0.83333333333333   1.008 > 1
   21 H   0.83333333333333  0.00000000000000  0.83333333333333   1.008 > 1
   22 H   0.16666666666667  0.33333333333333  0.83333333333333   1.008 > 1
   23 H   0.50000000000000  0.33333333333333  0.83333333333333   1.008 > 1
   24 H   0.83333333333333  0.33333333333333  0.83333333333333   1.008 > 1
   25 H   0.16666666666667  0.66666666666667  0.83333333333333   1.008 > 1
   26 H   0.50000000000000  0.66666666666667  0.83333333333333   1.008 > 1
   27 H   0.83333333333333  0.66666666666667  0.83333333333333   1.008 > 1
   28 H   0.16666666666667  0.16666666666667  0.00000000000000   1.008 > 28
   29 H   0.50000000000000  0.16666666666667  0.00000000000000   1.008 > 28
   30 H   0.83333333333333  0.16666666666667  0.00000000000000   1.008 > 28
   31 H   0.16666666666667  0.50000000000000  0.00000000000000   1.008 > 28
   32 H   0.50000000000000  0.50000000000000  0.00000000000000   1.008 > 28
   33 H   0.83333333333333  0.50000000000000  0.00000000000000   1.008 > 28
   34 H   0.16666666666667  0.83333333333333  0.00000000000000   1.008 > 28
   35 H   0.50000000000000  0.83333333333333  0.00000000000000   1.008 > 28
   36 H   0.83333333333333  0.83333333333333  0.00000000000000   1.008 > 28
   37 H   0.16666666666667  0.16666666666667  0.33333333333333   1.008 > 28
   38 H   0.50000000000000  0.16666666666667  0.33333333333333   1.008 > 28
   39 H   0.83333333333333  0.16666666666667  0.33333333333333   1.008 > 28
   40 H   0.16666666666667  0.50000000000000  0.33333333333333   1.008 > 28
   41 H   0.50000000000000  0.50000000000000  0.33333333333333   1.008 > 28
   42 H   0.83333333333333  0.50000000000000  0.33333333333333   1.008 > 28
   43 H   0.16666666666667  0.83333333333333  0.33333333333333   1.008 > 28
   44 H   0.50000000000000  0.83333333333333  0.33333333333333   1.008 > 28
   45 H   0.83333333333333  0.83333333333333  0.33333333333333   1.008 > 28
   46 H   0.16666666666667  0.16666666666667  0.66666666666667   1.008 > 28
   47 H   0.50000000000000  0.16666666666667  0.66666666666667   1.008 > 28
   48 H   0.83333333333333  0.16666666666667  0.66666666666667   1.008 > 28
   49 H   0.16666666666667  0.50000000000000  0.66666666666667   1.008 > 28
   50 H   0.50000000000000  0.50000000000000  0.66666666666667   1.008 > 28
   51 H   0.83333333333333  0.50000000000000  0.66666666666667   1.008 > 28
   52 H   0.16666666666667  0.83333333333333  0.66666666666667   1.008 > 28
   53 H   0.50000000000000  0.83333333333333  0.66666666666667   1.008 > 28
   54 H   0.83333333333333  0.83333333333333  0.66666666666667   1.008 > 28
   55 H   0.00000000000000  0.16666666666667  0.16666666666667   1.008 > 55
   56 H   0.33333333333333  0.16666666666667  0.16666666666667   1.008 > 55
   57 H   0.66666666666667  0.16666666666667  0.16666666666667   1.008 > 55
   58 H   0.00000000000000  0.50000000000000  0.16666666666667   1.008 > 55
   59 H   0.33333333333333  0.50000000000000  0.16666666666667   1.008 > 55
   60 H   0.66666666666667  0.50000000000000  0.16666666666667   1.008 > 55
   61 H   0.00000000000000  0.83333333333333  0.16666666666667   1.008 > 55
   62 H   0.33333333333333  0.83333333333333  0.16666666666667   1.008 > 55
   63 H   0.66666666666667  0.83333333333333  0.16666666666667   1.008 > 55
   64 H   0.00000000000000  0.16666666666667  0.50000000000000   1.008 > 55
   65 H   0.33333333333333  0.16666666666667  0.50000000000000   1.008 > 55
   66 H   0.66666666666667  0.16666666666667  0.50000000000000   1.008 > 55
   67 H   0.00000000000000  0.50000000000000  0.50000000000000   1.008 > 55
   68 H   0.33333333333333  0.50000000000000  0.50000000000000   1.008 > 55
   69 H   0.66666666666667  0.50000000000000  0.50000000000000   1.008 > 55
   70 H   0.00000000000000  0.83333333333333  0.50000000000000   1.008 > 55
   71 H   0.33333333333333  0.83333333333333  0.50000000000000   1.008 > 55
   72 H   0.66666666666667  0.83333333333333  0.50000000000000   1.008 > 55
   73 H   0.00000000000000  0.16666666666667  0.83333333333333   1.008 > 55
   74 H   0.33333333333333  0.16666666666667  0.83333333333333   1.008 > 55
   75 H   0.66666666666667  0.16666666666667  0.83333333333333   1.008 > 55
   76 H   0.00000000000000  0.50000000000000  0.83333333333333   1.008 > 55
   77 H   0.33333333333333  0.50000000000000  0.83333333333333   1.008 > 55
   78 H   0.66666666666667  0.50000000000000  0.83333333333333   1.008 > 55
   79 H   0.00000000000000  0.83333333333333  0.83333333333333   1.008 > 55
   80 H   0.33333333333333  0.83333333333333  0.83333333333333   1.008 > 55
   81 H   0.66666666666667  0.83333333333333  0.83333333333333   1.008 > 55
   82 Ca  0.16666666666667  0.16666666666667  0.16666666666667  40.078 > 82
   83 Ca  0.50000000000000  0.16666666666667  0.16666666666667  40.078 > 82
   84 Ca  0.83333333333333  0.16666666666667  0.16666666666667  40.078 > 82
   85 Ca  0.16666666666667  0.50000000000000  0.16666666666667  40.078 > 82
   86 Ca  0.50000000000000  0.50000000000000  0.16666666666667  40.078 > 82
   87 Ca  0.83333333333333  0.50000000000000  0.16666666666667  40.078 > 82
   88 Ca  0.16666666666667  0.83333333333333  0.16666666666667  40.078 > 82
   89 Ca  0.50000000000000  0.83333333333333  0.16666666666667  40.078 > 82
   90 Ca  0.83333333333333  0.83333333333333  0.16666666666667  40.078 > 82
   91 Ca  0.16666666666667  0.16666666666667  0.50000000000000  40.078 > 82
   92 Ca  0.50000000000000  0.16666666666667  0.50000000000000  40.078 > 82
   93 Ca  0.83333333333333  0.16666666666667  0.50000000000000  40.078 > 82
   94 Ca  0.16666666666667  0.50000000000000  0.50000000000000  40.078 > 82
   95 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 82
   96 Ca  0.83333333333333  0.50000000000000  0.50000000000000  40.078 > 82
   97 Ca  0.16666666666667  0.83333333333333  0.50000000000000  40.078 > 82
   98 Ca  0.50000000000000  0.83333333333333  0.50000000000000  40.078 > 82
   99 Ca  0.83333333333333  0.83333333333333  0.50000000000000  40.078 > 82
  100 Ca  0.16666666666667  0.16666666666667  0.83333333333333  40.078 > 82
  101 Ca  0.50000000000000  0.16666666666667  0.83333333333333  40.078 > 82
  102 Ca  0.83333333333333  0.16666666666667  0.83333333333333  40.078 > 82
  103 Ca  0.16666666666667  0.50000000000000  0.83333333333333  40.078 > 82
  104 Ca  0.50000000000000  0.50000000000000  0.83333333333333  40.078 > 82
  105 Ca  0.83333333333333  0.50000000000000  0.83333333333333  40.078 > 82
  106 Ca  0.16666666666667  0.83333333333333  0.83333333333333  40.078 > 82
  107 Ca  0.50000000000000  0.83333333333333  0.83333333333333  40.078 > 82
  108 Ca  0.83333333333333  0.83333333333333  0.83333333333333  40.078 > 82
  109 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468 > 109
  110 Rb  0.33333333333333  0.00000000000000  0.00000000000000  85.468 > 109
  111 Rb  0.66666666666667  0.00000000000000  0.00000000000000  85.468 > 109
  112 Rb  0.00000000000000  0.33333333333333  0.00000000000000  85.468 > 109
  113 Rb  0.33333333333333  0.33333333333333  0.00000000000000  85.468 > 109
  114 Rb  0.66666666666667  0.33333333333333  0.00000000000000  85.468 > 109
  115 Rb  0.00000000000000  0.66666666666667  0.00000000000000  85.468 > 109
  116 Rb  0.33333333333333  0.66666666666667  0.00000000000000  85.468 > 109
  117 Rb  0.66666666666667  0.66666666666667  0.00000000000000  85.468 > 109
  118 Rb  0.00000000000000  0.00000000000000  0.33333333333333  85.468 > 109
  119 Rb  0.33333333333333  0.00000000000000  0.33333333333333  85.468 > 109
  120 Rb  0.66666666666667  0.00000000000000  0.33333333333333  85.468 > 109
  121 Rb  0.00000000000000  0.33333333333333  0.33333333333333  85.468 > 109
  122 Rb  0.33333333333333  0.33333333333333  0.33333333333333  85.468 > 109
  123 Rb  0.66666666666667  0.33333333333333  0.33333333333333  85.468 > 109
  124 Rb  0.00000000000000  0.66666666666667  0.33333333333333  85.468 > 109
  125 Rb  0.33333333333333  0.66666666666667  0.33333333333333  85.468 > 109
  126 Rb  0.66666666666667  0.66666666666667  0.33333333333333  85.468 > 109
  127 Rb  0.00000000000000  0.00000000000000  0.66666666666667  85.468 > 109
  128 Rb  0.33333333333333  0.00000000000000  0.66666666666667  85.468 > 109
  129 Rb  0.66666666666667  0.00000000000000  0.66666666666667  85.468 > 109
  130 Rb  0.00000000000000  0.33333333333333  0.66666666666667  85.468 > 109
  131 Rb  0.33333333333333  0.33333333333333  0.66666666666667  85.468 > 109
  132 Rb  0.66666666666667  0.33333333333333  0.66666666666667  85.468 > 109
  133 Rb  0.00000000000000  0.66666666666667  0.66666666666667  85.468 > 109
  134 Rb  0.33333333333333  0.66666666666667  0.66666666666667  85.468 > 109
  135 Rb  0.66666666666667  0.66666666666667  0.66666666666667  85.468 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.3619700    0.0000000    0.0000000
            0.0000000    3.3619700    0.0000000
            0.0000000    0.0000000    3.3619700
-------------------------- Born effective charges --------------------------
    1 H    -0.8777473    0.0000000    0.0000000
            0.0000000   -1.0960474    0.0000000
            0.0000000    0.0000000   -0.8777473
    2 H    -0.8777473    0.0000000    0.0000000
            0.0000000   -0.8777473    0.0000000
            0.0000000    0.0000000   -1.0960474
    3 H    -1.0960474    0.0000000    0.0000000
            0.0000000   -0.8777473    0.0000000
            0.0000000    0.0000000   -0.8777473
    4 Ca    1.7535224    0.0000000    0.0000000
            0.0000000    1.7535224    0.0000000
            0.0000000    0.0000000    1.7535224
    5 Rb    1.0980195    0.0000000    0.0000000
            0.0000000    1.0980195    0.0000000
            0.0000000    0.0000000    1.0980195
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
    [ 0.0000  0.0100  0.0000]
    [ 0.0000 -0.0100  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-08 05:42:27]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 05:42:28]-------------------------
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.449197749999998    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.449197749999998    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.449197749999998
Atomic positions (fractional):
    1 H   0.50000000000000  0.00000000000000  0.50000000000000   1.008
    2 H   0.50000000000000  0.50000000000000  0.00000000000000   1.008
    3 H   0.00000000000000  0.50000000000000  0.50000000000000   1.008
    4 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078
    5 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   13.347593249999994    0.000000000000000    0.000000000000000
  b    0.000000000000000   13.347593249999994    0.000000000000000
  c    0.000000000000000    0.000000000000000   13.347593249999994
Atomic positions (fractional):
    1 H   0.16666666666667  0.00000000000000  0.16666666666667   1.008 > 1
    2 H   0.50000000000000  0.00000000000000  0.16666666666667   1.008 > 1
    3 H   0.83333333333333  0.00000000000000  0.16666666666667   1.008 > 1
    4 H   0.16666666666667  0.33333333333333  0.16666666666667   1.008 > 1
    5 H   0.50000000000000  0.33333333333333  0.16666666666667   1.008 > 1
    6 H   0.83333333333333  0.33333333333333  0.16666666666667   1.008 > 1
    7 H   0.16666666666667  0.66666666666667  0.16666666666667   1.008 > 1
    8 H   0.50000000000000  0.66666666666667  0.16666666666667   1.008 > 1
    9 H   0.83333333333333  0.66666666666667  0.16666666666667   1.008 > 1
   10 H   0.16666666666667  0.00000000000000  0.50000000000000   1.008 > 1
   11 H   0.50000000000000  0.00000000000000  0.50000000000000   1.008 > 1
   12 H   0.83333333333333  0.00000000000000  0.50000000000000   1.008 > 1
   13 H   0.16666666666667  0.33333333333333  0.50000000000000   1.008 > 1
   14 H   0.50000000000000  0.33333333333333  0.50000000000000   1.008 > 1
   15 H   0.83333333333333  0.33333333333333  0.50000000000000   1.008 > 1
   16 H   0.16666666666667  0.66666666666667  0.50000000000000   1.008 > 1
   17 H   0.50000000000000  0.66666666666667  0.50000000000000   1.008 > 1
   18 H   0.83333333333333  0.66666666666667  0.50000000000000   1.008 > 1
   19 H   0.16666666666667  0.00000000000000  0.83333333333333   1.008 > 1
   20 H   0.50000000000000  0.00000000000000  0.83333333333333   1.008 > 1
   21 H   0.83333333333333  0.00000000000000  0.83333333333333   1.008 > 1
   22 H   0.16666666666667  0.33333333333333  0.83333333333333   1.008 > 1
   23 H   0.50000000000000  0.33333333333333  0.83333333333333   1.008 > 1
   24 H   0.83333333333333  0.33333333333333  0.83333333333333   1.008 > 1
   25 H   0.16666666666667  0.66666666666667  0.83333333333333   1.008 > 1
   26 H   0.50000000000000  0.66666666666667  0.83333333333333   1.008 > 1
   27 H   0.83333333333333  0.66666666666667  0.83333333333333   1.008 > 1
   28 H   0.16666666666667  0.16666666666667  0.00000000000000   1.008 > 28
   29 H   0.50000000000000  0.16666666666667  0.00000000000000   1.008 > 28
   30 H   0.83333333333333  0.16666666666667  0.00000000000000   1.008 > 28
   31 H   0.16666666666667  0.50000000000000  0.00000000000000   1.008 > 28
   32 H   0.50000000000000  0.50000000000000  0.00000000000000   1.008 > 28
   33 H   0.83333333333333  0.50000000000000  0.00000000000000   1.008 > 28
   34 H   0.16666666666667  0.83333333333333  0.00000000000000   1.008 > 28
   35 H   0.50000000000000  0.83333333333333  0.00000000000000   1.008 > 28
   36 H   0.83333333333333  0.83333333333333  0.00000000000000   1.008 > 28
   37 H   0.16666666666667  0.16666666666667  0.33333333333333   1.008 > 28
   38 H   0.50000000000000  0.16666666666667  0.33333333333333   1.008 > 28
   39 H   0.83333333333333  0.16666666666667  0.33333333333333   1.008 > 28
   40 H   0.16666666666667  0.50000000000000  0.33333333333333   1.008 > 28
   41 H   0.50000000000000  0.50000000000000  0.33333333333333   1.008 > 28
   42 H   0.83333333333333  0.50000000000000  0.33333333333333   1.008 > 28
   43 H   0.16666666666667  0.83333333333333  0.33333333333333   1.008 > 28
   44 H   0.50000000000000  0.83333333333333  0.33333333333333   1.008 > 28
   45 H   0.83333333333333  0.83333333333333  0.33333333333333   1.008 > 28
   46 H   0.16666666666667  0.16666666666667  0.66666666666667   1.008 > 28
   47 H   0.50000000000000  0.16666666666667  0.66666666666667   1.008 > 28
   48 H   0.83333333333333  0.16666666666667  0.66666666666667   1.008 > 28
   49 H   0.16666666666667  0.50000000000000  0.66666666666667   1.008 > 28
   50 H   0.50000000000000  0.50000000000000  0.66666666666667   1.008 > 28
   51 H   0.83333333333333  0.50000000000000  0.66666666666667   1.008 > 28
   52 H   0.16666666666667  0.83333333333333  0.66666666666667   1.008 > 28
   53 H   0.50000000000000  0.83333333333333  0.66666666666667   1.008 > 28
   54 H   0.83333333333333  0.83333333333333  0.66666666666667   1.008 > 28
   55 H   0.00000000000000  0.16666666666667  0.16666666666667   1.008 > 55
   56 H   0.33333333333333  0.16666666666667  0.16666666666667   1.008 > 55
   57 H   0.66666666666667  0.16666666666667  0.16666666666667   1.008 > 55
   58 H   0.00000000000000  0.50000000000000  0.16666666666667   1.008 > 55
   59 H   0.33333333333333  0.50000000000000  0.16666666666667   1.008 > 55
   60 H   0.66666666666667  0.50000000000000  0.16666666666667   1.008 > 55
   61 H   0.00000000000000  0.83333333333333  0.16666666666667   1.008 > 55
   62 H   0.33333333333333  0.83333333333333  0.16666666666667   1.008 > 55
   63 H   0.66666666666667  0.83333333333333  0.16666666666667   1.008 > 55
   64 H   0.00000000000000  0.16666666666667  0.50000000000000   1.008 > 55
   65 H   0.33333333333333  0.16666666666667  0.50000000000000   1.008 > 55
   66 H   0.66666666666667  0.16666666666667  0.50000000000000   1.008 > 55
   67 H   0.00000000000000  0.50000000000000  0.50000000000000   1.008 > 55
   68 H   0.33333333333333  0.50000000000000  0.50000000000000   1.008 > 55
   69 H   0.66666666666667  0.50000000000000  0.50000000000000   1.008 > 55
   70 H   0.00000000000000  0.83333333333333  0.50000000000000   1.008 > 55
   71 H   0.33333333333333  0.83333333333333  0.50000000000000   1.008 > 55
   72 H   0.66666666666667  0.83333333333333  0.50000000000000   1.008 > 55
   73 H   0.00000000000000  0.16666666666667  0.83333333333333   1.008 > 55
   74 H   0.33333333333333  0.16666666666667  0.83333333333333   1.008 > 55
   75 H   0.66666666666667  0.16666666666667  0.83333333333333   1.008 > 55
   76 H   0.00000000000000  0.50000000000000  0.83333333333333   1.008 > 55
   77 H   0.33333333333333  0.50000000000000  0.83333333333333   1.008 > 55
   78 H   0.66666666666667  0.50000000000000  0.83333333333333   1.008 > 55
   79 H   0.00000000000000  0.83333333333333  0.83333333333333   1.008 > 55
   80 H   0.33333333333333  0.83333333333333  0.83333333333333   1.008 > 55
   81 H   0.66666666666667  0.83333333333333  0.83333333333333   1.008 > 55
   82 Ca  0.16666666666667  0.16666666666667  0.16666666666667  40.078 > 82
   83 Ca  0.50000000000000  0.16666666666667  0.16666666666667  40.078 > 82
   84 Ca  0.83333333333333  0.16666666666667  0.16666666666667  40.078 > 82
   85 Ca  0.16666666666667  0.50000000000000  0.16666666666667  40.078 > 82
   86 Ca  0.50000000000000  0.50000000000000  0.16666666666667  40.078 > 82
   87 Ca  0.83333333333333  0.50000000000000  0.16666666666667  40.078 > 82
   88 Ca  0.16666666666667  0.83333333333333  0.16666666666667  40.078 > 82
   89 Ca  0.50000000000000  0.83333333333333  0.16666666666667  40.078 > 82
   90 Ca  0.83333333333333  0.83333333333333  0.16666666666667  40.078 > 82
   91 Ca  0.16666666666667  0.16666666666667  0.50000000000000  40.078 > 82
   92 Ca  0.50000000000000  0.16666666666667  0.50000000000000  40.078 > 82
   93 Ca  0.83333333333333  0.16666666666667  0.50000000000000  40.078 > 82
   94 Ca  0.16666666666667  0.50000000000000  0.50000000000000  40.078 > 82
   95 Ca  0.50000000000000  0.50000000000000  0.50000000000000  40.078 > 82
   96 Ca  0.83333333333333  0.50000000000000  0.50000000000000  40.078 > 82
   97 Ca  0.16666666666667  0.83333333333333  0.50000000000000  40.078 > 82
   98 Ca  0.50000000000000  0.83333333333333  0.50000000000000  40.078 > 82
   99 Ca  0.83333333333333  0.83333333333333  0.50000000000000  40.078 > 82
  100 Ca  0.16666666666667  0.16666666666667  0.83333333333333  40.078 > 82
  101 Ca  0.50000000000000  0.16666666666667  0.83333333333333  40.078 > 82
  102 Ca  0.83333333333333  0.16666666666667  0.83333333333333  40.078 > 82
  103 Ca  0.16666666666667  0.50000000000000  0.83333333333333  40.078 > 82
  104 Ca  0.50000000000000  0.50000000000000  0.83333333333333  40.078 > 82
  105 Ca  0.83333333333333  0.50000000000000  0.83333333333333  40.078 > 82
  106 Ca  0.16666666666667  0.83333333333333  0.83333333333333  40.078 > 82
  107 Ca  0.50000000000000  0.83333333333333  0.83333333333333  40.078 > 82
  108 Ca  0.83333333333333  0.83333333333333  0.83333333333333  40.078 > 82
  109 Rb  0.00000000000000  0.00000000000000  0.00000000000000  85.468 > 109
  110 Rb  0.33333333333333  0.00000000000000  0.00000000000000  85.468 > 109
  111 Rb  0.66666666666667  0.00000000000000  0.00000000000000  85.468 > 109
  112 Rb  0.00000000000000  0.33333333333333  0.00000000000000  85.468 > 109
  113 Rb  0.33333333333333  0.33333333333333  0.00000000000000  85.468 > 109
  114 Rb  0.66666666666667  0.33333333333333  0.00000000000000  85.468 > 109
  115 Rb  0.00000000000000  0.66666666666667  0.00000000000000  85.468 > 109
  116 Rb  0.33333333333333  0.66666666666667  0.00000000000000  85.468 > 109
  117 Rb  0.66666666666667  0.66666666666667  0.00000000000000  85.468 > 109
  118 Rb  0.00000000000000  0.00000000000000  0.33333333333333  85.468 > 109
  119 Rb  0.33333333333333  0.00000000000000  0.33333333333333  85.468 > 109
  120 Rb  0.66666666666667  0.00000000000000  0.33333333333333  85.468 > 109
  121 Rb  0.00000000000000  0.33333333333333  0.33333333333333  85.468 > 109
  122 Rb  0.33333333333333  0.33333333333333  0.33333333333333  85.468 > 109
  123 Rb  0.66666666666667  0.33333333333333  0.33333333333333  85.468 > 109
  124 Rb  0.00000000000000  0.66666666666667  0.33333333333333  85.468 > 109
  125 Rb  0.33333333333333  0.66666666666667  0.33333333333333  85.468 > 109
  126 Rb  0.66666666666667  0.66666666666667  0.33333333333333  85.468 > 109
  127 Rb  0.00000000000000  0.00000000000000  0.66666666666667  85.468 > 109
  128 Rb  0.33333333333333  0.00000000000000  0.66666666666667  85.468 > 109
  129 Rb  0.66666666666667  0.00000000000000  0.66666666666667  85.468 > 109
  130 Rb  0.00000000000000  0.33333333333333  0.66666666666667  85.468 > 109
  131 Rb  0.33333333333333  0.33333333333333  0.66666666666667  85.468 > 109
  132 Rb  0.66666666666667  0.33333333333333  0.66666666666667  85.468 > 109
  133 Rb  0.00000000000000  0.66666666666667  0.66666666666667  85.468 > 109
  134 Rb  0.33333333333333  0.66666666666667  0.66666666666667  85.468 > 109
  135 Rb  0.66666666666667  0.66666666666667  0.66666666666667  85.468 > 109
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.3619700    0.0000000    0.0000000
            0.0000000    3.3619700    0.0000000
            0.0000000    0.0000000    3.3619700
-------------------------- Born effective charges --------------------------
    1 H    -0.8777473    0.0000000    0.0000000
            0.0000000   -1.0960474    0.0000000
            0.0000000    0.0000000   -0.8777473
    2 H    -0.8777473    0.0000000    0.0000000
            0.0000000   -0.8777473    0.0000000
            0.0000000    0.0000000   -1.0960474
    3 H    -1.0960474    0.0000000    0.0000000
            0.0000000   -0.8777473    0.0000000
            0.0000000    0.0000000   -0.8777473
    4 Ca    1.7535224    0.0000000    0.0000000
            0.0000000    1.7535224    0.0000000
            0.0000000    0.0000000    1.7535224
    5 Rb    1.0980195    0.0000000    0.0000000
            0.0000000    1.0980195    0.0000000
            0.0000000    0.0000000    1.0980195
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (yyy) 0.00000000 (yyy) 0.00000000 (yyy)
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
Length for sampling mesh generation: 50.00
Generating grid system ... [ 11 11 11 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 0.93, Number of G-points: 305, Lambda: 0.18
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/56) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 56
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
  -0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   0.000   (   0.000    0.000    0.000)    0.000
   3.262   (   0.000    0.000    0.000)    0.000
   3.262   (   0.000    0.000    0.000)    0.000
   3.262   (   0.000    0.000    0.000)    0.000
  11.227   (   0.000    0.000    0.000)    0.000
  11.227   (   0.000    0.000    0.000)    0.000
  11.227   (   0.000    0.000    0.000)    0.000
  12.823   (   0.000    0.000    0.000)    0.000
  12.823   (   0.000    0.000    0.000)    0.000
  12.823   (   0.000    0.000    0.000)    0.000
  34.955   (   0.000    0.000    0.000)    0.000
  34.955   (   0.000    0.000    0.000)    0.000
  34.955   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/56) =======================
q-point: ( 0.09  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.557   (  26.104    0.000    0.000)   26.104
   0.557   (  26.104    0.000    0.000)   26.104
   0.969   (  43.122    0.000    0.000)   43.122
   3.266   (   0.373    0.000    0.000)    0.373
   3.266   (   0.373    0.000    0.000)    0.373
   3.743   (  16.814    0.000    0.000)   16.814
  11.347   (  11.409    0.000    0.000)   11.409
  11.347   (  11.409    0.000    0.000)   11.409
  11.378   (  14.312    0.000    0.000)   14.312
  12.855   (   3.018    0.000    0.000)    3.018
  12.855   (   3.018    0.000    0.000)    3.018
  19.692   (   5.187    0.000    0.000)    5.187
  34.952   (  -0.240    0.000    0.000)    0.240
  34.952   (  -0.240    0.000    0.000)    0.240
  37.715   (   5.036    0.000    0.000)    5.036
======================= Grid point 2 (3/56) =======================
q-point: ( 0.18  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.049   (  21.702    0.000    0.000)   21.702
   1.049   (  21.702    0.000    0.000)   21.702
   1.711   (  29.238    0.000    0.000)   29.238
   3.273   (   0.200    0.000    0.000)    0.200
   3.273   (   0.200    0.000    0.000)    0.200
   4.223   (  28.553    0.000    0.000)   28.553
  11.667   (  19.080    0.000    0.000)   19.080
  11.667   (  19.080    0.000    0.000)   19.080
  11.774   (  23.233    0.000    0.000)   23.233
  12.938   (   4.896    0.000    0.000)    4.896
  12.938   (   4.896    0.000    0.000)    4.896
  19.837   (   8.537    0.000    0.000)    8.537
  34.945   (  -0.386    0.000    0.000)    0.386
  34.945   (  -0.386    0.000    0.000)    0.386
  37.858   (   8.646    0.000    0.000)    8.646
======================= Grid point 3 (4/56) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.436   (  16.149    0.000    0.000)   16.149
   1.436   (  16.149    0.000    0.000)   16.149
   2.184   (  17.811    0.000    0.000)   17.811
   3.271   (  -0.467    0.000    0.000)    0.467
   3.271   (  -0.467    0.000    0.000)    0.467
   4.829   (  28.831    0.000    0.000)   28.831
  12.086   (  20.887    0.000    0.000)   20.887
  12.086   (  20.887    0.000    0.000)   20.887
  12.273   (  24.204    0.000    0.000)   24.204
  13.041   (   4.772    0.000    0.000)    4.772
  13.041   (   4.772    0.000    0.000)    4.772
  20.021   (   9.007    0.000    0.000)    9.007
  34.937   (  -0.395    0.000    0.000)    0.395
  34.937   (  -0.395    0.000    0.000)    0.395
  38.050   (   9.648    0.000    0.000)    9.648
======================= Grid point 4 (5/56) =======================
q-point: ( 0.36  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.706   (  10.134    0.000    0.000)   10.134
   1.706   (  10.134    0.000    0.000)   10.134
   2.463   (   9.910    0.000    0.000)    9.910
   3.255   (  -0.928    0.000    0.000)    0.928
   3.255   (  -0.928    0.000    0.000)    0.928
   5.339   (  19.982    0.000    0.000)   19.982
  12.480   (  16.757    0.000    0.000)   16.757
  12.480   (  16.757    0.000    0.000)   16.757
  12.713   (  17.807    0.000    0.000)   17.807
  13.119   (   2.529    0.000    0.000)    2.529
  13.119   (   2.529    0.000    0.000)    2.529
  20.186   (   6.671    0.000    0.000)    6.671
  34.930   (  -0.284    0.000    0.000)    0.284
  34.930   (  -0.284    0.000    0.000)    0.284
  38.231   (   7.530    0.000    0.000)    7.530
======================= Grid point 5 (6/56) =======================
q-point: ( 0.45  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 126
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.846   (   3.511    0.000    0.000)    3.511
   1.846   (   3.511    0.000    0.000)    3.511
   2.597   (   3.236    0.000    0.000)    3.236
   3.239   (  -0.498    0.000    0.000)    0.498
   3.239   (  -0.498    0.000    0.000)    0.498
   5.618   (   7.019    0.000    0.000)    7.019
  12.734   (   7.044    0.000    0.000)    7.044
  12.734   (   7.044    0.000    0.000)    7.044
  12.967   (   6.501    0.000    0.000)    6.501
  13.142   (   0.079    0.000    0.000)    0.079
  13.142   (   0.079    0.000    0.000)    0.079
  20.281   (   2.440    0.000    0.000)    2.440
  34.926   (  -0.102    0.000    0.000)    0.102
  34.926   (  -0.102    0.000    0.000)    0.102
  38.340   (   2.842    0.000    0.000)    2.842
======================= Grid point 12 (7/56) =======================
q-point: ( 0.09  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.777   (  17.639   17.639    0.000)   24.945
   0.855   (  18.221   18.221    0.000)   25.769
   1.278   (  26.842   26.842    0.000)   37.960
   3.274   (   0.771    0.771    0.000)    1.090
   3.492   (  10.997   10.997    0.000)   15.552
   3.696   (   6.868    6.868    0.000)    9.713
  11.130   (  -6.027   -6.027    0.000)    8.523
  11.463   (  11.043   11.043    0.000)   15.618
  11.807   (  26.204   26.204    0.000)   37.058
  12.869   (   3.021    3.021    0.000)    4.273
  12.888   (   3.162    3.162    0.000)    4.472
  19.766   (   5.554    5.554    0.000)    7.855
  34.949   (  -0.251   -0.251    0.000)    0.355
  34.973   (   0.956    0.956    0.000)    1.352
  37.752   (   4.589    4.589    0.000)    6.490
======================= Grid point 13 (8/56) =======================
q-point: ( 0.18  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.158   (  18.189    9.645    0.000)   20.588
   1.209   (  16.428   11.843    0.000)   20.251
   1.821   (  23.048    9.909    0.000)   25.088
   3.293   (   0.915    1.878    0.000)    2.089
   3.526   (  -0.162   23.328    0.000)   23.328
   4.180   (  31.072   -2.883    0.000)   31.205
  11.122   (   4.713  -28.841    0.000)   29.223
  11.774   (  18.509   10.110    0.000)   21.091
  12.351   (  24.376   33.962    0.000)   41.804
  12.976   (   5.133    3.543    0.000)    6.237
  13.011   (  11.766    5.547    0.000)   13.008
  19.897   (   7.079    4.934    0.000)    8.629
  34.942   (  -0.404   -0.279    0.000)    0.491
  34.982   (   0.145    2.935    0.000)    2.939
  37.899   (   9.310    4.563    0.000)   10.369
======================= Grid point 14 (9/56) =======================
q-point: ( 0.27  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.490   (  14.021    4.829    0.000)   14.830
   1.514   (  13.143    6.045    0.000)   14.466
   2.193   (  13.777    0.716    0.000)   13.796
   3.307   (   0.401    3.432    0.000)    3.456
   3.538   (   1.101   24.815    0.000)   24.839
   4.822   (  29.941   -0.299    0.000)   29.942
  11.288   (  10.522  -49.189    0.000)   50.302
  12.180   (  20.317    8.979    0.000)   22.213
  12.715   (  11.355   23.895    0.000)   26.456
  13.084   (   5.023    4.031    0.000)    6.441
  13.373   (  22.157   22.993    0.000)   31.931
  20.046   (   7.148    2.188    0.000)    7.475
  34.934   (  -0.413   -0.314    0.000)    0.519
  34.983   (  -0.046    4.039    0.000)    4.040
  38.109   (  10.546    5.904    0.000)   12.086
======================= Grid point 15 (10/56) =======================
q-point: ( 0.36  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.724   (   8.763    1.629    0.000)    8.913
   1.736   (   8.402    2.405    0.000)    8.739
   2.403   (   7.185   -4.959    0.000)    8.730
   3.309   (  -0.159    5.009    0.000)    5.012
   3.561   (   1.056   27.488    0.000)   27.508
   5.346   (  20.380    0.803    0.000)   20.396
  11.507   (   9.864  -64.578    0.000)   65.327
  12.564   (  16.355    8.013    0.000)   18.212
  12.872   (   5.053   17.062    0.000)   17.795
  13.166   (   2.716    4.475    0.000)    5.235
  13.807   (  18.209   34.858    0.000)   39.328
  20.175   (   5.250   -0.682    0.000)    5.295
  34.926   (  -0.297   -0.344    0.000)    0.454
  34.982   (  -0.048    4.721    0.000)    4.722
  38.306   (   8.210    7.376    0.000)   11.037
======================= Grid point 16 (11/56) =======================
q-point: ( 0.45  0.09  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.845   (   3.003   -0.122    0.000)    3.006
   1.852   (   2.898    0.548    0.000)    2.949
   2.497   (   2.222   -7.685    0.000)    8.000
   3.305   (  -0.174    6.051    0.000)    6.054
   3.577   (   0.403   29.328    0.000)   29.331
   5.630   (   7.115    1.154    0.000)    7.208
  11.655   (   3.968  -73.209    0.000)   73.316
  12.812   (   6.923    7.492    0.000)   10.201
  12.939   (   1.569   16.215    0.000)   16.291
  13.192   (   0.125    4.712    0.000)    4.714
  14.068   (   6.697   38.682    0.000)   39.258
  20.250   (   1.918   -2.398    0.000)    3.071
  34.922   (  -0.107   -0.361    0.000)    0.377
  34.981   (  -0.016    5.079    0.000)    5.079
  38.425   (   3.090    8.259    0.000)    8.818
======================= Grid point 24 (12/56) =======================
q-point: ( 0.18  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.388   (  11.725   11.725    0.000)   16.582
   1.439   (  10.412   10.412    0.000)   14.724
   2.042   (  10.038   10.038    0.000)   14.196
   3.343   (   2.815    2.815    0.000)    3.980
   4.100   (  17.566   17.566    0.000)   24.842
   4.190   (  17.067   17.067    0.000)   24.136
  10.603   ( -20.793  -20.793    0.000)   29.406
  12.058   (  17.032   17.032    0.000)   24.087
  13.074   (   5.764    5.764    0.000)    8.151
  13.086   (  32.480   32.480    0.000)   45.934
  13.130   (   9.744    9.744    0.000)   13.780
  20.004   (   5.143    5.143    0.000)    7.274
  34.935   (  -0.449   -0.449    0.000)    0.635
  35.040   (   2.375    2.375    0.000)    3.359
  38.054   (  10.276   10.276    0.000)   14.532
======================= Grid point 25 (13/56) =======================
q-point: ( 0.27  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.607   (   9.343    5.993    0.000)   11.100
   1.635   (   8.686    5.497    0.000)   10.279
   2.203   (   5.775    0.093    0.000)    5.775
   3.399   (   2.519    5.228    0.000)    5.803
   4.183   (   3.254   35.357    0.000)   35.506
   4.836   (  31.226    2.026    0.000)   31.292
  10.283   ( -10.816  -46.947    0.000)   48.176
  12.434   (  18.813   15.194    0.000)   24.182
  13.196   (   5.711    6.596    0.000)    8.725
  13.226   (   4.715   23.958    0.000)   24.418
  13.801   (  31.194   18.871    0.000)   36.458
  20.101   (   4.281    3.043    0.000)    5.252
  34.925   (  -0.459   -0.504    0.000)    0.681
  35.078   (   1.446    4.749    0.000)    4.964
  38.294   (  12.383   11.969    0.000)   17.222
======================= Grid point 26 (14/56) =======================
q-point: ( 0.36  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.762   (   5.744    1.863    0.000)    6.038
   1.786   (   5.828    2.313    0.000)    6.270
   2.288   (   2.837   -5.882    0.000)    6.530
   3.441   (   1.522    7.362    0.000)    7.518
   4.251   (   3.007   37.020    0.000)   37.142
   5.375   (  20.740    2.187    0.000)   20.855
  10.138   (  -4.038  -67.168    0.000)   67.289
  12.791   (  15.259   13.560    0.000)   20.414
  13.291   (   3.251    7.395    0.000)    8.078
  13.311   (   3.499   23.747    0.000)   24.004
  14.333   (  20.317   16.694    0.000)   26.295
  20.176   (   2.980    1.141    0.000)    3.191
  34.917   (  -0.329   -0.551    0.000)    0.642
  35.101   (   0.828    6.359    0.000)    6.413
  38.527   (   9.774   13.949    0.000)   17.033
======================= Grid point 27 (15/56) =======================
q-point: ( 0.45  0.18  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.840   (   1.928   -0.320    0.000)    1.955
   1.867   (   2.027    0.854    0.000)    2.200
   2.326   (   0.879   -8.696    0.000)    8.741
   3.461   (   0.477    8.600    0.000)    8.614
   4.295   (   1.156   38.156    0.000)   38.173
   5.663   (   7.151    1.977    0.000)    7.419
  10.093   (  -0.848  -77.928    0.000)   77.933
  13.024   (   6.596   12.644    0.000)   14.261
  13.324   (   0.255    7.852    0.000)    7.856
  13.362   (   1.317   23.301    0.000)   23.338
  14.614   (   7.000   14.916    0.000)   16.477
  20.219   (   1.089    0.094    0.000)    1.093
  34.912   (  -0.118   -0.577    0.000)    0.589
  35.112   (   0.281    7.159    0.000)    7.165
  38.669   (   3.681   15.214    0.000)   15.653
======================= Grid point 36 (16/56) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.718   (   4.681    4.681    0.000)    6.620
   1.741   (   4.878    4.878    0.000)    6.899
   2.195   (  -0.849   -0.849    0.000)    1.201
   3.506   (   4.849    4.849    0.000)    6.858
   4.830   (  17.149   17.149    0.000)   24.252
   4.958   (  18.583   18.583    0.000)   26.280
   9.428   ( -35.917  -35.917    0.000)   50.794
  12.770   (  16.851   16.851    0.000)   23.832
  13.336   (   6.676    6.676    0.000)    9.442
  13.606   (  12.545   12.545    0.000)   17.741
  14.183   (  18.748   18.748    0.000)   26.514
  20.166   (   3.200    3.200    0.000)    4.526
  34.914   (  -0.514   -0.514    0.000)    0.726
  35.165   (   3.596    3.596    0.000)    5.086
  38.573   (  14.521   14.521    0.000)   20.535
======================= Grid point 37 (17/56) =======================
q-point: ( 0.36  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.795   (   2.786    1.297    0.000)    3.073
   1.831   (   3.794    2.191    0.000)    4.381
   2.173   (  -1.076   -5.199    0.000)    5.309
   3.588   (   3.122    6.620    0.000)    7.319
   4.972   (   4.705   31.520    0.000)   31.869
   5.441   (  19.729    4.566    0.000)   20.250
   8.836   ( -21.730  -58.650    0.000)   62.546
  13.090   (  13.721   14.961    0.000)   20.300
  13.450   (   4.090    7.717    0.000)    8.734
  13.770   (   4.629   19.797    0.000)   20.331
  14.526   (  13.705    3.689    0.000)   14.193
  20.227   (   2.661    3.672    0.000)    4.535
  34.905   (  -0.368   -0.560    0.000)    0.670
  35.223   (   2.172    5.352    0.000)    5.776
  38.849   (  11.587   16.594    0.000)   20.239
======================= Grid point 38 (18/56) =======================
q-point: ( 0.45  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.832   (   0.916   -0.440    0.000)    1.016
   1.887   (   1.475    1.094    0.000)    1.836
   2.156   (  -0.478   -7.572    0.000)    7.587
   3.631   (   1.041    7.538    0.000)    7.610
   5.045   (   1.998   33.497    0.000)   33.556
   5.710   (   6.484    2.616    0.000)    6.992
   8.543   (  -7.104  -72.001    0.000)   72.351
  13.301   (   6.128   13.776    0.000)   15.078
  13.495   (   0.474    8.425    0.000)    8.439
  13.825   (   1.185   20.874    0.000)   20.908
  14.718   (   4.781   -3.052    0.000)    5.672
  20.267   (   1.067    4.371    0.000)    4.499
  34.900   (  -0.132   -0.585    0.000)    0.600
  35.253   (   0.732    6.275    0.000)    6.317
  39.017   (   4.379   17.886    0.000)   18.414
======================= Grid point 48 (19/56) =======================
q-point: ( 0.36  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.815   (   0.688    0.688    0.000)    0.973
   1.875   (   2.064    2.064    0.000)    2.919
   2.081   (  -3.727   -3.727    0.000)    5.271
   3.700   (   4.204    4.204    0.000)    5.945
   5.442   (  12.351   12.351    0.000)   17.466
   5.605   (  12.378   12.378    0.000)   17.505
   7.805   ( -40.764  -40.764    0.000)   57.648
  13.372   (  11.997   11.997    0.000)   16.966
  13.587   (   5.203    5.203    0.000)    7.359
  14.079   (   9.850    9.850    0.000)   13.930
  14.553   (   0.508    0.508    0.000)    0.718
  20.310   (   3.992    3.992    0.000)    5.645
  34.895   (  -0.400   -0.400    0.000)    0.565
  35.313   (   3.372    3.372    0.000)    4.769
  39.164   (  13.234   13.234    0.000)   18.716
======================= Grid point 49 (20/56) =======================
q-point: ( 0.45  0.36  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.824   (   0.210   -0.374    0.000)    0.429
   1.910   (   1.088    1.080    0.000)    1.533
   2.026   (  -1.558   -5.069    0.000)    5.303
   3.758   (   1.399    4.728    0.000)    4.931
   5.606   (   5.363   19.773    0.000)   20.487
   5.783   (   2.991    5.657    0.000)    6.399
   7.231   ( -14.351  -54.356    0.000)   56.218
  13.557   (   5.484   10.565    0.000)   11.903
  13.649   (   0.899    6.287    0.000)    6.350
  14.182   (   1.482   13.558    0.000)   13.639
  14.576   (   1.216   -8.728    0.000)    8.812
  20.372   (   1.683    5.153    0.000)    5.420
  34.890   (  -0.143   -0.417    0.000)    0.441
  35.360   (   1.146    4.032    0.000)    4.192
  39.356   (   5.007   14.235    0.000)   15.090
======================= Grid point 60 (21/56) =======================
q-point: ( 0.45  0.45  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 216
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.818   (  -0.145   -0.145    0.000)    0.205
   1.928   (   0.613    0.613    0.000)    0.867
   1.954   (  -1.899   -1.899    0.000)    2.686
   3.822   (   1.568    1.568    0.000)    2.217
   5.838   (   7.904    7.904    0.000)   11.177
   5.943   (   4.117    4.117    0.000)    5.822
   6.407   ( -23.673  -23.673    0.000)   33.478
  13.712   (   4.323    4.323    0.000)    6.113
  13.737   (   1.933    1.933    0.000)    2.734
  14.364   (   3.682    3.682    0.000)    5.208
  14.437   (  -3.202   -3.202    0.000)    4.528
  20.451   (   2.168    2.168    0.000)    3.066
  34.884   (  -0.149   -0.149    0.000)    0.210
  35.415   (   1.377    1.377    0.000)    1.948
  39.563   (   5.386    5.386    0.000)    7.617
======================= Grid point 133 (22/56) =======================
q-point: ( 0.09  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.989   (  13.496   13.496   13.496)   23.376
   0.989   (  13.496   13.496   13.496)   23.376
   1.507   (  19.632   19.632   19.632)   34.004
   3.501   (   7.711    7.711    7.711)   13.357
   3.501   (   7.711    7.711    7.711)   13.357
   3.662   (   4.501    4.501    4.501)    7.796
  10.926   ( -10.089  -10.089  -10.089)   17.474
  11.784   (  11.059   11.059   11.059)   19.155
  11.784   (  11.059   11.059   11.059)   19.155
  12.994   (  10.081   10.081   10.081)   17.461
  12.994   (  10.081   10.081   10.081)   17.461
  19.829   (   4.967    4.967    4.967)    8.604
  34.972   (   0.688    0.688    0.688)    1.192
  34.972   (   0.688    0.688    0.688)    1.192
  37.792   (   4.567    4.567    4.567)    7.911
======================= Grid point 134 (23/56) =======================
q-point: ( 0.18  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.271   (  14.533    8.558    8.558)   18.913
   1.295   (  14.253    8.094    8.094)   18.281
   1.915   (  17.553    7.995    7.995)   20.880
   3.521   (   0.984   11.652   11.652)   16.508
   3.561   (   1.066   13.002   13.002)   18.418
   4.159   (  32.426   -1.156   -1.156)   32.467
  10.822   (   0.020  -21.627  -21.627)   30.586
  11.933   (  11.885    1.296    1.296)   12.026
  12.135   (  13.203   -3.106   -3.106)   13.914
  13.147   (   7.505   20.367   20.367)   29.765
  13.389   (  25.050   22.469   22.469)   40.462
  19.941   (   5.715    3.605    3.605)    7.658
  34.972   (  -0.033    1.299    1.299)    1.838
  34.990   (   0.339    1.655    1.655)    2.365
  37.943   (   9.768    4.671    4.671)   11.792
======================= Grid point 135 (24/56) =======================
q-point: ( 0.27  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.549   (  12.203    4.547    4.547)   13.794
   1.558   (  11.154    4.209    4.209)   12.643
   2.194   (  10.124   -0.045   -0.045)   10.124
   3.536   (   0.460   12.351   12.351)   17.473
   3.611   (   3.301   15.615   15.615)   22.328
   4.822   (  30.589    0.257    0.257)   30.591
  10.907   (   7.431  -32.875  -32.875)   47.082
  12.217   (  14.737   -9.413   -9.413)   19.859
  12.333   (   6.761  -10.394  -10.394)   16.179
  13.300   (   7.286   26.259   26.259)   37.843
  13.942   (  27.059   32.181   32.181)   52.947
  20.056   (   5.376    0.893    0.893)    5.522
  34.971   (  -0.005    1.631    1.631)    2.307
  34.991   (  -0.135    2.369    2.369)    3.353
  38.165   (  11.280    5.705    5.705)   13.868
======================= Grid point 136 (25/56) =======================
q-point: ( 0.36  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.743   (   6.871    1.312    1.312)    7.117
   1.757   (   7.933    2.147    2.147)    8.495
   2.347   (   5.183   -4.671   -4.671)    8.397
   3.539   (  -0.134   13.040   13.040)   18.441
   3.677   (   2.897   18.219   18.219)   25.928
   5.354   (  20.612    0.863    0.863)   20.648
  11.083   (   8.697  -41.185  -41.185)   58.890
  12.430   (   2.983  -13.711  -13.711)   19.619
  12.493   (  11.387  -17.217  -17.217)   26.879
  13.435   (   5.557   30.660   30.660)   43.714
  14.430   (  19.508   36.494   36.494)   55.174
  20.152   (   3.848   -1.937   -1.937)    4.723
  34.971   (   0.018    1.939    1.939)    2.743
  34.987   (  -0.180    2.643    2.643)    3.742
  38.376   (   8.797    6.939    6.939)   13.179
======================= Grid point 137 (26/56) =======================
q-point: ( 0.45  0.09  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.837   (   2.290   -0.352   -0.352)    2.343
   1.867   (   2.759    0.979    0.979)    3.087
   2.415   (   1.629   -6.596   -6.596)    9.469
   3.535   (  -0.168   13.476   13.476)   19.059
   3.719   (   1.088   19.691   19.691)   27.869
   5.641   (   7.169    1.042    1.042)    7.319
  11.219   (   3.817  -45.627  -45.627)   64.639
  12.465   (   0.727  -14.856  -14.856)   21.023
  12.657   (   4.221  -21.569  -21.569)   30.794
  13.515   (   2.098   33.261   33.261)   47.085
  14.706   (   7.009   38.061   38.061)   54.281
  20.207   (   1.394   -3.663   -3.663)    5.365
  34.972   (   0.011    2.124    2.124)    3.004
  34.985   (  -0.071    2.735    2.735)    3.868
  38.504   (   3.307    7.704    7.704)   11.386
======================= Grid point 145 (27/56) =======================
q-point: ( 0.18  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.461   (   9.751    9.751    6.116)   15.085
   1.484   (   9.303    9.303    4.241)   13.823
   2.079   (   6.919    6.919    2.891)   10.203
   3.591   (   3.311    3.311   22.875)   23.349
   4.112   (  17.539   17.539    1.186)   24.833
   4.203   (  18.229   18.229    1.820)   25.843
  10.396   ( -18.334  -18.334  -18.607)   31.913
  11.862   (   0.097    0.097  -11.394)   11.394
  12.258   (   6.948    6.948  -35.295)   36.637
  13.620   (  19.867   19.867   38.274)   47.479
  13.883   (  28.799   28.799   34.323)   53.262
  20.017   (   3.455    3.455    0.920)    4.972
  34.977   (  -0.071   -0.071    3.549)    3.551
  35.043   (   2.429    2.429    0.304)    3.449
  38.100   (  10.449   10.449    4.781)   15.532
======================= Grid point 146 (28/56) =======================
q-point: ( 0.27  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.648   (   7.993    5.130    3.489)   10.118
   1.662   (   8.044    5.179    2.585)    9.910
   2.185   (   3.626   -0.932   -1.633)    4.085
   3.660   (   3.115    5.593   23.466)   24.324
   4.215   (   4.158   34.408    3.319)   34.817
   4.848   (  31.094    2.631    1.331)   31.233
  10.107   (  -9.842  -41.756  -17.556)   46.353
  11.953   (   7.534   -6.117  -28.511)   30.117
  12.373   (   4.418    3.622  -45.646)   46.002
  13.924   (  11.472   30.867   44.740)   55.552
  14.515   (  30.021   24.223   40.930)   56.243
  20.077   (   2.444    1.083   -2.305)    3.530
  34.974   (  -0.153   -0.233    4.305)    4.314
  35.081   (   1.370    4.782    0.315)    4.984
  38.345   (  12.746   11.720    5.209)   18.082
======================= Grid point 147 (29/56) =======================
q-point: ( 0.36  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.775   (   4.450    1.754    1.129)    4.914
   1.807   (   5.758    2.395    1.946)    6.533
   2.240   (   1.962   -5.424   -4.005)    7.021
   3.710   (   1.729    7.508   23.620)   24.845
   4.300   (   3.679   35.454    5.073)   36.004
   5.385   (  20.630    2.229    0.952)   20.772
   9.979   (  -3.291  -61.299  -16.753)   63.632
  12.123   (   8.157   -7.454  -38.122)   39.691
  12.440   (   2.056    0.622  -52.576)   52.620
  14.113   (   7.085   33.221   47.082)   58.056
  15.040   (  20.417   22.425   44.021)   53.457
  20.117   (   1.525   -0.949   -5.374)    5.666
  34.971   (  -0.140   -0.391    4.864)    4.882
  35.102   (   0.713    6.374    0.109)    6.415
  38.587   (  10.118   13.323    5.860)   17.726
======================= Grid point 148 (30/56) =======================
q-point: ( 0.45  0.18  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.833   (   1.382    0.006   -0.252)    1.404
   1.888   (   2.043    0.974    1.552)    2.745
   2.267   (   0.713   -7.420   -4.713)    8.819
   3.732   (   0.493    8.658   23.490)   25.040
   4.353   (   1.411   36.154    6.139)   36.698
   5.671   (   7.110    1.780    0.733)    7.366
   9.946   (  -0.513  -72.099  -15.972)   73.849
  12.252   (   3.763   -5.448  -41.156)   41.686
  12.460   (   0.167   -2.902  -57.617)   57.691
  14.211   (   2.457   33.404   47.619)   58.219
  15.325   (   7.139   21.244   45.376)   50.609
  20.139   (   0.536   -2.055   -7.257)    7.562
  34.969   (  -0.056   -0.487    5.150)    5.173
  35.111   (   0.232    7.147   -0.067)    7.151
  38.733   (   3.816   14.386    6.312)   16.167
======================= Grid point 157 (31/56) =======================
q-point: ( 0.27  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.749   (   4.283    4.283    2.715)    6.638
   1.760   (   4.592    4.592    1.901)    6.767
   2.160   (  -1.338   -1.338   -2.893)    3.457
   3.773   (   5.109    5.109   23.962)   25.027
   4.837   (  16.967   16.967    0.744)   24.006
   4.987   (  18.443   18.443    3.063)   26.262
   9.312   ( -34.327  -34.327  -11.659)   49.926
  11.949   (   4.334    4.334  -50.015)   50.389
  12.435   (   2.245    2.245  -57.121)   57.209
  14.452   (  19.470   19.470   53.098)   59.812
  14.972   (  21.161   21.161   46.662)   55.434
  20.100   (   1.193    1.193   -6.195)    6.421
  34.968   (  -0.326   -0.326    4.818)    4.840
  35.167   (   3.503    3.503    0.240)    4.960
  38.621   (  14.420   14.420    4.812)   20.954
======================= Grid point 158 (32/56) =======================
q-point: ( 0.36  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.812   (   2.007    1.743    1.565)    3.085
   1.850   (   3.847    1.885    1.788)    4.642
   2.136   (  -0.867   -4.506   -3.018)    5.492
   3.859   (   3.187    6.718   24.199)   25.316
   4.993   (   5.106   30.270    2.133)   30.771
   5.452   (  19.147    4.741    1.152)   19.759
   8.743   ( -21.016  -57.065   -9.471)   61.545
  12.081   (   7.546    2.023  -61.143)   61.641
  12.457   (  -0.216    0.603  -65.316)   65.319
  14.728   (   8.524   25.869   55.529)   61.849
  15.360   (  15.711    9.912   51.290)   54.551
  20.125   (   1.188    1.673   -9.615)    9.831
  34.962   (  -0.274   -0.479    5.132)    5.162
  35.223   (   2.024    5.230   -0.018)    5.608
  38.896   (  11.605   16.083    4.707)   20.384
======================= Grid point 159 (33/56) =======================
q-point: ( 0.45  0.27  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.836   (   0.512    0.232    0.516)    0.763
   1.907   (   1.507    0.912    1.733)    2.471
   2.125   (  -0.250   -6.256   -2.482)    6.735
   3.902   (   1.034    7.635   24.231)   25.426
   5.071   (   2.130   32.220    2.751)   32.408
   5.713   (   6.293    2.376    0.304)    6.733
   8.459   (  -6.845  -70.491   -8.484)   71.329
  12.219   (   4.796    1.280  -63.481)   63.675
  12.431   (  -1.776   -0.892  -72.010)   72.037
  14.834   (   2.385   26.758   55.762)   61.895
  15.582   (   5.602    5.062   53.820)   54.347
  20.144   (   0.533    2.482  -11.581)   11.856
  34.958   (  -0.106   -0.571    5.268)    5.300
  35.251   (   0.665    6.126   -0.244)    6.167
  39.065   (   4.401   17.147    4.738)   18.325
======================= Grid point 169 (34/56) =======================
q-point: ( 0.36  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.837   (   0.589    0.589    2.035)    2.199
   1.889   (   2.003    2.003    1.333)    3.131
   2.060   (  -2.864   -2.864   -1.612)    4.359
   3.973   (   4.326    4.326   24.653)   25.401
   5.441   (  12.131   12.131   -0.053)   17.157
   5.618   (  11.860   11.860    1.479)   16.838
   7.731   ( -40.136  -40.136   -7.228)   57.219
  12.164   (   5.383    5.383  -73.937)   74.328
  12.452   (  -1.297   -1.297  -74.027)   74.050
  15.135   (  13.319   13.319   59.500)   62.410
  15.501   (   5.369    5.369   56.215)   56.726
  20.174   (   2.724    2.724  -12.927)   13.488
  34.953   (  -0.389   -0.389    5.236)    5.265
  35.310   (   3.195    3.195   -0.337)    4.531
  39.203   (  12.951   12.951    3.912)   18.729
======================= Grid point 170 (35/56) =======================
q-point: ( 0.45  0.36  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.840   (  -0.035    0.139    1.579)    1.585
   1.925   (   1.199    0.862    1.451)    2.070
   2.019   (  -1.159   -4.012   -0.401)    4.195
   4.032   (   1.449    4.977   24.961)   25.494
   5.609   (   5.655   18.688    0.363)   19.528
   5.782   (   2.359    5.560   -0.082)    6.040
   7.167   ( -14.073  -53.706   -6.263)   55.872
  12.272   (   4.326    3.282  -77.966)   78.155
  12.407   (  -2.603   -1.495  -80.101)   80.157
  15.283   (   2.554   16.783   59.991)   62.346
  15.598   (   3.149   -1.882   59.280)   59.394
  20.219   (   1.266    4.102  -14.587)   15.205
  34.947   (  -0.149   -0.456    5.247)    5.269
  35.353   (   1.065    3.836   -0.608)    4.027
  39.391   (   4.922   13.776    3.547)   15.053
======================= Grid point 181 (36/56) =======================
q-point: ( 0.45  0.45  0.09)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.841   (  -0.015   -0.015    2.176)    2.176
   1.940   (   0.603    0.603    1.218)    1.487
   1.963   (  -1.521   -1.521    0.909)    2.335
   4.101   (   1.691    1.691   25.565)   25.677
   5.825   (   7.480    7.480   -1.314)   10.659
   5.938   (   3.806    3.806   -0.500)    5.406
   6.354   ( -23.185  -23.185   -5.033)   33.173
  12.332   (   2.324    2.324  -84.385)   84.449
  12.374   (  -1.609   -1.609  -84.283)   84.314
  15.508   (   4.680    4.680   61.684)   62.038
  15.561   (  -0.405   -0.405   61.346)   61.349
  20.284   (   1.850    1.850  -15.920)   16.133
  34.940   (  -0.174   -0.174    5.183)    5.189
  35.406   (   1.289    1.289   -0.907)    2.036
  39.592   (   5.237    5.237    2.948)    7.971
======================= Grid point 266 (37/56) =======================
q-point: ( 0.18  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.599   (   6.722    6.722    6.722)   11.642
   1.599   (   6.722    6.722    6.722)   11.642
   2.128   (   1.465    1.465    1.465)    2.537
   4.145   (  12.351   12.351   12.351)   21.393
   4.145   (  12.351   12.351   12.351)   21.393
   4.326   (  16.232   16.232   16.232)   28.115
   9.958   ( -21.510  -21.510  -21.510)   37.257
  11.700   ( -12.085  -12.085  -12.085)   20.932
  11.700   ( -12.085  -12.085  -12.085)   20.932
  14.459   (  31.836   31.836   31.836)   55.142
  14.459   (  31.836   31.836   31.836)   55.142
  20.027   (  -0.106   -0.106   -0.106)    0.183
  35.052   (   1.894    1.894    1.894)    3.281
  35.052   (   1.894    1.894    1.894)    3.281
  38.258   (  10.606   10.606   10.606)   18.370
======================= Grid point 267 (38/56) =======================
q-point: ( 0.27  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.726   (   5.296    4.105    4.105)    7.858
   1.737   (   6.562    4.250    4.250)    8.899
   2.142   (   0.265   -2.406   -2.406)    3.413
   4.183   (   1.738   17.896   17.896)   25.368
   4.378   (   8.170   19.176   19.176)   28.323
   4.903   (  29.753    4.561    4.561)   30.444
   9.588   ( -13.852  -31.225  -31.225)   46.281
  11.452   (  -6.197  -23.050  -23.050)   33.181
  11.623   (  -2.990  -23.753  -23.753)   33.725
  14.822   (  13.677   40.698   40.698)   59.158
  15.240   (  34.071   29.988   29.988)   54.400
  20.012   (  -1.265   -3.860   -3.860)    5.604
  35.064   (   0.617    2.760    2.760)    3.951
  35.101   (   1.562    2.623    2.623)    4.026
  38.512   (  13.367   10.974   10.974)   20.483
======================= Grid point 268 (39/56) =======================
q-point: ( 0.36  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.807   (   2.572    2.021    2.021)    3.845
   1.854   (   4.672    2.461    2.461)    5.826
   2.151   (   0.678   -4.353   -4.353)    6.194
   4.211   (   0.970   18.410   18.410)   26.054
   4.528   (   6.146   21.284   21.284)   30.721
   5.418   (  19.768    2.446    2.446)   20.068
   9.394   (  -5.629  -40.667  -40.667)   57.786
  11.400   (  -0.002  -32.900  -32.900)   46.527
  11.568   (  -2.366  -28.525  -28.525)   40.410
  15.035   (   7.408   41.714   41.714)   59.455
  15.820   (  22.069   30.821   30.821)   48.857
  19.983   (  -1.337   -6.938   -6.938)    9.902
  35.075   (   0.494    2.928    2.928)    4.170
  35.120   (   0.497    3.429    3.429)    4.875
  38.767   (  10.820   11.677   11.677)   19.743
======================= Grid point 269 (40/56) =======================
q-point: ( 0.45  0.18  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.837   (   0.591    0.624    0.624)    1.063
   1.921   (   1.681    1.569    1.569)    2.784
   2.165   (   0.495   -4.911   -4.911)    6.963
   4.224   (   0.262   18.784   18.784)   26.566
   4.616   (   2.255   22.447   22.447)   31.825
   5.692   (   6.787    1.252    1.252)    7.014
   9.332   (  -1.138  -46.011  -46.011)   65.079
  11.414   (   0.708  -37.566  -37.566)   53.131
  11.532   (  -1.016  -30.970  -30.970)   43.810
  15.133   (   2.333   41.503   41.503)   58.740
  16.125   (   7.612   30.925   30.925)   44.393
  19.963   (  -0.538   -8.653   -8.653)   12.249
  35.083   (   0.189    3.026    3.026)    4.284
  35.125   (   0.093    3.720    3.720)    5.261
  38.925   (   4.110   12.229   12.229)   17.776
======================= Grid point 278 (41/56) =======================
q-point: ( 0.27  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.814   (   3.852    3.852    3.212)    6.325
   1.816   (   3.647    3.647    3.608)    6.294
   2.095   (  -1.949   -1.949   -3.237)    4.251
   4.360   (   5.438    5.438   30.321)   31.281
   4.860   (  16.507   16.507    1.504)   23.393
   5.098   (  17.747   17.747    8.317)   26.440
   8.938   ( -30.033  -30.033  -24.746)   49.156
  11.047   ( -13.987  -13.987  -36.381)   41.410
  11.326   (  -9.721   -9.721  -49.336)   51.216
  15.487   (  22.321   22.321   46.695)   56.364
  15.800   (  26.651   26.651   33.339)   50.319
  19.930   (  -3.856   -3.856   -9.835)   11.245
  35.084   (   0.346    0.346    5.776)    5.796
  35.174   (   3.250    3.250    0.416)    4.615
  38.774   (  13.969   13.969   10.166)   22.217
======================= Grid point 279 (42/56) =======================
q-point: ( 0.36  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.854   (   0.757    2.476    2.516)    3.610
   1.898   (   3.466    1.614    2.686)    4.672
   2.072   (  -0.230   -3.096   -2.974)    4.299
   4.447   (   3.101    8.243   29.920)   31.190
   5.068   (   7.284   25.361    5.556)   26.965
   5.497   (  16.591    5.888    3.666)   17.983
   8.427   ( -19.123  -49.756  -22.091)   57.700
  10.882   (  -3.003  -17.409  -54.104)   56.915
  11.158   (  -7.256  -14.552  -58.849)   61.054
  15.777   (   8.281   30.221   46.117)   55.755
  16.291   (  19.563   15.992   38.381)   45.951
  19.864   (  -2.485   -4.274  -14.930)   15.727
  35.088   (   0.109   -0.126    6.453)    6.456
  35.223   (   1.644    4.813    0.058)    5.087
  39.044   (  11.505   14.645    9.746)   21.020
======================= Grid point 280 (43/56) =======================
q-point: ( 0.45  0.27  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.857   (  -0.117    1.302    1.563)    2.037
   1.949   (   1.332    1.104    2.173)    2.778
   2.076   (   0.252   -3.702   -2.123)    4.275
   4.488   (   0.967    9.856   29.505)   31.122
   5.177   (   2.957   27.577    8.579)   29.032
   5.722   (   5.372    1.805    0.607)    5.700
   8.169   (  -6.199  -62.340  -21.022)   66.081
  10.885   (   2.020  -17.059  -66.973)   69.141
  11.029   (  -4.767  -18.692  -59.866)   62.898
  15.877   (   2.109   30.822   45.427)   54.937
  16.564   (   6.788   12.391   40.537)   42.929
  19.830   (  -0.791   -3.705  -17.781)   18.180
  35.089   (   0.016   -0.436    6.808)    6.822
  35.244   (   0.486    5.634   -0.293)    5.663
  39.213   (   4.408   15.139    9.643)   18.483
======================= Grid point 290 (44/56) =======================
q-point: ( 0.36  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.891   (   0.406    0.406    3.025)    3.079
   1.927   (   1.840    1.840    2.279)    3.459
   2.031   (  -1.056   -1.056   -1.141)    1.880
   4.587   (   5.344    5.344   32.131)   33.008
   5.440   (  11.443   11.443    0.019)   16.183
   5.677   (  10.098   10.098    4.822)   15.072
   7.507   ( -37.556  -37.556  -14.803)   55.137
  10.659   (  -5.223   -5.223  -71.231)   71.613
  10.915   (  -9.915   -9.915  -74.115)   75.429
  16.239   (  14.150   14.150   47.759)   51.782
  16.529   (   9.086    9.086   43.132)   45.005
  19.813   (  -0.961   -0.961  -21.376)   21.419
  35.083   (  -0.273   -0.273    6.813)    6.824
  35.300   (   2.717    2.717   -0.564)    3.884
  39.328   (  12.107   12.107    8.294)   19.025
======================= Grid point 291 (45/56) =======================
q-point: ( 0.45  0.36  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.885   (  -0.377    1.214    2.652)    2.941
   1.965   (   1.340    0.532    2.307)    2.720
   2.017   (  -0.423   -1.967    0.314)    2.036
   4.661   (   1.852    6.945   33.323)   34.090
   5.626   (   7.352   14.074    1.526)   15.952
   5.785   (  -0.423    5.876    0.778)    5.943
   6.977   ( -13.186  -50.802  -12.342)   53.917
  10.638   (   2.388   -7.788  -80.686)   81.096
  10.743   (  -6.807  -10.138  -80.051)   80.977
  16.382   (   2.144   18.551   47.016)   50.589
  16.682   (   4.606    0.920   45.872)   46.112
  19.808   (   0.111    0.858  -24.445)   24.460
  35.079   (  -0.130   -0.500    6.928)    6.947
  35.336   (   0.849    3.301   -1.018)    3.558
  39.505   (   4.659   12.494    7.585)   15.341
======================= Grid point 302 (46/56) =======================
q-point: ( 0.45  0.45  0.18)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.901   (   0.350    0.350    3.476)    3.511
   1.975   (   0.575    0.575    2.079)    2.232
   1.991   (  -0.845   -0.845    1.769)    2.135
   4.761   (   2.604    2.604   36.237)   36.423
   5.781   (   5.649    5.649   -3.069)    8.558
   5.924   (   2.490    2.490   -0.752)    3.601
   6.212   ( -21.055  -21.055   -8.523)   30.972
  10.553   (  -0.837   -0.837  -88.065)   88.073
  10.593   (  -4.657   -4.657  -88.212)   88.458
  16.628   (   4.798    4.798   47.451)   47.934
  16.675   (   0.378    0.378   47.111)   47.114
  19.833   (   0.959    0.959  -27.072)   27.106
  35.071   (  -0.214   -0.214    6.899)    6.906
  35.380   (   1.052    1.052   -1.527)    2.131
  39.688   (   4.812    4.812    6.516)    9.422
======================= Grid point 399 (47/56) =======================
q-point: ( 0.27  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.885   (   3.133    3.133    3.133)    5.427
   1.885   (   3.133    3.133    3.133)    5.427
   2.036   (  -2.191   -2.191   -2.191)    3.795
   4.898   (  11.435   11.435   11.435)   19.806
   4.898   (  11.435   11.435   11.435)   19.806
   5.344   (  15.313   15.313   15.313)   26.523
   8.362   ( -29.226  -29.226  -29.226)   50.621
  10.468   ( -27.763  -27.763  -27.763)   48.088
  10.468   ( -27.763  -27.763  -27.763)   48.088
  16.327   (  26.839   26.839   26.839)   46.486
  16.327   (  26.839   26.839   26.839)   46.486
  19.726   (  -9.498   -9.498   -9.498)   16.451
  35.183   (   2.097    2.097    2.097)    3.632
  35.183   (   2.097    2.097    2.097)    3.632
  39.021   (  13.275   13.275   13.275)   22.992
======================= Grid point 400 (48/56) =======================
q-point: ( 0.36  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.911   (  -0.269    2.831    2.831)    4.012
   1.950   (   2.712    2.156    2.156)    4.081
   2.021   (   0.767   -1.769   -1.769)    2.616
   4.941   (   1.972   16.245   16.245)   23.058
   5.235   (  13.506   10.366   10.366)   19.933
   5.633   (  10.544   10.202   10.202)   17.870
   7.835   ( -20.728  -34.601  -34.601)   53.142
   9.938   ( -18.036  -37.598  -37.598)   56.147
  10.128   ( -12.393  -40.842  -40.842)   59.073
  16.586   (   7.561   32.535   32.535)   46.629
  16.908   (  21.539   22.132   22.132)   37.995
  19.563   (  -6.067  -13.270  -13.270)   19.723
  35.192   (   0.373    2.589    2.589)    3.680
  35.236   (   1.545    1.964    1.964)    3.178
  39.281   (  11.203   12.768   12.768)   21.249
======================= Grid point 401 (49/56) =======================
q-point: ( 0.45  0.27  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.898   (  -0.505    2.259    2.259)    3.235
   1.990   (   1.017    1.718    1.718)    2.634
   2.042   (   0.713   -1.137   -1.137)    1.759
   4.971   (   0.813   16.775   16.775)   23.737
   5.454   (   6.810   16.766   16.766)   24.669
   5.745   (   1.327    2.354    2.354)    3.584
   7.549   (  -7.071  -39.843  -39.843)   56.788
   9.707   (  -5.198  -44.578  -44.578)   63.256
   9.956   (  -4.369  -45.262  -45.262)   64.159
  16.672   (   1.683   32.276   32.276)   45.676
  17.203   (   7.240   21.576   21.576)   31.360
  19.482   (  -1.940  -14.649  -14.649)   20.808
  35.198   (   0.142    2.394    2.394)    3.389
  35.253   (   0.357    2.241    2.241)    3.189
  39.446   (   4.347   12.568   12.568)   18.298
======================= Grid point 411 (50/56) =======================
q-point: ( 0.36  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.953   (   0.195    0.195    2.842)    2.856
   1.977   (   1.625    1.625    2.530)    3.418
   2.016   (   0.968    0.968   -0.264)    1.394
   5.157   (   6.819    6.819   21.321)   23.401
   5.444   (  10.209   10.209    0.465)   14.445
   5.849   (   7.723    7.723   12.481)   16.585
   7.124   ( -31.686  -31.686  -22.626)   50.199
   9.338   ( -21.239  -21.239  -56.064)   63.604
   9.499   ( -21.606  -21.606  -62.706)   69.754
  17.073   (  13.645   13.645   33.749)   38.876
  17.244   (  13.228   13.228   26.694)   32.597
  19.357   (  -6.684   -6.684  -21.695)   23.664
  35.206   (  -0.025   -0.025    4.842)    4.842
  35.288   (   2.088    2.088   -0.611)    3.015
  39.532   (  10.909   10.909   11.203)   19.067
======================= Grid point 412 (51/56) =======================
q-point: ( 0.45  0.36  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 666
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.944   (  -0.362    2.095    2.993)    3.671
   2.013   (   1.248    0.498    2.253)    2.623
   2.031   (   0.220    0.038    0.908)    0.935
   5.249   (   2.302   12.046   20.712)   24.071
   5.664   (  10.627    3.115    1.057)   11.124
   5.873   (  -5.106    7.967   10.530)   14.156
   6.658   ( -12.000  -42.255  -19.111)   47.903
   9.066   (  -5.968  -20.276  -69.957)   73.080
   9.207   (  -7.236  -27.311  -68.895)   74.463
  17.195   (   1.490   18.891   32.491)   37.613
  17.456   (   5.838    4.570   29.754)   30.663
  19.279   (  -1.519   -5.144  -25.662)   26.217
  35.204   (  -0.070   -0.397    4.951)    4.968
  35.314   (   0.584    2.500   -1.034)    2.768
  39.693   (   4.268   10.768   10.383)   15.556
======================= Grid point 423 (52/56) =======================
q-point: ( 0.45  0.45  0.27)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.975   (   0.799    0.799    3.568)    3.743
   2.021   (   0.533    0.533    2.307)    2.427
   2.031   (  -0.256   -0.256    2.058)    2.089
   5.459   (   6.923    6.923   27.699)   29.378
   5.695   (   0.818    0.818   -5.405)    5.527
   5.930   (  -2.689   -2.689    3.574)    5.219
   6.021   ( -15.526  -15.526   -9.997)   24.126
   8.810   (  -5.586   -5.586  -80.220)   80.608
   8.844   (  -8.757   -8.757  -80.777)   81.720
  17.442   (   4.485    4.485   32.349)   32.965
  17.481   (   0.940    0.940   31.917)   31.945
  19.235   (  -0.340   -0.340  -29.889)   29.892
  35.197   (  -0.225   -0.225    5.035)    5.045
  35.347   (   0.740    0.740   -1.663)    1.965
  39.853   (   4.225    4.225    9.198)   10.968
======================= Grid point 532 (53/56) =======================
q-point: ( 0.36  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.994   (   0.901    0.901    0.901)    1.561
   2.025   (   1.604    1.604    1.604)    2.778
   2.025   (   1.604    1.604    1.604)    2.778
   5.464   (   6.253    6.253    6.253)   10.830
   5.464   (   6.253    6.253    6.253)   10.830
   6.110   (   9.423    9.423    9.423)   16.321
   6.614   ( -25.510  -25.510  -25.510)   44.185
   8.430   ( -36.128  -36.128  -36.128)   62.576
   8.430   ( -36.128  -36.128  -36.128)   62.576
  17.619   (  15.081   15.081   15.081)   26.122
  17.619   (  15.081   15.081   15.081)   26.122
  18.991   ( -12.842  -12.842  -12.842)   22.244
  35.277   (   0.849    0.849    0.849)    1.470
  35.277   (   0.849    0.849    0.849)    1.470
  39.755   (   9.824    9.824    9.824)   17.015
======================= Grid point 533 (54/56) =======================
q-point: ( 0.45  0.36  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.001   (   0.112    2.503    2.503)    3.542
   2.049   (   0.719    1.037    1.037)    1.633
   2.053   (   0.739    1.300    1.300)    1.981
   5.508   (   2.099    8.059    8.059)   11.588
   5.652   (   6.756   -1.882   -1.882)    7.261
   6.161   (  -2.569   12.042   12.042)   17.222
   6.191   ( -12.741  -25.453  -25.453)   38.184
   7.912   ( -11.338  -40.664  -40.664)   58.614
   8.004   ( -10.037  -46.294  -46.294)   66.234
  17.715   (   1.003   18.660   18.660)   26.408
  17.887   (   6.504   12.051   12.051)   18.241
  18.839   (  -2.902  -15.338  -15.338)   21.885
  35.270   (  -0.173    1.157    1.157)    1.645
  35.302   (   0.514    0.281    0.281)    0.650
  39.902   (   3.897    9.219    9.219)   13.608
======================= Grid point 544 (55/56) =======================
q-point: ( 0.45  0.45  0.36)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 366
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.040   (   1.105    1.105    2.657)    3.082
   2.064   (   0.493    0.493    1.792)    1.923
   2.069   (   0.232    0.232    1.556)    1.590
   5.571   (  -5.723   -5.723   -6.058)   10.110
   5.657   (   5.350    5.350   -4.775)    8.947
   5.787   (  -8.636   -8.636  -13.789)   18.420
   6.225   (  -4.019   -4.019   19.864)   20.661
   7.401   ( -10.951  -10.951  -53.843)   56.026
   7.412   ( -12.177  -12.177  -55.656)   58.259
  17.958   (   4.057    4.057   18.440)   19.312
  17.987   (   1.834    1.834   17.809)   17.997
  18.679   (  -2.057   -2.057  -22.913)   23.097
  35.271   (  -0.220   -0.220    2.274)    2.295
  35.316   (   0.453    0.453   -1.271)    1.423
  40.040   (   3.683    3.683    8.345)    9.837
======================= Grid point 665 (56/56) =======================
q-point: ( 0.45  0.45  0.45)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 2.98e-04 1.10e-04 
Number of triplets: 146
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.080   (   1.151    1.151    1.151)    1.993
   2.090   (   0.538    0.538    0.538)    0.931
   2.090   (   0.538    0.538    0.538)    0.931
   5.474   (  -6.597   -6.597   -6.597)   11.426
   5.474   (  -6.597   -6.597   -6.597)   11.426
   5.486   (  -9.856   -9.856   -9.856)   17.072
   6.492   (   3.070    3.070    3.070)    5.318
   6.713   ( -12.872  -12.872  -12.872)   22.295
   6.713   ( -12.872  -12.872  -12.872)   22.295
  18.205   (   4.495    4.495    4.495)    7.785
  18.205   (   4.495    4.495    4.495)    7.785
  18.367   (  -5.968   -5.968   -5.968)   10.337
  35.298   (   0.030    0.030    0.030)    0.053
  35.298   (   0.030    0.030    0.030)    0.053
  40.165   (   3.370    3.370    3.370)    5.838
=================== End of collection of collisions ===================
----------- Thermal conductivity (W/m-k) with tetrahedron method -----------
#  T(K)        xx         yy         zz         yz         xz         xy        #ipm
    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/19965
   10.0    417.982    417.982    417.982      0.000      0.000      0.000 3/19965
   20.0    197.513    197.513    197.513      0.000      0.000      0.000 3/19965
   30.0    109.076    109.076    109.076      0.000      0.000      0.000 3/19965
   40.0     56.776     56.776     56.776      0.000      0.000      0.000 3/19965
   50.0     32.982     32.982     32.982      0.000      0.000      0.000 3/19965
   60.0     21.843     21.843     21.843      0.000      0.000      0.000 3/19965
   70.0     15.941     15.941     15.941      0.000      0.000      0.000 3/19965
   80.0     12.427     12.427     12.427      0.000      0.000      0.000 3/19965
   90.0     10.138     10.138     10.138     -0.000      0.000      0.000 3/19965
  100.0      8.547      8.547      8.547     -0.000      0.000      0.000 3/19965
  110.0      7.387      7.387      7.387     -0.000      0.000      0.000 3/19965
  120.0      6.509      6.509      6.509     -0.000      0.000      0.000 3/19965
  130.0      5.825      5.825      5.825     -0.000      0.000      0.000 3/19965
  140.0      5.280      5.280      5.280     -0.000      0.000      0.000 3/19965
  150.0      4.835      4.835      4.835     -0.000      0.000      0.000 3/19965
  160.0      4.467      4.467      4.467     -0.000      0.000      0.000 3/19965
  170.0      4.156      4.156      4.156     -0.000      0.000      0.000 3/19965
  180.0      3.891      3.891      3.891     -0.000      0.000      0.000 3/19965
  190.0      3.661      3.661      3.661     -0.000      0.000      0.000 3/19965
  200.0      3.460      3.460      3.460     -0.000      0.000      0.000 3/19965
  210.0      3.282      3.282      3.282     -0.000      0.000      0.000 3/19965
  220.0      3.123      3.123      3.123     -0.000      0.000      0.000 3/19965
  230.0      2.981      2.981      2.981     -0.000      0.000      0.000 3/19965
  240.0      2.852      2.852      2.852     -0.000      0.000      0.000 3/19965
  250.0      2.734      2.734      2.734     -0.000      0.000      0.000 3/19965
  260.0      2.626      2.626      2.626     -0.000      0.000      0.000 3/19965
  270.0      2.527      2.527      2.527     -0.000      0.000      0.000 3/19965
  280.0      2.436      2.436      2.436     -0.000      0.000      0.000 3/19965
  290.0      2.351      2.351      2.351     -0.000      0.000      0.000 3/19965
  300.0      2.272      2.272      2.272     -0.000      0.000      0.000 3/19965
  310.0      2.198      2.198      2.198     -0.000      0.000      0.000 3/19965
  320.0      2.129      2.129      2.129     -0.000      0.000      0.000 3/19965
  330.0      2.065      2.065      2.065     -0.000      0.000      0.000 3/19965
  340.0      2.004      2.004      2.004     -0.000      0.000      0.000 3/19965
  350.0      1.946      1.946      1.946     -0.000      0.000      0.000 3/19965
  360.0      1.892      1.892      1.892     -0.000      0.000      0.000 3/19965
  370.0      1.841      1.841      1.841     -0.000      0.000      0.000 3/19965
  380.0      1.793      1.793      1.793     -0.000      0.000      0.000 3/19965
  390.0      1.747      1.747      1.747     -0.000      0.000      0.000 3/19965
  400.0      1.703      1.703      1.703      0.000      0.000      0.000 3/19965
  410.0      1.662      1.662      1.662      0.000      0.000      0.000 3/19965
  420.0      1.622      1.622      1.622      0.000      0.000      0.000 3/19965
  430.0      1.584      1.584      1.584      0.000      0.000      0.000 3/19965
  440.0      1.548      1.548      1.548      0.000      0.000      0.000 3/19965
  450.0      1.514      1.514      1.514      0.000      0.000      0.000 3/19965
  460.0      1.481      1.481      1.481      0.000      0.000      0.000 3/19965
  470.0      1.449      1.449      1.449      0.000      0.000      0.000 3/19965
  480.0      1.419      1.419      1.419      0.000      0.000      0.000 3/19965
  490.0      1.390      1.390      1.390      0.000      0.000      0.000 3/19965
  500.0      1.362      1.362      1.362      0.000      0.000      0.000 3/19965
  510.0      1.336      1.336      1.336      0.000      0.000      0.000 3/19965
  520.0      1.310      1.310      1.310      0.000      0.000      0.000 3/19965
  530.0      1.285      1.285      1.285      0.000      0.000      0.000 3/19965
  540.0      1.262      1.262      1.262      0.000      0.000      0.000 3/19965
  550.0      1.239      1.239      1.239      0.000      0.000      0.000 3/19965
  560.0      1.217      1.217      1.217      0.000      0.000      0.000 3/19965
  570.0      1.195      1.195      1.195      0.000      0.000      0.000 3/19965
  580.0      1.175      1.175      1.175      0.000      0.000      0.000 3/19965
  590.0      1.155      1.155      1.155      0.000      0.000      0.000 3/19965
  600.0      1.136      1.136      1.136      0.000      0.000      0.000 3/19965
  610.0      1.117      1.117      1.117      0.000      0.000      0.000 3/19965
  620.0      1.099      1.099      1.099      0.000      0.000      0.000 3/19965
  630.0      1.082      1.082      1.082      0.000      0.000      0.000 3/19965
  640.0      1.065      1.065      1.065      0.000      0.000      0.000 3/19965
  650.0      1.048      1.048      1.048      0.000      0.000      0.000 3/19965
  660.0      1.032      1.032      1.032      0.000      0.000      0.000 3/19965
  670.0      1.017      1.017      1.017      0.000      0.000      0.000 3/19965
  680.0      1.002      1.002      1.002      0.000      0.000      0.000 3/19965
  690.0      0.988      0.988      0.988      0.000      0.000      0.000 3/19965
  700.0      0.974      0.974      0.974      0.000      0.000      0.000 3/19965
  710.0      0.960      0.960      0.960      0.000      0.000      0.000 3/19965
  720.0      0.947      0.947      0.947      0.000      0.000      0.000 3/19965
  730.0      0.934      0.934      0.934      0.000      0.000      0.000 3/19965
  740.0      0.921      0.921      0.921      0.000      0.000      0.000 3/19965
  750.0      0.909      0.909      0.909      0.000      0.000      0.000 3/19965
  760.0      0.897      0.897      0.897      0.000      0.000      0.000 3/19965
  770.0      0.885      0.885      0.885      0.000      0.000      0.000 3/19965
  780.0      0.874      0.874      0.874      0.000      0.000      0.000 3/19965
  790.0      0.863      0.863      0.863      0.000      0.000      0.000 3/19965
  800.0      0.852      0.852      0.852      0.000      0.000      0.000 3/19965
  810.0      0.842      0.842      0.842      0.000      0.000      0.000 3/19965
  820.0      0.831      0.831      0.831      0.000      0.000      0.000 3/19965
  830.0      0.822      0.822      0.822      0.000      0.000      0.000 3/19965
  840.0      0.812      0.812      0.812      0.000      0.000      0.000 3/19965
  850.0      0.802      0.802      0.802      0.000      0.000      0.000 3/19965
  860.0      0.793      0.793      0.793      0.000      0.000      0.000 3/19965
  870.0      0.784      0.784      0.784      0.000      0.000      0.000 3/19965
  880.0      0.775      0.775      0.775      0.000      0.000      0.000 3/19965
  890.0      0.766      0.766      0.766      0.000      0.000      0.000 3/19965
  900.0      0.758      0.758      0.758      0.000      0.000      0.000 3/19965
  910.0      0.750      0.750      0.750      0.000      0.000      0.000 3/19965
  920.0      0.741      0.741      0.741      0.000      0.000      0.000 3/19965
  930.0      0.733      0.733      0.733      0.000      0.000      0.000 3/19965
  940.0      0.726      0.726      0.726      0.000      0.000      0.000 3/19965
  950.0      0.718      0.718      0.718      0.000      0.000      0.000 3/19965
  960.0      0.711      0.711      0.711      0.000      0.000      0.000 3/19965
  970.0      0.703      0.703      0.703      0.000      0.000      0.000 3/19965
  980.0      0.696      0.696      0.696      0.000      0.000      0.000 3/19965
  990.0      0.689      0.689      0.689      0.000      0.000      0.000 3/19965
 1000.0      0.682      0.682      0.682      0.000      0.000      0.000 3/19965

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

