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

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

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

Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
Unit of length: angstrom
Settings:
  Supercell: [5 5 3]
  Primitive matrix:
    [1. 0. 0.]
    [0. 1. 0.]
    [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
Number of symmetry operations in supercell: 1800
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.977318140791256    0.000000000000000    0.000000000000000
  b   -1.988659070395628    3.444458548857920    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.514740099999999
Atomic positions (fractional):
   *1 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948
    2 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    3.977318140791256    0.000000000000000    0.000000000000000
  b   -1.988659070395628    3.444458548857920    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.514740099999999
Atomic positions (fractional):
   *1 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948 > 1
    2 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948 > 2
-------------------------------- super cell --------------------------------
Lattice vectors:
  a   19.886590703956280    0.000000000000000    0.000000000000000
  b   -9.943295351978140   17.222292744289600    0.000000000000000
  c    0.000000000000000    0.000000000000000   19.544220299999999
Atomic positions (fractional):
   *1 Ar  0.06666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    2 Ar  0.26666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    3 Ar  0.46666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    4 Ar  0.66666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    5 Ar  0.86666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    6 Ar  0.06666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    7 Ar  0.26666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    8 Ar  0.46666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    9 Ar  0.66666666666667  0.33333333333333  0.08333333333333  39.948 > 1
   10 Ar  0.86666666666667  0.33333333333333  0.08333333333333  39.948 > 1
   11 Ar  0.06666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   12 Ar  0.26666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   13 Ar  0.46666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   14 Ar  0.66666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   15 Ar  0.86666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   16 Ar  0.06666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   17 Ar  0.26666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   18 Ar  0.46666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   19 Ar  0.66666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   20 Ar  0.86666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   21 Ar  0.06666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   22 Ar  0.26666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   23 Ar  0.46666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   24 Ar  0.66666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   25 Ar  0.86666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   26 Ar  0.06666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   27 Ar  0.26666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   28 Ar  0.46666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   29 Ar  0.66666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   30 Ar  0.86666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   31 Ar  0.06666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   32 Ar  0.26666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   33 Ar  0.46666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   34 Ar  0.66666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   35 Ar  0.86666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   36 Ar  0.06666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   37 Ar  0.26666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   38 Ar  0.46666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   39 Ar  0.66666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   40 Ar  0.86666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   41 Ar  0.06666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   42 Ar  0.26666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   43 Ar  0.46666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   44 Ar  0.66666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   45 Ar  0.86666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   46 Ar  0.06666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   47 Ar  0.26666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   48 Ar  0.46666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   49 Ar  0.66666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   50 Ar  0.86666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   51 Ar  0.06666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   52 Ar  0.26666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   53 Ar  0.46666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   54 Ar  0.66666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   55 Ar  0.86666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   56 Ar  0.06666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   57 Ar  0.26666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   58 Ar  0.46666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   59 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   60 Ar  0.86666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   61 Ar  0.06666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   62 Ar  0.26666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   63 Ar  0.46666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   64 Ar  0.66666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   65 Ar  0.86666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   66 Ar  0.06666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   67 Ar  0.26666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   68 Ar  0.46666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   69 Ar  0.66666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   70 Ar  0.86666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   71 Ar  0.06666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   72 Ar  0.26666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   73 Ar  0.46666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   74 Ar  0.66666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   75 Ar  0.86666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   76 Ar  0.13333333333333  0.06666666666667  0.25000000000000  39.948 > 2
   77 Ar  0.33333333333333  0.06666666666667  0.25000000000000  39.948 > 2
   78 Ar  0.53333333333333  0.06666666666667  0.25000000000000  39.948 > 2
   79 Ar  0.73333333333333  0.06666666666667  0.25000000000000  39.948 > 2
   80 Ar  0.93333333333333  0.06666666666667  0.25000000000000  39.948 > 2
   81 Ar  0.13333333333333  0.26666666666667  0.25000000000000  39.948 > 2
   82 Ar  0.33333333333333  0.26666666666667  0.25000000000000  39.948 > 2
   83 Ar  0.53333333333333  0.26666666666667  0.25000000000000  39.948 > 2
   84 Ar  0.73333333333333  0.26666666666667  0.25000000000000  39.948 > 2
   85 Ar  0.93333333333333  0.26666666666667  0.25000000000000  39.948 > 2
   86 Ar  0.13333333333333  0.46666666666667  0.25000000000000  39.948 > 2
   87 Ar  0.33333333333333  0.46666666666667  0.25000000000000  39.948 > 2
   88 Ar  0.53333333333333  0.46666666666667  0.25000000000000  39.948 > 2
   89 Ar  0.73333333333333  0.46666666666667  0.25000000000000  39.948 > 2
   90 Ar  0.93333333333333  0.46666666666667  0.25000000000000  39.948 > 2
   91 Ar  0.13333333333333  0.66666666666667  0.25000000000000  39.948 > 2
   92 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948 > 2
   93 Ar  0.53333333333333  0.66666666666667  0.25000000000000  39.948 > 2
   94 Ar  0.73333333333333  0.66666666666667  0.25000000000000  39.948 > 2
   95 Ar  0.93333333333333  0.66666666666667  0.25000000000000  39.948 > 2
   96 Ar  0.13333333333333  0.86666666666667  0.25000000000000  39.948 > 2
   97 Ar  0.33333333333333  0.86666666666667  0.25000000000000  39.948 > 2
   98 Ar  0.53333333333333  0.86666666666667  0.25000000000000  39.948 > 2
   99 Ar  0.73333333333333  0.86666666666667  0.25000000000000  39.948 > 2
  100 Ar  0.93333333333333  0.86666666666667  0.25000000000000  39.948 > 2
  101 Ar  0.13333333333333  0.06666666666667  0.58333333333333  39.948 > 2
  102 Ar  0.33333333333333  0.06666666666667  0.58333333333333  39.948 > 2
  103 Ar  0.53333333333333  0.06666666666667  0.58333333333333  39.948 > 2
  104 Ar  0.73333333333333  0.06666666666667  0.58333333333333  39.948 > 2
  105 Ar  0.93333333333333  0.06666666666667  0.58333333333333  39.948 > 2
  106 Ar  0.13333333333333  0.26666666666667  0.58333333333333  39.948 > 2
  107 Ar  0.33333333333333  0.26666666666667  0.58333333333333  39.948 > 2
  108 Ar  0.53333333333333  0.26666666666667  0.58333333333333  39.948 > 2
  109 Ar  0.73333333333333  0.26666666666667  0.58333333333333  39.948 > 2
  110 Ar  0.93333333333333  0.26666666666667  0.58333333333333  39.948 > 2
  111 Ar  0.13333333333333  0.46666666666667  0.58333333333333  39.948 > 2
  112 Ar  0.33333333333333  0.46666666666667  0.58333333333333  39.948 > 2
  113 Ar  0.53333333333333  0.46666666666667  0.58333333333333  39.948 > 2
  114 Ar  0.73333333333333  0.46666666666667  0.58333333333333  39.948 > 2
  115 Ar  0.93333333333333  0.46666666666667  0.58333333333333  39.948 > 2
  116 Ar  0.13333333333333  0.66666666666667  0.58333333333333  39.948 > 2
  117 Ar  0.33333333333333  0.66666666666667  0.58333333333333  39.948 > 2
  118 Ar  0.53333333333333  0.66666666666667  0.58333333333333  39.948 > 2
  119 Ar  0.73333333333333  0.66666666666667  0.58333333333333  39.948 > 2
  120 Ar  0.93333333333333  0.66666666666667  0.58333333333333  39.948 > 2
  121 Ar  0.13333333333333  0.86666666666667  0.58333333333333  39.948 > 2
  122 Ar  0.33333333333333  0.86666666666667  0.58333333333333  39.948 > 2
  123 Ar  0.53333333333333  0.86666666666667  0.58333333333333  39.948 > 2
  124 Ar  0.73333333333333  0.86666666666667  0.58333333333333  39.948 > 2
  125 Ar  0.93333333333333  0.86666666666667  0.58333333333333  39.948 > 2
  126 Ar  0.13333333333333  0.06666666666667  0.91666666666667  39.948 > 2
  127 Ar  0.33333333333333  0.06666666666667  0.91666666666667  39.948 > 2
  128 Ar  0.53333333333333  0.06666666666667  0.91666666666667  39.948 > 2
  129 Ar  0.73333333333333  0.06666666666667  0.91666666666667  39.948 > 2
  130 Ar  0.93333333333333  0.06666666666667  0.91666666666667  39.948 > 2
  131 Ar  0.13333333333333  0.26666666666667  0.91666666666667  39.948 > 2
  132 Ar  0.33333333333333  0.26666666666667  0.91666666666667  39.948 > 2
  133 Ar  0.53333333333333  0.26666666666667  0.91666666666667  39.948 > 2
  134 Ar  0.73333333333333  0.26666666666667  0.91666666666667  39.948 > 2
  135 Ar  0.93333333333333  0.26666666666667  0.91666666666667  39.948 > 2
  136 Ar  0.13333333333333  0.46666666666667  0.91666666666667  39.948 > 2
  137 Ar  0.33333333333333  0.46666666666667  0.91666666666667  39.948 > 2
  138 Ar  0.53333333333333  0.46666666666667  0.91666666666667  39.948 > 2
  139 Ar  0.73333333333333  0.46666666666667  0.91666666666667  39.948 > 2
  140 Ar  0.93333333333333  0.46666666666667  0.91666666666667  39.948 > 2
  141 Ar  0.13333333333333  0.66666666666667  0.91666666666667  39.948 > 2
  142 Ar  0.33333333333333  0.66666666666667  0.91666666666667  39.948 > 2
  143 Ar  0.53333333333333  0.66666666666667  0.91666666666667  39.948 > 2
  144 Ar  0.73333333333333  0.66666666666667  0.91666666666667  39.948 > 2
  145 Ar  0.93333333333333  0.66666666666667  0.91666666666667  39.948 > 2
  146 Ar  0.13333333333333  0.86666666666667  0.91666666666667  39.948 > 2
  147 Ar  0.33333333333333  0.86666666666667  0.91666666666667  39.948 > 2
  148 Ar  0.53333333333333  0.86666666666667  0.91666666666667  39.948 > 2
  149 Ar  0.73333333333333  0.86666666666667  0.91666666666667  39.948 > 2
  150 Ar  0.93333333333333  0.86666666666667  0.91666666666667  39.948 > 2
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            1.5972440    0.0000000    0.0000000
           -0.0000000    1.5972440    0.0000000
            0.0000000    0.0000000    1.5969012
-------------------------- Born effective charges --------------------------
    1 Ar    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    2 Ar    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
----------------------------------------------------------------------------
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: 450/450
Permutation basis: 2013/2013
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 92
Number of blocks in projector: 61
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 52
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 24
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 16
Use standard eigh solver.
Tree of FC basis block matrices:
- (92, 90), data: False
|-- (16, 16), data: True
|-- (24, 23), data: True
|-- (52, 51), data: True
-----
Solver_atoms: 1 -- 150 / 150
Time (Solver_compr_matrix_reshape): 0.015
Solver_block: 80 / 80
 - Time: 0.046
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.062
--------------------------------- Symfc end --------------------------------
Max drift of force constants: 0.00000000 (yy) 0.00000000 (yy) 
Permutation basis: 450/450
Permutation basis: 2013/2013
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 92
Number of blocks in projector: 61
Finding block diagonal structure in projector.
Using scipy connected_components.
Number of blocks in projector (Sum rule): 3
--- Eigsh_solver_block: 1 / 3 ---
Block_size: 52
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 24
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 16
Use standard eigh solver.
Tree of FC basis block matrices:
- (92, 90), data: False
|-- (16, 16), data: True
|-- (24, 23), data: True
|-- (52, 51), 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 17:57:09]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 17:57:09]-------------------------
Compiled with OpenMP support (max 128 threads).
Running in phono3py.load mode.
Python version 3.14.2
Spglib version 2.6.1
----------------------------- General settings -----------------------------
Run mode: force constants
HDF5 data compression filter: gzip
Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
Supercell (dim): [5 5 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.977318140791256    0.000000000000000    0.000000000000000
  b   -1.988659070395628    3.444458548857920    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.514740099999999
Atomic positions (fractional):
    1 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948
    2 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   19.886590703956280    0.000000000000000    0.000000000000000
  b   -9.943295351978140   17.222292744289600    0.000000000000000
  c    0.000000000000000    0.000000000000000   19.544220299999999
Atomic positions (fractional):
    1 Ar  0.06666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    2 Ar  0.26666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    3 Ar  0.46666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    4 Ar  0.66666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    5 Ar  0.86666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    6 Ar  0.06666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    7 Ar  0.26666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    8 Ar  0.46666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    9 Ar  0.66666666666667  0.33333333333333  0.08333333333333  39.948 > 1
   10 Ar  0.86666666666667  0.33333333333333  0.08333333333333  39.948 > 1
   11 Ar  0.06666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   12 Ar  0.26666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   13 Ar  0.46666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   14 Ar  0.66666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   15 Ar  0.86666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   16 Ar  0.06666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   17 Ar  0.26666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   18 Ar  0.46666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   19 Ar  0.66666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   20 Ar  0.86666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   21 Ar  0.06666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   22 Ar  0.26666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   23 Ar  0.46666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   24 Ar  0.66666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   25 Ar  0.86666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   26 Ar  0.06666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   27 Ar  0.26666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   28 Ar  0.46666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   29 Ar  0.66666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   30 Ar  0.86666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   31 Ar  0.06666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   32 Ar  0.26666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   33 Ar  0.46666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   34 Ar  0.66666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   35 Ar  0.86666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   36 Ar  0.06666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   37 Ar  0.26666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   38 Ar  0.46666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   39 Ar  0.66666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   40 Ar  0.86666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   41 Ar  0.06666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   42 Ar  0.26666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   43 Ar  0.46666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   44 Ar  0.66666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   45 Ar  0.86666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   46 Ar  0.06666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   47 Ar  0.26666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   48 Ar  0.46666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   49 Ar  0.66666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   50 Ar  0.86666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   51 Ar  0.06666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   52 Ar  0.26666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   53 Ar  0.46666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   54 Ar  0.66666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   55 Ar  0.86666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   56 Ar  0.06666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   57 Ar  0.26666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   58 Ar  0.46666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   59 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   60 Ar  0.86666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   61 Ar  0.06666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   62 Ar  0.26666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   63 Ar  0.46666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   64 Ar  0.66666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   65 Ar  0.86666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   66 Ar  0.06666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   67 Ar  0.26666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   68 Ar  0.46666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   69 Ar  0.66666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   70 Ar  0.86666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   71 Ar  0.06666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   72 Ar  0.26666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   73 Ar  0.46666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   74 Ar  0.66666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   75 Ar  0.86666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   76 Ar  0.13333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   77 Ar  0.33333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   78 Ar  0.53333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   79 Ar  0.73333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   80 Ar  0.93333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   81 Ar  0.13333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   82 Ar  0.33333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   83 Ar  0.53333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   84 Ar  0.73333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   85 Ar  0.93333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   86 Ar  0.13333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   87 Ar  0.33333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   88 Ar  0.53333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   89 Ar  0.73333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   90 Ar  0.93333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   91 Ar  0.13333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   92 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   93 Ar  0.53333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   94 Ar  0.73333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   95 Ar  0.93333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   96 Ar  0.13333333333333  0.86666666666667  0.25000000000000  39.948 > 76
   97 Ar  0.33333333333333  0.86666666666667  0.25000000000000  39.948 > 76
   98 Ar  0.53333333333333  0.86666666666667  0.25000000000000  39.948 > 76
   99 Ar  0.73333333333333  0.86666666666667  0.25000000000000  39.948 > 76
  100 Ar  0.93333333333333  0.86666666666667  0.25000000000000  39.948 > 76
  101 Ar  0.13333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  102 Ar  0.33333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  103 Ar  0.53333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  104 Ar  0.73333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  105 Ar  0.93333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  106 Ar  0.13333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  107 Ar  0.33333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  108 Ar  0.53333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  109 Ar  0.73333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  110 Ar  0.93333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  111 Ar  0.13333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  112 Ar  0.33333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  113 Ar  0.53333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  114 Ar  0.73333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  115 Ar  0.93333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  116 Ar  0.13333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  117 Ar  0.33333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  118 Ar  0.53333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  119 Ar  0.73333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  120 Ar  0.93333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  121 Ar  0.13333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  122 Ar  0.33333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  123 Ar  0.53333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  124 Ar  0.73333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  125 Ar  0.93333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  126 Ar  0.13333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  127 Ar  0.33333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  128 Ar  0.53333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  129 Ar  0.73333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  130 Ar  0.93333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  131 Ar  0.13333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  132 Ar  0.33333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  133 Ar  0.53333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  134 Ar  0.73333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  135 Ar  0.93333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  136 Ar  0.13333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  137 Ar  0.33333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  138 Ar  0.53333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  139 Ar  0.73333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  140 Ar  0.93333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  141 Ar  0.13333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  142 Ar  0.33333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  143 Ar  0.53333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  144 Ar  0.73333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  145 Ar  0.93333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  146 Ar  0.13333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  147 Ar  0.33333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  148 Ar  0.53333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  149 Ar  0.73333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  150 Ar  0.93333333333333  0.86666666666667  0.91666666666667  39.948 > 76
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            1.5972440    0.0000000    0.0000000
           -0.0000000    1.5972440    0.0000000
            0.0000000    0.0000000    1.5969012
-------------------------- Born effective charges --------------------------
    1 Ar    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    2 Ar    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
----------------------------------------------------------------------------
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.0000  0.0100]
    [ 0.0000  0.0000 -0.0100]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000000 (xxx) 0.00000000 (xxx) 0.00000000 (xxx)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 17:57:12]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 17:57:12]-------------------------
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): [5 5 3]
Primitive matrix:
  [1. 0. 0.]
  [0. 1. 0.]
  [0. 0. 1.]
Spacegroup: P6_3/mmc (194)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    3.977318140791256    0.000000000000000    0.000000000000000
  b   -1.988659070395628    3.444458548857920    0.000000000000000
  c    0.000000000000000    0.000000000000000    6.514740099999999
Atomic positions (fractional):
    1 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948
    2 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a   19.886590703956280    0.000000000000000    0.000000000000000
  b   -9.943295351978140   17.222292744289600    0.000000000000000
  c    0.000000000000000    0.000000000000000   19.544220299999999
Atomic positions (fractional):
    1 Ar  0.06666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    2 Ar  0.26666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    3 Ar  0.46666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    4 Ar  0.66666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    5 Ar  0.86666666666667  0.13333333333333  0.08333333333333  39.948 > 1
    6 Ar  0.06666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    7 Ar  0.26666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    8 Ar  0.46666666666667  0.33333333333333  0.08333333333333  39.948 > 1
    9 Ar  0.66666666666667  0.33333333333333  0.08333333333333  39.948 > 1
   10 Ar  0.86666666666667  0.33333333333333  0.08333333333333  39.948 > 1
   11 Ar  0.06666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   12 Ar  0.26666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   13 Ar  0.46666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   14 Ar  0.66666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   15 Ar  0.86666666666667  0.53333333333333  0.08333333333333  39.948 > 1
   16 Ar  0.06666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   17 Ar  0.26666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   18 Ar  0.46666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   19 Ar  0.66666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   20 Ar  0.86666666666667  0.73333333333333  0.08333333333333  39.948 > 1
   21 Ar  0.06666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   22 Ar  0.26666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   23 Ar  0.46666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   24 Ar  0.66666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   25 Ar  0.86666666666667  0.93333333333333  0.08333333333333  39.948 > 1
   26 Ar  0.06666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   27 Ar  0.26666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   28 Ar  0.46666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   29 Ar  0.66666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   30 Ar  0.86666666666667  0.13333333333333  0.41666666666667  39.948 > 1
   31 Ar  0.06666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   32 Ar  0.26666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   33 Ar  0.46666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   34 Ar  0.66666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   35 Ar  0.86666666666667  0.33333333333333  0.41666666666667  39.948 > 1
   36 Ar  0.06666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   37 Ar  0.26666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   38 Ar  0.46666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   39 Ar  0.66666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   40 Ar  0.86666666666667  0.53333333333333  0.41666666666667  39.948 > 1
   41 Ar  0.06666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   42 Ar  0.26666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   43 Ar  0.46666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   44 Ar  0.66666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   45 Ar  0.86666666666667  0.73333333333333  0.41666666666667  39.948 > 1
   46 Ar  0.06666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   47 Ar  0.26666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   48 Ar  0.46666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   49 Ar  0.66666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   50 Ar  0.86666666666667  0.93333333333333  0.41666666666667  39.948 > 1
   51 Ar  0.06666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   52 Ar  0.26666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   53 Ar  0.46666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   54 Ar  0.66666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   55 Ar  0.86666666666667  0.13333333333333  0.75000000000000  39.948 > 1
   56 Ar  0.06666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   57 Ar  0.26666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   58 Ar  0.46666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   59 Ar  0.66666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   60 Ar  0.86666666666667  0.33333333333333  0.75000000000000  39.948 > 1
   61 Ar  0.06666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   62 Ar  0.26666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   63 Ar  0.46666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   64 Ar  0.66666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   65 Ar  0.86666666666667  0.53333333333333  0.75000000000000  39.948 > 1
   66 Ar  0.06666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   67 Ar  0.26666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   68 Ar  0.46666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   69 Ar  0.66666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   70 Ar  0.86666666666667  0.73333333333333  0.75000000000000  39.948 > 1
   71 Ar  0.06666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   72 Ar  0.26666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   73 Ar  0.46666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   74 Ar  0.66666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   75 Ar  0.86666666666667  0.93333333333333  0.75000000000000  39.948 > 1
   76 Ar  0.13333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   77 Ar  0.33333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   78 Ar  0.53333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   79 Ar  0.73333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   80 Ar  0.93333333333333  0.06666666666667  0.25000000000000  39.948 > 76
   81 Ar  0.13333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   82 Ar  0.33333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   83 Ar  0.53333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   84 Ar  0.73333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   85 Ar  0.93333333333333  0.26666666666667  0.25000000000000  39.948 > 76
   86 Ar  0.13333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   87 Ar  0.33333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   88 Ar  0.53333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   89 Ar  0.73333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   90 Ar  0.93333333333333  0.46666666666667  0.25000000000000  39.948 > 76
   91 Ar  0.13333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   92 Ar  0.33333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   93 Ar  0.53333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   94 Ar  0.73333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   95 Ar  0.93333333333333  0.66666666666667  0.25000000000000  39.948 > 76
   96 Ar  0.13333333333333  0.86666666666667  0.25000000000000  39.948 > 76
   97 Ar  0.33333333333333  0.86666666666667  0.25000000000000  39.948 > 76
   98 Ar  0.53333333333333  0.86666666666667  0.25000000000000  39.948 > 76
   99 Ar  0.73333333333333  0.86666666666667  0.25000000000000  39.948 > 76
  100 Ar  0.93333333333333  0.86666666666667  0.25000000000000  39.948 > 76
  101 Ar  0.13333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  102 Ar  0.33333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  103 Ar  0.53333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  104 Ar  0.73333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  105 Ar  0.93333333333333  0.06666666666667  0.58333333333333  39.948 > 76
  106 Ar  0.13333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  107 Ar  0.33333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  108 Ar  0.53333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  109 Ar  0.73333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  110 Ar  0.93333333333333  0.26666666666667  0.58333333333333  39.948 > 76
  111 Ar  0.13333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  112 Ar  0.33333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  113 Ar  0.53333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  114 Ar  0.73333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  115 Ar  0.93333333333333  0.46666666666667  0.58333333333333  39.948 > 76
  116 Ar  0.13333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  117 Ar  0.33333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  118 Ar  0.53333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  119 Ar  0.73333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  120 Ar  0.93333333333333  0.66666666666667  0.58333333333333  39.948 > 76
  121 Ar  0.13333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  122 Ar  0.33333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  123 Ar  0.53333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  124 Ar  0.73333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  125 Ar  0.93333333333333  0.86666666666667  0.58333333333333  39.948 > 76
  126 Ar  0.13333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  127 Ar  0.33333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  128 Ar  0.53333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  129 Ar  0.73333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  130 Ar  0.93333333333333  0.06666666666667  0.91666666666667  39.948 > 76
  131 Ar  0.13333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  132 Ar  0.33333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  133 Ar  0.53333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  134 Ar  0.73333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  135 Ar  0.93333333333333  0.26666666666667  0.91666666666667  39.948 > 76
  136 Ar  0.13333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  137 Ar  0.33333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  138 Ar  0.53333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  139 Ar  0.73333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  140 Ar  0.93333333333333  0.46666666666667  0.91666666666667  39.948 > 76
  141 Ar  0.13333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  142 Ar  0.33333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  143 Ar  0.53333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  144 Ar  0.73333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  145 Ar  0.93333333333333  0.66666666666667  0.91666666666667  39.948 > 76
  146 Ar  0.13333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  147 Ar  0.33333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  148 Ar  0.53333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  149 Ar  0.73333333333333  0.86666666666667  0.91666666666667  39.948 > 76
  150 Ar  0.93333333333333  0.86666666666667  0.91666666666667  39.948 > 76
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            1.5972440    0.0000000    0.0000000
           -0.0000000    1.5972440    0.0000000
            0.0000000    0.0000000    1.5969012
-------------------------- Born effective charges --------------------------
    1 Ar    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
    2 Ar    0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
            0.0000000    0.0000000    0.0000000
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000000 (xxx) 0.00000000 (xxx) 0.00000000 (xxx)
Max drift of fc2: 0.00000000 (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 ... [ 15 15 8 ]
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: 313, Lambda: 0.12
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/135) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 135
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
   0.452   (   0.000    0.000    0.000)    0.000
   0.452   (   0.000    0.000    0.000)    0.000
   1.110   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/135) =======================
q-point: ( 0.07  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.081   (   3.727    2.152    0.000)    4.303
   0.144   (   6.109    3.527    0.000)    7.054
   0.195   (   8.632    4.984    0.000)    9.967
   0.470   (   1.567    0.905    0.000)    1.809
   0.496   (   3.727    2.152    0.000)    4.303
   1.105   (  -0.488   -0.282    0.000)    0.564
======================= Grid point 2 (3/135) =======================
q-point: ( 0.13  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.167   (   3.955    2.284    0.000)    4.567
   0.267   (   4.855    2.803    0.000)    5.606
   0.384   (   8.218    4.744    0.000)    9.489
   0.515   (   2.232    1.289    0.000)    2.577
   0.601   (   5.309    3.065    0.000)    6.130
   1.089   (  -0.914   -0.528    0.000)    1.056
======================= Grid point 3 (4/135) =======================
q-point: ( 0.20  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.255   (   3.857    2.227    0.000)    4.453
   0.361   (   3.635    2.099    0.000)    4.197
   0.560   (   7.490    4.324    0.000)    8.649
   0.565   (   2.235    1.290    0.000)    2.581
   0.722   (   5.295    3.057    0.000)    6.114
   1.065   (  -1.261   -0.728    0.000)    1.456
======================= Grid point 4 (5/135) =======================
q-point: ( 0.27  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.336   (   3.380    1.951    0.000)    3.903
   0.436   (   3.160    1.825    0.000)    3.649
   0.613   (   2.081    1.202    0.000)    2.403
   0.716   (   6.438    3.717    0.000)    7.434
   0.833   (   4.556    2.630    0.000)    5.261
   1.033   (  -1.562   -0.902    0.000)    1.804
======================= Grid point 5 (6/135) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.404   (   2.598    1.500    0.000)    3.000
   0.508   (   3.350    1.934    0.000)    3.868
   0.658   (   1.877    1.084    0.000)    2.167
   0.845   (   5.039    2.909    0.000)    5.818
   0.923   (   3.476    2.007    0.000)    4.014
   0.995   (  -1.836   -1.060    0.000)    2.119
======================= Grid point 6 (7/135) =======================
q-point: ( 0.40  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.451   (   1.632    0.942    0.000)    1.884
   0.583   (   3.218    1.858    0.000)    3.716
   0.695   (   1.394    0.805    0.000)    1.610
   0.939   (   3.277    1.892    0.000)    3.784
   0.952   (  -1.890   -1.091    0.000)    2.183
   0.986   (   2.160    1.247    0.000)    2.494
======================= Grid point 7 (8/135) =======================
q-point: ( 0.47  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.476   (   0.558    0.322    0.000)    0.645
   0.640   (   1.543    0.891    0.000)    1.782
   0.717   (   0.523    0.302    0.000)    0.604
   0.917   (  -0.997   -0.575    0.000)    1.151
   0.989   (   1.167    0.674    0.000)    1.347
   1.018   (   0.709    0.409    0.000)    0.819
======================= Grid point 16 (9/135) =======================
q-point: ( 0.07  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.155   (   2.654    4.597    0.000)    5.309
   0.237   (   3.011    5.215    0.000)    6.021
   0.325   (   4.382    7.590    0.000)    8.764
   0.507   (   1.582    2.739    0.000)    3.163
   0.567   (   2.808    4.863    0.000)    5.616
   1.094   (  -0.469   -0.812    0.000)    0.938
======================= Grid point 17 (10/135) =======================
q-point: ( 0.13  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.244   (   2.316    5.076    0.000)    5.579
   0.331   (   3.032    3.460    0.000)    4.601
   0.470   (   5.766    3.770    0.000)    6.890
   0.572   (   2.033    4.621    0.000)    5.048
   0.673   (   4.298    4.170    0.000)    5.989
   1.074   (  -0.868   -1.026    0.000)    1.345
======================= Grid point 18 (11/135) =======================
q-point: ( 0.20  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.328   (   1.869    4.972    0.000)    5.312
   0.407   (   2.721    2.527    0.000)    3.713
   0.575   (   3.954    1.176    0.000)    4.125
   0.673   (   3.865    5.706    0.000)    6.892
   0.784   (   4.461    3.224    0.000)    5.504
   1.046   (  -1.193   -1.232    0.000)    1.715
======================= Grid point 19 (12/135) =======================
q-point: ( 0.27  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.401   (   1.292    4.461    0.000)    4.644
   0.477   (   2.880    2.409    0.000)    3.755
   0.639   (   2.341    1.681    0.000)    2.882
   0.794   (   4.582    3.999    0.000)    6.082
   0.881   (   3.891    2.293    0.000)    4.517
   1.010   (  -1.470   -1.451    0.000)    2.065
======================= Grid point 20 (13/135) =======================
q-point: ( 0.33  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.456   (   0.538    3.788    0.000)    3.826
   0.551   (   3.096    2.565    0.000)    4.020
   0.689   (   1.349    2.157    0.000)    2.544
   0.897   (   3.880    2.269    0.000)    4.495
   0.956   (   2.895    1.394    0.000)    3.213
   0.969   (  -1.655   -1.639    0.000)    2.329
======================= Grid point 21 (14/135) =======================
q-point: ( 0.40  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.491   (  -0.373    3.095    0.000)    3.118
   0.623   (   2.302    2.212    0.000)    3.192
   0.724   (   0.269    2.123    0.000)    2.140
   0.927   (  -1.330   -1.561    0.000)    2.051
   0.964   (   2.416    0.697    0.000)    2.514
   1.004   (   1.597    0.519    0.000)    1.679
======================= Grid point 22 (15/135) =======================
q-point: ( 0.47  0.07  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.503   (  -1.409    2.440    0.000)    2.817
   0.656   (  -0.421    0.729    0.000)    0.842
   0.736   (  -0.934    1.618    0.000)    1.868
   0.905   (   0.383   -0.663    0.000)    0.766
   0.988   (   0.414   -0.717    0.000)    0.827
   1.019   (   0.178   -0.308    0.000)    0.355
======================= Grid point 33 (16/135) =======================
q-point: ( 0.13  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.345   (   2.905    5.031    0.000)    5.809
   0.396   (   1.874    3.245    0.000)    3.747
   0.538   (   2.035    3.524    0.000)    4.069
   0.675   (   3.189    5.524    0.000)    6.379
   0.758   (   2.596    4.497    0.000)    5.193
   1.050   (  -0.854   -1.478    0.000)    1.707
======================= Grid point 34 (17/135) =======================
q-point: ( 0.20  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.433   (   1.437    5.513    0.000)    5.697
   0.459   (   2.233    2.920    0.000)    3.676
   0.609   (   3.029    2.310    0.000)    3.809
   0.777   (   2.868    5.002    0.000)    5.766
   0.844   (   3.237    2.899    0.000)    4.346
   1.017   (  -1.182   -1.720    0.000)    2.087
======================= Grid point 35 (18/135) =======================
q-point: ( 0.27  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.499   (   0.145    5.426    0.000)    5.428
   0.530   (   2.717    3.111    0.000)    4.130
   0.680   (   2.586    2.559    0.000)    3.638
   0.865   (   2.732    3.250    0.000)    4.246
   0.919   (   3.112    1.604    0.000)    3.501
   0.977   (  -1.437   -1.971    0.000)    2.439
======================= Grid point 36 (19/135) =======================
q-point: ( 0.33  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.543   (  -0.880    4.977    0.000)    5.054
   0.607   (   2.640    3.116    0.000)    4.084
   0.740   (   1.148    3.048    0.000)    3.257
   0.929   (   2.401    1.051    0.000)    2.621
   0.932   (  -1.468   -2.123    0.000)    2.581
   0.976   (   2.294    0.663    0.000)    2.388
======================= Grid point 37 (20/135) =======================
q-point: ( 0.40  0.13  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.565   (  -1.869    4.408    0.000)    4.787
   0.667   (   0.916    2.267    0.000)    2.445
   0.774   (  -0.640    2.941    0.000)    3.009
   0.894   (  -0.490   -1.764    0.000)    1.831
   0.964   (   1.502   -0.786    0.000)    1.695
   1.007   (   1.031   -0.154    0.000)    1.042
======================= Grid point 49 (21/135) =======================
q-point: ( 0.20  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.522   (   2.070    3.585    0.000)    4.139
   0.532   (   2.615    4.529    0.000)    5.229
   0.664   (   1.932    3.347    0.000)    3.865
   0.864   (   2.229    3.861    0.000)    4.459
   0.895   (   1.390    2.407    0.000)    2.780
   0.979   (  -1.237   -2.143    0.000)    2.475
======================= Grid point 50 (22/135) =======================
q-point: ( 0.27  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.596   (   2.481    3.704    0.000)    4.459
   0.603   (   0.275    5.159    0.000)    5.167
   0.735   (   2.685    2.980    0.000)    4.011
   0.918   (   0.622    2.123    0.000)    2.212
   0.935   (  -1.448   -2.363    0.000)    2.771
   0.942   (   2.462    0.859    0.000)    2.608
======================= Grid point 51 (23/135) =======================
q-point: ( 0.33  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.643   (  -1.460    5.183    0.000)    5.385
   0.670   (   1.902    3.319    0.000)    3.825
   0.801   (   1.383    3.096    0.000)    3.391
   0.888   (  -1.211   -2.399    0.000)    2.687
   0.935   (   0.657   -0.525    0.000)    0.841
   0.983   (   1.663    0.131    0.000)    1.668
======================= Grid point 52 (24/135) =======================
q-point: ( 0.40  0.20  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.655   (  -2.724    4.719    0.000)    5.449
   0.707   (  -1.038    1.797    0.000)    2.075
   0.830   (  -1.485    2.572    0.000)    2.970
   0.861   (   0.855   -1.481    0.000)    1.710
   0.936   (   1.090   -1.888    0.000)    2.180
   0.998   (   0.415   -0.719    0.000)    0.830
======================= Grid point 66 (25/135) =======================
q-point: ( 0.27  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 320
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.669   (   2.148    3.720    0.000)    4.295
   0.690   (   2.086    3.613    0.000)    4.171
   0.794   (   1.811    3.136    0.000)    3.622
   0.887   (  -1.479   -2.561    0.000)    2.957
   0.934   (  -0.276   -0.479    0.000)    0.553
   0.964   (   0.758    1.313    0.000)    1.517
======================= Grid point 67 (26/135) =======================
q-point: ( 0.33  0.27  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 565
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.732   (   0.955    3.109    0.000)    3.252
   0.737   (  -1.122    4.338    0.000)    4.480
   0.840   (  -0.750   -2.540    0.000)    2.648
   0.853   (   1.927    2.126    0.000)    2.870
   0.910   (  -0.927   -1.868    0.000)    2.085
   0.982   (   0.690   -0.185    0.000)    0.715
======================= Grid point 82 (27/135) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.789   (  -1.377   -2.385    0.000)    2.754
   0.789   (   1.377    2.385    0.000)    2.754
   0.790   (  -0.000   -0.000    0.000)    0.000
   0.879   (  -0.405   -0.702    0.000)    0.810
   0.879   (   0.405    0.702    0.000)    0.810
   0.980   (  -0.000   -0.000    0.000)    0.000
======================= Grid point 241 (28/135) =======================
q-point: ( 0.00  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.095   (   0.000   -0.000    4.812)    4.812
   0.095   (   0.000    0.000    4.812)    4.812
   0.273   (   0.000    0.000   13.679)   13.679
   0.442   (   0.000    0.000   -0.952)    0.952
   0.442   (  -0.000    0.000   -0.952)    0.952
   1.092   (   0.000   -0.000   -1.838)    1.838
======================= Grid point 242 (29/135) =======================
q-point: ( 0.07  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.127   (   2.576    1.487    3.845)    4.861
   0.192   (   5.961    3.441    2.739)    7.408
   0.316   (   3.762    2.172   10.755)   11.599
   0.461   (   1.546    0.893   -0.981)    2.037
   0.487   (   3.692    2.132   -0.990)    4.377
   1.088   (  -0.387   -0.224   -1.726)    1.783
======================= Grid point 243 (30/135) =======================
q-point: ( 0.13  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.197   (   3.498    2.020    2.817)    4.925
   0.308   (   4.234    2.444    3.147)    5.814
   0.436   (   6.580    3.799    5.424)    9.335
   0.505   (   2.241    1.294   -0.997)    2.773
   0.591   (   5.267    3.041   -1.069)    6.175
   1.075   (  -0.770   -0.445   -1.452)    1.703
======================= Grid point 244 (31/135) =======================
q-point: ( 0.20  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.277   (   3.570    2.061    2.162)    4.655
   0.389   (   3.204    1.850    2.630)    4.539
   0.556   (   2.293    1.324   -0.915)    2.802
   0.590   (   6.908    3.989    3.013)    8.527
   0.710   (   5.238    3.024   -1.171)    6.160
   1.054   (  -1.106   -0.638   -1.120)    1.699
======================= Grid point 245 (32/135) =======================
q-point: ( 0.27  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.353   (   3.149    1.818    1.649)    3.993
   0.457   (   2.944    1.700    2.077)    3.984
   0.606   (   2.155    1.244   -0.761)    2.603
   0.736   (   6.117    3.532    2.039)    7.352
   0.818   (   4.321    2.494   -1.447)    5.195
   1.027   (  -1.227   -0.708   -0.617)    1.545
======================= Grid point 246 (33/135) =======================
q-point: ( 0.33  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.416   (   2.422    1.398    1.233)    3.057
   0.523   (   2.994    1.728    1.507)    3.771
   0.652   (   1.918    1.107   -0.638)    2.304
   0.859   (   4.817    2.781    1.451)    5.748
   0.895   (   2.212    1.277   -2.541)    3.603
   1.006   (  -0.392   -0.227    0.777)    0.899
======================= Grid point 247 (34/135) =======================
q-point: ( 0.40  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.460   (   1.525    0.880    0.943)    1.998
   0.588   (   2.657    1.534    0.505)    3.109
   0.690   (   1.395    0.805   -0.606)    1.721
   0.908   (  -0.630   -0.364   -3.176)    3.258
   0.949   (   3.144    1.815    1.080)    3.787
   1.019   (   1.208    0.697    1.920)    2.373
======================= Grid point 248 (35/135) =======================
q-point: ( 0.47  0.00  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.483   (   0.523    0.302    0.797)    1.000
   0.633   (   1.187    0.685   -0.574)    1.486
   0.711   (   0.517    0.299   -0.625)    0.864
   0.890   (  -0.612   -0.353   -2.396)    2.498
   0.998   (   1.125    0.650    0.890)    1.575
   1.041   (   0.557    0.322    1.783)    1.896
======================= Grid point 257 (36/135) =======================
q-point: ( 0.07  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.186   (   2.308    3.997    2.915)    5.459
   0.277   (   2.624    4.544    2.770)    5.933
   0.392   (   3.071    5.319    6.758)    9.132
   0.499   (   1.673    2.897   -0.808)    3.442
   0.557   (   2.778    4.812   -1.071)    5.659
   1.079   (  -0.387   -0.671   -1.533)    1.718
======================= Grid point 258 (37/135) =======================
q-point: ( 0.13  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.267   (   2.125    4.668    2.191)    5.578
   0.358   (   2.897    2.626    2.470)    4.624
   0.500   (   4.335    2.766    2.590)    5.758
   0.572   (   2.627    5.338    0.500)    5.971
   0.662   (   4.281    4.114   -1.150)    6.047
   1.062   (  -0.757   -0.885   -1.235)    1.697
======================= Grid point 259 (38/135) =======================
q-point: ( 0.20  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.345   (   1.759    4.561    1.626)    5.152
   0.427   (   2.766    2.052    1.982)    3.974
   0.579   (   2.877    1.665    0.183)    3.329
   0.688   (   4.409    5.340    1.719)    7.135
   0.771   (   4.372    3.120   -1.290)    5.524
   1.038   (  -1.021   -1.039   -0.848)    1.686
======================= Grid point 260 (39/135) =======================
q-point: ( 0.27  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.412   (   1.279    4.036    1.158)    4.389
   0.493   (   2.791    2.129    1.643)    3.876
   0.637   (   1.930    2.094   -0.274)    2.861
   0.809   (   4.588    3.705    1.539)    6.094
   0.862   (   3.362    1.906   -1.867)    4.292
   1.011   (  -0.834   -0.923    0.002)    1.244
======================= Grid point 261 (40/135) =======================
q-point: ( 0.33  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.464   (   0.593    3.462    0.838)    3.611
   0.561   (   2.686    2.247    0.984)    3.637
   0.685   (   1.147    2.289   -0.471)    2.603
   0.904   (   1.902    0.525   -1.643)    2.567
   0.912   (   2.482    1.555   -0.303)    2.945
   1.004   (   0.839    0.052    1.711)    1.907
======================= Grid point 262 (41/135) =======================
q-point: ( 0.40  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.498   (  -0.309    2.894    0.677)    2.988
   0.621   (   1.788    1.934   -0.069)    2.635
   0.717   (   0.201    1.999   -0.694)    2.125
   0.895   (  -0.666   -0.959   -2.638)    2.885
   0.975   (   2.301    0.714    1.012)    2.613
   1.028   (   1.239    0.222    1.835)    2.226
======================= Grid point 263 (42/135) =======================
q-point: ( 0.47  0.07  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.509   (  -1.327    2.299    0.630)    2.729
   0.648   (  -0.461    0.798   -0.694)    1.154
   0.729   (  -0.807    1.398   -0.829)    1.814
   0.882   (   0.247   -0.428   -2.055)    2.113
   0.998   (   0.374   -0.648    0.923)    1.188
   1.040   (   0.216   -0.374    1.729)    1.782
======================= Grid point 274 (43/135) =======================
q-point: ( 0.13  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.361   (   2.731    4.730    1.579)    5.685
   0.409   (   1.597    2.766    1.468)    3.515
   0.557   (   1.860    3.222    1.547)    4.029
   0.683   (   3.295    5.708    1.080)    6.678
   0.745   (   2.533    4.387   -1.299)    5.230
   1.041   (  -0.733   -1.269   -0.900)    1.720
======================= Grid point 275 (44/135) =======================
q-point: ( 0.20  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.441   (   1.027    5.035    0.874)    5.212
   0.470   (   2.715    2.517    1.166)    3.881
   0.623   (   2.391    2.774    1.100)    3.824
   0.787   (   2.986    4.823    1.112)    5.780
   0.827   (   3.059    2.566   -1.627)    4.312
   1.015   (  -0.862   -1.289   -0.289)    1.577
======================= Grid point 276 (45/135) =======================
q-point: ( 0.27  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.498   (   0.446    4.530    0.185)    4.556
   0.543   (   2.511    3.116    1.155)    4.166
   0.687   (   1.821    3.031    0.474)    3.567
   0.871   (   2.573    2.177   -0.008)    3.371
   0.891   (   2.214    1.412   -1.573)    3.061
   0.994   (  -0.064   -0.754    1.132)    1.361
======================= Grid point 277 (46/135) =======================
q-point: ( 0.33  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.542   (  -0.424    4.333   -0.041)    4.353
   0.612   (   2.074    3.022    0.567)    3.709
   0.738   (   0.524    3.059   -0.394)    3.129
   0.891   (  -0.020   -1.189   -2.763)    3.008
   0.943   (   2.149    1.138    1.079)    2.660
   1.005   (   1.487   -0.040    1.895)    2.409
======================= Grid point 278 (47/135) =======================
q-point: ( 0.40  0.13  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.566   (  -1.541    3.919   -0.021)    4.212
   0.663   (   0.516    2.385   -0.184)    2.447
   0.761   (  -0.769    2.214   -1.380)    2.719
   0.876   (  -0.225   -0.765   -1.543)    1.737
   0.975   (   1.371   -0.649    1.121)    1.887
   1.026   (   0.984   -0.418    1.635)    1.953
======================= Grid point 290 (48/135) =======================
q-point: ( 0.20  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.518   (   1.774    3.072   -0.067)    3.548
   0.539   (   2.527    4.376    0.720)    5.104
   0.684   (   2.001    3.465    1.624)    4.318
   0.866   (   0.807    1.398   -2.585)    3.048
   0.872   (   2.189    3.792    0.911)    4.472
   0.991   (  -0.582   -1.009    0.829)    1.430
======================= Grid point 291 (49/135) =======================
q-point: ( 0.27  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.577   (   1.260    3.653   -1.071)    4.009
   0.616   (   1.662    4.220    0.738)    4.595
   0.751   (   1.721    3.500    1.224)    4.088
   0.875   (   0.227   -1.084   -3.130)    3.320
   0.932   (   1.259    2.205    0.745)    2.646
   0.983   (   0.635   -0.515    1.951)    2.115
======================= Grid point 292 (50/135) =======================
q-point: ( 0.33  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.623   (  -0.306    3.923   -1.567)    4.235
   0.678   (   1.092    3.773    0.663)    3.984
   0.793   (  -0.588    2.162   -1.221)    2.552
   0.868   (   0.730   -0.742   -1.055)    1.483
   0.953   (   0.265   -0.142    1.291)    1.326
   1.000   (   1.543   -0.519    1.490)    2.207
======================= Grid point 293 (51/135) =======================
q-point: ( 0.40  0.20  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.640   (  -2.014    3.488   -1.517)    4.304
   0.713   (  -1.533    2.655    0.685)    3.141
   0.788   (  -0.212    0.367   -3.115)    3.143
   0.877   (  -0.589    1.021    0.462)    1.266
   0.952   (   0.918   -1.589    1.477)    2.355
   1.012   (   0.631   -1.093    1.260)    1.783
======================= Grid point 307 (52/135) =======================
q-point: ( 0.27  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.643   (   1.820    3.153   -2.060)    4.183
   0.692   (   1.976    3.422    0.117)    3.953
   0.819   (   1.915    3.318    1.889)    4.271
   0.843   (  -1.078   -1.866   -3.000)    3.694
   0.967   (   0.649    1.124    0.433)    1.368
   0.969   (  -0.463   -0.801    2.214)    2.399
======================= Grid point 308 (53/135) =======================
q-point: ( 0.33  0.27  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.692   (  -0.078    3.146   -2.766)    4.189
   0.748   (   1.221    2.882   -0.285)    3.143
   0.797   (  -2.119   -1.170   -2.208)    3.276
   0.885   (   2.267    1.993    1.814)    3.522
   0.939   (  -1.637   -1.188    2.001)    2.845
   0.986   (   1.062   -0.892    0.485)    1.469
======================= Grid point 323 (54/135) =======================
q-point: ( 0.33  0.33  0.12)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.745   (  -1.160   -2.008   -3.398)    4.114
   0.745   (   1.160    2.008   -3.398)    4.114
   0.797   (  -0.000   -0.000    0.718)    0.718
   0.916   (  -0.524   -0.907    2.785)    2.975
   0.916   (   0.524    0.907    2.785)    2.975
   0.973   (  -0.000   -0.000   -0.659)    0.659
======================= Grid point 482 (55/135) =======================
q-point: ( 0.00  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 135
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.183   (   0.000    0.000    4.304)    4.304
   0.183   (   0.000   -0.000    4.304)    4.304
   0.416   (   0.000    0.000   -1.771)    1.771
   0.416   (   0.000    0.000   -1.771)    1.771
   0.516   (   0.000    0.000   11.526)   11.526
   1.042   (  -0.000    0.000   -3.370)    3.370
======================= Grid point 483 (56/135) =======================
q-point: ( 0.07  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.205   (   1.902    1.098    4.050)    4.607
   0.251   (   5.039    2.909    3.245)    6.662
   0.434   (   1.511    0.872   -1.796)    2.504
   0.459   (   3.574    2.063   -1.768)    4.489
   0.528   (   1.141    0.659   10.606)   10.687
   1.040   (  -0.185   -0.107   -3.259)    3.266
======================= Grid point 484 (57/135) =======================
q-point: ( 0.13  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.261   (   2.914    1.682    3.546)    4.888
   0.366   (   4.705    2.717    2.933)    6.174
   0.478   (   2.286    1.320   -1.780)    3.184
   0.557   (   4.647    2.683    0.216)    5.370
   0.578   (   3.672    2.120    5.973)    7.325
   1.033   (  -0.436   -0.252   -2.938)    2.981
======================= Grid point 485 (58/135) =======================
q-point: ( 0.20  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.329   (   3.085    1.781    3.034)    4.679
   0.450   (   2.960    1.709    3.493)    4.887
   0.531   (   2.441    1.409   -1.638)    3.260
   0.665   (   5.307    3.064    3.214)    6.919
   0.682   (   4.967    2.868   -0.252)    5.741
   1.021   (  -0.636   -0.367   -2.376)    2.487
======================= Grid point 486 (59/135) =======================
q-point: ( 0.27  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.395   (   2.733    1.578    2.573)    4.072
   0.509   (   2.474    1.428    3.269)    4.341
   0.585   (   2.307    1.332   -1.444)    3.031
   0.778   (   3.934    2.271   -2.076)    4.995
   0.790   (   5.059    2.921    2.790)    6.473
   1.009   (  -0.292   -0.169   -1.351)    1.392
======================= Grid point 487 (60/135) =======================
q-point: ( 0.33  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.449   (   2.104    1.215    2.198)    3.276
   0.562   (   2.204    1.273    2.505)    3.571
   0.633   (   1.972    1.139   -1.343)    2.644
   0.835   (   1.203    0.694   -3.565)    3.826
   0.898   (   4.297    2.481    2.378)    5.502
   1.016   (   1.050    0.606    0.197)    1.229
======================= Grid point 488 (61/135) =======================
q-point: ( 0.40  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.488   (   1.327    0.766    1.934)    2.467
   0.604   (   1.500    0.866    1.322)    2.179
   0.671   (   1.364    0.788   -1.372)    2.089
   0.842   (  -0.226   -0.130   -3.688)    3.698
   0.978   (   2.815    1.625    1.805)    3.718
   1.048   (   1.475    0.851    1.047)    1.999
======================= Grid point 489 (62/135) =======================
q-point: ( 0.47  0.00  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.508   (   0.455    0.263    1.794)    1.870
   0.627   (   0.507    0.292    0.195)    0.617
   0.692   (   0.491    0.284   -1.438)    1.546
   0.836   (  -0.184   -0.106   -3.180)    3.187
   1.021   (   0.989    0.571    1.462)    1.855
   1.072   (   0.594    0.343    1.298)    1.468
======================= Grid point 498 (63/135) =======================
q-point: ( 0.07  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.251   (   1.886    3.267    3.594)    5.211
   0.329   (   2.640    4.573    2.742)    5.950
   0.468   (   1.390    2.408   -1.396)    3.111
   0.528   (   2.720    4.711   -1.873)    5.754
   0.558   (   1.515    2.624    8.664)    9.179
   1.036   (  -0.210   -0.364   -3.045)    3.074
======================= Grid point 499 (64/135) =======================
q-point: ( 0.13  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.319   (   1.831    3.958    3.025)    5.307
   0.412   (   3.188    1.994    2.852)    4.719
   0.519   (   1.671    2.690   -0.674)    3.237
   0.626   (   3.537    4.385    1.010)    5.724
   0.638   (   4.005    3.890    2.570)    6.147
   1.026   (  -0.449   -0.520   -2.612)    2.701
======================= Grid point 500 (65/135) =======================
q-point: ( 0.20  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.386   (   1.575    3.840    2.534)    4.863
   0.476   (   2.873    1.169    2.891)    4.241
   0.574   (   1.620    2.841   -0.676)    3.339
   0.732   (   4.219    3.254   -0.675)    5.371
   0.745   (   4.243    3.834    2.087)    6.087
   1.012   (  -0.467   -0.493   -1.846)    1.967
======================= Grid point 501 (66/135) =======================
q-point: ( 0.27  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.444   (   1.197    3.315    2.119)    4.112
   0.534   (   2.624    1.339    2.465)    3.841
   0.626   (   1.275    2.825   -0.838)    3.210
   0.813   (   2.527    1.259   -2.970)    4.098
   0.852   (   4.184    3.253    2.623)    5.913
   1.007   (   0.381    0.037   -0.491)    0.623
======================= Grid point 502 (67/135) =======================
q-point: ( 0.33  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.489   (   0.649    2.745    1.748)    3.318
   0.588   (   1.965    1.543    1.725)    3.036
   0.669   (   0.729    2.438   -1.149)    2.792
   0.840   (   0.379   -0.046   -3.618)    3.639
   0.942   (   3.372    1.999    2.073)    4.434
   1.028   (   1.572    0.570    0.781)    1.845
======================= Grid point 503 (68/135) =======================
q-point: ( 0.40  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.518   (  -0.100    2.212    1.456)    2.650
   0.626   (   0.760    1.545    0.771)    1.887
   0.696   (   0.105    1.615   -1.555)    2.244
   0.838   (  -0.225   -0.170   -3.174)    3.186
   1.001   (   2.024    0.681    1.648)    2.698
   1.059   (   1.310    0.254    1.220)    1.808
======================= Grid point 504 (69/135) =======================
q-point: ( 0.47  0.07  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.528   (  -0.973    1.685    1.338)    2.361
   0.641   (  -0.654    1.133    0.232)    1.328
   0.704   (  -0.510    0.883   -1.723)    2.002
   0.835   (  -0.045    0.078   -2.784)    2.785
   1.021   (   0.300   -0.520    1.473)    1.591
   1.072   (   0.234   -0.405    1.329)    1.409
======================= Grid point 515 (70/135) =======================
q-point: ( 0.13  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.400   (   2.365    4.096    2.326)    5.271
   0.450   (   1.234    2.137    2.536)    3.538
   0.576   (   1.808    3.131    0.256)    3.624
   0.711   (   2.382    4.126   -2.163)    5.232
   0.727   (   2.988    5.176    3.532)    6.942
   1.014   (  -0.381   -0.659   -1.997)    2.138
======================= Grid point 516 (71/135) =======================
q-point: ( 0.20  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.467   (   0.684    4.266    1.814)    4.686
   0.502   (   2.887    1.635    2.015)    3.882
   0.637   (   1.463    3.572    0.254)    3.869
   0.783   (   2.439    2.008   -2.578)    4.078
   0.826   (   3.217    4.186    2.581)    5.876
   1.004   (  -0.102   -0.319   -0.920)    0.979
======================= Grid point 517 (72/135) =======================
q-point: ( 0.27  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.512   (   0.426    3.324    1.252)    3.577
   0.569   (   2.491    2.460    1.488)    3.804
   0.690   (   0.535    3.656   -0.272)    3.705
   0.825   (   1.396    0.074   -3.169)    3.464
   0.911   (   2.851    2.754    2.144)    4.506
   1.010   (   0.973    0.194    0.433)    1.083
======================= Grid point 518 (73/135) =======================
q-point: ( 0.33  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.546   (  -0.003    2.858    0.562)    2.913
   0.628   (   1.501    2.744    1.063)    3.304
   0.721   (  -0.595    2.680   -1.298)    3.037
   0.838   (   0.347   -0.042   -2.763)    2.785
   0.972   (   2.157    1.049    1.859)    3.035
   1.035   (   1.613    0.060    1.125)    1.968
======================= Grid point 519 (74/135) =======================
q-point: ( 0.40  0.13  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.565   (  -0.810    2.434    0.082)    2.567
   0.669   (  -0.136    2.890    0.881)    3.024
   0.724   (  -0.860    0.944   -2.322)    2.650
   0.842   (  -0.286    0.802   -1.984)    2.159
   1.003   (   1.187   -0.454    1.676)    2.103
   1.056   (   1.029   -0.530    1.294)    1.736
======================= Grid point 531 (75/135) =======================
q-point: ( 0.20  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.528   (   1.183    2.048    1.132)    2.622
   0.557   (   2.269    3.930    1.074)    4.663
   0.708   (   2.067    3.581    0.622)    4.181
   0.809   (   0.398    0.689   -3.134)    3.234
   0.901   (   2.012    3.485    2.053)    4.517
   1.002   (   0.119    0.207    0.239)    0.338
======================= Grid point 532 (76/135) =======================
q-point: ( 0.27  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.566   (   0.751    2.383   -0.009)    2.499
   0.629   (   1.939    3.401    0.517)    3.949
   0.763   (  -0.391    3.744   -0.220)    3.771
   0.819   (   1.447   -0.523   -2.384)    2.838
   0.957   (   1.367    1.915    1.700)    2.903
   1.012   (   0.910   -0.095    1.053)    1.395
======================= Grid point 533 (77/135) =======================
q-point: ( 0.33  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.596   (  -0.015    2.394   -1.038)    2.609
   0.686   (   1.362    2.707   -0.273)    3.042
   0.761   (  -2.144    1.342   -1.338)    2.861
   0.853   (   0.992    1.920   -0.883)    2.335
   0.982   (   0.657   -0.019    1.739)    1.860
   1.029   (   1.435   -0.772    1.267)    2.064
======================= Grid point 534 (78/135) =======================
q-point: ( 0.40  0.20  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.608   (  -1.157    2.005   -1.461)    2.737
   0.725   (   0.255   -0.441   -3.346)    3.384
   0.734   (  -1.906    3.302    1.530)    4.108
   0.872   (  -1.216    2.107   -0.904)    2.595
   0.986   (   0.669   -1.158    1.927)    2.346
   1.038   (   0.802   -1.390    1.240)    2.028
======================= Grid point 548 (79/135) =======================
q-point: ( 0.27  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 488
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.609   (   1.174    2.033   -1.340)    2.703
   0.686   (   1.238    2.145   -1.095)    2.708
   0.797   (  -0.225   -0.390   -1.357)    1.430
   0.845   (   1.953    3.382    0.504)    3.938
   0.984   (   0.413    0.715    1.239)    1.489
   1.005   (  -0.294   -0.510    1.474)    1.587
======================= Grid point 549 (80/135) =======================
q-point: ( 0.33  0.27  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 901
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.640   (  -0.058    2.194   -2.324)    3.196
   0.707   (   0.617   -0.625   -3.171)    3.290
   0.793   (  -1.172    2.051    1.110)    2.610
   0.904   (   0.964    3.019   -0.031)    3.169
   0.974   (  -1.023   -0.675    1.784)    2.165
   1.004   (   1.254   -1.735    1.244)    2.476
======================= Grid point 564 (81/135) =======================
q-point: ( 0.33  0.33  0.25)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 192
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.680   (  -0.927   -1.605   -3.168)    3.671
   0.680   (   0.927    1.605   -3.168)    3.671
   0.817   (  -0.000   -0.000    1.390)    1.390
   0.954   (  -0.000   -0.000   -1.293)    1.293
   0.966   (  -0.434   -0.752    2.258)    2.419
   0.966   (   0.434    0.752    2.258)    2.419
======================= Grid point 723 (82/135) =======================
q-point: ( 0.00  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 108
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.259   (   0.000    0.000    3.658)    3.658
   0.259   (   0.000   -0.000    3.658)    3.658
   0.375   (  -0.000    0.000   -2.431)    2.431
   0.375   (   0.000    0.000   -2.431)    2.431
   0.712   (   0.000    0.000    8.844)    8.844
   0.964   (  -0.000    0.000   -4.739)    4.739
======================= Grid point 724 (83/135) =======================
q-point: ( 0.07  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.278   (   1.651    0.953    3.558)    4.037
   0.312   (   4.247    2.452    3.130)    5.818
   0.393   (   1.505    0.869   -2.427)    2.985
   0.419   (   3.601    2.079   -2.368)    4.785
   0.713   (   0.198    0.114    8.587)    8.590
   0.964   (  -0.061   -0.035   -4.719)    4.720
======================= Grid point 725 (84/135) =======================
q-point: ( 0.13  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.327   (   2.590    1.495    3.332)    4.477
   0.421   (   4.949    2.858    2.769)    6.351
   0.438   (   2.364    1.365   -2.373)    3.617
   0.520   (   5.058    2.920   -2.322)    6.285
   0.727   (   1.200    0.693    7.594)    7.720
   0.962   (  -0.079   -0.046   -4.513)    4.514
======================= Grid point 726 (85/135) =======================
q-point: ( 0.20  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.389   (   2.782    1.606    3.093)    4.459
   0.494   (   2.587    1.493   -2.250)    3.740
   0.517   (   3.458    1.996    3.375)    5.228
   0.632   (   4.776    2.758   -2.590)    6.093
   0.774   (   3.063    1.768    5.567)    6.595
   0.962   (   0.160    0.092   -3.771)    3.776
======================= Grid point 727 (86/135) =======================
q-point: ( 0.27  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.448   (   2.485    1.435    2.907)    4.085
   0.550   (   2.408    1.390   -2.147)    3.512
   0.578   (   2.156    1.245    3.792)    4.537
   0.721   (   2.958    1.708   -3.366)    4.795
   0.857   (   4.081    2.356    3.578)    5.916
   0.974   (   1.077    0.622   -2.256)    2.576
======================= Grid point 728 (87/135) =======================
q-point: ( 0.33  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.498   (   1.921    1.109    2.781)    3.557
   0.600   (   1.957    1.130   -2.146)    3.116
   0.618   (   1.422    0.821    3.263)    3.653
   0.762   (   0.810    0.468   -3.957)    4.066
   0.945   (   3.649    2.107    2.441)    4.869
   1.011   (   2.078    1.200   -0.759)    2.516
======================= Grid point 729 (88/135) =======================
q-point: ( 0.40  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.533   (   1.210    0.699    2.698)    3.038
   0.636   (   1.284    0.741   -2.235)    2.682
   0.641   (   0.621    0.359    2.473)    2.574
   0.769   (   0.067    0.039   -3.832)    3.833
   1.014   (   2.418    1.396    1.842)    3.345
   1.057   (   1.820    1.051   -0.106)    2.104
======================= Grid point 730 (89/135) =======================
q-point: ( 0.47  0.00  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.551   (   0.413    0.239    2.653)    2.696
   0.648   (   0.093    0.054    2.003)    2.006
   0.656   (   0.449    0.259   -2.320)    2.378
   0.770   (   0.030    0.017   -3.629)    3.630
   1.051   (   0.835    0.482    1.557)    1.832
   1.086   (   0.691    0.399    0.099)    0.804
======================= Grid point 739 (90/135) =======================
q-point: ( 0.07  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.318   (   1.614    2.795    3.256)    4.585
   0.379   (   2.247    3.891    2.215)    5.009
   0.434   (   1.884    3.263   -1.737)    4.149
   0.487   (   2.705    4.685   -2.308)    5.881
   0.721   (   0.473    0.818    7.961)    8.017
   0.962   (  -0.051   -0.088   -4.605)    4.606
======================= Grid point 740 (91/135) =======================
q-point: ( 0.13  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.377   (   1.733    3.314    2.892)    4.728
   0.451   (   3.077    1.635    0.471)    3.517
   0.508   (   2.145    4.105    0.465)    4.655
   0.588   (   3.992    3.957   -2.356)    6.094
   0.752   (   1.794    2.044    6.369)    6.925
   0.961   (   0.012    0.000   -4.126)    4.126
======================= Grid point 741 (92/135) =======================
q-point: ( 0.20  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.437   (   1.605    3.203    2.686)    4.478
   0.514   (   2.876    1.072   -0.021)    3.069
   0.574   (   1.435    3.309    1.706)    3.990
   0.686   (   3.543    2.601   -2.848)    5.237
   0.820   (   3.232    2.929    4.255)    6.094
   0.966   (   0.505    0.384   -2.959)    3.026
======================= Grid point 742 (93/135) =======================
q-point: ( 0.27  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.490   (   1.245    2.724    2.568)    3.945
   0.569   (   2.495    0.877   -0.051)    2.645
   0.622   (   0.724    2.919    1.531)    3.375
   0.748   (   1.711    1.072   -3.598)    4.125
   0.907   (   3.557    2.651    2.831)    5.263
   0.990   (   1.602    0.979   -1.319)    2.294
======================= Grid point 743 (94/135) =======================
q-point: ( 0.33  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.530   (   0.753    2.030    2.439)    3.262
   0.614   (   1.851    0.846   -0.244)    2.050
   0.654   (  -0.052    2.424    0.829)    2.563
   0.770   (   0.261    0.461   -3.693)    3.731
   0.983   (   2.917    1.713    2.074)    3.968
   1.033   (   2.111    0.989   -0.297)    2.350
======================= Grid point 744 (95/135) =======================
q-point: ( 0.40  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.554   (   0.189    1.292    2.284)    2.631
   0.646   (   0.913    0.940   -0.376)    1.363
   0.666   (  -0.733    1.624    0.169)    1.790
   0.775   (  -0.186    0.668   -3.377)    3.447
   1.034   (   1.725    0.588    1.681)    2.479
   1.071   (   1.473    0.391    0.049)    1.525
======================= Grid point 745 (96/135) =======================
q-point: ( 0.47  0.07  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.562   (  -0.430    0.744    2.206)    2.367
   0.663   (  -0.880    1.523    2.047)    2.699
   0.663   (  -0.229    0.396   -2.524)    2.565
   0.777   (  -0.427    0.739   -3.242)    3.352
   1.051   (   0.248   -0.429    1.554)    1.631
   1.086   (   0.230   -0.399    0.125)    0.477
======================= Grid point 756 (97/135) =======================
q-point: ( 0.13  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.445   (   2.017    3.494    2.253)    4.621
   0.490   (   1.290    2.234    1.006)    2.769
   0.581   (   1.985    3.438    0.864)    4.062
   0.665   (   2.197    3.805   -2.542)    5.076
   0.807   (   2.049    3.548    4.520)    6.101
   0.964   (   0.205    0.356   -3.205)    3.231
======================= Grid point 757 (98/135) =======================
q-point: ( 0.20  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.505   (   0.860    3.606    1.959)    4.193
   0.537   (   2.398    1.267    1.198)    2.965
   0.641   (   1.101    3.616    0.652)    3.836
   0.728   (   1.988    1.702   -2.976)    3.963
   0.882   (   2.695    3.329    3.040)    5.253
   0.978   (   0.869    0.879   -1.748)    2.141
======================= Grid point 758 (99/135) =======================
q-point: ( 0.27  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.543   (   0.375    2.424    1.903)    3.104
   0.591   (   2.213    1.618    0.455)    2.779
   0.685   (  -0.082    3.558    0.274)    3.570
   0.765   (   1.120    0.714   -3.064)    3.340
   0.955   (   2.588    2.245    2.205)    4.074
   1.009   (   1.688    0.931   -0.488)    1.989
======================= Grid point 759 (100/135) =======================
q-point: ( 0.33  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.566   (   0.096    1.521    1.570)    2.188
   0.637   (   1.418    1.450   -0.783)    2.174
   0.705   (  -1.158    2.838    0.459)    3.099
   0.786   (   0.262    1.429   -2.634)    3.009
   1.009   (   1.957    0.840    1.795)    2.785
   1.046   (   1.808    0.327    0.038)    1.838
======================= Grid point 760 (101/135) =======================
q-point: ( 0.40  0.13  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.576   (  -0.287    0.991    1.263)    1.631
   0.665   (   0.572    0.547   -2.486)    2.609
   0.707   (  -1.684    2.775    1.742)    3.684
   0.800   (  -0.670    1.952   -2.411)    3.174
   1.036   (   1.034   -0.381    1.642)    1.978
   1.071   (   1.060   -0.478    0.145)    1.172
======================= Grid point 772 (102/135) =======================
q-point: ( 0.20  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.557   (   0.744    1.289    1.611)    2.193
   0.576   (   1.786    3.094    0.767)    3.654
   0.711   (   1.960    3.394   -0.054)    3.920
   0.752   (   0.545    0.944   -2.561)    2.783
   0.944   (   1.679    2.908    2.177)    4.002
   0.999   (   0.701    1.214   -0.585)    1.519
======================= Grid point 773 (103/135) =======================
q-point: ( 0.27  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.577   (   0.302    1.206    1.097)    1.658
   0.629   (   1.406    1.825   -0.738)    2.419
   0.753   (  -0.831    3.001   -0.283)    3.127
   0.785   (   1.945    1.836   -1.397)    3.017
   0.991   (   1.335    1.434    1.666)    2.572
   1.023   (   1.172    0.394    0.117)    1.242
======================= Grid point 774 (104/135) =======================
q-point: ( 0.33  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.590   (  -0.079    1.034    0.442)    1.128
   0.656   (   0.744    0.256   -2.395)    2.521
   0.764   (  -1.282    3.170    1.313)    3.663
   0.829   (   0.286    2.919   -1.553)    3.319
   1.015   (   0.879   -0.140    1.587)    1.820
   1.044   (   1.359   -0.613    0.280)    1.517
======================= Grid point 775 (105/135) =======================
q-point: ( 0.40  0.20  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.594   (  -0.501    0.868    0.143)    1.012
   0.663   (   0.328   -0.568   -2.986)    3.057
   0.769   (  -1.905    3.300    2.013)    4.309
   0.846   (  -1.520    2.633   -1.726)    3.496
   1.021   (   0.584   -1.011    1.663)    2.032
   1.053   (   0.772   -1.338    0.244)    1.564
======================= Grid point 789 (106/135) =======================
q-point: ( 0.27  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.596   (   0.495    0.857    0.095)    0.994
   0.650   (   0.209    0.361   -2.261)    2.299
   0.795   (   1.018    1.764    0.836)    2.201
   0.840   (   1.817    3.147   -0.865)    3.735
   1.009   (   0.233    0.403    1.230)    1.315
   1.024   (  -0.089   -0.154    0.558)    0.586
======================= Grid point 790 (107/135) =======================
q-point: ( 0.33  0.27  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 900
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.609   (  -0.101    1.046   -0.777)    1.307
   0.649   (   0.240   -0.802   -2.609)    2.740
   0.822   (  -0.533    2.536    1.769)    3.137
   0.889   (   0.120    2.937   -1.383)    3.248
   1.006   (  -0.123   -0.676    1.413)    1.571
   1.023   (   0.865   -1.428    0.586)    1.769
======================= Grid point 805 (108/135) =======================
q-point: ( 0.33  0.33  0.38)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 184
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.630   (  -0.542   -0.938   -1.886)    2.175
   0.630   (   0.542    0.938   -1.886)    2.175
   0.849   (  -0.000   -0.000    1.881)    1.881
   0.924   (  -0.000   -0.000   -1.803)    1.803
   1.000   (  -0.238   -0.412    1.208)    1.299
   1.000   (   0.238    0.412    1.208)    1.299
======================= Grid point 964 (109/135) =======================
q-point: ( 0.00  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 81
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.323   (  -0.000    0.000   -3.027)    3.027
   0.323   (   0.000   -0.000   -3.027)    3.027
   0.323   (   0.000    0.000    3.027)    3.027
   0.323   (  -0.000    0.000    3.027)    3.027
   0.858   (   0.000    0.000   -6.464)    6.464
   0.858   (   0.000    0.000    6.464)    6.464
======================= Grid point 966 (110/135) =======================
q-point: ( 0.07  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.341   (   1.544    0.892   -2.986)    3.478
   0.341   (   1.544    0.892    2.986)    3.478
   0.369   (   3.790    2.188   -2.786)    5.188
   0.369   (   3.790    2.188    2.786)    5.188
   0.857   (  -0.026   -0.015   -6.452)    6.452
   0.857   (  -0.026   -0.015    6.452)    6.452
======================= Grid point 968 (111/135) =======================
q-point: ( 0.13  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.387   (   2.448    1.413   -2.887)    4.040
   0.387   (   2.448    1.413    2.887)    4.040
   0.472   (   5.018    2.897   -2.571)    6.339
   0.472   (   5.018    2.897    2.571)    6.339
   0.859   (   0.333    0.192   -6.155)    6.167
   0.859   (   0.333    0.192    6.155)    6.167
======================= Grid point 970 (112/135) =======================
q-point: ( 0.20  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.445   (   2.669    1.541   -2.787)    4.155
   0.445   (   2.669    1.541    2.787)    4.155
   0.578   (   4.193    2.421   -3.017)    5.705
   0.578   (   4.193    2.421    3.017)    5.705
   0.877   (   1.351    0.780   -5.008)    5.245
   0.877   (   1.351    0.780    5.008)    5.245
======================= Grid point 972 (113/135) =======================
q-point: ( 0.27  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.503   (   2.428    1.402   -2.739)    3.920
   0.503   (   2.428    1.402    2.739)    3.920
   0.652   (   2.328    1.344   -3.816)    4.667
   0.652   (   2.328    1.344    3.816)    4.667
   0.922   (   2.629    1.518   -3.129)    4.360
   0.922   (   2.629    1.518    3.129)    4.360
======================= Grid point 974 (114/135) =======================
q-point: ( 0.33  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.552   (   1.905    1.100   -2.757)    3.527
   0.552   (   1.905    1.100    2.757)    3.527
   0.686   (   0.907    0.524   -3.825)    3.965
   0.686   (   0.907    0.524    3.825)    3.965
   0.987   (   2.914    1.682   -1.770)    3.801
   0.987   (   2.914    1.682    1.770)    3.801
======================= Grid point 976 (115/135) =======================
q-point: ( 0.40  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.587   (   1.210    0.698   -2.810)    3.138
   0.587   (   1.210    0.698    2.810)    3.138
   0.698   (   0.258    0.149   -3.443)    3.456
   0.698   (   0.258    0.149    3.443)    3.456
   1.044   (   2.101    1.213   -1.189)    2.702
   1.044   (   2.101    1.213    1.189)    2.702
======================= Grid point 978 (116/135) =======================
q-point: ( 0.47  0.00 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.605   (   0.414    0.239   -2.852)    2.892
   0.605   (   0.414    0.239    2.852)    2.892
   0.701   (   0.039    0.023   -3.380)    3.380
   0.701   (   0.039    0.023    3.380)    3.380
   1.077   (   0.750    0.433   -1.014)    1.334
   1.077   (   0.750    0.433    1.014)    1.334
======================= Grid point 996 (117/135) =======================
q-point: ( 0.07  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.374   (   1.362    2.359   -2.440)    3.657
   0.374   (   1.362    2.359    2.440)    3.657
   0.442   (   2.814    4.874   -2.248)    6.060
   0.442   (   2.814    4.874    2.248)    6.060
   0.858   (   0.103    0.178   -6.294)    6.297
   0.858   (   0.103    0.178    6.294)    6.297
======================= Grid point 998 (118/135) =======================
q-point: ( 0.13  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.427   (   1.970    2.584   -2.257)    3.956
   0.427   (   1.970    2.584    2.257)    3.956
   0.542   (   3.579    3.973   -2.389)    5.856
   0.542   (   3.579    3.973    2.389)    5.856
   0.868   (   0.681    0.764   -5.540)    5.634
   0.868   (   0.681    0.764    5.540)    5.634
======================= Grid point 1000 (119/135) =======================
q-point: ( 0.20  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.485   (   1.970    2.493   -2.269)    3.904
   0.485   (   1.970    2.493    2.269)    3.904
   0.627   (   2.672    2.507   -3.164)    4.841
   0.627   (   2.672    2.507    3.164)    4.841
   0.900   (   1.812    1.578   -3.896)    4.577
   0.900   (   1.812    1.578    3.896)    4.577
======================= Grid point 1002 (120/135) =======================
q-point: ( 0.27  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.538   (   1.604    2.130   -2.341)    3.548
   0.538   (   1.604    2.130    2.341)    3.548
   0.678   (   1.220    1.445   -3.586)    4.054
   0.678   (   1.220    1.445    3.586)    4.054
   0.956   (   2.678    1.874   -2.223)    3.953
   0.956   (   2.678    1.874    2.223)    3.953
======================= Grid point 1004 (121/135) =======================
q-point: ( 0.33  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.579   (   1.065    1.523   -2.501)    3.116
   0.579   (   1.065    1.523    2.501)    3.116
   0.702   (   0.192    1.190   -3.292)    3.506
   0.702   (   0.192    1.190    3.292)    3.506
   1.017   (   2.515    1.367   -1.353)    3.166
   1.017   (   2.515    1.367    1.353)    3.166
======================= Grid point 1006 (122/135) =======================
q-point: ( 0.40  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.603   (   0.479    0.808   -2.696)    2.855
   0.603   (   0.479    0.808    2.696)    2.855
   0.711   (  -0.427    1.325   -3.118)    3.414
   0.711   (  -0.427    1.325    3.118)    3.414
   1.061   (   1.569    0.495   -1.059)    1.956
   1.061   (   1.569    0.495    1.059)    1.956
======================= Grid point 1008 (123/135) =======================
q-point: ( 0.47  0.07 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.611   (  -0.151    0.261   -2.790)    2.806
   0.611   (  -0.151    0.261    2.790)    2.806
   0.714   (  -0.781    1.353   -3.142)    3.509
   0.714   (  -0.781    1.353    3.142)    3.509
   1.077   (   0.229   -0.396   -1.006)    1.105
   1.077   (   0.229   -0.396    1.006)    1.105
======================= Grid point 1030 (124/135) =======================
q-point: ( 0.13  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.482   (   1.701    2.947   -1.469)    3.706
   0.482   (   1.701    2.947    1.469)    3.706
   0.616   (   2.021    3.500   -2.458)    4.731
   0.616   (   2.021    3.500    2.458)    4.731
   0.892   (   1.039    1.799   -4.190)    4.676
   0.892   (   1.039    1.799    4.190)    4.676
======================= Grid point 1032 (125/135) =======================
q-point: ( 0.20  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.536   (   1.482    2.592   -1.127)    3.191
   0.536   (   1.482    2.592    1.127)    3.191
   0.673   (   1.480    2.258   -2.552)    3.715
   0.673   (   1.480    2.258    2.552)    3.715
   0.937   (   1.866    2.166   -2.547)    3.829
   0.937   (   1.866    2.166    2.547)    3.829
======================= Grid point 1034 (126/135) =======================
q-point: ( 0.27  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.579   (   0.960    1.823   -1.564)    2.587
   0.579   (   0.960    1.823    1.564)    2.587
   0.712   (   0.641    2.106   -2.308)    3.190
   0.712   (   0.641    2.106    2.308)    3.190
   0.991   (   2.200    1.622   -1.453)    3.096
   0.991   (   2.200    1.622    1.453)    3.096
======================= Grid point 1036 (127/135) =======================
q-point: ( 0.33  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.603   (   0.472    0.892   -2.176)    2.398
   0.603   (   0.472    0.892    2.176)    2.398
   0.737   (  -0.258    2.528   -2.353)    3.463
   0.737   (  -0.258    2.528    2.353)    3.463
   1.037   (   1.869    0.593   -1.043)    2.221
   1.037   (   1.869    0.593    1.043)    2.221
======================= Grid point 1038 (128/135) =======================
q-point: ( 0.40  0.13 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.614   (   0.120    0.255   -2.485)    2.501
   0.614   (   0.120    0.255    2.485)    2.501
   0.751   (  -1.127    2.705   -2.548)    3.883
   0.751   (  -1.127    2.705    2.548)    3.883
   1.062   (   1.022   -0.408   -0.994)    1.483
   1.062   (   1.022   -0.408    0.994)    1.483
======================= Grid point 1062 (129/135) =======================
q-point: ( 0.20  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.580   (   1.025    1.776   -0.555)    2.125
   0.580   (   1.025    1.776    0.555)    2.125
   0.717   (   1.327    2.298   -0.998)    2.835
   0.717   (   1.327    2.298    0.998)    2.835
   0.979   (   1.230    2.130   -1.463)    2.862
   0.979   (   1.230    2.130    1.463)    2.862
======================= Grid point 1064 (130/135) =======================
q-point: ( 0.27  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.604   (   0.443    0.789   -1.581)    1.821
   0.604   (   0.443    0.789    1.581)    1.821
   0.761   (   0.802    2.914   -1.015)    3.188
   0.761   (   0.802    2.914    1.015)    3.188
   1.016   (   1.284    0.915   -0.856)    1.794
   1.016   (   1.284    0.915    0.856)    1.794
======================= Grid point 1066 (131/135) =======================
q-point: ( 0.33  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.613   (   0.161    0.149   -1.886)    1.899
   0.613   (   0.161    0.149    1.886)    1.899
   0.796   (  -0.398    3.325   -1.837)    3.819
   0.796   (  -0.398    3.325    1.837)    3.819
   1.039   (   1.128   -0.365   -0.800)    1.430
   1.039   (   1.128   -0.365    0.800)    1.430
======================= Grid point 1068 (132/135) =======================
q-point: ( 0.40  0.20 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.615   (   0.046   -0.080   -1.921)    1.923
   0.615   (   0.046   -0.080    1.921)    1.923
   0.809   (  -1.747    3.026   -2.091)    4.072
   0.809   (  -1.747    3.026    2.091)    4.072
   1.046   (   0.657   -1.138   -0.884)    1.584
   1.046   (   0.657   -1.138    0.884)    1.584
======================= Grid point 1096 (133/135) =======================
q-point: ( 0.27  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 304
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.612   (   0.083    0.144   -1.509)    1.519
   0.612   (   0.083    0.144    1.509)    1.519
   0.817   (   1.556    2.696   -1.319)    3.381
   0.817   (   1.556    2.696    1.319)    3.381
   1.026   (   0.081    0.140   -0.414)    0.445
   1.026   (   0.081    0.140    0.414)    0.445
======================= Grid point 1098 (134/135) =======================
q-point: ( 0.33  0.27 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 563
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.613   (   0.032   -0.081   -1.107)    1.110
   0.613   (   0.032   -0.081    1.107)    1.110
   0.857   (  -0.228    2.724   -1.839)    3.295
   0.857   (  -0.228    2.724    1.839)    3.295
   1.024   (   0.410   -0.978   -0.478)    1.164
   1.024   (   0.410   -0.978    0.478)    1.164
======================= Grid point 1128 (135/135) =======================
q-point: ( 0.33  0.33 -0.50)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 3.51e-05 3.51e-05 
Number of triplets: 119
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.612   (  -0.000   -0.000    0.000)    0.000
   0.612   (  -0.000   -0.000    0.000)    0.000
   0.887   (  -0.000   -0.000   -2.033)    2.033
   0.887   (  -0.000   -0.000    2.033)    2.033
   1.012   (  -0.000   -0.000    0.000)    0.000
   1.012   (  -0.000   -0.000    0.000)    0.000
=================== End of collection of collisions ===================
----------- Thermal conductivity (W/m-k) with tetrahedron method -----------
#  T(K)        xx         yy         zz         yz         xz         xy        #ipm
    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/10800
   10.0      0.097      0.097      0.080      0.000     -0.000     -0.000 3/10800
   20.0      0.060      0.060      0.056      0.000     -0.000     -0.000 3/10800
   30.0      0.043      0.043      0.041      0.000     -0.000     -0.000 3/10800
   40.0      0.033      0.033      0.032      0.000     -0.000     -0.000 3/10800
   50.0      0.026      0.026      0.026      0.000     -0.000     -0.000 3/10800
   60.0      0.022      0.022      0.022      0.000     -0.000     -0.000 3/10800
   70.0      0.019      0.019      0.019      0.000     -0.000     -0.000 3/10800
   80.0      0.017      0.017      0.016      0.000     -0.000     -0.000 3/10800
   90.0      0.015      0.015      0.015      0.000     -0.000     -0.000 3/10800
  100.0      0.013      0.013      0.013      0.000     -0.000     -0.000 3/10800
  110.0      0.012      0.012      0.012      0.000     -0.000     -0.000 3/10800
  120.0      0.011      0.011      0.011      0.000     -0.000     -0.000 3/10800
  130.0      0.010      0.010      0.010      0.000     -0.000     -0.000 3/10800
  140.0      0.010      0.010      0.009      0.000     -0.000     -0.000 3/10800
  150.0      0.009      0.009      0.009      0.000     -0.000     -0.000 3/10800
  160.0      0.008      0.008      0.008      0.000     -0.000     -0.000 3/10800
  170.0      0.008      0.008      0.008      0.000     -0.000     -0.000 3/10800
  180.0      0.007      0.007      0.007      0.000     -0.000     -0.000 3/10800
  190.0      0.007      0.007      0.007      0.000     -0.000     -0.000 3/10800
  200.0      0.007      0.007      0.007      0.000     -0.000     -0.000 3/10800
  210.0      0.006      0.006      0.006      0.000     -0.000     -0.000 3/10800
  220.0      0.006      0.006      0.006      0.000     -0.000     -0.000 3/10800
  230.0      0.006      0.006      0.006      0.000     -0.000     -0.000 3/10800
  240.0      0.006      0.006      0.006      0.000     -0.000     -0.000 3/10800
  250.0      0.005      0.005      0.005      0.000     -0.000     -0.000 3/10800
  260.0      0.005      0.005      0.005      0.000     -0.000     -0.000 3/10800
  270.0      0.005      0.005      0.005      0.000     -0.000     -0.000 3/10800
  280.0      0.005      0.005      0.005      0.000     -0.000     -0.000 3/10800
  290.0      0.005      0.005      0.005      0.000     -0.000     -0.000 3/10800
  300.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  310.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  320.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  330.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  340.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  350.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  360.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  370.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  380.0      0.004      0.004      0.004      0.000     -0.000     -0.000 3/10800
  390.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  400.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  410.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  420.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  430.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  440.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  450.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  460.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  470.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  480.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  490.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  500.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  510.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  520.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  530.0      0.003      0.003      0.003      0.000     -0.000     -0.000 3/10800
  540.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  550.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  560.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  570.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  580.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  590.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  600.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  610.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  620.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  630.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  640.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  650.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  660.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  670.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  680.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  690.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  700.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  710.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  720.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  730.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  740.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  750.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  760.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  770.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  780.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  790.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  800.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  810.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  820.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  830.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  840.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  850.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  860.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  870.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  880.0      0.002      0.002      0.002      0.000     -0.000     -0.000 3/10800
  890.0      0.002      0.002      0.001      0.000     -0.000     -0.000 3/10800
  900.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  910.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  920.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  930.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  940.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  950.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  960.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  970.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  980.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
  990.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800
 1000.0      0.001      0.001      0.001      0.000     -0.000     -0.000 3/10800

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

