# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/b4905e7e-381e-41ac-a3b0-77418ffec1a7/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 Bi3BrO4 / P4/mmm (123) / materials id 866311](https://mdr.nims.go.jp/datasets/14cc9a89-ed11-4203-89d7-5457d273ae6e)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-09 04:15:03]-------------------------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.933297820000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.933297820000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.182162070000000Atomic positions (fractional):   *1 Bi  0.00000000000000  0.00000000000000  0.27912707883787 208.980    2 Bi  0.00000000000000  0.00000000000000  0.72087292116213 208.980   *3 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980   *4 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904   *5 O   0.00000000000000  0.50000000000000  0.16056373316296  15.999    6 O   0.50000000000000  0.00000000000000  0.83943626683704  15.999    7 O   0.00000000000000  0.50000000000000  0.83943626683704  15.999    8 O   0.50000000000000  0.00000000000000  0.16056373316296  15.999-------------------------------- unit cell ---------------------------------Lattice vectors:  a    3.933297820000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.933297820000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.182162070000000Atomic positions (fractional):   *1 Bi  0.00000000000000  0.00000000000000  0.27912707883787 208.980 > 1    2 Bi  0.00000000000000  0.00000000000000  0.72087292116213 208.980 > 2   *3 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980 > 3   *4 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 4   *5 O   0.00000000000000  0.50000000000000  0.16056373316296  15.999 > 5    6 O   0.50000000000000  0.00000000000000  0.83943626683704  15.999 > 6    7 O   0.00000000000000  0.50000000000000  0.83943626683704  15.999 > 7    8 O   0.50000000000000  0.00000000000000  0.16056373316296  15.999 > 8-------------------------------- super cell --------------------------------Lattice vectors:  a   11.799893460000000    0.000000000000000    0.000000000000000  b    0.000000000000000   11.799893460000000    0.000000000000000  c    0.000000000000000    0.000000000000000   18.364324140000001Atomic positions (fractional):   *1 Bi  0.00000000000000  0.00000000000000  0.13956353941893 208.980 > 1    2 Bi  0.33333333333333  0.00000000000000  0.13956353941893 208.980 > 1    3 Bi  0.66666666666667  0.00000000000000  0.13956353941893 208.980 > 1    4 Bi  0.00000000000000  0.33333333333333  0.13956353941893 208.980 > 1    5 Bi  0.33333333333333  0.33333333333333  0.13956353941893 208.980 > 1    6 Bi  0.66666666666667  0.33333333333333  0.13956353941893 208.980 > 1    7 Bi  0.00000000000000  0.66666666666667  0.13956353941893 208.980 > 1    8 Bi  0.33333333333333  0.66666666666667  0.13956353941893 208.980 > 1    9 Bi  0.66666666666667  0.66666666666667  0.13956353941893 208.980 > 1   10 Bi  0.00000000000000  0.00000000000000  0.63956353941893 208.980 > 1   11 Bi  0.33333333333333  0.00000000000000  0.63956353941893 208.980 > 1   12 Bi  0.66666666666667  0.00000000000000  0.63956353941893 208.980 > 1   13 Bi  0.00000000000000  0.33333333333333  0.63956353941893 208.980 > 1   14 Bi  0.33333333333333  0.33333333333333  0.63956353941893 208.980 > 1   15 Bi  0.66666666666667  0.33333333333333  0.63956353941893 208.980 > 1   16 Bi  0.00000000000000  0.66666666666667  0.63956353941893 208.980 > 1   17 Bi  0.33333333333333  0.66666666666667  0.63956353941893 208.980 > 1   18 Bi  0.66666666666667  0.66666666666667  0.63956353941893 208.980 > 1   19 Bi  0.00000000000000  0.00000000000000  0.36043646058107 208.980 > 2   20 Bi  0.33333333333333  0.00000000000000  0.36043646058107 208.980 > 2   21 Bi  0.66666666666667  0.00000000000000  0.36043646058107 208.980 > 2   22 Bi  0.00000000000000  0.33333333333333  0.36043646058107 208.980 > 2   23 Bi  0.33333333333333  0.33333333333333  0.36043646058107 208.980 > 2   24 Bi  0.66666666666667  0.33333333333333  0.36043646058107 208.980 > 2   25 Bi  0.00000000000000  0.66666666666667  0.36043646058107 208.980 > 2   26 Bi  0.33333333333333  0.66666666666667  0.36043646058107 208.980 > 2   27 Bi  0.66666666666667  0.66666666666667  0.36043646058107 208.980 > 2   28 Bi  0.00000000000000  0.00000000000000  0.86043646058107 208.980 > 2   29 Bi  0.33333333333333  0.00000000000000  0.86043646058107 208.980 > 2   30 Bi  0.66666666666667  0.00000000000000  0.86043646058107 208.980 > 2   31 Bi  0.00000000000000  0.33333333333333  0.86043646058107 208.980 > 2   32 Bi  0.33333333333333  0.33333333333333  0.86043646058107 208.980 > 2   33 Bi  0.66666666666667  0.33333333333333  0.86043646058107 208.980 > 2   34 Bi  0.00000000000000  0.66666666666667  0.86043646058107 208.980 > 2   35 Bi  0.33333333333333  0.66666666666667  0.86043646058107 208.980 > 2   36 Bi  0.66666666666667  0.66666666666667  0.86043646058107 208.980 > 2  *37 Bi  0.16666666666667  0.16666666666667  0.00000000000000 208.980 > 3   38 Bi  0.50000000000000  0.16666666666667  0.00000000000000 208.980 > 3   39 Bi  0.83333333333333  0.16666666666667  0.00000000000000 208.980 > 3   40 Bi  0.16666666666667  0.50000000000000  0.00000000000000 208.980 > 3   41 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980 > 3   42 Bi  0.83333333333333  0.50000000000000  0.00000000000000 208.980 > 3   43 Bi  0.16666666666667  0.83333333333333  0.00000000000000 208.980 > 3   44 Bi  0.50000000000000  0.83333333333333  0.00000000000000 208.980 > 3   45 Bi  0.83333333333333  0.83333333333333  0.00000000000000 208.980 > 3   46 Bi  0.16666666666667  0.16666666666667  0.50000000000000 208.980 > 3   47 Bi  0.50000000000000  0.16666666666667  0.50000000000000 208.980 > 3   48 Bi  0.83333333333333  0.16666666666667  0.50000000000000 208.980 > 3   49 Bi  0.16666666666667  0.50000000000000  0.50000000000000 208.980 > 3   50 Bi  0.50000000000000  0.50000000000000  0.50000000000000 208.980 > 3   51 Bi  0.83333333333333  0.50000000000000  0.50000000000000 208.980 > 3   52 Bi  0.16666666666667  0.83333333333333  0.50000000000000 208.980 > 3   53 Bi  0.50000000000000  0.83333333333333  0.50000000000000 208.980 > 3   54 Bi  0.83333333333333  0.83333333333333  0.50000000000000 208.980 > 3  *55 Br  0.16666666666667  0.16666666666667  0.25000000000000  79.904 > 4   56 Br  0.50000000000000  0.16666666666667  0.25000000000000  79.904 > 4   57 Br  0.83333333333333  0.16666666666667  0.25000000000000  79.904 > 4   58 Br  0.16666666666667  0.50000000000000  0.25000000000000  79.904 > 4   59 Br  0.50000000000000  0.50000000000000  0.25000000000000  79.904 > 4   60 Br  0.83333333333333  0.50000000000000  0.25000000000000  79.904 > 4   61 Br  0.16666666666667  0.83333333333333  0.25000000000000  79.904 > 4   62 Br  0.50000000000000  0.83333333333333  0.25000000000000  79.904 > 4   63 Br  0.83333333333333  0.83333333333333  0.25000000000000  79.904 > 4   64 Br  0.16666666666667  0.16666666666667  0.75000000000000  79.904 > 4   65 Br  0.50000000000000  0.16666666666667  0.75000000000000  79.904 > 4   66 Br  0.83333333333333  0.16666666666667  0.75000000000000  79.904 > 4   67 Br  0.16666666666667  0.50000000000000  0.75000000000000  79.904 > 4   68 Br  0.50000000000000  0.50000000000000  0.75000000000000  79.904 > 4   69 Br  0.83333333333333  0.50000000000000  0.75000000000000  79.904 > 4   70 Br  0.16666666666667  0.83333333333333  0.75000000000000  79.904 > 4   71 Br  0.50000000000000  0.83333333333333  0.75000000000000  79.904 > 4   72 Br  0.83333333333333  0.83333333333333  0.75000000000000  79.904 > 4  *73 O   0.00000000000000  0.16666666666667  0.08028186658148  15.999 > 5   74 O   0.33333333333333  0.16666666666667  0.08028186658148  15.999 > 5   75 O   0.66666666666667  0.16666666666667  0.08028186658148  15.999 > 5   76 O   0.00000000000000  0.50000000000000  0.08028186658148  15.999 > 5   77 O   0.33333333333333  0.50000000000000  0.08028186658148  15.999 > 5   78 O   0.66666666666667  0.50000000000000  0.08028186658148  15.999 > 5   79 O   0.00000000000000  0.83333333333333  0.08028186658148  15.999 > 5   80 O   0.33333333333333  0.83333333333333  0.08028186658148  15.999 > 5   81 O   0.66666666666667  0.83333333333333  0.08028186658148  15.999 > 5   82 O   0.00000000000000  0.16666666666667  0.58028186658148  15.999 > 5   83 O   0.33333333333333  0.16666666666667  0.58028186658148  15.999 > 5   84 O   0.66666666666667  0.16666666666667  0.58028186658148  15.999 > 5   85 O   0.00000000000000  0.50000000000000  0.58028186658148  15.999 > 5   86 O   0.33333333333333  0.50000000000000  0.58028186658148  15.999 > 5   87 O   0.66666666666667  0.50000000000000  0.58028186658148  15.999 > 5   88 O   0.00000000000000  0.83333333333333  0.58028186658148  15.999 > 5   89 O   0.33333333333333  0.83333333333333  0.58028186658148  15.999 > 5   90 O   0.66666666666667  0.83333333333333  0.58028186658148  15.999 > 5   91 O   0.16666666666667  0.00000000000000  0.41971813341852  15.999 > 6   92 O   0.50000000000000  0.00000000000000  0.41971813341852  15.999 > 6   93 O   0.83333333333333  0.00000000000000  0.41971813341852  15.999 > 6   94 O   0.16666666666667  0.33333333333333  0.41971813341852  15.999 > 6   95 O   0.50000000000000  0.33333333333333  0.41971813341852  15.999 > 6   96 O   0.83333333333333  0.33333333333333  0.41971813341852  15.999 > 6   97 O   0.16666666666667  0.66666666666667  0.41971813341852  15.999 > 6   98 O   0.50000000000000  0.66666666666667  0.41971813341852  15.999 > 6   99 O   0.83333333333333  0.66666666666667  0.41971813341852  15.999 > 6  100 O   0.16666666666667  0.00000000000000  0.91971813341852  15.999 > 6  101 O   0.50000000000000  0.00000000000000  0.91971813341852  15.999 > 6  102 O   0.83333333333333  0.00000000000000  0.91971813341852  15.999 > 6  103 O   0.16666666666667  0.33333333333333  0.91971813341852  15.999 > 6  104 O   0.50000000000000  0.33333333333333  0.91971813341852  15.999 > 6  105 O   0.83333333333333  0.33333333333333  0.91971813341852  15.999 > 6  106 O   0.16666666666667  0.66666666666667  0.91971813341852  15.999 > 6  107 O   0.50000000000000  0.66666666666667  0.91971813341852  15.999 > 6  108 O   0.83333333333333  0.66666666666667  0.91971813341852  15.999 > 6  109 O   0.00000000000000  0.16666666666667  0.41971813341852  15.999 > 7  110 O   0.33333333333333  0.16666666666667  0.41971813341852  15.999 > 7  111 O   0.66666666666667  0.16666666666667  0.41971813341852  15.999 > 7  112 O   0.00000000000000  0.50000000000000  0.41971813341852  15.999 > 7  113 O   0.33333333333333  0.50000000000000  0.41971813341852  15.999 > 7  114 O   0.66666666666667  0.50000000000000  0.41971813341852  15.999 > 7  115 O   0.00000000000000  0.83333333333333  0.41971813341852  15.999 > 7  116 O   0.33333333333333  0.83333333333333  0.41971813341852  15.999 > 7  117 O   0.66666666666667  0.83333333333333  0.41971813341852  15.999 > 7  118 O   0.00000000000000  0.16666666666667  0.91971813341852  15.999 > 7  119 O   0.33333333333333  0.16666666666667  0.91971813341852  15.999 > 7  120 O   0.66666666666667  0.16666666666667  0.91971813341852  15.999 > 7  121 O   0.00000000000000  0.50000000000000  0.91971813341852  15.999 > 7  122 O   0.33333333333333  0.50000000000000  0.91971813341852  15.999 > 7  123 O   0.66666666666667  0.50000000000000  0.91971813341852  15.999 > 7  124 O   0.00000000000000  0.83333333333333  0.91971813341852  15.999 > 7  125 O   0.33333333333333  0.83333333333333  0.91971813341852  15.999 > 7  126 O   0.66666666666667  0.83333333333333  0.91971813341852  15.999 > 7  127 O   0.16666666666667  0.00000000000000  0.08028186658148  15.999 > 8  128 O   0.50000000000000  0.00000000000000  0.08028186658148  15.999 > 8  129 O   0.83333333333333  0.00000000000000  0.08028186658148  15.999 > 8  130 O   0.16666666666667  0.33333333333333  0.08028186658148  15.999 > 8  131 O   0.50000000000000  0.33333333333333  0.08028186658148  15.999 > 8  132 O   0.83333333333333  0.33333333333333  0.08028186658148  15.999 > 8  133 O   0.16666666666667  0.66666666666667  0.08028186658148  15.999 > 8  134 O   0.50000000000000  0.66666666666667  0.08028186658148  15.999 > 8  135 O   0.83333333333333  0.66666666666667  0.08028186658148  15.999 > 8  136 O   0.16666666666667  0.00000000000000  0.58028186658148  15.999 > 8  137 O   0.50000000000000  0.00000000000000  0.58028186658148  15.999 > 8  138 O   0.83333333333333  0.00000000000000  0.58028186658148  15.999 > 8  139 O   0.16666666666667  0.33333333333333  0.58028186658148  15.999 > 8  140 O   0.50000000000000  0.33333333333333  0.58028186658148  15.999 > 8  141 O   0.83333333333333  0.33333333333333  0.58028186658148  15.999 > 8  142 O   0.16666666666667  0.66666666666667  0.58028186658148  15.999 > 8  143 O   0.50000000000000  0.66666666666667  0.58028186658148  15.999 > 8  144 O   0.83333333333333  0.66666666666667  0.58028186658148  15.999 > 8----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            7.3288954    0.0000000    0.0000000            0.0000000    7.3288954    0.0000000            0.0000000    0.0000000    6.1816784-------------------------- Born effective charges --------------------------    1 Bi    5.2106629    0.0000000    0.0000000            0.0000000    5.2106629    0.0000000            0.0000000    0.0000000    3.0859888    2 Bi    5.2106629    0.0000000    0.0000000            0.0000000    5.2106629    0.0000000            0.0000000    0.0000000    3.0859888    3 Bi    4.6668054    0.0000000    0.0000000            0.0000000    4.6668054    0.0000000            0.0000000    0.0000000    5.6419354    4 Br   -2.0787426    0.0000000    0.0000000            0.0000000   -2.0787426    0.0000000            0.0000000    0.0000000   -1.9263779    5 O    -2.6653376    0.0000000    0.0000000            0.0000000   -3.8393567    0.0000000            0.0000000    0.0000000   -2.4718838    6 O    -3.8393567    0.0000000    0.0000000            0.0000000   -2.6653376    0.0000000            0.0000000    0.0000000   -2.4718838    7 O    -2.6653376    0.0000000    0.0000000            0.0000000   -3.8393567    0.0000000            0.0000000    0.0000000   -2.4718838    8 O    -3.8393567    0.0000000    0.0000000            0.0000000   -2.6653376    0.0000000            0.0000000    0.0000000   -2.4718838----------------------------------------------------------------------------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.004Solver_block: 80 / 80 - Time: 0.364Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.371--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) 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 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-09 04:15:08]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-09 04:15:09]-------------------------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.933297820000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.933297820000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.182162070000000Atomic positions (fractional):    1 Bi  0.00000000000000  0.00000000000000  0.27912707883787 208.980    2 Bi  0.00000000000000  0.00000000000000  0.72087292116213 208.980    3 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980    4 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904    5 O   0.00000000000000  0.50000000000000  0.16056373316296  15.999    6 O   0.50000000000000  0.00000000000000  0.83943626683704  15.999    7 O   0.00000000000000  0.50000000000000  0.83943626683704  15.999    8 O   0.50000000000000  0.00000000000000  0.16056373316296  15.999-------------------------------- supercell ---------------------------------Lattice vectors:  a   11.799893460000000    0.000000000000000    0.000000000000000  b    0.000000000000000   11.799893460000000    0.000000000000000  c    0.000000000000000    0.000000000000000   18.364324140000001Atomic positions (fractional):    1 Bi  0.00000000000000  0.00000000000000  0.13956353941893 208.980 > 1    2 Bi  0.33333333333333  0.00000000000000  0.13956353941893 208.980 > 1    3 Bi  0.66666666666667  0.00000000000000  0.13956353941893 208.980 > 1    4 Bi  0.00000000000000  0.33333333333333  0.13956353941893 208.980 > 1    5 Bi  0.33333333333333  0.33333333333333  0.13956353941893 208.980 > 1    6 Bi  0.66666666666667  0.33333333333333  0.13956353941893 208.980 > 1    7 Bi  0.00000000000000  0.66666666666667  0.13956353941893 208.980 > 1    8 Bi  0.33333333333333  0.66666666666667  0.13956353941893 208.980 > 1    9 Bi  0.66666666666667  0.66666666666667  0.13956353941893 208.980 > 1   10 Bi  0.00000000000000  0.00000000000000  0.63956353941893 208.980 > 1   11 Bi  0.33333333333333  0.00000000000000  0.63956353941893 208.980 > 1   12 Bi  0.66666666666667  0.00000000000000  0.63956353941893 208.980 > 1   13 Bi  0.00000000000000  0.33333333333333  0.63956353941893 208.980 > 1   14 Bi  0.33333333333333  0.33333333333333  0.63956353941893 208.980 > 1   15 Bi  0.66666666666667  0.33333333333333  0.63956353941893 208.980 > 1   16 Bi  0.00000000000000  0.66666666666667  0.63956353941893 208.980 > 1   17 Bi  0.33333333333333  0.66666666666667  0.63956353941893 208.980 > 1   18 Bi  0.66666666666667  0.66666666666667  0.63956353941893 208.980 > 1   19 Bi  0.00000000000000  0.00000000000000  0.36043646058107 208.980 > 19   20 Bi  0.33333333333333  0.00000000000000  0.36043646058107 208.980 > 19   21 Bi  0.66666666666667  0.00000000000000  0.36043646058107 208.980 > 19   22 Bi  0.00000000000000  0.33333333333333  0.36043646058107 208.980 > 19   23 Bi  0.33333333333333  0.33333333333333  0.36043646058107 208.980 > 19   24 Bi  0.66666666666667  0.33333333333333  0.36043646058107 208.980 > 19   25 Bi  0.00000000000000  0.66666666666667  0.36043646058107 208.980 > 19   26 Bi  0.33333333333333  0.66666666666667  0.36043646058107 208.980 > 19   27 Bi  0.66666666666667  0.66666666666667  0.36043646058107 208.980 > 19   28 Bi  0.00000000000000  0.00000000000000  0.86043646058107 208.980 > 19   29 Bi  0.33333333333333  0.00000000000000  0.86043646058107 208.980 > 19   30 Bi  0.66666666666667  0.00000000000000  0.86043646058107 208.980 > 19   31 Bi  0.00000000000000  0.33333333333333  0.86043646058107 208.980 > 19   32 Bi  0.33333333333333  0.33333333333333  0.86043646058107 208.980 > 19   33 Bi  0.66666666666667  0.33333333333333  0.86043646058107 208.980 > 19   34 Bi  0.00000000000000  0.66666666666667  0.86043646058107 208.980 > 19   35 Bi  0.33333333333333  0.66666666666667  0.86043646058107 208.980 > 19   36 Bi  0.66666666666667  0.66666666666667  0.86043646058107 208.980 > 19   37 Bi  0.16666666666667  0.16666666666667  0.00000000000000 208.980 > 37   38 Bi  0.50000000000000  0.16666666666667  0.00000000000000 208.980 > 37   39 Bi  0.83333333333333  0.16666666666667  0.00000000000000 208.980 > 37   40 Bi  0.16666666666667  0.50000000000000  0.00000000000000 208.980 > 37   41 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980 > 37   42 Bi  0.83333333333333  0.50000000000000  0.00000000000000 208.980 > 37   43 Bi  0.16666666666667  0.83333333333333  0.00000000000000 208.980 > 37   44 Bi  0.50000000000000  0.83333333333333  0.00000000000000 208.980 > 37   45 Bi  0.83333333333333  0.83333333333333  0.00000000000000 208.980 > 37   46 Bi  0.16666666666667  0.16666666666667  0.50000000000000 208.980 > 37   47 Bi  0.50000000000000  0.16666666666667  0.50000000000000 208.980 > 37   48 Bi  0.83333333333333  0.16666666666667  0.50000000000000 208.980 > 37   49 Bi  0.16666666666667  0.50000000000000  0.50000000000000 208.980 > 37   50 Bi  0.50000000000000  0.50000000000000  0.50000000000000 208.980 > 37   51 Bi  0.83333333333333  0.50000000000000  0.50000000000000 208.980 > 37   52 Bi  0.16666666666667  0.83333333333333  0.50000000000000 208.980 > 37   53 Bi  0.50000000000000  0.83333333333333  0.50000000000000 208.980 > 37   54 Bi  0.83333333333333  0.83333333333333  0.50000000000000 208.980 > 37   55 Br  0.16666666666667  0.16666666666667  0.25000000000000  79.904 > 55   56 Br  0.50000000000000  0.16666666666667  0.25000000000000  79.904 > 55   57 Br  0.83333333333333  0.16666666666667  0.25000000000000  79.904 > 55   58 Br  0.16666666666667  0.50000000000000  0.25000000000000  79.904 > 55   59 Br  0.50000000000000  0.50000000000000  0.25000000000000  79.904 > 55   60 Br  0.83333333333333  0.50000000000000  0.25000000000000  79.904 > 55   61 Br  0.16666666666667  0.83333333333333  0.25000000000000  79.904 > 55   62 Br  0.50000000000000  0.83333333333333  0.25000000000000  79.904 > 55   63 Br  0.83333333333333  0.83333333333333  0.25000000000000  79.904 > 55   64 Br  0.16666666666667  0.16666666666667  0.75000000000000  79.904 > 55   65 Br  0.50000000000000  0.16666666666667  0.75000000000000  79.904 > 55   66 Br  0.83333333333333  0.16666666666667  0.75000000000000  79.904 > 55   67 Br  0.16666666666667  0.50000000000000  0.75000000000000  79.904 > 55   68 Br  0.50000000000000  0.50000000000000  0.75000000000000  79.904 > 55   69 Br  0.83333333333333  0.50000000000000  0.75000000000000  79.904 > 55   70 Br  0.16666666666667  0.83333333333333  0.75000000000000  79.904 > 55   71 Br  0.50000000000000  0.83333333333333  0.75000000000000  79.904 > 55   72 Br  0.83333333333333  0.83333333333333  0.75000000000000  79.904 > 55   73 O   0.00000000000000  0.16666666666667  0.08028186658148  15.999 > 73   74 O   0.33333333333333  0.16666666666667  0.08028186658148  15.999 > 73   75 O   0.66666666666667  0.16666666666667  0.08028186658148  15.999 > 73   76 O   0.00000000000000  0.50000000000000  0.08028186658148  15.999 > 73   77 O   0.33333333333333  0.50000000000000  0.08028186658148  15.999 > 73   78 O   0.66666666666667  0.50000000000000  0.08028186658148  15.999 > 73   79 O   0.00000000000000  0.83333333333333  0.08028186658148  15.999 > 73   80 O   0.33333333333333  0.83333333333333  0.08028186658148  15.999 > 73   81 O   0.66666666666667  0.83333333333333  0.08028186658148  15.999 > 73   82 O   0.00000000000000  0.16666666666667  0.58028186658148  15.999 > 73   83 O   0.33333333333333  0.16666666666667  0.58028186658148  15.999 > 73   84 O   0.66666666666667  0.16666666666667  0.58028186658148  15.999 > 73   85 O   0.00000000000000  0.50000000000000  0.58028186658148  15.999 > 73   86 O   0.33333333333333  0.50000000000000  0.58028186658148  15.999 > 73   87 O   0.66666666666667  0.50000000000000  0.58028186658148  15.999 > 73   88 O   0.00000000000000  0.83333333333333  0.58028186658148  15.999 > 73   89 O   0.33333333333333  0.83333333333333  0.58028186658148  15.999 > 73   90 O   0.66666666666667  0.83333333333333  0.58028186658148  15.999 > 73   91 O   0.16666666666667  0.00000000000000  0.41971813341852  15.999 > 91   92 O   0.50000000000000  0.00000000000000  0.41971813341852  15.999 > 91   93 O   0.83333333333333  0.00000000000000  0.41971813341852  15.999 > 91   94 O   0.16666666666667  0.33333333333333  0.41971813341852  15.999 > 91   95 O   0.50000000000000  0.33333333333333  0.41971813341852  15.999 > 91   96 O   0.83333333333333  0.33333333333333  0.41971813341852  15.999 > 91   97 O   0.16666666666667  0.66666666666667  0.41971813341852  15.999 > 91   98 O   0.50000000000000  0.66666666666667  0.41971813341852  15.999 > 91   99 O   0.83333333333333  0.66666666666667  0.41971813341852  15.999 > 91  100 O   0.16666666666667  0.00000000000000  0.91971813341852  15.999 > 91  101 O   0.50000000000000  0.00000000000000  0.91971813341852  15.999 > 91  102 O   0.83333333333333  0.00000000000000  0.91971813341852  15.999 > 91  103 O   0.16666666666667  0.33333333333333  0.91971813341852  15.999 > 91  104 O   0.50000000000000  0.33333333333333  0.91971813341852  15.999 > 91  105 O   0.83333333333333  0.33333333333333  0.91971813341852  15.999 > 91  106 O   0.16666666666667  0.66666666666667  0.91971813341852  15.999 > 91  107 O   0.50000000000000  0.66666666666667  0.91971813341852  15.999 > 91  108 O   0.83333333333333  0.66666666666667  0.91971813341852  15.999 > 91  109 O   0.00000000000000  0.16666666666667  0.41971813341852  15.999 > 109  110 O   0.33333333333333  0.16666666666667  0.41971813341852  15.999 > 109  111 O   0.66666666666667  0.16666666666667  0.41971813341852  15.999 > 109  112 O   0.00000000000000  0.50000000000000  0.41971813341852  15.999 > 109  113 O   0.33333333333333  0.50000000000000  0.41971813341852  15.999 > 109  114 O   0.66666666666667  0.50000000000000  0.41971813341852  15.999 > 109  115 O   0.00000000000000  0.83333333333333  0.41971813341852  15.999 > 109  116 O   0.33333333333333  0.83333333333333  0.41971813341852  15.999 > 109  117 O   0.66666666666667  0.83333333333333  0.41971813341852  15.999 > 109  118 O   0.00000000000000  0.16666666666667  0.91971813341852  15.999 > 109  119 O   0.33333333333333  0.16666666666667  0.91971813341852  15.999 > 109  120 O   0.66666666666667  0.16666666666667  0.91971813341852  15.999 > 109  121 O   0.00000000000000  0.50000000000000  0.91971813341852  15.999 > 109  122 O   0.33333333333333  0.50000000000000  0.91971813341852  15.999 > 109  123 O   0.66666666666667  0.50000000000000  0.91971813341852  15.999 > 109  124 O   0.00000000000000  0.83333333333333  0.91971813341852  15.999 > 109  125 O   0.33333333333333  0.83333333333333  0.91971813341852  15.999 > 109  126 O   0.66666666666667  0.83333333333333  0.91971813341852  15.999 > 109  127 O   0.16666666666667  0.00000000000000  0.08028186658148  15.999 > 127  128 O   0.50000000000000  0.00000000000000  0.08028186658148  15.999 > 127  129 O   0.83333333333333  0.00000000000000  0.08028186658148  15.999 > 127  130 O   0.16666666666667  0.33333333333333  0.08028186658148  15.999 > 127  131 O   0.50000000000000  0.33333333333333  0.08028186658148  15.999 > 127  132 O   0.83333333333333  0.33333333333333  0.08028186658148  15.999 > 127  133 O   0.16666666666667  0.66666666666667  0.08028186658148  15.999 > 127  134 O   0.50000000000000  0.66666666666667  0.08028186658148  15.999 > 127  135 O   0.83333333333333  0.66666666666667  0.08028186658148  15.999 > 127  136 O   0.16666666666667  0.00000000000000  0.58028186658148  15.999 > 127  137 O   0.50000000000000  0.00000000000000  0.58028186658148  15.999 > 127  138 O   0.83333333333333  0.00000000000000  0.58028186658148  15.999 > 127  139 O   0.16666666666667  0.33333333333333  0.58028186658148  15.999 > 127  140 O   0.50000000000000  0.33333333333333  0.58028186658148  15.999 > 127  141 O   0.83333333333333  0.33333333333333  0.58028186658148  15.999 > 127  142 O   0.16666666666667  0.66666666666667  0.58028186658148  15.999 > 127  143 O   0.50000000000000  0.66666666666667  0.58028186658148  15.999 > 127  144 O   0.83333333333333  0.66666666666667  0.58028186658148  15.999 > 127----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            7.3288954    0.0000000    0.0000000            0.0000000    7.3288954    0.0000000            0.0000000    0.0000000    6.1816784-------------------------- Born effective charges --------------------------    1 Bi    5.2106629    0.0000000    0.0000000            0.0000000    5.2106629    0.0000000            0.0000000    0.0000000    3.0859888    2 Bi    5.2106629    0.0000000    0.0000000            0.0000000    5.2106629    0.0000000            0.0000000    0.0000000    3.0859888    3 Bi    4.6668054    0.0000000    0.0000000            0.0000000    4.6668054    0.0000000            0.0000000    0.0000000    5.6419354    4 Br   -2.0787426    0.0000000    0.0000000            0.0000000   -2.0787426    0.0000000            0.0000000    0.0000000   -1.9263779    5 O    -2.6653376    0.0000000    0.0000000            0.0000000   -3.8393567    0.0000000            0.0000000    0.0000000   -2.4718838    6 O    -3.8393567    0.0000000    0.0000000            0.0000000   -2.6653376    0.0000000            0.0000000    0.0000000   -2.4718838    7 O    -2.6653376    0.0000000    0.0000000            0.0000000   -3.8393567    0.0000000            0.0000000    0.0000000   -2.4718838    8 O    -3.8393567    0.0000000    0.0000000            0.0000000   -2.6653376    0.0000000            0.0000000    0.0000000   -2.4718838----------------------------------------------------------------------------Sets of supercell forces were read from "FORCES_FC3.xz".Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".----------------------------- Force constants ------------------------------Computing fc3[ 1, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 37, 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[ 55, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 73, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0100  0.0000]    [ 0.0000 -0.0100  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Expanding fc3.Symmetrizing fc3 by traditional approach (N=3).Symmetrizing fc2 by traditional approach (N=3).Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz)fc3 was written into "fc3.hdf5".Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330----------- None of ph-ph interaction calculation was performed. -----------Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-09 04:15:18]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-09 04:15:18]-------------------------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.933297820000000    0.000000000000000    0.000000000000000  b    0.000000000000000    3.933297820000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.182162070000000Atomic positions (fractional):    1 Bi  0.00000000000000  0.00000000000000  0.27912707883787 208.980    2 Bi  0.00000000000000  0.00000000000000  0.72087292116213 208.980    3 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980    4 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904    5 O   0.00000000000000  0.50000000000000  0.16056373316296  15.999    6 O   0.50000000000000  0.00000000000000  0.83943626683704  15.999    7 O   0.00000000000000  0.50000000000000  0.83943626683704  15.999    8 O   0.50000000000000  0.00000000000000  0.16056373316296  15.999-------------------------------- supercell ---------------------------------Lattice vectors:  a   11.799893460000000    0.000000000000000    0.000000000000000  b    0.000000000000000   11.799893460000000    0.000000000000000  c    0.000000000000000    0.000000000000000   18.364324140000001Atomic positions (fractional):    1 Bi  0.00000000000000  0.00000000000000  0.13956353941893 208.980 > 1    2 Bi  0.33333333333333  0.00000000000000  0.13956353941893 208.980 > 1    3 Bi  0.66666666666667  0.00000000000000  0.13956353941893 208.980 > 1    4 Bi  0.00000000000000  0.33333333333333  0.13956353941893 208.980 > 1    5 Bi  0.33333333333333  0.33333333333333  0.13956353941893 208.980 > 1    6 Bi  0.66666666666667  0.33333333333333  0.13956353941893 208.980 > 1    7 Bi  0.00000000000000  0.66666666666667  0.13956353941893 208.980 > 1    8 Bi  0.33333333333333  0.66666666666667  0.13956353941893 208.980 > 1    9 Bi  0.66666666666667  0.66666666666667  0.13956353941893 208.980 > 1   10 Bi  0.00000000000000  0.00000000000000  0.63956353941893 208.980 > 1   11 Bi  0.33333333333333  0.00000000000000  0.63956353941893 208.980 > 1   12 Bi  0.66666666666667  0.00000000000000  0.63956353941893 208.980 > 1   13 Bi  0.00000000000000  0.33333333333333  0.63956353941893 208.980 > 1   14 Bi  0.33333333333333  0.33333333333333  0.63956353941893 208.980 > 1   15 Bi  0.66666666666667  0.33333333333333  0.63956353941893 208.980 > 1   16 Bi  0.00000000000000  0.66666666666667  0.63956353941893 208.980 > 1   17 Bi  0.33333333333333  0.66666666666667  0.63956353941893 208.980 > 1   18 Bi  0.66666666666667  0.66666666666667  0.63956353941893 208.980 > 1   19 Bi  0.00000000000000  0.00000000000000  0.36043646058107 208.980 > 19   20 Bi  0.33333333333333  0.00000000000000  0.36043646058107 208.980 > 19   21 Bi  0.66666666666667  0.00000000000000  0.36043646058107 208.980 > 19   22 Bi  0.00000000000000  0.33333333333333  0.36043646058107 208.980 > 19   23 Bi  0.33333333333333  0.33333333333333  0.36043646058107 208.980 > 19   24 Bi  0.66666666666667  0.33333333333333  0.36043646058107 208.980 > 19   25 Bi  0.00000000000000  0.66666666666667  0.36043646058107 208.980 > 19   26 Bi  0.33333333333333  0.66666666666667  0.36043646058107 208.980 > 19   27 Bi  0.66666666666667  0.66666666666667  0.36043646058107 208.980 > 19   28 Bi  0.00000000000000  0.00000000000000  0.86043646058107 208.980 > 19   29 Bi  0.33333333333333  0.00000000000000  0.86043646058107 208.980 > 19   30 Bi  0.66666666666667  0.00000000000000  0.86043646058107 208.980 > 19   31 Bi  0.00000000000000  0.33333333333333  0.86043646058107 208.980 > 19   32 Bi  0.33333333333333  0.33333333333333  0.86043646058107 208.980 > 19   33 Bi  0.66666666666667  0.33333333333333  0.86043646058107 208.980 > 19   34 Bi  0.00000000000000  0.66666666666667  0.86043646058107 208.980 > 19   35 Bi  0.33333333333333  0.66666666666667  0.86043646058107 208.980 > 19   36 Bi  0.66666666666667  0.66666666666667  0.86043646058107 208.980 > 19   37 Bi  0.16666666666667  0.16666666666667  0.00000000000000 208.980 > 37   38 Bi  0.50000000000000  0.16666666666667  0.00000000000000 208.980 > 37   39 Bi  0.83333333333333  0.16666666666667  0.00000000000000 208.980 > 37   40 Bi  0.16666666666667  0.50000000000000  0.00000000000000 208.980 > 37   41 Bi  0.50000000000000  0.50000000000000  0.00000000000000 208.980 > 37   42 Bi  0.83333333333333  0.50000000000000  0.00000000000000 208.980 > 37   43 Bi  0.16666666666667  0.83333333333333  0.00000000000000 208.980 > 37   44 Bi  0.50000000000000  0.83333333333333  0.00000000000000 208.980 > 37   45 Bi  0.83333333333333  0.83333333333333  0.00000000000000 208.980 > 37   46 Bi  0.16666666666667  0.16666666666667  0.50000000000000 208.980 > 37   47 Bi  0.50000000000000  0.16666666666667  0.50000000000000 208.980 > 37   48 Bi  0.83333333333333  0.16666666666667  0.50000000000000 208.980 > 37   49 Bi  0.16666666666667  0.50000000000000  0.50000000000000 208.980 > 37   50 Bi  0.50000000000000  0.50000000000000  0.50000000000000 208.980 > 37   51 Bi  0.83333333333333  0.50000000000000  0.50000000000000 208.980 > 37   52 Bi  0.16666666666667  0.83333333333333  0.50000000000000 208.980 > 37   53 Bi  0.50000000000000  0.83333333333333  0.50000000000000 208.980 > 37   54 Bi  0.83333333333333  0.83333333333333  0.50000000000000 208.980 > 37   55 Br  0.16666666666667  0.16666666666667  0.25000000000000  79.904 > 55   56 Br  0.50000000000000  0.16666666666667  0.25000000000000  79.904 > 55   57 Br  0.83333333333333  0.16666666666667  0.25000000000000  79.904 > 55   58 Br  0.16666666666667  0.50000000000000  0.25000000000000  79.904 > 55   59 Br  0.50000000000000  0.50000000000000  0.25000000000000  79.904 > 55   60 Br  0.83333333333333  0.50000000000000  0.25000000000000  79.904 > 55   61 Br  0.16666666666667  0.83333333333333  0.25000000000000  79.904 > 55   62 Br  0.50000000000000  0.83333333333333  0.25000000000000  79.904 > 55   63 Br  0.83333333333333  0.83333333333333  0.25000000000000  79.904 > 55   64 Br  0.16666666666667  0.16666666666667  0.75000000000000  79.904 > 55   65 Br  0.50000000000000  0.16666666666667  0.75000000000000  79.904 > 55   66 Br  0.83333333333333  0.16666666666667  0.75000000000000  79.904 > 55   67 Br  0.16666666666667  0.50000000000000  0.75000000000000  79.904 > 55   68 Br  0.50000000000000  0.50000000000000  0.75000000000000  79.904 > 55   69 Br  0.83333333333333  0.50000000000000  0.75000000000000  79.904 > 55   70 Br  0.16666666666667  0.83333333333333  0.75000000000000  79.904 > 55   71 Br  0.50000000000000  0.83333333333333  0.75000000000000  79.904 > 55   72 Br  0.83333333333333  0.83333333333333  0.75000000000000  79.904 > 55   73 O   0.00000000000000  0.16666666666667  0.08028186658148  15.999 > 73   74 O   0.33333333333333  0.16666666666667  0.08028186658148  15.999 > 73   75 O   0.66666666666667  0.16666666666667  0.08028186658148  15.999 > 73   76 O   0.00000000000000  0.50000000000000  0.08028186658148  15.999 > 73   77 O   0.33333333333333  0.50000000000000  0.08028186658148  15.999 > 73   78 O   0.66666666666667  0.50000000000000  0.08028186658148  15.999 > 73   79 O   0.00000000000000  0.83333333333333  0.08028186658148  15.999 > 73   80 O   0.33333333333333  0.83333333333333  0.08028186658148  15.999 > 73   81 O   0.66666666666667  0.83333333333333  0.08028186658148  15.999 > 73   82 O   0.00000000000000  0.16666666666667  0.58028186658148  15.999 > 73   83 O   0.33333333333333  0.16666666666667  0.58028186658148  15.999 > 73   84 O   0.66666666666667  0.16666666666667  0.58028186658148  15.999 > 73   85 O   0.00000000000000  0.50000000000000  0.58028186658148  15.999 > 73   86 O   0.33333333333333  0.50000000000000  0.58028186658148  15.999 > 73   87 O   0.66666666666667  0.50000000000000  0.58028186658148  15.999 > 73   88 O   0.00000000000000  0.83333333333333  0.58028186658148  15.999 > 73   89 O   0.33333333333333  0.83333333333333  0.58028186658148  15.999 > 73   90 O   0.66666666666667  0.83333333333333  0.58028186658148  15.999 > 73   91 O   0.16666666666667  0.00000000000000  0.41971813341852  15.999 > 91   92 O   0.50000000000000  0.00000000000000  0.41971813341852  15.999 > 91   93 O   0.83333333333333  0.00000000000000  0.41971813341852  15.999 > 91   94 O   0.16666666666667  0.33333333333333  0.41971813341852  15.999 > 91   95 O   0.50000000000000  0.33333333333333  0.41971813341852  15.999 > 91   96 O   0.83333333333333  0.33333333333333  0.41971813341852  15.999 > 91   97 O   0.16666666666667  0.66666666666667  0.41971813341852  15.999 > 91   98 O   0.50000000000000  0.66666666666667  0.41971813341852  15.999 > 91   99 O   0.83333333333333  0.66666666666667  0.41971813341852  15.999 > 91  100 O   0.16666666666667  0.00000000000000  0.91971813341852  15.999 > 91  101 O   0.50000000000000  0.00000000000000  0.91971813341852  15.999 > 91  102 O   0.83333333333333  0.00000000000000  0.91971813341852  15.999 > 91  103 O   0.16666666666667  0.33333333333333  0.91971813341852  15.999 > 91  104 O   0.50000000000000  0.33333333333333  0.91971813341852  15.999 > 91  105 O   0.83333333333333  0.33333333333333  0.91971813341852  15.999 > 91  106 O   0.16666666666667  0.66666666666667  0.91971813341852  15.999 > 91  107 O   0.50000000000000  0.66666666666667  0.91971813341852  15.999 > 91  108 O   0.83333333333333  0.66666666666667  0.91971813341852  15.999 > 91  109 O   0.00000000000000  0.16666666666667  0.41971813341852  15.999 > 109  110 O   0.33333333333333  0.16666666666667  0.41971813341852  15.999 > 109  111 O   0.66666666666667  0.16666666666667  0.41971813341852  15.999 > 109  112 O   0.00000000000000  0.50000000000000  0.41971813341852  15.999 > 109  113 O   0.33333333333333  0.50000000000000  0.41971813341852  15.999 > 109  114 O   0.66666666666667  0.50000000000000  0.41971813341852  15.999 > 109  115 O   0.00000000000000  0.83333333333333  0.41971813341852  15.999 > 109  116 O   0.33333333333333  0.83333333333333  0.41971813341852  15.999 > 109  117 O   0.66666666666667  0.83333333333333  0.41971813341852  15.999 > 109  118 O   0.00000000000000  0.16666666666667  0.91971813341852  15.999 > 109  119 O   0.33333333333333  0.16666666666667  0.91971813341852  15.999 > 109  120 O   0.66666666666667  0.16666666666667  0.91971813341852  15.999 > 109  121 O   0.00000000000000  0.50000000000000  0.91971813341852  15.999 > 109  122 O   0.33333333333333  0.50000000000000  0.91971813341852  15.999 > 109  123 O   0.66666666666667  0.50000000000000  0.91971813341852  15.999 > 109  124 O   0.00000000000000  0.83333333333333  0.91971813341852  15.999 > 109  125 O   0.33333333333333  0.83333333333333  0.91971813341852  15.999 > 109  126 O   0.66666666666667  0.83333333333333  0.91971813341852  15.999 > 109  127 O   0.16666666666667  0.00000000000000  0.08028186658148  15.999 > 127  128 O   0.50000000000000  0.00000000000000  0.08028186658148  15.999 > 127  129 O   0.83333333333333  0.00000000000000  0.08028186658148  15.999 > 127  130 O   0.16666666666667  0.33333333333333  0.08028186658148  15.999 > 127  131 O   0.50000000000000  0.33333333333333  0.08028186658148  15.999 > 127  132 O   0.83333333333333  0.33333333333333  0.08028186658148  15.999 > 127  133 O   0.16666666666667  0.66666666666667  0.08028186658148  15.999 > 127  134 O   0.50000000000000  0.66666666666667  0.08028186658148  15.999 > 127  135 O   0.83333333333333  0.66666666666667  0.08028186658148  15.999 > 127  136 O   0.16666666666667  0.00000000000000  0.58028186658148  15.999 > 127  137 O   0.50000000000000  0.00000000000000  0.58028186658148  15.999 > 127  138 O   0.83333333333333  0.00000000000000  0.58028186658148  15.999 > 127  139 O   0.16666666666667  0.33333333333333  0.58028186658148  15.999 > 127  140 O   0.50000000000000  0.33333333333333  0.58028186658148  15.999 > 127  141 O   0.83333333333333  0.33333333333333  0.58028186658148  15.999 > 127  142 O   0.16666666666667  0.66666666666667  0.58028186658148  15.999 > 127  143 O   0.50000000000000  0.66666666666667  0.58028186658148  15.999 > 127  144 O   0.83333333333333  0.66666666666667  0.58028186658148  15.999 > 127----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            7.3288954    0.0000000    0.0000000            0.0000000    7.3288954    0.0000000            0.0000000    0.0000000    6.1816784-------------------------- Born effective charges --------------------------    1 Bi    5.2106629    0.0000000    0.0000000            0.0000000    5.2106629    0.0000000            0.0000000    0.0000000    3.0859888    2 Bi    5.2106629    0.0000000    0.0000000            0.0000000    5.2106629    0.0000000            0.0000000    0.0000000    3.0859888    3 Bi    4.6668054    0.0000000    0.0000000            0.0000000    4.6668054    0.0000000            0.0000000    0.0000000    5.6419354    4 Br   -2.0787426    0.0000000    0.0000000            0.0000000   -2.0787426    0.0000000            0.0000000    0.0000000   -1.9263779    5 O    -2.6653376    0.0000000    0.0000000            0.0000000   -3.8393567    0.0000000            0.0000000    0.0000000   -2.4718838    6 O    -3.8393567    0.0000000    0.0000000            0.0000000   -2.6653376    0.0000000            0.0000000    0.0000000   -2.4718838    7 O    -2.6653376    0.0000000    0.0000000            0.0000000   -3.8393567    0.0000000            0.0000000    0.0000000   -2.4718838    8 O    -3.8393567    0.0000000    0.0000000            0.0000000   -2.6653376    0.0000000            0.0000000    0.0000000   -2.4718838----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz)Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330Length for sampling mesh generation: 50.00Generating grid system ... [ 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.80, Number of G-points: 299, Lambda: 0.22Running 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: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 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.671   (   0.000    0.000    0.000)    0.000   1.671   (   0.000    0.000    0.000)    0.000   2.001   (   0.000    0.000    0.000)    0.000   2.001   (   0.000    0.000    0.000)    0.000   2.205   (   0.000    0.000    0.000)    0.000   2.512   (   0.000    0.000    0.000)    0.000   2.512   (   0.000    0.000    0.000)    0.000   3.382   (   0.000    0.000    0.000)    0.000   5.138   (   0.000    0.000    0.000)    0.000   6.605   (   0.000    0.000    0.000)    0.000   6.605   (   0.000    0.000    0.000)    0.000   8.214   (   0.000    0.000    0.000)    0.000   8.214   (   0.000    0.000    0.000)    0.000   9.280   (   0.000    0.000    0.000)    0.000  11.486   (   0.000    0.000    0.000)    0.000  11.664   (   0.000    0.000    0.000)    0.000  11.946   (   0.000    0.000    0.000)    0.000  11.946   (   0.000    0.000    0.000)    0.000  12.885   (   0.000    0.000    0.000)    0.000  12.885   (   0.000    0.000    0.000)    0.000  14.918   (   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: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.277   (  15.556    0.000    0.000)   15.556   0.524   (  25.638    0.000    0.000)   25.638   0.760   (  37.808    0.000    0.000)   37.808   1.703   (   3.095    0.000    0.000)    3.095   1.797   (  11.745    0.000    0.000)   11.745   2.012   (   1.035    0.000    0.000)    1.035   2.222   (   2.010    0.000    0.000)    2.010   2.427   (   1.578    0.000    0.000)    1.578   2.503   (  -0.915    0.000    0.000)    0.915   2.608   (   8.458    0.000    0.000)    8.458   3.421   (   4.330    0.000    0.000)    4.330   4.996   ( -13.863    0.000    0.000)   13.863   6.714   (  10.885    0.000    0.000)   10.885   8.122   (  -9.041    0.000    0.000)    9.041   8.354   (  14.029    0.000    0.000)   14.029   9.266   (  -1.258    0.000    0.000)    1.258   9.925   (   6.123    0.000    0.000)    6.123  11.546   (   4.416    0.000    0.000)    4.416  11.613   (  -5.156    0.000    0.000)    5.156  11.955   (   0.955    0.000    0.000)    0.955  12.850   (  -3.562    0.000    0.000)    3.562  12.959   (   7.886    0.000    0.000)    7.886  14.809   ( -10.721    0.000    0.000)   10.721  14.842   (  -2.463    0.000    0.000)    2.463======================= Grid point 2 (3/84) =======================q-point: ( 0.15  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.608   (  17.949    0.000    0.000)   17.949   0.987   (  21.477    0.000    0.000)   21.477   1.454   (  32.375    0.000    0.000)   32.375   1.779   (   4.330    0.000    0.000)    4.330   2.039   (   1.710    0.000    0.000)    1.710   2.048   (  10.817    0.000    0.000)   10.817   2.320   (   9.737    0.000    0.000)    9.737   2.440   (  -0.179    0.000    0.000)    0.179   2.477   (  -1.789    0.000    0.000)    1.789   2.843   (  14.737    0.000    0.000)   14.737   3.564   (  10.482    0.000    0.000)   10.482   4.642   ( -20.841    0.000    0.000)   20.841   7.009   (  18.407    0.000    0.000)   18.407   7.878   ( -15.497    0.000    0.000)   15.497   8.742   (  24.926    0.000    0.000)   24.926   9.235   (  -1.909    0.000    0.000)    1.909  10.076   (   8.459    0.000    0.000)    8.459  11.471   (  -8.940    0.000    0.000)    8.940  11.582   (  -1.833    0.000    0.000)    1.833  11.984   (   1.994    0.000    0.000)    1.994  12.746   (  -7.041    0.000    0.000)    7.041  13.199   (  16.062    0.000    0.000)   16.062  14.523   ( -17.491    0.000    0.000)   17.491  14.775   (  -4.270    0.000    0.000)    4.270======================= Grid point 3 (4/84) =======================q-point: ( 0.23  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.962   (  17.784    0.000    0.000)   17.784   1.361   (  16.803    0.000    0.000)   16.803   1.862   (   3.865    0.000    0.000)    3.865   1.985   (  20.721    0.000    0.000)   20.721   2.075   (   1.877    0.000    0.000)    1.877   2.179   (   4.254    0.000    0.000)    4.254   2.429   (  -0.410    0.000    0.000)    0.410   2.434   (  -2.605    0.000    0.000)    2.605   2.594   (  15.302    0.000    0.000)   15.302   3.174   (  19.219    0.000    0.000)   19.219   3.821   (  14.999    0.000    0.000)   14.999   4.249   ( -17.413    0.000    0.000)   17.413   7.400   (  20.665    0.000    0.000)   20.665   7.529   ( -19.789    0.000    0.000)   19.789   9.182   (  -4.015    0.000    0.000)    4.015   9.297   (  30.997    0.000    0.000)   30.997  10.219   (   5.129    0.000    0.000)    5.129  11.282   (  -9.748    0.000    0.000)    9.748  11.481   (  -7.500    0.000    0.000)    7.500  12.033   (   2.966    0.000    0.000)    2.966  12.575   ( -10.450    0.000    0.000)   10.450  13.539   (  17.387    0.000    0.000)   17.387  14.166   ( -17.941    0.000    0.000)   17.941  14.679   (  -5.488    0.000    0.000)    5.488======================= Grid point 4 (5/84) =======================q-point: ( 0.31  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.285   (  14.642    0.000    0.000)   14.642   1.648   (  12.733    0.000    0.000)   12.733   1.926   (   2.720    0.000    0.000)    2.720   2.110   (   1.605    0.000    0.000)    1.605   2.217   (   3.109    0.000    0.000)    3.109   2.251   (   3.252    0.000    0.000)    3.252   2.375   (  -3.369    0.000    0.000)    3.369   2.475   (   5.780    0.000    0.000)    5.780   2.820   (   6.606    0.000    0.000)    6.606   3.582   (  21.518    0.000    0.000)   21.518   4.036   (  -3.119    0.000    0.000)    3.119   4.116   (  14.378    0.000    0.000)   14.378   7.125   ( -20.879    0.000    0.000)   20.879   7.784   (  17.912    0.000    0.000)   17.912   9.064   (  -8.102    0.000    0.000)    8.102   9.923   (  32.318    0.000    0.000)   32.318  10.247   (  -2.502    0.000    0.000)    2.502  11.119   (  -6.362    0.000    0.000)    6.362  11.334   (  -6.592    0.000    0.000)    6.592  12.096   (   3.409    0.000    0.000)    3.409  12.337   ( -13.896    0.000    0.000)   13.896  13.840   (  12.868    0.000    0.000)   12.868  13.859   ( -12.832    0.000    0.000)   12.832  14.564   (  -6.146    0.000    0.000)    6.146======================= Grid point 5 (6/84) =======================q-point: ( 0.38  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.509   (   7.775    0.000    0.000)    7.775   1.863   (   9.315    0.000    0.000)    9.315   1.968   (   1.539    0.000    0.000)    1.539   2.136   (   1.042    0.000    0.000)    1.042   2.210   (  -1.748    0.000    0.000)    1.748   2.302   (  -4.005    0.000    0.000)    4.005   2.305   (   2.255    0.000    0.000)    2.255   2.590   (   4.776    0.000    0.000)    4.776   2.881   (   1.474    0.000    0.000)    1.474   3.966   (  16.615    0.000    0.000)   16.615   4.095   (   6.679    0.000    0.000)    6.679   4.357   (   9.815    0.000    0.000)    9.815   6.745   ( -16.996    0.000    0.000)   16.996   8.078   (  11.800    0.000    0.000)   11.800   8.884   (  -9.340    0.000    0.000)    9.340  10.143   (  -6.806    0.000    0.000)    6.806  10.537   (  29.869    0.000    0.000)   29.869  11.037   (  -2.370    0.000    0.000)    2.370  11.241   (  -2.925    0.000    0.000)    2.925  12.032   ( -17.300    0.000    0.000)   17.300  12.159   (   2.808    0.000    0.000)    2.808  13.676   (  -5.932    0.000    0.000)    5.932  14.035   (   7.061    0.000    0.000)    7.061  14.446   (  -5.569    0.000    0.000)    5.569======================= Grid point 6 (7/84) =======================q-point: ( 0.46  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.594   (   1.597    0.000    0.000)    1.597   1.987   (   0.491    0.000    0.000)    0.491   2.008   (   4.933    0.000    0.000)    4.933   2.149   (   0.275    0.000    0.000)    0.275   2.181   (  -0.865    0.000    0.000)    0.865   2.226   (  -3.129    0.000    0.000)    3.129   2.336   (   0.852    0.000    0.000)    0.852   2.651   (   1.548    0.000    0.000)    1.548   2.913   (   1.392    0.000    0.000)    1.392   4.196   (   6.199    0.000    0.000)    6.199   4.209   (   3.596    0.000    0.000)    3.596   4.488   (   3.396    0.000    0.000)    3.396   6.502   (  -6.820    0.000    0.000)    6.820   8.235   (   4.077    0.000    0.000)    4.077   8.742   (  -4.165    0.000    0.000)    4.165  10.034   (  -3.354    0.000    0.000)    3.354  11.011   (  -0.564    0.000    0.000)    0.564  11.060   (  21.489    0.000    0.000)   21.489  11.209   (  -0.619    0.000    0.000)    0.619  11.678   ( -17.047    0.000    0.000)   17.047  12.199   (   1.097    0.000    0.000)    1.097  13.611   (  -1.303    0.000    0.000)    1.303  14.123   (   2.142    0.000    0.000)    2.142  14.363   (  -2.441    0.000    0.000)    2.441======================= Grid point 14 (8/84) =======================q-point: ( 0.08  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.394   (  11.048   11.048    0.000)   15.624   0.436   (  10.746   10.746    0.000)   15.197   1.191   (  27.532   27.532    0.000)   38.937   1.827   (   6.897    6.897    0.000)    9.754   1.827   (   7.472    7.472    0.000)   10.567   2.029   (   1.697    1.697    0.000)    2.400   2.237   (   1.782    1.782    0.000)    2.521   2.458   (  -1.072   -1.072    0.000)    1.516   2.522   (   3.806    3.806    0.000)    5.383   2.699   (   9.097    9.097    0.000)   12.866   3.451   (   3.377    3.377    0.000)    4.775   4.886   ( -11.033  -11.033    0.000)   15.604   6.648   (   1.888    1.888    0.000)    2.670   8.078   (  -6.492   -6.492    0.000)    9.181   8.395   (   6.885    6.885    0.000)    9.737   9.238   (  -2.677   -2.677    0.000)    3.786  10.002   (   6.541    6.541    0.000)    9.250  11.503   (  -7.493   -7.493    0.000)   10.597  11.605   (   4.328    4.328    0.000)    6.121  11.875   (  -4.692   -4.692    0.000)    6.635  12.767   (  -5.979   -5.979    0.000)    8.455  13.087   (  10.934   10.934    0.000)   15.463  14.708   (  -9.879   -9.879    0.000)   13.972  14.981   (   5.621    5.621    0.000)    7.950======================= Grid point 15 (9/84) =======================q-point: ( 0.15  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.667   (  16.007    6.230    0.000)   17.176   0.825   (  23.286   -7.333    0.000)   24.413   1.704   (  23.673   17.930    0.000)   29.697   1.898   (   4.201   11.076    0.000)   11.846   2.074   (  10.653    2.795    0.000)   11.013   2.077   (   3.284    3.199    0.000)    4.584   2.326   (   9.106    0.662    0.000)    9.130   2.408   (  -2.807   -1.934    0.000)    3.409   2.596   (   2.293    9.963    0.000)   10.224   2.936   (  14.537    9.022    0.000)   17.109   3.573   (   9.462    1.237    0.000)    9.542   4.590   ( -18.133   -5.473    0.000)   18.941   6.838   (  15.639  -12.458    0.000)   19.994   7.865   ( -14.287   -2.118    0.000)   14.443   8.651   (  18.437   -6.807    0.000)   19.653   9.174   (  -3.572   -6.063    0.000)    7.037  10.148   (   7.568    6.334    0.000)    9.869  11.352   (  -8.405  -10.855    0.000)   13.728  11.639   (  -1.916    4.179    0.000)    4.597  11.796   (  -3.252  -14.341    0.000)   14.705  12.632   (  -8.053   -9.547    0.000)   12.490  13.370   (  17.479   15.274    0.000)   23.212  14.449   ( -15.725   -7.483    0.000)   17.414  15.008   (  -1.725   17.176    0.000)   17.263======================= Grid point 16 (10/84) =======================q-point: ( 0.23  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.992   (  16.579    3.337    0.000)   16.912   1.257   (  20.097   -6.862    0.000)   21.236   1.982   (   4.152   11.247    0.000)   11.989   2.044   (   8.825    5.678    0.000)   10.494   2.180   (   7.567    6.940    0.000)   10.268   2.198   (   3.817    2.084    0.000)    4.349   2.377   (   1.047   -1.372    0.000)    1.726   2.591   (  14.931   -0.339    0.000)   14.935   2.599   (  -1.794   12.978    0.000)   13.101   3.257   (  18.302    8.196    0.000)   20.053   3.817   (  14.734   -0.086    0.000)   14.734   4.237   ( -16.049   -1.362    0.000)   16.107   7.194   (  19.533  -18.789    0.000)   27.103   7.543   ( -18.293    0.912    0.000)   18.316   9.009   (   5.956  -16.444    0.000)   17.489   9.172   (  14.537  -10.907    0.000)   18.174  10.261   (   3.033    3.669    0.000)    4.760  11.175   (  -9.094  -10.463    0.000)   13.863  11.535   (  -7.729    3.938    0.000)    8.675  11.752   (  -1.280  -23.488    0.000)   23.523  12.446   ( -11.104  -11.690    0.000)   16.123  13.730   (  18.089   17.371    0.000)   25.079  14.132   ( -15.777   -3.836    0.000)   16.237  14.945   (  -4.235   22.107    0.000)   22.509======================= Grid point 17 (11/84) =======================q-point: ( 0.31  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.294   (  13.722    1.327    0.000)   13.786   1.599   (  14.898   -3.273    0.000)   15.253   2.057   (   3.481   12.258    0.000)   12.742   2.132   (   2.523    3.339    0.000)    4.185   2.242   (  -0.379    2.108    0.000)    2.141   2.263   (   2.968    1.208    0.000)    3.204   2.482   (   5.025    0.775    0.000)    5.085   2.546   (  -1.476   14.511    0.000)   14.586   2.814   (   6.675   -0.556    0.000)    6.699   3.640   (  20.073    5.827    0.000)   20.902   4.038   (  -3.005    0.114    0.000)    3.008   4.111   (  14.467   -0.232    0.000)   14.469   7.168   ( -19.352    4.144    0.000)   19.790   7.566   (  17.792  -21.374    0.000)   27.810   8.921   (  -8.180  -14.157    0.000)   16.350   9.663   (  28.117  -23.540    0.000)   36.671  10.236   (  -5.809   -1.719    0.000)    6.058  11.027   (  -5.478   -9.197    0.000)   10.705  11.379   (  -7.391    2.836    0.000)    7.916  11.740   (  -0.191  -30.550    0.000)   30.551  12.196   ( -14.303  -13.072    0.000)   19.377  13.864   ( -10.991    0.105    0.000)   10.991  14.044   (  13.499   18.530    0.000)   22.925  14.854   (  -4.810   24.707    0.000)   25.171======================= Grid point 18 (12/84) =======================q-point: ( 0.38  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.505   (   7.548    0.090    0.000)    7.548   1.840   (   9.703   -1.115    0.000)    9.767   2.113   (   2.037   13.660    0.000)   13.811   2.166   (   1.016    3.111    0.000)    3.273   2.232   (  -0.036    3.320    0.000)    3.320   2.313   (   2.112    0.768    0.000)    2.247   2.428   (  -5.137    8.929    0.000)   10.301   2.629   (   4.537    4.657    0.000)    6.502   2.877   (   1.562   -0.258    0.000)    1.583   4.000   (  15.620    3.419    0.000)   15.990   4.094   (   6.490   -0.143    0.000)    6.491   4.355   (   9.964   -0.044    0.000)    9.964   6.816   ( -15.746    7.069    0.000)   17.260   7.871   (  12.940  -20.811    0.000)   24.506   8.744   (  -8.979  -13.776    0.000)   16.444  10.067   (  -9.941   -8.326    0.000)   12.967  10.203   (  26.232  -30.635    0.000)   40.331  10.959   (  -1.992   -8.255    0.000)    8.492  11.263   (  -4.558    0.384    0.000)    4.574  11.738   (  -0.042  -35.888    0.000)   35.888  11.896   ( -15.915  -12.782    0.000)   20.413  13.709   (  -5.061    3.034    0.000)    5.901  14.250   (   7.622   19.487    0.000)   20.925  14.766   (  -3.982   26.861    0.000)   27.155======================= Grid point 19 (13/84) =======================q-point: ( 0.46  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.593   (   1.883    0.049    0.000)    1.884   1.974   (   3.677   -1.143    0.000)    3.851   2.134   (   0.342   14.160    0.000)   14.164   2.176   (   0.137    2.577    0.000)    2.581   2.251   (   1.963    5.764    0.000)    6.089   2.331   (  -4.238    8.975    0.000)    9.925   2.342   (   0.814    0.598    0.000)    1.010   2.687   (   1.464    3.630    0.000)    3.914   2.909   (   1.381   -0.256    0.000)    1.404   4.207   (   3.566   -0.358    0.000)    3.584   4.216   (   5.854    2.159    0.000)    6.239   4.488   (   3.484    0.144    0.000)    3.487   6.591   (  -6.304    8.946    0.000)   10.944   8.052   (   5.026  -18.312    0.000)   18.989   8.609   (  -3.939  -13.055    0.000)   13.636   9.909   (  -4.829  -13.460    0.000)   14.300  10.638   (  15.559  -38.177    0.000)   41.226  10.936   (  -0.518   -8.106    0.000)    8.122  11.184   (  -3.751   -2.815    0.000)    4.690  11.624   (  -9.450   -7.416    0.000)   12.012  11.737   (  -0.060  -39.095    0.000)   39.095  13.652   (  -1.200    4.307    0.000)    4.471  14.346   (   2.378   20.118    0.000)   20.258  14.709   (  -1.602   28.466    0.000)   28.511======================= Grid point 28 (14/84) =======================q-point: ( 0.15  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.847   (  11.655   11.655    0.000)   16.482   0.872   (  12.082   12.082    0.000)   17.087   2.032   (  14.324   14.324    0.000)   20.256   2.100   (   5.989    5.989    0.000)    8.469   2.147   (   4.191    4.191    0.000)    5.927   2.205   (  11.269   11.269    0.000)   15.937   2.359   (   4.351    4.351    0.000)    6.153   2.367   (  -2.016   -2.016    0.000)    2.851   2.803   (   8.972    8.972    0.000)   12.689   3.166   (  13.893   13.893    0.000)   19.648   3.641   (   6.698    6.698    0.000)    9.472   4.425   ( -11.109  -11.109    0.000)   15.710   6.691   (  -1.064   -1.064    0.000)    1.504   7.757   (  -9.262   -9.262    0.000)   13.098   8.593   (   2.390    2.390    0.000)    3.381   9.018   (  -9.076   -9.076    0.000)   12.835  10.277   (   5.685    5.685    0.000)    8.040  11.119   ( -11.842  -11.842    0.000)   16.747  11.518   ( -13.687  -13.687    0.000)   19.357  11.674   (  -1.740   -1.740    0.000)    2.460  12.419   ( -11.926  -11.926    0.000)   16.866  13.733   (  20.622   20.622    0.000)   29.164  14.255   ( -11.507  -11.507    0.000)   16.273  15.269   (   8.530    8.530    0.000)   12.064======================= Grid point 29 (15/84) =======================q-point: ( 0.23  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.099   (  13.418    7.463    0.000)   15.354   1.209   (  19.871    3.360    0.000)   20.153   2.187   (   3.369    8.182    0.000)    8.848   2.199   (   4.923    6.299    0.000)    7.995   2.275   (   4.097    2.959    0.000)    5.054   2.313   (   1.784   12.765    0.000)   12.889   2.367   (   4.924    0.375    0.000)    4.938   2.584   (  14.692    0.053    0.000)   14.692   2.891   (   0.211   14.314    0.000)   14.316   3.467   (  16.531   12.336    0.000)   20.626   3.846   (  13.612    4.255    0.000)   14.261   4.188   ( -11.680   -3.777    0.000)   12.275   6.781   (   9.131  -20.799    0.000)   22.715   7.530   ( -13.144   -3.330    0.000)   13.560   8.681   (   0.356  -13.209    0.000)   13.214   8.897   (   2.395  -16.801    0.000)   16.971  10.322   (  -2.408    1.424    0.000)    2.797  10.907   (  -9.407  -15.894    0.000)   18.469  11.274   ( -10.793  -23.663    0.000)   26.008  11.569   (  -8.317   -1.547    0.000)    8.460  12.166   ( -14.016  -16.325    0.000)   21.516  14.034   ( -10.342   -5.395    0.000)   11.664  14.134   (  19.438   22.282    0.000)   29.568  15.357   (   1.233   17.574    0.000)   17.617======================= Grid point 30 (16/84) =======================q-point: ( 0.31  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.347   (  11.339    4.084    0.000)   12.052   1.583   (  17.297    2.429    0.000)   17.467   2.234   (   1.439    6.313    0.000)    6.475   2.288   (   3.495    2.897    0.000)    4.540   2.294   (   0.947    3.359    0.000)    3.489   2.355   (   3.461   14.787    0.000)   15.187   2.500   (   5.603    0.959    0.000)    5.684   2.813   (   7.508    1.329    0.000)    7.625   2.848   (  -3.486   14.625    0.000)   15.035   3.796   (  16.687    9.507    0.000)   19.205   4.036   (  -2.478   -0.424    0.000)    2.514   4.130   (  14.249    2.971    0.000)   14.555   7.010   (  13.667  -33.452    0.000)   36.136   7.261   ( -13.979    4.162    0.000)   14.585   8.550   (  -9.343  -21.828    0.000)   23.743   9.138   (  17.079  -27.513    0.000)   32.383  10.149   ( -14.989   -7.861    0.000)   16.925  10.767   (  -4.950  -17.390    0.000)   18.081  11.105   (  -6.623  -31.235    0.000)   31.929  11.385   (  -9.577   -2.918    0.000)   10.012  11.870   ( -15.996  -19.740    0.000)   25.408  13.873   (  -5.840    1.358    0.000)    5.996  14.467   (  14.317   22.904    0.000)   27.010  15.346   (  -1.793   22.837    0.000)   22.908======================= Grid point 31 (17/84) =======================q-point: ( 0.38  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.522   (   6.337    1.637    0.000)    6.545   1.861   (  10.850    3.706    0.000)   11.465   2.244   (  -0.030    4.471    0.000)    4.471   2.334   (   1.721    1.281    0.000)    2.145   2.335   (   2.800    6.683    0.000)    7.246   2.440   (   4.110   17.883    0.000)   18.349   2.551   (  -0.558    3.649    0.000)    3.691   2.795   (  -1.312   10.966    0.000)   11.044   2.888   (   1.598    2.416    0.000)    2.897   4.086   (   5.906   -0.782    0.000)    5.958   4.096   (  13.267    6.219    0.000)   14.653   4.371   (   9.895    2.255    0.000)   10.149   7.004   ( -11.609   11.357    0.000)   16.240   7.298   (  15.195  -36.347    0.000)   39.396   8.376   (  -7.872  -22.893    0.000)   24.208   9.480   (  16.702  -40.076    0.000)   43.417   9.789   ( -19.893  -20.572    0.000)   28.617  10.699   (  -2.286  -18.771    0.000)   18.910  11.004   (  -3.943  -35.134    0.000)   35.354  11.214   (  -7.730   -5.593    0.000)    9.542  11.564   ( -14.316  -20.965    0.000)   25.387  13.802   (  -1.808    6.693    0.000)    6.932  14.688   (   8.251   23.253    0.000)   24.673  15.304   (  -2.089   25.598    0.000)   25.683======================= Grid point 32 (18/84) =======================q-point: ( 0.46  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.598   (   1.745    0.550    0.000)    1.830   2.003   (   3.643    4.260    0.000)    5.606   2.243   (  -0.030    4.049    0.000)    4.049   2.359   (   0.725    1.028    0.000)    1.258   2.390   (   2.336    7.861    0.000)    8.200   2.490   (   1.000   20.473    0.000)   20.498   2.507   (  -2.595    7.941    0.000)    8.354   2.790   (   0.233    6.698    0.000)    6.702   2.915   (   1.117    1.629    0.000)    1.975   4.192   (   3.438   -1.265    0.000)    3.663   4.283   (   5.151    4.720    0.000)    6.987   4.503   (   3.446    1.743    0.000)    3.862   6.837   (  -4.687   15.595    0.000)   16.284   7.551   (   8.402  -31.674    0.000)   32.769   8.264   (  -3.091  -21.186    0.000)   21.410   9.463   ( -10.597  -32.846    0.000)   34.513   9.737   (   7.916  -50.677    0.000)   51.291  10.672   (  -0.629  -19.390    0.000)   19.401  10.950   (  -1.509  -37.046    0.000)   37.077  11.087   (  -4.485   -7.318    0.000)    8.583  11.362   (  -5.317  -19.267    0.000)   19.988  13.786   (  -0.195    9.376    0.000)    9.378  14.793   (   2.634   23.472    0.000)   23.619  15.274   (  -0.866   26.850    0.000)   26.863======================= Grid point 42 (19/84) =======================q-point: ( 0.23  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.269   (   9.459    9.459    0.000)   13.378   1.389   (  14.077   14.077    0.000)   19.908   2.280   (   3.313    3.313    0.000)    4.685   2.317   (   1.840    1.840    0.000)    2.602   2.326   (   4.497    4.497    0.000)    6.359   2.406   (   4.541    4.541    0.000)    6.423   2.526   (   3.874    3.874    0.000)    5.479   2.678   (  12.690   12.690    0.000)   17.946   3.075   (   3.714    3.714    0.000)    5.252   3.722   (  13.728   13.728    0.000)   19.415   4.006   (  11.467   11.467    0.000)   16.216   4.094   (  -5.304   -5.304    0.000)    7.501   6.451   ( -11.839  -11.839    0.000)   16.743   7.407   (  -7.980   -7.980    0.000)   11.286   8.521   ( -15.937  -15.937    0.000)   22.538   8.583   (  -2.441   -2.441    0.000)    3.452  10.273   (  -8.057   -8.057    0.000)   11.394  10.584   ( -16.445  -16.445    0.000)   23.256  10.858   ( -18.095  -18.095    0.000)   25.591  11.467   (  -8.513   -8.513    0.000)   12.039  11.818   ( -19.068  -19.068    0.000)   26.966  13.949   (  -2.774   -2.774    0.000)    3.923  14.554   (  19.853   19.853    0.000)   28.076  15.606   (   8.287    8.287    0.000)   11.720======================= Grid point 43 (20/84) =======================q-point: ( 0.31  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.446   (   8.022    5.689    0.000)    9.835   1.704   (  16.625    9.824    0.000)   19.311   2.331   (   1.869    1.935    0.000)    2.690   2.340   (   0.242    4.312    0.000)    4.319   2.389   (   3.543    5.823    0.000)    6.816   2.519   (   5.493    1.049    0.000)    5.592   2.593   (   4.755    7.668    0.000)    9.023   2.926   (  10.166   11.203    0.000)   15.128   3.068   (  -3.117    7.234    0.000)    7.877   3.989   (  13.135    9.896    0.000)   16.445   4.019   (  -1.246   -1.283    0.000)    1.788   4.249   (  12.326    8.976    0.000)   15.248   6.349   (   1.994  -31.382    0.000)   31.445   7.275   (  -4.755   -3.571    0.000)    5.947   8.161   ( -17.298  -16.805    0.000)   24.117   8.661   (   6.025  -19.266    0.000)   20.187   9.922   ( -27.455  -14.802    0.000)   31.191  10.331   (  -8.971  -27.773    0.000)   29.186  10.587   ( -10.125  -21.176    0.000)   23.472  11.262   ( -11.652   -9.402    0.000)   14.972  11.433   ( -19.749  -24.797    0.000)   31.701  13.932   (   0.919    4.757    0.000)    4.845  14.890   (  14.308   19.609    0.000)   24.275  15.704   (   2.383   13.745    0.000)   13.950======================= Grid point 44 (21/84) =======================q-point: ( 0.38  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.567   (   4.308    2.910    0.000)    5.199   1.987   (  11.379    9.004    0.000)   14.510   2.334   (  -0.566    4.650    0.000)    4.684   2.360   (   1.329    1.368    0.000)    1.908   2.469   (   4.279    6.417    0.000)    7.713   2.599   (   2.547    1.501    0.000)    2.957   2.711   (   6.578    8.343    0.000)   10.624   2.992   (  -3.914    8.548    0.000)    9.401   3.047   (   2.642   14.922    0.000)   15.154   4.062   (   4.969   -1.603    0.000)    5.221   4.228   (  11.053    6.704    0.000)   12.927   4.458   (   8.412    6.547    0.000)   10.659   6.532   (  16.022  -38.959    0.000)   42.125   7.212   (  -2.429    8.079    0.000)    8.436   7.897   (  -8.599  -24.480    0.000)   25.946   8.756   (   2.721  -31.011    0.000)   31.130   9.291   ( -34.799  -28.204    0.000)   44.793  10.206   (  -3.997  -31.622    0.000)   31.874  10.429   (  -6.536  -24.558    0.000)   25.413  11.040   ( -10.439  -14.991    0.000)   18.268  11.082   ( -14.886  -25.292    0.000)   29.348  13.969   (   2.359   10.099    0.000)   10.371  15.110   (   8.229   19.282    0.000)   20.965  15.723   (   0.009   16.985    0.000)   16.985======================= Grid point 45 (22/84) =======================q-point: ( 0.46  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.619   (   1.184    1.695    0.000)    2.067   2.135   (   3.731    8.867    0.000)    9.620   2.324   (  -0.334    4.034    0.000)    4.048   2.381   (   0.629    1.196    0.000)    1.351   2.537   (   2.198    6.839    0.000)    7.183   2.622   (   0.157    3.555    0.000)    3.559   2.813   (   3.026    8.968    0.000)    9.465   2.936   (  -1.519    8.052    0.000)    8.194   3.065   (   0.068   16.796    0.000)   16.796   4.156   (   3.188   -2.325    0.000)    3.946   4.395   (   5.031    6.360    0.000)    8.109   4.565   (   2.590    4.623    0.000)    5.299   6.883   (  14.871  -34.138    0.000)   37.236   7.166   (  -1.845   16.369    0.000)   16.473   7.826   (  -0.193  -22.109    0.000)   22.110   8.687   ( -21.504  -42.851    0.000)   47.944   8.766   (  -0.790  -45.936    0.000)   45.943  10.168   (  -0.355  -31.990    0.000)   31.992  10.330   (  -3.162  -27.922    0.000)   28.100  10.874   (  -4.758  -30.488    0.000)   30.857  10.895   (  -4.785  -12.190    0.000)   13.096  14.007   (   1.135   12.664    0.000)   12.715  15.216   (   2.637   19.101    0.000)   19.283  15.719   (  -0.224   18.420    0.000)   18.422======================= Grid point 56 (23/84) =======================q-point: ( 0.31  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.554   (   4.973    4.973    0.000)    7.033   1.946   (  13.901   13.901    0.000)   19.658   2.362   (   1.306    1.306    0.000)    1.846   2.418   (   2.955    2.955    0.000)    4.180   2.471   (   2.685    2.685    0.000)    3.798   2.561   (   3.568    3.568    0.000)    5.046   2.710   (   6.080    6.080    0.000)    8.598   3.142   (   1.105    1.105    0.000)    1.563   3.185   (  13.056   13.056    0.000)   18.463   3.998   (  -0.321   -0.321    0.000)    0.454   4.170   (   8.391    8.391    0.000)   11.867   4.444   (   9.897    9.897    0.000)   13.996   5.897   ( -12.182  -12.182    0.000)   17.228   7.179   (  -3.129   -3.129    0.000)    4.425   7.829   ( -18.438  -18.438    0.000)   26.075   8.447   (  -3.863   -3.863    0.000)    5.464   9.573   ( -27.007  -27.007    0.000)   38.194   9.729   ( -25.798  -25.798    0.000)   36.483  10.247   ( -14.749  -14.749    0.000)   20.859  10.916   ( -26.918  -26.918    0.000)   38.068  11.041   ( -12.713  -12.713    0.000)   17.979  14.049   (   6.777    6.777    0.000)    9.584  15.219   (  13.810   13.810    0.000)   19.531  15.896   (   6.371    6.371    0.000)    9.010======================= Grid point 57 (24/84) =======================q-point: ( 0.38  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.631   (   2.861    3.424    0.000)    4.462   2.194   (  10.073   11.175    0.000)   15.045   2.386   (   1.172    1.238    0.000)    1.705   2.421   (  -1.451    4.144    0.000)    4.391   2.552   (   4.984    0.937    0.000)    5.071   2.629   (   2.784    2.661    0.000)    3.851   2.841   (   6.703    6.713    0.000)    9.487   3.136   (  -1.196    6.759    0.000)    6.864   3.383   (   6.625   17.206    0.000)   18.437   4.028   (   3.606   -1.594    0.000)    3.943   4.337   (   9.414    4.209    0.000)   10.312   4.607   (   6.239    7.923    0.000)   10.084   5.918   (  13.683  -20.177    0.000)   24.379   7.239   (   9.980   -5.068    0.000)   11.193   7.483   ( -15.277  -17.425    0.000)   23.173   8.367   (  -7.807   -7.941    0.000)   11.136   8.776   ( -46.586  -22.955    0.000)   51.935   9.478   (  -6.022  -42.184    0.000)   42.611   9.985   ( -12.081  -21.723    0.000)   24.857  10.436   ( -20.169  -35.162    0.000)   40.536  10.797   ( -11.537  -14.134    0.000)   18.245  14.176   (   5.694   10.516    0.000)   11.959  15.430   (   7.859   13.308    0.000)   15.455  15.975   (   2.245    9.197    0.000)    9.467======================= Grid point 58 (25/84) =======================q-point: ( 0.46  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.667   (   0.893    3.163    0.000)    3.287   2.320   (   3.249    8.597    0.000)    9.190   2.392   (  -1.336    2.766    0.000)    3.072   2.404   (   0.533    1.096    0.000)    1.218   2.634   (   2.847    1.016    0.000)    3.023   2.661   (   0.439    3.184    0.000)    3.214   2.942   (   2.930    6.732    0.000)    7.342   3.114   (  -0.719   10.419    0.000)   10.444   3.453   (   1.306   19.897    0.000)   19.939   4.105   (   2.854   -2.683    0.000)    3.918   4.515   (   6.744    5.640    0.000)    8.792   4.674   (   0.782    6.018    0.000)    6.069   6.330   (  22.006  -19.430    0.000)   29.356   7.262   (  -5.786  -13.416    0.000)   14.610   7.564   (  20.082   -2.036    0.000)   20.185   7.923   ( -33.531  -32.656    0.000)   46.805   8.085   ( -19.385  -17.405    0.000)   26.052   9.433   (   0.204  -42.770    0.000)   42.771   9.794   (  -6.285  -27.711    0.000)   28.415  10.181   (  -5.836  -39.825    0.000)   40.250  10.624   (  -5.116  -14.953    0.000)   15.804  14.256   (   2.182   12.236    0.000)   12.429  15.530   (   2.509   13.029    0.000)   13.268  15.999   (   0.477   10.593    0.000)   10.604======================= Grid point 70 (26/84) =======================q-point: ( 0.38  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.694   (   2.829    2.829    0.000)    4.000   2.357   (   4.099    4.099    0.000)    5.797   2.408   (   0.995    0.995    0.000)    1.407   2.469   (  -0.921   -0.921    0.000)    1.303   2.575   (   4.646    4.646    0.000)    6.570   2.688   (   2.679    2.679    0.000)    3.789   2.983   (   7.149    7.149    0.000)   10.111   3.261   (   5.543    5.543    0.000)    7.839   3.670   (  11.164   11.164    0.000)   15.788   4.015   (   0.772    0.772    0.000)    1.091   4.436   (   7.854    7.854    0.000)   11.107   4.739   (   4.920    4.920    0.000)    6.958   5.813   (   8.985    8.985    0.000)   12.706   7.128   (  -9.480   -9.480    0.000)   13.406   7.240   (  -0.090   -0.090    0.000)    0.127   8.321   ( -33.747  -33.747    0.000)   47.725   8.321   (  -2.329   -2.329    0.000)    3.293   8.696   ( -26.018  -26.018    0.000)   36.795   9.572   ( -20.100  -20.100    0.000)   28.426   9.759   ( -31.165  -31.165    0.000)   44.074  10.526   ( -12.704  -12.704    0.000)   17.966  14.359   (   7.707    7.707    0.000)   10.900  15.632   (   7.495    7.495    0.000)   10.600  16.100   (   4.060    4.060    0.000)    5.742======================= Grid point 71 (27/84) =======================q-point: ( 0.46  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 255Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.734   (   1.112    3.343    0.000)    3.523   2.389   (   0.683   -0.419    0.000)    0.801   2.422   (   0.393    0.763    0.000)    0.858   2.429   (  -2.001   -0.902    0.000)    2.195   2.667   (   2.819    4.305    0.000)    5.146   2.726   (   1.107    2.581    0.000)    2.808   3.089   (   3.023    7.599    0.000)    8.178   3.335   (   1.890   11.319    0.000)   11.476   3.797   (   2.373   14.597    0.000)   14.789   4.059   (   2.410   -1.897    0.000)    3.067   4.630   (   9.027    6.761    0.000)   11.278   4.780   (  -0.876    4.295    0.000)    4.383   6.199   (  25.438    6.493    0.000)   26.254   6.950   (  -5.367  -14.312    0.000)   15.285   7.457   (  16.373   -2.209    0.000)   16.521   7.464   ( -40.195  -14.297    0.000)   42.662   8.125   ( -16.355   11.148    0.000)   19.793   8.527   (  -0.319  -48.403    0.000)   48.404   9.250   ( -10.519  -26.892    0.000)   28.876   9.352   (  -9.066  -43.495    0.000)   44.430  10.338   (  -5.522  -13.468    0.000)   14.556  14.464   (   2.802    8.603    0.000)    9.048  15.728   (   2.384    7.297    0.000)    7.677  16.148   (   1.089    5.047    0.000)    5.163======================= Grid point 84 (28/84) =======================q-point: ( 0.46  0.46  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 147Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.784   (   1.424    1.424    0.000)    2.013   2.369   (  -0.658   -0.658    0.000)    0.930   2.385   (  -2.121   -2.121    0.000)    3.000   2.433   (   0.271    0.271    0.000)    0.383   2.746   (   2.363    2.363    0.000)    3.342   2.760   (   0.941    0.941    0.000)    1.331   3.203   (   3.287    3.287    0.000)    4.649   3.506   (   4.987    4.987    0.000)    7.053   3.996   (   4.664    4.664    0.000)    6.595   4.039   (   0.341    0.341    0.000)    0.482   4.766   (   5.870    5.870    0.000)    8.302   4.827   (   0.119    0.119    0.000)    0.168   6.557   (  27.381   27.381    0.000)   38.722   6.752   (  -5.362   -5.362    0.000)    7.583   7.250   ( -16.496  -16.496    0.000)   23.329   7.605   (  14.439   14.439    0.000)   20.420   7.689   ( -26.994  -26.994    0.000)   38.176   8.276   (  -0.290   -0.290    0.000)    0.410   8.602   ( -25.465  -25.465    0.000)   36.013   8.824   ( -13.536  -13.536    0.000)   19.143  10.138   (  -5.843   -5.843    0.000)    8.263  14.580   (   3.071    3.071    0.000)    4.343  15.821   (   2.317    2.317    0.000)    3.277  16.209   (   1.440    1.440    0.000)    2.037======================= Grid point 169 (29/84) =======================q-point: ( 0.00  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 84Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.345   (  -0.000    0.000   14.609)   14.609   0.345   (   0.000    0.000   14.609)   14.609   0.618   (   0.000    0.000   26.395)   26.395   1.604   (  -0.000   -0.000   -5.416)    5.416   1.604   (  -0.000    0.000   -5.416)    5.416   2.021   (   0.000    0.000    1.343)    1.343   2.021   (   0.000    0.000    1.343)    1.343   2.514   (  -0.000   -0.000    0.167)    0.167   2.514   (   0.000    0.000    0.167)    0.167   2.929   (  -0.000   -0.000  -16.274)   16.274   3.573   (   0.000    0.000    9.099)    9.099   5.158   (  -0.000   -0.000   -0.011)    0.011   6.615   (   0.000    0.000    0.829)    0.829   6.615   (  -0.000    0.000    0.829)    0.829   8.202   (  -0.000    0.000   -0.986)    0.986   8.202   (   0.000    0.000   -0.986)    0.986   9.286   (   0.000    0.000    0.511)    0.511  11.662   (  -0.000   -0.000   -0.202)    0.202  11.955   (  -0.000    0.000    0.700)    0.700  11.955   (   0.000    0.000    0.700)    0.700  12.880   (   0.000    0.000   -0.381)    0.381  12.880   (   0.000    0.000   -0.381)    0.381  14.774   (  -0.000   -0.000    1.206)    1.206  14.912   (   0.000    0.000   -0.508)    0.508======================= Grid point 170 (30/84) =======================q-point: ( 0.08  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.498   (   6.641    0.000   14.700)   16.130   0.622   (  21.056    0.000    7.673)   22.411   0.882   (  26.963    0.000    8.643)   28.315   1.641   (   3.563    0.000   -5.021)    6.156   1.757   (  14.055    0.000   -3.575)   14.503   2.031   (   0.982    0.000    1.310)    1.637   2.111   (   6.070    0.000   -1.151)    6.178   2.505   (  -0.962    0.000    0.130)    0.971   2.585   ( -10.276    0.000    0.625)   10.295   2.623   (   7.498    0.000    1.525)    7.651   3.461   (   0.426    0.000    4.101)    4.124   5.009   ( -13.640    0.000    0.899)   13.670   6.725   (  10.881    0.000    0.826)   10.912   8.126   (  -6.941    0.000    0.130)    6.942   8.342   (  14.052    0.000   -0.970)   14.085   8.485   (  86.591    0.000  -60.965)  105.899   9.275   (  -0.910    0.000    0.682)    1.137  11.596   (  -5.860    0.000   -0.105)    5.861  11.963   (   0.825    0.000    0.596)    1.018  12.100   (   7.193    0.000   17.364)   18.795  12.847   (  -3.445    0.000   -0.287)    3.457  12.910   (   6.003    0.000   -3.040)    6.729  14.608   (  -9.453    0.000   -9.184)   13.180  14.802   ( -10.754    0.000   -0.535)   10.767======================= Grid point 171 (31/84) =======================q-point: ( 0.15  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.689   (  13.476    0.000    6.713)   15.055   1.033   (  19.916    0.000    3.859)   20.287   1.442   (  27.304    0.000   -1.733)   27.359   1.729   (   5.104    0.000   -4.002)    6.486   2.057   (   1.589    0.000    1.206)    1.995   2.069   (  14.548    0.000    1.181)   14.596   2.221   (   6.716    0.000   -3.578)    7.609   2.477   (  -1.876    0.000    0.025)    1.876   2.515   (   2.611    0.000    2.686)    3.746   2.859   (  13.918    0.000    1.021)   13.956   3.554   (   8.849    0.000   -0.686)    8.876   4.664   ( -20.248    0.000    1.699)   20.319   7.019   (  18.388    0.000    0.809)   18.406   7.909   ( -14.705    0.000    2.452)   14.908   8.731   (  24.991    0.000   -0.904)   25.008   9.248   (  -1.119    0.000    0.821)    1.388   9.556   (  33.006    0.000  -31.946)   45.934  11.445   (  -9.334    0.000   -0.807)    9.369  11.988   (   1.770    0.000    0.314)    1.798  12.038   ( -12.463    0.000   27.742)   30.413  12.746   (  -6.848    0.000   -0.037)    6.848  13.181   (  19.485    0.000   -1.981)   19.586  14.501   (  -7.887    0.000  -13.100)   15.291  14.528   ( -10.875    0.000   -6.619)   12.731======================= Grid point 172 (32/84) =======================q-point: ( 0.23  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.994   (  16.737    0.000    2.671)   16.948   1.385   (  15.934    0.000    2.018)   16.061   1.828   (   4.768    0.000   -2.666)    5.463   1.885   (  17.034    0.000   -6.857)   18.362   2.090   (   1.691    0.000    1.027)    1.979   2.230   (   4.612    0.000    2.895)    5.445   2.376   (   4.889    0.000   -3.456)    5.987   2.432   (  -2.712    0.000   -0.131)    2.715   2.683   (  14.047    0.000    6.361)   15.420   3.161   (  17.269    0.000   -1.174)   17.309   3.794   (  14.775    0.000   -2.120)   14.926   4.280   ( -17.109    0.000    2.401)   17.277   7.409   (  20.623    0.000    0.762)   20.637   7.563   ( -20.274    0.000    2.654)   20.447   9.204   (  -3.812    0.000    1.792)    4.212   9.288   (  31.118    0.000   -0.759)   31.127  10.014   (  15.519    0.000  -14.294)   21.099  11.252   (  -9.616    0.000   -1.457)    9.725  11.718   ( -17.973    0.000   16.360)   24.304  12.032   (   2.704    0.000   -0.070)    2.705  12.579   ( -10.244    0.000    0.280)   10.248  13.567   (  18.120    0.000    2.368)   18.274  14.165   ( -17.254    0.000   -0.127)   17.255  14.506   (  -0.250    0.000  -14.034)   14.036======================= Grid point 173 (33/84) =======================q-point: ( 0.31  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.312   (  15.075    0.000    2.209)   15.236   1.658   (  12.139    0.000    0.840)   12.168   1.910   (   3.572    0.000   -1.281)    3.795   2.112   (   7.046    0.000   -6.106)    9.323   2.120   (   1.407    0.000    0.795)    1.616   2.290   (   1.424    0.000    1.136)    1.822   2.371   (  -3.461    0.000   -0.297)    3.474   2.432   (   2.309    0.000   -2.765)    3.602   2.926   (   8.210    0.000    7.752)   11.292   3.524   (  18.749    0.000   -4.455)   19.271   4.073   (   0.616    0.000    2.188)    2.273   4.095   (  12.061    0.000   -1.248)   12.126   7.140   ( -22.151    0.000    1.192)   22.183   7.792   (  17.853    0.000    0.681)   17.866   9.079   (  -9.152    0.000    1.170)    9.226   9.917   (  32.485    0.000   -0.531)   32.489  10.181   (   1.742    0.000   -5.100)    5.389  11.104   (  -4.874    0.000   -0.893)    4.955  11.405   ( -12.793    0.000    5.357)   13.869  12.090   (   3.143    0.000   -0.460)    3.176  12.345   ( -13.709    0.000    0.566)   13.721  13.784   (  -3.151    0.000   -2.388)    3.953  13.951   (   2.876    0.000    5.423)    6.139  14.486   (  -1.781    0.000   -6.441)    6.683======================= Grid point 174 (34/84) =======================q-point: ( 0.38  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.551   (   8.659    0.000    3.362)    9.289   1.863   (   8.854    0.000    0.003)    8.854   1.966   (   2.186    0.000   -0.103)    2.189   2.143   (   0.904    0.000    0.569)    1.068   2.186   (   0.997    0.000   -1.847)    2.099   2.298   (  -4.025    0.000   -0.403)    4.045   2.308   (   1.011    0.000    0.052)    1.013   2.476   (   2.121    0.000   -6.991)    7.305   3.007   (   1.814    0.000    8.067)    8.269   3.841   (  12.881    0.000   -7.453)   14.882   4.181   (   9.755    0.000    4.233)   10.634   4.340   (  10.172    0.000   -1.381)   10.265   6.733   ( -18.323    0.000   -0.937)   18.347   8.086   (  11.746    0.000    0.588)   11.760   8.868   ( -11.176    0.000   -1.227)   11.243  10.124   (  -5.684    0.000   -1.560)    5.894  10.533   (  30.019    0.000   -0.272)   30.020  11.057   (  -0.635    0.000    1.560)    1.685  11.234   (  -5.054    0.000   -0.460)    5.074  12.039   ( -16.971    0.000    0.638)   16.983  12.151   (   2.408    0.000   -0.656)    2.496  13.654   (  -5.092    0.000   -1.666)    5.358  14.071   (   4.275    0.000    2.634)    5.021  14.448   (  -1.635    0.000    0.284)    1.659======================= Grid point 175 (35/84) =======================q-point: ( 0.46  0.00  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.648   (   1.903    0.000    4.380)    4.776   1.987   (   1.816    0.000   -0.020)    1.816   2.008   (   3.560    0.000    0.030)    3.560   2.155   (   0.236    0.000    0.439)    0.498   2.184   (  -0.374    0.000    0.425)    0.567   2.222   (  -3.025    0.000   -0.343)    3.045   2.327   (   0.673    0.000   -0.774)    1.025   2.511   (   1.169    0.000   -8.884)    8.961   3.038   (   1.193    0.000    7.492)    7.587   4.015   (   4.614    0.000   -8.285)    9.483   4.342   (   4.912    0.000    4.347)    6.559   4.476   (   3.545    0.000   -0.960)    3.672   6.471   (  -7.383    0.000   -2.463)    7.783   8.242   (   4.055    0.000    0.526)    4.089   8.697   (  -5.023    0.000   -3.409)    6.071  10.026   (  -3.128    0.000   -0.613)    3.187  11.056   (   0.187    0.000    3.496)    3.501  11.057   (  21.356    0.000   -0.190)   21.357  11.180   (  -1.107    0.000   -2.332)    2.581  11.688   ( -16.822    0.000    0.860)   16.844  12.188   (   1.036    0.000   -0.871)    1.353  13.601   (  -0.958    0.000   -0.813)    1.256  14.105   (   0.146    0.000   -1.070)    1.080  14.432   (  -0.206    0.000    5.098)    5.102======================= Grid point 183 (36/84) =======================q-point: ( 0.08  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.554   (   8.385    8.385    8.961)   14.863   0.570   (   4.662    4.662   12.594)   14.215   1.227   (  23.896   23.896    2.129)   33.861   1.744   (   6.684    6.684   -6.215)   11.313   1.822   (   9.206    9.206   -0.617)   13.034   2.062   (   2.419    2.419    2.071)    3.999   2.158   (   3.941    3.941   -1.864)    5.877   2.480   (  -3.403   -3.403   -0.032)    4.813   2.575   (   0.980    0.980    0.885)    1.644   2.704   (   9.212    9.212    0.373)   13.033   3.471   (   1.817    1.817    1.960)    3.232   4.902   ( -10.655  -10.655    1.231)   15.119   6.667   (   2.249    2.249    1.569)    3.547   8.070   (  -6.281   -6.281   -0.671)    8.908   8.410   (   7.997    7.997    1.189)   11.372   9.095   (  40.453   40.453  -48.223)   74.822   9.247   (  -2.590   -2.590    0.762)    3.740  11.510   (  -7.238   -7.238    0.496)   10.248  11.863   (  -5.415   -5.415   -0.911)    7.712  12.184   (   0.678    0.678   29.218)   29.233  12.769   (  -5.706   -5.706    0.099)    8.070  13.051   (  13.151   13.151   -3.592)   18.941  14.647   (   2.662    2.662  -19.040)   19.409  14.705   (  -8.969   -8.969   -1.304)   12.751======================= Grid point 184 (37/84) =======================q-point: ( 0.15  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.730   (  12.743    4.818    5.272)   14.607   0.890   (  20.507   -6.292    5.099)   22.048   1.677   (  20.582   16.594   -2.235)   26.533   1.846   (   4.944   10.574   -4.618)   12.553   2.082   (  11.860    2.523    1.883)   12.271   2.122   (   4.565    4.906    1.796)    6.938   2.249   (   6.511    2.722   -3.803)    8.017   2.432   (  -0.873   -2.325    1.442)    2.872   2.615   (   2.907    7.465    0.460)    8.024   2.944   (  14.132    8.613    0.549)   16.559   3.563   (   8.258    1.235   -0.692)    8.378   4.616   ( -17.618   -5.100    1.989)   18.449   6.857   (  15.374  -11.710    1.525)   19.386   7.881   ( -12.932   -3.397    1.199)   13.424   8.663   (  17.584   -4.865    0.961)   18.270   9.187   (  -3.177   -6.219    1.084)    7.067   9.730   (  24.578   14.801  -26.597)   39.122  11.347   (  -9.388   -9.268   -0.009)   13.192  11.775   (  -3.557  -14.982   -1.051)   15.435  12.056   ( -12.564   -1.497   26.380)   29.257  12.639   (  -7.777   -9.148    0.560)   12.020  13.394   (  20.318   18.084    2.293)   27.297  14.445   ( -15.513   -7.433   -0.351)   17.206  14.719   (   3.058   17.042  -22.594)   28.465======================= Grid point 185 (38/84) =======================q-point: ( 0.23  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.020   (  15.928    3.050    2.350)   16.387   1.283   (  18.710   -6.346    2.120)   19.871   1.942   (   4.886   10.327   -3.133)   11.847   1.990   (  10.904    7.706   -3.907)   13.912   2.169   (   2.082    6.529    0.666)    6.885   2.260   (   5.864    3.259    2.541)    7.174   2.360   (   2.194   -0.750   -1.528)    2.777   2.490   (   4.419    4.708   -3.457)    7.324   2.719   (   9.596    4.333    5.542)   11.898   3.242   (  16.603    8.183   -1.205)   18.549   3.794   (  14.589    0.340   -1.726)   14.695   4.271   ( -15.779   -1.098    2.587)   16.028   7.207   (  19.216  -18.470    1.035)   26.673   7.571   ( -18.364    0.435    2.282)   18.510   9.008   (   7.257  -17.048   -0.112)   18.529   9.190   (  12.871   -9.998    1.495)   16.367  10.084   (  11.630    5.989  -12.632)   18.185  11.155   (  -9.429   -9.543   -1.035)   13.455  11.702   (  -8.042  -12.364    5.341)   15.687  11.779   ( -10.064  -10.533    8.197)   16.716  12.456   ( -10.984  -11.264    0.849)   15.756  13.784   (  18.202   18.806    5.071)   26.659  14.133   ( -15.075   -3.643    0.147)   15.510  14.762   (   1.076   22.076  -15.382)   26.928======================= Grid point 186 (39/84) =======================q-point: ( 0.31  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.322   (  14.214    1.361    2.244)   14.455   1.605   (  14.118   -3.339    0.473)   14.515   2.030   (   4.159   10.414   -1.910)   11.375   2.135   (   4.880    4.574   -1.155)    6.787   2.193   (   0.550    4.159   -0.863)    4.283   2.309   (   0.433    1.691    1.859)    2.550   2.439   (   4.709    1.983   -3.198)    6.028   2.493   (  -2.587   10.436   -3.233)   11.228   2.919   (   7.479   -0.583    7.239)   10.425   3.586   (  17.635    6.241   -4.078)   19.146   4.074   (   2.094    0.047    1.863)    2.803   4.097   (  10.620    0.462   -0.412)   10.638   7.184   ( -20.487    4.171    1.235)   20.943   7.573   (  17.562  -21.446    0.609)   27.726   8.926   (  -8.892  -15.115    0.327)   17.540   9.663   (  27.824  -22.971    0.045)   36.081  10.175   (  -2.181   -1.253   -4.717)    5.346  11.015   (  -4.304   -8.980   -0.776)    9.989  11.443   ( -12.338    2.162    4.907)   13.453  11.728   (  -0.434  -30.443   -0.928)   30.461  12.209   ( -14.266  -12.805    0.929)   19.192  13.847   ( -10.161    1.724   -1.097)   10.364  14.121   (  12.437   18.123    6.406)   22.894  14.761   (  -1.024   23.936   -8.023)   25.266======================= Grid point 187 (40/84) =======================q-point: ( 0.38  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.546   (   8.194   -0.079    3.280)    8.827   1.833   (   9.177   -1.490   -0.540)    9.313   2.101   (   2.933   11.452   -0.758)   11.846   2.185   (   0.175    4.975    1.251)    5.133   2.211   (   1.730    3.059   -1.433)    3.796   2.322   (   1.000    1.552    0.580)    1.935   2.411   (  -4.518    8.041   -1.397)    9.329   2.520   (   2.890    5.673   -6.593)    9.165   2.995   (   1.784   -1.093    7.397)    7.688   3.886   (  12.210    4.452   -6.482)   14.524   4.177   (   9.370   -0.408    3.927)   10.168   4.342   (  10.239    0.344   -1.050)   10.298   6.806   ( -17.093    7.211   -0.827)   18.570   7.875   (  12.808  -21.105    0.337)   24.690   8.726   ( -10.380  -13.992   -1.444)   17.482  10.047   (  -8.917   -8.414   -1.620)   12.367  10.200   (  26.077  -30.559   -0.247)   40.173  10.976   (  -0.448   -8.503    1.367)    8.624  11.265   (  -6.141    1.089    0.218)    6.240  11.724   (  -0.198  -35.661   -1.063)   35.678  11.908   ( -15.936  -12.892    0.882)   20.517  13.698   (  -4.850    3.507   -0.835)    6.043  14.310   (   6.587   20.720    4.985)   22.306  14.733   (  -1.445   24.946   -2.844)   25.149======================= Grid point 188 (41/84) =======================q-point: ( 0.46  0.08  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.642   (   2.025   -0.519    3.927)    4.449   1.959   (   3.454   -1.527   -1.093)    3.932   2.140   (   1.003   13.352    0.383)   13.395   2.182   (  -0.178    2.928    0.490)    2.974   2.242   (   1.527    4.707   -0.797)    5.013   2.326   (  -3.201    8.372   -0.206)    8.965   2.341   (   0.377    1.803   -0.149)    1.848   2.562   (   1.269    5.281   -7.681)    9.407   3.025   (   1.143   -1.203    6.857)    7.055   4.051   (   4.425    3.599   -6.986)    9.018   4.334   (   4.808   -0.843    3.889)    6.241   4.479   (   3.601    0.416   -0.729)    3.697   6.561   (  -6.908    9.018   -2.433)   11.618   8.054   (   4.978  -18.702    0.200)   19.354   8.569   (  -4.600  -12.643   -3.112)   13.809   9.900   (  -4.586  -13.556   -0.706)   14.328  10.629   (  15.049  -37.900   -0.777)   40.786  10.976   (   0.104   -8.642    3.084)    9.176  11.178   (  -3.136   -1.979   -0.419)    3.731  11.625   (  -9.701   -8.425    0.234)   12.851  11.726   (  -0.323  -38.280   -0.930)   38.293  13.645   (  -1.083    4.490   -0.574)    4.654  14.386   (   1.596   22.841    3.066)   23.101  14.714   (  -0.438   25.097    0.488)   25.105======================= Grid point 197 (42/84) =======================q-point: ( 0.15  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.886   (  10.389   10.389    3.109)   15.018   0.934   (  11.122   11.122    5.073)   16.527   1.967   (  11.852   11.852   -2.903)   17.010   2.050   (   7.478    7.478   -6.293)   12.306   2.173   (   7.335    7.335    2.731)   10.727   2.251   (   8.025    8.025    2.050)   11.532   2.319   (   4.599    4.599   -2.553)    6.987   2.397   (  -0.318   -0.318    2.805)    2.841   2.771   (   7.364    7.364   -2.761)   10.774   3.174   (  13.951   13.951    0.671)   19.742   3.629   (   6.368    6.368   -0.984)    9.059   4.460   ( -10.658  -10.658    2.644)   15.303   6.717   (  -1.159   -1.159    2.161)    2.712   7.758   (  -9.014   -9.014    0.107)   12.749   8.627   (   2.233    2.233    2.782)    4.208   9.026   (  -9.236   -9.236    0.721)   13.082  10.022   (  13.186   13.186  -17.235)   25.393  11.130   ( -11.891  -11.891    0.748)   16.833  11.490   ( -13.656  -13.656   -2.137)   19.430  11.955   (  -9.074   -9.074   19.565)   23.398  12.430   ( -11.769  -11.769    0.866)   16.667  13.797   (  21.736   21.736    6.229)   31.364  14.253   ( -11.327  -11.327   -0.152)   16.019  15.025   (  12.449   12.449  -20.445)   26.980======================= Grid point 198 (43/84) =======================q-point: ( 0.23  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.123   (  13.286    7.390    2.003)   15.334   1.245   (  18.370    3.552    2.937)   18.939   2.129   (   5.360    6.205   -3.322)    8.847   2.149   (   4.140    8.472   -3.729)   10.140   2.318   (   4.345    6.007    0.306)    7.420   2.343   (   0.156    3.818    3.250)    5.016   2.364   (   5.600    5.999    0.357)    8.214   2.568   (  12.994    2.743   -1.349)   13.348   2.868   (   2.680    9.555   -1.775)   10.082   3.461   (  15.104   13.318   -0.405)   20.141   3.830   (  13.656    4.158   -1.193)   14.324   4.230   ( -11.440   -3.304    3.093)   12.303   6.799   (   8.449  -20.742    1.431)   22.443   7.545   ( -12.501   -4.097    1.204)   13.210   8.684   (  -1.973  -12.648    0.184)   12.802   8.916   (   3.269  -17.274    1.560)   17.650  10.193   (   2.997    3.903   -9.382)   10.594  10.905   ( -10.412  -15.234   -0.133)   18.453  11.254   ( -10.113  -23.279   -1.415)   25.420  11.718   ( -14.151   -4.803   10.889)   18.490  12.180   ( -13.953  -16.403    1.092)   21.562  14.027   (  -9.851   -4.664   -0.285)   10.903  14.213   (  18.960   21.741    7.032)   29.692  15.194   (   5.289   19.527  -13.869)   24.528======================= Grid point 199 (44/84) =======================q-point: ( 0.31  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.375   (  11.756    4.054    2.282)   12.643   1.592   (  16.121    2.738    0.764)   16.370   2.208   (   3.078    7.533   -2.049)    8.391   2.219   (   2.827    2.757   -3.297)    5.144   2.345   (   0.908    3.429    1.985)    4.065   2.374   (   2.954   11.066    1.312)   11.528   2.486   (   4.758    4.358   -0.805)    6.502   2.727   (   1.938   10.751   -4.983)   12.008   2.925   (   3.852    2.381    3.864)    5.953   3.758   (  14.791   10.781   -2.624)   18.490   4.073   (   0.696   -0.249    2.416)    2.526   4.127   (  11.969    3.166   -0.060)   12.380   7.015   (  13.057  -33.520    0.410)   35.975   7.277   ( -14.560    3.994    1.264)   15.150   8.537   (  -9.676  -22.429   -1.103)   24.452   9.147   (  16.271  -27.308    0.752)   31.796  10.098   ( -12.537   -7.338   -3.967)   15.058  10.757   (  -4.444  -17.450   -0.713)   18.022  11.097   (  -6.236  -30.355   -0.567)   30.994  11.445   ( -12.915   -3.469    4.537)   14.121  11.884   ( -16.096  -19.879    1.039)   25.600  13.867   (  -5.912    1.435   -0.475)    6.102  14.542   (  14.053   22.642    6.596)   27.453  15.251   (   0.978   23.362   -8.235)   24.790======================= Grid point 200 (45/84) =======================q-point: ( 0.38  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.557   (   6.548    1.256    2.882)    7.263   1.851   (  10.068    3.759   -0.797)   10.776   2.236   (   0.529    3.952   -0.806)    4.068   2.284   (   3.315    5.081   -2.263)    6.475   2.368   (   1.311    3.589    1.054)    3.964   2.442   (   3.184   13.810   -0.300)   14.175   2.533   (   0.233    7.697   -0.789)    7.740   2.717   (  -0.966   12.215   -4.269)   12.976   2.977   (   1.478    0.043    5.014)    5.227   4.006   (   9.962    7.173   -3.552)   12.779   4.165   (   8.754   -0.515    2.668)    9.166   4.366   (  10.041    2.579   -0.398)   10.375   6.997   ( -12.874   11.536   -0.500)   17.293   7.296   (  14.850  -36.536   -0.160)   39.439   8.355   (  -8.326  -22.919   -1.692)   24.443   9.478   (  16.304  -40.087   -0.170)   43.276   9.768   ( -18.949  -20.556   -1.684)   28.008  10.708   (  -1.177  -19.433    0.676)   19.480  10.998   (  -4.031  -34.423   -0.402)   34.661  11.232   (  -8.871   -4.924    1.351)   10.236  11.574   ( -14.559  -21.109    0.790)   25.655  13.795   (  -1.774    6.694   -0.528)    6.945  14.757   (   7.883   23.076    5.945)   25.100  15.251   (  -0.514   25.263   -4.718)   25.705======================= Grid point 201 (46/84) =======================q-point: ( 0.46  0.15  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.634   (   1.692   -0.081    2.970)    3.419   1.982   (   3.368    4.126   -1.608)    5.564   2.242   (   0.150    3.675   -0.199)    3.684   2.334   (   1.477    3.940   -1.793)    4.574   2.394   (   1.210    5.047    0.204)    5.194   2.466   (  -0.632   11.990   -1.731)   12.131   2.526   (  -0.200   15.795    1.359)   15.854   2.712   (   0.181    9.534   -4.336)   10.476   3.001   (   0.914   -0.582    4.939)    5.056   4.144   (   3.849    5.054   -2.930)    6.996   4.315   (   4.644   -0.465    1.633)    4.944   4.501   (   3.506    2.058   -0.219)    4.072   6.810   (  -5.353   15.860   -2.196)   16.883   7.544   (   8.252  -32.124   -0.511)   33.171   8.236   (  -3.308  -20.392   -2.212)   20.776   9.454   ( -10.220  -32.705   -0.696)   34.272   9.731   (   7.782  -50.746   -0.560)   51.342  10.698   (  -0.100  -20.205    2.043)   20.308  10.940   (  -1.699  -36.796   -0.733)   36.843  11.096   (  -4.465   -6.634    0.688)    8.026  11.366   (  -5.616  -19.181    0.294)   19.988  13.780   (  -0.150    9.311   -0.452)    9.323  14.855   (   2.402   23.559    5.294)   24.266  15.242   (  -0.299   25.998   -2.908)   26.162======================= Grid point 211 (47/84) =======================q-point: ( 0.23  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.293   (   9.530    9.530    1.974)   13.622   1.420   (  13.454   13.454    2.575)   19.200   2.214   (   3.400    3.400   -3.642)    6.032   2.287   (   5.065    5.065   -3.170)    7.833   2.361   (   1.166    1.166    1.335)    2.121   2.391   (   3.928    3.928   -0.401)    5.570   2.571   (   4.861    4.861    3.642)    7.780   2.675   (  12.556   12.556   -0.223)   17.758   3.008   (   3.627    3.627   -5.579)    7.579   3.731   (  13.693   13.693    0.706)   19.378   3.978   (  10.400   10.400   -1.531)   14.788   4.146   (  -4.347   -4.347    3.391)    7.021   6.465   ( -12.289  -12.289    1.089)   17.414   7.413   (  -8.027   -8.027    0.500)   11.363   8.520   ( -16.219  -16.219   -0.096)   22.937   8.600   (  -2.925   -2.925    1.367)    4.357  10.189   (  -5.866   -5.866   -6.308)   10.421  10.583   ( -17.061  -17.061   -0.293)   24.129  10.853   ( -16.999  -16.999   -0.167)   24.041  11.559   ( -11.088  -11.088    7.001)   17.173  11.829   ( -19.246  -19.246    0.827)   27.230  13.944   (  -2.879   -2.879   -0.365)    4.087  14.628   (  19.506   19.506    6.596)   28.363  15.488   (  10.564   10.564  -10.188)   18.083======================= Grid point 212 (48/84) =======================q-point: ( 0.31  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.472   (   8.145    5.557    2.163)   10.095   1.719   (  15.621    9.964    1.182)   18.566   2.273   (   2.573    3.003   -3.230)    5.106   2.345   (   1.667    6.171   -1.191)    6.503   2.397   (   2.465    2.690    0.821)    3.740   2.482   (   4.871    2.333   -2.091)    5.792   2.648   (   4.585    6.958    3.899)    9.200   2.918   (   9.355   11.946   -0.951)   15.203   3.023   (  -0.750    5.546   -3.515)    6.609   3.979   (  10.704   10.773   -0.304)   15.190   4.059   (   0.491   -0.957    2.746)    2.949   4.240   (  12.177    8.374   -0.801)   14.800   6.352   (   1.462  -31.431    0.262)   31.466   7.283   (  -4.750   -4.212    0.664)    6.383   8.148   ( -17.908  -16.309   -1.038)   24.244   8.670   (   5.524  -19.383    0.680)   20.166   9.881   ( -25.817  -14.385   -3.241)   29.731  10.317   (  -8.950  -27.790   -1.094)   29.216  10.595   ( -10.032  -20.669    0.710)   22.986  11.312   ( -13.389   -9.927    3.866)   17.110  11.442   ( -19.732  -25.019    0.720)   31.872  13.925   (   0.781    4.587   -0.597)    4.691  14.957   (  13.951   19.103    5.853)   24.369  15.625   (   4.084   14.717   -6.850)   16.739======================= Grid point 213 (49/84) =======================q-point: ( 0.38  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.595   (   4.256    2.514    2.235)    5.425   1.982   (  10.491    9.406   -0.409)   14.097   2.307   (   0.904    3.546   -1.869)    4.109   2.373   (   1.254    3.784    0.599)    4.031   2.457   (   3.489    5.197   -0.886)    6.322   2.565   (   3.265    2.990   -2.236)    4.960   2.748   (   5.167    7.796    2.707)    9.737   2.957   (  -2.078   10.454   -2.615)   10.974   3.063   (   2.570   11.861    0.987)   12.176   4.084   (   4.997   -0.704    1.854)    5.376   4.218   (  10.197    6.669   -0.952)   12.221   4.453   (   8.708    6.279   -0.379)   10.743   6.526   (  15.615  -38.916   -0.452)   41.934   7.212   (  -3.378    8.181   -0.040)    8.851   7.879   (  -8.371  -24.070   -1.405)   25.523   8.753   (   2.122  -30.914   -0.204)   30.988   9.272   ( -33.890  -28.050   -1.489)   44.018  10.201   (  -3.299  -32.013   -0.356)   32.185  10.433   (  -6.970  -24.504    0.356)   25.478  11.060   ( -12.074  -19.774    1.317)   23.206  11.097   ( -13.949  -20.191    1.499)   24.587  13.960   (   2.296    9.822   -0.763)   10.116  15.171   (   7.951   18.728    5.262)   21.015  15.670   (   1.008   17.235   -4.652)   17.881======================= Grid point 214 (50/84) =======================q-point: ( 0.46  0.23  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.645   (   1.102    1.269    2.082)    2.676   2.117   (   3.407    9.332   -1.365)   10.027   2.315   (   0.063    3.658   -0.815)    3.748   2.391   (   0.555    2.512    0.899)    2.725   2.514   (   1.791    6.272   -1.923)    6.800   2.603   (   0.719    4.272   -1.423)    4.560   2.825   (   2.040    8.881    0.820)    9.150   2.927   (  -0.649   10.241   -0.532)   10.275   3.085   (   0.268   13.870    1.420)   13.945   4.178   (   3.224   -1.992    1.921)    4.249   4.373   (   4.663    6.443   -1.923)    8.183   4.566   (   2.815    4.631    0.065)    5.420   6.869   (  14.424  -34.227   -1.119)   37.159   7.146   (  -2.482   16.631   -1.599)   16.891   7.817   (   0.160  -20.949   -0.707)   20.962   8.681   ( -20.295  -42.782   -0.374)   47.353   8.756   (  -1.681  -46.061   -0.841)   46.099  10.178   (   0.292  -32.735    0.729)   32.744  10.323   (  -3.744  -27.922   -0.491)   28.177  10.881   (  -4.716  -30.311    0.522)   30.680  10.917   (  -4.954  -11.551    1.727)   12.686  13.996   (   1.125   12.280   -0.818)   12.358  15.272   (   2.512   18.587    4.885)   19.382  15.679   (   0.117   18.337   -3.557)   18.680======================= Grid point 225 (51/84) =======================q-point: ( 0.31  0.31  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.577   (   4.823    4.823    1.915)    7.085   1.959   (  13.568   13.568    1.058)   19.217   2.332   (   2.673    2.673   -2.233)    4.390   2.431   (   2.194    2.194    0.797)    3.203   2.453   (   3.306    3.306   -1.487)    4.907   2.537   (   3.731    3.731   -1.690)    5.540   2.753   (   5.420    5.420    3.526)    8.438   3.093   (   1.894    1.894   -4.159)    4.947   3.183   (  13.038   13.038   -0.164)   18.439   4.027   (  -0.807   -0.807    2.336)    2.600   4.174   (   8.110    8.110    0.318)   11.473   4.429   (  10.023   10.023   -1.184)   14.224   5.900   ( -12.231  -12.231    0.219)   17.298   7.177   (  -3.516   -3.516   -0.210)    4.976   7.823   ( -18.328  -18.328   -0.498)   25.925   8.457   (  -3.751   -3.751    0.741)    5.356   9.532   ( -26.605  -26.605   -3.213)   37.762   9.728   ( -25.378  -25.378   -0.062)   35.890  10.255   ( -15.055  -15.055    0.664)   21.301  10.917   ( -27.129  -27.129    0.120)   38.367  11.085   ( -13.071  -13.071    3.447)   18.803  14.038   (   6.529    6.529   -0.944)    9.282  15.276   (  13.351   13.351    4.966)   19.524  15.835   (   7.244    7.244   -5.249)   11.511======================= Grid point 226 (52/84) =======================q-point: ( 0.38  0.31  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.652   (   2.763    3.204    1.694)    4.557   2.197   (   9.544   11.517    0.240)   14.959   2.371   (   1.253    2.714   -1.048)    3.168   2.440   (  -0.333    3.370    1.345)    3.643   2.538   (   4.851    1.981   -1.184)    5.372   2.609   (   3.031    2.564   -1.481)    4.238   2.870   (   6.006    6.462    2.526)    9.176   3.108   (  -0.037    6.835   -2.476)    7.269   3.381   (   6.660   16.810   -0.156)   18.082   4.045   (   2.990   -2.296    1.343)    4.003   4.332   (   8.980    4.110   -0.365)    9.882   4.599   (   6.596    7.951   -0.667)   10.353   5.917   (  13.288  -19.787   -0.089)   23.835   7.231   (   9.826   -5.753   -0.592)   11.401   7.478   ( -15.164  -16.722   -0.396)   22.577   8.374   (  -8.249   -7.094    0.558)   10.894   8.757   ( -45.737  -23.242   -1.540)   51.326   9.472   (  -5.940  -41.937   -0.422)   42.357   9.986   ( -12.466  -21.911    0.086)   25.210  10.437   ( -19.906  -35.403    0.126)   40.616  10.834   ( -11.838  -13.860    2.959)   18.466  14.160   (   5.490   10.120   -1.315)   11.588  15.479   (   7.554   12.773    4.300)   15.450  15.929   (   2.814    9.567   -4.042)   10.761======================= Grid point 227 (53/84) =======================q-point: ( 0.46  0.31  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.687   (   0.855    3.001    1.563)    3.490   2.317   (   3.045    9.566   -0.289)   10.043   2.380   (  -0.213    2.941   -0.387)    2.974   2.434   (  -0.132    1.773    2.065)    2.725   2.618   (   2.788    3.209   -1.690)    4.575   2.644   (   0.470    1.688   -1.101)    2.069   2.958   (   2.393    6.871    1.366)    7.403   3.106   (  -0.020    9.991   -0.727)   10.017   3.452   (   1.306   19.294   -0.145)   19.339   4.115   (   2.684   -3.750    0.720)    4.668   4.503   (   6.397    6.011   -0.953)    8.830   4.673   (   1.147    5.885   -0.090)    5.997   6.317   (  21.156  -19.235   -1.095)   28.614   7.256   (  -5.705  -11.846   -0.462)   13.156   7.562   (  20.769   -3.748   -0.154)   21.105   7.922   ( -32.601  -32.511   -0.128)   46.041   8.077   ( -20.504  -16.163   -0.614)   26.116   9.430   (   0.358  -43.333   -0.218)   43.335   9.788   (  -6.546  -27.476   -0.435)   28.248  10.191   (  -5.475  -39.728    0.746)   40.111  10.657   (  -5.286  -14.528    2.578)   15.673  14.237   (   2.106   11.733   -1.545)   12.020  15.576   (   2.399   12.468    3.946)   13.296  15.961   (   0.675   10.728   -3.410)   11.277======================= Grid point 239 (54/84) =======================q-point: ( 0.38  0.38  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.713   (   2.795    2.795    1.481)    4.221   2.369   (   4.834    4.834    0.970)    6.905   2.410   (   1.386    1.386    0.183)    1.969   2.485   (  -0.234   -0.234    1.328)    1.369   2.574   (   4.326    4.326   -0.035)    6.117   2.669   (   2.775    2.775   -1.612)    4.243   3.008   (   6.940    6.940    2.192)   10.056   3.237   (   5.859    5.859   -2.154)    8.561   3.661   (  10.815   10.815   -0.780)   15.315   4.017   (   0.000    0.000    0.148)    0.148   4.426   (   7.616    7.616   -0.859)   10.805   4.734   (   5.129    5.129   -0.404)    7.264   5.817   (   8.838    8.838    0.255)   12.501   7.131   (  -7.360   -7.360    0.152)   10.410   7.225   (  -2.145   -2.145   -1.119)    3.233   8.279   ( -29.458  -29.458   -2.214)   41.719   8.355   (  -5.844   -5.844    1.586)    8.415   8.701   ( -26.252  -26.252    0.407)   37.128   9.571   ( -20.045  -20.045   -0.046)   28.348   9.758   ( -31.021  -31.021   -0.059)   43.871  10.566   ( -12.739  -12.739    3.152)   18.290  14.335   (   7.379    7.379   -1.907)   10.608  15.673   (   7.135    7.135    3.530)   10.690  16.061   (   4.341    4.341   -3.453)    7.044======================= Grid point 240 (55/84) =======================q-point: ( 0.46  0.38  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.753   (   1.125    3.389    1.474)    3.863   2.410   (   0.689    1.010    1.660)    2.062   2.431   (   0.791    0.879    1.071)    1.595   2.457   (  -1.625   -0.048    2.161)    2.704   2.660   (   2.710    3.787   -0.659)    4.703   2.706   (   0.936    2.807   -1.780)    3.453   3.110   (   2.872    8.030    1.807)    8.718   3.317   (   2.114   10.796   -1.568)   11.112   3.783   (   2.289   13.972   -1.100)   14.201   4.048   (   2.037   -2.742   -0.835)    3.517   4.620   (   9.237    6.783   -0.777)   11.487   4.777   (  -0.852    4.278   -0.218)    4.368   6.189   (  24.258    6.518   -0.796)   25.131   6.964   (  -5.206  -13.678    1.143)   14.680   7.422   (   5.018   -4.742   -2.127)    7.224   7.462   ( -28.061  -13.086   -0.820)   30.974   8.149   ( -16.331   12.321    1.925)   20.548   8.516   (  -0.949  -48.573   -0.901)   48.591   9.253   ( -10.237  -26.192    0.271)   28.122   9.363   (  -8.343  -43.416    0.870)   44.219  10.375   (  -5.651  -13.400    2.983)   14.846  14.436   (   2.677    8.206   -2.276)    8.927  15.764   (   2.259    6.900    3.120)    7.902  16.113   (   1.191    5.183   -3.118)    6.165======================= Grid point 253 (56/84) =======================q-point: ( 0.46  0.46  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.804   (   1.481    1.481    1.588)    2.628   2.407   (  -0.321   -0.321    3.237)    3.268   2.421   (  -1.457   -1.457    2.974)    3.618   2.451   (   0.182    0.182    1.499)    1.521   2.730   (   2.192    2.192   -1.424)    3.412   2.742   (   0.955    0.955   -1.561)    2.064   3.232   (   3.527    3.527    2.504)    5.582   3.481   (   4.773    4.773   -2.270)    7.122   3.973   (   4.421    4.421   -1.859)    6.523   4.017   (  -0.023   -0.023   -1.765)    1.765   4.760   (   6.110    6.110   -0.548)    8.658   4.824   (   0.045    0.045   -0.262)    0.269   6.541   (  26.538   26.538   -1.220)   37.550   6.775   (  -4.935   -4.935    1.776)    7.201   7.230   ( -16.156  -16.156   -1.481)   22.896   7.564   (  13.311   13.311   -3.308)   19.113   7.672   ( -27.544  -27.544   -1.385)   38.978   8.320   (   0.032    0.032    3.400)    3.400   8.620   ( -24.327  -24.327    1.503)   34.437   8.844   ( -12.847  -12.847    1.466)   18.227  10.175   (  -5.944   -5.944    2.912)    8.897  14.546   (   2.923    2.923   -2.717)    4.947  15.851   (   2.177    2.177    2.661)    4.069  16.176   (   1.493    1.493   -2.931)    3.612======================= Grid point 338 (57/84) =======================q-point: ( 0.00  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 84Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.595   (   0.000    0.000    6.967)    6.967   0.595   (   0.000    0.000    6.967)    6.967   1.078   (  -0.000    0.000   13.164)   13.164   1.485   (   0.000    0.000   -3.991)    3.991   1.485   (   0.000    0.000   -3.991)    3.991   2.044   (  -0.000    0.000    0.591)    0.591   2.044   (   0.000    0.000    0.591)    0.591   2.518   (  -0.000    0.000    0.102)    0.102   2.518   (   0.000    0.000    0.102)    0.102   2.617   (   0.000    0.000   -9.755)    9.755   3.710   (  -0.000   -0.000    3.250)    3.250   5.158   (   0.000    0.000   -0.008)    0.008   6.632   (   0.000    0.000    0.523)    0.523   6.632   (   0.000   -0.000    0.523)    0.523   8.181   (  -0.000    0.000   -0.619)    0.619   8.181   (   0.000    0.000   -0.619)    0.619   9.296   (   0.000    0.000    0.317)    0.317  11.658   (  -0.000   -0.000   -0.126)    0.126  11.969   (   0.000    0.000    0.437)    0.437  11.969   (   0.000    0.000    0.437)    0.437  12.873   (  -0.000   -0.000   -0.241)    0.241  12.873   (   0.000    0.000   -0.241)    0.241  14.802   (  -0.000   -0.000    0.972)    0.972  14.899   (   0.000    0.000   -0.541)    0.541======================= Grid point 339 (58/84) =======================q-point: ( 0.08  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.781   (  15.969    0.000    5.040)   16.746   0.804   (   8.625    0.000   12.028)   14.801   0.987   (   5.478    0.000   -1.354)    5.643   1.531   (   4.447    0.000   -3.673)    5.768   1.699   (  16.722    0.000   -0.921)   16.748   2.053   (   0.931    0.000    0.580)    1.096   2.173   (   9.884    0.000    5.360)   11.244   2.469   (  -8.406    0.000   -8.351)   11.849   2.507   (  -1.036    0.000    0.080)    1.039   2.629   (  10.391    0.000    0.111)   10.392   3.552   (  -6.581    0.000    2.748)    7.132   5.024   ( -12.885    0.000    0.398)   12.891   6.742   (  10.874    0.000    0.521)   10.886   7.738   (  79.815    0.000  -13.506)   80.950   8.137   (  -5.015    0.000    0.202)    5.019   8.322   (  14.090    0.000   -0.608)   14.103   9.290   (  -0.543    0.000    0.446)    0.703  11.593   (  -6.170    0.000   -0.145)    6.172  11.975   (   0.610    0.000    0.371)    0.714  12.294   (  29.435    0.000    3.310)   29.621  12.841   (  -3.252    0.000   -0.181)    3.257  12.862   (  -1.681    0.000   -1.188)    2.058  14.514   ( -23.773    0.000   -1.219)   23.805  14.791   ( -10.620    0.000   -0.383)   10.627======================= Grid point 340 (59/84) =======================q-point: ( 0.15  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.847   (   5.636    0.000    5.784)    8.076   1.120   (  17.287    0.000    2.941)   17.535   1.371   (  22.537    0.000   -3.588)   22.821   1.643   (   6.509    0.000   -2.872)    7.114   2.051   (  17.681    0.000   -1.440)   17.740   2.078   (   1.489    0.000    0.550)    1.588   2.255   (   0.036    0.000    3.957)    3.957   2.477   (  -2.017    0.000    0.016)    2.017   2.521   (  11.806    0.000   -1.433)   11.893   2.871   (  13.047    0.000    0.087)   13.047   3.544   (   5.543    0.000   -0.196)    5.547   4.697   ( -19.330    0.000    0.930)   19.352   7.035   (  18.356    0.000    0.510)   18.363   7.959   ( -13.410    0.000    1.531)   13.497   8.712   (  25.097    0.000   -0.567)   25.104   9.033   (  49.699    0.000  -12.265)   51.190   9.280   (  -0.391    0.000    0.525)    0.655  11.430   ( -10.005    0.000   -0.448)   10.015  11.994   (   1.405    0.000    0.193)    1.419  12.543   ( -17.529    0.000   15.959)   23.706  12.745   (  -6.533    0.000   -0.023)    6.533  13.093   (  31.791    0.000   -5.546)   32.271  14.155   (  -8.189    0.000   -9.752)   12.734  14.512   ( -16.832    0.000   -0.247)   16.834======================= Grid point 341 (60/84) =======================q-point: ( 0.23  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.055   (  14.364    0.000    2.155)   14.524   1.431   (  14.279    0.000    1.584)   14.366   1.755   (  16.104    0.000   -3.848)   16.558   1.771   (   6.339    0.000   -1.857)    6.605   2.108   (   1.550    0.000    0.508)    1.631   2.276   (   1.904    0.000    1.151)    2.225   2.315   (   7.981    0.000   -1.623)    8.144   2.429   (  -2.896    0.000   -0.084)    2.897   2.809   (  16.643    0.000    3.956)   17.107   3.129   (  14.045    0.000   -1.307)   14.106   3.751   (  14.397    0.000   -1.322)   14.457   4.326   ( -16.751    0.000    1.315)   16.803   7.425   (  20.553    0.000    0.479)   20.558   7.617   ( -21.087    0.000    1.657)   21.152   9.242   (  -3.723    0.000    1.197)    3.911   9.272   (  31.317    0.000   -0.477)   31.320   9.761   (  25.950    0.000   -6.807)   26.828  11.228   (  -9.787    0.000   -0.676)    9.810  12.025   ( -29.754    0.000    8.962)   31.075  12.030   (   2.295    0.000   -0.041)    2.295  12.585   (  -9.922    0.000    0.171)    9.923  13.628   (  19.408    0.000    2.591)   19.580  14.144   (  -4.314    0.000   -3.745)    5.713  14.216   (  -3.506    0.000   -6.834)    7.681======================= Grid point 342 (61/84) =======================q-point: ( 0.31  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.359   (  15.525    0.000    1.456)   15.593   1.677   (  10.982    0.000    0.632)   11.000   1.883   (   5.023    0.000   -0.852)    5.094   2.004   (   9.779    0.000   -2.968)   10.220   2.136   (   1.235    0.000    0.454)    1.316   2.302   (   0.647    0.000    0.131)    0.660   2.365   (  -3.634    0.000   -0.200)    3.640   2.385   (   0.176    0.000   -1.224)    1.236   3.081   (   8.697    0.000    4.972)   10.018   3.428   (  15.813    0.000   -3.344)   16.163   4.052   (  15.195    0.000   -1.205)   15.242   4.128   (  -1.489    0.000    1.382)    2.032   7.164   ( -24.255    0.000    0.720)   24.266   7.806   (  17.754    0.000    0.426)   17.760   9.102   ( -10.874    0.000    0.664)   10.895   9.906   (  32.756    0.000   -0.334)   32.758  10.081   (   7.147    0.000   -2.901)    7.713  11.090   (  -3.364    0.000   -0.321)    3.379  11.511   ( -20.920    0.000    3.199)   21.163  12.081   (   2.811    0.000   -0.256)    2.823  12.356   ( -13.502    0.000    0.322)   13.506  13.743   (  -6.297    0.000   -1.203)    6.410  14.071   (   3.921    0.000    4.687)    6.111  14.334   (   7.199    0.000   -5.849)    9.275======================= Grid point 343 (62/84) =======================q-point: ( 0.38  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.620   (  10.174    0.000    2.087)   10.386   1.863   (   8.096    0.000    0.003)    8.096   1.964   (   3.215    0.000   -0.061)    3.216   2.150   (   5.233    0.000   -1.050)    5.337   2.155   (   0.745    0.000    0.389)    0.840   2.289   (  -4.093    0.000   -0.292)    4.103   2.307   (   0.030    0.000   -0.019)    0.036   2.360   (  -1.758    0.000   -3.040)    3.512   3.149   (   0.565    0.000    3.967)    4.007   3.710   (  12.167    0.000   -3.680)   12.712   4.241   (   9.722    0.000    1.338)    9.814   4.311   (  10.741    0.000   -0.890)   10.778   6.714   ( -20.419    0.000   -0.599)   20.428   8.098   (  11.656    0.000    0.367)   11.662   8.844   ( -13.875    0.000   -0.716)   13.894  10.091   (  -3.850    0.000   -1.056)    3.992  10.528   (  30.260    0.000   -0.173)   30.261  11.088   (   1.957    0.000    0.946)    2.174  11.227   (  -8.686    0.000   -0.116)    8.686  12.055   ( -16.912    0.000    0.574)   16.922  12.136   (   2.235    0.000   -0.583)    2.310  13.621   (  -4.082    0.000   -0.977)    4.197  14.117   (  -0.267    0.000    1.189)    1.219  14.459   (   5.242    0.000    0.480)    5.264======================= Grid point 344 (63/84) =======================q-point: ( 0.46  0.00  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.739   (   2.461    0.000    2.826)    3.747   1.987   (   3.793    0.000   -0.011)    3.793   2.009   (   1.506    0.000    0.017)    1.506   2.165   (   0.216    0.000    0.358)    0.418   2.210   (   1.169    0.000    1.914)    2.243   2.214   (  -2.902    0.000   -0.300)    2.917   2.309   (   0.306    0.000   -0.632)    0.702   2.343   (  -0.126    0.000   -5.640)    5.641   3.162   (   0.732    0.000    3.213)    3.296   3.878   (   4.533    0.000   -3.569)    5.770   4.395   (   4.465    0.000    0.894)    4.553   4.456   (   3.807    0.000   -0.624)    3.857   6.422   (  -8.198    0.000   -1.467)    8.328   8.252   (   4.018    0.000    0.327)    4.032   8.632   (  -6.233    0.000   -1.882)    6.511  10.013   (  -2.676    0.000   -0.401)    2.706  11.053   (  21.142    0.000   -0.124)   21.142  11.124   (   0.955    0.000    1.756)    1.999  11.133   (  -1.727    0.000   -1.178)    2.090  11.707   ( -16.481    0.000    0.584)   16.491  12.170   (   0.961    0.000   -0.585)    1.125  13.584   (  -0.430    0.000   -0.505)    0.663  14.089   (  -1.442    0.000   -0.341)    1.481  14.529   (   1.828    0.000    2.805)    3.348======================= Grid point 352 (64/84) =======================q-point: ( 0.08  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.735   (   6.262    6.262    5.623)   10.490   0.822   (  -2.667   -2.667    8.256)    9.077   1.230   (  20.370   20.370   -1.745)   28.861   1.615   (   6.277    6.277   -4.223)    9.831   1.806   (  10.906   10.906   -0.488)   15.432   2.095   (   2.773    2.773    0.821)    4.006   2.229   (   3.860    3.860    7.686)    9.427   2.419   (  -0.950   -0.950   -7.193)    7.317   2.536   (   0.821    0.821   -1.441)    1.850   2.711   (   9.386    9.386    0.219)   13.275   3.518   (  -1.603   -1.603    1.553)    2.748   4.925   (  -9.987   -9.987    0.616)   14.137   6.700   (   2.851    2.851    0.993)    4.152   8.056   (  -5.980   -5.980   -0.455)    8.469   8.400   (  43.333   43.333  -11.046)   62.269   8.446   (  16.383   16.383   -2.853)   23.343   9.263   (  -2.395   -2.395    0.516)    3.426  11.517   (  -6.945   -6.945    0.092)    9.822  11.847   (  -6.478   -6.478   -0.363)    9.168  12.652   (  17.716   17.716   14.901)   29.151  12.770   (  -5.284   -5.284    0.045)    7.473  12.922   (  11.530   11.530   -8.093)   18.204  14.359   ( -11.850  -11.850   -6.242)   17.883  14.693   (  -9.631   -9.631   -0.325)   13.625======================= Grid point 353 (65/84) =======================q-point: ( 0.15  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.851   (   7.225    1.790    4.229)    8.561   0.992   (  16.771   -4.433    3.004)   17.605   1.627   (  18.569   17.568   -1.676)   25.618   1.742   (   6.744    8.800   -3.478)   11.619   2.077   (  14.870    2.904   -1.660)   15.242   2.145   (   2.476    6.082    0.732)    6.608   2.249   (   0.594   -0.087    1.940)    2.031   2.480   (   1.057    0.654    3.161)    3.397   2.583   (   8.057    4.977   -3.145)    9.978   2.952   (  13.527    8.606    0.175)   16.034   3.551   (   5.998    1.030   -0.304)    6.093   4.654   ( -16.817   -4.574    1.083)   17.461   6.889   (  14.937  -10.471    0.976)   18.267   7.904   ( -10.621   -5.568    0.637)   12.009   8.684   (  16.094   -1.712    0.732)   16.201   9.210   (  -2.106   -6.123    0.785)    6.522   9.294   (  36.934   21.214  -10.970)   43.983  11.345   ( -10.296   -8.227   -0.148)   13.180  11.755   (  -3.213  -16.999   -0.516)   17.307  12.513   ( -20.795   -5.212   11.781)   24.462  12.657   (  -7.642   -8.538    1.332)   11.536  13.474   (  31.011   30.341    4.409)   43.609  14.235   (   1.495    8.716  -16.541)   18.757  14.439   ( -15.441   -7.164   -0.231)   17.023======================= Grid point 354 (66/84) =======================q-point: ( 0.23  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.073   (  14.345    2.450    1.854)   14.670   1.327   (  16.454   -5.455    1.358)   17.388   1.874   (   6.556    9.039   -2.216)   11.384   1.921   (  11.403   11.348   -1.830)   16.191   2.186   (   1.705    7.621    0.541)    7.828   2.288   (   3.739    1.955    0.436)    4.242   2.317   (   7.026    0.508   -1.731)    7.254   2.459   (  -1.624    2.681   -0.027)    3.135   2.816   (  14.077    0.641    2.726)   14.353   3.215   (  13.827    8.768   -0.917)   16.398   3.760   (  14.401    1.112   -1.045)   14.481   4.320   ( -15.490   -0.755    1.410)   15.573   7.229   (  18.686  -17.933    0.657)   25.907   7.618   ( -18.424   -0.388    1.408)   18.482   9.006   (   8.523  -17.633   -0.065)   19.584   9.221   (  10.795   -8.914    0.978)   14.034   9.857   (  20.408    8.425   -6.125)   22.912  11.138   (  -9.928   -8.895   -0.516)   13.340  11.714   (  -1.230  -24.864   -0.241)   24.895  12.016   ( -26.951   -2.090    7.231)   27.982  12.474   ( -10.999  -10.542    0.593)   15.247  13.948   (  13.895   23.225    8.926)   28.498  14.145   (  -9.886   -1.555    1.683)   10.148  14.376   (   9.887   18.550  -17.119)   27.109======================= Grid point 355 (67/84) =======================q-point: ( 0.31  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.369   (  14.850    1.458    1.488)   14.995   1.615   (  12.743   -3.454    0.277)   13.205   1.994   (   5.586    9.292   -1.086)   10.896   2.091   (   6.410    5.390   -1.564)    8.520   2.209   (   0.643    7.661    1.026)    7.756   2.330   (   0.695    2.902    0.080)    2.985   2.391   (   1.427    1.850   -0.652)    2.425   2.429   (  -1.776    5.998   -2.191)    6.628   3.055   (   8.022   -2.660    3.999)    9.350   3.504   (  15.150    7.607   -2.610)   17.152   4.061   (  15.139    1.062   -0.861)   15.201   4.133   (  -1.715    0.413    1.370)    2.233   7.209   ( -22.379    4.272    0.763)   22.796   7.586   (  17.197  -21.597    0.373)   27.610   8.931   (  -9.987  -16.714    0.110)   19.471   9.665   (  27.295  -21.968    0.087)   35.037  10.082   (   2.591   -0.640   -2.737)    3.822  11.003   (  -3.028   -8.884   -0.305)    9.391  11.543   ( -19.178    1.567    3.155)   19.498  11.707   (  -0.436  -30.613   -0.802)   30.627  12.227   ( -14.274  -12.378    0.524)   18.901  13.829   ( -10.466    2.859   -0.493)   10.861  14.297   (  10.654   21.379    9.137)   25.574  14.546   (   6.946   20.421  -10.461)   23.973======================= Grid point 356 (68/84) =======================q-point: ( 0.38  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.615   (   9.405   -0.175    2.132)    9.645   1.822   (   8.377   -2.098   -0.374)    8.644   2.089   (   4.023   10.202   -0.313)   10.971   2.179   (   2.537    1.438   -0.817)    3.028   2.215   (   0.133    7.158    0.838)    7.208   2.333   (   0.041    2.402    0.436)    2.442   2.380   (  -2.654    5.544   -0.934)    6.217   2.408   (   0.443    8.232   -3.179)    8.836   3.123   (   0.831   -2.610    3.458)    4.411   3.775   (  11.734    6.456   -3.005)   13.726   4.233   (   9.142   -0.894    1.233)    9.268   4.320   (  10.688    1.005   -0.655)   10.755   6.789   ( -19.252    7.449   -0.530)   20.650   7.882   (  12.613  -21.587    0.215)   25.002   8.697   ( -12.416  -14.382   -0.874)   19.020  10.013   (  -7.232   -8.482   -1.067)   11.198  10.195   (  25.756  -30.422   -0.129)   39.861  11.003   (   1.819   -8.948    0.807)    9.166  11.271   (  -8.824    2.148    0.253)    9.085  11.704   (  -0.098  -35.679   -0.546)   35.684  11.924   ( -16.210  -12.676    0.444)   20.582  13.682   (  -4.590    4.205   -0.502)    6.245  14.427   (   2.881   25.878    4.413)   26.409  14.659   (   4.836   18.807   -3.165)   19.675======================= Grid point 357 (69/84) =======================q-point: ( 0.46  0.08  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.725   (   2.330   -1.229    2.684)    3.761   1.938   (   3.261   -2.322   -0.624)    4.051   2.146   (   1.712   12.435    0.138)   12.553   2.194   (  -0.238    2.071    0.607)    2.171   2.218   (  -0.020    3.546   -1.191)    3.740   2.337   (   0.155    2.959    0.192)    2.969   2.339   (  -1.082    7.473    1.316)    7.665   2.419   (   0.514    8.675   -4.940)    9.996   3.138   (   0.738   -2.356    2.897)    3.806   3.937   (   4.398    5.890   -2.921)    7.910   4.382   (   4.465   -1.409    0.850)    4.758   4.464   (   3.784    0.891   -0.460)    3.915   6.512   (  -7.821    9.078   -1.476)   12.073   8.059   (   4.909  -19.317    0.144)   19.932   8.508   (  -5.527  -12.097   -1.758)   13.415   9.886   (  -4.133  -13.683   -0.449)   14.301  10.612   (  14.048  -37.435   -0.559)   39.988  11.036   (   1.046   -9.434    1.772)    9.656  11.170   (  -2.214   -0.685   -0.241)    2.330  11.635   ( -10.431   -8.719    0.480)   13.603  11.704   (  -0.316  -38.284   -0.821)   38.294  13.633   (  -0.910    4.784   -0.358)    4.883  14.445   (  -0.060   28.743    1.680)   28.792  14.727   (   1.866   18.173    0.486)   18.275======================= Grid point 366 (70/84) =======================q-point: ( 0.15  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.951   (   8.443    8.443    2.043)   12.114   1.045   (   9.562    9.562    3.687)   14.016   1.925   (   8.615    8.615   -3.896)   12.792   1.930   (  11.780   11.780   -0.603)   16.670   2.172   (   6.763    6.763   -1.510)    9.683   2.281   (   7.122    7.122    0.670)   10.094   2.283   (   4.269    4.269   -0.550)    6.062   2.492   (   0.582    0.582    4.525)    4.599   2.695   (   6.356    6.356   -3.525)    9.655   3.187   (  14.045   14.045    0.394)   19.867   3.609   (   5.896    5.896   -0.558)    8.357   4.510   ( -10.042  -10.042    1.425)   14.274   6.762   (  -1.309   -1.309    1.386)    2.312   7.759   (  -8.621   -8.621   -0.027)   12.192   8.685   (   1.962    1.962    1.774)    3.294   9.043   (  -9.489   -9.489    0.553)   13.431   9.725   (  20.148   20.148   -7.775)   29.536  11.139   ( -12.084  -12.084    0.058)   17.089  11.452   ( -13.490  -13.490   -0.943)   19.101  12.306   ( -16.484  -16.484    9.587)   25.206  12.447   ( -11.523  -11.523    0.491)   16.303  14.032   (  24.092   24.092   15.798)   37.556  14.250   ( -11.021  -11.021   -0.087)   15.586  14.504   (  17.358   17.358  -24.353)   34.578======================= Grid point 367 (71/84) =======================q-point: ( 0.23  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.167   (  12.936    7.067    1.458)   14.813   1.308   (  15.681    4.095    2.065)   16.338   2.062   (   5.820    9.106   -2.344)   11.059   2.083   (   5.059    4.957   -1.564)    7.253   2.313   (   6.410    1.626   -0.248)    6.617   2.367   (   2.569    4.276    0.380)    5.003   2.402   (   4.496   11.754    1.062)   12.630   2.543   (   5.176    6.449   -0.532)    8.287   2.833   (   7.995    1.198   -1.034)    8.150   3.455   (  13.081   14.779   -0.129)   19.737   3.808   (  13.977    4.379   -0.610)   14.660   4.287   ( -11.370   -2.885    1.618)   11.841   6.828   (   7.316  -20.640    0.917)   21.918   7.569   ( -11.333   -5.394    0.688)   12.570   8.687   (  -5.459  -11.861    0.119)   13.057   8.949   (   4.203  -17.872    1.048)   18.390  10.021   (   8.761    6.749   -4.758)   12.039  10.899   ( -11.441  -14.977   -0.328)   18.850  11.232   (  -9.074  -22.665   -0.482)   24.419  11.919   ( -21.212   -8.247    5.596)   23.437  12.203   ( -13.873  -16.609    0.730)   21.653  14.024   ( -10.580   -5.131   -0.003)   11.759  14.437   (  17.556   21.289   13.066)   30.531  14.828   (  13.888   23.924  -17.528)   32.749======================= Grid point 368 (72/84) =======================q-point: ( 0.31  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.424   (  12.419    4.065    1.562)   13.160   1.608   (  14.025    3.293    0.482)   14.414   2.160   (   3.532    2.757   -1.568)    4.747   2.169   (   4.660    7.775   -1.255)    9.151   2.379   (   0.080    2.172    1.049)    2.413   2.393   (   1.903    4.762    0.226)    5.133   2.488   (   3.973   12.900    1.102)   13.543   2.622   (   2.048   11.852   -3.854)   12.630   2.989   (   6.481   -3.428    1.680)    7.521   3.712   (  13.163   13.016   -1.232)   18.553   4.101   (  13.646    3.145    0.200)   14.005   4.140   (  -0.973    0.332    0.867)    1.345   7.024   (  12.231  -33.770    0.302)   35.918   7.302   ( -15.696    3.872    0.746)   16.183   8.513   ( -10.232  -23.300   -0.748)   25.459   9.163   (  14.969  -27.052    0.515)   30.922  10.019   (  -9.116   -6.592   -2.348)   11.492  10.743   (  -3.841  -17.872   -0.438)   18.286  11.089   (  -5.819  -29.116   -0.107)   29.692  11.532   ( -17.368   -3.935    2.537)   17.988  11.905   ( -16.300  -20.028    0.625)   25.830  13.858   (  -6.118    1.488   -0.266)    6.302  14.733   (  12.635   21.551   10.015)   26.914  15.025   (   6.593   25.007  -11.115)   28.149======================= Grid point 369 (73/84) =======================q-point: ( 0.38  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.620   (   7.038    0.698    2.011)    7.353   1.834   (   8.831    3.838   -0.534)    9.644   2.210   (   1.723    2.259   -1.115)    3.052   2.253   (   3.610    6.867   -0.532)    7.776   2.385   (   0.588    3.139    0.559)    3.242   2.418   (   0.928    7.399   -1.282)    7.567   2.543   (   1.572   15.497    1.213)   15.624   2.638   (   0.172   12.703   -2.626)   12.973   3.060   (   1.533   -3.123    2.115)    4.072   3.949   (  10.435   10.827   -1.421)   15.104   4.204   (   7.440   -2.079    0.902)    7.778   4.358   (  10.347    3.063   -0.273)   10.794   6.987   ( -15.093   12.080   -0.285)   19.334   7.292   (  14.438  -37.116   -0.148)   39.826   8.320   (  -8.890  -22.941   -1.029)   24.624   9.475   (  15.627  -40.136   -0.081)   43.071   9.733   ( -17.473  -20.472   -1.065)   26.936  10.720   (   0.471  -20.635    0.308)   20.643  10.993   (  -4.225  -33.256   -0.087)   33.524  11.259   ( -10.615   -3.842    0.836)   11.320  11.590   ( -14.926  -21.284    0.461)   26.000  13.784   (  -1.739    6.686   -0.323)    6.916  14.922   (   6.542   22.901    8.613)   25.327  15.110   (   2.829   24.650   -7.885)   26.034======================= Grid point 370 (74/84) =======================q-point: ( 0.46  0.15  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.700   (   1.646   -1.172    2.159)    2.957   1.950   (   2.982    3.892   -0.983)    5.001   2.230   (   0.438    2.584   -0.718)    2.718   2.302   (   1.325    5.929   -0.731)    6.119   2.397   (   0.444    3.465    0.040)    3.494   2.431   (   0.291    9.884   -1.071)    9.946   2.557   (   0.167   17.764    1.217)   17.807   2.640   (   0.066   12.169   -2.078)   12.345   3.082   (   0.797   -2.734    2.034)    3.500   4.096   (   4.038    9.805   -1.219)   10.673   4.336   (   4.233   -3.308    0.373)    5.385   4.496   (   3.616    2.571   -0.132)    4.439   6.765   (  -6.421   16.279   -1.364)   17.553   7.534   (   8.032  -32.928   -0.320)   33.895   8.193   (  -3.585  -19.158   -1.257)   19.531   9.440   (  -9.624  -32.385   -0.406)   33.787   9.718   (   7.511  -50.909   -0.431)   51.462  10.738   (   0.876  -21.533    1.158)   21.582  10.927   (  -2.136  -36.340   -0.349)   36.404  11.110   (  -4.448   -5.542    0.421)    7.119  11.372   (  -6.091  -19.037    0.176)   19.989  13.771   (  -0.087    9.203   -0.283)    9.208  14.995   (   1.439   25.065    6.665)   25.976  15.150   (   1.217   23.276   -5.201)   23.881======================= Grid point 380 (75/84) =======================q-point: ( 0.23  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.336   (   9.643    9.643    1.358)   13.705   1.477   (  12.296   12.296    1.879)   17.490   2.149   (   3.108    3.108   -1.825)    4.760   2.218   (   6.153    6.153   -2.266)    8.991   2.359   (   3.477    3.477   -0.373)    4.931   2.404   (   0.794    0.794    0.606)    1.276   2.667   (   6.382    6.382    4.507)   10.088   2.671   (  12.316   12.316   -0.105)   17.418   2.870   (   2.698    2.698   -5.795)    6.938   3.745   (  13.644   13.644    0.422)   19.300   3.955   (  11.018   11.018   -0.540)   15.591   4.205   (  -5.119   -5.119    1.558)    7.404   6.487   ( -13.039  -13.039    0.699)   18.453   7.422   (  -8.064   -8.064    0.253)   11.407   8.519   ( -16.721  -16.721    0.006)   23.647   8.628   (  -3.726   -3.726    0.855)    5.338  10.069   (  -3.237   -3.237   -3.388)    5.696  10.570   ( -17.703  -17.703   -0.543)   25.042  10.856   ( -15.571  -15.571    0.244)   22.022  11.694   ( -14.398  -14.398    3.885)   20.730  11.845   ( -19.506  -19.506    0.486)   27.591  13.937   (  -3.131   -3.131   -0.185)    4.432  14.819   (  17.629   17.629    9.546)   26.696  15.222   (  15.661   15.661  -11.861)   25.123======================= Grid point 381 (76/84) =======================q-point: ( 0.31  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.519   (   8.383    5.328    1.500)   10.045   1.744   (  13.749   10.263    0.810)   17.176   2.213   (   3.318    3.270   -1.752)    4.977   2.313   (   4.044    6.617   -1.111)    7.834   2.414   (   1.175    1.587    0.633)    2.074   2.446   (   4.199    2.366   -1.068)    4.937   2.746   (   3.619    8.671    5.019)   10.653   2.858   (   3.759   12.092   -4.677)   13.499   2.964   (   7.673    3.195   -1.097)    8.383   3.977   (  10.561   12.953    0.028)   16.712   4.115   (  -1.408   -1.515    1.906)    2.813   4.218   (  12.977    7.575   -0.986)   15.059   6.358   (   0.596  -31.534    0.176)   31.540   7.297   (  -4.693   -5.261    0.393)    7.061   8.127   ( -18.981  -15.442   -0.637)   24.477   8.684   (   4.723  -19.595    0.428)   20.161   9.815   ( -23.420  -13.805   -1.959)   27.256  10.295   (  -8.817  -27.611   -0.641)   28.991  10.611   ( -10.060  -20.141    0.523)   22.520  11.385   ( -16.318  -11.578    2.024)   20.110  11.459   ( -19.264  -24.417    0.669)   31.109  13.913   (   0.534    4.294   -0.359)    4.342  15.118   (  12.760   17.460    7.497)   22.888  15.443   (   7.502   17.150   -8.157)   20.419======================= Grid point 382 (77/84) =======================q-point: ( 0.38  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.643   (   4.158    1.739    1.563)    4.770   1.973   (   9.039   10.089   -0.290)   13.549   2.269   (   2.147    3.793   -1.232)    4.529   2.381   (   2.782    6.045    0.214)    6.658   2.441   (   1.634    2.876   -0.442)    3.337   2.525   (   3.587    4.104   -1.090)    5.559   2.812   (   2.663    8.691    2.893)    9.539   2.893   (   1.273   10.590   -2.916)   11.058   3.079   (   3.209    9.496    0.394)   10.031   4.133   (   5.276    0.179    2.365)    5.784   4.186   (   8.482    7.065   -1.838)   11.191   4.445   (   9.298    5.898   -0.226)   11.013   6.518   (  14.966  -38.904   -0.265)   41.684   7.211   (  -5.023    8.443   -0.061)    9.824   7.851   (  -7.875  -23.443   -0.851)   24.745   8.749   (   1.156  -30.762   -0.116)   30.784   9.242   ( -32.453  -27.822   -0.934)   42.757  10.194   (  -2.311  -32.383   -0.228)   32.467  10.442   (  -7.595  -24.623    0.286)   25.770  11.079   ( -13.987  -26.217    0.417)   29.718  11.135   ( -13.113  -13.258    1.283)   18.692  13.944   (   2.180    9.362   -0.473)    9.625  15.315   (   7.310   17.194    6.831)   19.893  15.539   (   2.842   18.289   -6.470)   19.607======================= Grid point 383 (78/84) =======================q-point: ( 0.46  0.23  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.689   (   0.940    0.346    1.437)    1.752   2.090   (   2.928   10.052   -0.845)   10.504   2.295   (   0.600    3.979   -0.722)    4.088   2.416   (   0.844    5.103    1.241)    5.320   2.469   (   0.910    4.242   -1.803)    4.698   2.576   (   1.340    5.156   -0.774)    5.383   2.839   (   0.393    9.183    0.392)    9.199   2.920   (   1.018   12.057   -0.091)   12.100   3.110   (   0.592   10.711    0.665)   10.748   4.230   (   3.353   -3.105    2.410)    5.167   4.320   (   4.032    8.181   -2.458)    9.446   4.568   (   3.130    4.677    0.056)    5.628   6.847   (  13.813  -34.640   -0.640)   37.298   7.112   (  -3.593   17.345   -1.085)   17.746   7.804   (   0.746  -19.154   -0.364)   19.172   8.678   ( -19.512  -42.308    0.040)   46.591   8.735   (  -1.985  -46.585   -0.822)   46.634  10.191   (   1.459  -33.672    0.371)   33.706  10.315   (  -4.798  -28.127   -0.176)   28.534  10.891   (  -4.671  -30.041    0.300)   30.404  10.952   (  -5.220  -10.540    1.050)   11.808  13.980   (   1.102   11.648   -0.513)   11.711  15.406   (   2.225   17.313    6.683)   18.691  15.572   (   0.764   18.650   -5.865)   19.566======================= Grid point 394 (79/84) =======================q-point: ( 0.31  0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.618   (   4.513    4.513    1.308)    6.515   1.982   (  12.942   12.942    0.741)   18.318   2.286   (   3.887    3.887   -1.416)    5.677   2.421   (   4.392    4.392   -1.004)    6.292   2.445   (   1.823    1.823    0.457)    2.618   2.505   (   3.790    3.790   -0.917)    5.437   2.852   (   4.479    4.479    5.378)    8.310   2.979   (   3.373    3.373   -5.887)    7.578   3.179   (  13.010   13.010   -0.102)   18.399   4.076   (  -1.778   -1.778    1.547)    2.953   4.181   (   7.659    7.659    0.194)   10.833   4.403   (  10.489   10.489   -0.908)   14.862   5.905   ( -12.315  -12.315    0.139)   17.417   7.172   (  -4.121   -4.121   -0.136)    5.829   7.813   ( -18.172  -18.172   -0.302)   25.701   8.472   (  -3.552   -3.552    0.468)    5.044   9.468   ( -26.049  -26.049   -1.902)   36.888   9.727   ( -24.677  -24.677   -0.042)   34.899  10.269   ( -15.566  -15.566    0.413)   22.018  10.920   ( -27.467  -27.467    0.072)   38.844  11.154   ( -13.580  -13.580    2.059)   19.315  14.018   (   6.112    6.112   -0.584)    8.663  15.408   (  12.163   12.163    5.894)   18.183  15.697   (   9.127    9.127   -6.093)   14.273======================= Grid point 395 (80/84) =======================q-point: ( 0.38  0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.687   (   2.513    2.718    1.100)    3.862   2.202   (   8.629   12.097    0.144)   14.860   2.350   (   2.324    4.099   -0.648)    4.756   2.472   (   1.316    3.165    1.461)    3.726   2.507   (   3.837    3.675   -1.479)    5.515   2.581   (   3.421    2.259   -0.789)    4.175   2.943   (   4.172    6.346    3.873)    8.525   3.036   (   2.725    6.791   -3.872)    8.279   3.378   (   6.680   16.213   -0.108)   17.536   4.073   (   1.957   -3.563    0.850)    4.153   4.324   (   8.241    4.113   -0.267)    9.215   4.585   (   7.276    7.988   -0.421)   10.813   5.915   (  12.652  -19.149   -0.051)   22.952   7.219   (   9.558   -6.870   -0.371)   11.777   7.470   ( -14.969  -15.558   -0.242)   21.591   8.386   (  -8.935   -5.635    0.376)   10.571   8.725   ( -44.331  -23.821   -1.006)   50.336   9.465   (  -5.911  -41.523   -0.180)   41.942   9.988   ( -13.082  -22.235    0.055)   25.798  10.440   ( -19.483  -35.777    0.075)   40.738  10.894   ( -12.282  -13.410    1.792)   18.273  14.133   (   5.145    9.461   -0.825)   10.802  15.594   (   6.903   11.454    5.143)   14.328  15.819   (   3.908   10.628   -4.993)   12.375======================= Grid point 396 (81/84) =======================q-point: ( 0.46  0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.719   (   0.754    2.625    0.961)    2.896   2.310   (   2.682   11.074   -0.217)   11.396   2.376   (   0.515    4.023   -0.020)    4.056   2.486   (   0.256    2.276    2.464)    3.365   2.564   (   1.635    4.897   -2.690)    5.821   2.627   (   1.216    1.010   -0.405)    1.632   2.992   (   0.932    7.421    1.285)    7.589   3.086   (   1.797    8.971   -0.832)    9.187   3.448   (   1.270   18.373   -0.133)   18.418   4.129   (   2.394   -5.590    0.438)    6.097   4.483   (   5.890    6.769   -0.629)    8.995   4.671   (   1.696    5.638   -0.038)    5.887   6.295   (  19.832  -18.902   -0.662)   27.405   7.246   (  -5.568   -9.357   -0.291)   10.893   7.559   (  21.858   -6.368   -0.100)   22.767   7.919   ( -31.064  -32.961   -0.126)   45.293   8.065   ( -22.362  -13.542   -0.333)   26.145   9.426   (   0.562  -44.232   -0.105)   44.236   9.780   (  -6.917  -27.162   -0.243)   28.030  10.206   (  -4.934  -39.541    0.438)   39.850  10.709   (  -5.540  -13.818    1.579)   14.970  14.205   (   1.976   10.898   -0.974)   11.119  15.682   (   2.189   11.074    4.909)   12.310  15.865   (   1.033   11.433   -4.579)   12.359======================= Grid point 408 (82/84) =======================q-point: ( 0.38  0.38  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.743   (   2.690    2.690    0.896)    3.909   2.390   (   6.193    6.193    0.633)    8.782   2.414   (   2.323    2.323    0.137)    3.288   2.516   (   0.976    0.976    1.118)    1.777   2.574   (   3.668    3.668   -0.037)    5.187   2.632   (   2.530    2.530   -1.307)    3.808   3.070   (   6.131    6.131    3.080)    9.201   3.175   (   6.878    6.878   -3.062)   10.197   3.645   (  10.234   10.234   -0.485)   14.481   4.020   (  -1.245   -1.245    0.086)    1.763   4.408   (   7.220    7.220   -0.548)   10.225   4.725   (   5.468    5.468   -0.248)    7.738   5.822   (   8.617    8.617    0.171)   12.187   7.129   (  -2.655   -2.655   -0.239)    3.762   7.207   (  -6.744   -6.744   -0.352)    9.544   8.238   ( -32.846  -32.846   -1.165)   46.466   8.384   (  -1.236   -1.236    0.821)    1.931   8.709   ( -26.670  -26.670    0.223)   37.717   9.570   ( -19.949  -19.949   -0.021)   28.213   9.757   ( -30.777  -30.777   -0.032)   43.526  10.630   ( -12.752  -12.752    1.932)   18.137  14.296   (   6.831    6.831   -1.205)    9.735  15.767   (   6.324    6.324    4.131)    9.852  15.969   (   5.031    5.031   -4.088)    8.206======================= Grid point 409 (83/84) =======================q-point: ( 0.46  0.38  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 423Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.782   (   1.140    3.467    0.873)    3.752   2.440   (   0.603    2.567    0.752)    2.742   2.462   (   1.792    3.364    1.165)    3.985   2.511   (  -0.885    0.306    2.283)    2.468   2.640   (   2.514    3.022   -1.165)    4.100   2.662   (   0.303    2.656   -1.469)    3.051   3.156   (   2.293    9.155    1.932)    9.634   3.276   (   2.840    9.464   -1.766)   10.038   3.760   (   2.129   12.973   -0.706)   13.165   4.031   (   1.428   -4.092   -0.524)    4.366   4.604   (   9.581    6.805   -0.500)   11.762   4.772   (  -0.827    4.219   -0.140)    4.302   6.173   (  22.435    6.614   -0.453)   23.394   6.988   (  -4.963  -12.639    0.710)   13.597   7.383   (  13.413   -4.183   -1.082)   14.092   7.441   ( -35.332  -15.631   -0.718)   38.642   8.188   ( -16.243   14.299    1.205)   21.674   8.497   (  -1.891  -49.202   -0.606)   49.242   9.258   (  -9.860  -25.147    0.133)   27.012   9.381   (  -7.169  -43.190    0.570)   43.785  10.436   (  -5.833  -13.241    1.856)   14.587  14.389   (   2.468    7.544   -1.442)    8.067  15.847   (   1.992    5.965    3.680)    7.286  16.029   (   1.423    5.698   -3.679)    6.930======================= Grid point 422 (84/84) =======================q-point: ( 0.46  0.46  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 3.36e-05 3.36e-05 3.36e-05 Number of triplets: 231Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.836   (   1.587    1.587    0.954)    2.439   2.472   (   0.773    0.773    1.013)    1.490   2.484   (   0.375    0.375    2.814)    2.863   2.501   (  -1.043   -1.043    2.590)    2.981   2.691   (   1.775    1.775   -1.645)    3.001   2.700   (   0.908    0.908   -1.769)    2.186   3.300   (   4.249    4.249    3.213)    6.814   3.417   (   4.091    4.091   -3.068)    6.549   3.934   (   4.033    4.033   -1.170)    5.823   3.981   (  -0.619   -0.619   -1.104)    1.409   4.748   (   6.492    6.492   -0.364)    9.188   4.818   (  -0.089   -0.089   -0.186)    0.224   6.518   (  25.860   25.860   -0.606)   36.577   6.810   (  -4.815   -4.815    1.062)    6.892   7.202   ( -15.641  -15.641   -0.774)   22.134   7.496   (  11.607   11.607   -2.092)   16.548   7.643   ( -28.538  -28.538   -0.904)   40.369   8.387   (   0.570    0.570    1.965)    2.124   8.653   ( -22.423  -22.423    1.092)   31.730   8.870   ( -11.940  -11.940    0.717)   16.901  10.235   (  -6.082   -6.082    1.843)    8.797  14.490   (   2.675    2.675   -1.725)    4.158  15.922   (   1.854    1.854    3.084)    4.048  16.099   (   1.672    1.672   -3.251)    4.020=================== 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   1047.582   1047.582    393.435     -0.000     -0.000      0.000 3/20280   20.0    199.414    199.414     44.862     -0.000      0.000      0.000 3/20280   30.0     96.535     96.535     16.516     -0.000     -0.000      0.000 3/20280   40.0     62.579     62.579      9.789     -0.000     -0.000      0.000 3/20280   50.0     45.945     45.945      6.976     -0.000     -0.000      0.000 3/20280   60.0     35.955     35.955      5.416     -0.000     -0.000      0.000 3/20280   70.0     29.258     29.258      4.409     -0.000     -0.000      0.000 3/20280   80.0     24.481     24.481      3.704     -0.000     -0.000      0.000 3/20280   90.0     20.939     20.939      3.184     -0.000     -0.000      0.000 3/20280  100.0     18.235     18.235      2.785     -0.000     -0.000      0.000 3/20280  110.0     16.121     16.121      2.473     -0.000     -0.000      0.000 3/20280  120.0     14.434     14.434      2.221     -0.000     -0.000      0.000 3/20280  130.0     13.062     13.062      2.016     -0.000     -0.000      0.000 3/20280  140.0     11.928     11.928      1.845     -0.000     -0.000      0.000 3/20280  150.0     10.977     10.977      1.701     -0.000     -0.000      0.000 3/20280  160.0     10.169     10.169      1.579     -0.000     -0.000      0.000 3/20280  170.0      9.474      9.474      1.473     -0.000     -0.000      0.000 3/20280  180.0      8.871      8.871      1.381     -0.000     -0.000      0.000 3/20280  190.0      8.343      8.343      1.300     -0.000     -0.000      0.000 3/20280  200.0      7.875      7.875      1.228     -0.000     -0.000      0.000 3/20280  210.0      7.459      7.459      1.164     -0.000     -0.000      0.000 3/20280  220.0      7.086      7.086      1.107     -0.000     -0.000      0.000 3/20280  230.0      6.750      6.750      1.055     -0.000     -0.000      0.000 3/20280  240.0      6.446      6.446      1.008     -0.000     -0.000      0.000 3/20280  250.0      6.168      6.168      0.965     -0.000     -0.000      0.000 3/20280  260.0      5.914      5.914      0.925     -0.000     -0.000      0.000 3/20280  270.0      5.681      5.681      0.889     -0.000     -0.000      0.000 3/20280  280.0      5.466      5.466      0.856     -0.000     -0.000      0.000 3/20280  290.0      5.267      5.267      0.825     -0.000     -0.000      0.000 3/20280  300.0      5.082      5.082      0.796     -0.000     -0.000      0.000 3/20280  310.0      4.911      4.911      0.770     -0.000     -0.000      0.000 3/20280  320.0      4.750      4.750      0.745     -0.000     -0.000      0.000 3/20280  330.0      4.600      4.600      0.721     -0.000     -0.000      0.000 3/20280  340.0      4.460      4.460      0.699     -0.000     -0.000      0.000 3/20280  350.0      4.328      4.328      0.679     -0.000     -0.000      0.000 3/20280  360.0      4.203      4.203      0.659     -0.000     -0.000      0.000 3/20280  370.0      4.086      4.086      0.641     -0.000     -0.000      0.000 3/20280  380.0      3.975      3.975      0.624     -0.000     -0.000      0.000 3/20280  390.0      3.871      3.871      0.608     -0.000     -0.000      0.000 3/20280  400.0      3.771      3.771      0.592     -0.000     -0.000      0.000 3/20280  410.0      3.677      3.677      0.577     -0.000     -0.000      0.000 3/20280  420.0      3.587      3.587      0.563     -0.000     -0.000      0.000 3/20280  430.0      3.502      3.502      0.550     -0.000     -0.000      0.000 3/20280  440.0      3.421      3.421      0.537     -0.000     -0.000      0.000 3/20280  450.0      3.343      3.343      0.525     -0.000     -0.000      0.000 3/20280  460.0      3.269      3.269      0.514     -0.000     -0.000      0.000 3/20280  470.0      3.198      3.198      0.502     -0.000     -0.000      0.000 3/20280  480.0      3.130      3.130      0.492     -0.000     -0.000      0.000 3/20280  490.0      3.066      3.066      0.482     -0.000     -0.000      0.000 3/20280  500.0      3.003      3.003      0.472     -0.000     -0.000      0.000 3/20280  510.0      2.943      2.943      0.463     -0.000     -0.000      0.000 3/20280  520.0      2.886      2.886      0.454     -0.000     -0.000      0.000 3/20280  530.0      2.831      2.831      0.445     -0.000     -0.000      0.000 3/20280  540.0      2.778      2.778      0.437     -0.000     -0.000      0.000 3/20280  550.0      2.727      2.727      0.429     -0.000     -0.000      0.000 3/20280  560.0      2.677      2.677      0.421     -0.000     -0.000      0.000 3/20280  570.0      2.630      2.630      0.413     -0.000     -0.000      0.000 3/20280  580.0      2.584      2.584      0.406     -0.000     -0.000      0.000 3/20280  590.0      2.540      2.540      0.399     -0.000     -0.000      0.000 3/20280  600.0      2.497      2.497      0.393     -0.000     -0.000      0.000 3/20280  610.0      2.456      2.456      0.386     -0.000     -0.000      0.000 3/20280  620.0      2.416      2.416      0.380     -0.000     -0.000      0.000 3/20280  630.0      2.377      2.377      0.374     -0.000     -0.000      0.000 3/20280  640.0      2.340      2.340      0.368     -0.000     -0.000      0.000 3/20280  650.0      2.303      2.303      0.362     -0.000     -0.000      0.000 3/20280  660.0      2.268      2.268      0.357     -0.000     -0.000      0.000 3/20280  670.0      2.234      2.234      0.351     -0.000     -0.000      0.000 3/20280  680.0      2.201      2.201      0.346     -0.000     -0.000      0.000 3/20280  690.0      2.169      2.169      0.341     -0.000     -0.000      0.000 3/20280  700.0      2.138      2.138      0.336     -0.000     -0.000      0.000 3/20280  710.0      2.107      2.107      0.331     -0.000     -0.000      0.000 3/20280  720.0      2.078      2.078      0.327     -0.000     -0.000      0.000 3/20280  730.0      2.049      2.049      0.322     -0.000     -0.000      0.000 3/20280  740.0      2.021      2.021      0.318     -0.000     -0.000      0.000 3/20280  750.0      1.994      1.994      0.314     -0.000     -0.000      0.000 3/20280  760.0      1.968      1.968      0.310     -0.000     -0.000      0.000 3/20280  770.0      1.942      1.942      0.306     -0.000     -0.000      0.000 3/20280  780.0      1.917      1.917      0.302     -0.000     -0.000      0.000 3/20280  790.0      1.893      1.893      0.298     -0.000     -0.000      0.000 3/20280  800.0      1.869      1.869      0.294     -0.000     -0.000      0.000 3/20280  810.0      1.846      1.846      0.290     -0.000     -0.000      0.000 3/20280  820.0      1.823      1.823      0.287     -0.000     -0.000      0.000 3/20280  830.0      1.801      1.801      0.283     -0.000     -0.000      0.000 3/20280  840.0      1.780      1.780      0.280     -0.000     -0.000      0.000 3/20280  850.0      1.758      1.758      0.277     -0.000     -0.000      0.000 3/20280  860.0      1.738      1.738      0.273     -0.000     -0.000      0.000 3/20280  870.0      1.718      1.718      0.270     -0.000     -0.000      0.000 3/20280  880.0      1.698      1.698      0.267     -0.000     -0.000      0.000 3/20280  890.0      1.679      1.679      0.264     -0.000     -0.000      0.000 3/20280  900.0      1.660      1.660      0.261     -0.000     -0.000      0.000 3/20280  910.0      1.642      1.642      0.258     -0.000     -0.000      0.000 3/20280  920.0      1.624      1.624      0.256     -0.000     -0.000      0.000 3/20280  930.0      1.607      1.607      0.253     -0.000     -0.000      0.000 3/20280  940.0      1.590      1.590      0.250     -0.000     -0.000      0.000 3/20280  950.0      1.573      1.573      0.248     -0.000     -0.000      0.000 3/20280  960.0      1.556      1.556      0.245     -0.000     -0.000      0.000 3/20280  970.0      1.540      1.540      0.242     -0.000     -0.000      0.000 3/20280  980.0      1.524      1.524      0.240     -0.000     -0.000      0.000 3/20280  990.0      1.509      1.509      0.238     -0.000     -0.000      0.000 3/20280 1000.0      1.494      1.494      0.235     -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-09 04:15:42]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|