# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/16508778-7389-4d21-b0b1-fe72cabf2f28/download)

## Creator

[Atsushi Togo](https://orcid.org/0000-0001-8393-9766)

## Rights

Creative Commons Attribution 4.0 International[Creative Commons BY Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/)

## Other metadata

[First-principles lattice thermal conductivity calculation for K4I2O / I4/mmm (139) / materials id 768368](https://mdr.nims.go.jp/datasets/67a695b8-2762-4c56-8569-3a433abf09fa)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-09 00:53:34]-------------------------Compiled with OpenMP support (max 128 threads).Running in phonopy.load mode.Python version 3.14.2Spglib version 2.6.1Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".Unit of length: angstromSettings:  Supercell: [3 3 1]  Primitive matrix:    [-0.5  0.5  0.5]    [ 0.5 -0.5  0.5]    [ 0.5  0.5 -0.5]Spacegroup: I4/mmm (139)Number of symmetry operations in supercell: 288------------------------------ primitive cell ------------------------------Lattice vectors:  a   -2.642727035000000    2.642727035000000    8.727389710000001  b    2.642727035000000   -2.642727035000000    8.727389710000001  c    2.642727035000000    2.642727035000000   -8.727389710000001Atomic positions (fractional):   *1 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098   *2 K   0.14793473746599  0.14793473746599  0.00000000000000  39.098    3 K   0.85206526253401  0.85206526253401  0.00000000000000  39.098    4 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098   *5 I   0.35340590890074  0.35340590890074  0.00000000000000 126.904    6 I   0.64659409109926  0.64659409109926  0.00000000000000 126.904   *7 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999-------------------------------- unit cell ---------------------------------Lattice vectors:  a    5.285454070000000    0.000000000000000    0.000000000000000  b    0.000000000000000    5.285454070000000    0.000000000000000  c    0.000000000000000    0.000000000000000   17.454779420000001Atomic positions (fractional):   *1 K   0.50000000000000  0.00000000000000  0.00000000000000  39.098 > 1   *2 K   0.00000000000000  0.00000000000000  0.14793473746599  39.098 > 2    3 K   0.50000000000000  0.50000000000000  0.35206526253401  39.098 > 3    4 K   0.00000000000000  0.50000000000000  0.00000000000000  39.098 > 4    5 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098 > 1    6 K   0.50000000000000  0.50000000000000  0.64793473746599  39.098 > 2    7 K   0.00000000000000  0.00000000000000  0.85206526253401  39.098 > 3    8 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098 > 4   *9 I   0.00000000000000  0.00000000000000  0.35340590890074 126.904 > 5   10 I   0.50000000000000  0.50000000000000  0.14659409109926 126.904 > 6   11 I   0.50000000000000  0.50000000000000  0.85340590890074 126.904 > 5   12 I   0.00000000000000  0.00000000000000  0.64659409109926 126.904 > 6  *13 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 7   14 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 7-------------------------------- super cell --------------------------------Lattice vectors:  a   15.856362210000000    0.000000000000000    0.000000000000000  b    0.000000000000000   15.856362210000000    0.000000000000000  c    0.000000000000000    0.000000000000000   17.454779420000001Atomic positions (fractional):   *1 K   0.16666666666667  0.00000000000000  0.00000000000000  39.098 > 1    2 K   0.50000000000000  0.00000000000000  0.00000000000000  39.098 > 1    3 K   0.83333333333333  0.00000000000000  0.00000000000000  39.098 > 1    4 K   0.16666666666667  0.33333333333333  0.00000000000000  39.098 > 1    5 K   0.50000000000000  0.33333333333333  0.00000000000000  39.098 > 1    6 K   0.83333333333333  0.33333333333333  0.00000000000000  39.098 > 1    7 K   0.16666666666667  0.66666666666667  0.00000000000000  39.098 > 1    8 K   0.50000000000000  0.66666666666667  0.00000000000000  39.098 > 1    9 K   0.83333333333333  0.66666666666667  0.00000000000000  39.098 > 1  *10 K   0.00000000000000  0.00000000000000  0.14793473746599  39.098 > 2   11 K   0.33333333333333  0.00000000000000  0.14793473746599  39.098 > 2   12 K   0.66666666666667  0.00000000000000  0.14793473746599  39.098 > 2   13 K   0.00000000000000  0.33333333333333  0.14793473746599  39.098 > 2   14 K   0.33333333333333  0.33333333333333  0.14793473746599  39.098 > 2   15 K   0.66666666666667  0.33333333333333  0.14793473746599  39.098 > 2   16 K   0.00000000000000  0.66666666666667  0.14793473746599  39.098 > 2   17 K   0.33333333333333  0.66666666666667  0.14793473746599  39.098 > 2   18 K   0.66666666666667  0.66666666666667  0.14793473746599  39.098 > 2   19 K   0.16666666666667  0.16666666666667  0.35206526253401  39.098 > 3   20 K   0.50000000000000  0.16666666666667  0.35206526253401  39.098 > 3   21 K   0.83333333333333  0.16666666666667  0.35206526253401  39.098 > 3   22 K   0.16666666666667  0.50000000000000  0.35206526253401  39.098 > 3   23 K   0.50000000000000  0.50000000000000  0.35206526253401  39.098 > 3   24 K   0.83333333333333  0.50000000000000  0.35206526253401  39.098 > 3   25 K   0.16666666666667  0.83333333333333  0.35206526253401  39.098 > 3   26 K   0.50000000000000  0.83333333333333  0.35206526253401  39.098 > 3   27 K   0.83333333333333  0.83333333333333  0.35206526253401  39.098 > 3   28 K   0.00000000000000  0.16666666666667  0.00000000000000  39.098 > 4   29 K   0.33333333333333  0.16666666666667  0.00000000000000  39.098 > 4   30 K   0.66666666666667  0.16666666666667  0.00000000000000  39.098 > 4   31 K   0.00000000000000  0.50000000000000  0.00000000000000  39.098 > 4   32 K   0.33333333333333  0.50000000000000  0.00000000000000  39.098 > 4   33 K   0.66666666666667  0.50000000000000  0.00000000000000  39.098 > 4   34 K   0.00000000000000  0.83333333333333  0.00000000000000  39.098 > 4   35 K   0.33333333333333  0.83333333333333  0.00000000000000  39.098 > 4   36 K   0.66666666666667  0.83333333333333  0.00000000000000  39.098 > 4   37 K   0.00000000000000  0.16666666666667  0.50000000000000  39.098 > 1   38 K   0.33333333333333  0.16666666666667  0.50000000000000  39.098 > 1   39 K   0.66666666666667  0.16666666666667  0.50000000000000  39.098 > 1   40 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098 > 1   41 K   0.33333333333333  0.50000000000000  0.50000000000000  39.098 > 1   42 K   0.66666666666667  0.50000000000000  0.50000000000000  39.098 > 1   43 K   0.00000000000000  0.83333333333333  0.50000000000000  39.098 > 1   44 K   0.33333333333333  0.83333333333333  0.50000000000000  39.098 > 1   45 K   0.66666666666667  0.83333333333333  0.50000000000000  39.098 > 1   46 K   0.16666666666667  0.16666666666667  0.64793473746599  39.098 > 2   47 K   0.50000000000000  0.16666666666667  0.64793473746599  39.098 > 2   48 K   0.83333333333333  0.16666666666667  0.64793473746599  39.098 > 2   49 K   0.16666666666667  0.50000000000000  0.64793473746599  39.098 > 2   50 K   0.50000000000000  0.50000000000000  0.64793473746599  39.098 > 2   51 K   0.83333333333333  0.50000000000000  0.64793473746599  39.098 > 2   52 K   0.16666666666667  0.83333333333333  0.64793473746599  39.098 > 2   53 K   0.50000000000000  0.83333333333333  0.64793473746599  39.098 > 2   54 K   0.83333333333333  0.83333333333333  0.64793473746599  39.098 > 2   55 K   0.00000000000000  0.00000000000000  0.85206526253401  39.098 > 3   56 K   0.33333333333333  0.00000000000000  0.85206526253401  39.098 > 3   57 K   0.66666666666667  0.00000000000000  0.85206526253401  39.098 > 3   58 K   0.00000000000000  0.33333333333333  0.85206526253401  39.098 > 3   59 K   0.33333333333333  0.33333333333333  0.85206526253401  39.098 > 3   60 K   0.66666666666667  0.33333333333333  0.85206526253401  39.098 > 3   61 K   0.00000000000000  0.66666666666667  0.85206526253401  39.098 > 3   62 K   0.33333333333333  0.66666666666667  0.85206526253401  39.098 > 3   63 K   0.66666666666667  0.66666666666667  0.85206526253401  39.098 > 3   64 K   0.16666666666667  0.00000000000000  0.50000000000000  39.098 > 4   65 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098 > 4   66 K   0.83333333333333  0.00000000000000  0.50000000000000  39.098 > 4   67 K   0.16666666666667  0.33333333333333  0.50000000000000  39.098 > 4   68 K   0.50000000000000  0.33333333333333  0.50000000000000  39.098 > 4   69 K   0.83333333333333  0.33333333333333  0.50000000000000  39.098 > 4   70 K   0.16666666666667  0.66666666666667  0.50000000000000  39.098 > 4   71 K   0.50000000000000  0.66666666666667  0.50000000000000  39.098 > 4   72 K   0.83333333333333  0.66666666666667  0.50000000000000  39.098 > 4  *73 I   0.00000000000000  0.00000000000000  0.35340590890074 126.904 > 5   74 I   0.33333333333333  0.00000000000000  0.35340590890074 126.904 > 5   75 I   0.66666666666667  0.00000000000000  0.35340590890074 126.904 > 5   76 I   0.00000000000000  0.33333333333333  0.35340590890074 126.904 > 5   77 I   0.33333333333333  0.33333333333333  0.35340590890074 126.904 > 5   78 I   0.66666666666667  0.33333333333333  0.35340590890074 126.904 > 5   79 I   0.00000000000000  0.66666666666667  0.35340590890074 126.904 > 5   80 I   0.33333333333333  0.66666666666667  0.35340590890074 126.904 > 5   81 I   0.66666666666667  0.66666666666667  0.35340590890074 126.904 > 5   82 I   0.16666666666667  0.16666666666667  0.14659409109926 126.904 > 6   83 I   0.50000000000000  0.16666666666667  0.14659409109926 126.904 > 6   84 I   0.83333333333333  0.16666666666667  0.14659409109926 126.904 > 6   85 I   0.16666666666667  0.50000000000000  0.14659409109926 126.904 > 6   86 I   0.50000000000000  0.50000000000000  0.14659409109926 126.904 > 6   87 I   0.83333333333333  0.50000000000000  0.14659409109926 126.904 > 6   88 I   0.16666666666667  0.83333333333333  0.14659409109926 126.904 > 6   89 I   0.50000000000000  0.83333333333333  0.14659409109926 126.904 > 6   90 I   0.83333333333333  0.83333333333333  0.14659409109926 126.904 > 6   91 I   0.16666666666667  0.16666666666667  0.85340590890074 126.904 > 5   92 I   0.50000000000000  0.16666666666667  0.85340590890074 126.904 > 5   93 I   0.83333333333333  0.16666666666667  0.85340590890074 126.904 > 5   94 I   0.16666666666667  0.50000000000000  0.85340590890074 126.904 > 5   95 I   0.50000000000000  0.50000000000000  0.85340590890074 126.904 > 5   96 I   0.83333333333333  0.50000000000000  0.85340590890074 126.904 > 5   97 I   0.16666666666667  0.83333333333333  0.85340590890074 126.904 > 5   98 I   0.50000000000000  0.83333333333333  0.85340590890074 126.904 > 5   99 I   0.83333333333333  0.83333333333333  0.85340590890074 126.904 > 5  100 I   0.00000000000000  0.00000000000000  0.64659409109926 126.904 > 6  101 I   0.33333333333333  0.00000000000000  0.64659409109926 126.904 > 6  102 I   0.66666666666667  0.00000000000000  0.64659409109926 126.904 > 6  103 I   0.00000000000000  0.33333333333333  0.64659409109926 126.904 > 6  104 I   0.33333333333333  0.33333333333333  0.64659409109926 126.904 > 6  105 I   0.66666666666667  0.33333333333333  0.64659409109926 126.904 > 6  106 I   0.00000000000000  0.66666666666667  0.64659409109926 126.904 > 6  107 I   0.33333333333333  0.66666666666667  0.64659409109926 126.904 > 6  108 I   0.66666666666667  0.66666666666667  0.64659409109926 126.904 > 6 *109 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 7  110 O   0.33333333333333  0.00000000000000  0.00000000000000  15.999 > 7  111 O   0.66666666666667  0.00000000000000  0.00000000000000  15.999 > 7  112 O   0.00000000000000  0.33333333333333  0.00000000000000  15.999 > 7  113 O   0.33333333333333  0.33333333333333  0.00000000000000  15.999 > 7  114 O   0.66666666666667  0.33333333333333  0.00000000000000  15.999 > 7  115 O   0.00000000000000  0.66666666666667  0.00000000000000  15.999 > 7  116 O   0.33333333333333  0.66666666666667  0.00000000000000  15.999 > 7  117 O   0.66666666666667  0.66666666666667  0.00000000000000  15.999 > 7  118 O   0.16666666666667  0.16666666666667  0.50000000000000  15.999 > 7  119 O   0.50000000000000  0.16666666666667  0.50000000000000  15.999 > 7  120 O   0.83333333333333  0.16666666666667  0.50000000000000  15.999 > 7  121 O   0.16666666666667  0.50000000000000  0.50000000000000  15.999 > 7  122 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 7  123 O   0.83333333333333  0.50000000000000  0.50000000000000  15.999 > 7  124 O   0.16666666666667  0.83333333333333  0.50000000000000  15.999 > 7  125 O   0.50000000000000  0.83333333333333  0.50000000000000  15.999 > 7  126 O   0.83333333333333  0.83333333333333  0.50000000000000  15.999 > 7----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            3.0759119    0.0000000    0.0000000            0.0000000    3.0759119    0.0000000            0.0000000    0.0000000    3.0891469-------------------------- Born effective charges --------------------------    1 K     1.1710065    0.0000000    0.0000000            0.0000000    1.1424984    0.0000000            0.0000000    0.0000000    0.9864279    2 K     1.1290615    0.0000000    0.0000000            0.0000000    1.1290615    0.0000000            0.0000000    0.0000000    1.2349842    3 K     1.1290615    0.0000000    0.0000000            0.0000000    1.1290615    0.0000000            0.0000000    0.0000000    1.2349842    4 K     1.1424984    0.0000000    0.0000000            0.0000000    1.1710065    0.0000000            0.0000000    0.0000000    0.9864279    5 I    -1.1615912    0.0000000    0.0000000            0.0000000   -1.1615912    0.0000000            0.0000000    0.0000000   -1.1299427    6 I    -1.1615912    0.0000000    0.0000000            0.0000000   -1.1615912    0.0000000            0.0000000    0.0000000   -1.1299427    7 O    -2.2484455    0.0000000    0.0000000            0.0000000   -2.2484455    0.0000000            0.0000000    0.0000000   -2.1829389----------------------------------------------------------------------------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: 378/378Permutation basis: 5061/5061Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 352Number of blocks in projector: 358Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 6--- Eigsh_solver_block: 1 / 6 ---Block_size: 144Use standard eigh solver.--- Eigsh_solver_block: 2 / 6 ---Block_size: 92Use standard eigh solver.--- Eigsh_solver_block: 3 / 6 ---Block_size: 80Use standard eigh solver.--- Eigsh_solver_block: 4 / 6 ---Block_size: 34Use standard eigh solver.--- Eigsh_solver_block: 5 / 6 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 6 / 6 ---Block_size: 1Use standard eigh solver.Tree of FC basis block matrices:- (352, 343), data: False|-- (1, 1), data: True|-- (1, 1), data: True|-- (34, 34), data: True|-- (80, 80), data: True|-- (92, 88), data: True|-- (144, 139), data: True-----Solver_atoms: 1 -- 126 / 126Time (Solver_compr_matrix_reshape): 0.005Solver_block: 80 / 80 - Time: 0.264Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.271--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 378/378Permutation basis: 5061/5061Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 352Number of blocks in projector: 358Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 6--- Eigsh_solver_block: 1 / 6 ---Block_size: 144Use standard eigh solver.--- Eigsh_solver_block: 2 / 6 ---Block_size: 92Use standard eigh solver.--- Eigsh_solver_block: 3 / 6 ---Block_size: 80Use standard eigh solver.--- Eigsh_solver_block: 4 / 6 ---Block_size: 34Use standard eigh solver.--- Eigsh_solver_block: 5 / 6 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 6 / 6 ---Block_size: 1Use standard eigh solver.Tree of FC basis block matrices:- (352, 343), data: False|-- (1, 1), data: True|-- (1, 1), data: True|-- (34, 34), data: True|-- (80, 80), data: True|-- (92, 88), data: True|-- (144, 139), data: TrueMax drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-09 00:53:38]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-09 00:53:39]-------------------------Compiled with OpenMP support (max 128 threads).Running in phono3py.load mode.Python version 3.14.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: force constantsHDF5 data compression filter: gzipCrystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".Supercell (dim): [3 3 1]Primitive matrix:  [-0.5  0.5  0.5]  [ 0.5 -0.5  0.5]  [ 0.5  0.5 -0.5]Spacegroup: I4/mmm (139)------------------------------ primitive cell ------------------------------Lattice vectors:  a   -2.642727035000000    2.642727035000000    8.727389710000001  b    2.642727035000000   -2.642727035000000    8.727389710000001  c    2.642727035000000    2.642727035000000   -8.727389710000001Atomic positions (fractional):    1 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098    2 K   0.14793473746599  0.14793473746599  0.00000000000000  39.098    3 K   0.85206526253401  0.85206526253401  0.00000000000000  39.098    4 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098    5 I   0.35340590890074  0.35340590890074  0.00000000000000 126.904    6 I   0.64659409109926  0.64659409109926  0.00000000000000 126.904    7 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999-------------------------------- supercell ---------------------------------Lattice vectors:  a   15.856362210000000    0.000000000000000    0.000000000000000  b    0.000000000000000   15.856362210000000    0.000000000000000  c    0.000000000000000    0.000000000000000   17.454779420000001Atomic positions (fractional):    1 K   0.16666666666667  0.00000000000000  0.00000000000000  39.098 > 1    2 K   0.50000000000000  0.00000000000000  0.00000000000000  39.098 > 1    3 K   0.83333333333333  0.00000000000000  0.00000000000000  39.098 > 1    4 K   0.16666666666667  0.33333333333333  0.00000000000000  39.098 > 1    5 K   0.50000000000000  0.33333333333333  0.00000000000000  39.098 > 1    6 K   0.83333333333333  0.33333333333333  0.00000000000000  39.098 > 1    7 K   0.16666666666667  0.66666666666667  0.00000000000000  39.098 > 1    8 K   0.50000000000000  0.66666666666667  0.00000000000000  39.098 > 1    9 K   0.83333333333333  0.66666666666667  0.00000000000000  39.098 > 1   10 K   0.00000000000000  0.00000000000000  0.14793473746599  39.098 > 10   11 K   0.33333333333333  0.00000000000000  0.14793473746599  39.098 > 10   12 K   0.66666666666667  0.00000000000000  0.14793473746599  39.098 > 10   13 K   0.00000000000000  0.33333333333333  0.14793473746599  39.098 > 10   14 K   0.33333333333333  0.33333333333333  0.14793473746599  39.098 > 10   15 K   0.66666666666667  0.33333333333333  0.14793473746599  39.098 > 10   16 K   0.00000000000000  0.66666666666667  0.14793473746599  39.098 > 10   17 K   0.33333333333333  0.66666666666667  0.14793473746599  39.098 > 10   18 K   0.66666666666667  0.66666666666667  0.14793473746599  39.098 > 10   19 K   0.16666666666667  0.16666666666667  0.35206526253401  39.098 > 19   20 K   0.50000000000000  0.16666666666667  0.35206526253401  39.098 > 19   21 K   0.83333333333333  0.16666666666667  0.35206526253401  39.098 > 19   22 K   0.16666666666667  0.50000000000000  0.35206526253401  39.098 > 19   23 K   0.50000000000000  0.50000000000000  0.35206526253401  39.098 > 19   24 K   0.83333333333333  0.50000000000000  0.35206526253401  39.098 > 19   25 K   0.16666666666667  0.83333333333333  0.35206526253401  39.098 > 19   26 K   0.50000000000000  0.83333333333333  0.35206526253401  39.098 > 19   27 K   0.83333333333333  0.83333333333333  0.35206526253401  39.098 > 19   28 K   0.00000000000000  0.16666666666667  0.00000000000000  39.098 > 28   29 K   0.33333333333333  0.16666666666667  0.00000000000000  39.098 > 28   30 K   0.66666666666667  0.16666666666667  0.00000000000000  39.098 > 28   31 K   0.00000000000000  0.50000000000000  0.00000000000000  39.098 > 28   32 K   0.33333333333333  0.50000000000000  0.00000000000000  39.098 > 28   33 K   0.66666666666667  0.50000000000000  0.00000000000000  39.098 > 28   34 K   0.00000000000000  0.83333333333333  0.00000000000000  39.098 > 28   35 K   0.33333333333333  0.83333333333333  0.00000000000000  39.098 > 28   36 K   0.66666666666667  0.83333333333333  0.00000000000000  39.098 > 28   37 K   0.00000000000000  0.16666666666667  0.50000000000000  39.098 > 1   38 K   0.33333333333333  0.16666666666667  0.50000000000000  39.098 > 1   39 K   0.66666666666667  0.16666666666667  0.50000000000000  39.098 > 1   40 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098 > 1   41 K   0.33333333333333  0.50000000000000  0.50000000000000  39.098 > 1   42 K   0.66666666666667  0.50000000000000  0.50000000000000  39.098 > 1   43 K   0.00000000000000  0.83333333333333  0.50000000000000  39.098 > 1   44 K   0.33333333333333  0.83333333333333  0.50000000000000  39.098 > 1   45 K   0.66666666666667  0.83333333333333  0.50000000000000  39.098 > 1   46 K   0.16666666666667  0.16666666666667  0.64793473746599  39.098 > 10   47 K   0.50000000000000  0.16666666666667  0.64793473746599  39.098 > 10   48 K   0.83333333333333  0.16666666666667  0.64793473746599  39.098 > 10   49 K   0.16666666666667  0.50000000000000  0.64793473746599  39.098 > 10   50 K   0.50000000000000  0.50000000000000  0.64793473746599  39.098 > 10   51 K   0.83333333333333  0.50000000000000  0.64793473746599  39.098 > 10   52 K   0.16666666666667  0.83333333333333  0.64793473746599  39.098 > 10   53 K   0.50000000000000  0.83333333333333  0.64793473746599  39.098 > 10   54 K   0.83333333333333  0.83333333333333  0.64793473746599  39.098 > 10   55 K   0.00000000000000  0.00000000000000  0.85206526253401  39.098 > 19   56 K   0.33333333333333  0.00000000000000  0.85206526253401  39.098 > 19   57 K   0.66666666666667  0.00000000000000  0.85206526253401  39.098 > 19   58 K   0.00000000000000  0.33333333333333  0.85206526253401  39.098 > 19   59 K   0.33333333333333  0.33333333333333  0.85206526253401  39.098 > 19   60 K   0.66666666666667  0.33333333333333  0.85206526253401  39.098 > 19   61 K   0.00000000000000  0.66666666666667  0.85206526253401  39.098 > 19   62 K   0.33333333333333  0.66666666666667  0.85206526253401  39.098 > 19   63 K   0.66666666666667  0.66666666666667  0.85206526253401  39.098 > 19   64 K   0.16666666666667  0.00000000000000  0.50000000000000  39.098 > 28   65 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098 > 28   66 K   0.83333333333333  0.00000000000000  0.50000000000000  39.098 > 28   67 K   0.16666666666667  0.33333333333333  0.50000000000000  39.098 > 28   68 K   0.50000000000000  0.33333333333333  0.50000000000000  39.098 > 28   69 K   0.83333333333333  0.33333333333333  0.50000000000000  39.098 > 28   70 K   0.16666666666667  0.66666666666667  0.50000000000000  39.098 > 28   71 K   0.50000000000000  0.66666666666667  0.50000000000000  39.098 > 28   72 K   0.83333333333333  0.66666666666667  0.50000000000000  39.098 > 28   73 I   0.00000000000000  0.00000000000000  0.35340590890074 126.904 > 73   74 I   0.33333333333333  0.00000000000000  0.35340590890074 126.904 > 73   75 I   0.66666666666667  0.00000000000000  0.35340590890074 126.904 > 73   76 I   0.00000000000000  0.33333333333333  0.35340590890074 126.904 > 73   77 I   0.33333333333333  0.33333333333333  0.35340590890074 126.904 > 73   78 I   0.66666666666667  0.33333333333333  0.35340590890074 126.904 > 73   79 I   0.00000000000000  0.66666666666667  0.35340590890074 126.904 > 73   80 I   0.33333333333333  0.66666666666667  0.35340590890074 126.904 > 73   81 I   0.66666666666667  0.66666666666667  0.35340590890074 126.904 > 73   82 I   0.16666666666667  0.16666666666667  0.14659409109926 126.904 > 82   83 I   0.50000000000000  0.16666666666667  0.14659409109926 126.904 > 82   84 I   0.83333333333333  0.16666666666667  0.14659409109926 126.904 > 82   85 I   0.16666666666667  0.50000000000000  0.14659409109926 126.904 > 82   86 I   0.50000000000000  0.50000000000000  0.14659409109926 126.904 > 82   87 I   0.83333333333333  0.50000000000000  0.14659409109926 126.904 > 82   88 I   0.16666666666667  0.83333333333333  0.14659409109926 126.904 > 82   89 I   0.50000000000000  0.83333333333333  0.14659409109926 126.904 > 82   90 I   0.83333333333333  0.83333333333333  0.14659409109926 126.904 > 82   91 I   0.16666666666667  0.16666666666667  0.85340590890074 126.904 > 73   92 I   0.50000000000000  0.16666666666667  0.85340590890074 126.904 > 73   93 I   0.83333333333333  0.16666666666667  0.85340590890074 126.904 > 73   94 I   0.16666666666667  0.50000000000000  0.85340590890074 126.904 > 73   95 I   0.50000000000000  0.50000000000000  0.85340590890074 126.904 > 73   96 I   0.83333333333333  0.50000000000000  0.85340590890074 126.904 > 73   97 I   0.16666666666667  0.83333333333333  0.85340590890074 126.904 > 73   98 I   0.50000000000000  0.83333333333333  0.85340590890074 126.904 > 73   99 I   0.83333333333333  0.83333333333333  0.85340590890074 126.904 > 73  100 I   0.00000000000000  0.00000000000000  0.64659409109926 126.904 > 82  101 I   0.33333333333333  0.00000000000000  0.64659409109926 126.904 > 82  102 I   0.66666666666667  0.00000000000000  0.64659409109926 126.904 > 82  103 I   0.00000000000000  0.33333333333333  0.64659409109926 126.904 > 82  104 I   0.33333333333333  0.33333333333333  0.64659409109926 126.904 > 82  105 I   0.66666666666667  0.33333333333333  0.64659409109926 126.904 > 82  106 I   0.00000000000000  0.66666666666667  0.64659409109926 126.904 > 82  107 I   0.33333333333333  0.66666666666667  0.64659409109926 126.904 > 82  108 I   0.66666666666667  0.66666666666667  0.64659409109926 126.904 > 82  109 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 109  110 O   0.33333333333333  0.00000000000000  0.00000000000000  15.999 > 109  111 O   0.66666666666667  0.00000000000000  0.00000000000000  15.999 > 109  112 O   0.00000000000000  0.33333333333333  0.00000000000000  15.999 > 109  113 O   0.33333333333333  0.33333333333333  0.00000000000000  15.999 > 109  114 O   0.66666666666667  0.33333333333333  0.00000000000000  15.999 > 109  115 O   0.00000000000000  0.66666666666667  0.00000000000000  15.999 > 109  116 O   0.33333333333333  0.66666666666667  0.00000000000000  15.999 > 109  117 O   0.66666666666667  0.66666666666667  0.00000000000000  15.999 > 109  118 O   0.16666666666667  0.16666666666667  0.50000000000000  15.999 > 109  119 O   0.50000000000000  0.16666666666667  0.50000000000000  15.999 > 109  120 O   0.83333333333333  0.16666666666667  0.50000000000000  15.999 > 109  121 O   0.16666666666667  0.50000000000000  0.50000000000000  15.999 > 109  122 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 109  123 O   0.83333333333333  0.50000000000000  0.50000000000000  15.999 > 109  124 O   0.16666666666667  0.83333333333333  0.50000000000000  15.999 > 109  125 O   0.50000000000000  0.83333333333333  0.50000000000000  15.999 > 109  126 O   0.83333333333333  0.83333333333333  0.50000000000000  15.999 > 109----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            3.0759119    0.0000000    0.0000000            0.0000000    3.0759119    0.0000000            0.0000000    0.0000000    3.0891469-------------------------- Born effective charges --------------------------    1 K     1.1710065    0.0000000    0.0000000            0.0000000    1.1424984    0.0000000            0.0000000    0.0000000    0.9864279    2 K     1.1290615    0.0000000    0.0000000            0.0000000    1.1290615    0.0000000            0.0000000    0.0000000    1.2349842    3 K     1.1290615    0.0000000    0.0000000            0.0000000    1.1290615    0.0000000            0.0000000    0.0000000    1.2349842    4 K     1.1424984    0.0000000    0.0000000            0.0000000    1.1710065    0.0000000            0.0000000    0.0000000    0.9864279    5 I    -1.1615912    0.0000000    0.0000000            0.0000000   -1.1615912    0.0000000            0.0000000    0.0000000   -1.1299427    6 I    -1.1615912    0.0000000    0.0000000            0.0000000   -1.1615912    0.0000000            0.0000000    0.0000000   -1.1299427    7 O    -2.2484455    0.0000000    0.0000000            0.0000000   -2.2484455    0.0000000            0.0000000    0.0000000   -2.1829389----------------------------------------------------------------------------Sets of supercell forces were read from "FORCES_FC3.xz".Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".----------------------------- Force constants ------------------------------Computing fc3[ 1, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0100  0.0000]    [ 0.0000 -0.0100  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 10, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 73, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 109, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Expanding fc3.Symmetrizing fc3 by traditional approach (N=3).Symmetrizing fc2 by traditional approach (N=3).Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)fc3 was written into "fc3.hdf5".Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency 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-09 00:53:46]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-09 00:53:46]-------------------------Compiled with OpenMP support (max 128 threads).Running in phono3py.load mode.Python version 3.14.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: conductivity-RTAHDF5 data compression filter: gzipCrystal structure was read from "phono3py.yaml".Supercell (dim): [3 3 1]Primitive matrix:  [-0.5  0.5  0.5]  [ 0.5 -0.5  0.5]  [ 0.5  0.5 -0.5]Spacegroup: I4/mmm (139)------------------------------ primitive cell ------------------------------Lattice vectors:  a   -2.642727035000000    2.642727035000000    8.727389710000001  b    2.642727035000000   -2.642727035000000    8.727389710000001  c    2.642727035000000    2.642727035000000   -8.727389710000001Atomic positions (fractional):    1 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098    2 K   0.14793473746599  0.14793473746599  0.00000000000000  39.098    3 K   0.85206526253401  0.85206526253401  0.00000000000000  39.098    4 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098    5 I   0.35340590890074  0.35340590890074  0.00000000000000 126.904    6 I   0.64659409109926  0.64659409109926  0.00000000000000 126.904    7 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999-------------------------------- supercell ---------------------------------Lattice vectors:  a   15.856362210000000    0.000000000000000    0.000000000000000  b    0.000000000000000   15.856362210000000    0.000000000000000  c    0.000000000000000    0.000000000000000   17.454779420000001Atomic positions (fractional):    1 K   0.16666666666667  0.00000000000000  0.00000000000000  39.098 > 1    2 K   0.50000000000000  0.00000000000000  0.00000000000000  39.098 > 1    3 K   0.83333333333333  0.00000000000000  0.00000000000000  39.098 > 1    4 K   0.16666666666667  0.33333333333333  0.00000000000000  39.098 > 1    5 K   0.50000000000000  0.33333333333333  0.00000000000000  39.098 > 1    6 K   0.83333333333333  0.33333333333333  0.00000000000000  39.098 > 1    7 K   0.16666666666667  0.66666666666667  0.00000000000000  39.098 > 1    8 K   0.50000000000000  0.66666666666667  0.00000000000000  39.098 > 1    9 K   0.83333333333333  0.66666666666667  0.00000000000000  39.098 > 1   10 K   0.00000000000000  0.00000000000000  0.14793473746599  39.098 > 10   11 K   0.33333333333333  0.00000000000000  0.14793473746599  39.098 > 10   12 K   0.66666666666667  0.00000000000000  0.14793473746599  39.098 > 10   13 K   0.00000000000000  0.33333333333333  0.14793473746599  39.098 > 10   14 K   0.33333333333333  0.33333333333333  0.14793473746599  39.098 > 10   15 K   0.66666666666667  0.33333333333333  0.14793473746599  39.098 > 10   16 K   0.00000000000000  0.66666666666667  0.14793473746599  39.098 > 10   17 K   0.33333333333333  0.66666666666667  0.14793473746599  39.098 > 10   18 K   0.66666666666667  0.66666666666667  0.14793473746599  39.098 > 10   19 K   0.16666666666667  0.16666666666667  0.35206526253401  39.098 > 19   20 K   0.50000000000000  0.16666666666667  0.35206526253401  39.098 > 19   21 K   0.83333333333333  0.16666666666667  0.35206526253401  39.098 > 19   22 K   0.16666666666667  0.50000000000000  0.35206526253401  39.098 > 19   23 K   0.50000000000000  0.50000000000000  0.35206526253401  39.098 > 19   24 K   0.83333333333333  0.50000000000000  0.35206526253401  39.098 > 19   25 K   0.16666666666667  0.83333333333333  0.35206526253401  39.098 > 19   26 K   0.50000000000000  0.83333333333333  0.35206526253401  39.098 > 19   27 K   0.83333333333333  0.83333333333333  0.35206526253401  39.098 > 19   28 K   0.00000000000000  0.16666666666667  0.00000000000000  39.098 > 28   29 K   0.33333333333333  0.16666666666667  0.00000000000000  39.098 > 28   30 K   0.66666666666667  0.16666666666667  0.00000000000000  39.098 > 28   31 K   0.00000000000000  0.50000000000000  0.00000000000000  39.098 > 28   32 K   0.33333333333333  0.50000000000000  0.00000000000000  39.098 > 28   33 K   0.66666666666667  0.50000000000000  0.00000000000000  39.098 > 28   34 K   0.00000000000000  0.83333333333333  0.00000000000000  39.098 > 28   35 K   0.33333333333333  0.83333333333333  0.00000000000000  39.098 > 28   36 K   0.66666666666667  0.83333333333333  0.00000000000000  39.098 > 28   37 K   0.00000000000000  0.16666666666667  0.50000000000000  39.098 > 1   38 K   0.33333333333333  0.16666666666667  0.50000000000000  39.098 > 1   39 K   0.66666666666667  0.16666666666667  0.50000000000000  39.098 > 1   40 K   0.00000000000000  0.50000000000000  0.50000000000000  39.098 > 1   41 K   0.33333333333333  0.50000000000000  0.50000000000000  39.098 > 1   42 K   0.66666666666667  0.50000000000000  0.50000000000000  39.098 > 1   43 K   0.00000000000000  0.83333333333333  0.50000000000000  39.098 > 1   44 K   0.33333333333333  0.83333333333333  0.50000000000000  39.098 > 1   45 K   0.66666666666667  0.83333333333333  0.50000000000000  39.098 > 1   46 K   0.16666666666667  0.16666666666667  0.64793473746599  39.098 > 10   47 K   0.50000000000000  0.16666666666667  0.64793473746599  39.098 > 10   48 K   0.83333333333333  0.16666666666667  0.64793473746599  39.098 > 10   49 K   0.16666666666667  0.50000000000000  0.64793473746599  39.098 > 10   50 K   0.50000000000000  0.50000000000000  0.64793473746599  39.098 > 10   51 K   0.83333333333333  0.50000000000000  0.64793473746599  39.098 > 10   52 K   0.16666666666667  0.83333333333333  0.64793473746599  39.098 > 10   53 K   0.50000000000000  0.83333333333333  0.64793473746599  39.098 > 10   54 K   0.83333333333333  0.83333333333333  0.64793473746599  39.098 > 10   55 K   0.00000000000000  0.00000000000000  0.85206526253401  39.098 > 19   56 K   0.33333333333333  0.00000000000000  0.85206526253401  39.098 > 19   57 K   0.66666666666667  0.00000000000000  0.85206526253401  39.098 > 19   58 K   0.00000000000000  0.33333333333333  0.85206526253401  39.098 > 19   59 K   0.33333333333333  0.33333333333333  0.85206526253401  39.098 > 19   60 K   0.66666666666667  0.33333333333333  0.85206526253401  39.098 > 19   61 K   0.00000000000000  0.66666666666667  0.85206526253401  39.098 > 19   62 K   0.33333333333333  0.66666666666667  0.85206526253401  39.098 > 19   63 K   0.66666666666667  0.66666666666667  0.85206526253401  39.098 > 19   64 K   0.16666666666667  0.00000000000000  0.50000000000000  39.098 > 28   65 K   0.50000000000000  0.00000000000000  0.50000000000000  39.098 > 28   66 K   0.83333333333333  0.00000000000000  0.50000000000000  39.098 > 28   67 K   0.16666666666667  0.33333333333333  0.50000000000000  39.098 > 28   68 K   0.50000000000000  0.33333333333333  0.50000000000000  39.098 > 28   69 K   0.83333333333333  0.33333333333333  0.50000000000000  39.098 > 28   70 K   0.16666666666667  0.66666666666667  0.50000000000000  39.098 > 28   71 K   0.50000000000000  0.66666666666667  0.50000000000000  39.098 > 28   72 K   0.83333333333333  0.66666666666667  0.50000000000000  39.098 > 28   73 I   0.00000000000000  0.00000000000000  0.35340590890074 126.904 > 73   74 I   0.33333333333333  0.00000000000000  0.35340590890074 126.904 > 73   75 I   0.66666666666667  0.00000000000000  0.35340590890074 126.904 > 73   76 I   0.00000000000000  0.33333333333333  0.35340590890074 126.904 > 73   77 I   0.33333333333333  0.33333333333333  0.35340590890074 126.904 > 73   78 I   0.66666666666667  0.33333333333333  0.35340590890074 126.904 > 73   79 I   0.00000000000000  0.66666666666667  0.35340590890074 126.904 > 73   80 I   0.33333333333333  0.66666666666667  0.35340590890074 126.904 > 73   81 I   0.66666666666667  0.66666666666667  0.35340590890074 126.904 > 73   82 I   0.16666666666667  0.16666666666667  0.14659409109926 126.904 > 82   83 I   0.50000000000000  0.16666666666667  0.14659409109926 126.904 > 82   84 I   0.83333333333333  0.16666666666667  0.14659409109926 126.904 > 82   85 I   0.16666666666667  0.50000000000000  0.14659409109926 126.904 > 82   86 I   0.50000000000000  0.50000000000000  0.14659409109926 126.904 > 82   87 I   0.83333333333333  0.50000000000000  0.14659409109926 126.904 > 82   88 I   0.16666666666667  0.83333333333333  0.14659409109926 126.904 > 82   89 I   0.50000000000000  0.83333333333333  0.14659409109926 126.904 > 82   90 I   0.83333333333333  0.83333333333333  0.14659409109926 126.904 > 82   91 I   0.16666666666667  0.16666666666667  0.85340590890074 126.904 > 73   92 I   0.50000000000000  0.16666666666667  0.85340590890074 126.904 > 73   93 I   0.83333333333333  0.16666666666667  0.85340590890074 126.904 > 73   94 I   0.16666666666667  0.50000000000000  0.85340590890074 126.904 > 73   95 I   0.50000000000000  0.50000000000000  0.85340590890074 126.904 > 73   96 I   0.83333333333333  0.50000000000000  0.85340590890074 126.904 > 73   97 I   0.16666666666667  0.83333333333333  0.85340590890074 126.904 > 73   98 I   0.50000000000000  0.83333333333333  0.85340590890074 126.904 > 73   99 I   0.83333333333333  0.83333333333333  0.85340590890074 126.904 > 73  100 I   0.00000000000000  0.00000000000000  0.64659409109926 126.904 > 82  101 I   0.33333333333333  0.00000000000000  0.64659409109926 126.904 > 82  102 I   0.66666666666667  0.00000000000000  0.64659409109926 126.904 > 82  103 I   0.00000000000000  0.33333333333333  0.64659409109926 126.904 > 82  104 I   0.33333333333333  0.33333333333333  0.64659409109926 126.904 > 82  105 I   0.66666666666667  0.33333333333333  0.64659409109926 126.904 > 82  106 I   0.00000000000000  0.66666666666667  0.64659409109926 126.904 > 82  107 I   0.33333333333333  0.66666666666667  0.64659409109926 126.904 > 82  108 I   0.66666666666667  0.66666666666667  0.64659409109926 126.904 > 82  109 O   0.00000000000000  0.00000000000000  0.00000000000000  15.999 > 109  110 O   0.33333333333333  0.00000000000000  0.00000000000000  15.999 > 109  111 O   0.66666666666667  0.00000000000000  0.00000000000000  15.999 > 109  112 O   0.00000000000000  0.33333333333333  0.00000000000000  15.999 > 109  113 O   0.33333333333333  0.33333333333333  0.00000000000000  15.999 > 109  114 O   0.66666666666667  0.33333333333333  0.00000000000000  15.999 > 109  115 O   0.00000000000000  0.66666666666667  0.00000000000000  15.999 > 109  116 O   0.33333333333333  0.66666666666667  0.00000000000000  15.999 > 109  117 O   0.66666666666667  0.66666666666667  0.00000000000000  15.999 > 109  118 O   0.16666666666667  0.16666666666667  0.50000000000000  15.999 > 109  119 O   0.50000000000000  0.16666666666667  0.50000000000000  15.999 > 109  120 O   0.83333333333333  0.16666666666667  0.50000000000000  15.999 > 109  121 O   0.16666666666667  0.50000000000000  0.50000000000000  15.999 > 109  122 O   0.50000000000000  0.50000000000000  0.50000000000000  15.999 > 109  123 O   0.83333333333333  0.50000000000000  0.50000000000000  15.999 > 109  124 O   0.16666666666667  0.83333333333333  0.50000000000000  15.999 > 109  125 O   0.50000000000000  0.83333333333333  0.50000000000000  15.999 > 109  126 O   0.83333333333333  0.83333333333333  0.50000000000000  15.999 > 109----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            3.0759119    0.0000000    0.0000000            0.0000000    3.0759119    0.0000000            0.0000000    0.0000000    3.0891469-------------------------- Born effective charges --------------------------    1 K     1.1710065    0.0000000    0.0000000            0.0000000    1.1424984    0.0000000            0.0000000    0.0000000    0.9864279    2 K     1.1290615    0.0000000    0.0000000            0.0000000    1.1290615    0.0000000            0.0000000    0.0000000    1.2349842    3 K     1.1290615    0.0000000    0.0000000            0.0000000    1.1290615    0.0000000            0.0000000    0.0000000    1.2349842    4 K     1.1424984    0.0000000    0.0000000            0.0000000    1.1710065    0.0000000            0.0000000    0.0000000    0.9864279    5 I    -1.1615912    0.0000000    0.0000000            0.0000000   -1.1615912    0.0000000            0.0000000    0.0000000   -1.1299427    6 I    -1.1615912    0.0000000    0.0000000            0.0000000   -1.1615912    0.0000000            0.0000000    0.0000000   -1.1299427    7 O    -2.2484455    0.0000000    0.0000000            0.0000000   -2.2484455    0.0000000            0.0000000    0.0000000   -2.1829389----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz)Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330Length for sampling mesh generation: 50.00Generating grid system ... Generalized regular grid: [ 3 9 18 ]Grid generation matrix:  [ 0 9 9 ]  [ 9 0 9 ]  [ 3 3 0 ]fc3-r2q-transformation over three atoms: True--------------------------- Phonon calculations ----------------------------Use NAC by Gonze et al. (no real space sum in current implementation)  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)  G-cutoff distance: 0.66, Number of G-points: 295, Lambda: 0.12Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/60) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 60Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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   1.219   (   0.000    0.000    0.000)    0.000   1.219   (   0.000    0.000    0.000)    0.000   1.840   (   0.000    0.000    0.000)    0.000   1.921   (   0.000    0.000    0.000)    0.000   1.921   (   0.000    0.000    0.000)    0.000   2.275   (   0.000    0.000    0.000)    0.000   2.275   (   0.000    0.000    0.000)    0.000   2.361   (   0.000    0.000    0.000)    0.000   2.361   (   0.000    0.000    0.000)    0.000   2.413   (   0.000    0.000    0.000)    0.000   2.539   (   0.000    0.000    0.000)    0.000   4.119   (   0.000    0.000    0.000)    0.000   4.119   (   0.000    0.000    0.000)    0.000   4.329   (   0.000    0.000    0.000)    0.000   4.510   (   0.000    0.000    0.000)    0.000   7.256   (   0.000    0.000    0.000)    0.000   7.256   (   0.000    0.000    0.000)    0.000   8.488   (   0.000    0.000    0.000)    0.000======================= Grid point 1 (2/60) =======================q-point: ( 0.00  0.33  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.931   (  12.071    0.000    6.191)   13.567   0.992   (  11.121    0.000    1.391)   11.207   1.275   (   1.421    0.000   -0.477)    1.499   1.402   (   2.924    0.000   -0.150)    2.928   1.462   (   8.914    0.000    4.125)    9.822   1.640   (  -4.114   -0.000   -2.964)    5.071   1.781   (  -4.045   -0.000    0.745)    4.113   1.895   (  -4.063   -0.000    0.707)    4.124   2.295   (   0.478    0.000   -2.995)    3.033   2.309   (  -0.840   -0.000   -0.647)    1.060   2.410   (   0.725    0.000    0.991)    1.228   2.575   (   6.651    0.000    0.376)    6.661   3.104   (  11.386    0.000   -7.455)   13.610   3.360   (  10.502    0.000    4.881)   11.581   3.950   (  -7.021   -0.000   -3.301)    7.758   3.951   (  -3.958   -0.000    0.039)    3.958   4.174   (  -5.643   -0.000    3.156)    6.466   5.376   (  15.977    0.000   -1.402)   16.038   7.308   (   1.024    0.000    0.031)    1.025   8.331   (  -2.724   -0.000    5.207)    5.877   9.013   (   6.927    0.000   -2.878)    7.501======================= Grid point 3 (3/60) =======================q-point: ( 0.00  0.00  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.354   (   9.123    9.123    0.000)   12.902   0.569   (  11.688   11.688    0.000)   16.529   0.849   (  16.788   16.788    0.000)   23.742   1.256   (   1.710    1.710    0.000)    2.418   1.281   (   2.355    2.355    0.000)    3.330   1.783   (  -2.134   -2.134    0.000)    3.019   2.002   (   3.184    3.184    0.000)    4.502   2.148   (  -6.159   -6.159    0.000)    8.710   2.248   (  -1.352   -1.352    0.000)    1.912   2.254   (  -3.452   -3.452    0.000)    4.882   2.302   (   2.917    2.917    0.000)    4.125   2.533   (   4.078    4.078    0.000)    5.768   2.680   (  11.118   11.118    0.000)   15.724   2.765   (   6.781    6.781    0.000)    9.589   4.223   (   4.034    4.034    0.000)    5.704   4.254   (  -5.576   -5.576    0.000)    7.886   4.273   (  -2.700   -2.700    0.000)    3.818   4.826   (   8.046    8.046    0.000)   11.378   7.353   (   5.096    5.096    0.000)    7.207   8.416   (  -2.817   -2.817    0.000)    3.983   8.689   (   1.397    1.397    0.000)    1.976======================= Grid point 4 (4/60) =======================q-point: ( 0.00  0.33  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.173   (   7.577    0.753    0.747)    7.651   1.211   (   6.963    1.859    1.968)    7.471   1.322   (   0.737    2.159   -1.009)    2.494   1.382   (  -4.095    0.654    1.129)    4.298   1.532   (  -1.567   -0.858   -0.012)    1.787   1.614   (   2.003   -0.935   -1.869)    2.895   1.817   (  -1.002    1.870    0.591)    2.202   1.856   (  -1.070    6.951    1.311)    7.154   2.078   (  -1.597   -9.614   -1.185)    9.817   2.276   (   0.965   -9.590   -0.023)    9.639   2.533   (  -0.141    8.729   -0.842)    8.771   2.685   (   1.527   -0.125    0.455)    1.598   3.340   (   3.726    6.833   -2.514)    8.179   3.531   (   7.796   -1.438    5.282)    9.526   3.791   (  -6.336   -2.675   -6.269)    9.306   4.032   (  -2.258    3.700    2.319)    4.916   4.061   (  -2.379    6.235    0.627)    6.703   5.566   (   5.342   -2.445   -0.337)    5.885   7.537   (   0.993   18.350   -0.017)   18.377   8.307   (   0.196    1.275    1.498)    1.977   9.069   (   2.824   -4.883   -0.513)    5.664======================= Grid point 5 (5/60) =======================q-point: ( 0.00 -0.33  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.741   (  -8.506    3.702  -12.501)   15.567   0.837   ( -10.711    7.943   -2.962)   13.660   1.208   ( -12.303    4.829   -1.034)   13.257   1.276   (  -2.702    2.704    0.977)    3.945   1.368   (  -3.986    1.227   -2.201)    4.716   1.702   (   2.752   -1.182    2.488)    3.894   1.961   (   4.994    2.255    1.071)    5.583   2.007   (   5.080    4.656    0.099)    6.891   2.159   (   1.548   -6.355    3.618)    7.475   2.223   (  -0.318   -4.590   -0.201)    4.605   2.392   (  -8.352   -2.968    3.713)    9.610   2.609   (  -1.132   10.435   -2.090)   10.703   2.899   ( -12.018    7.891    4.604)   15.096   3.063   ( -15.121   -0.266   -2.306)   15.298   4.077   (   6.972   -3.010    1.664)    7.774   4.209   (   3.380   11.061    0.233)   11.568   4.267   (   5.620   -2.447   -4.201)    7.431   5.010   ( -16.875    2.785    5.129)   17.855   7.438   (  -3.266   12.819   -0.031)   13.229   8.308   (  -7.475    1.058    1.818)    7.765   8.862   (   2.519   -6.635   -4.164)    8.228======================= Grid point 6 (6/60) =======================q-point: ( 0.00  0.00  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.754   (   9.548    9.548    0.000)   13.503   0.969   (   7.508    7.508    0.000)   10.617   1.347   (   2.327    2.327    0.000)    3.291   1.349   (   7.076    7.076    0.000)   10.007   1.370   (   1.659    1.659    0.000)    2.347   1.697   (  -1.738   -1.738    0.000)    2.458   1.915   (  -4.711   -4.711    0.000)    6.662   2.084   (  -1.344   -1.344    0.000)    1.900   2.084   (  -4.624   -4.624    0.000)    6.539   2.199   (   0.537    0.537    0.000)    0.760   2.416   (   1.973    1.973    0.000)    2.790   2.731   (   4.991    4.991    0.000)    7.058   3.031   (   5.532    5.532    0.000)    7.823   3.240   (  13.246   13.246    0.000)   18.732   4.057   (  -3.723   -3.723    0.000)    5.265   4.090   (  -6.139   -6.139    0.000)    8.681   4.350   (   1.003    1.003    0.000)    1.418   5.220   (  11.071   11.071    0.000)   15.657   7.721   (  12.234   12.234    0.000)   17.301   8.319   (  -1.107   -1.107    0.000)    1.566   8.748   (   1.159    1.159    0.000)    1.639======================= Grid point 7 (7/60) =======================q-point: ( 0.00 -0.67  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.243   (  -2.816    2.896   -1.938)    4.480   1.299   (  -4.976    2.728   -1.552)    5.883   1.330   (   0.483    1.801    2.074)    2.789   1.421   (   2.049    0.282   -1.698)    2.676   1.503   (   1.479   -0.040    1.850)    2.369   1.592   (  -0.565   -1.391   -2.065)    2.553   1.779   (  -0.484   -3.243    0.931)    3.409   1.870   (   0.505   -7.321    1.584)    7.507   2.053   (   0.853  -10.056    0.078)   10.093   2.075   (   1.536    9.089   -0.056)    9.218   2.549   (  -2.998   -4.576   -0.375)    5.483   2.806   (   1.652    9.937    0.326)   10.079   3.433   (  -8.565    1.700   -1.519)    8.863   3.571   (  -6.569    0.533    1.304)    6.719   3.631   (   7.944   -4.539    2.555)    9.499   4.010   (   3.685   -2.334   -1.118)    4.503   4.259   (   1.606    4.448   -0.062)    4.729   5.628   (  -5.509   10.598    0.441)   11.953   7.997   (  -1.859   23.360    0.023)   23.434   8.377   (   2.601    1.444   -2.077)    3.628   8.933   (  -4.372   -5.757    0.912)    7.286======================= Grid point 9 (8/60) =======================q-point: ( 0.00  0.00  0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.124   (   7.633    7.633    0.000)   10.795   1.212   (   4.178    4.178    0.000)    5.909   1.417   (   0.594    0.594    0.000)    0.840   1.424   (   1.093    1.093    0.000)    1.545   1.490   (   0.561    0.561    0.000)    0.793   1.642   (  -1.056   -1.056    0.000)    1.494   1.763   (  -2.432   -2.432    0.000)    3.439   1.869   (  -5.442   -5.442    0.000)    7.696   1.922   (  -5.155   -5.155    0.000)    7.290   2.256   (   1.403    1.403    0.000)    1.985   2.444   (  -0.610   -0.610    0.000)    0.862   2.928   (   4.366    4.366    0.000)    6.174   3.210   (   2.859    2.859    0.000)    4.044   3.650   (   5.511    5.511    0.000)    7.794   3.768   (  -8.764   -8.764    0.000)   12.394   3.952   (  -0.996   -0.996    0.000)    1.408   4.300   (  -2.731   -2.731    0.000)    3.862   5.743   (  12.811   12.811    0.000)   18.118   8.293   (  13.387   13.387    0.000)   18.932   8.341   (   1.803    1.803    0.000)    2.550   8.774   (   0.065    0.065    0.000)    0.092======================= Grid point 10 (9/60) =======================q-point: ( 0.00 -0.67  0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.247   (  -5.039    5.039   -3.958)    8.152   1.262   (  -3.287    3.287   -2.443)    5.251   1.370   (  -2.059    2.059    2.558)    3.877   1.391   (  -0.370    0.370    0.593)    0.791   1.543   (   0.699   -0.699   -2.204)    2.415   1.599   (   0.522   -0.522    0.410)    0.844   1.688   (   2.786   -2.786    2.058)    4.445   1.818   (   3.390   -3.390    2.858)    5.581   1.931   (   5.719   -5.719   -0.072)    8.088   2.316   (  -1.719    1.719   -1.786)    3.017   2.342   (  -0.851    0.851    2.654)    2.914   3.045   (  -2.705    2.705   -3.943)    5.494   3.187   (  -2.767    2.767    0.738)    3.983   3.551   (  -7.786    7.786    5.663)   12.382   3.669   (   7.205   -7.205    1.602)   10.314   4.041   (   3.852   -3.852   -2.619)    6.044   4.298   (   2.735   -2.735    0.074)    3.868   5.713   ( -14.094   14.094    0.884)   19.951   8.294   ( -13.119   13.119   -0.039)   18.553   8.425   (   0.771   -0.771   -3.459)    3.627   8.757   (  -1.448    1.448    1.311)    2.432======================= Grid point 12 (10/60) =======================q-point: ( 0.00  0.00  0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.334   (   1.674    1.674    0.000)    2.368   1.345   (   1.805    1.805    0.000)    2.552   1.441   (  -0.007   -0.007    0.000)    0.010   1.460   (   2.751    2.751    0.000)    3.891   1.477   (  -0.509   -0.509    0.000)    0.720   1.601   (  -0.679   -0.679    0.000)    0.960   1.652   (  -4.497   -4.497    0.000)    6.360   1.704   (  -0.586   -0.586    0.000)    0.829   1.719   (  -3.453   -3.453    0.000)    4.883   2.304   (   0.719    0.719    0.000)    1.017   2.397   (  -0.984   -0.984    0.000)    1.391   3.122   (   5.274    5.274    0.000)    7.459   3.280   (   0.683    0.683    0.000)    0.965   3.404   (  -7.924   -7.924    0.000)   11.207   3.736   (   0.083    0.083    0.000)    0.118   3.967   (   0.796    0.796    0.000)    1.125   4.189   (  -1.760   -1.760    0.000)    2.489   6.169   (   5.933    5.933    0.000)    8.390   8.419   (   1.283    1.283    0.000)    1.815   8.718   (   5.674    5.674    0.000)    8.024   8.764   (  -0.323   -0.323    0.000)    0.456======================= Grid point 27 (11/60) =======================q-point: ( 0.06 -0.06 -0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.249   ( -12.874   -0.000    0.000)   12.874   0.361   ( -16.756   -0.000    0.000)   16.756   0.649   ( -27.010   -0.000    0.000)   27.010   1.228   (  -0.835   -0.000    0.000)    0.835   1.263   (  -3.746   -0.000    0.000)    3.746   1.809   (   2.499    0.000   -0.000)    2.499   1.903   (   1.728    0.000   -0.000)    1.728   2.208   (   9.000    0.000   -0.000)    9.000   2.249   (  -2.011   -0.000    0.000)    2.011   2.316   (  -3.782   -0.000    0.000)    3.782   2.355   (   0.594    0.000   -0.000)    0.594   2.423   (   0.565    0.000   -0.000)    0.565   2.625   ( -16.208   -0.000    0.000)   16.208   2.701   ( -11.309   -0.000    0.000)   11.309   4.097   (   2.108    0.000   -0.000)    2.108   4.302   (   2.466    0.000   -0.000)    2.466   4.325   (   8.759    0.000   -0.000)    8.759   4.743   ( -12.358   -0.000    0.000)   12.358   7.264   (  -0.724   -0.000    0.000)    0.724   8.448   (   3.498    0.000   -0.000)    3.498   8.718   (  -5.679   -0.000    0.000)    5.679======================= Grid point 28 (12/60) =======================q-point: ( 0.06  0.28 -0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.703   (   9.123    0.000   13.500)   16.293   0.723   (  14.105    0.000    2.600)   14.343   1.156   (  16.877    0.000    3.528)   17.242   1.241   (   1.600    0.000   -0.976)    1.874   1.356   (   4.320    0.000    1.226)    4.491   1.718   (  -3.107   -0.000   -2.627)    4.069   1.860   (  -3.256   -0.000    0.940)    3.389   2.013   (  -6.861   -0.000   -1.722)    7.073   2.257   (   3.021    0.000   -2.873)    4.169   2.327   (  -0.797   -0.000   -1.003)    1.281   2.407   (  -0.089   -0.000   -0.041)    0.098   2.431   (   6.502    0.000    0.353)    6.511   2.819   (  14.856    0.000   -8.553)   17.142   3.090   (  14.658    0.000    3.184)   15.000   4.034   (  -3.679   -0.000    0.020)    3.679   4.110   (  -7.741   -0.000   -1.403)    7.867   4.296   (  -5.879   -0.000    4.369)    7.325   4.996   (  19.044    0.000   -4.427)   19.551   7.284   (   1.140    0.000    0.029)    1.140   8.297   (  10.163    0.000   -3.342)   10.699   8.950   (  -5.195   -0.000    4.884)    7.130======================= Grid point 29 (13/60) =======================q-point: ( 0.06 -0.39 -0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.171   (  -5.170   -0.000   -0.445)    5.189   1.187   ( -11.582   -0.000   -1.505)   11.679   1.296   (  -0.471   -0.000    0.138)    0.491   1.375   (   4.736    0.000   -1.356)    4.926   1.545   (   2.304    0.000    0.952)    2.493   1.631   (  -3.891   -0.000    0.927)    4.000   1.708   (   2.245    0.000   -0.298)    2.264   1.847   (   0.737    0.000   -1.806)    1.951   2.288   (   0.736   -0.000    2.244)    2.361   2.294   (   0.490    0.000    0.229)    0.541   2.429   (  -0.691   -0.000   -1.048)    1.255   2.684   (  -3.073   -0.000   -0.170)    3.077   3.256   (  -1.961   -0.000    2.935)    3.530   3.549   (  -8.416   -0.000   -5.624)   10.122   3.819   (   5.816    0.000    6.790)    8.940   3.885   (   1.845    0.000   -0.015)    1.845   4.077   (   2.730    0.000   -3.232)    4.231   5.618   (  -6.089   -0.000    0.242)    6.094   7.324   (  -0.419   -0.000   -0.014)    0.419   8.293   (   0.014    0.000   -1.614)    1.614   9.124   (  -2.495   -0.000    0.701)    2.591======================= Grid point 31 (14/60) =======================q-point: ( 0.06  0.28  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.952   (  11.089    2.315    5.803)   12.728   1.051   (   9.173    4.292    2.077)   10.338   1.299   (   1.506    2.176   -1.015)    2.834   1.403   (   2.537    0.070   -0.344)    2.561   1.473   (   6.840    0.911    3.642)    7.803   1.634   (  -3.521   -0.344   -2.988)    4.631   1.863   (  -3.544   -0.396   -0.101)    3.567   1.912   (  -3.962    8.239    0.967)    9.194   2.118   (  -2.071   -8.257   -2.594)    8.899   2.247   (   1.584   -8.042    0.445)    8.209   2.518   (   2.766    3.743   -2.138)    5.122   2.631   (   2.305    5.622    1.122)    6.179   3.159   (  11.606    4.886   -4.495)   13.371   3.346   (  10.997    0.269    3.219)   11.462   3.927   (  -6.928   -2.124   -3.449)    8.025   4.119   (  -4.822    9.176    1.172)   10.432   4.146   (  -5.430    0.402    1.998)    5.800   5.351   (  14.323   -0.473   -1.728)   14.435   7.499   (   2.481   16.591   -0.022)   16.775   8.334   (  -2.038    0.271    4.803)    5.225   8.950   (   7.368   -5.348   -2.237)    9.375======================= Grid point 33 (15/60) =======================q-point: ( 0.06 -0.06  0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.597   (   5.105   13.289    0.000)   14.236   0.807   (   7.866   11.232    0.000)   13.712   1.194   (   6.674   14.327    0.000)   15.805   1.284   (   2.794    1.249   -0.000)    3.060   1.354   (   0.372    3.638    0.000)    3.657   1.729   (  -1.312   -2.862   -0.000)    3.148   1.978   (   9.758   -3.554   -0.000)   10.385   1.997   (  -2.775   -7.240   -0.000)    7.753   2.199   (  -6.598   -2.451    0.000)    7.038   2.220   (  -4.504    0.130    0.000)    4.505   2.453   (  -1.540    7.999    0.000)    8.146   2.571   (   9.279    0.265   -0.000)    9.283   2.969   (   1.194   11.530    0.000)   11.592   3.010   (   5.662   16.794    0.000)   17.723   4.118   (  -1.581   -6.527   -0.000)    6.715   4.191   (  -3.357   -5.073   -0.000)    6.083   4.222   (  10.527   -3.070   -0.000)   10.965   5.067   (   3.114   14.102    0.000)   14.441   7.438   (  12.749    3.316   -0.000)   13.173   8.345   (  -1.144   -3.499   -0.000)    3.681   8.799   (  -4.959    7.820    0.000)    9.260======================= Grid point 34 (16/60) =======================q-point: ( 0.06  0.28  0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.211   (   6.436    2.891    1.235)    7.163   1.258   (   4.977    2.470    2.077)    5.932   1.371   (   0.390    2.200   -1.920)    2.946   1.400   (  -2.537    1.009    0.639)    2.804   1.514   (  -0.906   -0.784    0.839)    1.462   1.597   (  -0.276   -0.670   -2.570)    2.671   1.781   (  -0.083   -2.756    0.667)    2.836   1.898   (  -2.854   -7.455   -1.091)    8.057   2.046   (  -0.254   -7.306    0.783)    7.352   2.077   (  -1.613    6.542   -0.196)    6.741   2.562   (   1.743   -4.498   -1.696)    5.113   2.796   (  -0.547    9.301    1.388)    9.420   3.414   (   9.568   -2.387    0.558)    9.877   3.562   (  -0.133    5.322    1.011)    5.419   3.709   (  -6.355   -4.318   -4.004)    8.664   3.982   (  -1.473   -2.706    1.638)    3.490   4.258   (  -1.560    4.429    0.027)    4.696   5.637   (   4.561   10.499   -0.432)   11.454   7.997   (   1.827   23.361   -0.036)   23.432   8.343   (   0.665    2.096    1.197)    2.503   8.944   (   3.202   -6.303   -0.155)    7.071======================= Grid point 35 (17/60) =======================q-point: ( 0.06 -0.39  0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.857   (  -7.226    7.226   -9.432)   13.907   0.996   (  -7.154    7.154   -2.617)   10.450   1.325   (  -4.760    4.760    1.325)    6.861   1.327   (  -2.339    2.339    2.002)    3.866   1.418   (  -3.547    3.547   -3.869)    6.335   1.672   (   1.639   -1.639    2.192)    3.190   1.890   (   4.691   -4.691    2.086)    6.954   2.050   (   4.092   -4.092    3.528)    6.777   2.107   (   1.943   -1.943   -1.879)    3.329   2.205   (  -1.051    1.051   -0.652)    1.623   2.343   (  -2.740    2.740    5.612)    6.820   2.789   (  -5.075    5.075   -4.498)    8.470   3.029   (  -5.376    5.376    0.232)    7.606   3.211   ( -12.764   12.764    3.229)   18.338   3.986   (   5.536   -5.536    2.767)    8.303   4.180   (   5.106   -5.106   -4.682)    8.606   4.349   (  -1.020    1.020    0.084)    1.445   5.173   ( -12.386   12.386    4.288)   18.033   7.723   ( -12.236   12.236   -0.204)   17.306   8.338   (  -1.306    1.306   -2.412)    3.038   8.761   (   1.171   -1.171   -0.553)    1.746======================= Grid point 36 (18/60) =======================q-point: ( 0.06 -0.06  0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.978   (   5.823   11.158    0.000)   12.586   1.121   (   4.398    6.675    0.000)    7.994   1.373   (   2.975    0.823   -0.000)    3.086   1.413   (   0.334    1.174    0.000)    1.221   1.467   (   1.440    3.835    0.000)    4.097   1.659   (  -0.597   -1.673   -0.000)    1.777   1.819   (  -2.618   -3.878   -0.000)    4.679   1.992   (  -6.408   -4.206    0.000)    7.665   2.049   (  -4.880   -1.196    0.000)    5.024   2.139   (   6.796   -3.941   -0.000)    7.856   2.521   (  -5.212    5.154    0.000)    7.330   2.777   (   9.311    0.212   -0.000)    9.313   3.199   (  -0.948    9.613    0.000)    9.659   3.477   (  10.016    8.059   -0.000)   12.856   3.939   (  -7.280   -8.125   -0.000)   10.909   3.981   (  -1.373   -2.868   -0.000)    3.180   4.314   (   3.169   -3.094   -0.000)    4.429   5.472   (  10.569   11.262    0.000)   15.445   7.915   (  20.812    6.302   -0.000)   21.745   8.308   (   1.267    0.224   -0.000)    1.287   8.838   (  -5.475    6.106    0.000)    8.202======================= Grid point 37 (19/60) =======================q-point: ( 0.06 -0.72  0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.307   (  -1.179    2.854   -1.388)    3.385   1.349   (  -2.300    1.633   -1.284)    3.100   1.376   (  -0.978    2.529    1.241)    2.982   1.424   (  -0.735    0.262   -0.388)    0.872   1.509   (   2.237    0.706    0.831)    2.488   1.574   (   0.399   -0.284   -1.256)    1.348   1.691   (  -2.212   -4.634    1.346)    5.308   1.756   (   2.379   -3.809    0.821)    4.566   1.848   (   2.093   -8.288    0.011)    8.548   2.253   (   2.109    6.711    0.309)    7.041   2.437   (  -2.848   -5.319   -0.078)    6.034   3.008   (   2.328    9.012   -0.768)    9.339   3.345   (  -8.729   -6.151   -1.446)   10.776   3.502   (   7.669   -5.922    1.044)    9.746   3.689   (  -4.524    5.528    5.193)    8.832   3.967   (   2.885   -1.820   -3.012)    4.551   4.249   (   1.349   -3.278    0.011)    3.545   5.930   (  -5.667   15.058    0.247)   16.091   8.408   (   0.798    1.522   -1.159)    2.073   8.448   (  -2.590   17.939    0.083)   18.125   8.826   (  -3.073   -4.077    0.222)    5.111======================= Grid point 39 (20/60) =======================q-point: ( 0.06 -0.06  0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.275   (   4.131    5.954    0.000)    7.246   1.291   (   2.432    3.170    0.000)    3.996   1.420   (   1.171   -0.801   -0.000)    1.419   1.433   (   1.225    0.611   -0.000)    1.369   1.484   (  -0.460   -1.324   -0.000)    1.401   1.617   (  -0.851   -1.170   -0.000)    1.447   1.726   (  -1.760   -1.249    0.000)    2.159   1.765   (  -5.928   -4.052    0.000)    7.181   1.846   (  -8.195   -1.919    0.000)    8.417   2.236   (   5.937   -1.133   -0.000)    6.044   2.465   (  -5.308    0.764    0.000)    5.362   2.972   (   8.531    1.062   -0.000)    8.596   3.324   (  -3.882    7.136    0.000)    8.124   3.577   (  -8.412   -9.302   -0.000)   12.541   3.691   (   4.466   -1.217   -0.000)    4.629   3.956   (   0.088    1.412    0.000)    1.415   4.250   (  -3.187   -1.426    0.000)    3.492   5.939   (  14.757    4.833   -0.000)   15.528   8.381   (   2.203    1.948   -0.000)    2.941   8.449   (  18.275    2.540   -0.000)   18.451   8.825   (  -4.548    2.809    0.000)    5.345======================= Grid point 42 (21/60) =======================q-point: ( 0.06  0.94 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.353   (   0.539   -0.539   -0.000)    0.762   1.361   (  -0.751    0.751    0.000)    1.063   1.422   (   1.752   -1.752   -0.000)    2.478   1.447   (   2.040   -2.040   -0.000)    2.885   1.522   (  -0.886    0.886    0.000)    1.253   1.592   (   0.078   -0.078   -0.000)    0.111   1.593   (  -0.730    0.730    0.000)    1.032   1.706   (  -0.858    0.858    0.000)    1.213   1.719   (  -3.451    3.451    0.000)    4.881   2.317   (  -0.442    0.442    0.000)    0.625   2.381   (   0.512   -0.512   -0.000)    0.724   3.179   (   2.714   -2.714   -0.000)    3.838   3.257   (   0.951   -0.951   -0.000)    1.345   3.369   (  -5.584    5.584    0.000)    7.897   3.710   (   2.501   -2.501   -0.000)    3.537   3.983   (  -0.730    0.730    0.000)    1.032   4.189   (  -1.682    1.682    0.000)    2.378   6.166   (   6.257   -6.257   -0.000)    8.849   8.432   (   0.009   -0.009   -0.000)    0.013   8.703   (   5.373   -5.373   -0.000)    7.599   8.779   (   0.007   -0.007   -0.000)    0.010======================= Grid point 45 (22/60) =======================q-point: ( 0.06  0.61 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.196   (   1.492   -4.017    1.424)    4.516   1.246   (   2.149   -6.512    1.170)    6.957   1.302   (   0.894   -0.125   -0.903)    1.276   1.411   (   0.520    5.110    1.822)    5.450   1.506   (  -0.236    0.339   -1.967)    2.010   1.627   (  -1.707   -1.175    2.555)    3.289   1.820   (   1.557    1.083   -0.436)    1.946   1.861   (   6.966    1.426   -0.769)    7.152   2.056   ( -10.158   -0.419   -0.863)   10.203   2.282   (  -9.547   -0.435    0.568)    9.574   2.534   (   8.771    0.170    0.950)    8.824   2.675   (  -0.284   -2.335   -1.428)    2.751   3.288   (   8.066   -7.405   -2.533)   11.239   3.617   (  -2.406   -8.495    1.695)    8.990   3.722   (  -3.270    8.032    0.730)    8.703   4.048   (  12.434    2.410    0.155)   12.667   4.061   (  -1.811    3.312   -1.170)    3.952   5.557   (  -2.558   -6.221   -0.461)    6.742   7.537   (  18.315   -1.015   -0.012)   18.343   8.353   (   0.826    4.224    2.940)    5.212   9.046   (  -4.397   -5.125   -1.679)    6.958======================= Grid point 48 (23/60) =======================q-point: ( 0.06  0.61 -0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.141   (   4.539   -5.326    5.073)    8.643   1.191   (   3.792   -5.720    2.876)    7.441   1.316   (   2.444   -1.488   -2.428)    3.753   1.392   (  -0.043   -3.447   -0.639)    3.506   1.558   (  -0.783   -0.651    3.879)    4.010   1.589   (   0.728    2.486   -1.039)    2.791   1.781   (  -4.815    1.464   -1.280)    5.193   1.903   (  -5.188    2.162   -4.279)    7.064   2.081   (  -8.289    1.954    0.560)    8.535   2.137   (  10.370    4.019   -0.337)   11.126   2.436   (  -4.863   -7.643   -0.934)    9.107   2.862   (  11.380    3.234    1.589)   11.936   3.220   (   0.102  -10.984    1.630)   11.104   3.377   (   6.031  -11.268   -5.541)   13.930   3.810   (  -5.874    8.168   -0.568)   10.076   4.117   (  -3.211    6.031    1.515)    6.999   4.312   (   3.389    2.897    0.031)    4.458   5.415   (  11.390  -14.089   -1.653)   18.193   7.916   (  20.982   -6.013    0.061)   21.827   8.444   (  -1.029    2.941    6.096)    6.847   8.786   (  -3.536   -7.869   -3.420)    9.280======================= Grid point 53 (24/60) =======================q-point: ( 0.06  0.28 -0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.603   (   4.847   -4.847   18.020)   19.280   0.612   (  11.095  -11.095    4.087)   16.214   0.898   (  14.330  -14.330    5.456)   20.987   1.238   (   1.659   -1.659   -1.806)    2.961   1.289   (   2.967   -2.967    1.084)    4.333   1.750   (  -1.826    1.826   -3.085)    4.023   2.010   (   2.760   -2.760    0.778)    3.980   2.136   (  -5.595    5.595   -1.147)    7.995   2.143   (   2.603   -2.603   -8.197)    8.986   2.253   (  -1.319    1.319    0.470)    1.923   2.257   (  -1.677    1.677   -3.501)    4.229   2.501   (   8.905   -8.905   -1.122)   12.644   2.758   (   5.464   -5.464    2.019)    7.986   2.766   (   6.843   -6.843    0.033)    9.678   4.200   (  -4.290    4.290   -1.715)    6.304   4.221   (   4.024   -4.024   -0.136)    5.692   4.365   (  -3.156    3.156    4.967)    6.678   4.704   (  10.015  -10.015   -8.778)   16.663   7.354   (   5.133   -5.133    0.104)    7.260   8.031   (  15.318  -15.318  -16.193)   27.046   9.072   ( -14.764   14.764   16.082)   26.355======================= Grid point 55 (25/60) =======================q-point: ( 0.11 -0.11 -0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.544   ( -14.954   -0.000    0.000)   14.954   0.696   ( -14.868   -0.000    0.000)   14.868   1.117   ( -17.909   -0.000    0.000)   17.909   1.251   (  -1.338   -0.000    0.000)    1.338   1.350   (  -3.947   -0.000    0.000)    3.947   1.747   (   3.299    0.000   -0.000)    3.299   1.850   (   3.221    0.000   -0.000)    3.221   2.026   (   8.100    0.000   -0.000)    8.100   2.303   (  -2.960   -0.000    0.000)    2.960   2.337   (   0.985    0.000   -0.000)    0.985   2.401   (   1.434    0.000   -0.000)    1.434   2.427   (  -6.483   -0.000    0.000)    6.483   2.958   ( -12.879   -0.000    0.000)   12.879   3.000   ( -16.713   -0.000    0.000)   16.713   4.034   (   3.702    0.000   -0.000)    3.702   4.150   (   7.570    0.000   -0.000)    7.570   4.227   (   4.646    0.000   -0.000)    4.646   5.047   ( -16.294   -0.000    0.000)   16.294   7.284   (  -1.134   -0.000    0.000)    1.134   8.358   (   4.494    0.000   -0.000)    4.494   8.872   (  -8.308   -0.000    0.000)    8.308======================= Grid point 56 (26/60) =======================q-point: ( 0.11  0.22 -0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.415   (  14.539    0.000    4.899)   15.342   0.543   (   9.080    0.000   16.066)   18.455   0.746   (  19.114    0.000   12.733)   22.966   1.213   (   1.026    0.000   -1.465)    1.789   1.259   (   4.412    0.000   -0.227)    4.417   1.773   (  -2.093   -0.000   -3.399)    3.992   1.913   (  -1.706   -0.000    0.946)    1.951   2.122   (   5.612    0.000   -7.049)    9.010   2.205   (  -5.305   -0.000   -1.964)    5.657   2.314   (   3.465    0.000   -0.165)    3.469   2.345   (  -0.130    0.000   -0.867)    0.876   2.349   (   6.532    0.000   -1.584)    6.721   2.596   (   2.553    0.000   -3.260)    4.141   2.800   (  10.802    0.000    1.544)   10.911   4.096   (  -2.080   -0.000   -0.037)    2.081   4.246   (  -4.580   -0.000   -1.973)    4.987   4.395   (   1.503    0.000    1.319)    1.999   4.623   (  12.238    0.000   -3.939)   12.856   7.264   (   0.732    0.000    0.014)    0.732   7.839   (  32.588    0.000  -19.302)   37.875   9.289   ( -25.834   -0.000   16.932)   30.889======================= Grid point 57 (27/60) =======================q-point: ( 0.11  0.56 -0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.180   (   4.242    0.000    0.473)    4.268   1.220   (   7.979    0.000    1.573)    8.133   1.293   (   0.786    0.000   -0.164)    0.803   1.406   (  -6.674   -0.000    1.519)    6.845   1.510   (   0.176    0.000   -1.752)    1.761   1.655   (   1.722    0.000    2.798)    3.285   1.714   (  -2.814   -0.000    0.297)    2.830   1.851   (  -0.863   -0.000   -1.496)    1.727   2.276   (   0.286   -0.000    1.037)    1.076   2.290   (  -0.036    0.000   -0.229)    0.232   2.429   (   0.751    0.000   -0.974)    1.230   2.687   (   2.761    0.000    0.168)    2.766   3.188   (   7.551    0.000   -3.568)    8.352   3.644   (   7.269    0.000    1.545)    7.431   3.757   (  -7.139   -0.000    2.327)    7.509   3.885   (  -1.876   -0.000    0.015)    1.876   4.086   (  -3.386   -0.000   -2.001)    3.933   5.611   (   6.821    0.000   -0.422)    6.834   7.324   (   0.393    0.000    0.014)    0.393   8.344   (  -4.917   -0.000    3.310)    5.927   9.096   (   5.359    0.000   -2.041)    5.734======================= Grid point 62 (28/60) =======================q-point: ( 0.11  0.22  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.033   (   9.070    5.266    5.240)   11.724   1.140   (   6.440    4.115    1.998)    7.899   1.356   (   1.013    2.900   -1.610)    3.468   1.405   (   2.060    0.146   -0.677)    2.173   1.494   (   2.929    0.920    2.542)    3.986   1.629   (  -2.399   -0.101   -2.550)    3.503   1.806   (  -2.705   -3.603   -1.213)    4.666   1.970   (  -3.713   -5.952   -2.399)    7.414   2.064   (  -1.485   -6.357    1.171)    6.632   2.140   (  -4.084    8.481    0.030)    9.413   2.477   (   6.190   -4.641   -3.175)    8.363   2.816   (  -0.900    9.433    2.949)    9.924   3.203   (  10.000   -0.977    0.420)   10.056   3.457   (   9.277    9.162   -2.439)   13.264   3.853   (  -6.871   -5.004   -4.173)    9.469   4.063   (  -5.439   -3.965    4.141)    7.902   4.313   (  -3.041    3.317   -0.087)    4.501   5.452   (  12.290   10.842   -1.837)   16.491   7.915   (   6.200   20.872    0.028)   21.774   8.344   (  -0.417    0.744    3.781)    3.876   8.830   (   6.251   -5.060   -1.117)    8.119======================= Grid point 64 (29/60) =======================q-point: ( 0.11 -0.11  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.890   (   2.525   14.098    0.000)   14.322   1.030   (   4.058    9.544    0.000)   10.371   1.310   (   2.610    1.171   -0.000)    2.861   1.407   (   0.232    0.811    0.000)    0.844   1.436   (   1.335    8.505    0.000)    8.609   1.670   (  -0.457   -2.495   -0.000)    2.537   1.864   (  -1.530   -4.956   -0.000)    5.187   1.896   (   9.998   -3.781   -0.000)   10.688   2.145   (  -8.195   -2.644    0.000)    8.611   2.244   (  -8.331    1.751    0.000)    8.513   2.563   (  10.624   -0.757   -0.000)   10.650   2.599   (  -0.807    5.523    0.000)    5.582   3.222   (  -0.878   12.164    0.000)   12.196   3.292   (   5.623    8.963    0.000)   10.581   4.001   (  -0.282   -4.083   -0.000)    4.093   4.063   (  -4.233   -7.120   -0.000)    8.283   4.135   (  11.858   -4.521   -0.000)   12.690   5.369   (  -0.187   13.198    0.000)   13.199   7.499   (  16.617    2.487   -0.000)   16.802   8.288   (   0.586   -1.502   -0.000)    1.612   8.968   (  -5.845    7.277    0.000)    9.334======================= Grid point 65 (30/60) =======================q-point: ( 0.11  0.22  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.285   (   3.486    3.633    0.912)    5.117   1.310   (   4.376    2.399    1.942)    5.355   1.407   (   0.347    1.212   -1.671)    2.094   1.422   (  -0.917    1.166   -0.137)    1.490   1.503   (  -0.699   -0.221    1.294)    1.487   1.588   (  -1.618   -0.175   -2.231)    2.761   1.717   (   0.549   -3.002   -1.077)    3.236   1.764   (  -3.917   -5.376   -0.112)    6.653   1.847   (  -2.049   -8.191    0.046)    8.444   2.247   (  -1.494    6.364    0.819)    6.588   2.449   (   1.559   -5.274   -1.225)    5.635   2.987   (  -0.068    8.471    1.284)    8.568   3.327   (   7.604   -4.453    0.435)    8.822   3.536   (  -8.157   -6.470   -2.522)   10.713   3.737   (   0.459    4.722    1.605)    5.008   3.935   (  -0.372   -1.687   -0.851)    1.926   4.250   (  -1.403   -3.215   -0.023)    3.508   5.936   (   5.121   14.858   -0.270)   15.718   8.388   (   1.098    2.103    0.768)    2.494   8.449   (   2.546   18.143   -0.055)   18.321   8.827   (   2.841   -4.479    0.110)    5.305======================= Grid point 67 (31/60) =======================q-point: ( 0.11 -0.11  0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.199   (   2.702    9.111    0.000)    9.503   1.236   (   2.625    3.929    0.000)    4.725   1.391   (   2.215    0.722   -0.000)    2.330   1.396   (   1.616   -2.362   -0.000)    2.861   1.501   (  -1.087   -0.623    0.000)    1.253   1.635   (  -0.802   -0.586    0.000)    0.993   1.767   (  -2.011   -1.036    0.000)    2.262   1.908   (  -7.714   -3.623    0.000)    8.522   2.036   (  -3.782   -0.051    0.000)    3.782   2.079   (   2.983   -1.464   -0.000)    3.323   2.588   (  -5.599    0.976    0.000)    5.684   2.774   (  10.191    0.081   -0.000)   10.192   3.409   (  -3.697    9.990    0.000)   10.652   3.551   (   8.432   -0.976   -0.000)    8.488   3.752   (  -7.976   -9.567   -0.000)   12.456   3.962   (  -0.515    1.295    0.000)    1.394   4.258   (   4.394   -1.581   -0.000)    4.670   5.642   (  10.482    4.058   -0.000)   11.240   7.998   (  23.366    1.798   -0.000)   23.435   8.333   (   2.208    2.191   -0.000)    3.110   8.944   (  -6.420    2.752    0.000)    6.985======================= Grid point 68 (32/60) =======================q-point: ( 0.11 -0.78  0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.348   (  -0.900    0.900   -0.495)    1.365   1.359   (   0.362   -0.362   -0.246)    0.568   1.427   (  -1.878    1.878    0.413)    2.688   1.445   (  -1.368    1.368   -0.186)    1.944   1.522   (   0.250   -0.250   -0.198)    0.405   1.577   (   0.013   -0.013   -0.943)    0.943   1.620   (   0.581   -0.581    2.536)    2.666   1.694   (   2.075   -2.075   -1.186)    3.165   1.720   (   3.480   -3.480    0.038)    4.922   2.338   (   0.498   -0.498    1.942)    2.065   2.363   (  -0.565    0.565   -1.733)    1.909   3.161   (  -3.348    3.348   -1.511)    4.970   3.263   (  -0.881    0.881    0.553)    1.363   3.373   (   5.985   -5.985    0.472)    8.478   3.773   (  -2.442    2.442    5.461)    6.462   3.935   (   1.110   -1.110   -4.383)    4.656   4.189   (   1.701   -1.701    0.022)    2.406   6.166   (  -6.178    6.178    0.071)    8.737   8.428   (  -0.390    0.390   -0.345)    0.651   8.707   (  -5.448    5.448    0.360)    7.714   8.776   (   0.131   -0.131   -0.331)    0.379======================= Grid point 70 (33/60) =======================q-point: (-0.89 -0.11  0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.332   (   1.489   -0.318   -0.000)    1.523   1.341   (   2.739    1.452   -0.000)    3.100   1.371   (   2.879   -3.377   -0.000)    4.438   1.433   (  -0.882   -3.279   -0.000)    3.396   1.497   (   1.585    4.808    0.000)    5.062   1.599   (  -0.800   -0.494    0.000)    0.940   1.675   (  -4.998   -2.580    0.000)    5.625   1.743   (  -2.853    1.016    0.000)    3.029   1.847   (  -8.349    2.055    0.000)    8.598   2.247   (   6.508    2.264   -0.000)    6.890   2.441   (  -5.225   -3.284    0.000)    6.171   3.015   (   9.553    3.509   -0.000)   10.177   3.366   (  -7.863   -9.449   -0.000)   12.293   3.493   (  -5.865    7.430    0.000)    9.466   3.622   (   5.987   -5.080   -0.000)    7.852   4.001   (  -0.935    2.739    0.000)    2.894   4.249   (  -3.316    1.319    0.000)    3.569   5.928   (  15.157   -5.927   -0.000)   16.275   8.422   (   0.863    1.953    0.000)    2.136   8.446   (  18.026   -2.704   -0.000)   18.227   8.822   (  -3.723   -3.308    0.000)    4.980======================= Grid point 77 (34/60) =======================q-point: ( 0.11  0.56 -0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.064   (   2.882   -6.806    4.561)    8.685   1.103   (   4.151   -7.525    3.050)    9.119   1.276   (   1.248   -1.913   -1.178)    2.570   1.393   (  -0.061   -5.028   -0.699)    5.077   1.561   (   0.536   -3.407    3.791)    5.125   1.592   (  -1.012    3.604    0.417)    3.767   1.854   (  -1.601    2.385   -1.181)    3.106   1.917   (   8.809    3.568   -0.266)    9.508   2.060   (  -9.608    0.699   -3.020)   10.095   2.261   (  -7.817   -1.435    0.849)    7.993   2.506   (   0.674   -4.612    0.845)    4.737   2.625   (   8.657   -0.594   -1.676)    8.838   3.074   (   8.994  -12.355   -3.999)   15.796   3.384   (  -2.733  -12.949    0.802)   13.258   3.906   (  -3.105    8.594    0.875)    9.179   4.117   (  12.025    4.105   -0.611)   12.721   4.176   (  -1.518    6.791    0.601)    6.984   5.315   (  -0.907  -16.240   -1.656)   16.350   7.499   (  16.562   -2.449    0.000)   16.742   8.466   (  -0.919    4.600    8.444)    9.659   8.872   (  -3.884   -9.251   -5.784)   11.581======================= Grid point 80 (35/60) =======================q-point: ( 0.11  0.56 -0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.044   (   4.018   -4.018    8.292)   10.052   1.059   (   6.585   -6.585    3.646)   10.001   1.276   (   2.253   -2.253   -3.052)    4.412   1.302   (   4.558   -4.558   -1.132)    6.544   1.516   (   2.448   -2.448    6.333)    7.217   1.633   (  -1.303    1.303   -1.510)    2.383   1.852   (  -4.888    4.888   -1.531)    7.080   1.954   (  -2.952    2.952   -6.095)    7.388   2.138   (  -3.040    3.040    1.150)    4.451   2.221   (   2.223   -2.223    0.829)    3.252   2.251   (   3.193   -3.193   -3.548)    5.743   2.875   (   6.145   -6.145    4.620)    9.842   3.024   (   5.054   -5.054   -0.243)    7.152   3.115   (  11.734  -11.734   -6.359)   17.771   3.966   (  -6.443    6.443    0.202)    9.114   4.244   (  -5.521    5.521    2.006)    8.062   4.347   (   1.054   -1.054   -0.083)    1.493   5.090   (  14.563  -14.563   -3.627)   20.912   7.727   (  12.242  -12.242    0.204)   17.314   8.416   (   5.716   -5.716    5.228)    9.627   8.744   (  -5.285    5.285   -2.373)    7.841======================= Grid point 84 (36/60) =======================q-point: ( 0.17 -0.17 -0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.865   ( -15.269   -0.000    0.000)   15.269   0.977   ( -11.635   -0.000    0.000)   11.635   1.279   (  -1.214   -0.000    0.000)    1.214   1.404   (  -0.625   -0.000    0.000)    0.625   1.420   ( -11.024   -0.000    0.000)   11.024   1.677   (   3.039    0.000   -0.000)    3.039   1.773   (   3.889    0.000   -0.000)    3.889   1.881   (   5.308    0.000   -0.000)    5.308   2.316   (   1.030    0.000   -0.000)    1.030   2.364   (  -2.654   -0.000    0.000)    2.654   2.365   (   1.872    0.000   -0.000)    1.872   2.571   (  -6.693   -0.000    0.000)    6.693   3.231   ( -12.855   -0.000    0.000)   12.855   3.254   (  -6.764   -0.000    0.000)    6.764   3.950   (   3.954    0.000   -0.000)    3.954   4.021   (   4.256    0.000   -0.000)    4.256   4.108   (   6.643    0.000   -0.000)    6.643   5.391   ( -15.018   -0.000    0.000)   15.018   7.308   (  -1.028   -0.000    0.000)    1.028   8.283   (   2.244    0.000   -0.000)    2.244   9.038   (  -6.743   -0.000    0.000)    6.743======================= Grid point 85 (37/60) =======================q-point: ( 0.17  0.17 -0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 45Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.205   (  -0.000   -0.000    9.979)    9.979   0.205   (  -0.000   -0.000    9.979)    9.979   0.602   (  -0.000   -0.000   30.644)   30.644   1.202   (   0.000    0.000   -1.673)    1.673   1.202   (   0.000    0.000   -1.673)    1.673   1.800   (   0.000    0.000   -3.851)    3.851   1.931   (  -0.000   -0.000    0.914)    0.914   1.931   (  -0.000   -0.000    0.914)    0.914   2.275   (   0.000    0.000   -0.032)    0.032   2.275   (   0.000    0.000   -0.032)    0.032   2.349   (   0.000    0.000   -1.154)    1.154   2.349   (   0.000    0.000   -1.154)    1.154   2.535   (   0.000    0.000   -0.340)    0.340   2.776   (   0.000    0.000  -11.697)   11.697   4.119   (   0.000    0.000   -0.068)    0.068   4.119   (   0.000    0.000   -0.068)    0.068   4.294   (   0.000    0.000   -2.298)    2.298   4.622   (  -0.000   -0.000    8.906)    8.906   7.256   (  -0.000   -0.000    0.005)    0.005   7.256   (  -0.000   -0.000    0.005)    0.005   9.751   (   0.000    0.000   -2.540)    2.540======================= Grid point 86 (38/60) =======================q-point: ( 0.17  0.50 -0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.022   (  10.061    0.000    1.501)   10.172   1.057   (   6.808    0.000    5.726)    8.896   1.263   (   1.890    0.000   -0.602)    1.983   1.394   (   5.695    0.000   -0.760)    5.746   1.553   (   3.253    0.000    3.080)    4.480   1.608   (  -2.412   -0.000    1.680)    2.939   1.797   (  -4.389   -0.000    0.782)    4.459   1.883   (  -2.320   -0.000   -2.087)    3.121   2.274   (   0.231   -0.000    0.895)    0.925   2.295   (  -0.435   -0.000   -0.677)    0.805   2.402   (   1.504    0.000   -1.731)    2.293   2.583   (   6.576    0.000    0.367)    6.586   2.966   (  13.168    0.000   -6.330)   14.610   3.418   (  12.884    0.000    1.330)   12.953   3.939   (  -8.804   -0.000    1.689)    8.965   3.952   (  -3.966   -0.000    0.039)    3.966   4.198   (  -6.807   -0.000   -0.438)    6.821   5.346   (  17.611    0.000   -1.435)   17.669   7.309   (   1.016    0.000    0.031)    1.017   8.477   (  -5.236   -0.000    9.434)   10.789   8.918   (   9.061    0.000   -6.941)   11.414======================= Grid point 94 (39/60) =======================q-point: ( 0.17  0.17  0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.163   (   6.740    6.740    3.873)   10.288   1.225   (   3.938    3.938    1.365)    5.734   1.406   (   0.511    0.511   -0.873)    1.133   1.410   (   1.356    1.356   -1.469)    2.416   1.508   (   0.296    0.296    1.566)    1.621   1.622   (  -0.953   -0.953   -1.719)    2.184   1.735   (  -2.532   -2.532   -2.479)    4.355   1.858   (  -4.867   -4.867   -1.264)    6.998   1.927   (  -5.385   -5.385    0.323)    7.623   2.277   (   1.534    1.534    1.912)    2.893   2.405   (   0.088    0.088   -3.383)    3.385   2.968   (   3.696    3.696    3.605)    6.349   3.202   (   2.830    2.830   -0.703)    4.063   3.628   (   4.851    4.851   -2.193)    7.202   3.729   (  -5.337   -5.337   -3.781)    8.442   3.982   (  -3.518   -3.518    2.918)    5.768   4.299   (  -2.732   -2.732   -0.075)    3.865   5.733   (  13.267   13.267   -0.983)   18.788   8.293   (  13.298   13.298    0.038)   18.806   8.364   (   1.147    1.147    2.429)    2.920   8.773   (   0.325    0.325   -0.306)    0.552======================= Grid point 96 (40/60) =======================q-point: ( 0.17 -0.17  0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.169   (   0.283   11.678    0.000)   11.681   1.187   (   1.847    4.821    0.000)    5.163   1.333   (   2.771    1.256   -0.000)    3.042   1.371   (   0.873   -4.086   -0.000)    4.179   1.526   (  -1.281   -0.175    0.000)    1.293   1.650   (  -0.594    0.865    0.000)    1.049   1.803   (  -1.246   -0.973    0.000)    1.581   1.839   (  10.340   -1.184   -0.000)   10.407   2.091   (  -9.662   -2.381    0.000)    9.951   2.278   (  -9.975    1.261    0.000)   10.054   2.548   (  11.166   -0.464   -0.000)   11.176   2.674   (  -2.065    1.401    0.000)    2.495   3.371   (   7.483   -1.287   -0.000)    7.592   3.469   (  -1.909   11.051    0.000)   11.215   3.894   (  -5.048   -8.789   -0.000)   10.136   3.969   (   0.441    1.317    0.000)    1.389   4.056   (  11.626   -2.509   -0.000)   11.893   5.569   (  -2.347    4.885    0.000)    5.420   7.537   (  18.370    0.978   -0.000)   18.396   8.294   (   1.335    2.240    0.000)    2.608   9.071   (  -5.016    1.868    0.000)    5.353======================= Grid point 97 (41/60) =======================q-point: ( 0.17  0.17  0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.338   (   1.454    1.454    0.427)    2.100   1.351   (   1.006    1.006    0.517)    1.514   1.440   (   1.950    1.950   -1.016)    2.939   1.442   (   0.386    0.386    0.089)    0.553   1.504   (   0.366    0.366    1.576)    1.659   1.582   (  -0.583   -0.583   -1.324)    1.560   1.668   (  -3.788   -3.788    1.359)    5.526   1.674   (  -0.547   -0.547   -2.640)    2.751   1.720   (  -3.471   -3.471    0.018)    4.908   2.325   (   0.668    0.668    1.922)    2.142   2.378   (  -0.768   -0.768   -1.858)    2.153   3.134   (   4.598    4.598    1.128)    6.599   3.274   (   0.747    0.747   -0.529)    1.181   3.389   (  -7.061   -7.061   -1.114)   10.048   3.795   (   0.539    0.539    4.675)    4.737   3.919   (   0.279    0.279   -4.096)    4.115   4.189   (  -1.741   -1.741   -0.023)    2.462   6.168   (   6.016    6.016   -0.087)    8.508   8.422   (   1.026    1.026    0.288)    1.480   8.714   (   5.599    5.599   -0.360)    7.926   8.768   (  -0.295   -0.295    0.401)    0.579======================= Grid point 99 (42/60) =======================q-point: (-0.83 -0.17  0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.274   (   3.215   -0.417   -0.000)    3.242   1.293   (   2.119   -2.539   -0.000)    3.307   1.316   (   1.938   -1.994   -0.000)    2.780   1.458   (  -1.235   -2.950   -0.000)    3.198   1.467   (   1.388    5.939    0.000)    6.100   1.627   (  -1.755   -0.360    0.000)    1.792   1.760   (  -2.843   -0.153    0.000)    2.847   1.852   (  -7.676   -1.436    0.000)    7.809   2.051   (  -5.996    1.358    0.000)    6.147   2.076   (   4.795    1.222   -0.000)    4.948   2.558   (  -5.360   -3.790    0.000)    6.565   2.798   (  11.075    2.448   -0.000)   11.342   3.460   (   9.140   -7.037   -0.000)   11.535   3.529   (  -6.326  -10.126   -0.000)   11.940   3.621   (  -6.702    8.570    0.000)   10.879   4.021   (  -0.980    4.030    0.000)    4.148   4.260   (   4.435    1.671   -0.000)    4.740   5.623   (  10.678   -5.952   -0.000)   12.225   7.996   (  23.363   -1.864   -0.000)   23.437   8.402   (   0.872    4.460    0.000)    4.544   8.921   (  -5.291   -5.213    0.000)    7.428======================= Grid point 103 (43/60) =======================q-point: ( 0.17  0.83 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.286   (   4.206   -4.206   -0.000)    5.948   1.309   (   1.519   -1.519   -0.000)    2.149   1.322   (   3.853   -3.853   -0.000)    5.449   1.385   (   0.300   -0.300   -0.000)    0.425   1.576   (  -1.403    1.403    0.000)    1.985   1.599   (  -0.509    0.509    0.000)    0.719   1.667   (  -2.929    2.929    0.000)    4.142   1.777   (  -2.168    2.168    0.000)    3.065   1.931   (  -5.845    5.845    0.000)    8.266   2.316   (   1.054   -1.054   -0.000)    1.491   2.334   (   1.788   -1.788   -0.000)    2.529   3.091   (   2.015   -2.015   -0.000)    2.850   3.179   (   2.735   -2.735   -0.000)    3.868   3.479   (   8.506   -8.506   -0.000)   12.029   3.657   (  -7.682    7.682    0.000)   10.864   4.068   (  -3.365    3.365    0.000)    4.759   4.297   (  -2.736    2.736    0.000)    3.869   5.704   (  14.466  -14.466   -0.000)   20.458   8.294   (  13.030  -13.030   -0.000)   18.427   8.466   (  -2.380    2.380    0.000)    3.367   8.739   (   2.661   -2.661   -0.000)    3.763======================= Grid point 110 (44/60) =======================q-point: ( 0.17  0.50 -0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.899   (   8.400   -9.177    3.330)   12.879   0.977   (   2.588   -3.864    9.921)   10.957   1.212   (   2.369   -8.223   -0.697)    8.585   1.242   (   1.774   -4.646   -2.389)    5.517   1.467   (   1.939   -4.642    7.956)    9.413   1.655   (  -0.908    1.979   -1.979)    2.943   1.933   (  -1.786    5.298   -2.538)    6.140   1.989   (   7.873    2.893   -0.737)    8.420   2.084   (  -7.541    1.485   -3.895)    8.616   2.228   (  -3.812   -1.902    0.312)    4.272   2.345   (  -5.241   -9.076   -1.318)   10.563   2.638   (  13.214   -0.458    1.288)   13.284   2.805   (  11.759  -11.142   -4.984)   16.948   3.082   (  -2.457  -15.560    0.019)   15.753   4.073   (  -3.547    6.973    0.705)    7.855   4.201   (  11.701    3.363   -0.517)   12.186   4.333   (  -2.806    7.800    2.433)    8.639   4.914   (   1.833  -20.496   -3.988)   20.960   7.439   (  12.955   -3.181    0.063)   13.340   8.312   (   3.168  -19.332    1.797)   19.672   8.905   (  -7.725   12.880    0.336)   15.023======================= Grid point 114 (45/60) =======================q-point: ( 0.22 -0.22 -0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.166   (  -5.602   -0.000    0.000)    5.602   1.171   ( -13.483   -0.000    0.000)   13.483   1.297   (  -0.349   -0.000    0.000)    0.349   1.362   (   4.477    0.000   -0.000)    4.477   1.546   (   0.200    0.000   -0.000)    0.200   1.659   (  -2.432   -0.000    0.000)    2.432   1.704   (   1.968    0.000   -0.000)    1.968   1.816   (   0.938    0.000   -0.000)    0.938   2.297   (   0.711    0.000   -0.000)    0.711   2.326   (   1.745    0.000   -0.000)    1.745   2.407   (  -1.230   -0.000    0.000)    1.230   2.682   (  -3.230   -0.000    0.000)    3.230   3.286   (   3.150    0.000   -0.000)    3.150   3.489   ( -11.304   -0.000    0.000)   11.304   3.885   (   1.829    0.000   -0.000)    1.829   3.949   (   8.338    0.000   -0.000)    8.338   3.985   (  -0.814   -0.000    0.000)    0.814   5.620   (  -5.717   -0.000    0.000)    5.717   7.324   (  -0.431   -0.000    0.000)    0.431   8.280   (  -2.244   -0.000    0.000)    2.244   9.128   (  -1.285   -0.000    0.000)    1.285======================= Grid point 116 (46/60) =======================q-point: ( 0.22  0.44 -0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.779   (  12.672    0.000    2.709)   12.959   0.948   (   4.125    0.000   10.457)   11.242   1.194   (  12.424    0.000   -0.813)   12.451   1.219   (   2.173    0.000   -1.178)    2.472   1.445   (   5.631    0.000    8.413)   10.124   1.667   (  -2.041   -0.000   -2.276)    3.057   1.881   (  -3.358   -0.000    1.026)    3.511   1.954   (  -4.593   -0.000   -4.119)    6.170   2.239   (   3.719    0.000    1.026)    3.858   2.305   (  -0.497   -0.000   -1.046)    1.158   2.389   (  -0.315    0.000   -2.057)    2.081   2.437   (   6.644    0.000    0.307)    6.651   2.657   (  14.971    0.000   -7.664)   16.819   3.114   (  15.775    0.000    0.014)   15.775   4.035   (  -3.633   -0.000    0.020)    3.633   4.111   (  -7.183   -0.000    0.909)    7.241   4.363   (  -8.542   -0.000    2.418)    8.878   4.916   (  21.854    0.000   -3.283)   22.100   7.285   (   1.152    0.000    0.029)    1.152   8.281   (  23.990    0.000    1.101)   24.015   8.998   ( -16.689   -0.000    0.346)   16.693======================= Grid point 129 (47/60) =======================q-point: (-0.78 -0.22  0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.214   (   2.381   -2.336   -0.000)    3.336   1.252   (   1.521   -5.768   -0.000)    5.965   1.293   (   0.413   -0.347   -0.000)    0.539   1.437   (   1.263    8.089    0.000)    8.187   1.481   (  -0.887   -3.338   -0.000)    3.454   1.666   (  -1.746    0.081    0.000)    1.748   1.805   (  -1.502    0.747    0.000)    1.678   1.851   (   9.898    2.038   -0.000)   10.106   2.048   ( -10.715   -1.628    0.000)   10.838   2.290   (  -9.935   -0.235    0.000)    9.938   2.552   (  12.207    1.021   -0.000)   12.249   2.651   (  -3.294   -3.547   -0.000)    4.840   3.261   (   8.944   -8.580   -0.000)   12.394   3.616   (  -2.389   -4.188   -0.000)    4.822   3.759   (  -5.396    3.077    0.000)    6.211   4.038   (   1.407    3.223    0.000)    3.517   4.050   (  10.728    3.505   -0.000)   11.286   5.552   (  -2.575   -6.633   -0.000)    7.115   7.537   (  18.300   -1.022   -0.000)   18.329   8.388   (   0.407    6.824    0.000)    6.836   9.025   (  -4.009   -6.707   -0.000)    7.814======================= Grid point 133 (48/60) =======================q-point: (-0.78 -0.22  0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.200   (   3.712   -4.225   -0.000)    5.624   1.222   (   5.273   -4.268   -0.000)    6.784   1.281   (   1.108   -2.215   -0.000)    2.476   1.385   (   0.039   -4.051   -0.000)    4.051   1.578   (   0.406    3.614    0.000)    3.637   1.619   (  -1.355   -0.260    0.000)    1.380   1.765   (  -4.658    1.344    0.000)    4.848   1.846   (  -5.292    0.514    0.000)    5.317   2.086   (  -8.531    2.202    0.000)    8.811   2.132   (  10.341    3.775   -0.000)   11.008   2.431   (  -5.049   -8.234   -0.000)    9.659   2.877   (  12.705    4.762   -0.000)   13.568   3.269   (  11.235  -10.942   -0.000)   15.683   3.282   (  -6.240  -12.728   -0.000)   14.176   3.809   (  -6.604    8.547    0.000)   10.801   4.129   (  -2.308    6.123    0.000)    6.543   4.313   (   3.288    2.818   -0.000)    4.331   5.398   (  11.661  -14.865   -0.000)   18.893   7.917   (  21.031   -5.925   -0.000)   21.850   8.524   (  -3.082    7.235    0.000)    7.864   8.734   (  -1.616  -11.753   -0.000)   11.864======================= Grid point 144 (49/60) =======================q-point: ( 0.22  0.44 -0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.699   (  10.297  -10.297    4.282)   15.178   0.843   (  10.842  -10.842    6.582)   16.686   1.048   (   1.971   -1.971    6.556)    7.124   1.195   (   1.455   -1.455   -2.276)    3.069   1.369   (   4.489   -4.489    9.097)   11.092   1.686   (  -1.150    1.150   -3.196)    3.586   2.031   (   1.958   -1.958    1.471)    3.135   2.032   (   1.148   -1.148   -3.243)    3.626   2.107   (  -4.627    4.627   -1.781)    6.782   2.159   (  -1.562    1.562   -6.671)    7.027   2.262   (  -1.141    1.141    0.409)    1.665   2.490   (   9.631   -9.631   -0.187)   13.622   2.723   (   7.098   -7.098   -3.984)   10.799   2.766   (   6.969   -6.969    0.033)    9.855   4.187   (  -3.764    3.764   -0.036)    5.323   4.219   (   4.005   -4.005   -0.135)    5.665   4.413   (   4.925   -4.925   -0.380)    6.976   4.611   (   0.568   -0.568   -0.482)    0.936   7.356   (   5.208   -5.208    0.104)    7.365   7.863   (  18.981  -18.981   -3.569)   27.079   9.238   ( -15.746   15.746    3.445)   22.533======================= Grid point 145 (50/60) =======================q-point: ( 0.28  0.72 -0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.185   (   3.735    0.000   -0.000)    3.735   1.235   (   6.363    0.000   -0.000)    6.363   1.291   (   0.990    0.000   -0.000)    0.990   1.422   (  -9.548   -0.000    0.000)    9.548   1.493   (   3.755    0.000   -0.000)    3.755   1.697   (  -0.053   -0.000    0.000)    0.053   1.717   (  -3.108   -0.000    0.000)    3.108   1.821   (  -0.931   -0.000    0.000)    0.931   2.287   (   0.198    0.000   -0.000)    0.198   2.298   (   0.761    0.000   -0.000)    0.761   2.411   (   0.891    0.000   -0.000)    0.891   2.689   (   2.606    0.000   -0.000)    2.606   3.148   (   9.475    0.000   -0.000)    9.475   3.640   (   2.508    0.000   -0.000)    2.508   3.819   (  -1.543   -0.000    0.000)    1.543   3.886   (  -1.891   -0.000    0.000)    1.891   4.047   (  -4.867   -0.000    0.000)    4.867   5.606   (   7.169    0.000   -0.000)    7.169   7.324   (   0.380    0.000   -0.000)    0.380   8.384   (  -7.839   -0.000    0.000)    7.839   9.070   (   7.306    0.000   -0.000)    7.306======================= Grid point 148 (51/60) =======================q-point: ( 0.28  0.39 -0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.513   (  11.660    0.000    4.484)   12.493   0.715   (  20.341    0.000    5.295)   21.019   1.051   (  -1.488   -0.000    9.893)   10.004   1.179   (   1.441    0.000   -1.731)    2.253   1.311   (   7.249    0.000    9.301)   11.792   1.701   (  -1.292    0.000   -3.681)    3.901   1.934   (  -1.670   -0.000    1.062)    1.979   2.044   (   2.373    0.000   -1.563)    2.841   2.152   (  -2.810   -0.000   -3.857)    4.772   2.308   (   1.317    0.000   -0.518)    1.415   2.329   (   2.320    0.000   -0.629)    2.403   2.337   (   5.448    0.000   -0.106)    5.449   2.486   (  -0.378    0.000   -7.626)    7.636   2.774   (  16.025    0.000   -2.022)   16.152   4.096   (  -2.026   -0.000   -0.037)    2.026   4.227   (  -3.727   -0.000   -0.346)    3.743   4.361   (  16.358    0.000   -2.338)   16.524   4.662   (  -6.152   -0.000    4.241)    7.472   7.264   (   0.747    0.000    0.014)    0.747   7.659   (  31.005    0.000   -3.443)   31.195   9.432   ( -21.081   -0.000    1.953)   21.171======================= Grid point 165 (52/60) =======================q-point: (-0.72 -0.28  0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.107   (   5.316   -7.155   -0.000)    8.914   1.138   (   2.071   -4.329   -0.000)    4.799   1.263   (   0.604   -2.211   -0.000)    2.292   1.384   (   0.045   -6.362   -0.000)    6.362   1.568   (   0.455    3.591    0.000)    3.620   1.650   (  -1.391   -1.280    0.000)    1.890   1.828   (  -1.731    1.576    0.000)    2.341   1.915   (   9.155    3.766   -0.000)    9.900   2.026   ( -11.231   -0.435    0.000)   11.240   2.271   (  -7.995   -1.421    0.000)    8.120   2.529   (  -4.211   -7.889   -0.000)    8.942   2.593   (  14.106    2.746   -0.000)   14.370   3.031   (  10.457  -12.877   -0.000)   16.588   3.389   (  -3.615  -13.433   -0.000)   13.911   3.923   (  -4.000    8.476    0.000)    9.373   4.108   (  10.736    4.216   -0.000)   11.534   4.179   (   0.537    7.450    0.000)    7.469   5.298   (  -1.061  -17.045   -0.000)   17.077   7.499   (  16.559   -2.426   -0.000)   16.735   8.587   (  -2.576   12.232    0.000)   12.500   8.780   (  -2.127  -16.642   -0.000)   16.778======================= Grid point 169 (53/60) =======================q-point: ( 0.28  0.72 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.106   (   6.726   -6.726   -0.000)    9.513   1.128   (   2.682   -2.682   -0.000)    3.793   1.236   (   1.798   -1.798   -0.000)    2.543   1.288   (   4.987   -4.987   -0.000)    7.052   1.607   (   0.364   -0.364   -0.000)    0.515   1.619   (  -0.283    0.283    0.000)    0.400   1.837   (  -5.017    5.017    0.000)    7.095   1.869   (  -2.336    2.336    0.000)    3.304   2.149   (  -3.521    3.521    0.000)    4.980   2.216   (   3.419   -3.419   -0.000)    4.835   2.230   (   2.779   -2.779   -0.000)    3.930   2.956   (   4.850   -4.850   -0.000)    6.860   3.012   (  13.121  -13.121   -0.000)   18.556   3.021   (   4.888   -4.888   -0.000)    6.913   3.971   (  -6.712    6.712    0.000)    9.492   4.262   (  -5.800    5.800    0.000)    8.203   4.346   (   1.071   -1.071   -0.000)    1.514   5.053   (  15.456  -15.456   -0.000)   21.859   7.729   (  12.245  -12.245   -0.000)   17.317   8.485   (  10.281  -10.281   -0.000)   14.540   8.704   (  -9.649    9.649    0.000)   13.645======================= Grid point 177 (54/60) =======================q-point: ( 0.33  0.67 -0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.038   (   9.504    0.000   -0.000)    9.504   1.116   (   4.392    0.000   -0.000)    4.392   1.257   (   2.161    0.000   -0.000)    2.161   1.383   (   7.298    0.000   -0.000)    7.298   1.562   (  -3.460   -0.000    0.000)    3.460   1.669   (   2.050    0.000   -0.000)    2.050   1.805   (  -4.581   -0.000    0.000)    4.581   1.846   (  -1.465   -0.000    0.000)    1.465   2.288   (  -0.217   -0.000    0.000)    0.217   2.301   (  -1.225   -0.000    0.000)    1.225   2.368   (   3.256    0.000   -0.000)    3.256   2.586   (   6.543    0.000   -0.000)    6.543   2.898   (  14.130    0.000   -0.000)   14.130   3.428   (  13.521    0.000   -0.000)   13.521   3.952   (  -3.970   -0.000    0.000)    3.970   3.967   (  -8.248   -0.000    0.000)    8.248   4.185   (  -8.229   -0.000    0.000)    8.229   5.331   (  18.297    0.000   -0.000)   18.297   7.309   (   1.012    0.000   -0.000)    1.012   8.617   ( -14.481   -0.000    0.000)   14.481   8.805   (  18.220    0.000   -0.000)   18.220======================= Grid point 180 (55/60) =======================q-point: ( 0.33  0.33 -0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 60Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.365   (  -0.000   -0.000    6.351)    6.351   0.365   (  -0.000   -0.000    6.351)    6.351   1.143   (  -0.000   -0.000   24.779)   24.779   1.163   (   0.000    0.000   -1.967)    1.967   1.163   (   0.000    0.000   -1.967)    1.967   1.717   (   0.000    0.000   -4.095)    4.095   1.951   (  -0.000   -0.000    1.031)    1.031   1.951   (  -0.000   -0.000    1.031)    1.031   2.276   (  -0.000   -0.000    0.201)    0.201   2.276   (  -0.000   -0.000    0.201)    0.201   2.322   (   0.000    0.000   -1.502)    1.502   2.322   (   0.000    0.000   -1.502)    1.502   2.495   (   0.000    0.000  -15.924)   15.924   2.528   (   0.000    0.000   -0.341)    0.341   4.117   (   0.000    0.000   -0.067)    0.067   4.117   (   0.000    0.000   -0.067)    0.067   4.266   (   0.000    0.000   -0.738)    0.738   4.782   (  -0.000   -0.000    6.755)    6.755   7.256   (  -0.000   -0.000    0.005)    0.005   7.256   (  -0.000   -0.000    0.005)    0.005   9.697   (   0.000    0.000   -2.594)    2.594======================= Grid point 202 (56/60) =======================q-point: (-0.67 -0.33  0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.935   (   9.130   -8.748   -0.000)   12.645   1.084   (   1.481   -0.757   -0.000)    1.663   1.194   (   3.368  -11.891   -0.000)   12.359   1.211   (   0.607   -2.483   -0.000)    2.556   1.585   (   1.261   -1.546   -0.000)    1.995   1.630   (  -0.664   -0.315    0.000)    0.735   1.885   (  -1.044    4.315    0.000)    4.440   1.995   (   9.037    3.607   -0.000)    9.731   2.031   ( -11.373    0.904    0.000)   11.409   2.232   (  -3.356   -3.036    0.000)    4.525   2.334   (  -5.987   -9.380   -0.000)   11.128   2.658   (  16.073    3.099   -0.000)   16.368   2.743   (  11.737  -13.715   -0.000)   18.052   3.080   (  -2.950  -15.851   -0.000)   16.123   4.083   (  -3.719    6.464    0.000)    7.457   4.195   (  11.707    3.530   -0.000)   12.228   4.357   (  -2.840    9.029    0.000)    9.465   4.875   (   1.327  -21.671   -0.000)   21.711   7.440   (  13.021   -3.143   -0.000)   13.396   8.338   (   3.946  -24.519   -0.000)   24.834   8.901   (  -8.012   17.243    0.000)   19.013======================= Grid point 210 (57/60) =======================q-point: ( 0.39  0.61 -0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.808   (  11.984    0.000   -0.000)   11.984   1.068   (   0.032    0.000   -0.000)    0.032   1.159   (  14.498    0.000   -0.000)   14.498   1.206   (   2.491    0.000   -0.000)    2.491   1.570   (   2.363    0.000   -0.000)    2.363   1.639   (   0.422    0.000   -0.000)    0.422   1.892   (  -3.429   -0.000    0.000)    3.429   1.896   (  -3.941   -0.000    0.000)    3.941   2.270   (   6.351    0.000   -0.000)    6.351   2.294   (  -0.352   -0.000    0.000)    0.352   2.345   (  -2.585   -0.000    0.000)    2.585   2.441   (   6.740    0.000   -0.000)    6.740   2.575   (  15.335    0.000   -0.000)   15.335   3.112   (  16.330    0.000   -0.000)   16.330   4.035   (  -3.610   -0.000    0.000)    3.610   4.123   (  -6.299   -0.000    0.000)    6.299   4.386   ( -10.166   -0.000    0.000)   10.166   4.883   (  22.461    0.000   -0.000)   22.461   7.285   (   1.157    0.000   -0.000)    1.157   8.300   (  29.456    0.000   -0.000)   29.456   8.994   ( -20.978   -0.000    0.000)   20.978======================= Grid point 240 (58/60) =======================q-point: ( 0.39  0.61 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.746   (  10.113  -10.113   -0.000)   14.302   0.903   (  13.337  -13.337   -0.000)   18.862   1.093   (  -0.937    0.937    0.000)    1.326   1.169   (   1.240   -1.240   -0.000)    1.753   1.528   (   3.086   -3.086   -0.000)    4.365   1.638   (  -0.827    0.827    0.000)    1.170   2.009   (  -4.834    4.834    0.000)    6.836   2.056   (   3.785   -3.785   -0.000)    5.352   2.061   (   3.633   -3.633   -0.000)    5.137   2.079   (  -6.399    6.399    0.000)    9.050   2.266   (  -1.026    1.026    0.000)    1.451   2.489   (   9.874   -9.874   -0.000)   13.964   2.674   (   8.332   -8.332   -0.000)   11.783   2.767   (   7.032   -7.032   -0.000)    9.944   4.188   (  -3.452    3.452    0.000)    4.882   4.217   (   3.995   -3.995   -0.000)    5.650   4.399   (   8.996   -8.996   -0.000)   12.722   4.620   (  -5.050    5.050    0.000)    7.142   7.357   (   5.245   -5.245   -0.000)    7.417   7.836   (  19.599  -19.599   -0.000)   27.717   9.264   ( -15.221   15.221    0.000)   21.525======================= Grid point 244 (59/60) =======================q-point: ( 0.44  0.56 -0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.560   (  10.620    0.000   -0.000)   10.620   0.770   (  21.798    0.000   -0.000)   21.798   1.103   (  -2.687   -0.000    0.000)    2.687   1.160   (   1.670    0.000   -0.000)    1.670   1.489   (   4.407    0.000   -0.000)    4.407   1.650   (  -1.218   -0.000    0.000)    1.218   1.945   (  -1.656   -0.000    0.000)    1.656   2.040   ( -10.090   -0.000    0.000)   10.090   2.084   (  10.006    0.000   -0.000)   10.006   2.301   (  -0.237   -0.000    0.000)    0.237   2.324   (   4.028    0.000   -0.000)    4.028   2.340   (   6.198    0.000   -0.000)    6.198   2.390   (  -1.061   -0.000    0.000)    1.061   2.754   (  16.923    0.000   -0.000)   16.923   4.095   (  -1.998   -0.000    0.000)    1.998   4.225   (  -3.335   -0.000    0.000)    3.335   4.338   (  17.271    0.000   -0.000)   17.271   4.709   (  -9.499   -0.000    0.000)    9.499   7.264   (   0.754    0.000   -0.000)    0.754   7.633   (  30.396    0.000   -0.000)   30.396   9.442   ( -19.183   -0.000    0.000)   19.183======================= Grid point 279 (60/60) =======================q-point: (-0.50 -0.50  0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 30Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.429   (  -0.000   -0.000    0.000)    0.000   0.429   (  -0.000   -0.000    0.000)    0.000   1.142   (  -0.000   -0.000    0.000)    0.000   1.142   (  -0.000   -0.000    0.000)    0.000   1.431   (  -0.000   -0.000    0.000)    0.000   1.669   (   0.000    0.000   -0.000)    0.000   1.962   (  -0.000   -0.000    0.000)    0.000   1.962   (  -0.000   -0.000    0.000)    0.000   2.280   (  -0.000   -0.000    0.000)    0.000   2.280   (  -0.000   -0.000    0.000)    0.000   2.295   (  -0.000   -0.000    0.000)    0.000   2.303   (  -0.000   -0.000    0.000)    0.000   2.303   (  -0.000   -0.000    0.000)    0.000   2.524   (  -0.000   -0.000    0.000)    0.000   4.116   (  -0.000   -0.000    0.000)    0.000   4.116   (  -0.000   -0.000    0.000)    0.000   4.260   (  -0.000   -0.000    0.000)    0.000   4.850   (  -0.000   -0.000    0.000)    0.000   7.256   (  -0.000   -0.000    0.000)    0.000   7.256   (  -0.000   -0.000    0.000)    0.000   9.670   (  -0.000   -0.000    0.000)    0.000=================== End of collection of collisions ===================----------- Thermal conductivity (W/m-k) with tetrahedron method -----------#  T(K)        xx         yy         zz         yz         xz         xy        #ipm    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/10206   10.0     86.617     86.617     54.051      0.000     -0.000      0.000 3/10206   20.0     13.944     13.944      8.943      0.000     -0.000      0.000 3/10206   30.0      6.312      6.312      3.954      0.000     -0.000      0.000 3/10206   40.0      4.109      4.109      2.465      0.000     -0.000      0.000 3/10206   50.0      3.127      3.127      1.791      0.000     -0.000      0.000 3/10206   60.0      2.570      2.570      1.412      0.000     -0.000      0.000 3/10206   70.0      2.204      2.204      1.169      0.000     -0.000      0.000 3/10206   80.0      1.941      1.941      0.999      0.000     -0.000      0.000 3/10206   90.0      1.739      1.739      0.874      0.000     -0.000      0.000 3/10206  100.0      1.578      1.578      0.777      0.000     -0.000      0.000 3/10206  110.0      1.446      1.446      0.700      0.000     -0.000      0.000 3/10206  120.0      1.335      1.335      0.638      0.000     -0.000      0.000 3/10206  130.0      1.240      1.240      0.586      0.000     -0.000      0.000 3/10206  140.0      1.158      1.158      0.542      0.000     -0.000      0.000 3/10206  150.0      1.087      1.087      0.504      0.000     -0.000      0.000 3/10206  160.0      1.023      1.023      0.471      0.000     -0.000      0.000 3/10206  170.0      0.967      0.967      0.442      0.000     -0.000      0.000 3/10206  180.0      0.917      0.917      0.417      0.000     -0.000      0.000 3/10206  190.0      0.871      0.871      0.395      0.000     -0.000      0.000 3/10206  200.0      0.830      0.830      0.374      0.000     -0.000      0.000 3/10206  210.0      0.792      0.792      0.356      0.000     -0.000      0.000 3/10206  220.0      0.758      0.758      0.340      0.000     -0.000      0.000 3/10206  230.0      0.726      0.726      0.325      0.000     -0.000      0.000 3/10206  240.0      0.697      0.697      0.311      0.000     -0.000      0.000 3/10206  250.0      0.671      0.671      0.298      0.000     -0.000      0.000 3/10206  260.0      0.646      0.646      0.287      0.000     -0.000      0.000 3/10206  270.0      0.623      0.623      0.276      0.000     -0.000      0.000 3/10206  280.0      0.601      0.601      0.266      0.000     -0.000      0.000 3/10206  290.0      0.581      0.581      0.257      0.000     -0.000      0.000 3/10206  300.0      0.562      0.562      0.248      0.000     -0.000      0.000 3/10206  310.0      0.545      0.545      0.240      0.000     -0.000      0.000 3/10206  320.0      0.528      0.528      0.233      0.000     -0.000      0.000 3/10206  330.0      0.513      0.513      0.225      0.000     -0.000      0.000 3/10206  340.0      0.498      0.498      0.219      0.000     -0.000      0.000 3/10206  350.0      0.484      0.484      0.212      0.000     -0.000      0.000 3/10206  360.0      0.471      0.471      0.207      0.000     -0.000      0.000 3/10206  370.0      0.458      0.458      0.201      0.000     -0.000      0.000 3/10206  380.0      0.447      0.447      0.196      0.000     -0.000      0.000 3/10206  390.0      0.435      0.435      0.191      0.000     -0.000      0.000 3/10206  400.0      0.425      0.425      0.186      0.000     -0.000      0.000 3/10206  410.0      0.414      0.414      0.181      0.000     -0.000      0.000 3/10206  420.0      0.405      0.405      0.177      0.000     -0.000      0.000 3/10206  430.0      0.396      0.396      0.173      0.000     -0.000      0.000 3/10206  440.0      0.387      0.387      0.169      0.000     -0.000      0.000 3/10206  450.0      0.378      0.378      0.165      0.000     -0.000      0.000 3/10206  460.0      0.370      0.370      0.161      0.000     -0.000      0.000 3/10206  470.0      0.362      0.362      0.158      0.000     -0.000      0.000 3/10206  480.0      0.355      0.355      0.155      0.000     -0.000      0.000 3/10206  490.0      0.348      0.348      0.152      0.000     -0.000      0.000 3/10206  500.0      0.341      0.341      0.149      0.000     -0.000      0.000 3/10206  510.0      0.334      0.334      0.146      0.000     -0.000      0.000 3/10206  520.0      0.328      0.328      0.143      0.000     -0.000      0.000 3/10206  530.0      0.322      0.322      0.140      0.000     -0.000      0.000 3/10206  540.0      0.316      0.316      0.137      0.000     -0.000      0.000 3/10206  550.0      0.310      0.310      0.135      0.000     -0.000      0.000 3/10206  560.0      0.305      0.305      0.133      0.000     -0.000      0.000 3/10206  570.0      0.299      0.299      0.130      0.000     -0.000      0.000 3/10206  580.0      0.294      0.294      0.128      0.000     -0.000      0.000 3/10206  590.0      0.289      0.289      0.126      0.000     -0.000      0.000 3/10206  600.0      0.285      0.285      0.124      0.000     -0.000      0.000 3/10206  610.0      0.280      0.280      0.122      0.000     -0.000      0.000 3/10206  620.0      0.276      0.276      0.120      0.000     -0.000      0.000 3/10206  630.0      0.271      0.271      0.118      0.000     -0.000      0.000 3/10206  640.0      0.267      0.267      0.116      0.000     -0.000      0.000 3/10206  650.0      0.263      0.263      0.114      0.000     -0.000      0.000 3/10206  660.0      0.259      0.259      0.112      0.000     -0.000      0.000 3/10206  670.0      0.255      0.255      0.111      0.000     -0.000      0.000 3/10206  680.0      0.251      0.251      0.109      0.000     -0.000      0.000 3/10206  690.0      0.248      0.248      0.108      0.000     -0.000      0.000 3/10206  700.0      0.244      0.244      0.106      0.000     -0.000      0.000 3/10206  710.0      0.241      0.241      0.105      0.000     -0.000      0.000 3/10206  720.0      0.238      0.238      0.103      0.000     -0.000      0.000 3/10206  730.0      0.234      0.234      0.102      0.000     -0.000      0.000 3/10206  740.0      0.231      0.231      0.100      0.000     -0.000      0.000 3/10206  750.0      0.228      0.228      0.099      0.000     -0.000      0.000 3/10206  760.0      0.225      0.225      0.098      0.000     -0.000      0.000 3/10206  770.0      0.222      0.222      0.096      0.000     -0.000      0.000 3/10206  780.0      0.219      0.219      0.095      0.000     -0.000      0.000 3/10206  790.0      0.217      0.217      0.094      0.000     -0.000      0.000 3/10206  800.0      0.214      0.214      0.093      0.000     -0.000      0.000 3/10206  810.0      0.211      0.211      0.092      0.000     -0.000      0.000 3/10206  820.0      0.209      0.209      0.090      0.000     -0.000      0.000 3/10206  830.0      0.206      0.206      0.089      0.000     -0.000      0.000 3/10206  840.0      0.204      0.204      0.088      0.000     -0.000      0.000 3/10206  850.0      0.201      0.201      0.087      0.000     -0.000      0.000 3/10206  860.0      0.199      0.199      0.086      0.000     -0.000      0.000 3/10206  870.0      0.197      0.197      0.085      0.000     -0.000      0.000 3/10206  880.0      0.195      0.195      0.084      0.000     -0.000      0.000 3/10206  890.0      0.192      0.192      0.083      0.000     -0.000      0.000 3/10206  900.0      0.190      0.190      0.082      0.000     -0.000      0.000 3/10206  910.0      0.188      0.188      0.082      0.000     -0.000      0.000 3/10206  920.0      0.186      0.186      0.081      0.000     -0.000      0.000 3/10206  930.0      0.184      0.184      0.080      0.000     -0.000      0.000 3/10206  940.0      0.182      0.182      0.079      0.000     -0.000      0.000 3/10206  950.0      0.180      0.180      0.078      0.000     -0.000      0.000 3/10206  960.0      0.178      0.178      0.077      0.000     -0.000      0.000 3/10206  970.0      0.177      0.177      0.076      0.000     -0.000      0.000 3/10206  980.0      0.175      0.175      0.076      0.000     -0.000      0.000 3/10206  990.0      0.173      0.173      0.075      0.000     -0.000      0.000 3/10206 1000.0      0.171      0.171      0.074      0.000     -0.000      0.000 3/10206Thermal conductivity related properties were written into "kappa-m3918.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-09 00:53:59]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|