# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/d68d33a6-c083-44a9-977c-c97c12015d8d/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 PrBi2BrO4 / P4/mmm (123) / materials id 551900](https://mdr.nims.go.jp/datasets/45182e16-56ee-4298-b170-14506dd581b2)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-08 10:51:18]-------------------------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 2]  Primitive matrix:    [1. 0. 0.]    [0. 1. 0.]    [0. 0. 1.]Spacegroup: P4/mmm (123)Number of symmetry operations in supercell: 288------------------------------ primitive cell ------------------------------Lattice vectors:  a    3.941866730000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.941866730000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.188490310000001Atomic positions (fractional):   *1 O   0.50000000000000  0.00000000000000  0.16017525025169  15.999    2 O   0.50000000000000  0.00000000000000  0.83982474974831  15.999    3 O   0.00000000000000  0.50000000000000  0.83982474974831  15.999    4 O   0.00000000000000  0.50000000000000  0.16017525025169  15.999   *5 Bi  0.50000000000000  0.50000000000000  0.27816677469145 208.980    6 Bi  0.50000000000000  0.50000000000000  0.72183322530855 208.980   *7 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908   *8 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904-------------------------------- unit cell ---------------------------------Lattice vectors:  a    3.941866730000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.941866730000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.188490310000001Atomic positions (fractional):   *1 O   0.50000000000000  0.00000000000000  0.16017525025169  15.999 > 1    2 O   0.50000000000000  0.00000000000000  0.83982474974831  15.999 > 2    3 O   0.00000000000000  0.50000000000000  0.83982474974831  15.999 > 3    4 O   0.00000000000000  0.50000000000000  0.16017525025169  15.999 > 4   *5 Bi  0.50000000000000  0.50000000000000  0.27816677469145 208.980 > 5    6 Bi  0.50000000000000  0.50000000000000  0.72183322530855 208.980 > 6   *7 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908 > 7   *8 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 8-------------------------------- super cell --------------------------------Lattice vectors:  a   11.825600189999999    0.000000000000000    0.000000000000000  b    0.000000000000000   11.825600189999999    0.000000000000000  c    0.000000000000000    0.000000000000000   18.376980620000001Atomic positions (fractional):   *1 O   0.16666666666667  0.00000000000000  0.08008762512584  15.999 > 1    2 O   0.50000000000000  0.00000000000000  0.08008762512584  15.999 > 1    3 O   0.83333333333333  0.00000000000000  0.08008762512584  15.999 > 1    4 O   0.16666666666667  0.33333333333333  0.08008762512584  15.999 > 1    5 O   0.50000000000000  0.33333333333333  0.08008762512584  15.999 > 1    6 O   0.83333333333333  0.33333333333333  0.08008762512584  15.999 > 1    7 O   0.16666666666667  0.66666666666667  0.08008762512584  15.999 > 1    8 O   0.50000000000000  0.66666666666667  0.08008762512584  15.999 > 1    9 O   0.83333333333333  0.66666666666667  0.08008762512584  15.999 > 1   10 O   0.16666666666667  0.00000000000000  0.58008762512584  15.999 > 1   11 O   0.50000000000000  0.00000000000000  0.58008762512584  15.999 > 1   12 O   0.83333333333333  0.00000000000000  0.58008762512584  15.999 > 1   13 O   0.16666666666667  0.33333333333333  0.58008762512584  15.999 > 1   14 O   0.50000000000000  0.33333333333333  0.58008762512584  15.999 > 1   15 O   0.83333333333333  0.33333333333333  0.58008762512584  15.999 > 1   16 O   0.16666666666667  0.66666666666667  0.58008762512584  15.999 > 1   17 O   0.50000000000000  0.66666666666667  0.58008762512584  15.999 > 1   18 O   0.83333333333333  0.66666666666667  0.58008762512584  15.999 > 1   19 O   0.16666666666667  0.00000000000000  0.41991237487416  15.999 > 2   20 O   0.50000000000000  0.00000000000000  0.41991237487416  15.999 > 2   21 O   0.83333333333333  0.00000000000000  0.41991237487416  15.999 > 2   22 O   0.16666666666667  0.33333333333333  0.41991237487416  15.999 > 2   23 O   0.50000000000000  0.33333333333333  0.41991237487416  15.999 > 2   24 O   0.83333333333333  0.33333333333333  0.41991237487416  15.999 > 2   25 O   0.16666666666667  0.66666666666667  0.41991237487416  15.999 > 2   26 O   0.50000000000000  0.66666666666667  0.41991237487416  15.999 > 2   27 O   0.83333333333333  0.66666666666667  0.41991237487416  15.999 > 2   28 O   0.16666666666667  0.00000000000000  0.91991237487416  15.999 > 2   29 O   0.50000000000000  0.00000000000000  0.91991237487416  15.999 > 2   30 O   0.83333333333333  0.00000000000000  0.91991237487416  15.999 > 2   31 O   0.16666666666667  0.33333333333333  0.91991237487416  15.999 > 2   32 O   0.50000000000000  0.33333333333333  0.91991237487416  15.999 > 2   33 O   0.83333333333333  0.33333333333333  0.91991237487416  15.999 > 2   34 O   0.16666666666667  0.66666666666667  0.91991237487416  15.999 > 2   35 O   0.50000000000000  0.66666666666667  0.91991237487416  15.999 > 2   36 O   0.83333333333333  0.66666666666667  0.91991237487416  15.999 > 2   37 O   0.00000000000000  0.16666666666667  0.41991237487416  15.999 > 3   38 O   0.33333333333333  0.16666666666667  0.41991237487416  15.999 > 3   39 O   0.66666666666667  0.16666666666667  0.41991237487416  15.999 > 3   40 O   0.00000000000000  0.50000000000000  0.41991237487416  15.999 > 3   41 O   0.33333333333333  0.50000000000000  0.41991237487416  15.999 > 3   42 O   0.66666666666667  0.50000000000000  0.41991237487416  15.999 > 3   43 O   0.00000000000000  0.83333333333333  0.41991237487416  15.999 > 3   44 O   0.33333333333333  0.83333333333333  0.41991237487416  15.999 > 3   45 O   0.66666666666667  0.83333333333333  0.41991237487416  15.999 > 3   46 O   0.00000000000000  0.16666666666667  0.91991237487416  15.999 > 3   47 O   0.33333333333333  0.16666666666667  0.91991237487416  15.999 > 3   48 O   0.66666666666667  0.16666666666667  0.91991237487416  15.999 > 3   49 O   0.00000000000000  0.50000000000000  0.91991237487416  15.999 > 3   50 O   0.33333333333333  0.50000000000000  0.91991237487416  15.999 > 3   51 O   0.66666666666667  0.50000000000000  0.91991237487416  15.999 > 3   52 O   0.00000000000000  0.83333333333333  0.91991237487416  15.999 > 3   53 O   0.33333333333333  0.83333333333333  0.91991237487416  15.999 > 3   54 O   0.66666666666667  0.83333333333333  0.91991237487416  15.999 > 3   55 O   0.00000000000000  0.16666666666667  0.08008762512584  15.999 > 4   56 O   0.33333333333333  0.16666666666667  0.08008762512584  15.999 > 4   57 O   0.66666666666667  0.16666666666667  0.08008762512584  15.999 > 4   58 O   0.00000000000000  0.50000000000000  0.08008762512584  15.999 > 4   59 O   0.33333333333333  0.50000000000000  0.08008762512584  15.999 > 4   60 O   0.66666666666667  0.50000000000000  0.08008762512584  15.999 > 4   61 O   0.00000000000000  0.83333333333333  0.08008762512584  15.999 > 4   62 O   0.33333333333333  0.83333333333333  0.08008762512584  15.999 > 4   63 O   0.66666666666667  0.83333333333333  0.08008762512584  15.999 > 4   64 O   0.00000000000000  0.16666666666667  0.58008762512584  15.999 > 4   65 O   0.33333333333333  0.16666666666667  0.58008762512584  15.999 > 4   66 O   0.66666666666667  0.16666666666667  0.58008762512584  15.999 > 4   67 O   0.00000000000000  0.50000000000000  0.58008762512584  15.999 > 4   68 O   0.33333333333333  0.50000000000000  0.58008762512584  15.999 > 4   69 O   0.66666666666667  0.50000000000000  0.58008762512584  15.999 > 4   70 O   0.00000000000000  0.83333333333333  0.58008762512584  15.999 > 4   71 O   0.33333333333333  0.83333333333333  0.58008762512584  15.999 > 4   72 O   0.66666666666667  0.83333333333333  0.58008762512584  15.999 > 4  *73 Bi  0.16666666666667  0.16666666666667  0.13908338734572 208.980 > 5   74 Bi  0.50000000000000  0.16666666666667  0.13908338734572 208.980 > 5   75 Bi  0.83333333333333  0.16666666666667  0.13908338734572 208.980 > 5   76 Bi  0.16666666666667  0.50000000000000  0.13908338734572 208.980 > 5   77 Bi  0.50000000000000  0.50000000000000  0.13908338734572 208.980 > 5   78 Bi  0.83333333333333  0.50000000000000  0.13908338734572 208.980 > 5   79 Bi  0.16666666666667  0.83333333333333  0.13908338734572 208.980 > 5   80 Bi  0.50000000000000  0.83333333333333  0.13908338734572 208.980 > 5   81 Bi  0.83333333333333  0.83333333333333  0.13908338734572 208.980 > 5   82 Bi  0.16666666666667  0.16666666666667  0.63908338734572 208.980 > 5   83 Bi  0.50000000000000  0.16666666666667  0.63908338734572 208.980 > 5   84 Bi  0.83333333333333  0.16666666666667  0.63908338734572 208.980 > 5   85 Bi  0.16666666666667  0.50000000000000  0.63908338734572 208.980 > 5   86 Bi  0.50000000000000  0.50000000000000  0.63908338734572 208.980 > 5   87 Bi  0.83333333333333  0.50000000000000  0.63908338734572 208.980 > 5   88 Bi  0.16666666666667  0.83333333333333  0.63908338734572 208.980 > 5   89 Bi  0.50000000000000  0.83333333333333  0.63908338734572 208.980 > 5   90 Bi  0.83333333333333  0.83333333333333  0.63908338734572 208.980 > 5   91 Bi  0.16666666666667  0.16666666666667  0.36091661265428 208.980 > 6   92 Bi  0.50000000000000  0.16666666666667  0.36091661265428 208.980 > 6   93 Bi  0.83333333333333  0.16666666666667  0.36091661265428 208.980 > 6   94 Bi  0.16666666666667  0.50000000000000  0.36091661265428 208.980 > 6   95 Bi  0.50000000000000  0.50000000000000  0.36091661265428 208.980 > 6   96 Bi  0.83333333333333  0.50000000000000  0.36091661265428 208.980 > 6   97 Bi  0.16666666666667  0.83333333333333  0.36091661265428 208.980 > 6   98 Bi  0.50000000000000  0.83333333333333  0.36091661265428 208.980 > 6   99 Bi  0.83333333333333  0.83333333333333  0.36091661265428 208.980 > 6  100 Bi  0.16666666666667  0.16666666666667  0.86091661265428 208.980 > 6  101 Bi  0.50000000000000  0.16666666666667  0.86091661265428 208.980 > 6  102 Bi  0.83333333333333  0.16666666666667  0.86091661265428 208.980 > 6  103 Bi  0.16666666666667  0.50000000000000  0.86091661265428 208.980 > 6  104 Bi  0.50000000000000  0.50000000000000  0.86091661265428 208.980 > 6  105 Bi  0.83333333333333  0.50000000000000  0.86091661265428 208.980 > 6  106 Bi  0.16666666666667  0.83333333333333  0.86091661265428 208.980 > 6  107 Bi  0.50000000000000  0.83333333333333  0.86091661265428 208.980 > 6  108 Bi  0.83333333333333  0.83333333333333  0.86091661265428 208.980 > 6 *109 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908 > 7  110 Pr  0.33333333333333  0.00000000000000  0.00000000000000 140.908 > 7  111 Pr  0.66666666666667  0.00000000000000  0.00000000000000 140.908 > 7  112 Pr  0.00000000000000  0.33333333333333  0.00000000000000 140.908 > 7  113 Pr  0.33333333333333  0.33333333333333  0.00000000000000 140.908 > 7  114 Pr  0.66666666666667  0.33333333333333  0.00000000000000 140.908 > 7  115 Pr  0.00000000000000  0.66666666666667  0.00000000000000 140.908 > 7  116 Pr  0.33333333333333  0.66666666666667  0.00000000000000 140.908 > 7  117 Pr  0.66666666666667  0.66666666666667  0.00000000000000 140.908 > 7  118 Pr  0.00000000000000  0.00000000000000  0.50000000000000 140.908 > 7  119 Pr  0.33333333333333  0.00000000000000  0.50000000000000 140.908 > 7  120 Pr  0.66666666666667  0.00000000000000  0.50000000000000 140.908 > 7  121 Pr  0.00000000000000  0.33333333333333  0.50000000000000 140.908 > 7  122 Pr  0.33333333333333  0.33333333333333  0.50000000000000 140.908 > 7  123 Pr  0.66666666666667  0.33333333333333  0.50000000000000 140.908 > 7  124 Pr  0.00000000000000  0.66666666666667  0.50000000000000 140.908 > 7  125 Pr  0.33333333333333  0.66666666666667  0.50000000000000 140.908 > 7  126 Pr  0.66666666666667  0.66666666666667  0.50000000000000 140.908 > 7 *127 Br  0.00000000000000  0.00000000000000  0.25000000000000  79.904 > 8  128 Br  0.33333333333333  0.00000000000000  0.25000000000000  79.904 > 8  129 Br  0.66666666666667  0.00000000000000  0.25000000000000  79.904 > 8  130 Br  0.00000000000000  0.33333333333333  0.25000000000000  79.904 > 8  131 Br  0.33333333333333  0.33333333333333  0.25000000000000  79.904 > 8  132 Br  0.66666666666667  0.33333333333333  0.25000000000000  79.904 > 8  133 Br  0.00000000000000  0.66666666666667  0.25000000000000  79.904 > 8  134 Br  0.33333333333333  0.66666666666667  0.25000000000000  79.904 > 8  135 Br  0.66666666666667  0.66666666666667  0.25000000000000  79.904 > 8  136 Br  0.00000000000000  0.00000000000000  0.75000000000000  79.904 > 8  137 Br  0.33333333333333  0.00000000000000  0.75000000000000  79.904 > 8  138 Br  0.66666666666667  0.00000000000000  0.75000000000000  79.904 > 8  139 Br  0.00000000000000  0.33333333333333  0.75000000000000  79.904 > 8  140 Br  0.33333333333333  0.33333333333333  0.75000000000000  79.904 > 8  141 Br  0.66666666666667  0.33333333333333  0.75000000000000  79.904 > 8  142 Br  0.00000000000000  0.66666666666667  0.75000000000000  79.904 > 8  143 Br  0.33333333333333  0.66666666666667  0.75000000000000  79.904 > 8  144 Br  0.66666666666667  0.66666666666667  0.75000000000000  79.904 > 8----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            6.8146433    0.0000000    0.0000000            0.0000000    6.8146433    0.0000000            0.0000000    0.0000000    5.6056542-------------------------- Born effective charges --------------------------    1 O    -2.1895334    0.0000000    0.0000000            0.0000000   -4.0265966    0.0000000            0.0000000    0.0000000   -2.3026890    2 O    -2.1895334    0.0000000    0.0000000            0.0000000   -4.0265966    0.0000000            0.0000000    0.0000000   -2.3026890    3 O    -4.0265966    0.0000000    0.0000000            0.0000000   -2.1895334    0.0000000            0.0000000    0.0000000   -2.3026890    4 O    -4.0265966    0.0000000    0.0000000            0.0000000   -2.1895334    0.0000000            0.0000000    0.0000000   -2.3026890    5 Bi    5.3589907    0.0000000    0.0000000            0.0000000    5.3589907    0.0000000            0.0000000    0.0000000    2.6408645    6 Bi    5.3589907    0.0000000    0.0000000            0.0000000    5.3589907    0.0000000            0.0000000    0.0000000    2.6408645    7 Pr    3.7895244    0.0000000    0.0000000            0.0000000    3.7895244    0.0000000            0.0000000    0.0000000    5.6329384    8 Br   -2.0752457    0.0000000    0.0000000            0.0000000   -2.0752457    0.0000000            0.0000000    0.0000000   -1.7039114----------------------------------------------------------------------------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: 432/432Permutation basis: 5784/5784Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 428Number of blocks in projector: 428Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 8--- Eigsh_solver_block: 1 / 8 ---Block_size: 164Use standard eigh solver.--- Eigsh_solver_block: 2 / 8 ---Block_size: 124Use standard eigh solver.--- Eigsh_solver_block: 3 / 8 ---Block_size: 99Use standard eigh solver.--- Eigsh_solver_block: 4 / 8 ---Block_size: 37Use standard eigh solver.--- Eigsh_solver_block: 5 / 8 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 6 / 8 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 7 / 8 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 8 / 8 ---Block_size: 1Use standard eigh solver.Tree of FC basis block matrices:- (428, 419), data: False|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (37, 37), data: True|-- (99, 95), data: True|-- (124, 124), data: True|-- (164, 159), data: True-----Solver_atoms: 1 -- 144 / 144Time (Solver_compr_matrix_reshape): 0.008Solver_block: 80 / 80 - Time: 0.349Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.361--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 432/432Permutation basis: 5784/5784Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 428Number of blocks in projector: 428Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 8--- Eigsh_solver_block: 1 / 8 ---Block_size: 164Use standard eigh solver.--- Eigsh_solver_block: 2 / 8 ---Block_size: 124Use standard eigh solver.--- Eigsh_solver_block: 3 / 8 ---Block_size: 99Use standard eigh solver.--- Eigsh_solver_block: 4 / 8 ---Block_size: 37Use standard eigh solver.--- Eigsh_solver_block: 5 / 8 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 6 / 8 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 7 / 8 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 8 / 8 ---Block_size: 1Use standard eigh solver.Tree of FC basis block matrices:- (428, 419), data: False|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (37, 37), data: True|-- (99, 95), data: True|-- (124, 124), data: True|-- (164, 159), data: TrueMax drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-08 10:51:23]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 10:51:24]-------------------------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 2]Primitive matrix:  [1. 0. 0.]  [0. 1. 0.]  [0. 0. 1.]Spacegroup: P4/mmm (123)------------------------------ primitive cell ------------------------------Lattice vectors:  a    3.941866730000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.941866730000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.188490310000001Atomic positions (fractional):    1 O   0.50000000000000  0.00000000000000  0.16017525025169  15.999    2 O   0.50000000000000  0.00000000000000  0.83982474974831  15.999    3 O   0.00000000000000  0.50000000000000  0.83982474974831  15.999    4 O   0.00000000000000  0.50000000000000  0.16017525025169  15.999    5 Bi  0.50000000000000  0.50000000000000  0.27816677469145 208.980    6 Bi  0.50000000000000  0.50000000000000  0.72183322530855 208.980    7 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908    8 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904-------------------------------- supercell ---------------------------------Lattice vectors:  a   11.825600189999999    0.000000000000000    0.000000000000000  b    0.000000000000000   11.825600189999999    0.000000000000000  c    0.000000000000000    0.000000000000000   18.376980620000001Atomic positions (fractional):    1 O   0.16666666666667  0.00000000000000  0.08008762512584  15.999 > 1    2 O   0.50000000000000  0.00000000000000  0.08008762512584  15.999 > 1    3 O   0.83333333333333  0.00000000000000  0.08008762512584  15.999 > 1    4 O   0.16666666666667  0.33333333333333  0.08008762512584  15.999 > 1    5 O   0.50000000000000  0.33333333333333  0.08008762512584  15.999 > 1    6 O   0.83333333333333  0.33333333333333  0.08008762512584  15.999 > 1    7 O   0.16666666666667  0.66666666666667  0.08008762512584  15.999 > 1    8 O   0.50000000000000  0.66666666666667  0.08008762512584  15.999 > 1    9 O   0.83333333333333  0.66666666666667  0.08008762512584  15.999 > 1   10 O   0.16666666666667  0.00000000000000  0.58008762512584  15.999 > 1   11 O   0.50000000000000  0.00000000000000  0.58008762512584  15.999 > 1   12 O   0.83333333333333  0.00000000000000  0.58008762512584  15.999 > 1   13 O   0.16666666666667  0.33333333333333  0.58008762512584  15.999 > 1   14 O   0.50000000000000  0.33333333333333  0.58008762512584  15.999 > 1   15 O   0.83333333333333  0.33333333333333  0.58008762512584  15.999 > 1   16 O   0.16666666666667  0.66666666666667  0.58008762512584  15.999 > 1   17 O   0.50000000000000  0.66666666666667  0.58008762512584  15.999 > 1   18 O   0.83333333333333  0.66666666666667  0.58008762512584  15.999 > 1   19 O   0.16666666666667  0.00000000000000  0.41991237487416  15.999 > 19   20 O   0.50000000000000  0.00000000000000  0.41991237487416  15.999 > 19   21 O   0.83333333333333  0.00000000000000  0.41991237487416  15.999 > 19   22 O   0.16666666666667  0.33333333333333  0.41991237487416  15.999 > 19   23 O   0.50000000000000  0.33333333333333  0.41991237487416  15.999 > 19   24 O   0.83333333333333  0.33333333333333  0.41991237487416  15.999 > 19   25 O   0.16666666666667  0.66666666666667  0.41991237487416  15.999 > 19   26 O   0.50000000000000  0.66666666666667  0.41991237487416  15.999 > 19   27 O   0.83333333333333  0.66666666666667  0.41991237487416  15.999 > 19   28 O   0.16666666666667  0.00000000000000  0.91991237487416  15.999 > 19   29 O   0.50000000000000  0.00000000000000  0.91991237487416  15.999 > 19   30 O   0.83333333333333  0.00000000000000  0.91991237487416  15.999 > 19   31 O   0.16666666666667  0.33333333333333  0.91991237487416  15.999 > 19   32 O   0.50000000000000  0.33333333333333  0.91991237487416  15.999 > 19   33 O   0.83333333333333  0.33333333333333  0.91991237487416  15.999 > 19   34 O   0.16666666666667  0.66666666666667  0.91991237487416  15.999 > 19   35 O   0.50000000000000  0.66666666666667  0.91991237487416  15.999 > 19   36 O   0.83333333333333  0.66666666666667  0.91991237487416  15.999 > 19   37 O   0.00000000000000  0.16666666666667  0.41991237487416  15.999 > 37   38 O   0.33333333333333  0.16666666666667  0.41991237487416  15.999 > 37   39 O   0.66666666666667  0.16666666666667  0.41991237487416  15.999 > 37   40 O   0.00000000000000  0.50000000000000  0.41991237487416  15.999 > 37   41 O   0.33333333333333  0.50000000000000  0.41991237487416  15.999 > 37   42 O   0.66666666666667  0.50000000000000  0.41991237487416  15.999 > 37   43 O   0.00000000000000  0.83333333333333  0.41991237487416  15.999 > 37   44 O   0.33333333333333  0.83333333333333  0.41991237487416  15.999 > 37   45 O   0.66666666666667  0.83333333333333  0.41991237487416  15.999 > 37   46 O   0.00000000000000  0.16666666666667  0.91991237487416  15.999 > 37   47 O   0.33333333333333  0.16666666666667  0.91991237487416  15.999 > 37   48 O   0.66666666666667  0.16666666666667  0.91991237487416  15.999 > 37   49 O   0.00000000000000  0.50000000000000  0.91991237487416  15.999 > 37   50 O   0.33333333333333  0.50000000000000  0.91991237487416  15.999 > 37   51 O   0.66666666666667  0.50000000000000  0.91991237487416  15.999 > 37   52 O   0.00000000000000  0.83333333333333  0.91991237487416  15.999 > 37   53 O   0.33333333333333  0.83333333333333  0.91991237487416  15.999 > 37   54 O   0.66666666666667  0.83333333333333  0.91991237487416  15.999 > 37   55 O   0.00000000000000  0.16666666666667  0.08008762512584  15.999 > 55   56 O   0.33333333333333  0.16666666666667  0.08008762512584  15.999 > 55   57 O   0.66666666666667  0.16666666666667  0.08008762512584  15.999 > 55   58 O   0.00000000000000  0.50000000000000  0.08008762512584  15.999 > 55   59 O   0.33333333333333  0.50000000000000  0.08008762512584  15.999 > 55   60 O   0.66666666666667  0.50000000000000  0.08008762512584  15.999 > 55   61 O   0.00000000000000  0.83333333333333  0.08008762512584  15.999 > 55   62 O   0.33333333333333  0.83333333333333  0.08008762512584  15.999 > 55   63 O   0.66666666666667  0.83333333333333  0.08008762512584  15.999 > 55   64 O   0.00000000000000  0.16666666666667  0.58008762512584  15.999 > 55   65 O   0.33333333333333  0.16666666666667  0.58008762512584  15.999 > 55   66 O   0.66666666666667  0.16666666666667  0.58008762512584  15.999 > 55   67 O   0.00000000000000  0.50000000000000  0.58008762512584  15.999 > 55   68 O   0.33333333333333  0.50000000000000  0.58008762512584  15.999 > 55   69 O   0.66666666666667  0.50000000000000  0.58008762512584  15.999 > 55   70 O   0.00000000000000  0.83333333333333  0.58008762512584  15.999 > 55   71 O   0.33333333333333  0.83333333333333  0.58008762512584  15.999 > 55   72 O   0.66666666666667  0.83333333333333  0.58008762512584  15.999 > 55   73 Bi  0.16666666666667  0.16666666666667  0.13908338734572 208.980 > 73   74 Bi  0.50000000000000  0.16666666666667  0.13908338734572 208.980 > 73   75 Bi  0.83333333333333  0.16666666666667  0.13908338734572 208.980 > 73   76 Bi  0.16666666666667  0.50000000000000  0.13908338734572 208.980 > 73   77 Bi  0.50000000000000  0.50000000000000  0.13908338734572 208.980 > 73   78 Bi  0.83333333333333  0.50000000000000  0.13908338734572 208.980 > 73   79 Bi  0.16666666666667  0.83333333333333  0.13908338734572 208.980 > 73   80 Bi  0.50000000000000  0.83333333333333  0.13908338734572 208.980 > 73   81 Bi  0.83333333333333  0.83333333333333  0.13908338734572 208.980 > 73   82 Bi  0.16666666666667  0.16666666666667  0.63908338734572 208.980 > 73   83 Bi  0.50000000000000  0.16666666666667  0.63908338734572 208.980 > 73   84 Bi  0.83333333333333  0.16666666666667  0.63908338734572 208.980 > 73   85 Bi  0.16666666666667  0.50000000000000  0.63908338734572 208.980 > 73   86 Bi  0.50000000000000  0.50000000000000  0.63908338734572 208.980 > 73   87 Bi  0.83333333333333  0.50000000000000  0.63908338734572 208.980 > 73   88 Bi  0.16666666666667  0.83333333333333  0.63908338734572 208.980 > 73   89 Bi  0.50000000000000  0.83333333333333  0.63908338734572 208.980 > 73   90 Bi  0.83333333333333  0.83333333333333  0.63908338734572 208.980 > 73   91 Bi  0.16666666666667  0.16666666666667  0.36091661265428 208.980 > 91   92 Bi  0.50000000000000  0.16666666666667  0.36091661265428 208.980 > 91   93 Bi  0.83333333333333  0.16666666666667  0.36091661265428 208.980 > 91   94 Bi  0.16666666666667  0.50000000000000  0.36091661265428 208.980 > 91   95 Bi  0.50000000000000  0.50000000000000  0.36091661265428 208.980 > 91   96 Bi  0.83333333333333  0.50000000000000  0.36091661265428 208.980 > 91   97 Bi  0.16666666666667  0.83333333333333  0.36091661265428 208.980 > 91   98 Bi  0.50000000000000  0.83333333333333  0.36091661265428 208.980 > 91   99 Bi  0.83333333333333  0.83333333333333  0.36091661265428 208.980 > 91  100 Bi  0.16666666666667  0.16666666666667  0.86091661265428 208.980 > 91  101 Bi  0.50000000000000  0.16666666666667  0.86091661265428 208.980 > 91  102 Bi  0.83333333333333  0.16666666666667  0.86091661265428 208.980 > 91  103 Bi  0.16666666666667  0.50000000000000  0.86091661265428 208.980 > 91  104 Bi  0.50000000000000  0.50000000000000  0.86091661265428 208.980 > 91  105 Bi  0.83333333333333  0.50000000000000  0.86091661265428 208.980 > 91  106 Bi  0.16666666666667  0.83333333333333  0.86091661265428 208.980 > 91  107 Bi  0.50000000000000  0.83333333333333  0.86091661265428 208.980 > 91  108 Bi  0.83333333333333  0.83333333333333  0.86091661265428 208.980 > 91  109 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908 > 109  110 Pr  0.33333333333333  0.00000000000000  0.00000000000000 140.908 > 109  111 Pr  0.66666666666667  0.00000000000000  0.00000000000000 140.908 > 109  112 Pr  0.00000000000000  0.33333333333333  0.00000000000000 140.908 > 109  113 Pr  0.33333333333333  0.33333333333333  0.00000000000000 140.908 > 109  114 Pr  0.66666666666667  0.33333333333333  0.00000000000000 140.908 > 109  115 Pr  0.00000000000000  0.66666666666667  0.00000000000000 140.908 > 109  116 Pr  0.33333333333333  0.66666666666667  0.00000000000000 140.908 > 109  117 Pr  0.66666666666667  0.66666666666667  0.00000000000000 140.908 > 109  118 Pr  0.00000000000000  0.00000000000000  0.50000000000000 140.908 > 109  119 Pr  0.33333333333333  0.00000000000000  0.50000000000000 140.908 > 109  120 Pr  0.66666666666667  0.00000000000000  0.50000000000000 140.908 > 109  121 Pr  0.00000000000000  0.33333333333333  0.50000000000000 140.908 > 109  122 Pr  0.33333333333333  0.33333333333333  0.50000000000000 140.908 > 109  123 Pr  0.66666666666667  0.33333333333333  0.50000000000000 140.908 > 109  124 Pr  0.00000000000000  0.66666666666667  0.50000000000000 140.908 > 109  125 Pr  0.33333333333333  0.66666666666667  0.50000000000000 140.908 > 109  126 Pr  0.66666666666667  0.66666666666667  0.50000000000000 140.908 > 109  127 Br  0.00000000000000  0.00000000000000  0.25000000000000  79.904 > 127  128 Br  0.33333333333333  0.00000000000000  0.25000000000000  79.904 > 127  129 Br  0.66666666666667  0.00000000000000  0.25000000000000  79.904 > 127  130 Br  0.00000000000000  0.33333333333333  0.25000000000000  79.904 > 127  131 Br  0.33333333333333  0.33333333333333  0.25000000000000  79.904 > 127  132 Br  0.66666666666667  0.33333333333333  0.25000000000000  79.904 > 127  133 Br  0.00000000000000  0.66666666666667  0.25000000000000  79.904 > 127  134 Br  0.33333333333333  0.66666666666667  0.25000000000000  79.904 > 127  135 Br  0.66666666666667  0.66666666666667  0.25000000000000  79.904 > 127  136 Br  0.00000000000000  0.00000000000000  0.75000000000000  79.904 > 127  137 Br  0.33333333333333  0.00000000000000  0.75000000000000  79.904 > 127  138 Br  0.66666666666667  0.00000000000000  0.75000000000000  79.904 > 127  139 Br  0.00000000000000  0.33333333333333  0.75000000000000  79.904 > 127  140 Br  0.33333333333333  0.33333333333333  0.75000000000000  79.904 > 127  141 Br  0.66666666666667  0.33333333333333  0.75000000000000  79.904 > 127  142 Br  0.00000000000000  0.66666666666667  0.75000000000000  79.904 > 127  143 Br  0.33333333333333  0.66666666666667  0.75000000000000  79.904 > 127  144 Br  0.66666666666667  0.66666666666667  0.75000000000000  79.904 > 127----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            6.8146433    0.0000000    0.0000000            0.0000000    6.8146433    0.0000000            0.0000000    0.0000000    5.6056542-------------------------- Born effective charges --------------------------    1 O    -2.1895334    0.0000000    0.0000000            0.0000000   -4.0265966    0.0000000            0.0000000    0.0000000   -2.3026890    2 O    -2.1895334    0.0000000    0.0000000            0.0000000   -4.0265966    0.0000000            0.0000000    0.0000000   -2.3026890    3 O    -4.0265966    0.0000000    0.0000000            0.0000000   -2.1895334    0.0000000            0.0000000    0.0000000   -2.3026890    4 O    -4.0265966    0.0000000    0.0000000            0.0000000   -2.1895334    0.0000000            0.0000000    0.0000000   -2.3026890    5 Bi    5.3589907    0.0000000    0.0000000            0.0000000    5.3589907    0.0000000            0.0000000    0.0000000    2.6408645    6 Bi    5.3589907    0.0000000    0.0000000            0.0000000    5.3589907    0.0000000            0.0000000    0.0000000    2.6408645    7 Pr    3.7895244    0.0000000    0.0000000            0.0000000    3.7895244    0.0000000            0.0000000    0.0000000    5.6329384    8 Br   -2.0752457    0.0000000    0.0000000            0.0000000   -2.0752457    0.0000000            0.0000000    0.0000000   -1.7039114----------------------------------------------------------------------------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[ 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]Computing fc3[ 127, 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 (zxx) 0.00000000 (zxx) 0.00000000 (zxx)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-08 10:51:33]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 10:51:33]-------------------------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 2]Primitive matrix:  [1. 0. 0.]  [0. 1. 0.]  [0. 0. 1.]Spacegroup: P4/mmm (123)------------------------------ primitive cell ------------------------------Lattice vectors:  a    3.941866730000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.941866730000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.188490310000001Atomic positions (fractional):    1 O   0.50000000000000  0.00000000000000  0.16017525025169  15.999    2 O   0.50000000000000  0.00000000000000  0.83982474974831  15.999    3 O   0.00000000000000  0.50000000000000  0.83982474974831  15.999    4 O   0.00000000000000  0.50000000000000  0.16017525025169  15.999    5 Bi  0.50000000000000  0.50000000000000  0.27816677469145 208.980    6 Bi  0.50000000000000  0.50000000000000  0.72183322530855 208.980    7 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908    8 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904-------------------------------- supercell ---------------------------------Lattice vectors:  a   11.825600189999999    0.000000000000000    0.000000000000000  b    0.000000000000000   11.825600189999999    0.000000000000000  c    0.000000000000000    0.000000000000000   18.376980620000001Atomic positions (fractional):    1 O   0.16666666666667  0.00000000000000  0.08008762512584  15.999 > 1    2 O   0.50000000000000  0.00000000000000  0.08008762512584  15.999 > 1    3 O   0.83333333333333  0.00000000000000  0.08008762512584  15.999 > 1    4 O   0.16666666666667  0.33333333333333  0.08008762512584  15.999 > 1    5 O   0.50000000000000  0.33333333333333  0.08008762512584  15.999 > 1    6 O   0.83333333333333  0.33333333333333  0.08008762512584  15.999 > 1    7 O   0.16666666666667  0.66666666666667  0.08008762512584  15.999 > 1    8 O   0.50000000000000  0.66666666666667  0.08008762512584  15.999 > 1    9 O   0.83333333333333  0.66666666666667  0.08008762512584  15.999 > 1   10 O   0.16666666666667  0.00000000000000  0.58008762512584  15.999 > 1   11 O   0.50000000000000  0.00000000000000  0.58008762512584  15.999 > 1   12 O   0.83333333333333  0.00000000000000  0.58008762512584  15.999 > 1   13 O   0.16666666666667  0.33333333333333  0.58008762512584  15.999 > 1   14 O   0.50000000000000  0.33333333333333  0.58008762512584  15.999 > 1   15 O   0.83333333333333  0.33333333333333  0.58008762512584  15.999 > 1   16 O   0.16666666666667  0.66666666666667  0.58008762512584  15.999 > 1   17 O   0.50000000000000  0.66666666666667  0.58008762512584  15.999 > 1   18 O   0.83333333333333  0.66666666666667  0.58008762512584  15.999 > 1   19 O   0.16666666666667  0.00000000000000  0.41991237487416  15.999 > 19   20 O   0.50000000000000  0.00000000000000  0.41991237487416  15.999 > 19   21 O   0.83333333333333  0.00000000000000  0.41991237487416  15.999 > 19   22 O   0.16666666666667  0.33333333333333  0.41991237487416  15.999 > 19   23 O   0.50000000000000  0.33333333333333  0.41991237487416  15.999 > 19   24 O   0.83333333333333  0.33333333333333  0.41991237487416  15.999 > 19   25 O   0.16666666666667  0.66666666666667  0.41991237487416  15.999 > 19   26 O   0.50000000000000  0.66666666666667  0.41991237487416  15.999 > 19   27 O   0.83333333333333  0.66666666666667  0.41991237487416  15.999 > 19   28 O   0.16666666666667  0.00000000000000  0.91991237487416  15.999 > 19   29 O   0.50000000000000  0.00000000000000  0.91991237487416  15.999 > 19   30 O   0.83333333333333  0.00000000000000  0.91991237487416  15.999 > 19   31 O   0.16666666666667  0.33333333333333  0.91991237487416  15.999 > 19   32 O   0.50000000000000  0.33333333333333  0.91991237487416  15.999 > 19   33 O   0.83333333333333  0.33333333333333  0.91991237487416  15.999 > 19   34 O   0.16666666666667  0.66666666666667  0.91991237487416  15.999 > 19   35 O   0.50000000000000  0.66666666666667  0.91991237487416  15.999 > 19   36 O   0.83333333333333  0.66666666666667  0.91991237487416  15.999 > 19   37 O   0.00000000000000  0.16666666666667  0.41991237487416  15.999 > 37   38 O   0.33333333333333  0.16666666666667  0.41991237487416  15.999 > 37   39 O   0.66666666666667  0.16666666666667  0.41991237487416  15.999 > 37   40 O   0.00000000000000  0.50000000000000  0.41991237487416  15.999 > 37   41 O   0.33333333333333  0.50000000000000  0.41991237487416  15.999 > 37   42 O   0.66666666666667  0.50000000000000  0.41991237487416  15.999 > 37   43 O   0.00000000000000  0.83333333333333  0.41991237487416  15.999 > 37   44 O   0.33333333333333  0.83333333333333  0.41991237487416  15.999 > 37   45 O   0.66666666666667  0.83333333333333  0.41991237487416  15.999 > 37   46 O   0.00000000000000  0.16666666666667  0.91991237487416  15.999 > 37   47 O   0.33333333333333  0.16666666666667  0.91991237487416  15.999 > 37   48 O   0.66666666666667  0.16666666666667  0.91991237487416  15.999 > 37   49 O   0.00000000000000  0.50000000000000  0.91991237487416  15.999 > 37   50 O   0.33333333333333  0.50000000000000  0.91991237487416  15.999 > 37   51 O   0.66666666666667  0.50000000000000  0.91991237487416  15.999 > 37   52 O   0.00000000000000  0.83333333333333  0.91991237487416  15.999 > 37   53 O   0.33333333333333  0.83333333333333  0.91991237487416  15.999 > 37   54 O   0.66666666666667  0.83333333333333  0.91991237487416  15.999 > 37   55 O   0.00000000000000  0.16666666666667  0.08008762512584  15.999 > 55   56 O   0.33333333333333  0.16666666666667  0.08008762512584  15.999 > 55   57 O   0.66666666666667  0.16666666666667  0.08008762512584  15.999 > 55   58 O   0.00000000000000  0.50000000000000  0.08008762512584  15.999 > 55   59 O   0.33333333333333  0.50000000000000  0.08008762512584  15.999 > 55   60 O   0.66666666666667  0.50000000000000  0.08008762512584  15.999 > 55   61 O   0.00000000000000  0.83333333333333  0.08008762512584  15.999 > 55   62 O   0.33333333333333  0.83333333333333  0.08008762512584  15.999 > 55   63 O   0.66666666666667  0.83333333333333  0.08008762512584  15.999 > 55   64 O   0.00000000000000  0.16666666666667  0.58008762512584  15.999 > 55   65 O   0.33333333333333  0.16666666666667  0.58008762512584  15.999 > 55   66 O   0.66666666666667  0.16666666666667  0.58008762512584  15.999 > 55   67 O   0.00000000000000  0.50000000000000  0.58008762512584  15.999 > 55   68 O   0.33333333333333  0.50000000000000  0.58008762512584  15.999 > 55   69 O   0.66666666666667  0.50000000000000  0.58008762512584  15.999 > 55   70 O   0.00000000000000  0.83333333333333  0.58008762512584  15.999 > 55   71 O   0.33333333333333  0.83333333333333  0.58008762512584  15.999 > 55   72 O   0.66666666666667  0.83333333333333  0.58008762512584  15.999 > 55   73 Bi  0.16666666666667  0.16666666666667  0.13908338734572 208.980 > 73   74 Bi  0.50000000000000  0.16666666666667  0.13908338734572 208.980 > 73   75 Bi  0.83333333333333  0.16666666666667  0.13908338734572 208.980 > 73   76 Bi  0.16666666666667  0.50000000000000  0.13908338734572 208.980 > 73   77 Bi  0.50000000000000  0.50000000000000  0.13908338734572 208.980 > 73   78 Bi  0.83333333333333  0.50000000000000  0.13908338734572 208.980 > 73   79 Bi  0.16666666666667  0.83333333333333  0.13908338734572 208.980 > 73   80 Bi  0.50000000000000  0.83333333333333  0.13908338734572 208.980 > 73   81 Bi  0.83333333333333  0.83333333333333  0.13908338734572 208.980 > 73   82 Bi  0.16666666666667  0.16666666666667  0.63908338734572 208.980 > 73   83 Bi  0.50000000000000  0.16666666666667  0.63908338734572 208.980 > 73   84 Bi  0.83333333333333  0.16666666666667  0.63908338734572 208.980 > 73   85 Bi  0.16666666666667  0.50000000000000  0.63908338734572 208.980 > 73   86 Bi  0.50000000000000  0.50000000000000  0.63908338734572 208.980 > 73   87 Bi  0.83333333333333  0.50000000000000  0.63908338734572 208.980 > 73   88 Bi  0.16666666666667  0.83333333333333  0.63908338734572 208.980 > 73   89 Bi  0.50000000000000  0.83333333333333  0.63908338734572 208.980 > 73   90 Bi  0.83333333333333  0.83333333333333  0.63908338734572 208.980 > 73   91 Bi  0.16666666666667  0.16666666666667  0.36091661265428 208.980 > 91   92 Bi  0.50000000000000  0.16666666666667  0.36091661265428 208.980 > 91   93 Bi  0.83333333333333  0.16666666666667  0.36091661265428 208.980 > 91   94 Bi  0.16666666666667  0.50000000000000  0.36091661265428 208.980 > 91   95 Bi  0.50000000000000  0.50000000000000  0.36091661265428 208.980 > 91   96 Bi  0.83333333333333  0.50000000000000  0.36091661265428 208.980 > 91   97 Bi  0.16666666666667  0.83333333333333  0.36091661265428 208.980 > 91   98 Bi  0.50000000000000  0.83333333333333  0.36091661265428 208.980 > 91   99 Bi  0.83333333333333  0.83333333333333  0.36091661265428 208.980 > 91  100 Bi  0.16666666666667  0.16666666666667  0.86091661265428 208.980 > 91  101 Bi  0.50000000000000  0.16666666666667  0.86091661265428 208.980 > 91  102 Bi  0.83333333333333  0.16666666666667  0.86091661265428 208.980 > 91  103 Bi  0.16666666666667  0.50000000000000  0.86091661265428 208.980 > 91  104 Bi  0.50000000000000  0.50000000000000  0.86091661265428 208.980 > 91  105 Bi  0.83333333333333  0.50000000000000  0.86091661265428 208.980 > 91  106 Bi  0.16666666666667  0.83333333333333  0.86091661265428 208.980 > 91  107 Bi  0.50000000000000  0.83333333333333  0.86091661265428 208.980 > 91  108 Bi  0.83333333333333  0.83333333333333  0.86091661265428 208.980 > 91  109 Pr  0.00000000000000  0.00000000000000  0.00000000000000 140.908 > 109  110 Pr  0.33333333333333  0.00000000000000  0.00000000000000 140.908 > 109  111 Pr  0.66666666666667  0.00000000000000  0.00000000000000 140.908 > 109  112 Pr  0.00000000000000  0.33333333333333  0.00000000000000 140.908 > 109  113 Pr  0.33333333333333  0.33333333333333  0.00000000000000 140.908 > 109  114 Pr  0.66666666666667  0.33333333333333  0.00000000000000 140.908 > 109  115 Pr  0.00000000000000  0.66666666666667  0.00000000000000 140.908 > 109  116 Pr  0.33333333333333  0.66666666666667  0.00000000000000 140.908 > 109  117 Pr  0.66666666666667  0.66666666666667  0.00000000000000 140.908 > 109  118 Pr  0.00000000000000  0.00000000000000  0.50000000000000 140.908 > 109  119 Pr  0.33333333333333  0.00000000000000  0.50000000000000 140.908 > 109  120 Pr  0.66666666666667  0.00000000000000  0.50000000000000 140.908 > 109  121 Pr  0.00000000000000  0.33333333333333  0.50000000000000 140.908 > 109  122 Pr  0.33333333333333  0.33333333333333  0.50000000000000 140.908 > 109  123 Pr  0.66666666666667  0.33333333333333  0.50000000000000 140.908 > 109  124 Pr  0.00000000000000  0.66666666666667  0.50000000000000 140.908 > 109  125 Pr  0.33333333333333  0.66666666666667  0.50000000000000 140.908 > 109  126 Pr  0.66666666666667  0.66666666666667  0.50000000000000 140.908 > 109  127 Br  0.00000000000000  0.00000000000000  0.25000000000000  79.904 > 127  128 Br  0.33333333333333  0.00000000000000  0.25000000000000  79.904 > 127  129 Br  0.66666666666667  0.00000000000000  0.25000000000000  79.904 > 127  130 Br  0.00000000000000  0.33333333333333  0.25000000000000  79.904 > 127  131 Br  0.33333333333333  0.33333333333333  0.25000000000000  79.904 > 127  132 Br  0.66666666666667  0.33333333333333  0.25000000000000  79.904 > 127  133 Br  0.00000000000000  0.66666666666667  0.25000000000000  79.904 > 127  134 Br  0.33333333333333  0.66666666666667  0.25000000000000  79.904 > 127  135 Br  0.66666666666667  0.66666666666667  0.25000000000000  79.904 > 127  136 Br  0.00000000000000  0.00000000000000  0.75000000000000  79.904 > 127  137 Br  0.33333333333333  0.00000000000000  0.75000000000000  79.904 > 127  138 Br  0.66666666666667  0.00000000000000  0.75000000000000  79.904 > 127  139 Br  0.00000000000000  0.33333333333333  0.75000000000000  79.904 > 127  140 Br  0.33333333333333  0.33333333333333  0.75000000000000  79.904 > 127  141 Br  0.66666666666667  0.33333333333333  0.75000000000000  79.904 > 127  142 Br  0.00000000000000  0.66666666666667  0.75000000000000  79.904 > 127  143 Br  0.33333333333333  0.66666666666667  0.75000000000000  79.904 > 127  144 Br  0.66666666666667  0.66666666666667  0.75000000000000  79.904 > 127----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            6.8146433    0.0000000    0.0000000            0.0000000    6.8146433    0.0000000            0.0000000    0.0000000    5.6056542-------------------------- Born effective charges --------------------------    1 O    -2.1895334    0.0000000    0.0000000            0.0000000   -4.0265966    0.0000000            0.0000000    0.0000000   -2.3026890    2 O    -2.1895334    0.0000000    0.0000000            0.0000000   -4.0265966    0.0000000            0.0000000    0.0000000   -2.3026890    3 O    -4.0265966    0.0000000    0.0000000            0.0000000   -2.1895334    0.0000000            0.0000000    0.0000000   -2.3026890    4 O    -4.0265966    0.0000000    0.0000000            0.0000000   -2.1895334    0.0000000            0.0000000    0.0000000   -2.3026890    5 Bi    5.3589907    0.0000000    0.0000000            0.0000000    5.3589907    0.0000000            0.0000000    0.0000000    2.6408645    6 Bi    5.3589907    0.0000000    0.0000000            0.0000000    5.3589907    0.0000000            0.0000000    0.0000000    2.6408645    7 Pr    3.7895244    0.0000000    0.0000000            0.0000000    3.7895244    0.0000000            0.0000000    0.0000000    5.6329384    8 Br   -2.0752457    0.0000000    0.0000000            0.0000000   -2.0752457    0.0000000            0.0000000    0.0000000   -1.7039114----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: 0.00000000 (zxx) 0.00000000 (zxx) 0.00000000 (zxx)Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- 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 ... [ 13 13 5 ]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.79, Number of G-points: 299, Lambda: 0.21Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/84) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 84Calculating 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.707   (   0.000    0.000    0.000)    0.000   1.707   (   0.000    0.000    0.000)    0.000   1.874   (   0.000    0.000    0.000)    0.000   1.874   (   0.000    0.000    0.000)    0.000   2.451   (   0.000    0.000    0.000)    0.000   3.250   (   0.000    0.000    0.000)    0.000   3.250   (   0.000    0.000    0.000)    0.000   3.457   (   0.000    0.000    0.000)    0.000   5.961   (   0.000    0.000    0.000)    0.000   8.387   (   0.000    0.000    0.000)    0.000   8.387   (   0.000    0.000    0.000)    0.000   8.778   (   0.000    0.000    0.000)    0.000   8.778   (   0.000    0.000    0.000)    0.000   9.371   (   0.000    0.000    0.000)    0.000  12.044   (   0.000    0.000    0.000)    0.000  12.110   (   0.000    0.000    0.000)    0.000  12.927   (   0.000    0.000    0.000)    0.000  12.927   (   0.000    0.000    0.000)    0.000  13.254   (   0.000    0.000    0.000)    0.000  13.254   (   0.000    0.000    0.000)    0.000  15.396   (   0.000    0.000    0.000)    0.000======================= Grid point 1 (2/84) =======================q-point: ( 0.08  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.262   (  16.357    0.000    0.000)   16.357   0.566   (  27.770    0.000    0.000)   27.770   0.860   (  42.880    0.000    0.000)   42.880   1.733   (   2.522    0.000    0.000)    2.522   1.847   (  13.095    0.000    0.000)   13.095   1.887   (   1.257    0.000    0.000)    1.257   2.426   (   1.692    0.000    0.000)    1.692   2.447   (  -0.255    0.000    0.000)    0.255   3.245   (  -0.527    0.000    0.000)    0.527   3.269   (  -3.355    0.000    0.000)    3.355   3.663   (  16.624    0.000    0.000)   16.624   5.858   (  -9.632    0.000    0.000)    9.632   8.506   (  12.013    0.000    0.000)   12.013   8.887   (  11.007    0.000    0.000)   11.007   8.892   (   8.592    0.000    0.000)    8.592   9.372   (   1.831    0.000    0.000)    1.831  11.856   (  -1.383    0.000    0.000)    1.383  12.045   (   0.646    0.000    0.000)    0.646  12.099   (  -1.210    0.000    0.000)    1.210  12.902   (  -2.588    0.000    0.000)    2.588  13.216   (  -3.859    0.000    0.000)    3.859  13.319   (   6.528    0.000    0.000)    6.528  15.191   (  -6.242    0.000    0.000)    6.242  15.358   (  -4.085    0.000    0.000)    4.085======================= Grid point 2 (3/84) =======================q-point: ( 0.15  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.630   (  20.714    0.000    0.000)   20.714   1.065   (  23.029    0.000    0.000)   23.029   1.646   (  36.741    0.000    0.000)   36.741   1.796   (   3.719    0.000    0.000)    3.719   1.921   (   2.166    0.000    0.000)    2.166   2.145   (  15.046    0.000    0.000)   15.046   2.461   (   2.846    0.000    0.000)    2.846   2.463   (   1.714    0.000    0.000)    1.714   3.230   (  -1.029    0.000    0.000)    1.029   3.243   (   1.238    0.000    0.000)    1.238   4.057   (  23.012    0.000    0.000)   23.012   5.633   ( -12.209    0.000    0.000)   12.209   8.841   (  21.768    0.000    0.000)   21.768   8.924   (  -6.976    0.000    0.000)    6.976   9.195   (  20.126    0.000    0.000)   20.126   9.560   (  17.575    0.000    0.000)   17.575  11.794   (  -4.721    0.000    0.000)    4.721  12.055   (  -3.531    0.000    0.000)    3.531  12.085   (   3.303    0.000    0.000)    3.303  12.826   (  -5.215    0.000    0.000)    5.215  13.104   (  -7.711    0.000    0.000)    7.711  13.498   (  11.341    0.000    0.000)   11.341  15.032   (  -9.557    0.000    0.000)    9.557  15.235   (  -8.276    0.000    0.000)    8.276======================= Grid point 3 (4/84) =======================q-point: ( 0.23  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.048   (  21.668    0.000    0.000)   21.668   1.458   (  17.241    0.000    0.000)   17.241   1.868   (   3.510    0.000    0.000)    3.510   1.967   (   2.504    0.000    0.000)    2.504   2.243   (  22.896    0.000    0.000)   22.896   2.321   (   3.186    0.000    0.000)    3.186   2.482   (   0.317    0.000    0.000)    0.317   2.630   (  13.427    0.000    0.000)   13.427   3.206   (  -1.413    0.000    0.000)    1.413   3.354   (  11.241    0.000    0.000)   11.241   4.528   (  24.337    0.000    0.000)   24.337   5.416   (  -9.460    0.000    0.000)    9.460   8.688   ( -15.618    0.000    0.000)   15.618   9.329   (  27.506    0.000    0.000)   27.506   9.652   (  26.105    0.000    0.000)   26.105   9.953   (  20.726    0.000    0.000)   20.726  11.690   (  -5.431    0.000    0.000)    5.431  11.961   (  -5.972    0.000    0.000)    5.972  12.157   (   3.563    0.000    0.000)    3.563  12.698   (  -7.869    0.000    0.000)    7.869  12.916   ( -11.539    0.000    0.000)   11.539  13.735   (  12.250    0.000    0.000)   12.250  14.833   ( -10.534    0.000    0.000)   10.534  15.048   ( -10.564    0.000    0.000)   10.564======================= Grid point 4 (5/84) =======================q-point: ( 0.31  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.464   (  20.626    0.000    0.000)   20.626   1.740   (  11.780    0.000    0.000)   11.780   1.928   (   2.574    0.000    0.000)    2.574   2.014   (   2.197    0.000    0.000)    2.197   2.355   (   1.270    0.000    0.000)    1.270   2.425   (  -1.616    0.000    0.000)    1.616   2.580   (   8.612    0.000    0.000)    8.612   2.869   (   9.486    0.000    0.000)    9.486   3.177   (  -1.524    0.000    0.000)    1.524   3.675   (  19.832    0.000    0.000)   19.832   4.977   (  21.091    0.000    0.000)   21.091   5.280   (  -4.262    0.000    0.000)    4.262   8.358   ( -17.303    0.000    0.000)   17.303   9.882   (  28.382    0.000    0.000)   28.382  10.191   (  28.561    0.000    0.000)   28.561  10.329   (  17.235    0.000    0.000)   17.235  11.598   (  -3.808    0.000    0.000)    3.808  11.833   (  -6.781    0.000    0.000)    6.781  12.211   (   1.837    0.000    0.000)    1.837  12.520   ( -10.286    0.000    0.000)   10.286  12.654   ( -15.239    0.000    0.000)   15.239  13.952   (   9.548    0.000    0.000)    9.548  14.634   (  -9.462    0.000    0.000)    9.462  14.840   ( -10.295    0.000    0.000)   10.295======================= Grid point 5 (6/84) =======================q-point: ( 0.38  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.831   (  16.087    0.000    0.000)   16.087   1.923   (   7.094    0.000    0.000)    7.094   1.968   (   1.489    0.000    0.000)    1.489   2.049   (   1.247    0.000    0.000)    1.247   2.375   (  -2.738    0.000    0.000)    2.738   2.385   (   1.990    0.000    0.000)    1.990   2.694   (   3.593    0.000    0.000)    3.593   2.982   (   2.258    0.000    0.000)    2.258   3.149   (  -1.214    0.000    0.000)    1.214   4.045   (  16.399    0.000    0.000)   16.399   5.241   (  -0.221    0.000    0.000)    0.221   5.327   (  14.265    0.000    0.000)   14.265   8.048   ( -13.605    0.000    0.000)   13.605  10.401   (  23.764    0.000    0.000)   23.764  10.611   (  11.367    0.000    0.000)   11.367  10.744   (  27.435    0.000    0.000)   27.435  11.546   (  -1.541    0.000    0.000)    1.541  11.711   (  -5.351    0.000    0.000)    5.351  12.227   (  -0.088    0.000    0.000)    0.088  12.306   ( -11.183    0.000    0.000)   11.183  12.325   ( -18.211    0.000    0.000)   18.211  14.099   (   5.504    0.000    0.000)    5.504  14.478   (  -6.235    0.000    0.000)    6.235  14.660   (  -7.788    0.000    0.000)    7.788======================= Grid point 6 (7/84) =======================q-point: ( 0.46  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.987   (   0.479    0.000    0.000)    0.479   2.022   (   3.193    0.000    0.000)    3.193   2.052   (   5.816    0.000    0.000)    5.816   2.058   (  -0.412    0.000    0.000)    0.412   2.334   (  -1.155    0.000    0.000)    1.155   2.431   (   2.092    0.000    0.000)    2.092   2.737   (   1.019    0.000    0.000)    1.019   2.985   (  -0.832    0.000    0.000)    0.832   3.132   (  -0.470    0.000    0.000)    0.470   4.272   (   6.130    0.000    0.000)    6.130   5.249   (   0.458    0.000    0.000)    0.458   5.519   (   5.042    0.000    0.000)    5.042   7.858   (  -5.204    0.000    0.000)    5.204  10.755   (  10.619    0.000    0.000)   10.619  10.763   (   4.024    0.000    0.000)    4.024  11.217   (  18.230    0.000    0.000)   18.230  11.531   (  -0.233    0.000    0.000)    0.233  11.636   (  -2.061    0.000    0.000)    2.061  11.976   ( -14.945    0.000    0.000)   14.945  12.122   (  -6.132    0.000    0.000)    6.132  12.219   (  -0.427    0.000    0.000)    0.427  14.169   (   1.724    0.000    0.000)    1.724  14.398   (  -2.012    0.000    0.000)    2.012  14.551   (  -3.052    0.000    0.000)    3.052======================= Grid point 14 (8/84) =======================q-point: ( 0.08  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.387   (  12.505   12.505    0.000)   17.685   0.563   (  13.930   13.930    0.000)   19.699   1.288   (  29.043   29.043    0.000)   41.073   1.871   (   7.644    7.644    0.000)   10.810   1.871   (   7.712    7.712    0.000)   10.906   1.908   (   1.971    1.971    0.000)    2.787   2.442   (  -0.377   -0.377    0.000)    0.534   2.482   (   3.954    3.954    0.000)    5.591   3.154   (  -4.815   -4.815    0.000)    6.809   3.433   (   8.980    8.980    0.000)   12.699   3.712   (   8.459    8.459    0.000)   11.963   5.779   (  -7.603   -7.603    0.000)   10.753   8.396   (   0.393    0.393    0.000)    0.555   8.667   (  -5.633   -5.633    0.000)    7.966   9.287   (  21.257   21.257    0.000)   30.062   9.339   (  -1.887   -1.887    0.000)    2.668  11.960   (  -4.594   -4.594    0.000)    6.497  11.968   (   4.889    4.889    0.000)    6.914  12.105   (   0.341    0.341    0.000)    0.483  12.801   (  -6.374   -6.374    0.000)    9.014  13.111   (  -7.293   -7.293    0.000)   10.314  13.447   (   9.606    9.606    0.000)   13.585  15.191   (  -5.435   -5.435    0.000)    7.686  15.411   (   3.235    3.235    0.000)    4.575======================= Grid point 15 (9/84) =======================q-point: ( 0.15  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.709   (  19.262    8.364    0.000)   21.000   0.983   (  24.214   -2.320    0.000)   24.325   1.821   (  21.780   13.256    0.000)   25.496   1.932   (   3.679   12.798    0.000)   13.316   1.974   (   7.071    4.234    0.000)    8.242   2.169   (  15.123    2.481    0.000)   15.325   2.450   (   2.502   -0.967    0.000)    2.683   2.554   (   3.256    8.425    0.000)    9.032   3.088   (  -1.587   -4.140    0.000)    4.433   3.529   (   3.069   18.641    0.000)   18.892   4.042   (  22.627   -0.974    0.000)   22.648   5.590   ( -10.950   -4.466    0.000)   11.825   8.595   (  17.944  -17.606    0.000)   25.139   8.643   (   1.192  -18.522    0.000)   18.560   9.253   (  -3.375    0.169    0.000)    3.379   9.746   (  23.003   14.739    0.000)   27.320  11.808   (  -8.335   -2.431    0.000)    8.683  12.090   (   5.405    2.445    0.000)    5.933  12.099   (  -1.291    4.268    0.000)    4.459  12.677   (  -6.856  -11.759    0.000)   13.611  12.952   (  -9.092  -12.138    0.000)   15.165  13.667   (  12.611   13.799    0.000)   18.694  15.028   (  -9.967   -3.355    0.000)   10.516  15.400   (  -3.998   13.933    0.000)   14.495======================= Grid point 16 (10/84) =======================q-point: ( 0.23  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.107   (  21.002    6.373    0.000)   21.948   1.422   (  20.001   -1.983    0.000)   20.099   1.987   (   3.412    3.373    0.000)    4.798   2.007   (   3.770   13.011    0.000)   13.546   2.282   (  15.068    1.445    0.000)   15.137   2.335   (   2.444    1.359    0.000)    2.796   2.603   (   1.945   10.815    0.000)   10.989   2.627   (  14.097   -0.206    0.000)   14.099   3.093   (   1.685   -3.482    0.000)    3.868   3.627   (   8.251   20.113    0.000)   21.740   4.515   (  24.683   -1.044    0.000)   24.705   5.385   (  -9.253   -3.112    0.000)    9.762   8.593   (  -8.415  -11.265    0.000)   14.062   9.034   (  25.929  -25.731    0.000)   36.529   9.304   (  10.380  -19.902    0.000)   22.446  10.170   (  20.034    9.016    0.000)   21.969  11.649   (  -7.616   -6.152    0.000)    9.791  12.046   (  -4.142    8.224    0.000)    9.208  12.173   (   2.781    2.067    0.000)    3.465  12.522   (  -8.980  -14.502    0.000)   17.057  12.753   ( -11.354  -14.151    0.000)   18.143  13.922   (  12.976   16.326    0.000)   20.854  14.829   (  -9.982   -1.249    0.000)   10.059  15.279   (  -7.901   19.802    0.000)   21.320======================= Grid point 17 (11/84) =======================q-point: ( 0.31  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.513   (  20.292    5.496    0.000)   21.023   1.752   (  13.711    1.546    0.000)   13.798   2.044   (   2.538    3.423    0.000)    4.261   2.075   (   3.192   13.734    0.000)   14.100   2.361   (   1.109    0.610    0.000)    1.265   2.446   (   2.335    0.584    0.000)    2.407   2.654   (   3.574    7.335    0.000)    8.159   2.873   (   9.696    0.601    0.000)    9.714   3.125   (   0.981   -2.492    0.000)    2.678   3.867   (  15.424   16.979    0.000)   22.939   4.971   (  21.404   -0.460    0.000)   21.409   5.249   (  -4.458   -3.191    0.000)    5.482   8.347   ( -14.909   -1.746    0.000)   15.010   9.570   (  28.290  -29.091    0.000)   40.578   9.618   (  19.971  -36.952    0.000)   42.003  10.515   (  15.041    0.253    0.000)   15.043  11.522   (  -5.234   -9.574    0.000)   10.911  11.944   (  -6.111    9.227    0.000)   11.068  12.162   (  -6.183   -6.459    0.000)    8.941  12.354   (  -6.137  -11.036    0.000)   12.628  12.514   ( -12.722  -11.417    0.000)   17.094  14.152   (  10.168   17.924    0.000)   20.607  14.653   (  -7.816    1.872    0.000)    8.037  15.113   (  -8.628   23.266    0.000)   24.815======================= Grid point 18 (12/84) =======================q-point: ( 0.38  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.875   (  15.930    5.429    0.000)   16.830   1.961   (   7.921    3.937    0.000)    8.845   2.084   (   1.537    3.715    0.000)    4.020   2.129   (   2.226   14.677    0.000)   14.845   2.389   (   1.965    0.447    0.000)    2.016   2.440   (  -1.600    5.217    0.000)    5.457   2.721   (   2.791    3.225    0.000)    4.266   2.987   (   2.156    0.707    0.000)    2.269   3.130   (  -0.193   -1.297    0.000)    1.311   4.168   (  13.888   11.950    0.000)   18.321   5.207   (  -0.342   -3.563    0.000)    3.580   5.326   (  14.475    0.002    0.000)   14.475   8.068   ( -12.592    1.803    0.000)   12.721  10.029   (  21.037  -41.430    0.000)   46.465  10.110   (  26.201  -28.156    0.000)   38.461  10.737   (   7.109  -13.597    0.000)   15.343  11.448   (  -2.270  -13.954    0.000)   14.137  11.818   (  -6.399    4.627    0.000)    7.897  11.921   ( -16.194  -19.909    0.000)   25.663  12.279   ( -10.655   -0.340    0.000)   10.660  12.301   (  -1.024   -3.683    0.000)    3.823  14.310   (   5.982   19.038    0.000)   19.956  14.531   (  -4.567    5.340    0.000)    7.026  14.960   (  -6.629   25.507    0.000)   26.354======================= Grid point 19 (13/84) =======================q-point: ( 0.46  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.068   (   3.235    4.591    0.000)    5.616   2.098   (   6.513    7.307    0.000)    9.789   2.099   (  -0.139    4.217    0.000)    4.219   2.154   (  -0.084   13.492    0.000)   13.492   2.412   (  -0.874    6.902    0.000)    6.957   2.437   (   2.310    0.738    0.000)    2.425   2.758   (   0.972    2.546    0.000)    2.725   2.987   (  -1.010    0.322    0.000)    1.060   3.125   (  -0.218   -0.680    0.000)    0.715   4.364   (   5.384    9.249    0.000)   10.702   5.213   (   0.447   -3.688    0.000)    3.715   5.520   (   5.116    0.253    0.000)    5.123   7.891   (  -4.889    3.161    0.000)    5.821  10.386   (  13.145  -33.431    0.000)   35.923  10.535   (  14.416  -21.309    0.000)   25.727  10.779   (  -2.101  -36.230    0.000)   36.291  11.438   (   1.724  -27.284    0.000)   27.338  11.606   ( -13.556  -19.660    0.000)   23.881  11.718   (  -3.098   -1.171    0.000)    3.312  12.130   (  -4.056   14.632    0.000)   15.183  12.288   (  -0.399   -2.979    0.000)    3.005  14.386   (   1.909   19.660    0.000)   19.753  14.474   (  -1.403    7.514    0.000)    7.644  14.868   (  -2.513   26.833    0.000)   26.950======================= Grid point 28 (14/84) =======================q-point: ( 0.15  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.952   (  15.891   15.891    0.000)   22.473   1.101   (  14.044   14.044    0.000)   19.861   2.037   (   4.472    4.472    0.000)    6.324   2.086   (  10.972   10.972    0.000)   15.516   2.200   (   8.392    8.392    0.000)   11.867   2.273   (  12.086   12.086    0.000)   17.092   2.433   (   0.284    0.284    0.000)    0.402   2.708   (   5.238    5.238    0.000)    7.407   3.080   (   3.134    3.134    0.000)    4.432   3.904   (  14.461   14.461    0.000)   20.451   4.091   (  11.472   11.472    0.000)   16.223   5.464   (  -8.079   -8.079    0.000)   11.425   8.351   ( -10.274  -10.274    0.000)   14.530   8.390   (  -1.383   -1.383    0.000)    1.955   9.192   (  -6.120   -6.120    0.000)    8.655  10.018   (  12.430   12.430    0.000)   17.579  11.659   ( -11.015  -11.015    0.000)   15.578  12.175   (   4.321    4.321    0.000)    6.111  12.198   (   4.960    4.960    0.000)    7.014  12.450   ( -11.192  -11.192    0.000)   15.829  12.690   ( -14.237  -14.237    0.000)   20.134  13.969   (  16.459   16.459    0.000)   23.277  14.921   (  -6.616   -6.616    0.000)    9.356  15.611   (   6.127    6.127    0.000)    8.665======================= Grid point 29 (15/84) =======================q-point: ( 0.23  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.300   (  19.235   13.191    0.000)   23.323   1.455   (  20.080    6.454    0.000)   21.092   2.098   (   2.644    7.159    0.000)    7.631   2.293   (   9.283    0.375    0.000)    9.290   2.305   (   4.844   13.601    0.000)   14.438   2.379   (  -0.256    5.174    0.000)    5.181   2.637   (  16.193    2.030    0.000)   16.320   2.758   (   1.098    3.257    0.000)    3.437   3.153   (   2.657    8.726    0.000)    9.122   4.067   (   8.054   24.119    0.000)   25.429   4.508   (  24.827    1.208    0.000)   24.856   5.295   (  -8.537   -6.044    0.000)   10.460   8.246   (  -1.346  -21.074    0.000)   21.117   8.526   (  14.107  -22.473    0.000)   26.534   9.036   (  -8.607   -9.474    0.000)   12.800  10.230   (   9.158   -1.476    0.000)    9.276  11.429   ( -11.947  -15.305    0.000)   19.416  12.200   (  -6.637   -4.224    0.000)    7.868  12.228   (  -5.640    6.996    0.000)    8.987  12.249   (  -5.070   -9.034    0.000)   10.360  12.442   (  -7.604  -14.916    0.000)   16.742  14.286   (  15.448   19.629    0.000)   24.979  14.797   (  -5.619   -1.456    0.000)    5.804  15.646   (  -1.757   15.351    0.000)   15.451======================= Grid point 30 (16/84) =======================q-point: ( 0.31  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.684   (  19.633   11.755    0.000)   22.883   1.823   (  16.572    6.233    0.000)   17.705   2.146   (   2.330    6.618    0.000)    7.017   2.376   (   0.733    0.949    0.000)    1.199   2.400   (   4.347   17.801    0.000)   18.324   2.439   (   5.553   -0.898    0.000)    5.626   2.783   (   1.814    4.068    0.000)    4.454   2.906   (  10.349    3.385    0.000)   10.888   3.165   (  -0.996    7.069    0.000)    7.138   4.262   (  11.302   22.359    0.000)   25.053   4.972   (  21.843    1.006    0.000)   21.866   5.158   (  -5.047   -6.098    0.000)    7.916   8.224   (  -3.364  -13.712    0.000)   14.118   8.896   (  22.944  -36.369    0.000)   43.001   8.944   (   2.048  -27.968    0.000)   28.043  10.359   (   3.317  -13.450    0.000)   13.853  11.203   ( -11.375  -21.561    0.000)   24.378  11.951   ( -15.222  -12.512    0.000)   19.704  11.968   ( -16.430  -12.433    0.000)   20.604  12.210   (  -2.280   -6.671    0.000)    7.050  12.413   (   0.890    5.582    0.000)    5.653  14.554   (  11.683   21.504    0.000)   24.473  14.713   (  -2.936    4.435    0.000)    5.319  15.573   (  -5.040   21.200    0.000)   21.791======================= Grid point 31 (17/84) =======================q-point: ( 0.38  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.043   (  16.194   11.299    0.000)   19.746   2.079   (   9.461    8.106    0.000)   12.458   2.187   (   1.850    6.460    0.000)    6.719   2.402   (   2.062    0.844    0.000)    2.228   2.477   (   3.270   19.323    0.000)   19.598   2.512   (   2.169    1.622    0.000)    2.708   2.828   (   2.503    6.613    0.000)    7.070   3.024   (   1.843    3.759    0.000)    4.186   3.136   (  -1.602    2.929    0.000)    3.338   4.477   (   9.880   18.782    0.000)   21.222   5.103   (  -0.800   -6.965    0.000)    7.011   5.335   (  14.797    1.134    0.000)   14.840   8.093   (  -8.390   -0.378    0.000)    8.399   9.151   (  17.769  -46.159    0.000)   49.461   9.398   (  27.993  -41.911    0.000)   50.400  10.318   (  -8.460  -25.659    0.000)   27.018  10.956   ( -15.256  -31.460    0.000)   34.964  11.645   ( -14.657  -11.049    0.000)   18.355  11.679   ( -11.823  -16.989    0.000)   20.698  12.163   (  -2.182  -10.696    0.000)   10.916  12.430   (   0.671   12.396    0.000)   12.414  14.679   (  -0.822    9.439    0.000)    9.475  14.735   (   6.867   22.498    0.000)   23.522  15.474   (  -4.595   24.472    0.000)   24.900======================= Grid point 32 (18/84) =======================q-point: ( 0.46  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.190   (   2.068    7.143    0.000)    7.436   2.221   (   1.558    8.439    0.000)    8.582   2.264   (   5.363    9.448    0.000)   10.864   2.464   (   3.812    3.303    0.000)    5.044   2.518   (   0.661   20.563    0.000)   20.574   2.534   (   0.392    4.069    0.000)    4.088   2.868   (   1.234    9.158    0.000)    9.241   3.010   (  -1.721    2.825    0.000)    3.308   3.112   (  -0.673   -0.295    0.000)    0.734   4.618   (   3.934   16.205    0.000)   16.676   5.105   (   0.375   -7.348    0.000)    7.358   5.533   (   5.231    1.240    0.000)    5.376   7.962   (  -3.886    3.323    0.000)    5.112   9.570   (  23.563  -48.515    0.000)   53.935   9.958   (  27.604  -37.586    0.000)   46.633  10.026   ( -20.351  -37.201    0.000)   42.403  10.571   ( -23.140  -45.006    0.000)   50.607  11.443   (  -5.192   -9.151    0.000)   10.521  11.532   (  -3.416  -16.656    0.000)   17.002  12.133   (  -0.821  -12.392    0.000)   12.419  12.439   (   0.237   14.993    0.000)   14.995  14.672   (  -0.066   12.172    0.000)   12.172  14.823   (   2.212   22.941    0.000)   23.047  15.409   (  -1.797   25.951    0.000)   26.013======================= Grid point 42 (19/84) =======================q-point: ( 0.23  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.603   (  17.269   17.269    0.000)   24.422   1.675   (  15.288   15.288    0.000)   21.620   2.227   (   4.885    4.885    0.000)    6.908   2.323   (   3.513    3.513    0.000)    4.968   2.381   (   0.927    0.927    0.000)    1.311   2.584   (   7.327    7.327    0.000)   10.361   2.779   (  13.518   13.518    0.000)   19.117   2.788   (   1.328    1.328    0.000)    1.878   3.307   (   4.854    4.854    0.000)    6.864   4.508   (  15.839   15.839    0.000)   22.399   4.619   (  14.862   14.862    0.000)   21.019   5.157   (  -7.698   -7.698    0.000)   10.886   7.903   ( -11.999  -11.999    0.000)   16.970   8.235   (  -6.767   -6.767    0.000)    9.570   8.837   ( -12.106  -12.106    0.000)   17.121  10.172   (  -3.476   -3.476    0.000)    4.916  11.099   ( -17.678  -17.678    0.000)   25.000  11.972   ( -12.622  -12.622    0.000)   17.851  12.000   ( -21.073  -21.073    0.000)   29.801  12.171   (  -6.024   -6.024    0.000)    8.519  12.463   (   7.550    7.550    0.000)   10.677  14.653   (  17.382   17.382    0.000)   24.582  14.788   (   0.733    0.733    0.000)    1.037  15.849   (   5.779    5.779    0.000)    8.173======================= Grid point 43 (20/84) =======================q-point: ( 0.31  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.957   (  18.553   15.749    0.000)   24.336   2.002   (  16.687   11.960    0.000)   20.531   2.283   (   1.922    7.029    0.000)    7.287   2.395   (   0.773    0.906    0.000)    1.190   2.425   (   5.711   -0.297    0.000)    5.719   2.700   (   5.596   11.482    0.000)   12.773   2.834   (   3.524    2.281    0.000)    4.197   3.045   (  11.399   11.540    0.000)   16.221   3.328   (  -2.074    7.452    0.000)    7.736   4.707   (   9.004   22.500    0.000)   24.234   5.017   (  -6.055   -8.107    0.000)   10.119   5.026   (  20.989    4.733    0.000)   21.516   7.814   (   2.829  -23.427    0.000)   23.598   8.265   (   9.624  -26.136    0.000)   27.851   8.567   ( -14.510  -14.570    0.000)   20.563  10.085   (  -6.704  -12.974    0.000)   14.604  10.748   ( -18.543  -23.934    0.000)   30.277  11.576   ( -21.092  -24.769    0.000)   32.532  11.714   ( -14.486  -11.731    0.000)   18.641  12.013   (  -8.781  -14.075    0.000)   16.590  12.584   (   4.797    8.816    0.000)   10.036  14.825   (   2.768    6.795    0.000)    7.337  14.948   (  12.572   18.211    0.000)   22.129  15.896   (  -0.244   11.829    0.000)   11.832======================= Grid point 44 (21/84) =======================q-point: ( 0.38  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.268   (   8.535   10.538    0.000)   13.561   2.297   (  14.314   13.919    0.000)   19.965   2.327   (   3.823    8.170    0.000)    9.021   2.424   (   3.527    1.877    0.000)    3.995   2.521   (   3.871   -0.182    0.000)    3.875   2.803   (   4.622   12.267    0.000)   13.109   2.926   (   5.798    3.517    0.000)    6.782   3.183   (   2.573   13.707    0.000)   13.946   3.256   (  -4.678    8.296    0.000)    9.524   4.868   (   7.133   20.575    0.000)   21.776   4.939   (  -1.852   -9.704    0.000)    9.879   5.379   (  14.506    3.441    0.000)   14.909   7.959   (   7.788  -17.761    0.000)   19.394   8.357   (  -2.345  -29.984    0.000)   30.076   8.594   (  23.493  -37.551    0.000)   44.295   9.849   ( -18.558  -20.273    0.000)   27.484  10.348   ( -23.391  -29.702    0.000)   37.807  11.237   ( -12.464  -27.594    0.000)   30.278  11.414   ( -14.758  -12.829    0.000)   19.555  11.873   (  -5.178  -19.140    0.000)   19.828  12.654   (   2.471    9.723    0.000)   10.032  14.883   (   2.854   11.042    0.000)   11.404  15.140   (   7.224   18.416    0.000)   19.782  15.868   (  -2.061   15.624    0.000)   15.760======================= Grid point 45 (22/84) =======================q-point: ( 0.46  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.325   (   0.743    6.637    0.000)    6.679   2.402   (   0.882    3.586    0.000)    3.693   2.429   (   3.216   12.431    0.000)   12.841   2.571   (   1.293    0.526    0.000)    1.396   2.600   (   8.255   11.936    0.000)   14.513   2.866   (   1.624   11.536    0.000)   11.649   3.046   (   6.088    4.932    0.000)    7.835   3.158   (  -5.300    7.968    0.000)    9.570   3.171   (  -2.107   16.120    0.000)   16.257   4.929   (   0.146  -10.604    0.000)   10.605   4.969   (   2.769   18.984    0.000)   19.185   5.574   (   5.150    2.939    0.000)    5.930   7.968   (  -2.654   -4.518    0.000)    5.240   8.593   (  18.939  -48.764    0.000)   52.313   9.138   (  30.014  -42.335    0.000)   51.896   9.407   ( -22.092  -24.503    0.000)   32.991   9.822   ( -29.273  -32.319    0.000)   43.605  11.093   (  -3.056  -27.938    0.000)   28.105  11.203   (  -5.825  -15.245    0.000)   16.320  11.809   (  -1.562  -20.771    0.000)   20.829  12.685   (   0.767   10.111    0.000)   10.140  14.924   (   1.166   13.126    0.000)   13.178  15.233   (   2.308   18.426    0.000)   18.570  15.834   (  -1.015   17.375    0.000)   17.405======================= Grid point 56 (23/84) =======================q-point: ( 0.31  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.270   (  14.508   14.508    0.000)   20.517   2.283   (  17.138   17.138    0.000)   24.237   2.403   (   4.281    4.281    0.000)    6.054   2.413   (   1.167    1.167    0.000)    1.651   2.439   (   2.734    2.734    0.000)    3.867   2.867   (   6.884    6.884    0.000)    9.736   2.894   (   4.013    4.013    0.000)    5.676   3.309   (  13.363   13.363    0.000)   18.898   3.417   (   1.599    1.599    0.000)    2.262   4.856   (  -7.850   -7.850    0.000)   11.102   5.096   (  13.762   13.762    0.000)   19.463   5.170   (  12.590   12.590    0.000)   17.804   7.475   (  -9.158   -9.158    0.000)   12.952   7.938   (  -6.333   -6.333    0.000)    8.956   8.266   ( -16.633  -16.633    0.000)   23.523   9.865   ( -10.895  -10.895    0.000)   15.408  10.286   ( -23.456  -23.456    0.000)   33.172  11.059   ( -26.677  -26.677    0.000)   37.728  11.485   ( -12.745  -12.745    0.000)   18.025  11.679   ( -18.278  -18.278    0.000)   25.848  12.724   (   5.474    5.474    0.000)    7.742  14.965   (   7.228    7.228    0.000)   10.222  15.250   (  12.673   12.673    0.000)   17.923  16.048   (   4.335    4.335    0.000)    6.131======================= Grid point 57 (24/84) =======================q-point: ( 0.38  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.427   (   1.548    6.055    0.000)    6.250   2.429   (   1.003    1.435    0.000)    1.751   2.519   (   3.973    0.238    0.000)    3.980   2.550   (  10.793   13.611    0.000)   17.371   2.618   (  15.950   15.734    0.000)   22.404   2.982   (   4.589    7.708    0.000)    8.970   2.989   (   5.519    3.243    0.000)    6.402   3.407   (  -2.024    6.355    0.000)    6.669   3.508   (   6.268   17.239    0.000)   18.343   4.735   (  -3.999  -10.806    0.000)   11.522   5.252   (   5.607   18.087    0.000)   18.936   5.465   (  13.273    5.103    0.000)   14.220   7.483   (   9.515  -22.235    0.000)   24.186   7.951   ( -14.357  -18.097    0.000)   23.101   8.029   (  15.018  -18.565    0.000)   23.879   9.570   ( -21.147   -8.800    0.000)   22.905   9.814   ( -25.684  -24.942    0.000)   35.802  10.627   ( -15.349  -33.966    0.000)   37.273  11.152   ( -19.267  -13.656    0.000)   23.616  11.424   (  -7.773  -26.268    0.000)   27.394  12.803   (   2.760    5.672    0.000)    6.308  15.092   (   5.378    9.931    0.000)   11.294  15.442   (   7.104   12.461    0.000)   14.344  16.090   (   0.615    7.653    0.000)    7.678======================= Grid point 58 (25/84) =======================q-point: ( 0.46  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.443   (   0.393    1.308    0.000)    1.365   2.445   (   0.491    5.568    0.000)    5.590   2.572   (   1.383    0.259    0.000)    1.407   2.694   (   3.662   13.880    0.000)   14.355   2.878   (   9.326   15.255    0.000)   17.880   3.027   (  -0.297    7.647    0.000)    7.652   3.086   (   3.446    1.314    0.000)    3.688   3.363   (  -1.784   10.325    0.000)   10.478   3.559   (   0.128   20.469    0.000)   20.470   4.698   (  -0.423  -12.797    0.000)   12.804   5.328   (   2.056   17.164    0.000)   17.287   5.645   (   4.766    4.126    0.000)    6.304   7.741   (  -0.115  -19.105    0.000)   19.106   7.766   (   8.115  -32.850    0.000)   33.837   8.450   (  23.693  -25.721    0.000)   34.971   9.116   ( -19.093   -6.080    0.000)   20.038   9.285   ( -25.643  -22.642    0.000)   34.208  10.462   (  -3.159  -36.006    0.000)   36.144  10.866   (  -8.126  -18.694    0.000)   20.384  11.328   (  -2.392  -28.101    0.000)   28.202  12.837   (   0.831    5.688    0.000)    5.748  15.165   (   1.970   11.156    0.000)   11.329  15.532   (   2.242   12.284    0.000)   12.487  16.091   (  -0.155    9.335    0.000)    9.336======================= Grid point 70 (26/84) =======================q-point: ( 0.38  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.450   (   0.871    0.871    0.000)    1.232   2.518   (   2.564    2.564    0.000)    3.626   2.535   (   2.037    2.037    0.000)    2.881   2.802   (  11.123   11.123    0.000)   15.730   2.919   (  14.261   14.261    0.000)   20.169   3.055   (   3.426    3.426    0.000)    4.846   3.124   (   6.650    6.650    0.000)    9.404   3.496   (   2.960    2.960    0.000)    4.186   3.795   (  11.050   11.050    0.000)   15.628   4.541   (  -8.280   -8.280    0.000)   11.709   5.543   (   8.773    8.773    0.000)   12.407   5.570   (   7.763    7.763    0.000)   10.978   7.235   (  -2.856   -2.856    0.000)    4.039   7.605   ( -16.049  -16.049    0.000)   22.697   7.906   (   5.808    5.808    0.000)    8.213   9.363   ( -22.011  -22.011    0.000)   31.129   9.418   ( -10.872  -10.872    0.000)   15.376   9.996   ( -26.244  -26.244    0.000)   37.115  10.875   ( -20.771  -20.771    0.000)   29.375  10.924   ( -16.351  -16.351    0.000)   23.124  12.883   (   2.733    2.733    0.000)    3.864  15.257   (   6.764    6.764    0.000)    9.565  15.629   (   6.858    6.858    0.000)    9.698  16.187   (   2.781    2.781    0.000)    3.933======================= Grid point 71 (27/84) =======================q-point: ( 0.46  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.462   (   0.329    0.734    0.000)    0.804   2.538   (   0.397    3.772    0.000)    3.793   2.580   (   1.338    0.532    0.000)    1.440   2.949   (   3.617   11.286    0.000)   11.852   3.113   (   2.017    1.341    0.000)    2.422   3.127   (   5.641    8.539    0.000)   10.234   3.230   (   3.725   13.793    0.000)   14.287   3.531   (   0.783    6.719    0.000)    6.765   3.913   (   1.660   15.077    0.000)   15.168   4.442   (  -1.834  -12.924    0.000)   13.054   5.614   (   1.614   11.685    0.000)   11.796   5.722   (   4.312    3.487    0.000)    5.546   7.303   (   5.778  -12.997    0.000)   14.223   7.377   (  -6.360  -17.011    0.000)   18.161   8.178   (  17.113   -2.472    0.000)   17.290   8.940   ( -19.331  -13.022    0.000)   23.308   9.089   ( -15.485    1.313    0.000)   15.540   9.737   (  -3.568  -36.335    0.000)   36.510  10.505   ( -11.023  -17.178    0.000)   20.410  10.749   (  -4.599  -29.607    0.000)   29.962  12.916   (   0.790    2.646    0.000)    2.762  15.348   (   2.387    7.342    0.000)    7.720  15.716   (   2.145    6.698    0.000)    7.033  16.216   (   0.587    3.969    0.000)    4.013======================= Grid point 84 (28/84) =======================q-point: ( 0.46  0.46  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.472   (   0.248    0.248    0.000)    0.351   2.588   (   0.933    0.933    0.000)    1.319   2.590   (   0.710    0.710    0.000)    1.004   3.104   (   3.966    3.966    0.000)    5.608   3.135   (   0.842    0.842    0.000)    1.191   3.237   (   2.903    2.903    0.000)    4.106   3.472   (   8.501    8.501    0.000)   12.022   3.620   (   2.345    2.345    0.000)    3.316   4.116   (   4.611    4.611    0.000)    6.521   4.233   (  -6.390   -6.390    0.000)    9.037   5.769   (   2.862    2.862    0.000)    4.047   5.772   (   2.578    2.578    0.000)    3.646   7.110   (  -6.741   -6.741    0.000)    9.534   7.240   (   0.876    0.876    0.000)    1.239   8.273   (   7.926    7.926    0.000)   11.209   8.751   (  -6.612   -6.612    0.000)    9.351   9.099   (  -4.592   -4.592    0.000)    6.495   9.181   ( -12.706  -12.706    0.000)   17.969  10.241   (  -9.481   -9.481    0.000)   13.408  10.287   ( -13.020  -13.020    0.000)   18.413  12.948   (   0.741    0.741    0.000)    1.048  15.445   (   2.549    2.549    0.000)    3.604  15.800   (   2.086    2.086    0.000)    2.951  16.262   (   1.012    1.012    0.000)    1.432======================= Grid point 169 (29/84) =======================q-point: ( 0.00  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 84Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.344   (  -0.000    0.000   14.392)   14.392   0.344   (   0.000    0.000   14.392)   14.392   0.657   (   0.000    0.000   28.057)   28.057   1.644   (  -0.000   -0.000   -4.775)    4.775   1.644   (   0.000    0.000   -4.775)    4.775   1.897   (   0.000    0.000    1.441)    1.441   1.897   (   0.000    0.000    1.441)    1.441   2.987   (   0.000   -0.000  -15.964)   15.964   3.252   (   0.000    0.000    0.156)    0.156   3.252   (  -0.000    0.000    0.156)    0.156   3.623   (   0.000    0.000    8.259)    8.259   6.017   (   0.000    0.000    0.107)    0.107   8.399   (  -0.000   -0.000    1.020)    1.020   8.399   (   0.000    0.000    1.020)    1.020   8.763   (  -0.000    0.000   -1.194)    1.194   8.763   (   0.000    0.000   -1.194)    1.194   9.377   (   0.000    0.000    0.445)    0.445  12.041   (   0.000    0.000   -0.281)    0.281  12.932   (   0.000    0.000    0.408)    0.408  12.932   (   0.000    0.000    0.408)    0.408  13.251   (   0.000    0.000   -0.252)    0.252  13.251   (   0.000    0.000   -0.252)    0.252  15.287   (   0.000    0.000    2.217)    2.217  15.398   (   0.000    0.000   -0.005)    0.005======================= Grid point 170 (30/84) =======================q-point: ( 0.08  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.542   (   5.170    0.000   17.077)   17.843   0.657   (  23.404    0.000    7.111)   24.461   0.955   (  33.257    0.000    6.961)   33.978   1.675   (   2.986    0.000   -4.411)    5.327   1.815   (  15.071    0.000   -3.009)   15.369   1.908   (   1.112    0.000    1.352)    1.751   2.124   (  13.682    0.000   -6.127)   14.991   2.657   ( -13.376    0.000   -0.670)   13.393   3.247   (  -0.566    0.000    0.124)    0.580   3.347   (  -3.791    0.000    4.272)    5.711   3.626   (  14.568    0.000   -0.998)   14.602   5.887   ( -10.199    0.000    1.843)   10.364   8.518   (  11.968    0.000    0.982)   12.008   8.873   (  11.092    0.000   -1.125)   11.149   8.899   (  10.683    0.000    0.532)   10.696   9.376   (   1.552    0.000    0.315)    1.583  10.203   (  88.908    0.000  -58.602)  106.484  12.002   (  -3.689    0.000   -0.563)    3.732  12.906   (  -2.767    0.000    0.637)    2.839  12.906   (  -7.468    0.000   16.303)   17.932  13.214   (  -3.787    0.000   -0.193)    3.792  13.282   (   7.527    0.000   -2.732)    8.008  15.053   ( -14.983    0.000   -5.073)   15.818  15.324   (  -6.301    0.000   -0.780)    6.349======================= Grid point 171 (31/84) =======================q-point: ( 0.15  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.734   (  15.160    0.000    7.998)   17.140   1.108   (  21.614    0.000    3.541)   21.902   1.628   (  32.484    0.000   -1.940)   32.541   1.751   (   4.495    0.000   -3.472)    5.680   1.938   (   1.878    0.000    1.097)    2.175   2.141   (  15.682    0.000   -0.967)   15.712   2.321   (   6.675    0.000   -4.153)    7.862   2.558   (   2.800    0.000    1.965)    3.421   3.231   (  -1.088    0.000    0.045)    1.089   3.281   (  -1.182    0.000    2.779)    3.019   4.019   (  23.393    0.000   -3.019)   23.587   5.659   ( -12.306    0.000    1.906)   12.452   8.851   (  21.683    0.000    0.872)   21.700   8.966   (  -6.104    0.000    3.309)    6.943   9.184   (  20.287    0.000   -0.926)   20.309   9.566   (  18.198    0.000    0.408)   18.202  11.262   (  29.865    0.000  -28.218)   41.087  11.918   (  -4.041    0.000   -1.495)    4.309  12.640   ( -16.769    0.000   21.413)   27.197  12.828   (  -5.374    0.000    0.139)    5.375  13.103   (  -7.593    0.000   -0.040)    7.593  13.526   (  14.634    0.000    2.568)   14.858  14.858   (  -5.866    0.000  -15.887)   16.935  15.156   ( -10.571    0.000   -2.490)   10.860======================= Grid point 172 (32/84) =======================q-point: ( 0.23  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.091   (  20.359    0.000    3.464)   20.652   1.481   (  16.506    0.000    1.944)   16.620   1.840   (   4.398    0.000   -2.255)    4.943   1.978   (   2.118    0.000    0.707)    2.233   2.145   (  18.588    0.000   -6.592)   19.722   2.329   (   5.286    0.000    0.020)    5.286   2.423   (   3.560    0.000   -3.243)    4.816   2.748   (  14.209    0.000    7.364)   16.004   3.205   (  -1.466    0.000   -0.045)    1.467   3.346   (   8.927    0.000   -0.673)    8.952   4.497   (  24.810    0.000   -2.434)   24.929   5.438   (  -9.668    0.000    1.718)    9.820   8.732   ( -16.240    0.000    3.512)   16.616   9.338   (  27.397    0.000    0.707)   27.406   9.644   (  26.329    0.000   -0.615)   26.336   9.966   (  20.738    0.000    0.993)   20.761  11.582   (   4.916    0.000   -7.507)    8.974  11.892   (   2.025    0.000   -3.009)    3.627  12.326   ( -13.899    0.000   10.028)   17.139  12.696   (  -8.067    0.000   -0.140)    8.068  12.918   ( -11.420    0.000    0.149)   11.421  13.801   (  12.998    0.000    6.086)   14.352  14.766   (  -4.797    0.000   -9.445)   10.593  14.931   ( -11.705    0.000   -5.774)   13.052======================= Grid point 173 (33/84) =======================q-point: ( 0.31  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.500   (  20.974    0.000    2.845)   21.166   1.752   (  11.313    0.000    0.994)   11.356   1.916   (   3.344    0.000   -1.020)    3.496   2.018   (   1.855    0.000    0.245)    1.871   2.314   (   2.483    0.000   -3.258)    4.096   2.428   (   0.091    0.000   -0.083)    0.123   2.488   (   6.457    0.000   -5.135)    8.249   2.993   (   9.322    0.000    8.145)   12.379   3.175   (  -1.554    0.000   -0.111)    1.558   3.624   (  17.756    0.000   -4.205)   18.247   4.956   (  21.552    0.000   -1.651)   21.615   5.302   (  -3.927    0.000    1.699)    4.279   8.374   ( -19.511    0.000    1.263)   19.552   9.889   (  28.283    0.000    0.530)   28.288  10.188   (  28.834    0.000   -0.218)   28.835  10.337   (  16.914    0.000    0.612)   16.925  11.583   (  -1.788    0.000   -0.957)    2.028  11.874   (  -6.611    0.000    3.460)    7.461  12.179   (  -0.480    0.000   -2.915)    2.954  12.514   ( -10.499    0.000   -0.446)   10.508  12.658   ( -15.167    0.000    0.290)   15.170  14.019   (   9.162    0.000    6.070)   10.991  14.634   (  -8.780    0.000   -0.200)    8.782  14.741   (  -7.235    0.000   -8.416)   11.098======================= Grid point 174 (34/84) =======================q-point: ( 0.38  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.883   (  17.400    0.000    4.084)   17.873   1.927   (   6.755    0.000    0.342)    6.764   1.968   (   1.989    0.000   -0.008)    1.989   2.047   (   1.089    0.000   -0.135)    1.097   2.338   (   0.114    0.000   -1.841)    1.845   2.408   (  -0.549    0.000    0.878)    1.035   2.584   (   3.069    0.000   -6.889)    7.542   3.097   (   1.633    0.000    7.498)    7.674   3.147   (  -1.220    0.000   -0.138)    1.228   3.955   (  14.509    0.000   -7.342)   16.261   5.274   (   2.093    0.000    2.369)    3.161   5.319   (  13.484    0.000   -0.532)   13.495   8.014   ( -16.244    0.000   -2.712)   16.469  10.407   (  23.726    0.000    0.410)   23.730  10.615   (  11.324    0.000    0.289)   11.328  10.747   (  27.763    0.000    0.250)   27.764  11.559   (  -0.478    0.000    1.023)    1.129  11.716   (  -7.606    0.000    0.527)    7.624  12.204   (   1.277    0.000   -1.909)    2.297  12.295   ( -11.825    0.000   -0.735)   11.848  12.331   ( -17.837    0.000    0.274)   17.839  14.152   (   4.462    0.000    4.477)    6.321  14.475   (  -6.608    0.000   -0.248)    6.612  14.627   (  -4.709    0.000   -2.966)    5.565======================= Grid point 175 (35/84) =======================q-point: ( 0.46  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.992   (   0.654    0.000    0.444)    0.791   2.023   (   3.046    0.000   -0.002)    3.046   2.055   (  -0.437    0.000   -0.237)    0.498   2.132   (   6.981    0.000    6.413)    9.479   2.321   (  -1.120    0.000   -0.891)    1.432   2.429   (   1.932    0.000   -0.075)    1.934   2.615   (   0.378    0.000   -7.645)    7.654   3.090   (  -1.121    0.000    6.318)    6.417   3.130   (  -0.468    0.000   -0.140)    0.489   4.153   (   5.292    0.000   -9.325)   10.722   5.306   (   1.204    0.000    4.167)    4.337   5.512   (   5.159    0.000   -0.501)    5.183   7.784   (  -6.405    0.000   -5.876)    8.692  10.760   (  10.651    0.000    0.417)   10.659  10.769   (   4.149    0.000    0.470)    4.175  11.227   (  18.717    0.000    0.876)   18.737  11.559   (   0.252    0.000    2.231)    2.245  11.608   (  -2.987    0.000   -2.132)    3.670  11.974   ( -15.306    0.000   -0.146)   15.307  12.110   (  -6.282    0.000   -0.924)    6.349  12.215   (   0.061    0.000   -0.346)    0.352  14.204   (   1.168    0.000    2.700)    2.941  14.385   (  -2.497    0.000   -0.954)    2.673  14.565   (  -1.516    0.000    1.136)    1.895======================= Grid point 183 (36/84) =======================q-point: ( 0.08  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.594   (   4.707    4.707   14.064)   15.559   0.654   (  11.780   11.780    7.153)   18.130   1.320   (  26.223   26.223    2.013)   37.139   1.787   (   6.406    6.406   -5.191)   10.442   1.868   (   8.988    8.988   -0.527)   12.722   1.951   (   3.440    3.440    1.795)    5.185   2.230   (   8.400    8.400   -6.398)   13.493   2.579   (  -4.258   -4.258   -1.142)    6.130   3.218   (  -5.829   -5.829    3.807)    9.079   3.436   (   9.016    9.016    0.241)   12.753   3.682   (   8.285    8.285   -1.671)   11.835   5.807   (  -7.696   -7.696    1.926)   11.053   8.411   (   0.475    0.475    1.250)    1.419   8.665   (  -5.080   -5.080   -0.220)    7.188   9.296   (  22.181   22.181    0.767)   31.378   9.344   (  -1.943   -1.943    0.378)    2.773  10.875   (  42.885   42.885  -44.239)   75.068  11.950   (  -4.857   -4.857   -0.783)    6.913  12.804   (  -6.437   -6.437    0.260)    9.107  12.864   (  -7.784   -7.784   23.562)   26.007  13.111   (  -7.128   -7.128    0.037)   10.081  13.442   (  12.758   12.758   -0.557)   18.052  15.027   (  -2.376   -2.376  -14.194)   14.587  15.279   (  -3.964   -3.964   -2.823)    6.276======================= Grid point 184 (37/84) =======================q-point: ( 0.15  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.788   (  15.191    6.374    6.205)   17.604   1.032   (  22.360   -1.987    3.956)   22.794   1.814   (  20.874   13.973   -0.647)   25.127   1.873   (   5.239   10.629   -3.783)   12.440   2.012   (   5.366    6.926    1.745)    8.934   2.160   (  14.434    2.164   -1.087)   14.635   2.361   (   4.906    2.843   -4.910)    7.501   2.584   (   5.297    3.889    1.164)    6.674   3.116   (  -3.758   -4.689    1.856)    6.289   3.534   (   2.849   18.148    0.390)   18.374   4.012   (  22.848   -0.210   -2.335)   22.968   5.616   ( -11.050   -4.414    1.973)   12.062   8.606   (  17.616  -17.330    0.945)   24.729   8.661   (   2.696  -18.969    1.366)   19.208   9.258   (  -3.757    0.049    0.382)    3.777   9.765   (  23.167   15.920    1.523)   28.151  11.494   (  20.229   18.328  -21.092)   34.496  11.842   (  -4.718   -7.427   -1.999)    9.023  12.619   ( -13.488   -2.504   17.393)   22.152  12.686   (  -7.935  -10.814    2.918)   13.727  12.962   (  -8.491  -11.905    1.019)   14.658  13.720   (  14.787   15.601    5.227)   22.122  14.956   (  -5.305    3.426  -10.578)   12.319  15.191   (  -4.646    6.635  -11.334)   13.931======================= Grid point 185 (38/84) =======================q-point: ( 0.23  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.141   (  19.987    5.632    2.790)   20.952   1.445   (  19.085   -1.843    1.878)   19.266   1.943   (   3.957    7.343   -2.638)    8.748   2.040   (   4.111    9.826    1.848)   10.811   2.244   (  12.734    5.834   -3.222)   14.372   2.322   (   4.179   -0.506   -0.682)    4.265   2.487   (   6.117    6.259   -5.359)   10.262   2.760   (  11.781    1.464    6.451)   13.511   3.086   (   0.747   -4.367   -0.647)    4.477   3.620   (   7.170   20.195   -0.507)   21.436   4.490   (  25.024   -0.452   -1.936)   25.103   5.408   (  -9.459   -3.056    1.768)   10.096   8.633   (  -8.388  -11.756    3.190)   14.790   9.040   (  25.698  -25.801    0.488)   36.419   9.301   (  10.413  -19.869   -0.335)   22.435  10.183   (  19.377    9.526    1.085)   21.619  11.611   (  -4.461   -2.515   -3.011)    5.941  11.896   (   7.742    3.669   -6.786)   10.929  12.354   ( -12.642    2.474    9.126)   15.787  12.522   (  -9.180  -14.236   -0.035)   16.939  12.771   ( -11.116  -12.893    1.553)   17.094  13.998   (  13.345   16.958    7.163)   22.737  14.819   (  -8.328    0.327   -1.462)    8.461  15.102   (  -4.699   18.444  -14.389)   23.860======================= Grid point 186 (39/84) =======================q-point: ( 0.31  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.542   (  20.474    4.811    2.260)   21.153   1.761   (  13.272    1.435    0.802)   13.373   2.016   (   3.425    6.019   -1.708)    7.133   2.103   (   2.618   11.878    1.719)   12.284   2.335   (   1.458    1.166   -2.098)    2.808   2.439   (   3.284    0.344   -0.664)    3.369   2.566   (   3.416    8.157   -5.252)   10.286   2.969   (   8.016   -2.328    5.664)   10.088   3.122   (   2.100   -2.323    0.006)    3.132   3.830   (  13.643   17.933   -2.888)   22.717   4.954   (  21.765   -0.090   -1.347)   21.807   5.271   (  -4.180   -3.166    1.718)    5.517   8.364   ( -16.972   -1.600    1.354)   17.101   9.573   (  28.190  -29.365    0.207)   40.707   9.622   (  20.404  -37.167    0.270)   42.400  10.514   (  14.330    0.184   -0.109)   14.331  11.513   (  -4.382   -9.471   -0.630)   10.455  11.953   (  -3.361    5.821    0.742)    6.763  12.148   (  -9.209   -3.653   -1.002)    9.957  12.349   (  -6.836  -11.415   -0.403)   13.311  12.534   ( -12.823   -9.782    1.485)   16.196  14.227   (   9.849   18.243    6.845)   21.833  14.655   (  -7.826    1.973    0.186)    8.073  15.004   (  -5.224   22.967   -9.467)   25.385======================= Grid point 187 (40/84) =======================q-point: ( 0.38  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.907   (  15.297    3.865    2.332)   15.949   1.968   (   8.384    3.973    0.488)    9.291   2.077   (   3.014    5.228   -0.246)    6.040   2.144   (   1.691   13.820    1.068)   13.964   2.361   (   1.323    1.217   -2.089)    2.756   2.447   (  -0.856    3.708    0.440)    3.831   2.629   (   2.521    5.135   -5.787)    8.137   3.059   (   1.926   -2.804    2.898)    4.469   3.144   (  -0.287   -0.312    2.451)    2.488   4.097   (  12.235   13.596   -5.605)   19.130   5.240   (   0.957   -3.587    2.549)    4.504   5.318   (  14.514    0.073   -0.674)   14.530   8.038   ( -15.066    2.220   -2.390)   15.415  10.038   (  21.119  -41.469    0.721)   46.542  10.112   (  26.255  -28.420    0.164)   38.692  10.729   (   7.186  -13.901   -0.689)   15.664  11.454   (  -1.766  -14.834    0.397)   14.944  11.824   (  -7.649    3.883    0.580)    8.597  11.907   ( -15.393  -18.388   -1.046)   24.003  12.286   (  -4.867   -2.670   -0.166)    5.554  12.296   (  -6.929   -0.431    0.322)    6.950  14.375   (   5.272   19.163    5.534)   20.631  14.531   (  -4.649    5.905    0.180)    7.518  14.908   (  -4.315   24.480   -4.687)   25.295======================= Grid point 188 (41/84) =======================q-point: ( 0.46  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.058   (   3.133    3.922   -0.382)    5.034   2.078   (   1.135    4.264   -0.801)    4.485   2.151   (   1.553   12.672    1.629)   12.870   2.194   (   4.225    6.672    4.510)    9.094   2.383   (   0.539    4.904   -2.238)    5.418   2.445   (   1.039    2.130    0.613)    2.448   2.657   (   0.432    4.659   -6.165)    7.740   3.074   (   0.124   -1.295    4.388)    4.577   3.118   (  -1.447   -1.025    0.192)    1.783   4.268   (   4.649   11.298   -7.349)   14.257   5.267   (   1.113   -4.029    3.911)    5.724   5.515   (   5.220    0.344   -0.427)    5.248   7.823   (  -6.008    3.836   -5.347)    8.911  10.391   (  12.682  -33.382    0.426)   35.712  10.540   (  14.624  -21.451    0.379)   25.965  10.780   (  -1.298  -36.471    0.128)   36.494  11.441   (   0.461  -31.712    0.082)   31.715  11.613   ( -11.887  -15.289    0.681)   19.378  11.704   (  -3.707   -0.979   -1.061)    3.978  12.134   (  -4.418   15.121    0.329)   15.757  12.283   (  -0.212   -3.045   -0.369)    3.074  14.440   (   1.833   16.455    2.628)   16.765  14.472   (  -1.804   12.093    1.743)   12.350  14.848   (  -1.572   24.950   -1.852)   25.068======================= Grid point 197 (42/84) =======================q-point: ( 0.15  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.993   (  14.216   14.216    3.296)   20.373   1.146   (  13.352   13.352    3.730)   19.248   2.014   (   5.062    5.062   -2.067)    7.451   2.094   (  11.759   11.759    0.967)   16.658   2.208   (   7.453    7.453    0.195)   10.541   2.255   (  11.530   11.530   -1.400)   16.366   2.396   (   2.245    2.245   -3.015)    4.378   2.709   (   6.557    6.557    0.165)    9.274   3.065   (  -0.137   -0.137   -1.121)    1.137   3.908   (  14.442   14.442    0.277)   20.426   4.073   (  11.950   11.950   -1.466)   16.964   5.491   (  -8.134   -8.134    2.023)   11.679   8.369   (  -9.865   -9.865    1.468)   14.028   8.403   (  -1.608   -1.608    1.111)    2.531   9.192   (  -6.282   -6.282    0.014)    8.884  10.053   (  12.577   12.577    2.734)   17.996  11.642   ( -11.067  -11.067   -1.403)   15.714  11.837   (  13.995   13.995  -13.149)   23.761  12.453   ( -11.102  -11.102    0.288)   15.703  12.570   (  -4.606   -4.606   15.796)   17.087  12.697   ( -14.128  -14.128    0.528)   19.987  14.047   (  17.113   17.113    7.569)   25.357  14.911   (  -5.393   -5.393   -1.348)    7.745  15.378   (   8.676    8.676  -18.901)   22.535======================= Grid point 198 (43/84) =======================q-point: ( 0.23  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.320   (  18.603   12.438    1.591)   22.434   1.482   (  19.064    6.575    2.236)   20.289   2.082   (   2.838    7.402   -1.361)    8.043   2.298   (   7.575    0.710    0.229)    7.612   2.311   (   5.701    8.366    0.057)   10.123   2.364   (   1.556   11.367   -0.352)   11.478   2.604   (  13.337    4.899   -2.731)   14.469   2.792   (   3.641    0.772    2.809)    4.663   3.103   (   2.815    5.880   -4.223)    7.768   4.064   (   7.474   24.355   -0.268)   25.478   4.498   (  24.972    1.847   -0.850)   25.054   5.319   (  -8.699   -5.999    1.864)   10.730   8.275   (  -0.894  -21.542    2.294)   21.682   8.532   (  13.651  -22.534    0.484)   26.350   9.032   (  -8.960   -9.583   -0.343)   13.124  10.253   (   8.046   -1.233    1.876)    8.353  11.413   ( -11.679  -15.169   -1.257)   19.185  12.031   (   4.533    6.518   -7.441)   10.881  12.242   ( -10.101  -13.593    0.365)   16.939  12.349   ( -16.329   -9.111    2.597)   18.878  12.531   (  -3.878   -5.177    7.163)    9.651  14.367   (  15.209   19.487    7.527)   25.840  14.800   (  -5.333   -1.188    0.206)    5.467  15.480   (   1.857   17.137  -14.207)   22.337======================= Grid point 199 (44/84) =======================q-point: ( 0.31  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.696   (  19.425   10.885    1.036)   22.291   1.832   (  15.868    6.307    0.762)   17.092   2.137   (   2.813    6.938   -0.795)    7.529   2.352   (   1.188    0.734   -1.897)    2.356   2.425   (   4.002   18.198    2.221)   18.765   2.442   (   5.142    0.286    0.175)    5.153   2.737   (   2.539    7.196   -3.387)    8.349   2.934   (   8.197    0.223    1.670)    8.369   3.143   (   1.270    4.546   -1.794)    5.050   4.242   (  10.348   23.126   -1.489)   25.379   4.964   (  21.978    1.486   -0.644)   22.038   5.181   (  -4.900   -6.052    1.772)    7.986   8.247   (  -4.810  -13.423    1.843)   14.377   8.894   (  22.001  -36.308   -0.151)   42.454   8.939   (   3.150  -28.786   -0.536)   28.963  10.360   (   2.315  -13.278    0.058)   13.479  11.194   ( -10.999  -21.496   -0.717)   24.157  11.922   ( -13.260   -9.320   -1.674)   16.294  11.990   ( -17.318  -14.645    1.172)   22.710  12.175   (  -1.064   -7.830   -2.163)    8.193  12.477   (  -1.740    6.095    4.465)    7.753  14.621   (   9.796   19.447    5.813)   22.538  14.723   (  -1.616    6.061    1.084)    6.366  15.469   (  -2.320   21.972   -8.956)   23.841======================= Grid point 200 (45/84) =======================q-point: ( 0.38  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.038   (  13.273    9.192   -0.062)   16.145   2.084   (  10.262    8.042   -0.112)   13.038   2.197   (   4.009    7.227    1.000)    8.325   2.380   (   2.103    1.088   -1.882)    3.025   2.486   (   1.877   14.646    0.704)   14.782   2.515   (   2.690    7.165    0.539)    7.673   2.782   (   2.217    8.894   -3.258)    9.728   3.030   (   2.091    1.035   -0.609)    2.411   3.151   (  -0.884    1.770    1.698)    2.607   4.438   (   8.937   20.272   -2.908)   22.344   5.135   (  -0.069   -7.140    2.357)    7.520   5.330   (  14.936    1.382   -0.395)   15.005   8.075   ( -10.543    0.449   -1.436)   10.649   9.155   (  18.200  -46.652    0.335)   50.078   9.396   (  28.040  -41.976   -0.162)   50.481  10.309   (  -8.464  -25.563   -0.766)   26.938  10.952   ( -15.098  -31.346   -0.274)   34.793  11.643   ( -13.592  -10.916   -0.059)   17.433  11.684   ( -12.602  -16.959    0.325)   21.131  12.150   (  -1.563  -10.772   -1.005)   10.931  12.457   (  -0.511   12.666    2.133)   12.855  14.665   (  -0.688    9.878   -0.711)    9.928  14.815   (   6.235   21.590    6.532)   23.403  15.412   (  -2.990   24.503   -5.427)   25.274======================= Grid point 201 (46/84) =======================q-point: ( 0.46  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.168   (   2.214    6.940   -1.165)    7.377   2.200   (   2.202    8.568   -1.508)    8.974   2.304   (   3.327    6.081    2.307)    7.306   2.460   (   5.669    4.125   -1.202)    7.113   2.504   (   0.434   11.125   -0.172)   11.135   2.543   (  -0.266   14.500    1.083)   14.543   2.813   (   0.716   10.994   -3.276)   11.494   3.056   (   1.417   -0.049    0.644)    1.558   3.105   (  -3.138    0.343    1.070)    3.333   4.564   (   3.480   18.210   -3.882)   18.941   5.150   (   0.864   -7.838    3.255)    8.531   5.530   (   5.291    1.394   -0.229)    5.476   7.913   (  -4.769    4.524   -3.879)    7.632   9.576   (  23.353  -48.737    0.426)   54.045   9.959   (  27.724  -37.724    0.035)   46.815  10.023   ( -19.781  -37.404   -0.256)   42.313  10.568   ( -23.233  -44.468   -0.196)   50.172  11.456   (  -4.835   -9.618    0.969)   10.808  11.526   (  -3.740  -16.251   -0.503)   16.683  12.125   (  -0.733  -12.591   -0.593)   12.626  12.452   (  -0.105   15.425    1.067)   15.462  14.658   (  -0.073   12.139   -1.076)   12.187  14.896   (   2.017   22.611    6.241)   23.543  15.369   (  -1.244   25.557   -3.671)   25.849======================= Grid point 211 (47/84) =======================q-point: ( 0.23  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.609   (  16.728   16.728    0.592)   23.664   1.699   (  14.859   14.859    1.983)   21.106   2.211   (   4.880    4.880   -1.402)    7.042   2.325   (   2.565    2.565   -0.309)    3.640   2.370   (   2.368    2.368   -0.482)    3.384   2.625   (   7.536    7.536    3.366)   11.176   2.761   (  13.800   13.800   -1.482)   19.572   2.795   (   0.688    0.688    0.953)    1.362   3.228   (   4.817    4.817   -6.811)    9.632   4.510   (  15.832   15.832    0.206)   22.390   4.612   (  14.953   14.953   -0.526)   21.154   5.182   (  -7.722   -7.722    1.880)   11.082   7.926   ( -12.219  -12.219    1.805)   17.374   8.237   (  -7.009   -7.009    0.211)    9.914   8.831   ( -12.240  -12.240   -0.540)   17.318  10.191   (  -4.127   -4.127    1.506)    6.027  11.087   ( -17.379  -17.379   -0.930)   24.596  11.980   ( -12.566  -12.566    0.614)   17.782  12.007   ( -21.189  -21.189    0.516)   29.970  12.084   (  -2.734   -2.734   -5.574)    6.783  12.564   (   3.997    3.997    6.664)    8.739  14.719   (  14.506   14.506    5.612)   21.268  14.801   (   3.000    3.000    1.397)    4.466  15.726   (   8.064    8.064  -10.625)   15.587======================= Grid point 212 (48/84) =======================q-point: ( 0.31  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.953   (  17.845   14.932   -0.242)   23.269   2.013   (  15.967   12.086    0.977)   20.049   2.276   (   2.691    6.861   -0.729)    7.406   2.369   (   1.651    1.106   -1.678)    2.601   2.445   (   4.947    0.515    1.329)    5.148   2.735   (   5.175   11.196    2.796)   12.646   2.823   (   2.851    3.113   -0.514)    4.253   3.037   (  11.499   10.928   -0.634)   15.876   3.267   (  -0.466    6.577   -5.265)    8.438   4.700   (   8.497   23.042   -0.478)   24.563   5.023   (  19.594    3.870   -0.045)   19.972   5.042   (  -4.483   -7.054    1.752)    8.540   7.831   (   2.416  -24.116    1.378)   24.276   8.263   (   9.518  -26.044   -0.093)   27.729   8.557   ( -14.705  -14.371   -0.797)   20.577  10.089   (  -7.407  -12.780    0.335)   14.775  10.743   ( -18.207  -23.677   -0.415)   29.871  11.581   ( -21.176  -24.994    0.363)   32.761  11.715   ( -14.878  -11.309    0.072)   18.688  11.978   (  -6.931  -13.723   -2.609)   15.594  12.637   (   3.175    8.021    3.977)    9.499  14.815   (   2.638    6.777   -0.445)    7.285  15.021   (  11.974   17.396    6.098)   21.981  15.813   (   1.540   12.981   -7.176)   14.912======================= Grid point 213 (49/84) =======================q-point: ( 0.38  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.242   (   8.372   10.729   -1.740)   13.720   2.282   (  10.751   12.065   -0.742)   16.177   2.348   (   5.247    7.544    0.850)    9.229   2.413   (   4.829    3.005   -1.074)    5.788   2.535   (   4.085    1.398    1.353)    4.525   2.830   (   4.209   12.445    2.260)   13.330   2.905   (   5.483    3.930   -1.577)    6.928   3.161   (   0.644   12.671   -2.436)   12.919   3.242   (  -1.215    7.815   -0.748)    7.944   4.850   (   6.310   20.407   -1.046)   21.386   4.969   (  -1.102   -8.861    2.011)    9.153   5.378   (  14.551    3.634   -0.028)   14.998   7.960   (   6.180  -17.201    0.014)   18.277   8.347   (  -1.678  -30.487   -0.765)   30.542   8.595   (  23.620  -37.285    0.024)   44.137   9.844   ( -18.746  -19.847   -0.396)   27.304  10.349   ( -23.186  -29.544    0.018)   37.556  11.240   ( -12.539  -27.585    0.251)   30.302  11.414   ( -14.610  -13.012   -0.058)   19.565  11.859   (  -4.724  -19.186   -1.136)   19.792  12.684   (   1.706    9.782    2.392)   10.214  14.869   (   2.653   10.849   -1.027)   11.216  15.206   (   6.920   17.691    5.569)   19.796  15.812   (  -0.986   16.098   -4.835)   16.837======================= Grid point 214 (50/84) =======================q-point: ( 0.46  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.311   (   1.116    7.253   -1.062)    7.415   2.393   (   0.777    3.339   -0.742)    3.508   2.409   (   2.900   12.454   -1.546)   12.881   2.545   (   2.810    2.386   -1.687)    4.054   2.638   (   6.286   10.098    2.685)   12.194   2.887   (   1.435   12.275    1.718)   12.477   3.020   (   5.715    5.499   -2.084)    8.200   3.114   (  -3.912    9.650   -2.729)   10.765   3.204   (  -2.480   12.572    1.545)   12.907   4.924   (   1.413    7.349   -0.220)    7.487   4.986   (   1.559    1.627    1.224)    2.564   5.574   (   5.167    3.095    0.021)    6.023   7.946   (  -3.065   -3.141   -1.815)    4.749   8.588   (  18.503  -49.559   -0.414)   52.902   9.142   (  30.299  -41.818    0.367)   51.642   9.405   ( -21.702  -24.322   -0.236)   32.598   9.824   ( -29.379  -32.376    0.150)   43.719  11.095   (  -3.101  -27.516    0.139)   27.690  11.206   (  -5.695  -15.716    0.235)   16.717  11.798   (  -1.530  -20.929   -0.887)   21.004  12.706   (   0.526   10.431    1.658)   10.575  14.908   (   1.097   12.889   -1.322)   13.003  15.294   (   2.196   17.706    5.280)   18.606  15.793   (  -0.653   17.536   -3.682)   17.931======================= Grid point 225 (51/84) =======================q-point: ( 0.31  0.31  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.258   (  15.358   15.358   -1.693)   21.786   2.281   (  14.153   14.153    0.963)   20.038   2.388   (   4.577    4.577   -1.235)    6.589   2.395   (   1.574    1.574   -0.990)    2.436   2.483   (   4.190    4.190    2.773)    6.542   2.859   (   4.322    4.322   -1.770)    6.363   2.919   (   5.921    5.921    3.549)    9.095   3.300   (  13.467   13.467   -0.721)   19.059   3.357   (   2.344    2.344   -5.213)    6.177   4.878   (  -7.874   -7.874    1.769)   11.275   5.099   (  13.811   13.811    0.285)   19.533   5.167   (  12.652   12.652   -0.247)   17.894   7.478   (  -9.898   -9.898    0.237)   14.000   7.940   (  -6.114   -6.114    0.130)    8.647   8.260   ( -16.427  -16.427   -0.428)   23.235   9.873   ( -10.731  -10.731    0.646)   15.189  10.286   ( -23.256  -23.256   -0.047)   32.890  11.059   ( -26.796  -26.796    0.032)   37.895  11.490   ( -12.966  -12.966    0.351)   18.340  11.650   ( -17.736  -17.736   -2.202)   25.178  12.764   (   4.917    4.917    3.090)    7.609  14.952   (   6.950    6.950   -0.982)    9.878  15.311   (  12.115   12.115    5.127)   17.884  15.986   (   5.346    5.346   -5.296)    9.230======================= Grid point 226 (52/84) =======================q-point: ( 0.38  0.31  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.409   (   1.681    5.246   -0.965)    5.592   2.428   (   1.935    3.461   -0.468)    3.992   2.498   (   6.056    2.407   -1.459)    6.679   2.550   (  10.025   14.239    0.076)   17.414   2.652   (  13.260   12.636    2.464)   18.481   2.959   (   5.893    3.546   -1.968)    7.153   3.021   (   4.085    7.386    2.936)    8.936   3.363   (  -1.170    6.764   -3.787)    7.840   3.503   (   6.443   16.625   -0.455)   17.836   4.755   (  -4.155  -11.047    1.596)   11.910   5.252   (   5.509   18.450   -0.043)   19.255   5.467   (  13.368    5.151    0.157)   14.327   7.472   (   8.781  -23.241   -0.901)   24.861   7.949   ( -14.203  -17.547   -0.141)   22.575   8.036   (  15.353  -18.155    0.550)   23.783   9.578   ( -21.344   -7.963    0.617)   22.789   9.816   ( -25.644  -24.924    0.153)   35.761  10.629   ( -14.972  -33.985    0.233)   37.138  11.146   ( -19.438  -14.007   -0.489)   23.964  11.410   (  -7.477  -26.134   -1.088)   27.204  12.833   (   2.402    5.639    2.430)    6.593  15.073   (   5.145    9.635   -1.478)   11.023  15.494   (   6.747   11.815    4.398)   14.299  16.045   (   1.273    8.176   -3.900)    9.148======================= Grid point 227 (53/84) =======================q-point: ( 0.46  0.31  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.428   (   0.432    3.480   -0.446)    3.535   2.456   (   0.816    4.561    0.104)    4.635   2.564   (   1.504    0.540   -0.683)    1.738   2.681   (   3.240   14.594   -0.991)   14.982   2.889   (   8.653   14.646    0.937)   17.037   3.054   (   1.117    6.579    1.743)    6.897   3.071   (   2.492    3.101   -0.593)    4.022   3.331   (  -1.442   10.197   -2.734)   10.655   3.557   (   0.225   19.641   -0.222)   19.643   4.716   (  -0.506  -13.417    1.390)   13.498   5.327   (   2.045   17.603   -0.104)   17.721   5.647   (   4.773    4.212    0.214)    6.370   7.728   (   8.961  -27.008   -1.480)   28.494   7.754   (  -1.486  -25.241   -0.434)   25.288   8.466   (  24.345  -25.196    1.190)   35.056   9.121   ( -19.015   -5.337    0.471)   19.756   9.285   ( -26.084  -22.806   -0.067)   34.648  10.472   (  -2.937  -35.578    0.800)   35.708  10.860   (  -8.060  -19.066   -0.433)   20.704  11.315   (  -2.419  -28.122   -1.038)   28.245  12.863   (   0.702    5.879    2.064)    6.270  15.143   (   1.883   10.844   -1.777)   11.149  15.580   (   2.117   11.582    4.030)   12.445  16.055   (   0.074    9.628   -3.176)   10.139======================= Grid point 239 (54/84) =======================q-point: ( 0.38  0.38  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.452   (   1.266    1.266    0.186)    1.800   2.513   (   2.940    2.940   -0.450)    4.182   2.536   (   2.530    2.530   -0.018)    3.578   2.811   (  11.126   11.126    0.739)   15.752   2.929   (  14.012   14.012    0.865)   19.834   3.035   (   4.454    4.454   -1.551)    6.488   3.154   (   5.926    5.926    2.619)    8.780   3.458   (   3.203    3.203   -3.338)    5.627   3.782   (  10.745   10.745   -1.090)   15.234   4.554   (  -8.755   -8.755    1.054)   12.427   5.548   (   8.775    8.775    0.373)   12.416   5.572   (   7.885    7.885    0.131)   11.152   7.205   (  -3.746   -3.746   -2.363)    5.800   7.610   ( -15.804  -15.804    0.460)   22.355   7.921   (   6.173    6.173    1.123)    8.802   9.364   ( -22.136  -22.136    0.077)   31.305   9.441   ( -10.433  -10.433    1.833)   14.868  10.000   ( -25.840  -25.840    0.357)   36.544  10.857   ( -20.612  -20.612   -1.386)   29.182  10.918   ( -16.576  -16.576   -0.437)   23.446  12.912   (   2.664    2.664    2.361)    4.446  15.233   (   6.519    6.519   -1.947)    9.422  15.671   (   6.422    6.422    3.491)    9.730  16.151   (   3.157    3.157   -3.105)    5.438======================= Grid point 240 (55/84) =======================q-point: ( 0.46  0.38  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.468   (   0.432    1.398    0.572)    1.571   2.547   (   0.869    3.858    0.623)    4.003   2.579   (   1.278    1.161   -0.017)    1.726   2.950   (   3.165   11.950    0.125)   12.363   3.105   (   3.233    2.292   -0.460)    3.990   3.160   (   6.461    9.717    2.405)   11.914   3.237   (   2.368   12.371    0.900)   12.627   3.496   (   0.857    6.581   -3.082)    7.317   3.897   (   1.631   14.413   -1.335)   14.567   4.445   (  -2.254  -13.798    0.259)   13.984   5.618   (   1.638   11.811    0.309)   11.928   5.726   (   4.326    3.540    0.310)    5.598   7.260   (   5.652  -13.493   -3.476)   15.037   7.383   (  -6.735  -17.378    0.639)   18.648   8.201   (  17.795   -2.161    1.772)   18.013   8.935   ( -19.976  -13.229   -0.395)   23.963   9.112   ( -15.433    2.123    1.838)   15.686   9.756   (  -3.044  -35.812    1.531)   35.974  10.494   ( -10.894  -17.353   -0.881)   20.508  10.737   (  -4.715  -29.517   -0.979)   29.908  12.944   (   0.756    2.730    2.283)    3.638  15.320   (   2.294    7.091   -2.240)    7.782  15.751   (   1.991    6.214    3.012)    7.187  16.185   (   0.722    4.200   -2.669)    5.028======================= Grid point 253 (56/84) =======================q-point: ( 0.46  0.46  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.487   (   0.470    0.470    1.180)    1.354   2.598   (   1.230    1.230    0.787)    1.909   2.601   (   0.858    0.858    0.908)    1.515   3.113   (   3.965    3.965    0.772)    5.661   3.138   (   1.376    1.376    0.357)    1.978   3.274   (   2.796    2.796    3.241)    5.112   3.482   (   9.147    9.147    0.699)   12.955   3.584   (   2.314    2.314   -3.240)    4.606   4.090   (   4.372    4.372   -2.125)    6.538   4.217   (  -7.353   -7.353   -1.255)   10.475   5.775   (   2.882    2.882    0.400)    4.095   5.777   (   2.618    2.618    0.386)    3.722   7.122   (  -5.848   -5.848    0.657)    8.296   7.175   (  -0.711   -0.711   -4.799)    4.903   8.304   (   8.371    8.371    2.230)   12.047   8.740   (  -7.051   -7.051   -0.872)   10.010   9.134   (  -4.437   -4.437    2.815)    6.878   9.211   ( -12.064  -12.064    2.377)   17.227  10.228   (  -9.395   -9.395   -1.011)   13.325  10.276   ( -13.079  -13.079   -0.887)   18.518  12.977   (   0.761    0.761    2.364)    2.598  15.414   (   2.454    2.454   -2.523)    4.291  15.829   (   1.916    1.916    2.458)    3.658  16.234   (   1.096    1.096   -2.379)    2.839======================= Grid point 338 (57/84) =======================q-point: ( 0.00  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 84Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.585   (   0.000    0.000    6.529)    6.529   0.585   (  -0.000   -0.000    6.529)    6.529   1.146   (   0.000    0.000   14.081)   14.081   1.544   (  -0.000    0.000   -3.231)    3.231   1.544   (   0.000    0.000   -3.231)    3.231   1.920   (   0.000    0.000    0.556)    0.556   1.920   (  -0.000   -0.000    0.556)    0.556   2.674   (  -0.000    0.000  -10.010)   10.010   3.256   (  -0.000    0.000    0.096)    0.096   3.256   (   0.000    0.000    0.096)    0.096   3.747   (   0.000    0.000    2.949)    2.949   6.019   (   0.000    0.000    0.065)    0.065   8.421   (  -0.000   -0.000    0.731)    0.731   8.421   (   0.000    0.000    0.731)    0.731   8.738   (   0.000    0.000   -0.838)    0.838   8.738   (  -0.000   -0.000   -0.838)    0.838   9.386   (   0.000    0.000    0.276)    0.276  12.035   (   0.000    0.000   -0.175)    0.175  12.941   (   0.000    0.000    0.260)    0.260  12.941   (   0.000    0.000    0.260)    0.260  13.246   (   0.000    0.000   -0.164)    0.164  13.246   (  -0.000   -0.000   -0.164)    0.164  15.342   (  -0.000    0.000    2.423)    2.423  15.387   (   0.000    0.000   -1.058)    1.058======================= Grid point 339 (58/84) =======================q-point: ( 0.08  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.803   (  18.436    0.000    4.532)   18.985   0.871   (   7.170    0.000   12.273)   14.214   1.038   (  12.339    0.000   -1.699)   12.456   1.583   (   3.731    0.000   -2.988)    4.780   1.756   (  16.327    0.000   -1.391)   16.386   1.930   (   0.972    0.000    0.523)    1.104   2.138   (  16.646    0.000    4.685)   17.293   2.518   (  -9.165    0.000   -8.161)   12.272   3.249   (  -0.629    0.000    0.077)    0.634   3.383   (   4.574    0.000    0.039)    4.574   3.650   (   1.818    0.000    1.580)    2.409   5.916   (  -9.499    0.000    0.659)    9.521   8.539   (  11.878    0.000    0.701)   11.899   8.849   (  11.245    0.000   -0.788)   11.273   8.910   (  14.122    0.000    0.303)   14.125   9.383   (   1.111    0.000    0.212)    1.131   9.506   (  81.840    0.000  -12.475)   82.785  11.990   (  -4.231    0.000   -0.351)    4.245  12.913   (  -2.827    0.000    0.213)    2.835  13.103   (  11.555    0.000    4.488)   12.396  13.210   (  -3.665    0.000   -0.125)    3.667  13.220   (   0.931    0.000   -2.577)    2.740  15.019   ( -27.957    0.000    0.193)   27.957  15.313   (  -6.638    0.000   -0.300)    6.644======================= Grid point 340 (59/84) =======================q-point: ( 0.15  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.899   (   7.390    0.000    5.292)    9.089   1.186   (  19.316    0.000    2.638)   19.495   1.566   (  28.275    0.000   -2.691)   28.403   1.678   (   5.667    0.000   -2.380)    6.146   1.956   (   1.623    0.000    0.431)    1.679   2.092   (  16.284    0.000   -2.057)   16.413   2.327   (   1.259    0.000    2.548)    2.842   2.557   (  14.208    0.000   -1.266)   14.264   3.232   (  -1.183    0.000    0.028)    1.184   3.331   (  -3.797    0.000    1.320)    4.020   3.959   (  23.217    0.000   -1.732)   23.282   5.694   ( -12.208    0.000    0.974)   12.247   8.870   (  21.514    0.000    0.611)   21.522   9.033   (  -4.820    0.000    2.042)    5.235   9.164   (  20.579    0.000   -0.644)   20.589   9.574   (  19.293    0.000    0.255)   19.294  10.830   (  48.899    0.000  -10.242)   49.960  11.892   (  -5.027    0.000   -0.742)    5.082  12.830   (  -5.634    0.000    0.086)    5.635  12.945   ( -18.963    0.000    6.677)   20.104  13.102   (  -7.399    0.000   -0.025)    7.399  13.596   (  24.629    0.000    2.865)   24.795  14.544   ( -15.455    0.000   -9.172)   17.972  15.133   ( -11.500    0.000   -0.358)   11.505======================= Grid point 341 (60/84) =======================q-point: ( 0.23  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.164   (  18.148    0.000    2.311)   18.294   1.526   (  15.184    0.000    1.565)   15.264   1.791   (   5.691    0.000   -1.622)    5.917   1.990   (   1.780    0.000    0.298)    1.805   2.024   (  17.399    0.000   -3.447)   17.737   2.326   (   7.073    0.000   -0.140)    7.075   2.369   (   3.696    0.000   -1.445)    3.969   2.876   (  16.598    0.000    3.462)   16.955   3.204   (  -1.553    0.000   -0.027)    1.553   3.332   (   5.029    0.000   -0.448)    5.049   4.447   (  25.610    0.000   -1.534)   25.656   5.471   (  -9.892    0.000    0.958)    9.938   8.803   ( -17.214    0.000    2.141)   17.346   9.353   (  27.185    0.000    0.479)   27.189   9.631   (  26.727    0.000   -0.422)   26.730   9.986   (  20.641    0.000    0.594)   20.650  11.454   (  16.389    0.000   -3.255)   16.709  11.846   (   2.696    0.000   -1.174)    2.941  12.492   ( -24.526    0.000    4.272)   24.895  12.694   (  -8.373    0.000   -0.083)    8.373  12.921   ( -11.239    0.000    0.089)   11.240  13.996   (  15.118    0.000   10.743)   18.546  14.439   (   1.815    0.000  -15.564)   15.670  14.884   ( -13.294    0.000   -0.469)   13.302======================= Grid point 342 (61/84) =======================q-point: ( 0.31  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.558   (  21.299    0.000    1.773)   21.373   1.775   (  10.400    0.000    0.845)   10.434   1.893   (   4.589    0.000   -0.838)    4.665   2.022   (   1.461    0.000    0.127)    1.466   2.243   (   6.196    0.000   -2.285)    6.604   2.416   (   3.317    0.000   -0.857)    3.425   2.423   (   1.025    0.000   -0.849)    1.331   3.144   (   9.672    0.000    4.547)   10.687   3.173   (  -1.603    0.000   -0.068)    1.604   3.529   (  13.714    0.000   -3.370)   14.122   4.923   (  22.433    0.000   -0.991)   22.455   5.336   (  -3.517    0.000    0.967)    3.648   8.399   ( -23.100    0.000    0.755)   23.112   9.900   (  28.096    0.000    0.343)   28.098  10.184   (  29.301    0.000   -0.150)   29.302  10.348   (  16.326    0.000    0.282)   16.329  11.568   (   0.573    0.000   -0.348)    0.671  11.968   (  -0.719    0.000    4.738)    4.793  12.094   (  -8.966    0.000   -4.481)   10.023  12.506   ( -10.788    0.000   -0.253)   10.791  12.663   ( -15.107    0.000    0.157)   15.108  14.203   (   6.665    0.000    9.973)   11.996  14.498   (   2.643    0.000  -11.406)   11.708  14.642   ( -10.786    0.000   -0.145)   10.787======================= Grid point 343 (62/84) =======================q-point: ( 0.38  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.935   (   6.031    0.000    0.263)    6.037   1.961   (  18.972    0.000    2.114)   19.089   1.967   (   3.022    0.000   -0.057)    3.023   2.044   (   0.783    0.000   -0.082)    0.787   2.312   (   1.671    0.000   -0.431)    1.726   2.425   (  -0.216    0.000    0.879)    0.905   2.466   (   1.776    0.000   -3.493)    3.918   3.145   (  -1.232    0.000   -0.086)    1.234   3.236   (   0.122    0.000    4.212)    4.214   3.794   (  11.992    0.000   -5.405)   13.154   5.296   (  14.801    0.000   -0.255)   14.803   5.332   (   2.755    0.000    1.190)    3.001   7.958   ( -20.426    0.000   -1.681)   20.495  10.415   (  23.659    0.000    0.257)   23.661  10.620   (  11.330    0.000    0.079)   11.330  10.752   (  28.304    0.000    0.152)   28.305  11.579   (   1.091    0.000    0.627)    1.258  11.731   ( -11.971    0.000    0.653)   11.989  12.162   (   4.193    0.000   -1.394)    4.419  12.282   ( -11.715    0.000   -0.348)   11.720  12.334   ( -18.386    0.000    0.061)   18.386  14.260   (  -0.047    0.000    4.082)    4.082  14.473   (  -6.576    0.000    0.247)    6.581  14.546   (   2.448    0.000   -3.562)    4.322======================= Grid point 344 (63/84) =======================q-point: ( 0.46  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.005   (   1.198    0.000    0.676)    1.375   2.018   (   2.175    0.000   -0.463)    2.224   2.051   (  -0.106    0.000   -0.087)    0.137   2.245   (   8.725    0.000    1.954)    8.941   2.323   (  -0.656    0.000    1.372)    1.521   2.436   (   1.630    0.000    0.928)    1.876   2.483   (  -0.121    0.000   -3.802)    3.804   3.128   (  -0.466    0.000   -0.087)    0.474   3.195   (  -2.508    0.000    2.770)    3.737   3.960   (   4.505    0.000   -6.051)    7.543   5.382   (   1.914    0.000    2.053)    2.807   5.502   (   5.395    0.000   -0.318)    5.405   7.665   (  -8.266    0.000   -3.559)    9.000  10.769   (  10.705    0.000    0.261)   10.708  10.779   (   4.492    0.000    0.287)    4.501  11.246   (  19.566    0.000    0.557)   19.574  11.566   (  -4.568    0.000   -1.233)    4.731  11.604   (   0.939    0.000    1.297)    1.601  11.972   ( -16.039    0.000   -0.050)   16.039  12.091   (  -6.440    0.000   -0.631)    6.470  12.207   (   0.828    0.000   -0.264)    0.869  14.256   (   0.342    0.000    1.471)    1.511  14.368   (  -3.712    0.000   -0.426)    3.737  14.588   (   1.308    0.000    0.707)    1.487======================= Grid point 352 (64/84) =======================q-point: ( 0.08  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.801   (   9.375    9.375    4.575)   14.026   0.856   (  -2.847   -2.847    7.867)    8.837   1.332   (  23.947   23.947   -0.855)   33.877   1.684   (   6.576    6.576   -3.224)    9.843   1.841   (  10.000   10.000   -1.295)   14.201   1.977   (   3.226    3.226    0.550)    4.596   2.248   (   8.721    8.721    6.670)   14.021   2.458   (  -2.138   -2.138   -8.457)    8.982   3.263   (  -4.719   -4.719    0.574)    6.698   3.441   (   9.076    9.076    0.149)   12.837   3.667   (   5.011    5.011    0.091)    7.088   5.840   (  -7.377   -7.377    0.840)   10.467   8.439   (   0.637    0.637    0.906)    1.277   8.658   (  -4.203   -4.203   -0.277)    5.950   9.312   (  23.649   23.649    0.479)   33.448   9.352   (  -2.037   -2.037    0.249)    2.892  10.265   (  54.858   54.858  -12.893)   78.645  11.934   (  -5.301   -5.301   -0.503)    7.514  12.809   (  -6.523   -6.523    0.179)    9.226  13.112   (  -6.865   -6.865    0.017)    9.708  13.168   (  -3.588   -3.588    6.338)    8.119  13.422   (  21.448   21.448   -1.038)   30.349  14.818   ( -16.421  -16.421   -4.272)   23.613  15.255   (  -5.759   -5.759   -0.349)    8.152======================= Grid point 353 (65/84) =======================q-point: ( 0.15  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.914   (   9.717    3.155    3.897)   10.934   1.116   (  19.518   -1.182    2.664)   19.734   1.789   (   7.275   10.516   -1.926)   12.931   1.806   (  20.244   17.339   -1.017)   26.674   2.035   (   3.592    7.173    0.408)    8.032   2.112   (  14.390    2.476   -2.279)   14.778   2.315   (   0.790   -0.854    0.069)    1.166   2.596   (  13.204    4.301    0.144)   13.888   3.145   (  -5.799   -5.341    0.627)    7.909   3.542   (   2.427   17.244    0.256)   17.416   3.966   (  22.784    0.926   -1.352)   22.842   5.653   ( -11.012   -4.279    1.027)   11.858   8.627   (  17.071  -16.800    0.705)   23.961   8.687   (   5.147  -19.540    0.696)   20.218   9.265   (  -4.415   -0.306    0.236)    4.431   9.797   (  23.444   17.721    0.957)   29.403  11.169   (  34.606   26.420   -7.755)   44.224  11.811   (  -6.344   -7.897   -0.849)   10.165  12.678   (  -7.259  -11.298    0.173)   13.431  12.891   ( -17.153   -5.791    3.637)   18.466  13.015   (  -9.902   -8.472    3.435)   13.477  13.894   (  23.062   23.697    8.922)   34.249  14.579   (  -4.655    4.705  -15.788)   17.119  15.098   ( -10.057   -3.361   -0.874)   10.640======================= Grid point 354 (66/84) =======================q-point: ( 0.23  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.200   (  18.288    4.450    1.862)   18.913   1.487   (  17.493   -1.596    1.403)   17.622   1.890   (   5.038    8.021   -1.686)    9.621   2.082   (   3.927   10.657    1.891)   11.514   2.169   (  10.886   10.330   -2.950)   15.294   2.323   (   6.659   -0.299    0.939)    6.731   2.385   (   6.275    2.897   -3.627)    7.805   2.874   (  14.371   -0.091    3.355)   14.757   3.067   (  -1.619   -6.431   -0.901)    6.693   3.611   (   5.537   20.351   -0.246)   21.092   4.451   (  25.595    0.560   -1.197)   25.629   5.442   (  -9.690   -2.955    0.992)   10.179   8.699   (  -8.268  -12.426    2.025)   15.062   9.050   (  25.237  -25.895    0.343)   36.161   9.292   (  10.478  -19.925   -0.338)   22.515  10.205   (  18.288   10.368    0.683)   21.034  11.551   (   2.102    2.804   -1.779)    3.930  11.799   (   9.539   -0.634   -2.178)    9.804  12.493   ( -20.261   -0.901    3.303)   20.548  12.521   (  -9.783  -13.391   -0.066)   16.584  12.807   ( -11.630   -9.949    1.259)   15.356  14.241   (  12.907   19.121   14.999)   27.516  14.625   (   6.306   17.152  -17.856)   25.550  14.881   ( -11.526    0.632   -2.208)   11.753======================= Grid point 355 (67/84) =======================q-point: ( 0.31  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.588   (  20.636    3.807    1.441)   21.033   1.779   (  12.257    1.109    0.565)   12.320   1.984   (   4.562    6.617   -0.922)    8.090   2.136   (   2.065   12.846    1.066)   13.054   2.288   (   2.905    2.477   -1.616)    4.146   2.429   (   3.367    1.103    0.109)    3.545   2.478   (   3.339    6.866   -2.555)    8.051   3.045   (   1.956   -6.932    1.122)    7.290   3.139   (   8.254   -1.825    1.276)    8.549   3.772   (  10.515   19.920   -1.792)   22.596   4.927   (  22.443    0.510   -0.798)   22.463   5.305   (  -3.856   -3.129    0.977)    5.061   8.392   ( -20.388   -1.319    0.827)   20.447   9.577   (  27.991  -29.829    0.115)   40.906   9.627   (  21.199  -37.503    0.140)   43.080  10.512   (  13.166    0.139   -0.073)   13.167  11.501   (  -3.151   -9.534   -0.348)   10.047  11.968   (   2.962   -1.195    0.462)    3.227  12.130   ( -15.699    2.687   -0.524)   15.936  12.341   (  -7.932  -12.040   -0.243)   14.420  12.563   ( -13.166   -7.362    0.835)   15.108  14.440   (   7.754   18.903   12.219)   23.807  14.659   (  -8.022    2.054    0.094)    8.281  14.739   (   1.971   22.464  -13.802)   26.439======================= Grid point 356 (68/84) =======================q-point: ( 0.38  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.950   (  11.751    2.385    1.394)   12.071   1.970   (   9.038    3.006   -0.452)    9.535   2.083   (   7.512    6.362    0.503)    9.857   2.164   (   0.881   14.106    0.591)   14.146   2.323   (   1.060    0.661   -0.959)    1.575   2.457   (   0.456    2.624    0.445)    2.700   2.531   (   2.061    7.289   -2.595)    8.007   3.069   (   0.723   -5.341   -0.462)    5.410   3.223   (   0.038   -1.570    2.704)    3.127   3.983   (   9.917   17.038   -3.529)   20.027   5.289   (   4.513   -2.995    1.436)    5.603   5.304   (  12.658   -0.444   -0.475)   12.675   7.989   ( -18.992    2.882   -1.486)   19.267  10.053   (  21.258  -41.448    0.446)   46.584  10.115   (  26.354  -28.889    0.083)   39.104  10.715   (   7.282  -14.348   -0.439)   16.096  11.461   (  -1.122  -16.341    0.197)   16.381  11.842   ( -10.647    2.920    0.820)   11.071  11.884   ( -13.095  -16.213   -0.907)   20.860  12.272   (  -0.190   -3.773   -0.552)    3.818  12.312   ( -11.693    2.196    0.535)   11.910  14.484   (  -3.459   10.676    1.894)   11.381  14.584   (   1.150   17.103    6.546)   18.350  14.766   (   1.048   21.243   -7.815)   22.659======================= Grid point 357 (69/84) =======================q-point: ( 0.46  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.050   (   2.239    3.841   -0.323)    4.457   2.069   (   1.548    2.938   -0.079)    3.322   2.171   (   0.114   14.552    0.396)   14.558   2.287   (   7.090    3.552    2.348)    8.271   2.338   (   0.587    1.707   -1.256)    2.199   2.471   (   1.106    3.362    1.508)    3.847   2.556   (   0.418    7.555   -2.640)    8.014   3.078   (   0.315   -3.823   -1.005)    3.965   3.180   (  -2.687   -1.470    2.292)    3.825   4.123   (   3.885   15.219   -4.302)   16.286   5.337   (   1.762   -4.491    1.924)    5.194   5.506   (   5.411    0.502   -0.264)    5.441   7.715   (  -7.742    4.901   -3.243)    9.719  10.400   (  11.992  -33.343    0.270)   35.435  10.547   (  14.908  -21.658    0.215)   26.294  10.782   (  -0.043  -36.845    0.030)   36.845  11.439   (  -1.683  -36.872   -0.150)   36.911  11.631   (  -9.128  -10.641    0.681)   14.037  11.683   (  -4.721   -0.290   -0.651)    4.775  12.140   (  -4.951   15.942    0.202)   16.694  12.275   (   0.092   -3.125   -0.268)    3.138  14.435   (  -1.110    8.228   -0.598)    8.324  14.586   (  -0.637   25.383    5.119)   25.902  14.790   (   0.936   19.078   -2.995)   19.334======================= Grid point 366 (70/84) =======================q-point: ( 0.15  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.061   (  11.822   11.822    2.085)   16.849   1.229   (  12.184   12.184    2.739)   17.447   1.962   (   6.381    6.381   -1.918)    9.225   2.134   (  12.483   12.483    2.606)   17.845   2.191   (   7.583    7.583   -1.787)   10.872   2.229   (  10.365   10.365   -0.670)   14.673   2.328   (   4.409    4.409   -2.341)    6.661   2.719   (   7.678    7.678    0.542)   10.872   3.043   (  -4.682   -4.682   -0.665)    6.654   3.913   (  14.417   14.417    0.177)   20.390   4.043   (  12.688   12.688   -0.902)   17.967   5.529   (  -8.130   -8.130    1.104)   11.551   8.398   (  -9.165   -9.165    0.817)   12.986   8.428   (  -1.954   -1.954    0.836)    2.887   9.192   (  -6.610   -6.610   -0.044)    9.348  10.108   (  12.870   12.870    1.663)   18.276  11.613   ( -11.133  -11.133   -0.891)   15.769  11.639   (  19.352   19.352   -4.675)   27.765  12.459   ( -10.960  -10.960    0.211)   15.502  12.707   ( -13.958  -13.958    0.307)   19.741  12.818   (  -9.515   -9.515    6.067)   14.760  14.318   (  17.640   17.640   17.317)   30.368  14.749   (  10.562   10.562  -18.213)   23.554  15.025   (  -2.422   -2.422   -6.088)    6.986======================= Grid point 367 (71/84) =======================q-point: ( 0.23  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.354   (  17.552   11.256    1.133)   20.881   1.532   (  17.208    6.795    1.667)   18.576   2.049   (   3.500    7.866   -1.204)    8.693   2.276   (   2.354    1.730   -1.213)    3.163   2.316   (   7.522    3.661    0.132)    8.367   2.397   (   6.836   15.312    2.588)   16.967   2.533   (  10.744   10.483   -3.169)   15.342   2.869   (   6.864   -0.706    3.968)    7.960   2.995   (   1.259   -0.426   -4.911)    5.087   4.059   (   6.689   24.665   -0.121)   25.556   4.481   (  25.205    2.875   -0.506)   25.373   5.356   (  -8.900   -5.911    1.060)   10.736   8.322   (   0.036  -22.365    1.454)   22.412   8.543   (  12.819  -22.604    0.381)   25.989   9.023   (  -9.659   -9.699   -0.344)   13.693  10.292   (   6.270   -0.825    1.183)    6.433  11.387   ( -11.221  -14.953   -0.779)   18.711  11.916   (   8.202    7.361   -2.815)   11.374  12.250   ( -10.159  -13.393    0.242)   16.812  12.379   ( -18.408  -13.874    0.580)   23.058  12.663   (  -6.666   -1.798    3.585)    7.780  14.604   (  12.691   17.519   12.792)   25.132  14.803   (  -3.984    0.034   -0.179)    3.988  15.111   (   8.585   20.788  -16.462)   27.872======================= Grid point 368 (72/84) =======================q-point: ( 0.31  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.719   (  19.095    9.551    0.760)   21.364   1.848   (  14.317    6.482    0.515)   15.725   2.119   (   3.859    7.264   -0.615)    8.248   2.313   (   1.547    0.848   -1.179)    2.122   2.442   (   5.379    4.726   -0.014)    7.161   2.488   (   2.571   16.295    2.712)   16.718   2.662   (   3.596    9.707   -2.810)   10.726   2.957   (   2.910   -1.307    0.341)    3.208   3.107   (   7.297   -0.872   -1.000)    7.417   4.214   (   8.894   24.368   -0.793)   25.953   4.951   (  22.228    2.222   -0.375)   22.342   5.215   (  -4.769   -5.989    1.005)    7.721   8.285   (  -7.418  -12.830    1.204)   14.869   8.893   (  21.836  -36.460    0.072)   42.499   8.924   (   3.828  -29.958   -0.589)   30.207  10.361   (   0.677  -13.000    0.048)   13.017  11.179   ( -10.403  -21.386   -0.440)   23.786  11.893   ( -10.567   -7.100   -0.800)   12.756  12.011   ( -18.365  -16.110    0.584)   24.436  12.140   (  -0.855   -9.020   -0.888)    9.103  12.558   (  -4.097    7.023    2.202)    8.424  14.709   (  -3.878    4.553    0.876)    6.044  14.832   (   9.818   18.876    8.530)   22.923  15.232   (   3.348   24.598  -10.647)   27.011======================= Grid point 369 (73/84) =======================q-point: ( 0.38  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.044   (  11.315    7.104    0.428)   13.367   2.071   (   8.426    7.532   -0.729)   11.326   2.222   (   8.144    7.769    0.828)   11.286   2.338   (   1.159    1.125   -1.399)    2.137   2.504   (   0.673    7.784    0.733)    7.847   2.535   (   3.091   16.089    0.888)   16.407   2.720   (   2.619   10.323   -1.880)   10.815   2.999   (   1.476   -0.777   -1.352)    2.147   3.190   (   0.558   -0.959    1.205)    1.638   4.383   (   7.758   22.762   -1.558)   24.098   5.179   (   0.715   -7.500    1.254)    7.638   5.322   (  15.175    1.817   -0.242)   15.285   8.046   ( -13.957    1.744   -0.916)   14.096   9.163   (  18.820  -47.408    0.249)   51.007   9.392   (  28.167  -42.071   -0.131)   50.630  10.293   (  -8.499  -25.399   -0.473)   26.788  10.947   ( -14.836  -31.126   -0.174)   34.481  11.644   ( -12.372  -10.560    0.068)   16.266  11.690   ( -13.430  -17.188    0.167)   21.813  12.130   (  -0.861  -10.928   -0.582)   10.977  12.499   (  -2.043   13.225    1.220)   13.437  14.653   (  -1.688    9.339   -0.296)    9.495  14.989   (   6.029   20.431    8.203)   22.827  15.258   (  -0.030   25.493   -7.663)   26.620======================= Grid point 370 (74/84) =======================q-point: ( 0.46  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.153   (   1.929    5.989   -0.278)    6.298   2.169   (   2.276    8.214   -0.951)    8.576   2.337   (   0.826    1.783    0.656)    2.072   2.415   (   4.805    5.919   -1.696)    7.811   2.530   (   2.984   15.665    1.197)   15.992   2.574   (   0.478   11.548    1.369)   11.638   2.758   (   1.095   12.311   -1.425)   12.441   3.019   (   0.660   -1.662   -1.866)    2.585   3.153   (  -2.673   -0.421    1.570)    3.128   4.493   (   3.059   21.486   -1.955)   21.791   5.209   (   1.332   -8.615    1.583)    8.860   5.526   (   5.387    1.646   -0.135)    5.635   7.835   (  -6.135    6.424   -2.363)    9.192   9.584   (  22.955  -49.127    0.248)   54.226   9.959   (  28.097  -37.971    0.019)   47.236  10.017   ( -19.055  -37.615   -0.185)   42.167  10.565   ( -23.329  -43.637   -0.102)   49.481  11.475   (  -4.301  -10.348    0.590)   11.222  11.516   (  -4.277  -15.628   -0.286)   16.206  12.113   (  -0.602  -12.894   -0.384)   12.914  12.474   (  -0.615   16.129    0.661)   16.154  14.636   (  -0.279   12.017   -0.660)   12.038  15.064   (   1.735   21.693    8.232)   23.267  15.252   (  -0.208   25.382   -6.667)   26.244======================= Grid point 380 (75/84) =======================q-point: ( 0.23  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.624   (  15.834   15.834    0.519)   22.398   1.743   (  14.076   14.076    1.484)   19.962   2.178   (   4.962    4.962   -1.109)    7.105   2.302   (   1.619    1.619   -1.046)    2.517   2.374   (   4.025    4.025    0.438)    5.709   2.691   (   6.169    6.169    1.694)    8.887   2.730   (  14.250   14.250   -0.930)   20.174   2.863   (   1.946    1.946    6.127)    6.717   3.043   (   4.219    4.219   -9.308)   11.057   4.515   (  15.816   15.816    0.135)   22.368   4.602   (  15.142   15.142   -0.301)   21.416   5.218   (  -7.772   -7.772    1.076)   11.044   7.963   ( -12.564  -12.564    1.135)   17.804   8.242   (  -7.450   -7.450    0.190)   10.537   8.819   ( -12.418  -12.418   -0.421)   17.567  10.222   (  -5.190   -5.190    0.957)    7.403  11.069   ( -16.889  -16.889   -0.569)   23.892  11.989   (  -1.007   -1.007   -2.471)    2.852  11.992   ( -12.484  -12.484    0.388)   17.659  12.017   ( -21.372  -21.372    0.307)   30.226  12.679   (   1.882    1.882    3.080)    4.071  14.793   (   0.140    0.140    0.753)    0.779  14.923   (  14.748   14.748    8.281)   22.441  15.456   (  13.316   13.316  -11.211)   21.916======================= Grid point 381 (76/84) =======================q-point: ( 0.31  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.950   (  16.821   13.673   -0.001)   21.677   2.035   (  14.535   12.426    0.742)   19.137   2.257   (   3.938    6.382   -0.721)    7.534   2.337   (   1.830    1.697   -0.943)    2.668   2.471   (   5.239    1.543    0.782)    5.517   2.765   (   3.315    3.445   -0.206)    4.786   2.863   (   2.797   14.520    3.633)   15.227   3.026   (   8.257    8.660   -0.019)   11.966   3.132   (   5.468    3.775   -5.686)    8.745   4.692   (   7.904   23.885   -0.210)   25.160   5.022   (  21.208    5.038   -0.049)   21.798   5.076   (  -5.993   -7.986    1.002)   10.035   7.859   (   1.734  -25.332    0.878)   25.407   8.262   (   9.396  -25.948   -0.014)   27.596   8.540   ( -15.059  -13.877   -0.573)   20.486  10.096   (  -8.565  -12.468    0.217)   15.128  10.735   ( -17.674  -23.265   -0.253)   29.218  11.589   ( -21.289  -25.493    0.260)   33.214  11.715   ( -15.356  -10.429   -0.091)   18.563  11.931   (  -4.883  -13.771   -1.251)   14.665  12.713   (   1.382    7.302    2.172)    7.742  14.811   (   1.654    6.023   -0.011)    6.246  15.176   (  10.950   15.763    6.442)   20.246  15.631   (   5.218   15.853   -7.356)   18.239======================= Grid point 382 (77/84) =======================q-point: ( 0.38  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.214   (   7.900    9.788   -0.719)   12.599   2.269   (   8.756   12.595   -0.365)   15.344   2.359   (   3.107    4.088    0.421)    5.152   2.385   (   7.323    5.349   -1.146)    9.141   2.569   (   5.151    3.443    1.210)    6.313   2.853   (   5.639    4.390   -0.777)    7.189   2.909   (   2.007   14.274    2.395)   14.612   3.087   (  -0.145    9.551   -3.317)   10.112   3.235   (   2.691    7.389   -0.082)    7.865   4.835   (   6.304   22.663   -0.355)   23.526   5.005   (  -1.216  -10.102    0.940)   10.218   5.378   (  14.622    3.937   -0.012)   15.143   7.959   (   3.454  -16.161   -0.026)   16.526   8.332   (  -0.578  -31.434   -0.432)   31.443   8.595   (  23.936  -36.835    0.001)   43.929   9.836   ( -19.057  -19.164   -0.243)   27.027  10.349   ( -22.854  -29.290    0.011)   37.151  11.245   ( -12.679  -27.609    0.185)   30.382  11.412   ( -14.305  -13.314   -0.049)   19.542  11.836   (  -4.175  -19.269   -0.661)   19.728  12.732   (   0.634    9.948    1.431)   10.070  14.850   (   2.052   10.397   -0.553)   10.612  15.346   (   6.463   15.813    5.758)   18.027  15.685   (   0.966   17.689   -5.396)   18.519======================= Grid point 383 (78/84) =======================q-point: ( 0.46  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.291   (   1.541    7.866   -0.567)    8.036   2.375   (   2.093    9.233   -0.809)    9.502   2.381   (   0.866    6.345   -0.602)    6.432   2.514   (   2.800    3.458   -0.857)    4.531   2.694   (   6.065    8.950    1.883)   10.974   2.930   (   0.880   14.240    1.857)   14.387   2.973   (   5.140    5.685   -1.770)    7.866   3.054   (  -2.943    8.805   -2.167)    9.533   3.228   (  -1.869   10.783    0.618)   10.962   4.922   (   2.349   20.612   -0.015)   20.745   5.008   (   0.714  -10.731    0.623)   10.772   5.574   (   5.194    3.342    0.017)    6.177   7.909   (  -3.710   -0.937   -1.115)    3.986   8.579   (  17.801  -50.784   -0.268)   53.814   9.150   (  30.787  -41.027    0.237)   51.294   9.400   ( -21.117  -24.036   -0.154)   31.995   9.827   ( -29.552  -32.465    0.093)   43.901  11.098   (  -3.181  -26.826    0.104)   27.014  11.211   (  -5.482  -16.463    0.152)   17.353  11.780   (  -1.498  -21.194   -0.562)   21.254  12.739   (   0.166   10.967    1.034)   11.017  14.881   (   0.904   12.477   -0.829)   12.537  15.429   (   2.077   15.733    5.694)   16.860  15.690   (  -0.056   18.618   -4.718)   19.207======================= Grid point 394 (79/84) =======================q-point: ( 0.31  0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.229   (  13.784   13.784   -0.738)   19.507   2.306   (  13.859   13.859    1.107)   19.631   2.358   (   4.775    4.775   -1.270)    6.871   2.376   (   2.154    2.154   -0.531)    3.092   2.535   (   5.518    5.518    1.532)    7.952   2.817   (   3.363    3.363   -1.424)    4.965   3.031   (   4.923    4.923    6.017)    9.202   3.215   (   4.232    4.232   -6.874)    9.114   3.285   (  13.634   13.634   -0.453)   19.287   4.914   (  -8.002   -8.002    1.081)   11.369   5.105   (  13.886   13.886    0.176)   19.638   5.162   (  12.849   12.849   -0.189)   18.172   7.482   ( -11.089  -11.089    0.143)   15.683   7.943   (  -5.794   -5.794    0.093)    8.194   8.252   ( -16.062  -16.062   -0.275)   22.716   9.886   ( -10.464  -10.464    0.411)   14.804  10.285   ( -22.937  -22.937   -0.028)   32.438  11.060   ( -26.989  -26.989    0.021)   38.169  11.497   ( -13.320  -13.320    0.216)   18.839  11.607   ( -17.027  -17.027   -1.278)   24.114  12.825   (   4.212    4.212    1.822)    6.229  14.934   (   6.275    6.275   -0.471)    8.886  15.437   (  10.860   10.860    4.881)   16.115  15.855   (   7.528    7.528   -5.128)   11.817======================= Grid point 395 (80/84) =======================q-point: ( 0.38  0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.393   (   3.405    7.529   -0.392)    8.272   2.415   (   1.778    3.011   -0.481)    3.529   2.472   (   6.667    2.834   -0.687)    7.276   2.552   (   9.109   15.375    0.078)   17.871   2.701   (  11.441   10.876    1.578)   15.864   2.915   (   6.510    3.278   -1.503)    7.443   3.105   (   2.698    7.438    4.152)    8.935   3.262   (   0.801    7.219   -4.636)    8.616   3.493   (   6.728   15.863   -0.320)   17.234   4.787   (  -4.399  -11.488    0.959)   12.339   5.250   (   5.375   19.088   -0.039)   19.830   5.470   (  13.519    5.239    0.102)   14.499   7.454   (   7.650  -24.777   -0.538)   25.936   7.946   ( -13.954  -16.645   -0.088)   21.720   8.047   (  15.836  -17.595    0.320)   23.675   9.590   ( -21.685   -6.591    0.384)   22.668   9.819   ( -25.577  -24.892    0.095)   35.690  10.634   ( -14.330  -34.025    0.160)   36.919  11.136   ( -19.600  -14.634   -0.336)   24.463  11.388   (  -7.185  -25.899   -0.629)   26.885  12.883   (   1.898    5.644    1.499)    6.141  15.043   (   4.625    9.073   -0.884)   10.222  15.601   (   6.089   10.325    4.089)   12.665  15.948   (   2.542    9.525   -3.838)   10.579======================= Grid point 396 (81/84) =======================q-point: ( 0.46  0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.429   (   0.646    3.944    0.475)    4.025   2.448   (   1.193    6.060   -0.718)    6.218   2.550   (   1.733    0.977   -0.428)    2.035   2.662   (   2.568   15.702   -0.542)   15.920   2.911   (   7.585   13.241    0.829)   15.282   3.033   (   4.170    3.122   -1.141)    5.332   3.130   (  -0.062    8.504    3.294)    9.120   3.258   (  -0.687    9.466   -3.330)   10.058   3.551   (   0.343   18.425   -0.202)   18.429   4.744   (  -0.632  -14.437    0.839)   14.475   5.325   (   2.035   18.360   -0.067)   18.472   5.652   (   4.785    4.359    0.137)    6.474   7.693   (  11.716  -33.578   -1.080)   35.579   7.751   (  -5.001  -18.999   -0.037)   19.647   8.489   (  25.296  -24.515    0.662)   35.232   9.131   ( -18.816   -4.070    0.296)   19.253   9.283   ( -26.856  -23.115   -0.045)   35.434  10.489   (  -2.572  -34.859    0.527)   34.957  10.852   (  -7.937  -19.689   -0.263)   21.230  11.294   (  -2.499  -28.177   -0.667)   28.296  12.905   (   0.514    6.221    1.305)    6.377  15.106   (   1.689   10.310   -1.129)   10.508  15.678   (   1.926    9.919    3.811)   10.800  15.973   (   0.478   10.638   -3.291)   11.146======================= Grid point 408 (82/84) =======================q-point: ( 0.38  0.38  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.458   (   1.871    1.871    0.263)    2.659   2.503   (   3.527    3.527   -0.296)    4.997   2.533   (   3.466    3.466   -0.193)    4.905   2.827   (  11.180   11.180    0.502)   15.818   2.954   (  13.172   13.172    1.301)   18.673   2.998   (   5.993    5.993   -1.638)    8.632   3.230   (   5.155    5.155    3.866)    8.252   3.367   (   3.599    3.599   -4.352)    6.697   3.759   (  10.221   10.221   -0.706)   14.472   4.575   (  -9.471   -9.471    0.637)   13.410   5.556   (   8.787    8.787    0.237)   12.429   5.575   (   8.082    8.082    0.083)   11.430   7.158   (  -5.132   -5.132   -1.387)    7.389   7.620   ( -15.373  -15.373    0.312)   21.743   7.942   (   6.657    6.657    0.582)    9.432   9.366   ( -22.334  -22.334    0.051)   31.585   9.479   (  -9.695   -9.695    1.167)   13.760  10.008   ( -25.163  -25.163    0.229)   35.587  10.829   ( -20.401  -20.401   -0.881)   28.865  10.909   ( -16.956  -16.956   -0.276)   23.982  12.961   (   2.593    2.593    1.493)    3.959  15.193   (   6.052    6.052   -1.222)    8.646  15.754   (   5.530    5.530    3.098)    8.412  16.075   (   4.001    4.001   -2.886)    6.352======================= Grid point 409 (83/84) =======================q-point: ( 0.46  0.38  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.483   (   0.645    2.116    0.511)    2.270   2.556   (   1.703    4.226    0.213)    4.561   2.579   (   1.108    2.266    0.008)    2.522   2.954   (   2.548   13.012    0.125)   13.260   3.097   (   3.967    3.213   -0.214)    5.109   3.198   (   8.621   15.015    0.663)   17.327   3.289   (   0.472    7.666    3.769)    8.555   3.409   (   0.917    5.972   -4.386)    7.466   3.869   (   1.559   13.319   -0.883)   13.438   4.451   (  -2.903  -15.126    0.183)   15.403   5.625   (   1.670   12.006    0.190)   12.124   5.733   (   4.350    3.619    0.201)    5.662   7.190   (   4.713  -15.069   -2.090)   15.926   7.398   (  -6.533  -17.081    0.487)   18.294   8.234   (  18.822   -1.799    0.917)   18.930   8.927   ( -20.958  -13.610   -0.221)   24.990   9.150   ( -15.399    3.537    1.144)   15.841   9.788   (  -2.182  -34.938    1.002)   35.020  10.476   ( -10.658  -17.704   -0.574)   20.673  10.716   (  -4.962  -29.394   -0.647)   29.817  12.992   (   0.714    2.893    1.464)    3.320  15.274   (   2.113    6.658   -1.437)    7.132  15.823   (   1.703    5.160    2.655)    6.048  16.120   (   1.002    4.850   -2.446)    5.524======================= Grid point 422 (84/84) =======================q-point: ( 0.46  0.46  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 0.00e+00 1.56e-04 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.511   (   0.735    0.735    0.734)    1.273   2.613   (   1.716    1.716    0.461)    2.469   2.619   (   1.127    1.127    0.529)    1.680   3.130   (   3.987    3.987    0.584)    5.669   3.147   (   2.237    2.237    0.318)    3.180   3.368   (   2.813    2.813    5.087)    6.458   3.487   (   4.576    4.576   -3.904)    7.557   3.497   (   7.831    7.831   -0.876)   11.110   4.045   (   3.971    3.971   -1.389)    5.785   4.192   (  -8.995   -8.995   -0.730)   12.742   5.783   (   2.897    2.897    0.246)    4.105   5.785   (   2.676    2.676    0.232)    3.792   7.065   (  -0.433   -0.433   -2.936)    2.999   7.151   (  -7.119   -7.119    0.608)   10.087   8.345   (   9.022    9.022    1.117)   12.807   8.723   (  -7.810   -7.810   -0.525)   11.057   9.192   (  -4.168   -4.168    1.818)    6.169   9.260   ( -10.989  -10.989    1.529)   15.616  10.206   (  -9.278   -9.278   -0.708)   13.141  10.257   ( -13.226  -13.226   -0.613)   18.714  13.026   (   0.804    0.804    1.529)    1.906  15.362   (   2.290    2.290   -1.627)    3.624  15.887   (   1.560    1.560    2.106)    3.050  16.178   (   1.314    1.314   -2.057)    2.772=================== 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/20280   10.0   1216.058   1216.058    650.665     -0.000     -0.000      0.000 3/20280   20.0    332.228    332.228     57.641     -0.000      0.000      0.000 3/20280   30.0    206.170    206.170     20.977     -0.000      0.000      0.000 3/20280   40.0    147.438    147.438     12.949     -0.000      0.000      0.000 3/20280   50.0    112.263    112.263      9.496     -0.000      0.000      0.000 3/20280   60.0     88.933     88.933      7.460     -0.000      0.000      0.000 3/20280   70.0     72.397     72.397      6.082     -0.000      0.000      0.000 3/20280   80.0     60.204     60.204      5.085     -0.000      0.000      0.000 3/20280   90.0     50.988     50.988      4.339     -0.000     -0.000      0.000 3/20280  100.0     43.893     43.893      3.766     -0.000     -0.000      0.000 3/20280  110.0     38.340     38.340      3.316     -0.000     -0.000      0.000 3/20280  120.0     33.925     33.925      2.957     -0.000     -0.000      0.000 3/20280  130.0     30.362     30.362      2.666     -0.000     -0.000      0.000 3/20280  140.0     27.444     27.444      2.425     -0.000     -0.000      0.000 3/20280  150.0     25.023     25.023      2.225     -0.000     -0.000      0.000 3/20280  160.0     22.987     22.987      2.055     -0.000     -0.000      0.000 3/20280  170.0     21.257     21.257      1.909     -0.000     -0.000      0.000 3/20280  180.0     19.770     19.770      1.784     -0.000     -0.000      0.000 3/20280  190.0     18.479     18.479      1.674     -0.000     -0.000      0.000 3/20280  200.0     17.350     17.350      1.577     -0.000     -0.000      0.000 3/20280  210.0     16.354     16.354      1.491     -0.000     -0.000      0.000 3/20280  220.0     15.469     15.469      1.415     -0.000     -0.000      0.000 3/20280  230.0     14.678     14.678      1.346     -0.000     -0.000      0.000 3/20280  240.0     13.966     13.966      1.284     -0.000     -0.000      0.000 3/20280  250.0     13.322     13.322      1.228     -0.000     -0.000      0.000 3/20280  260.0     12.737     12.737      1.176     -0.000     -0.000      0.000 3/20280  270.0     12.202     12.202      1.129     -0.000     -0.000      0.000 3/20280  280.0     11.713     11.713      1.085     -0.000     -0.000      0.000 3/20280  290.0     11.262     11.262      1.045     -0.000     -0.000      0.000 3/20280  300.0     10.846     10.846      1.008     -0.000     -0.000      0.000 3/20280  310.0     10.460     10.460      0.974     -0.000     -0.000      0.000 3/20280  320.0     10.102     10.102      0.941     -0.000     -0.000      0.000 3/20280  330.0      9.769      9.769      0.911     -0.000     -0.000      0.000 3/20280  340.0      9.457      9.457      0.883     -0.000     -0.000      0.000 3/20280  350.0      9.165      9.165      0.857     -0.000     -0.000      0.000 3/20280  360.0      8.891      8.891      0.832     -0.000     -0.000      0.000 3/20280  370.0      8.634      8.634      0.809     -0.000     -0.000      0.000 3/20280  380.0      8.391      8.391      0.787     -0.000     -0.000      0.000 3/20280  390.0      8.162      8.162      0.766     -0.000     -0.000      0.000 3/20280  400.0      7.946      7.946      0.746     -0.000     -0.000      0.000 3/20280  410.0      7.741      7.741      0.727     -0.000     -0.000      0.000 3/20280  420.0      7.547      7.547      0.709     -0.000     -0.000      0.000 3/20280  430.0      7.362      7.362      0.692     -0.000     -0.000      0.000 3/20280  440.0      7.186      7.186      0.676     -0.000     -0.000      0.000 3/20280  450.0      7.019      7.019      0.661     -0.000     -0.000      0.000 3/20280  460.0      6.860      6.860      0.646     -0.000     -0.000      0.000 3/20280  470.0      6.707      6.707      0.632     -0.000     -0.000      0.000 3/20280  480.0      6.562      6.562      0.619     -0.000     -0.000      0.000 3/20280  490.0      6.423      6.423      0.606     -0.000     -0.000      0.000 3/20280  500.0      6.289      6.289      0.593     -0.000     -0.000      0.000 3/20280  510.0      6.162      6.162      0.582     -0.000     -0.000      0.000 3/20280  520.0      6.039      6.039      0.570     -0.000     -0.000      0.000 3/20280  530.0      5.921      5.921      0.559     -0.000     -0.000      0.000 3/20280  540.0      5.808      5.808      0.549     -0.000     -0.000      0.000 3/20280  550.0      5.699      5.699      0.539     -0.000     -0.000      0.000 3/20280  560.0      5.595      5.595      0.529     -0.000     -0.000      0.000 3/20280  570.0      5.494      5.494      0.519     -0.000     -0.000      0.000 3/20280  580.0      5.396      5.396      0.510     -0.000     -0.000      0.000 3/20280  590.0      5.303      5.303      0.502     -0.000     -0.000      0.000 3/20280  600.0      5.212      5.212      0.493     -0.000     -0.000      0.000 3/20280  610.0      5.124      5.124      0.485     -0.000     -0.000      0.000 3/20280  620.0      5.040      5.040      0.477     -0.000     -0.000      0.000 3/20280  630.0      4.958      4.958      0.469     -0.000     -0.000      0.000 3/20280  640.0      4.879      4.879      0.462     -0.000     -0.000      0.000 3/20280  650.0      4.802      4.802      0.455     -0.000     -0.000      0.000 3/20280  660.0      4.728      4.728      0.448     -0.000     -0.000      0.000 3/20280  670.0      4.656      4.656      0.441     -0.000     -0.000      0.000 3/20280  680.0      4.586      4.586      0.435     -0.000     -0.000      0.000 3/20280  690.0      4.519      4.519      0.428     -0.000     -0.000      0.000 3/20280  700.0      4.453      4.453      0.422     -0.000     -0.000      0.000 3/20280  710.0      4.389      4.389      0.416     -0.000     -0.000      0.000 3/20280  720.0      4.327      4.327      0.410     -0.000     -0.000      0.000 3/20280  730.0      4.267      4.267      0.405     -0.000     -0.000      0.000 3/20280  740.0      4.209      4.209      0.399     -0.000     -0.000      0.000 3/20280  750.0      4.152      4.152      0.394     -0.000     -0.000      0.000 3/20280  760.0      4.096      4.096      0.389     -0.000     -0.000      0.000 3/20280  770.0      4.042      4.042      0.383     -0.000     -0.000      0.000 3/20280  780.0      3.990      3.990      0.379     -0.000     -0.000      0.000 3/20280  790.0      3.939      3.939      0.374     -0.000     -0.000      0.000 3/20280  800.0      3.889      3.889      0.369     -0.000     -0.000      0.000 3/20280  810.0      3.840      3.840      0.364     -0.000     -0.000      0.000 3/20280  820.0      3.793      3.793      0.360     -0.000     -0.000      0.000 3/20280  830.0      3.746      3.746      0.356     -0.000     -0.000      0.000 3/20280  840.0      3.701      3.701      0.351     -0.000     -0.000      0.000 3/20280  850.0      3.657      3.657      0.347     -0.000     -0.000      0.000 3/20280  860.0      3.614      3.614      0.343     -0.000     -0.000      0.000 3/20280  870.0      3.572      3.572      0.339     -0.000     -0.000      0.000 3/20280  880.0      3.531      3.531      0.335     -0.000     -0.000      0.000 3/20280  890.0      3.491      3.491      0.332     -0.000     -0.000      0.000 3/20280  900.0      3.452      3.452      0.328     -0.000     -0.000      0.000 3/20280  910.0      3.414      3.414      0.324     -0.000     -0.000      0.000 3/20280  920.0      3.376      3.376      0.321     -0.000     -0.000      0.000 3/20280  930.0      3.340      3.340      0.317     -0.000     -0.000      0.000 3/20280  940.0      3.304      3.304      0.314     -0.000     -0.000      0.000 3/20280  950.0      3.269      3.269      0.311     -0.000     -0.000      0.000 3/20280  960.0      3.235      3.235      0.307     -0.000     -0.000      0.000 3/20280  970.0      3.201      3.201      0.304     -0.000     -0.000      0.000 3/20280  980.0      3.168      3.168      0.301     -0.000     -0.000      0.000 3/20280  990.0      3.136      3.136      0.298     -0.000     -0.000      0.000 3/20280 1000.0      3.104      3.104      0.295     -0.000     -0.000      0.000 3/20280Thermal conductivity related properties were written into "kappa-m13135.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 10:51:56]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|