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

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

-------------------------[time 2026-01-08 05:41:49]-------------------------
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: [2 2 2]
  Primitive matrix:
    [0.  0.5 0.5]
    [0.5 0.  0.5]
    [0.5 0.5 0. ]
Spacegroup: Fm-3m (225)
Number of symmetry operations in supercell: 1536
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    0.000000000000000    2.377529000000000    2.377529000000000
  b    2.377529000000000    0.000000000000000    2.377529000000000
  c    2.377529000000000    2.377529000000000    0.000000000000000
Atomic positions (fractional):
   *1 H   0.50000000000000  0.50000000000000  0.50000000000000   1.008
   *2 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990
-------------------------------- unit cell ---------------------------------
Lattice vectors:
  a    4.755058000000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    4.755058000000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    4.755058000000000
Atomic positions (fractional):
   *1 H   0.50000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    2 H   0.50000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    3 H   0.00000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    4 H   0.00000000000000  0.50000000000000  0.00000000000000   1.008 > 1
   *5 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990 > 2
    6 Na  0.00000000000000  0.50000000000000  0.50000000000000  22.990 > 2
    7 Na  0.50000000000000  0.00000000000000  0.50000000000000  22.990 > 2
    8 Na  0.50000000000000  0.50000000000000  0.00000000000000  22.990 > 2
-------------------------------- super cell --------------------------------
Lattice vectors:
  a    9.510116000000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    9.510116000000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.510116000000000
Atomic positions (fractional):
   *1 H   0.25000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    2 H   0.75000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    3 H   0.25000000000000  0.50000000000000  0.00000000000000   1.008 > 1
    4 H   0.75000000000000  0.50000000000000  0.00000000000000   1.008 > 1
    5 H   0.25000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    6 H   0.75000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    7 H   0.25000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    8 H   0.75000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    9 H   0.25000000000000  0.25000000000000  0.25000000000000   1.008 > 1
   10 H   0.75000000000000  0.25000000000000  0.25000000000000   1.008 > 1
   11 H   0.25000000000000  0.75000000000000  0.25000000000000   1.008 > 1
   12 H   0.75000000000000  0.75000000000000  0.25000000000000   1.008 > 1
   13 H   0.25000000000000  0.25000000000000  0.75000000000000   1.008 > 1
   14 H   0.75000000000000  0.25000000000000  0.75000000000000   1.008 > 1
   15 H   0.25000000000000  0.75000000000000  0.75000000000000   1.008 > 1
   16 H   0.75000000000000  0.75000000000000  0.75000000000000   1.008 > 1
   17 H   0.00000000000000  0.00000000000000  0.25000000000000   1.008 > 1
   18 H   0.50000000000000  0.00000000000000  0.25000000000000   1.008 > 1
   19 H   0.00000000000000  0.50000000000000  0.25000000000000   1.008 > 1
   20 H   0.50000000000000  0.50000000000000  0.25000000000000   1.008 > 1
   21 H   0.00000000000000  0.00000000000000  0.75000000000000   1.008 > 1
   22 H   0.50000000000000  0.00000000000000  0.75000000000000   1.008 > 1
   23 H   0.00000000000000  0.50000000000000  0.75000000000000   1.008 > 1
   24 H   0.50000000000000  0.50000000000000  0.75000000000000   1.008 > 1
   25 H   0.00000000000000  0.25000000000000  0.00000000000000   1.008 > 1
   26 H   0.50000000000000  0.25000000000000  0.00000000000000   1.008 > 1
   27 H   0.00000000000000  0.75000000000000  0.00000000000000   1.008 > 1
   28 H   0.50000000000000  0.75000000000000  0.00000000000000   1.008 > 1
   29 H   0.00000000000000  0.25000000000000  0.50000000000000   1.008 > 1
   30 H   0.50000000000000  0.25000000000000  0.50000000000000   1.008 > 1
   31 H   0.00000000000000  0.75000000000000  0.50000000000000   1.008 > 1
   32 H   0.50000000000000  0.75000000000000  0.50000000000000   1.008 > 1
  *33 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990 > 2
   34 Na  0.50000000000000  0.00000000000000  0.00000000000000  22.990 > 2
   35 Na  0.00000000000000  0.50000000000000  0.00000000000000  22.990 > 2
   36 Na  0.50000000000000  0.50000000000000  0.00000000000000  22.990 > 2
   37 Na  0.00000000000000  0.00000000000000  0.50000000000000  22.990 > 2
   38 Na  0.50000000000000  0.00000000000000  0.50000000000000  22.990 > 2
   39 Na  0.00000000000000  0.50000000000000  0.50000000000000  22.990 > 2
   40 Na  0.50000000000000  0.50000000000000  0.50000000000000  22.990 > 2
   41 Na  0.00000000000000  0.25000000000000  0.25000000000000  22.990 > 2
   42 Na  0.50000000000000  0.25000000000000  0.25000000000000  22.990 > 2
   43 Na  0.00000000000000  0.75000000000000  0.25000000000000  22.990 > 2
   44 Na  0.50000000000000  0.75000000000000  0.25000000000000  22.990 > 2
   45 Na  0.00000000000000  0.25000000000000  0.75000000000000  22.990 > 2
   46 Na  0.50000000000000  0.25000000000000  0.75000000000000  22.990 > 2
   47 Na  0.00000000000000  0.75000000000000  0.75000000000000  22.990 > 2
   48 Na  0.50000000000000  0.75000000000000  0.75000000000000  22.990 > 2
   49 Na  0.25000000000000  0.00000000000000  0.25000000000000  22.990 > 2
   50 Na  0.75000000000000  0.00000000000000  0.25000000000000  22.990 > 2
   51 Na  0.25000000000000  0.50000000000000  0.25000000000000  22.990 > 2
   52 Na  0.75000000000000  0.50000000000000  0.25000000000000  22.990 > 2
   53 Na  0.25000000000000  0.00000000000000  0.75000000000000  22.990 > 2
   54 Na  0.75000000000000  0.00000000000000  0.75000000000000  22.990 > 2
   55 Na  0.25000000000000  0.50000000000000  0.75000000000000  22.990 > 2
   56 Na  0.75000000000000  0.50000000000000  0.75000000000000  22.990 > 2
   57 Na  0.25000000000000  0.25000000000000  0.00000000000000  22.990 > 2
   58 Na  0.75000000000000  0.25000000000000  0.00000000000000  22.990 > 2
   59 Na  0.25000000000000  0.75000000000000  0.00000000000000  22.990 > 2
   60 Na  0.75000000000000  0.75000000000000  0.00000000000000  22.990 > 2
   61 Na  0.25000000000000  0.25000000000000  0.50000000000000  22.990 > 2
   62 Na  0.75000000000000  0.25000000000000  0.50000000000000  22.990 > 2
   63 Na  0.25000000000000  0.75000000000000  0.50000000000000  22.990 > 2
   64 Na  0.75000000000000  0.75000000000000  0.50000000000000  22.990 > 2
----------------------------------------------------------------------------
NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.2506430    0.0000000    0.0000000
            0.0000000    3.2506430    0.0000000
            0.0000000    0.0000000    3.2506430
-------------------------- Born effective charges --------------------------
    1 H    -0.9497898    0.0000000    0.0000000
            0.0000000   -0.9497898    0.0000000
            0.0000000    0.0000000   -0.9497898
    2 Na    0.9497898    0.0000000    0.0000000
            0.0000000    0.9497898    0.0000000
            0.0000000    0.0000000    0.9497898
----------------------------------------------------------------------------
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: 192/192
Permutation basis: 852/852
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 33
Number of blocks in projector: 33
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: 28
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 3
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 2
Use standard eigh solver.
Tree of FC basis block matrices:
- (33, 31), data: False
|-- (2, 2), data: True
|-- (3, 3), data: True
|-- (28, 26), data: True
-----
Solver_atoms: 1 -- 64 / 64
Time (Solver_compr_matrix_reshape): 0.001
Solver_block: 80 / 80
 - Time: 0.006
Solver: Calculate X.T @ X and X.T @ y
 (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.007
--------------------------------- Symfc end --------------------------------
Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) 
Permutation basis: 192/192
Permutation basis: 852/852
Construct permutation basis matrix.
Finding block diagonal structure in projector.
Using scipy connected_components.
Rank of projector: 33
Number of blocks in projector: 33
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: 28
Use standard eigh solver.
--- Eigsh_solver_block: 2 / 3 ---
Block_size: 3
Use standard eigh solver.
--- Eigsh_solver_block: 3 / 3 ---
Block_size: 2
Use standard eigh solver.
Tree of FC basis block matrices:
- (33, 31), data: False
|-- (2, 2), data: True
|-- (3, 3), data: True
|-- (28, 26), data: True
Max drift after symmetrization by symfc projector: 0.00000000 (xx) 0.00000000 (xx) 
Force constants are written into "force_constants.hdf5".

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

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


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

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

-------------------------[time 2026-01-08 05:41:52]-------------------------
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): [2 2 2]
Primitive matrix:
  [0.  0.5 0.5]
  [0.5 0.  0.5]
  [0.5 0.5 0. ]
Spacegroup: Fm-3m (225)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    0.000000000000000    2.377529000000000    2.377529000000000
  b    2.377529000000000    0.000000000000000    2.377529000000000
  c    2.377529000000000    2.377529000000000    0.000000000000000
Atomic positions (fractional):
    1 H   0.50000000000000  0.50000000000000  0.50000000000000   1.008
    2 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a    9.510116000000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    9.510116000000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.510116000000000
Atomic positions (fractional):
    1 H   0.25000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    2 H   0.75000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    3 H   0.25000000000000  0.50000000000000  0.00000000000000   1.008 > 1
    4 H   0.75000000000000  0.50000000000000  0.00000000000000   1.008 > 1
    5 H   0.25000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    6 H   0.75000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    7 H   0.25000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    8 H   0.75000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    9 H   0.25000000000000  0.25000000000000  0.25000000000000   1.008 > 1
   10 H   0.75000000000000  0.25000000000000  0.25000000000000   1.008 > 1
   11 H   0.25000000000000  0.75000000000000  0.25000000000000   1.008 > 1
   12 H   0.75000000000000  0.75000000000000  0.25000000000000   1.008 > 1
   13 H   0.25000000000000  0.25000000000000  0.75000000000000   1.008 > 1
   14 H   0.75000000000000  0.25000000000000  0.75000000000000   1.008 > 1
   15 H   0.25000000000000  0.75000000000000  0.75000000000000   1.008 > 1
   16 H   0.75000000000000  0.75000000000000  0.75000000000000   1.008 > 1
   17 H   0.00000000000000  0.00000000000000  0.25000000000000   1.008 > 1
   18 H   0.50000000000000  0.00000000000000  0.25000000000000   1.008 > 1
   19 H   0.00000000000000  0.50000000000000  0.25000000000000   1.008 > 1
   20 H   0.50000000000000  0.50000000000000  0.25000000000000   1.008 > 1
   21 H   0.00000000000000  0.00000000000000  0.75000000000000   1.008 > 1
   22 H   0.50000000000000  0.00000000000000  0.75000000000000   1.008 > 1
   23 H   0.00000000000000  0.50000000000000  0.75000000000000   1.008 > 1
   24 H   0.50000000000000  0.50000000000000  0.75000000000000   1.008 > 1
   25 H   0.00000000000000  0.25000000000000  0.00000000000000   1.008 > 1
   26 H   0.50000000000000  0.25000000000000  0.00000000000000   1.008 > 1
   27 H   0.00000000000000  0.75000000000000  0.00000000000000   1.008 > 1
   28 H   0.50000000000000  0.75000000000000  0.00000000000000   1.008 > 1
   29 H   0.00000000000000  0.25000000000000  0.50000000000000   1.008 > 1
   30 H   0.50000000000000  0.25000000000000  0.50000000000000   1.008 > 1
   31 H   0.00000000000000  0.75000000000000  0.50000000000000   1.008 > 1
   32 H   0.50000000000000  0.75000000000000  0.50000000000000   1.008 > 1
   33 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990 > 33
   34 Na  0.50000000000000  0.00000000000000  0.00000000000000  22.990 > 33
   35 Na  0.00000000000000  0.50000000000000  0.00000000000000  22.990 > 33
   36 Na  0.50000000000000  0.50000000000000  0.00000000000000  22.990 > 33
   37 Na  0.00000000000000  0.00000000000000  0.50000000000000  22.990 > 33
   38 Na  0.50000000000000  0.00000000000000  0.50000000000000  22.990 > 33
   39 Na  0.00000000000000  0.50000000000000  0.50000000000000  22.990 > 33
   40 Na  0.50000000000000  0.50000000000000  0.50000000000000  22.990 > 33
   41 Na  0.00000000000000  0.25000000000000  0.25000000000000  22.990 > 33
   42 Na  0.50000000000000  0.25000000000000  0.25000000000000  22.990 > 33
   43 Na  0.00000000000000  0.75000000000000  0.25000000000000  22.990 > 33
   44 Na  0.50000000000000  0.75000000000000  0.25000000000000  22.990 > 33
   45 Na  0.00000000000000  0.25000000000000  0.75000000000000  22.990 > 33
   46 Na  0.50000000000000  0.25000000000000  0.75000000000000  22.990 > 33
   47 Na  0.00000000000000  0.75000000000000  0.75000000000000  22.990 > 33
   48 Na  0.50000000000000  0.75000000000000  0.75000000000000  22.990 > 33
   49 Na  0.25000000000000  0.00000000000000  0.25000000000000  22.990 > 33
   50 Na  0.75000000000000  0.00000000000000  0.25000000000000  22.990 > 33
   51 Na  0.25000000000000  0.50000000000000  0.25000000000000  22.990 > 33
   52 Na  0.75000000000000  0.50000000000000  0.25000000000000  22.990 > 33
   53 Na  0.25000000000000  0.00000000000000  0.75000000000000  22.990 > 33
   54 Na  0.75000000000000  0.00000000000000  0.75000000000000  22.990 > 33
   55 Na  0.25000000000000  0.50000000000000  0.75000000000000  22.990 > 33
   56 Na  0.75000000000000  0.50000000000000  0.75000000000000  22.990 > 33
   57 Na  0.25000000000000  0.25000000000000  0.00000000000000  22.990 > 33
   58 Na  0.75000000000000  0.25000000000000  0.00000000000000  22.990 > 33
   59 Na  0.25000000000000  0.75000000000000  0.00000000000000  22.990 > 33
   60 Na  0.75000000000000  0.75000000000000  0.00000000000000  22.990 > 33
   61 Na  0.25000000000000  0.25000000000000  0.50000000000000  22.990 > 33
   62 Na  0.75000000000000  0.25000000000000  0.50000000000000  22.990 > 33
   63 Na  0.25000000000000  0.75000000000000  0.50000000000000  22.990 > 33
   64 Na  0.75000000000000  0.75000000000000  0.50000000000000  22.990 > 33
----------------------------------------------------------------------------
NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
--------------------------- Dielectric constant ----------------------------
            3.2506430    0.0000000    0.0000000
            0.0000000    3.2506430    0.0000000
            0.0000000    0.0000000    3.2506430
-------------------------- Born effective charges --------------------------
    1 H    -0.9497898    0.0000000    0.0000000
            0.0000000   -0.9497898    0.0000000
            0.0000000    0.0000000   -0.9497898
    2 Na    0.9497898    0.0000000    0.0000000
            0.0000000    0.9497898    0.0000000
            0.0000000    0.0000000    0.9497898
----------------------------------------------------------------------------
Sets of supercell forces were read from "FORCES_FC3.xz".
Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".
----------------------------- Force constants ------------------------------
Computing fc3[ 1, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Computing fc3[ 33, x, x ] using numpy.linalg.pinv.
Displacements (in Angstrom):
    [ 0.0100  0.0000  0.0000]
    [-0.0100  0.0000  0.0000]
Expanding fc3.
Symmetrizing fc3 by traditional approach (N=3).
Symmetrizing fc2 by traditional approach (N=3).
Max drift of fc3: 0.00000005 (yyy) 0.00000005 (yyy) 0.00000005 (yyy)
fc3 was written into "fc3.hdf5".
Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 
fc2 was written into "fc2.hdf5".
--------------------------- Calculation settings ---------------------------
Non-analytical term correction (NAC): True
NAC unit conversion factor:  14.39965
BZ integration: Tetrahedron-method
Temperatures: 0.0  300.0 
Cutoff frequency: 0.01
Frequency conversion factor to THz:  15.63330
----------- None of ph-ph interaction calculation was performed. -----------
Summary of calculation was written in "phono3py.yaml".
-------------------------[time 2026-01-08 05:41:53]-------------------------
                 _
   ___ _ __   __| |
  / _ \ '_ \ / _` |
 |  __/ | | | (_| |
  \___|_| |_|\__,_|


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

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

-------------------------[time 2026-01-08 05:41:53]-------------------------
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): [2 2 2]
Primitive matrix:
  [0.  0.5 0.5]
  [0.5 0.  0.5]
  [0.5 0.5 0. ]
Spacegroup: Fm-3m (225)
------------------------------ primitive cell ------------------------------
Lattice vectors:
  a    0.000000000000000    2.377529000000000    2.377529000000000
  b    2.377529000000000    0.000000000000000    2.377529000000000
  c    2.377529000000000    2.377529000000000    0.000000000000000
Atomic positions (fractional):
    1 H   0.50000000000000  0.50000000000000  0.50000000000000   1.008
    2 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990
-------------------------------- supercell ---------------------------------
Lattice vectors:
  a    9.510116000000000    0.000000000000000    0.000000000000000
  b    0.000000000000000    9.510116000000000    0.000000000000000
  c    0.000000000000000    0.000000000000000    9.510116000000000
Atomic positions (fractional):
    1 H   0.25000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    2 H   0.75000000000000  0.00000000000000  0.00000000000000   1.008 > 1
    3 H   0.25000000000000  0.50000000000000  0.00000000000000   1.008 > 1
    4 H   0.75000000000000  0.50000000000000  0.00000000000000   1.008 > 1
    5 H   0.25000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    6 H   0.75000000000000  0.00000000000000  0.50000000000000   1.008 > 1
    7 H   0.25000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    8 H   0.75000000000000  0.50000000000000  0.50000000000000   1.008 > 1
    9 H   0.25000000000000  0.25000000000000  0.25000000000000   1.008 > 1
   10 H   0.75000000000000  0.25000000000000  0.25000000000000   1.008 > 1
   11 H   0.25000000000000  0.75000000000000  0.25000000000000   1.008 > 1
   12 H   0.75000000000000  0.75000000000000  0.25000000000000   1.008 > 1
   13 H   0.25000000000000  0.25000000000000  0.75000000000000   1.008 > 1
   14 H   0.75000000000000  0.25000000000000  0.75000000000000   1.008 > 1
   15 H   0.25000000000000  0.75000000000000  0.75000000000000   1.008 > 1
   16 H   0.75000000000000  0.75000000000000  0.75000000000000   1.008 > 1
   17 H   0.00000000000000  0.00000000000000  0.25000000000000   1.008 > 1
   18 H   0.50000000000000  0.00000000000000  0.25000000000000   1.008 > 1
   19 H   0.00000000000000  0.50000000000000  0.25000000000000   1.008 > 1
   20 H   0.50000000000000  0.50000000000000  0.25000000000000   1.008 > 1
   21 H   0.00000000000000  0.00000000000000  0.75000000000000   1.008 > 1
   22 H   0.50000000000000  0.00000000000000  0.75000000000000   1.008 > 1
   23 H   0.00000000000000  0.50000000000000  0.75000000000000   1.008 > 1
   24 H   0.50000000000000  0.50000000000000  0.75000000000000   1.008 > 1
   25 H   0.00000000000000  0.25000000000000  0.00000000000000   1.008 > 1
   26 H   0.50000000000000  0.25000000000000  0.00000000000000   1.008 > 1
   27 H   0.00000000000000  0.75000000000000  0.00000000000000   1.008 > 1
   28 H   0.50000000000000  0.75000000000000  0.00000000000000   1.008 > 1
   29 H   0.00000000000000  0.25000000000000  0.50000000000000   1.008 > 1
   30 H   0.50000000000000  0.25000000000000  0.50000000000000   1.008 > 1
   31 H   0.00000000000000  0.75000000000000  0.50000000000000   1.008 > 1
   32 H   0.50000000000000  0.75000000000000  0.50000000000000   1.008 > 1
   33 Na  0.00000000000000  0.00000000000000  0.00000000000000  22.990 > 33
   34 Na  0.50000000000000  0.00000000000000  0.00000000000000  22.990 > 33
   35 Na  0.00000000000000  0.50000000000000  0.00000000000000  22.990 > 33
   36 Na  0.50000000000000  0.50000000000000  0.00000000000000  22.990 > 33
   37 Na  0.00000000000000  0.00000000000000  0.50000000000000  22.990 > 33
   38 Na  0.50000000000000  0.00000000000000  0.50000000000000  22.990 > 33
   39 Na  0.00000000000000  0.50000000000000  0.50000000000000  22.990 > 33
   40 Na  0.50000000000000  0.50000000000000  0.50000000000000  22.990 > 33
   41 Na  0.00000000000000  0.25000000000000  0.25000000000000  22.990 > 33
   42 Na  0.50000000000000  0.25000000000000  0.25000000000000  22.990 > 33
   43 Na  0.00000000000000  0.75000000000000  0.25000000000000  22.990 > 33
   44 Na  0.50000000000000  0.75000000000000  0.25000000000000  22.990 > 33
   45 Na  0.00000000000000  0.25000000000000  0.75000000000000  22.990 > 33
   46 Na  0.50000000000000  0.25000000000000  0.75000000000000  22.990 > 33
   47 Na  0.00000000000000  0.75000000000000  0.75000000000000  22.990 > 33
   48 Na  0.50000000000000  0.75000000000000  0.75000000000000  22.990 > 33
   49 Na  0.25000000000000  0.00000000000000  0.25000000000000  22.990 > 33
   50 Na  0.75000000000000  0.00000000000000  0.25000000000000  22.990 > 33
   51 Na  0.25000000000000  0.50000000000000  0.25000000000000  22.990 > 33
   52 Na  0.75000000000000  0.50000000000000  0.25000000000000  22.990 > 33
   53 Na  0.25000000000000  0.00000000000000  0.75000000000000  22.990 > 33
   54 Na  0.75000000000000  0.00000000000000  0.75000000000000  22.990 > 33
   55 Na  0.25000000000000  0.50000000000000  0.75000000000000  22.990 > 33
   56 Na  0.75000000000000  0.50000000000000  0.75000000000000  22.990 > 33
   57 Na  0.25000000000000  0.25000000000000  0.00000000000000  22.990 > 33
   58 Na  0.75000000000000  0.25000000000000  0.00000000000000  22.990 > 33
   59 Na  0.25000000000000  0.75000000000000  0.00000000000000  22.990 > 33
   60 Na  0.75000000000000  0.75000000000000  0.00000000000000  22.990 > 33
   61 Na  0.25000000000000  0.25000000000000  0.50000000000000  22.990 > 33
   62 Na  0.75000000000000  0.25000000000000  0.50000000000000  22.990 > 33
   63 Na  0.25000000000000  0.75000000000000  0.50000000000000  22.990 > 33
   64 Na  0.75000000000000  0.75000000000000  0.50000000000000  22.990 > 33
----------------------------------------------------------------------------
NAC parameters were read from "phono3py.yaml".
--------------------------- Dielectric constant ----------------------------
            3.2506430    0.0000000    0.0000000
            0.0000000    3.2506430    0.0000000
            0.0000000    0.0000000    3.2506430
-------------------------- Born effective charges --------------------------
    1 H    -0.9497898    0.0000000    0.0000000
            0.0000000   -0.9497898    0.0000000
            0.0000000    0.0000000   -0.9497898
    2 Na    0.9497898    0.0000000    0.0000000
            0.0000000    0.9497898    0.0000000
            0.0000000    0.0000000    0.9497898
----------------------------------------------------------------------------
fc3 was read from "fc3.hdf5".
fc2 was read from "fc2.hdf5".
----------------------------- Force constants ------------------------------
Max drift of fc3: 0.00000005 (yyy) 0.00000005 (yyy) 0.00000005 (yyy)
Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) 
--------------------------- 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 ... [ 18 18 18 ]
fc3-r2q-transformation over three atoms: True
--------------------------- Phonon calculations ----------------------------
Use NAC by Gonze et al. (no real space sum in current implementation)
  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)
  G-cutoff distance: 1.39, Number of G-points: 307, Lambda: 0.26
Running harmonic phonon calculations...
-------------------- Lattice thermal conductivity (RTA) --------------------
======================= Grid point 0 (1/195) =======================
q-point: ( 0.00  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 195
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
  16.382   (   0.000    0.000    0.000)    0.000
  16.382   (   0.000    0.000    0.000)    0.000
  16.382   (   0.000    0.000    0.000)    0.000
======================= Grid point 1 (2/195) =======================
q-point: ( 0.06  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 570
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.746   ( -20.802   20.802   20.802)   36.030
   0.746   ( -20.802   20.802   20.802)   36.030
   1.231   ( -34.749   34.749   34.749)   60.187
  16.467   (  -4.763    4.763    4.763)    8.250
  16.467   (  -4.763    4.763    4.763)    8.250
  27.257   (  -1.718    1.718    1.718)    2.976
======================= Grid point 2 (3/195) =======================
q-point: ( 0.11  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 580
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.445   ( -18.934   18.934   18.934)   32.794
   1.445   ( -18.934   18.934   18.934)   32.794
   2.424   ( -33.163   33.163   33.163)   57.440
  16.702   (  -8.353    8.353    8.353)   14.468
  16.702   (  -8.353    8.353    8.353)   14.468
  27.334   (  -2.410    2.410    2.410)    4.173
======================= Grid point 3 (4/195) =======================
q-point: ( 0.17  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 570
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.066   ( -16.485   16.485   16.485)   28.554
   2.066   ( -16.485   16.485   16.485)   28.554
   3.547   ( -30.782   30.782   30.782)   53.315
  17.031   ( -10.106   10.106   10.106)   17.504
  17.031   ( -10.106   10.106   10.106)   17.504
  27.408   (  -1.589    1.589    1.589)    2.752
======================= Grid point 4 (5/195) =======================
q-point: ( 0.22  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 580
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.599   ( -13.923   13.923   13.923)   24.116
   2.599   ( -13.923   13.923   13.923)   24.116
   4.574   ( -27.748   27.748   27.748)   48.062
  17.389   ( -10.060   10.060   10.060)   17.424
  17.389   ( -10.060   10.060   10.060)   17.424
  27.431   (   0.454   -0.454   -0.454)    0.786
======================= Grid point 5 (6/195) =======================
q-point: ( 0.28  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 570
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.043   ( -11.397   11.397   11.397)   19.741
   3.043   ( -11.397   11.397   11.397)   19.741
   5.484   ( -24.022   24.022   24.022)   41.607
  17.721   (  -8.735    8.735    8.735)   15.130
  17.721   (  -8.735    8.735    8.735)   15.130
  27.373   (   2.813   -2.813   -2.813)    4.873
======================= Grid point 6 (7/195) =======================
q-point: ( 0.33  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 580
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.397   (  -8.826    8.826    8.826)   15.287
   3.397   (  -8.826    8.826    8.826)   15.287
   6.248   ( -19.426   19.426   19.426)   33.647
  17.994   (  -6.747    6.747    6.747)   11.687
  17.994   (  -6.747    6.747    6.747)   11.687
  27.243   (   4.404   -4.404   -4.404)    7.629
======================= Grid point 7 (8/195) =======================
q-point: ( 0.39  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 570
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.659   (  -6.080    6.080    6.080)   10.531
   3.659   (  -6.080    6.080    6.080)   10.531
   6.833   ( -13.806   13.806   13.806)   23.912
  18.192   (  -4.525    4.525    4.525)    7.837
  18.192   (  -4.525    4.525    4.525)    7.837
  27.082   (   4.437   -4.437   -4.437)    7.686
======================= Grid point 8 (9/195) =======================
q-point: ( 0.44  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 580
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.821   (  -3.114    3.114    3.114)    5.394
   3.821   (  -3.114    3.114    3.114)    5.394
   7.204   (  -7.209    7.209    7.209)   12.486
  18.311   (  -2.260    2.260    2.260)    3.915
  18.311   (  -2.260    2.260    2.260)    3.915
  26.952   (   2.770   -2.770   -2.770)    4.797
======================= Grid point 9 (10/195) =======================
q-point: (-0.50  0.00  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 290
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.876   (  -0.000    0.000    0.000)    0.000
   3.876   (  -0.000    0.000    0.000)    0.000
   7.332   (  -0.000    0.000    0.000)    0.000
  18.350   (  -0.000    0.000    0.000)    0.000
  18.350   (  -0.000    0.000    0.000)    0.000
  26.901   (  -0.000    0.000    0.000)    0.000
======================= Grid point 20 (11/195) =======================
q-point: ( 0.06  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 475
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   0.887   (   0.000   -0.000   36.909)   36.909
   0.887   (   0.000   -0.000   36.909)   36.909
   1.380   (   0.000   -0.000   57.905)   57.905
  16.568   (   0.000   -0.000   15.674)   15.674
  16.568   (   0.000   -0.000   15.674)   15.674
  27.182   (  -0.000    0.000   -3.706)    3.706
======================= Grid point 21 (12/195) =======================
q-point: ( 0.11  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.416   (  -8.875    8.875   32.263)   34.618
   1.421   (  -9.225    9.225   31.372)   33.977
   2.279   ( -18.797   18.797   49.081)   55.818
  16.747   (   1.874   -1.874   20.944)   21.111
  16.823   (  -2.570    2.570   21.989)   22.287
  27.210   (  -5.685    5.685   -5.167)    9.557
======================= Grid point 22 (13/195) =======================
q-point: ( 0.17  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.995   ( -11.041   11.041   25.198)   29.644
   2.001   ( -10.784   10.784   27.260)   31.236
   3.308   ( -23.828   23.828   40.005)   52.306
  16.988   (  -1.091    1.091   21.579)   21.634
  17.181   (  -3.658    3.658   25.895)   26.406
  27.289   (  -6.107    6.107   -6.034)   10.536
======================= Grid point 23 (14/195) =======================
q-point: ( 0.22  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.513   ( -10.595   10.595   20.022)   25.008
   2.541   (  -9.899    9.899   23.421)   27.286
   4.303   ( -24.730   24.730   32.524)   47.759
  17.289   (  -3.366    3.366   20.312)   20.862
  17.575   (  -3.057    3.057   27.098)   27.441
  27.335   (  -4.401    4.401   -7.641)    9.855
======================= Grid point 24 (15/195) =======================
q-point: ( 0.28  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.956   (  -9.402    9.402   15.662)   20.545
   3.006   (  -8.043    8.043   20.227)   23.206
   5.211   ( -23.594   23.594   25.799)   42.178
  17.607   (  -3.987    3.987   18.992)   19.812
  17.938   (  -1.191    1.191   25.828)   25.883
  27.308   (  -1.917    1.917   -9.754)   10.123
======================= Grid point 25 (16/195) =======================
q-point: ( 0.33  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.318   (  -7.841    7.841   11.885)   16.255
   3.387   (  -5.788    5.788   17.188)   19.038
   5.997   ( -20.918   20.918   19.218)   35.277
  17.911   (  -3.066    3.066   18.722)   19.218
  18.222   (   1.195   -1.195   22.528)   22.591
  27.202   (   0.176   -0.176  -11.427)   11.429
======================= Grid point 26 (17/195) =======================
q-point: ( 0.39  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.598   (  -5.893    5.893    8.494)   11.900
   3.676   (  -3.271    3.271   14.012)   14.756
   6.625   ( -16.783   16.783   12.369)   26.765
  18.183   (  -0.836    0.836   19.881)   19.916
  18.404   (   3.397   -3.397   17.643)   18.285
  27.048   (   0.915   -0.915  -11.820)   11.891
======================= Grid point 27 (18/195) =======================
q-point: ( 0.44  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.789   (  -3.474    3.474    5.360)    7.271
   3.865   (  -0.530    0.530   10.670)   10.696
   7.058   ( -11.271   11.271    5.125)   16.744
  18.410   (   2.498   -2.498   22.088)   22.369
  18.476   (   4.956   -4.956   11.617)   13.568
  26.903   (  -0.122    0.122  -10.640)   10.641
======================= Grid point 28 (19/195) =======================
q-point: ( 0.50  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.883   (  -0.593    0.593    2.425)    2.566
   3.948   (   2.383   -2.383    7.361)    8.096
   7.262   (  -4.680    4.680   -2.277)    6.999
  18.447   (   5.715   -5.715    5.000)    9.504
  18.577   (   6.663   -6.663   24.450)   26.203
  26.822   (  -2.617    2.617   -8.278)    9.068
======================= Grid point 29 (20/195) =======================
q-point: (-0.44  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.872   (   2.608   -2.608   -0.319)    3.702
   3.926   (   5.389   -5.389    4.386)    8.793
   7.221   (   2.435   -2.435   -9.285)    9.903
  18.332   (   5.797   -5.797   -1.509)    8.336
  18.664   (  11.236  -11.236   25.970)   30.445
  26.837   (  -5.595    5.595   -5.649)    9.722
======================= Grid point 30 (21/195) =======================
q-point: (-0.39  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.754   (   5.902   -5.902   -2.853)    8.821
   3.801   (   8.414   -8.414    1.985)   12.064
   6.937   (   9.430   -9.430  -15.261)   20.267
  18.147   (   5.556   -5.556   -7.159)   10.630
  18.653   (  15.716  -15.716   25.929)   34.151
  26.942   (  -7.853    7.853   -3.770)   11.728
======================= Grid point 31 (22/195) =======================
q-point: (-0.33  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.531   (   9.079   -9.079   -5.135)   13.829
   3.582   (  11.420  -11.420    0.223)   16.152
   6.434   (  15.834  -15.834  -19.724)   29.841
  17.910   (   5.455   -5.455  -11.314)   13.694
  18.532   (  19.644  -19.644   24.080)   36.765
  27.097   (  -8.508    8.508   -3.272)   12.470
======================= Grid point 32 (23/195) =======================
q-point: (-0.28  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.212   (  12.059  -12.059   -7.090)   18.469
   3.275   (  14.428  -14.428   -1.032)   20.430
   5.749   (  21.500  -21.500  -22.421)   37.779
  17.631   (   5.835   -5.835  -13.554)   15.868
  18.299   (  22.617  -22.617   20.657)   38.076
  27.244   (  -7.366    7.366   -4.094)   11.193
======================= Grid point 33 (24/195) =======================
q-point: (-0.22  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.805   (  14.941  -14.941   -8.547)   22.793
   2.884   (  17.502  -17.502   -2.005)   24.833
   4.918   (  26.578  -26.578  -23.155)   44.147
  17.325   (   6.660   -6.660  -13.621)   16.560
  17.965   (  24.196  -24.196   16.245)   37.878
  27.334   (  -4.940    4.940   -5.470)    8.872
======================= Grid point 34 (25/195) =======================
q-point: (-0.17  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.317   (  17.934  -17.934   -9.062)   26.933
   2.410   (  20.670  -20.670   -2.767)   29.363
   3.977   (  31.368  -31.368  -21.410)   49.258
  17.011   (   7.421   -7.421  -11.346)   15.456
  17.562   (  23.820  -23.820   11.479)   35.588
  27.347   (  -2.180    2.180   -6.162)    6.890
======================= Grid point 35 (26/195) =======================
q-point: (-0.11  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.761   (  21.144  -21.144   -7.416)   30.808
   1.856   (  23.738  -23.738   -2.825)   33.690
   2.967   (  36.130  -36.130  -15.492)   53.393
  16.729   (   7.356   -7.356   -6.699)   12.373
  17.139   (  20.907  -20.907    6.421)   30.257
  27.305   (  -0.022    0.022   -4.748)    4.748
======================= Grid point 36 (27/195) =======================
q-point: (-0.06  0.06  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.179   (  23.968  -23.968    0.000)   33.896
   1.251   (  25.947  -25.947    0.000)   36.695
   1.966   (  40.287  -40.287    0.000)   56.974
  16.530   (   5.916   -5.916    0.000)    8.367
  16.750   (  15.289  -15.289    0.000)   21.622
  27.260   (   1.000   -1.000    0.000)    1.415
======================= Grid point 39 (28/195) =======================
q-point: ( 0.11  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.705   (   0.000   -0.000   32.849)   32.849
   1.705   (   0.000   -0.000   32.849)   32.849
   2.680   (   0.000   -0.000   52.842)   52.842
  17.088   (   0.000   -0.000   28.051)   28.051
  17.088   (   0.000   -0.000   28.051)   28.051
  27.056   (  -0.000    0.000   -6.980)    6.980
======================= Grid point 40 (29/195) =======================
q-point: ( 0.17  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.131   (  -4.770    4.770   28.740)   29.522
   2.148   (  -5.870    5.870   29.642)   30.782
   3.406   ( -10.941   10.941   46.202)   48.724
  17.379   (   4.549   -4.549   32.046)   32.686
  17.460   (  -1.478    1.478   31.661)   31.730
  27.041   (  -5.835    5.835   -8.956)   12.179
======================= Grid point 41 (30/195) =======================
q-point: ( 0.22  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.572   (  -6.401    6.401   23.590)   25.268
   2.639   (  -7.578    7.578   26.474)   28.561
   4.221   ( -16.443   16.443   37.122)   43.804
  17.648   (   4.982   -4.982   33.598)   34.329
  17.874   (  -1.354    1.354   32.497)   32.554
  27.079   (  -6.884    6.884  -11.350)   14.954
======================= Grid point 42 (31/195) =======================
q-point: ( 0.28  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.972   (  -6.517    6.517   18.792)   20.931
   3.103   (  -6.672    6.672   23.860)   25.659
   5.019   ( -18.938   18.938   28.094)   38.815
  17.937   (   4.123   -4.123   34.020)   34.516
  18.261   (   0.298   -0.298   30.816)   30.819
  27.077   (  -5.290    5.290  -13.829)   15.723
======================= Grid point 43 (32/195) =======================
q-point: ( 0.33  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.312   (  -5.903    5.903   14.573)   16.795
   3.500   (  -4.591    4.591   21.385)   22.349
   5.747   ( -19.180   19.180   19.555)   33.439
  18.249   (   3.751   -3.751   35.053)   35.452
  18.564   (   2.963   -2.963   27.150)   27.471
  27.000   (  -3.153    3.153  -15.961)   16.572
======================= Grid point 44 (33/195) =======================
q-point: ( 0.39  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.586   (  -4.785    4.785   10.849)   12.787
   3.813   (  -2.003    2.003   18.750)   18.963
   6.358   ( -17.387   17.387   11.376)   27.093
  18.577   (   4.570   -4.570   37.392)   37.946
  18.749   (   5.949   -5.949   22.038)   23.590
  26.862   (  -1.916    1.916  -17.152)   17.365
======================= Grid point 45 (34/195) =======================
q-point: ( 0.44  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.787   (  -3.145    3.145    7.511)    8.730
   4.030   (   0.837   -0.837   15.989)   16.033
   6.811   ( -13.713   13.713    3.419)   19.693
  18.800   (   8.606   -8.606   15.961)   20.072
  18.901   (   6.881   -6.881   40.575)   41.726
  26.708   (  -2.387    2.387  -17.081)   17.412
======================= Grid point 46 (35/195) =======================
q-point: ( 0.50  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.905   (  -0.937    0.937    4.501)    4.692
   4.146   (   3.823   -3.823   13.354)   14.407
   7.068   (  -8.464    8.464   -4.184)   12.680
  18.724   (  10.439  -10.439    9.434)   17.519
  19.190   (  10.667  -10.667   43.356)   45.905
  26.593   (  -4.612    4.612  -15.868)   17.156
======================= Grid point 47 (36/195) =======================
q-point: (-0.44  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.933   (   1.804   -1.804    1.819)    3.133
   4.164   (   6.883   -6.883   11.126)   14.783
   7.104   (  -2.162    2.162  -10.970)   11.388
  18.542   (  11.190  -11.190    3.098)   16.125
  19.400   (  15.482  -15.482   44.299)   49.415
  26.563   (  -7.841    7.841  -13.965)   17.832
======================= Grid point 48 (37/195) =======================
q-point: (-0.39  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.862   (   4.923   -4.923   -0.494)    6.980
   4.087   (   9.956   -9.956    9.427)   16.945
   6.915   (   4.547   -4.547  -16.296)   17.519
  18.285   (  10.922  -10.922   -2.318)   15.619
  19.489   (  20.510  -20.510   42.397)   51.369
  26.631   ( -10.782   10.782  -11.946)   19.371
======================= Grid point 49 (38/195) =======================
q-point: (-0.33  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.692   (   8.191   -8.191   -2.361)   11.821
   3.921   (  12.989  -12.989    8.126)   20.086
   6.520   (  11.102  -11.102  -19.546)   25.071
  17.988   (  10.040  -10.040   -6.153)   15.475
  19.426   (  24.879  -24.879   37.426)   51.367
  26.775   ( -12.207   12.207  -10.254)   20.080
======================= Grid point 50 (39/195) =======================
q-point: (-0.28  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.427   (  11.394  -11.394   -3.647)   16.521
   3.670   (  15.937  -15.937    6.828)   23.550
   5.953   (  17.216  -17.216  -20.198)   31.634
  17.680   (   9.135   -9.135   -7.967)   15.178
  19.200   (  27.959  -27.959   29.864)   49.550
  26.944   ( -11.535   11.535   -8.918)   18.591
======================= Grid point 51 (40/195) =======================
q-point: (-0.22  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.079   (  14.429  -14.429   -4.106)   20.815
   3.335   (  18.736  -18.736    5.004)   26.966
   5.259   (  22.903  -22.903  -17.684)   36.903
  17.381   (   8.646   -8.646   -7.554)   14.373
  18.823   (  29.350  -29.350   20.517)   46.301
  27.091   (  -9.012    9.012   -7.371)   14.723
======================= Grid point 52 (41/195) =======================
q-point: (-0.17  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.663   (  17.291  -17.291   -3.234)   24.666
   2.910   (  21.250  -21.250    2.363)   30.145
   4.487   (  28.301  -28.301  -11.210)   41.564
  17.105   (   8.551   -8.551   -4.869)   13.036
  18.320   (  28.736  -28.736   10.213)   41.903
  27.190   (  -5.521    5.521   -4.624)    9.074
======================= Grid point 53 (42/195) =======================
q-point: (-0.11  0.11  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 865
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.206   (  19.839  -19.839    0.000)   28.056
   2.402   (  23.191  -23.191    0.000)   32.797
   3.697   (  33.205  -33.205    0.000)   46.959
  16.877   (   8.339   -8.339    0.000)   11.792
  17.736   (  25.930  -25.930    0.000)   36.670
  27.251   (  -2.224    2.224    0.000)    3.146
======================= Grid point 58 (43/195) =======================
q-point: ( 0.17  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 475
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.413   (   0.000   -0.000   27.618)   27.618
   2.413   (   0.000   -0.000   27.618)   27.618
   3.828   (   0.000   -0.000   44.990)   44.990
  17.837   (   0.000   -0.000   35.062)   35.062
  17.837   (   0.000   -0.000   35.062)   35.062
  26.862   (  -0.000    0.000   -9.408)    9.408
======================= Grid point 59 (44/195) =======================
q-point: ( 0.22  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.750   (  -2.586    2.586   24.113)   24.389
   2.786   (  -5.284    5.284   24.846)   25.945
   4.402   (  -6.563    6.563   38.540)   39.642
  18.196   (   4.738   -4.738   36.896)   37.499
  18.258   (  -0.182    0.182   35.727)   35.728
  26.804   (  -4.520    4.520  -11.139)   12.843
======================= Grid point 60 (45/195) =======================
q-point: ( 0.28  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.083   (  -3.612    3.612   19.979)   20.622
   3.216   (  -7.226    7.226   22.640)   24.839
   5.009   ( -10.807   10.807   29.865)   33.549
  18.504   (   6.189   -6.189   38.664)   39.642
  18.659   (   1.121   -1.121   33.886)   33.923
  26.780   (  -5.382    5.382  -13.885)   15.835
======================= Grid point 61 (46/195) =======================
q-point: ( 0.33  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.380   (  -3.740    3.740   15.964)   16.818
   3.632   (  -6.149    6.149   21.136)   22.855
   5.591   ( -13.557   13.557   20.480)   28.054
  18.827   (   5.862   -5.862   41.178)   42.005
  18.980   (   3.682   -3.682   30.100)   30.547
  26.714   (  -4.191    4.191  -16.794)   17.809
======================= Grid point 62 (47/195) =======================
q-point: ( 0.39  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.627   (  -3.307    3.307   12.275)   13.136
   3.985   (  -3.393    3.393   19.847)   20.419
   6.108   ( -14.583   14.583   11.190)   23.464
  19.179   (   7.001   -7.001   24.971)   26.862
  19.197   (   5.632   -5.632   44.996)   45.696
  26.585   (  -2.943    2.943  -19.280)   19.724
======================= Grid point 63 (48/195) =======================
q-point: ( 0.44  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.818   (  -2.385    2.385    8.933)    9.549
   4.249   (  -0.012    0.012   18.478)   18.478
   6.520   ( -13.545   13.545    2.423)   19.308
  19.232   (  10.457  -10.457   19.014)   24.088
  19.608   (   6.819   -6.819   49.648)   50.576
  26.415   (  -3.023    3.023  -21.050)   21.480
======================= Grid point 64 (49/195) =======================
q-point: ( 0.50  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.945   (  -0.923    0.923    5.939)    6.081
   4.416   (   3.467   -3.467   17.095)   17.784
   6.785   ( -10.471   10.471   -5.501)   15.797
  19.137   (  13.401  -13.401   12.719)   22.825
  20.019   (  10.156  -10.156   53.628)   55.517
  26.253   (  -5.066    5.066  -21.979)   23.117
======================= Grid point 65 (50/195) =======================
q-point: (-0.44  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.997   (   1.120   -1.120    3.328)    3.685
   4.488   (   6.860   -6.860   15.838)   18.573
   6.873   (  -5.784    5.784  -12.110)   14.614
  18.913   (  15.258  -15.258    6.662)   22.583
  20.359   (  15.567  -15.567   54.802)   59.058
  26.156   (  -8.814    8.814  -21.758)   25.075
======================= Grid point 67 (51/195) =======================
q-point: (-0.39  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.968   (   3.703   -3.703    1.168)    5.365
   4.467   (  10.123  -10.123   14.619)   20.461
   6.773   (  -0.115    0.115  -16.705)   16.706
  18.595   (  15.685  -15.685    1.523)   22.235
  20.545   (  21.932  -21.932   51.164)   59.831
  26.166   ( -12.961   12.961  -19.797)   26.980
======================= Grid point 68 (52/195) =======================
q-point: (-0.33  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.851   (   6.666   -6.666   -0.440)    9.437
   4.356   (  13.214  -13.214   12.972)   22.748
   6.497   (   5.927   -5.927  -18.442)   20.257
  18.232   (  14.756  -14.756   -2.043)   20.969
  20.514   (  27.477  -27.477   42.083)   57.280
  26.295   ( -15.627   15.627  -15.912)   27.233
======================= Grid point 69 (53/195) =======================
q-point: (-0.28  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.648   (   9.764   -9.764   -1.348)   13.873
   4.150   (  16.038  -16.038   10.104)   24.830
   6.080   (  11.950  -11.950  -16.431)   23.571
  17.872   (  12.990  -12.990   -3.572)   18.714
  20.245   (  30.942  -30.942   28.928)   52.456
  26.509   ( -15.633   15.633  -10.826)   24.617
======================= Grid point 70 (54/195) =======================
q-point: (-0.22  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.368   (  12.732  -12.732   -1.328)   18.055
   3.839   (  18.436  -18.436    5.428)   26.632
   5.572   (  17.821  -17.821  -10.060)   27.136
  17.551   (  11.113  -11.113   -2.877)   15.977
  19.756   (  32.126  -32.126   14.076)   47.564
  26.754   ( -13.190   13.190   -5.389)   19.416
======================= Grid point 71 (55/195) =======================
q-point: (-0.17  0.17  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.029   (  15.340  -15.340    0.000)   21.694
   3.416   (  20.242  -20.242    0.000)   28.626
   5.028   (  23.432  -23.432    0.000)   33.137
  17.291   (   9.636   -9.636    0.000)   13.627
  19.091   (  31.329  -31.329    0.000)   44.306
  26.988   (  -9.399    9.399    0.000)   13.292
======================= Grid point 78 (56/195) =======================
q-point: ( 0.22  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.995   (   0.000   -0.000   22.290)   22.290
   2.995   (   0.000   -0.000   22.290)   22.290
   4.768   (   0.000   -0.000   35.247)   35.247
  18.682   (   0.000   -0.000   36.371)   36.371
  18.682   (   0.000   -0.000   36.371)   36.371
  26.625   (  -0.000    0.000  -10.643)   10.643
======================= Grid point 79 (57/195) =======================
q-point: ( 0.28  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.257   (  -1.372    1.372   19.286)   19.383
   3.308   (  -5.422    5.422   19.884)   21.311
   5.192   (  -3.741    3.741   28.870)   29.351
  19.062   (   3.221   -3.221   36.420)   36.704
  19.089   (   1.071   -1.071   34.585)   34.618
  26.532   (  -2.635    2.635  -11.844)   12.417
======================= Grid point 80 (58/195) =======================
q-point: ( 0.33  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.502   (  -1.897    1.897   15.916)   16.141
   3.693   (  -7.980    7.980   18.238)   21.447
   5.603   (  -6.371    6.371   20.812)   22.679
  19.414   (   3.483   -3.483   38.456)   38.770
  19.423   (   3.375   -3.375   30.907)   31.273
  26.445   (  -2.662    2.662  -14.449)   14.931
======================= Grid point 81 (59/195) =======================
q-point: ( 0.39  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.714   (  -1.882    1.882   12.578)   12.857
   4.085   (  -7.165    7.165   17.571)   20.284
   5.969   (  -8.489    8.489   11.851)   16.869
  19.641   (   6.657   -6.657   25.962)   27.616
  19.818   (   2.361   -2.361   42.704)   42.835
  26.308   (  -1.585    1.585  -17.728)   17.869
======================= Grid point 82 (60/195) =======================
q-point: ( 0.44  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.882   (  -1.462    1.462    9.433)    9.657
   4.423   (  -3.964    3.964   17.652)   18.521
   6.269   (  -9.737    9.737    2.688)   14.030
  19.712   (  10.478  -10.478   20.289)   25.125
  20.302   (   2.003   -2.003   48.549)   48.632
  26.109   (  -1.129    1.129  -21.344)   21.404
======================= Grid point 83 (61/195) =======================
q-point: ( 0.50  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.000   (  -0.626    0.626    6.567)    6.626
   4.675   (   0.299   -0.299   18.054)   18.059
   6.477   (  -9.476    9.476   -5.831)   14.615
  19.624   (  14.247  -14.247   14.357)   24.740
  20.838   (   4.241   -4.241   54.466)   54.795
  25.876   (  -2.642    2.642  -25.198)   25.473
======================= Grid point 84 (62/195) =======================
q-point: ( 0.56  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.062   (   0.700   -0.700    4.062)    4.181
   4.829   (   4.699   -4.699   18.260)   19.431
   6.568   (  -7.451    7.451  -12.697)   16.500
  19.388   (  17.257  -17.257    8.677)   25.901
  21.337   (  10.242  -10.242   57.733)   59.522
  25.668   (  -6.999    6.999  -28.369)   30.046
======================= Grid point 85 (63/195) =======================
q-point: (-0.39  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.060   (   2.564   -2.564    2.019)    4.150
   4.881   (   8.819   -8.819   17.469)   21.464
   6.528   (  -3.914    3.914  -16.653)   17.549
  19.035   (  18.844  -18.844    3.862)   26.928
  21.655   (  19.438  -19.438   53.803)   60.419
  25.574   ( -13.644   13.644  -27.607)   33.682
======================= Grid point 86 (64/195) =======================
q-point: (-0.33  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.988   (   4.922   -4.922    0.559)    6.983
   4.820   (  12.441  -12.441   14.405)   22.739
   6.371   (   0.709   -0.709  -16.147)   16.178
  18.618   (  18.643  -18.643    0.542)   26.371
  21.648   (  28.018  -28.018   39.374)   55.860
  25.679   ( -18.999   18.999  -19.857)   33.411
======================= Grid point 87 (65/195) =======================
q-point: (-0.28  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.843   (   7.605   -7.605   -0.175)   10.757
   4.629   (  15.306  -15.306    8.056)   23.096
   6.134   (   6.118   -6.118  -10.083)   13.286
  18.197   (  16.811  -16.811   -0.808)   23.788
  21.276   (  32.426  -32.426   19.053)   49.657
  25.975   ( -19.840   19.840   -8.859)   29.424
======================= Grid point 88 (66/195) =======================
q-point: (-0.22  0.22  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 865
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.631   (  10.327  -10.327    0.000)   14.604
   4.293   (  17.262  -17.262    0.000)   24.413
   5.861   (  12.014  -12.014    0.000)   16.990
  17.829   (  14.019  -14.019    0.000)   19.826
  20.611   (  33.282  -33.282    0.000)   47.068
  26.363   ( -17.200   17.200    0.000)   24.324
======================= Grid point 97 (67/195) =======================
q-point: ( 0.28  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 475
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.457   (   0.000   -0.000   17.316)   17.316
   3.457   (   0.000   -0.000   17.316)   17.316
   5.470   (   0.000   -0.000   24.870)   24.870
  19.500   (   0.000   -0.000   32.984)   32.984
  19.500   (   0.000   -0.000   32.984)   32.984
  26.376   (  -0.000    0.000  -10.499)   10.499
======================= Grid point 98 (68/195) =======================
q-point: ( 0.33  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.653   (  -0.651    0.651   14.697)   14.725
   3.717   (  -5.832    5.832   15.198)   17.292
   5.750   (  -1.696    1.696   19.019)   19.170
  19.844   (   2.113   -2.113   29.559)   29.710
  19.867   (   0.224   -0.224   31.883)   31.885
  26.261   (  -0.459    0.459  -11.104)   11.123
======================= Grid point 99 (69/195) =======================
q-point: ( 0.39  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.827   (  -0.815    0.815   11.851)   11.907
   4.069   (  -9.227    9.227   13.902)   19.067
   5.984   (  -2.684    2.684   12.070)   12.653
  20.080   (   5.175   -5.175   24.923)   25.975
  20.269   (  -2.248    2.248   34.090)   34.238
  26.116   (   0.904   -0.904  -13.482)   13.543
======================= Grid point 100 (70/195) =======================
q-point: ( 0.44  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.966   (  -0.656    0.656    9.022)    9.069
   4.452   (  -9.282    9.282   13.857)   19.087
   6.154   (  -3.542    3.542    4.323)    6.617
  20.176   (   8.971   -8.971   19.622)   23.366
  20.783   (  -5.210    5.210   39.064)   39.753
  25.896   (   2.644   -2.644  -17.329)   17.727
======================= Grid point 101 (71/195) =======================
q-point: ( 0.50  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.067   (  -0.224    0.224    6.360)    6.368
   4.809   (  -6.372    6.372   15.296)   17.753
   6.248   (  -4.346    4.346   -3.967)    7.315
  20.116   (  13.110  -13.110   14.110)   23.298
  21.416   (  -6.123    6.123   45.629)   46.443
  25.590   (   3.218   -3.218  -22.911)   23.359
======================= Grid point 102 (72/195) =======================
q-point: ( 0.56  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.124   (   0.529   -0.529    3.993)    4.063
   5.094   (  -1.234    1.234   17.778)   17.863
   6.264   (  -5.055    5.055  -11.958)   13.931
  19.897   (  16.966  -16.966    8.842)   25.571
  22.099   (  -2.206    2.206   52.306)   52.399
  25.226   (   0.499   -0.499  -30.581)   30.590
======================= Grid point 103 (73/195) =======================
q-point: (-0.39  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.134   (   1.688   -1.688    2.052)    3.149
   5.265   (   5.344   -5.344   18.732)   20.199
   6.216   (  -5.584    5.584  -16.752)   18.520
  19.541   (  19.753  -19.753    4.360)   28.273
  22.653   (  10.577  -10.577   53.483)   55.535
  24.916   (  -9.540    9.540  -35.992)   38.438
======================= Grid point 104 (74/195) =======================
q-point: (-0.33  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.091   (   3.308   -3.308    0.674)    4.726
   5.257   (  11.253  -11.253   12.174)   20.036
   6.159   (  -4.479    4.479  -12.163)   13.714
  19.095   (  20.759  -20.759    1.253)   29.385
  22.713   (  27.509  -27.509   30.646)   49.525
  24.952   ( -22.701   22.701  -20.753)   38.228
======================= Grid point 105 (75/195) =======================
q-point: (-0.28  0.28  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.994   (   5.343   -5.343    0.000)    7.556
   5.026   (  14.059  -14.059    0.000)   19.883
   6.143   (   0.112   -0.112    0.000)    0.158
  18.624   (  19.676  -19.676    0.000)   27.826
  22.161   (  32.538  -32.538    0.000)   46.016
  25.403   ( -23.654   23.654    0.000)   33.451
======================= Grid point 116 (76/195) =======================
q-point: ( 0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.808   (   0.000   -0.000   12.722)   12.722
   3.808   (   0.000   -0.000   12.722)   12.722
   5.937   (   0.000   -0.000   15.326)   15.326
  20.199   (   0.000   -0.000   26.450)   26.450
  20.199   (   0.000   -0.000   26.450)   26.450
  26.145   (  -0.000    0.000   -9.035)    9.035
======================= Grid point 117 (77/195) =======================
q-point: ( 0.39  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.946   (  -0.207    0.207   10.356)   10.360
   4.021   (  -6.354    6.354   10.777)   14.031
   6.092   (  -0.109    0.109   10.610)   10.612
  20.452   (   2.875   -2.875   22.157)   22.527
  20.531   (  -3.750    3.750   24.582)   25.147
  26.023   (   1.875   -1.875   -9.078)    9.457
======================= Grid point 118 (78/195) =======================
q-point: ( 0.44  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.057   (  -0.139    0.139    7.859)    7.862
   4.343   ( -10.719   10.719    9.589)   17.938
   6.183   (   0.429   -0.429    5.381)    5.415
  20.575   (   6.437   -6.437   17.227)   19.484
  20.979   (  -9.688    9.688   26.090)   29.468
  25.827   (   5.174   -5.174  -10.998)   13.210
======================= Grid point 119 (79/195) =======================
q-point: ( 0.50  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.135   (   0.100   -0.100    5.420)    5.422
   4.728   ( -12.259   12.259    9.559)   19.798
   6.194   (   1.262   -1.262   -0.234)    1.799
  20.549   (  10.524  -10.524   12.102)   19.182
  21.598   ( -15.277   15.277   29.508)   36.572
  25.515   (   8.859   -8.859  -14.660)   19.285
======================= Grid point 120 (80/195) =======================
q-point: ( 0.56  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.178   (   0.494   -0.494    3.195)    3.271
   5.125   ( -11.140   11.140   11.161)   19.307
   6.116   (   2.178   -2.178   -6.358)    7.065
  20.365   (  14.737  -14.737    7.192)   22.047
  22.366   ( -18.453   18.453   33.514)   42.476
  25.065   (  12.099  -12.099  -20.594)   26.775
======================= Grid point 121 (81/195) =======================
q-point: (-0.39  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.186   (   1.103   -1.103    1.338)    2.055
   5.495   (  -6.588    6.588   15.611)   18.180
   5.947   (   1.717   -1.717  -13.969)   14.179
  20.035   (  18.384  -18.384    2.976)   26.168
  23.212   ( -16.495   16.495   39.374)   45.766
  24.470   (  12.512  -12.512  -31.372)   36.018
======================= Grid point 123 (82/195) =======================
q-point: (-0.33  0.33  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 865
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.158   (   2.020   -2.020    0.000)    2.857
   5.605   (  10.668  -10.668    0.000)   15.087
   5.898   ( -10.088   10.088    0.000)   14.267
  19.592   (  20.617  -20.617    0.000)   29.157
  23.602   (  28.550  -28.550    0.000)   40.375
  24.177   ( -28.448   28.448    0.000)   40.231
======================= Grid point 136 (83/195) =======================
q-point: ( 0.39  0.39  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 475
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.054   (   0.000   -0.000    8.380)    8.380
   4.054   (   0.000   -0.000    8.380)    8.380
   6.203   (   0.000   -0.000    7.892)    7.892
  20.723   (   0.000   -0.000   18.196)   18.196
  20.723   (   0.000   -0.000   18.196)   18.196
  25.962   (  -0.000    0.000   -6.541)    6.541
======================= Grid point 137 (84/195) =======================
q-point: ( 0.44  0.39  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.139   (   0.056   -0.056    6.168)    6.168
   4.222   (  -6.871    6.871    6.474)   11.676
   6.266   (   1.070   -1.070    4.708)    4.945
  20.871   (   3.355   -3.355   13.579)   14.384
  21.002   (  -7.696    7.696   15.489)   18.931
  25.846   (   4.071   -4.071   -5.982)    8.302
======================= Grid point 138 (85/195) =======================
q-point: ( 0.50  0.39  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.194   (   0.238   -0.238    3.919)    3.933
   4.515   ( -12.086   12.086    4.994)   17.807
   6.259   (   2.731   -2.731    1.591)    4.177
  20.879   (   7.169   -7.169    8.746)   13.389
  21.460   ( -16.553   16.553   14.481)   27.527
  25.618   (   9.331   -9.331   -6.496)   14.708
======================= Grid point 139 (86/195) =======================
q-point: ( 0.56  0.39  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.219   (   0.464   -0.464    1.807)    1.923
   4.886   ( -14.958   14.958    3.605)   21.459
   6.174   (   4.854   -4.854   -0.813)    6.913
  20.738   (  11.296  -11.296    4.077)   16.487
  22.097   ( -24.178   24.178   11.727)   36.148
  25.254   (  15.270  -15.270   -6.447)   22.537
======================= Grid point 140 (87/195) =======================
q-point: (-0.39 -0.61  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 864
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.216   (   0.730   -0.730    0.000)    1.032
   5.273   ( -15.568   15.568    0.000)   22.016
   6.022   (   7.184   -7.184    0.000)   10.160
  20.450   (  15.286  -15.286    0.000)   21.618
  22.810   ( -28.882   28.882    0.000)   40.845
  24.775   (  21.162  -21.162    0.000)   29.928
======================= Grid point 157 (88/195) =======================
q-point: ( 0.44  0.44  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 480
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.200   (   0.000   -0.000    4.166)    4.166
   4.200   (   0.000   -0.000    4.166)    4.166
   6.326   (   0.000   -0.000    3.046)    3.046
  21.044   (   0.000   -0.000    9.208)    9.208
  21.044   (   0.000   -0.000    9.208)    9.208
  25.844   (  -0.000    0.000   -3.410)    3.410
======================= Grid point 158 (89/195) =======================
q-point: ( 0.50  0.44  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.235   (   0.180   -0.180    2.050)    2.065
   4.323   (  -7.212    7.212    2.166)   10.426
   6.331   (   1.729   -1.729    1.231)    2.738
  21.083   (   3.582   -3.582    4.559)    6.816
  21.246   ( -10.229   10.229    5.296)   15.405
  25.750   (   5.491   -5.491   -2.098)    8.044
======================= Grid point 159 (90/195) =======================
q-point: (-0.44 -0.56  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 865
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.240   (   0.360   -0.360    0.000)    0.509
   4.573   ( -12.673   12.673    0.000)   17.922
   6.275   (   3.623   -3.623    0.000)    5.124
  20.981   (   7.407   -7.407    0.000)   10.475
  21.633   ( -19.499   19.499    0.000)   27.575
  25.540   (  11.221  -11.221    0.000)   15.869
======================= Grid point 179 (91/195) =======================
q-point: (-0.50 -0.50  0.00)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 275
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.249   (   0.000   -0.000    0.000)    0.000
   4.249   (   0.000   -0.000    0.000)    0.000
   6.360   (   0.000   -0.000    0.000)    0.000
  21.152   (   0.000   -0.000    0.000)    0.000
  21.152   (   0.000   -0.000    0.000)    0.000
  25.804   (   0.000   -0.000    0.000)    0.000
======================= Grid point 378 (92/195) =======================
q-point: ( 0.17  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   1.830   (   0.000   11.150   29.161)   31.220
   1.917   (   0.000   16.709   29.124)   33.577
   2.997   (   0.000   24.534   45.324)   51.538
  16.908   (   0.000   -9.413   24.587)   26.327
  17.258   (   0.000   14.309   27.454)   30.959
  27.183   (   0.000    7.110   -6.703)    9.772
======================= Grid point 379 (93/195) =======================
q-point: ( 0.22  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.314   (  -4.779   11.790   24.571)   27.669
   2.446   (  -0.527   17.560   25.585)   31.036
   3.849   ( -11.652   24.922   37.785)   46.740
  17.121   (  -5.221   -9.136   24.922)   27.052
  17.750   (   8.332   19.392   30.542)   37.125
  27.198   (  -4.856    4.673   -8.446)   10.805
======================= Grid point 380 (94/195) =======================
q-point: ( 0.28  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.773   (  -5.831   11.158   19.947)   23.588
   2.925   (   0.374   15.713   22.727)   27.632
   4.706   ( -17.903   23.491   29.958)   42.069
  17.388   (  -8.872   -8.659   24.048)   27.055
  18.214   (  15.939   21.997   31.161)   41.339
  27.189   (  -6.679    0.863  -10.815)   12.740
======================= Grid point 381 (95/195) =======================
q-point: ( 0.33  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.170   (  -5.103   10.194   15.555)   19.285
   3.332   (   1.324   12.556   20.278)   23.887
   5.499   ( -20.443   20.928   22.477)   36.893
  17.680   ( -10.562   -8.500   23.269)   26.931
  18.605   (  22.925   22.298   29.787)   43.704
  27.115   (  -6.712   -2.757  -13.174)   15.041
======================= Grid point 382 (96/195) =======================
q-point: ( 0.39  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.491   (  -3.666    8.774   11.531)   14.947
   3.659   (   2.391    9.032   17.755)   20.063
   6.183   ( -20.174   17.244   15.168)   30.568
  17.974   ( -10.185   -8.493   23.008)   26.556
  18.887   (  29.016   20.739   27.249)   44.884
  26.981   (  -6.412   -5.032  -14.655)   16.769
======================= Grid point 383 (97/195) =======================
q-point: ( 0.44  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.728   (  -1.993    6.586    7.951)   10.515
   3.899   (   3.941    5.684   14.912)   16.438
   6.714   ( -17.546   12.392    7.796)   22.852
  18.248   (  -7.537   -8.172   22.373)   24.983
  19.048   (  33.386   17.355   24.793)   45.061
  26.822   (  -6.838   -5.229  -14.660)   17.001
======================= Grid point 384 (98/195) =======================
q-point: ( 0.50  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.874   (  -0.108    3.569    4.814)    5.994
   4.046   (   6.118    2.699   11.850)   13.606
   7.053   ( -12.908    6.507    0.376)   14.461
  18.464   (  -1.776   -6.128   18.628)   19.690
  19.103   (  34.455   11.074   24.826)   43.888
  26.695   (  -8.397   -3.296  -13.175)   15.967
======================= Grid point 385 (99/195) =======================
q-point: (-0.44  0.11  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1467
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.924   (   2.130    0.000    2.130)    3.012
   4.096   (   8.828   -0.000    8.828)   12.485
   7.169   (  -6.727    0.000   -6.727)    9.514
  18.554   (   8.671    0.000    8.671)   12.262
  19.106   (  29.993    0.000   29.993)   42.417
  26.646   ( -10.782    0.000  -10.782)   15.247
======================= Grid point 398 (100/195) =======================
q-point: ( 0.22  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.484   (   0.000    6.221   26.283)   27.010
   2.567   (   0.000   12.490   26.027)   28.869
   4.014   (   0.000   15.097   40.852)   43.552
  17.606   (   0.000  -16.214   33.964)   37.636
  17.993   (   0.000   13.154   34.495)   36.918
  26.975   (  -0.000    7.429  -10.656)   12.989
======================= Grid point 399 (101/195) =======================
q-point: ( 0.28  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.866   (  -2.400    7.340   22.254)   23.556
   3.020   (  -1.066   14.038   23.188)   27.127
   4.683   (  -8.293   16.473   32.793)   37.623
  17.828   (  -4.541  -18.494   34.296)   39.229
  18.526   (  10.099   17.952   35.034)   40.641
  26.940   (  -4.356    4.945  -13.060)   14.628
======================= Grid point 400 (102/195) =======================
q-point: ( 0.33  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.218   (  -2.406    7.659   17.869)   19.590
   3.445   (  -0.308   12.274   21.413)   24.683
   5.342   ( -14.035   16.627   23.925)   32.339
  18.083   (  -7.934  -18.920   34.160)   39.848
  18.983   (  19.746   20.654   33.976)   44.394
  26.872   (  -6.387    1.010  -15.726)   17.004
======================= Grid point 401 (103/195) =======================
q-point: ( 0.39  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.514   (  -1.605    7.259   13.750)   15.631
   3.805   (   1.406    8.906   19.840)   21.793
   5.943   ( -17.002   15.568   15.135)   27.577
  18.364   (  -9.473  -18.763   33.860)   39.854
  19.334   (  28.321   21.115   32.484)   47.991
  26.746   (  -7.245   -2.242  -17.963)   19.498
======================= Grid point 402 (104/195) =======================
q-point: ( 0.44  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.744   (  -0.560    6.041   10.076)   11.761
   4.080   (   3.568    5.136   17.996)   19.051
   6.441   ( -17.140   13.036    6.718)   22.558
  18.640   (  -7.757  -17.943   31.694)   37.237
  19.578   (  34.358   18.700   32.704)   50.987
  26.581   (  -8.165   -3.461  -19.225)   21.172
======================= Grid point 403 (105/195) =======================
q-point: ( 0.50  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.902   (   0.716    4.018    6.864)    7.986
   4.261   (   6.002    1.479   15.932)   17.089
   6.791   ( -14.683    9.025   -1.198)   17.276
  18.833   (  -0.107  -14.376   24.311)   28.244
  19.753   (  34.745   11.212   37.422)   52.281
  26.424   (  -9.870   -2.084  -19.317)   21.793
======================= Grid point 404 (106/195) =======================
q-point: ( 0.56  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.980   (   2.399    1.342    4.122)    4.955
   4.347   (   8.637   -1.958   13.906)   16.487
   6.957   ( -10.100    3.874   -8.306)   13.638
  18.851   (  12.148   -7.717   12.932)   19.349
  19.914   (  30.092   -1.176   44.399)   53.648
  26.328   ( -12.393    1.501  -18.334)   22.180
======================= Grid point 405 (107/195) =======================
q-point: (-0.39  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.972   (   4.651   -1.684    1.856)    5.283
   4.339   (  11.371   -5.286   12.101)   17.426
   6.921   (  -4.024   -1.826  -14.042)   14.721
  18.688   (  21.120   -2.356    3.570)   21.549
  20.036   (  27.473  -12.729   46.838)   55.772
  26.328   ( -15.119    5.903  -16.479)   23.130
======================= Grid point 406 (108/195) =======================
q-point: (-0.33  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.875   (   7.515   -4.666    0.043)    8.846
   4.241   (  14.041   -8.680   10.489)   19.558
   6.686   (   2.896   -7.418  -17.702)   19.411
  18.412   (  24.794    0.873   -2.374)   24.923
  20.045   (  28.350  -20.888   43.328)   55.833
  26.429   ( -17.071    9.224  -14.024)   23.941
======================= Grid point 407 (109/195) =======================
q-point: (-0.28  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.689   (  10.888   -7.284   -1.330)   13.167
   4.054   (  16.470  -12.189    8.771)   22.288
   6.278   (  10.184  -12.394  -18.592)   24.555
  18.094   (  25.347    3.217   -5.407)   26.116
  19.888   (  30.587  -25.728   35.263)   53.301
  26.600   ( -17.450   10.037  -11.269)   23.070
======================= Grid point 408 (110/195) =======================
q-point: (-0.22  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.422   (  14.545   -9.420   -2.120)   17.458
   3.776   (  18.594  -15.553    6.413)   25.075
   5.740   (  17.629  -16.548  -16.078)   29.036
  17.782   (  24.899    5.241   -5.851)   26.108
  19.554   (  32.610  -27.620   24.300)   49.161
  26.791   ( -16.090    8.046   -8.271)   19.800
======================= Grid point 409 (111/195) =======================
q-point: (-0.17  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.085   (  18.217  -11.087   -1.921)   21.412
   3.405   (  20.595  -18.157    3.116)   27.632
   5.121   (  25.122  -19.995   -9.754)   33.557
  17.505   (  24.708    6.986   -3.948)   25.978
  19.057   (  33.593  -26.962   11.989)   44.712
  26.960   ( -13.522    3.974   -4.651)   14.842
======================= Grid point 410 (112/195) =======================
q-point: (-0.11  0.17  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   2.702   (  21.671  -12.138   -0.000)   24.839
   2.949   (  22.902  -19.106   -0.000)   29.825
   4.474   (  32.218  -22.933   -0.000)   39.547
  17.280   (  24.988    8.409   -0.000)   26.365
  18.435   (  33.047  -24.127   -0.000)   40.917
  27.091   ( -10.705   -0.825    0.000)   10.736
======================= Grid point 417 (113/195) =======================
q-point: ( 0.28  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.047   (   0.000    4.346   21.831)   22.259
   3.123   (   0.000   10.493   21.539)   23.959
   4.877   (   0.000    9.019   32.580)   33.805
  18.445   (   0.000  -18.078   36.889)   41.080
  18.827   (   0.000   12.255   36.023)   38.051
  26.701   (  -0.000    5.130  -12.490)   13.503
======================= Grid point 418 (114/195) =======================
q-point: ( 0.33  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.343   (   0.026    5.505   18.412)   19.217
   3.519   (  -3.098   12.251   19.393)   23.147
   5.357   (  -5.534   10.017   24.633)   27.162
  18.673   (  -4.765  -21.258   37.052)   42.982
  19.350   (  10.982   17.035   34.779)   40.254
  26.609   (  -3.265    2.447  -14.912)   15.460
======================= Grid point 419 (115/195) =======================
q-point: ( 0.39  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.600   (   0.600    5.885   14.681)   15.828
   3.913   (  -2.750   10.645   18.417)   21.449
   5.807   ( -10.123   10.548   15.700)   21.453
  18.927   (  -7.764  -22.020   36.838)   43.614
  19.779   (  20.829   19.809   33.912)   44.455
  26.476   (  -5.165   -1.087  -17.759)   18.527
======================= Grid point 420 (116/195) =======================
q-point: ( 0.44  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.807   (   1.258    5.451   11.177)   12.499
   4.248   (  -0.233    6.990   17.914)   19.231
   6.197   ( -13.035   10.364    6.598)   17.913
  19.179   (  -6.154  -20.951   33.873)   40.302
  20.131   (  26.887   19.012   36.214)   48.947
  26.293   (  -6.673   -3.331  -20.507)   21.821
======================= Grid point 421 (117/195) =======================
q-point: ( 0.50  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.957   (   2.042    4.291    8.029)    9.330
   4.494   (   3.047    2.563   17.456)   17.904
   6.495   ( -13.542    8.901   -1.952)   16.323
  19.332   (   3.049  -16.897   25.602)   30.826
  20.465   (  26.345   12.080   43.586)   52.342
  26.089   (  -8.733   -2.888  -22.905)   24.683
======================= Grid point 422 (118/195) =======================
q-point: ( 0.56  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.045   (   3.133    2.506    5.290)    6.639
   4.644   (   6.340   -1.881   16.879)   18.129
   6.665   ( -11.514    5.996   -9.289)   15.963
  19.296   (  15.251  -11.650   15.315)   24.554
  20.808   (  23.530    0.404   51.183)   56.334
  25.915   ( -11.843    0.782  -24.563)   27.280
======================= Grid point 423 (119/195) =======================
q-point: (-0.39  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.064   (   4.734    0.234    3.005)    5.612
   4.697   (   9.401   -6.116   16.020)   19.555
   6.683   (  -7.377    2.004  -14.661)   16.534
  19.078   (  23.767   -7.992    7.027)   26.040
  21.079   (  24.189  -12.016   52.814)   59.320
  25.831   ( -15.740    6.867  -24.269)   29.730
======================= Grid point 424 (120/195) =======================
q-point: (-0.33  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.008   (   6.945   -2.337    1.218)    7.428
   4.655   (  12.149  -10.139   14.402)   21.397
   6.546   (  -1.761   -2.483  -17.102)   17.370
  18.743   (  27.636   -5.360    1.406)   28.186
  21.159   (  27.601  -22.548   46.153)   58.312
  25.886   ( -19.081   12.678  -20.497)   30.740
======================= Grid point 425 (121/195) =======================
q-point: (-0.28  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.875   (   9.694   -4.988    0.001)   10.902
   4.509   (  14.510  -13.862   11.209)   22.985
   6.278   (   4.784   -6.973  -15.539)   17.691
  18.364   (  28.175   -2.605   -1.572)   28.339
  20.975   (  31.048  -28.993   32.513)   53.495
  26.084   ( -20.278   15.148  -13.716)   28.789
======================= Grid point 426 (122/195) =======================
q-point: (-0.22  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.667   (  12.727   -7.499   -0.497)   14.781
   4.250   (  16.488  -16.877    5.947)   24.333
   5.924   (  11.833  -11.327   -9.419)   18.896
  18.006   (  27.021    0.616   -1.972)   27.100
  20.527   (  33.226  -31.055   15.872)   48.169
  26.370   ( -19.089   13.553   -6.328)   24.251
======================= Grid point 427 (123/195) =======================
q-point: (-0.17  0.22  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.396   (  15.690   -9.643   -0.000)   18.416
   3.872   (  18.386  -18.498    0.000)   26.081
   5.530   (  18.854  -15.695   -0.000)   24.531
  17.712   (  25.607    3.989   -0.000)   25.916
  19.864   (  34.085  -29.991   -0.000)   45.401
  26.678   ( -16.468    9.251    0.000)   18.889
======================= Grid point 436 (124/195) =======================
q-point: ( 0.33  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.504   (   0.000    3.768   17.267)   17.673
   3.573   (   0.000    9.509   16.946)   19.431
   5.526   (   0.000    4.708   22.935)   23.413
  19.293   (   0.000  -16.396   34.948)   38.603
  19.640   (   0.000   11.859   32.918)   34.989
  26.405   (  -0.000    1.710  -12.632)   12.748
======================= Grid point 437 (125/195) =======================
q-point: ( 0.39  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.725   (   2.230    4.901   14.337)   15.315
   3.924   (  -5.585   11.533   15.304)   19.961
   5.827   (  -3.208    4.829   15.654)   16.693
  19.519   (  -4.273  -18.708   34.549)   39.521
  20.124   (   9.285   16.724   31.082)   36.496
  26.258   (  -1.682   -1.546  -14.827)   15.002
======================= Grid point 438 (126/195) =======================
q-point: ( 0.44  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.902   (   3.006    5.074   11.194)   12.653
   4.303   (  -5.990   10.475   14.982)   19.236
   6.075   (  -6.163    4.949    7.478)   10.881
  19.743   (  -3.111  -18.131   31.782)   36.722
  20.563   (  13.851   18.611   33.107)   40.426
  26.055   (  -2.770   -4.858  -18.046)   18.893
======================= Grid point 439 (127/195) =======================
q-point: ( 0.50  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.034   (   3.546    4.534    8.238)   10.050
   4.641   (  -3.324    6.787   15.777)   17.493
   6.258   (  -8.470    5.078   -1.107)    9.938
  19.869   (   6.009  -14.937   24.390)   29.225
  21.037   (  11.930   15.047   40.798)   45.091
  25.793   (  -4.177   -6.223  -22.261)   23.489
======================= Grid point 440 (128/195) =======================
q-point: ( 0.56  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.117   (   4.158    3.434    5.586)    7.764
   4.898   (   0.893    1.511   17.185)   17.275
   6.361   (  -9.445    4.817   -9.154)   14.008
  19.811   (  17.203  -11.822   15.640)   26.083
  21.573   (  10.225    6.629   49.660)   51.133
  25.503   (  -7.070   -3.857  -27.467)   28.623
======================= Grid point 441 (129/195) =======================
q-point: (-0.39  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.146   (   5.061    1.841    3.332)    6.333
   5.051   (   5.477   -4.348   17.943)   19.258
   6.376   (  -8.600    3.838  -14.960)   17.677
  19.569   (  25.337  -10.273    8.358)   28.590
  22.053   (  14.488   -6.371   53.662)   55.947
  25.263   ( -12.808    3.910  -31.471)   34.201
======================= Grid point 443 (130/195) =======================
q-point: (-0.33  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.118   (   6.425   -0.185    1.572)    6.617
   5.080   (   9.567   -9.952   15.829)   21.002
   6.310   (  -5.738    1.953  -15.966)   17.078
  19.195   (  29.623   -9.162    3.178)   31.170
  22.256   (  23.462  -22.053   44.531)   54.953
  25.216   ( -20.268   15.096  -26.693)   36.759
======================= Grid point 444 (131/195) =======================
q-point: (-0.28  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.027   (   8.287   -2.549    0.416)    8.680
   4.956   (  12.530  -14.271    9.007)   21.019
   6.197   (  -0.833   -1.241  -10.162)   10.272
  18.762   (  30.601   -7.181    0.355)   31.434
  22.012   (  29.982  -31.460   22.080)   48.747
  25.464   ( -23.489   20.494  -12.190)   33.472
======================= Grid point 445 (132/195) =======================
q-point: (-0.22  0.28  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.875   (  10.497   -5.077   -0.000)   11.661
   4.667   (  14.549  -16.523    0.000)   22.016
   6.069   (   5.442   -5.927    0.000)    8.047
  18.345   (  29.301   -3.813   -0.000)   29.548
  21.387   (  32.475  -32.983    0.000)   46.287
  25.900   ( -22.011   18.461    0.000)   28.728
======================= Grid point 456 (133/195) =======================
q-point: ( 0.39  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.856   (   0.000    3.855   12.897)   13.461
   3.917   (   0.000    9.010   12.535)   15.437
   5.953   (   0.000    1.390   13.826)   13.896
  20.051   (   0.000  -11.968   29.449)   31.788
  20.341   (   0.000   11.961   26.602)   29.167
  26.123   (  -0.000   -2.268  -11.296)   11.522
======================= Grid point 457 (134/195) =======================
q-point: ( 0.44  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.013   (   4.134    4.949   10.289)   12.142
   4.234   (  -8.167   11.519   11.222)   18.036
   6.099   (  -1.160    0.481    7.892)    7.991
  20.252   (  -0.276  -12.292   27.464)   30.091
  20.787   (   2.455   16.810   25.410)   30.566
  25.928   (   0.487   -6.508  -13.157)   14.687
======================= Grid point 458 (135/195) =======================
q-point: ( 0.50  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.122   (   4.882    4.876    7.576)   10.247
   4.611   (  -9.834   11.585   11.289)   18.931
   6.173   (  -2.067   -0.370    1.336)    2.489
  20.370   (   7.733  -10.344   21.392)   24.988
  21.312   (  -2.566   18.275   30.232)   35.419
  25.645   (   1.448  -10.364  -16.697)   19.706
======================= Grid point 459 (136/195) =======================
q-point: ( 0.56  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.190   (   5.247    4.226    5.057)    8.425
   4.982   (  -8.162    8.810   13.108)   17.778
   6.168   (  -2.867   -0.776   -5.917)    6.621
  20.314   (  17.598   -9.092   13.920)   24.210
  21.977   (  -7.362   16.561   37.852)   41.967
  25.262   (   2.014  -12.164  -22.637)   25.777
======================= Grid point 460 (137/195) =======================
q-point: ( 0.61  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.216   (   5.598    3.137    2.866)    7.028
   5.297   (  -3.246    2.835   16.769)   17.314
   6.089   (  -4.327    0.257  -13.631)   14.304
  20.078   (  25.162   -9.468    7.463)   27.901
  22.712   (  -5.335    9.331   46.188)   47.422
  24.800   (  -0.560   -8.621  -32.410)   33.541
======================= Grid point 461 (138/195) =======================
q-point: (-0.33  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.199   (   6.146    1.629    1.136)    6.460
   5.465   (   5.910   -7.860   16.091)   18.858
   6.002   (  -8.144    4.951  -15.397)   18.108
  19.697   (  29.769  -10.166    2.688)   31.572
  23.230   (  14.266  -16.305   43.798)   48.864
  24.456   ( -17.444   13.034  -35.330)   41.502
======================= Grid point 462 (139/195) =======================
q-point: (-0.28  0.33  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.136   (   7.041   -0.287   -0.000)    7.046
   5.319   (  10.878  -13.734    0.000)   17.520
   6.068   (  -6.511    4.415    0.000)    7.866
  19.232   (  31.438   -9.700   -0.000)   32.901
  22.878   (  28.306  -32.767    0.000)   43.300
  24.809   ( -26.878   25.451    0.000)   37.016
======================= Grid point 475 (140/195) =======================
q-point: ( 0.44  0.39  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.108   (   0.000    4.316    8.657)    9.673
   4.160   (   0.000    8.756    8.287)   12.056
   6.189   (   0.000   -1.204    6.780)    6.887
  20.649   (   0.000   -6.100   21.353)   22.207
  20.869   (   0.000   12.355   18.427)   22.185
  25.887   (  -0.000   -6.396   -8.636)   10.747
======================= Grid point 476 (141/195) =======================
q-point: ( 0.50  0.39  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.206   (   5.639    5.292    6.225)    9.927
   4.448   ( -10.540   11.979    7.013)   17.429
   6.217   (   0.576   -3.012    2.783)    4.141
  20.776   (   7.491   -4.906   17.024)   19.236
  21.297   (  -9.039   17.616   17.783)   26.613
  25.658   (   3.164  -11.861   -9.612)   15.592
======================= Grid point 477 (142/195) =======================
q-point: ( 0.56  0.39  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.255   (   6.141    4.962    3.834)    8.776
   4.822   ( -13.709   13.615    6.534)   20.396
   6.167   (   1.796   -5.029   -1.136)    5.459
  20.740   (  16.386   -4.784   10.434)   20.006
  21.915   ( -18.782   22.083   19.788)   35.099
  25.304   (   7.504  -17.226  -11.570)   22.066
======================= Grid point 478 (143/195) =======================
q-point: ( 0.61  0.39  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.269   (   6.183    4.188    1.707)    7.660
   5.225   ( -14.577   13.370    6.473)   20.812
   6.029   (   3.832   -7.319   -4.398)    9.359
  20.529   (  23.091   -6.575    4.533)   24.433
  22.688   ( -25.286   25.694   19.076)   40.786
  24.806   (  13.234  -22.414  -13.293)   29.227
======================= Grid point 479 (144/195) =======================
q-point: (-0.33  0.39  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.249   (   6.113    3.079   -0.000)    6.845
   5.606   ( -14.542   12.532    0.000)   19.197
   5.814   (   7.253  -10.542    0.000)   12.796
  20.165   (  27.597   -8.892   -0.000)   28.994
  23.491   ( -29.939   29.152    0.000)   41.787
  24.187   (  21.081  -28.717    0.000)   35.624
======================= Grid point 494 (145/195) =======================
q-point: ( 0.50  0.44  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.260   (   0.000    4.858    4.382)    6.542
   4.305   (   0.000    8.635    4.127)    9.571
   6.291   (  -0.000   -2.973    2.426)    3.837
  21.033   (   0.000   -0.902   11.289)   11.325
  21.195   (   0.000   12.755    9.370)   15.827
  25.729   (  -0.000   -9.824   -4.752)   10.913
======================= Grid point 495 (146/195) =======================
q-point: ( 0.56  0.44  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.303   (   6.497    5.584    2.088)    8.818
   4.559   ( -12.069   12.493    2.422)   17.539
   6.251   (   1.667   -5.113    0.507)    5.402
  21.041   (  13.226   -0.470    5.649)   14.389
  21.590   ( -17.707   19.216    6.632)   26.958
  25.498   (   5.355  -15.831   -3.655)   17.107
======================= Grid point 496 (147/195) =======================
q-point: (-0.39 -0.56  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.300   (   6.592    5.031   -0.000)    8.293
   4.902   ( -15.590   14.838    0.000)   21.522
   6.149   (   3.636   -7.151    0.000)    8.022
  20.861   (  19.239   -3.033   -0.000)   19.477
  22.160   ( -25.991   25.371    0.000)   36.321
  25.158   (  11.227  -21.179    0.000)   23.971
======================= Grid point 515 (148/195) =======================
q-point: (-0.44 -0.50  0.06)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 819
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.311   (   0.000    5.104    0.000)    5.104
   4.353   (   0.000    8.599    0.000)    8.599
   6.317   (  -0.000   -3.620    0.000)    3.620
  21.167   (   0.000    1.216    0.000)    1.216
  21.305   (   0.000   12.919    0.000)   12.919
  25.672   (  -0.000  -11.214    0.000)   11.214
======================= Grid point 754 (149/195) =======================
q-point: ( 0.33  0.22  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.191   (   0.000    7.711   19.798)   21.247
   3.452   (   0.000   16.785   19.928)   26.055
   5.164   (   0.000   14.805   26.197)   30.091
  17.986   (   0.000  -18.049   34.292)   38.751
  19.241   (   0.000   22.722   35.135)   41.842
  26.782   (  -0.000    0.274  -15.186)   15.188
======================= Grid point 755 (150/195) =======================
q-point: ( 0.39  0.22  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.494   (  -0.339    7.169   16.547)   18.036
   3.858   (   2.037   16.181   17.945)   24.249
   5.633   (  -7.827   12.878   18.276)   23.688
  18.210   (  -3.414  -15.947   34.605)   38.255
  19.856   (  13.076   25.769   34.051)   44.660
  26.589   (  -4.492   -4.917  -17.687)   18.900
======================= Grid point 756 (151/195) =======================
q-point: ( 0.44  0.22  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2920
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.746   (   0.342    6.110   13.036)   14.401
   4.190   (   5.217   13.050   16.732)   21.851
   6.064   ( -13.587   10.764    9.913)   19.969
  18.495   (  -4.395  -13.213   34.891)   37.567
  20.323   (  26.005   26.178   32.806)   49.374
  26.365   (  -8.841   -8.559  -19.887)   23.386
======================= Grid point 757 (152/195) =======================
q-point: ( 0.50  0.22  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.931   (   1.461    4.516    9.726)   10.822
   4.429   (   8.133    8.796   15.536)   19.618
   6.419   ( -16.050    7.974    1.850)   18.017
  18.800   (  -1.851  -10.087   33.397)   34.936
  20.625   (  36.275   22.602   32.630)   53.772
  26.146   ( -13.096   -9.085  -21.165)   26.495
======================= Grid point 758 (153/195) =======================
q-point: ( 0.56  0.22  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2920
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.043   (   2.856    2.414    6.888)    7.837
   4.570   (  10.462    4.353   14.079)   18.072
   6.655   ( -15.081    4.287   -5.388)   16.579
  19.051   (   5.351   -5.890   27.542)   28.668
  20.776   (  41.098   13.681   34.956)   55.661
  25.983   ( -16.933   -5.879  -21.168)   27.737
======================= Grid point 759 (154/195) =======================
q-point: (-0.39  0.22  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1467
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.081   (   4.617    0.000    4.617)    6.529
   4.616   (  12.388   -0.000   12.388)   17.520
   6.738   ( -11.265    0.000  -11.265)   15.932
  19.151   (  16.711    0.000   16.711)   23.633
  20.817   (  39.361    0.000   39.361)   55.665
  25.922   ( -19.747    0.000  -19.747)   27.926
======================= Grid point 774 (155/195) =======================
q-point: ( 0.39  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.618   (   0.000    5.599   16.647)   17.563
   3.874   (   0.000   15.527   16.108)   22.373
   5.681   (   0.000    8.227   17.903)   19.703
  18.835   (   0.000  -20.337   37.464)   42.628
  20.042   (   0.000   22.107   32.753)   39.515
  26.405   (  -0.000   -2.756  -16.722)   16.947
======================= Grid point 775 (156/195) =======================
q-point: ( 0.44  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.851   (   1.658    5.308   13.897)   14.968
   4.239   (  -0.378   14.939   14.607)   20.897
   5.965   (  -5.675    6.598   10.137)   13.361
  19.072   (  -2.500  -17.661   38.044)   42.018
  20.621   (  12.125   25.314   30.990)   41.812
  26.154   (  -3.877   -7.836  -19.126)   21.030
======================= Grid point 776 (157/195) =======================
q-point: ( 0.50  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.026   (   3.210    4.718   10.844)   12.254
   4.554   (   2.008   11.135   14.439)   18.344
   6.201   ( -10.413    5.316    1.902)   11.845
  19.347   (  -0.545  -14.263   36.535)   39.224
  21.077   (  22.203   24.661   31.650)   45.857
  25.876   (  -8.365  -10.522  -21.851)   25.655
======================= Grid point 777 (158/195) =======================
q-point: ( 0.56  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.139   (   4.574    3.634    8.006)    9.910
   4.781   (   5.355    5.661   14.747)   16.679
   6.369   ( -12.804    3.943   -5.866)   14.625
  19.564   (   7.589   -9.919   30.083)   32.573
  21.429   (  27.416   17.425   36.483)   48.849
  25.612   ( -13.362   -8.822  -24.612)   29.362
======================= Grid point 778 (159/195) =======================
q-point: ( 0.61  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.191   (   5.962    2.111    5.553)    8.417
   4.904   (   8.380   -0.318   14.729)   16.949
   6.447   ( -12.125    2.090  -11.956)   17.156
  19.609   (  19.915   -4.674   19.178)   28.040
  21.692   (  28.106    2.886   42.884)   51.355
  25.432   ( -18.282   -1.808  -26.232)   32.025
======================= Grid point 779 (160/195) =======================
q-point: (-0.33  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.179   (   7.609    0.290    3.510)    8.384
   4.921   (  10.680   -6.213   13.543)   18.332
   6.421   (  -8.580   -0.187  -15.031)   17.308
  19.444   (  30.455   -0.063    8.638)   31.657
  21.808   (  28.130  -14.133   42.873)   53.189
  25.407   ( -21.953    8.079  -23.843)   33.402
======================= Grid point 780 (161/195) =======================
q-point: (-0.28  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.101   (   9.656   -1.682    1.858)    9.976
   4.823   (  12.295  -11.606   10.327)   19.812
   6.300   (  -2.833   -2.675  -13.830)   14.368
  19.130   (  35.857    3.585    1.828)   36.082
  21.680   (  29.139  -26.177   32.483)   50.887
  25.575   ( -23.153   14.769  -16.241)   31.905
======================= Grid point 781 (162/195) =======================
q-point: (-0.22  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.957   (  12.050   -3.679    0.634)   12.615
   4.602   (  13.570  -15.711    5.106)   21.379
   6.113   (   4.240   -5.480   -8.156)   10.702
  18.767   (  36.930    7.306   -0.831)   37.655
  21.280   (  30.735  -30.462   16.291)   46.238
  25.885   ( -22.084   15.096   -7.143)   27.688
======================= Grid point 782 (163/195) =======================
q-point: (-0.17  0.28  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.751   (  14.501   -5.586   -0.000)   15.540
   4.266   (  15.240  -17.451    0.000)   23.169
   5.888   (  11.568   -8.921   -0.000)   14.608
  18.442   (  35.874   11.595   -0.000)   37.701
  20.650   (  32.571  -29.421   -0.000)   43.891
  26.247   ( -19.952   10.932    0.000)   22.751
======================= Grid point 793 (164/195) =======================
q-point: ( 0.44  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   3.966   (   0.000    5.071   13.185)   14.126
   4.203   (   0.000   14.846   12.094)   19.148
   6.001   (   0.000    2.658    9.671)   10.029
  19.698   (   0.000  -16.748   35.571)   39.317
  20.745   (   0.000   22.277   26.998)   35.002
  26.016   (  -0.000   -7.336  -16.226)   17.808
======================= Grid point 794 (165/195) =======================
q-point: ( 0.50  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.141   (   4.200    5.202   10.815)   12.715
   4.540   (  -4.061   14.236   11.157)   18.538
   6.115   (  -3.293    0.824    2.995)    4.528
  19.937   (   0.869  -13.467   34.902)   37.420
  21.287   (   6.811   25.382   25.702)   36.759
  25.709   (  -2.286  -12.395  -18.710)   22.560
======================= Grid point 795 (166/195) =======================
q-point: ( 0.56  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2920
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.248   (   6.392    5.029    8.021)   11.423
   4.867   (  -3.221   10.197   12.317)   16.311
   6.167   (  -6.349   -0.148   -4.292)    7.665
  20.130   (   9.349   -9.414   28.991)   31.883
  21.803   (   8.557   22.980   30.480)   39.119
  25.351   (  -5.371  -14.372  -22.780)   27.465
======================= Grid point 796 (167/195) =======================
q-point: ( 0.61  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.298   (   7.670    4.234    5.404)   10.293
   5.124   (   0.593    3.491   14.741)   15.160
   6.162   (  -8.729    0.188  -11.429)   14.382
  20.142   (  21.985   -5.314   18.692)   29.343
  22.328   (   9.017   12.215   39.902)   42.693
  24.989   ( -10.217   -9.982  -28.858)   32.199
======================= Grid point 797 (168/195) =======================
q-point: (-0.33  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2920
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.294   (   8.653    2.892    3.129)    9.645
   5.259   (   5.628   -5.300   15.426)   17.255
   6.122   (  -9.665    1.987  -15.156)   18.085
  19.936   (  32.142   -2.515    8.980)   33.468
  22.720   (  15.770   -9.219   43.264)   46.962
  24.753   ( -18.821    5.344  -31.541)   37.116
======================= Grid point 798 (169/195) =======================
q-point: (-0.28  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.236   (   9.664    1.073    1.296)    9.810
   5.211   (   9.247  -12.994    9.168)   18.396
   6.099   (  -7.229    2.739  -10.049)   12.679
  19.572   (  37.354   -0.396    2.598)   37.446
  22.652   (  24.314  -29.463   23.941)   45.082
  24.891   ( -25.330   20.743  -16.248)   36.550
======================= Grid point 799 (170/195) =======================
q-point: (-0.22  0.33  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.124   (  10.833   -1.158   -0.000)   10.895
   4.968   (  11.188  -16.054    0.000)   19.568
   6.111   (  -1.736   -0.406    0.000)    1.783
  19.150   (  38.368    2.621   -0.000)   38.458
  22.100   (  27.955  -33.194    0.000)   43.397
  25.353   ( -24.478   20.740    0.000)   32.083
======================= Grid point 812 (171/195) =======================
q-point: ( 0.50  0.39  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.232   (   0.000    5.720    9.494)   11.084
   4.439   (   0.000   14.479    8.060)   16.571
   6.149   (   0.000   -2.112    3.484)    4.074
  20.460   (   0.000   -9.463   28.813)   30.327
  21.286   (   0.000   22.909   19.026)   29.779
  25.661   (  -0.000  -12.980  -13.798)   18.944
======================= Grid point 813 (172/195) =======================
q-point: ( 0.56  0.39  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.352   (   7.152    6.334    7.085)   11.894
   4.759   (  -8.807   14.194    7.526)   18.321
   6.129   (  -0.685   -4.469   -1.141)    4.663
  20.647   (   8.853   -5.827   24.667)   26.848
  21.811   (  -5.132   25.844   18.853)   32.399
  25.300   (   0.946  -18.498  -15.621)   24.230
======================= Grid point 814 (173/195) =======================
q-point: ( 0.61  0.39  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.394   (   9.173    6.115    4.334)   11.846
   5.123   ( -10.336   11.553    9.112)   17.981
   6.033   (  -0.808   -6.154   -6.236)    8.798
  20.652   (  21.981   -2.821   15.150)   26.844
  22.469   ( -13.035   24.458   24.743)   37.152
  24.840   (   2.590  -21.585  -19.623)   29.286
======================= Grid point 815 (174/195) =======================
q-point: (-0.33  0.39  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.385   (   9.700    5.051    1.917)   11.103
   5.465   (  -7.086    5.378   13.725)   16.356
   5.863   (  -1.876   -5.136  -13.018)   14.120
  20.430   (  31.537   -1.916    6.078)   32.174
  23.229   ( -14.001   16.793   34.059)   40.473
  24.282   (   2.460  -18.337  -29.401)   34.738
======================= Grid point 817 (175/195) =======================
q-point: (-0.28  0.39  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.331   (   9.730    3.346   -0.000)   10.290
   5.532   (   7.373  -13.532    0.000)   15.410
   5.850   ( -11.693    8.344    0.000)   14.365
  20.040   (  36.352   -1.813   -0.000)   36.397
  23.458   (  20.858  -33.086    0.000)   39.112
  24.159   ( -28.476   27.677    0.000)   39.710
======================= Grid point 832 (176/195) =======================
q-point: ( 0.56  0.44  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.405   (   0.000    6.933    5.154)    8.639
   4.580   (  -0.000   14.296    4.029)   14.853
   6.189   (  -0.000   -5.716    0.459)    5.735
  21.002   (   0.000   -1.720   16.719)   16.808
  21.625   (  -0.000   23.548    9.790)   25.502
  25.394   (  -0.000  -18.575   -8.496)   20.426
======================= Grid point 833 (177/195) =======================
q-point: ( 0.61  0.44  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.465   (   9.372    7.509    2.494)   12.265
   4.883   ( -12.783   14.766    2.843)   19.736
   6.094   (   1.476   -8.259   -1.151)    8.468
  21.045   (  17.663    0.581    8.770)   19.729
  22.135   ( -18.457   27.043    7.781)   33.653
  25.024   (   4.926  -24.332   -6.760)   25.730
======================= Grid point 834 (178/195) =======================
q-point: (-0.33 -0.56  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.445   (  10.304    6.673   -0.000)   12.276
   5.244   ( -15.537   14.151    0.000)   21.016
   5.942   (   3.664  -10.477    0.000)   11.100
  20.829   (  27.469   -0.004   -0.000)   27.469
  22.805   ( -27.559   29.402   -0.000)   40.299
  24.568   (  11.193  -28.824    0.000)   30.921
======================= Grid point 853 (179/195) =======================
q-point: (-0.39 -0.50  0.11)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 819
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.466   (   0.000    7.578    0.000)    7.578
   4.627   (   0.000   14.240    0.000)   14.240
   6.191   (  -0.000   -7.142    0.000)    7.142
  21.204   (   0.000    1.820    0.000)    1.820
  21.739   (   0.000   23.811    0.000)   23.811
  25.290   (  -0.000  -21.181    0.000)   21.181
======================= Grid point 1133 (180/195) =======================
q-point: ( 0.50  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.070   (   0.000    3.484   12.804)   13.270
   4.580   (   0.000   16.837   11.605)   20.449
   6.068   (   0.000    2.775    3.980)    4.852
  19.338   (   0.000  -12.753   39.619)   41.621
  21.358   (   0.000   29.620   27.496)   40.415
  25.765   (  -0.000  -14.282  -20.983)   25.383
======================= Grid point 1134 (181/195) =======================
q-point: ( 0.56  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.228   (   2.338    2.021   11.079)   11.501
   4.882   (   2.735   14.721   10.459)   18.264
   6.128   (  -6.310    0.002   -2.657)    6.846
  19.671   (   2.030   -8.451   39.065)   40.020
  21.943   (  13.263   30.093   25.702)   41.738
  25.351   (  -7.064  -17.961  -23.332)   30.280
======================= Grid point 1135 (182/195) =======================
q-point: ( 0.61  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.322   (   4.664    0.932    9.044)   10.218
   5.088   (   7.060    8.716   10.624)   15.449
   6.158   ( -11.394   -0.889   -8.941)   14.510
  19.961   (   9.894   -4.483   33.462)   35.181
  22.338   (  24.681   21.413   27.394)   42.639
  25.004   ( -16.871  -14.331  -25.692)   33.913
======================= Grid point 1136 (183/195) =======================
q-point: (-0.33  0.33  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1467
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.353   (   6.896    0.000    6.896)    9.753
   5.161   (  10.086    0.000   10.086)   14.264
   6.167   ( -12.565    0.000  -12.565)   17.770
  20.078   (  22.215    0.000   22.215)   31.416
  22.483   (  29.600    0.000   29.600)   41.860
  24.860   ( -24.742    0.000  -24.742)   34.990
======================= Grid point 1153 (184/195) =======================
q-point: ( 0.56  0.39  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1548
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.346   (   0.000    3.484   10.627)   11.183
   4.807   (  -0.000   16.482    7.794)   18.232
   6.090   (  -0.000   -2.995   -1.570)    3.382
  20.238   (   0.000   -8.833   36.143)   37.207
  21.915   (  -0.000   30.273   19.766)   36.155
  25.283   (  -0.000  -19.258  -19.791)   27.615
======================= Grid point 1154 (185/195) =======================
q-point: ( 0.61  0.39  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.466   (   5.164    2.890    8.896)   10.684
   5.093   (  -2.146   13.980    7.574)   16.044
   6.012   (  -3.460   -5.470   -6.553)    9.210
  20.518   (   8.829   -4.857   31.669)   33.233
  22.472   (   3.191   29.830   19.418)   35.736
  24.812   (  -4.245  -22.687  -22.153)   31.992
======================= Grid point 1155 (186/195) =======================
q-point: ( 0.67  0.39  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.505   (   8.561    2.747    6.057)   10.841
   5.333   (   0.540    5.461   10.592)   11.929
   5.905   (  -7.573   -3.758  -12.051)   14.721
  20.608   (  23.060   -0.951   20.222)   30.685
  22.977   (   6.179   15.507   27.546)   32.209
  24.377   ( -12.216  -14.539  -27.735)   33.614
======================= Grid point 1156 (187/195) =======================
q-point: (-0.28  0.39  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.475   (  10.384    2.079    2.934)   10.989
   5.387   (   5.938  -10.265    8.153)   14.391
   5.881   ( -10.705    5.117   -9.881)   15.440
  20.434   (  34.512    2.091    7.768)   35.437
  23.133   (  16.675  -23.415   22.170)   36.302
  24.298   ( -24.749   18.300  -19.985)   36.699
======================= Grid point 1157 (188/195) =======================
q-point: (-0.22  0.39  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.381   (  11.094    0.503   -0.000)   11.106
   5.188   (   7.555  -15.770    0.000)   17.486
   6.001   (  -7.127    4.862    0.000)    8.627
  20.075   (  39.318    4.618   -0.000)   39.588
  22.670   (  20.420  -33.705    0.000)   39.408
  24.770   ( -24.945   24.353    0.000)   34.862
======================= Grid point 1172 (189/195) =======================
q-point: ( 0.61  0.44  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.554   (   0.000    5.183    6.737)    8.500
   4.944   (  -0.000   16.297    3.911)   16.759
   6.027   (  -0.000   -8.082   -3.037)    8.634
  20.955   (   0.000   -2.280   23.409)   23.520
  22.269   (  -0.000   30.925   10.308)   32.598
  24.876   (  -0.000  -25.371  -13.883)   28.921
======================= Grid point 1173 (190/195) =======================
q-point: ( 0.67  0.44  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 2916
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.620   (   8.487    5.106    3.673)   10.564
   5.229   (  -9.286   14.499    3.614)   17.593
   5.869   (   0.756  -10.829   -4.263)   11.662
  21.059   (  18.196    0.526   12.726)   22.211
  22.826   ( -14.121   31.316    9.423)   35.621
  24.385   (   3.339  -29.707  -11.769)   32.128
======================= Grid point 1174 (191/195) =======================
q-point: (-0.28  0.44  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1539
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.580   (  10.825    4.263   -0.000)   11.634
   5.553   ( -12.615   12.189    0.000)   17.542
   5.662   (   3.718  -13.424    0.000)   13.930
  20.851   (  30.599    1.685   -0.000)   30.645
  23.515   ( -24.942   30.944   -0.000)   39.745
  23.835   (  11.056  -33.363    0.000)   35.147
======================= Grid point 1191 (192/195) =======================
q-point: (-0.33 -0.50  0.17)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 819
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.637   (   0.000    6.447    0.000)    6.447
   4.989   (   0.000   16.239    0.000)   16.239
   5.985   (  -0.000  -10.458    0.000)   10.458
  21.245   (   0.000    1.557    0.000)    1.557
  22.390   (   0.000   31.192    0.000)   31.192
  24.700   (  -0.000  -28.862    0.000)   28.862
======================= Grid point 1508 (193/195) =======================
q-point: ( 0.67  0.44  0.22)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1550
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.625   (   0.000    0.616    9.023)    9.044
   5.321   (  -0.000   15.616    3.802)   16.073
   5.814   (  -0.000  -10.137   -7.828)   12.808
  20.900   (   0.000   -2.350   28.785)   28.881
  23.038   (  -0.000   34.197   10.760)   35.850
  24.228   (  -0.000  -29.602  -19.824)   35.627
======================= Grid point 1509 (194/195) =======================
q-point: (-0.28  0.44  0.22)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 1467
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.682   (   5.891    0.000    5.891)    8.331
   5.497   (   3.823    0.000    3.823)    5.406
   5.646   (  -8.108   -0.000   -8.108)   11.466
  21.066   (  16.222    0.000   16.222)   22.941
  23.443   (  10.631    0.000   10.631)   15.035
  23.788   ( -17.195   -0.000  -17.195)   24.317
======================= Grid point 1528 (195/195) =======================
q-point: (-0.28 -0.50  0.22)
Boundary mean free path (millimeter): 1000.000
Mass variance parameters: 1.15e-04 0.00e+00 
Number of triplets: 819
Calculating Gamma of ph-isotope with tetrahedron method
Frequency     group velocity (x, y, z)     |gv|
   4.745   (   0.000    2.495    0.000)    2.495
   5.365   (   0.000   15.568    0.000)   15.568
   5.705   (  -0.000  -13.391    0.000)   13.391
  21.271   (   0.000    0.604    0.000)    0.604
  23.164   (   0.000   34.397    0.000)   34.397
  23.965   (  -0.000  -33.473    0.000)   33.473
=================== 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/34992
   10.0   8258.695   8258.695   8258.695     -0.000     -0.000      0.000 3/34992
   20.0   1595.852   1595.852   1595.852     -0.000     -0.000      0.000 3/34992
   30.0    472.330    472.330    472.330     -0.000      0.000      0.000 3/34992
   40.0    238.966    238.966    238.966     -0.000      0.000      0.000 3/34992
   50.0    167.995    167.995    167.995     -0.000      0.000      0.000 3/34992
   60.0    137.393    137.393    137.393     -0.000      0.000      0.000 3/34992
   70.0    119.904    119.904    119.904     -0.000      0.000      0.000 3/34992
   80.0    107.692    107.692    107.692     -0.000     -0.000      0.000 3/34992
   90.0     98.000     98.000     98.000     -0.000     -0.000      0.000 3/34992
  100.0     89.676     89.676     89.676     -0.000     -0.000      0.000 3/34992
  110.0     82.187     82.187     82.187     -0.000     -0.000      0.000 3/34992
  120.0     75.293     75.293     75.293     -0.000     -0.000      0.000 3/34992
  130.0     68.905     68.905     68.905     -0.000     -0.000      0.000 3/34992
  140.0     63.007     63.007     63.007     -0.000     -0.000      0.000 3/34992
  150.0     57.607     57.607     57.607     -0.000     -0.000      0.000 3/34992
  160.0     52.711     52.711     52.711     -0.000     -0.000      0.000 3/34992
  170.0     48.309     48.309     48.309     -0.000     -0.000      0.000 3/34992
  180.0     44.376     44.376     44.376     -0.000     -0.000      0.000 3/34992
  190.0     40.877     40.877     40.877     -0.000     -0.000      0.000 3/34992
  200.0     37.771     37.771     37.771     -0.000     -0.000      0.000 3/34992
  210.0     35.015     35.015     35.015     -0.000     -0.000      0.000 3/34992
  220.0     32.569     32.569     32.569     -0.000     -0.000      0.000 3/34992
  230.0     30.394     30.394     30.394     -0.000     -0.000      0.000 3/34992
  240.0     28.455     28.455     28.455     -0.000     -0.000      0.000 3/34992
  250.0     26.723     26.723     26.723     -0.000     -0.000      0.000 3/34992
  260.0     25.170     25.170     25.170     -0.000     -0.000      0.000 3/34992
  270.0     23.774     23.774     23.774     -0.000     -0.000      0.000 3/34992
  280.0     22.515     22.515     22.515     -0.000     -0.000      0.000 3/34992
  290.0     21.375     21.375     21.375     -0.000     -0.000      0.000 3/34992
  300.0     20.339     20.339     20.339     -0.000     -0.000      0.000 3/34992
  310.0     19.396     19.396     19.396     -0.000     -0.000      0.000 3/34992
  320.0     18.534     18.534     18.534     -0.000     -0.000      0.000 3/34992
  330.0     17.743     17.743     17.743     -0.000     -0.000      0.000 3/34992
  340.0     17.017     17.017     17.017     -0.000     -0.000      0.000 3/34992
  350.0     16.347     16.347     16.347     -0.000     -0.000      0.000 3/34992
  360.0     15.727     15.727     15.727     -0.000     -0.000      0.000 3/34992
  370.0     15.153     15.153     15.153     -0.000     -0.000      0.000 3/34992
  380.0     14.620     14.620     14.620     -0.000     -0.000      0.000 3/34992
  390.0     14.123     14.123     14.123     -0.000     -0.000      0.000 3/34992
  400.0     13.660     13.660     13.660     -0.000     -0.000      0.000 3/34992
  410.0     13.226     13.226     13.226     -0.000     -0.000      0.000 3/34992
  420.0     12.820     12.820     12.820     -0.000     -0.000      0.000 3/34992
  430.0     12.438     12.438     12.438     -0.000     -0.000      0.000 3/34992
  440.0     12.079     12.079     12.079     -0.000     -0.000      0.000 3/34992
  450.0     11.741     11.741     11.741     -0.000     -0.000      0.000 3/34992
  460.0     11.421     11.421     11.421     -0.000     -0.000      0.000 3/34992
  470.0     11.119     11.119     11.119     -0.000     -0.000      0.000 3/34992
  480.0     10.833     10.833     10.833     -0.000     -0.000      0.000 3/34992
  490.0     10.562     10.562     10.562     -0.000     -0.000      0.000 3/34992
  500.0     10.305     10.305     10.305     -0.000     -0.000      0.000 3/34992
  510.0     10.060     10.060     10.060     -0.000     -0.000      0.000 3/34992
  520.0      9.827      9.827      9.827     -0.000     -0.000      0.000 3/34992
  530.0      9.605      9.605      9.605     -0.000     -0.000      0.000 3/34992
  540.0      9.393      9.393      9.393     -0.000     -0.000      0.000 3/34992
  550.0      9.191      9.191      9.191     -0.000     -0.000      0.000 3/34992
  560.0      8.997      8.997      8.997     -0.000     -0.000      0.000 3/34992
  570.0      8.812      8.812      8.812     -0.000     -0.000      0.000 3/34992
  580.0      8.634      8.634      8.634     -0.000     -0.000      0.000 3/34992
  590.0      8.464      8.464      8.464     -0.000     -0.000      0.000 3/34992
  600.0      8.300      8.300      8.300     -0.000     -0.000      0.000 3/34992
  610.0      8.143      8.143      8.143     -0.000     -0.000      0.000 3/34992
  620.0      7.992      7.992      7.992     -0.000     -0.000      0.000 3/34992
  630.0      7.847      7.847      7.847     -0.000     -0.000      0.000 3/34992
  640.0      7.707      7.707      7.707     -0.000     -0.000      0.000 3/34992
  650.0      7.572      7.572      7.572     -0.000     -0.000      0.000 3/34992
  660.0      7.442      7.442      7.442     -0.000     -0.000      0.000 3/34992
  670.0      7.317      7.317      7.317     -0.000     -0.000      0.000 3/34992
  680.0      7.195      7.195      7.195     -0.000     -0.000      0.000 3/34992
  690.0      7.078      7.078      7.078     -0.000     -0.000      0.000 3/34992
  700.0      6.965      6.965      6.965     -0.000     -0.000      0.000 3/34992
  710.0      6.856      6.856      6.856     -0.000     -0.000      0.000 3/34992
  720.0      6.749      6.749      6.749     -0.000     -0.000      0.000 3/34992
  730.0      6.647      6.647      6.647     -0.000     -0.000      0.000 3/34992
  740.0      6.547      6.547      6.547     -0.000     -0.000      0.000 3/34992
  750.0      6.451      6.451      6.451     -0.000     -0.000      0.000 3/34992
  760.0      6.357      6.357      6.357     -0.000     -0.000      0.000 3/34992
  770.0      6.266      6.266      6.266     -0.000     -0.000      0.000 3/34992
  780.0      6.178      6.178      6.178     -0.000     -0.000      0.000 3/34992
  790.0      6.092      6.092      6.092     -0.000     -0.000      0.000 3/34992
  800.0      6.009      6.009      6.009     -0.000     -0.000      0.000 3/34992
  810.0      5.928      5.928      5.928     -0.000     -0.000      0.000 3/34992
  820.0      5.850      5.850      5.850     -0.000     -0.000      0.000 3/34992
  830.0      5.773      5.773      5.773     -0.000     -0.000      0.000 3/34992
  840.0      5.698      5.698      5.698     -0.000     -0.000      0.000 3/34992
  850.0      5.626      5.626      5.626     -0.000     -0.000      0.000 3/34992
  860.0      5.555      5.555      5.555     -0.000     -0.000      0.000 3/34992
  870.0      5.486      5.486      5.486     -0.000     -0.000      0.000 3/34992
  880.0      5.419      5.419      5.419     -0.000     -0.000      0.000 3/34992
  890.0      5.353      5.353      5.353     -0.000     -0.000      0.000 3/34992
  900.0      5.290      5.290      5.290     -0.000     -0.000      0.000 3/34992
  910.0      5.227      5.227      5.227     -0.000     -0.000      0.000 3/34992
  920.0      5.166      5.166      5.166     -0.000     -0.000      0.000 3/34992
  930.0      5.107      5.107      5.107     -0.000     -0.000      0.000 3/34992
  940.0      5.049      5.049      5.049     -0.000     -0.000      0.000 3/34992
  950.0      4.992      4.992      4.992     -0.000     -0.000      0.000 3/34992
  960.0      4.937      4.937      4.937     -0.000     -0.000      0.000 3/34992
  970.0      4.882      4.882      4.882     -0.000     -0.000      0.000 3/34992
  980.0      4.830      4.830      4.830     -0.000     -0.000      0.000 3/34992
  990.0      4.778      4.778      4.778     -0.000     -0.000      0.000 3/34992
 1000.0      4.727      4.727      4.727     -0.000     -0.000      0.000 3/34992

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

