# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/e8b5096d-00d9-46b6-959a-7152582b7e32/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 BaLiAs / P-6m2 (187) / materials id 10616](https://mdr.nims.go.jp/datasets/d0c69f06-691e-45a4-a968-2234e735898d)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-07 23:44:33]-------------------------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: [4 4 3]  Primitive matrix:    [1. 0. 0.]    [0. 1. 0.]    [0. 0. 1.]Spacegroup: P-6m2 (187)Number of symmetry operations in supercell: 576------------------------------ primitive cell ------------------------------Lattice vectors:  a    4.488603559793336    0.000000000000000    0.000000000000000  b   -2.244301779896668    3.887244710298293    0.000000000000000  c    0.000000000000000    0.000000000000000    4.539274940000000Atomic positions (fractional):   *1 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941   *2 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922   *3 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327-------------------------------- unit cell ---------------------------------Lattice vectors:  a    4.488603559793336    0.000000000000000    0.000000000000000  b   -2.244301779896668    3.887244710298293    0.000000000000000  c    0.000000000000000    0.000000000000000    4.539274940000000Atomic positions (fractional):   *1 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1   *2 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922 > 2   *3 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327 > 3-------------------------------- super cell --------------------------------Lattice vectors:  a   17.954414239173342    0.000000000000000    0.000000000000000  b   -8.977207119586671   15.548978841193172    0.000000000000000  c    0.000000000000000    0.000000000000000   13.617824820000001Atomic positions (fractional):   *1 Li  0.00000000000000  0.00000000000000  0.16666666666667   6.941 > 1    2 Li  0.25000000000000  0.00000000000000  0.16666666666667   6.941 > 1    3 Li  0.50000000000000  0.00000000000000  0.16666666666667   6.941 > 1    4 Li  0.75000000000000  0.00000000000000  0.16666666666667   6.941 > 1    5 Li  0.00000000000000  0.25000000000000  0.16666666666667   6.941 > 1    6 Li  0.25000000000000  0.25000000000000  0.16666666666667   6.941 > 1    7 Li  0.50000000000000  0.25000000000000  0.16666666666667   6.941 > 1    8 Li  0.75000000000000  0.25000000000000  0.16666666666667   6.941 > 1    9 Li  0.00000000000000  0.50000000000000  0.16666666666667   6.941 > 1   10 Li  0.25000000000000  0.50000000000000  0.16666666666667   6.941 > 1   11 Li  0.50000000000000  0.50000000000000  0.16666666666667   6.941 > 1   12 Li  0.75000000000000  0.50000000000000  0.16666666666667   6.941 > 1   13 Li  0.00000000000000  0.75000000000000  0.16666666666667   6.941 > 1   14 Li  0.25000000000000  0.75000000000000  0.16666666666667   6.941 > 1   15 Li  0.50000000000000  0.75000000000000  0.16666666666667   6.941 > 1   16 Li  0.75000000000000  0.75000000000000  0.16666666666667   6.941 > 1   17 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1   18 Li  0.25000000000000  0.00000000000000  0.50000000000000   6.941 > 1   19 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1   20 Li  0.75000000000000  0.00000000000000  0.50000000000000   6.941 > 1   21 Li  0.00000000000000  0.25000000000000  0.50000000000000   6.941 > 1   22 Li  0.25000000000000  0.25000000000000  0.50000000000000   6.941 > 1   23 Li  0.50000000000000  0.25000000000000  0.50000000000000   6.941 > 1   24 Li  0.75000000000000  0.25000000000000  0.50000000000000   6.941 > 1   25 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1   26 Li  0.25000000000000  0.50000000000000  0.50000000000000   6.941 > 1   27 Li  0.50000000000000  0.50000000000000  0.50000000000000   6.941 > 1   28 Li  0.75000000000000  0.50000000000000  0.50000000000000   6.941 > 1   29 Li  0.00000000000000  0.75000000000000  0.50000000000000   6.941 > 1   30 Li  0.25000000000000  0.75000000000000  0.50000000000000   6.941 > 1   31 Li  0.50000000000000  0.75000000000000  0.50000000000000   6.941 > 1   32 Li  0.75000000000000  0.75000000000000  0.50000000000000   6.941 > 1   33 Li  0.00000000000000  0.00000000000000  0.83333333333333   6.941 > 1   34 Li  0.25000000000000  0.00000000000000  0.83333333333333   6.941 > 1   35 Li  0.50000000000000  0.00000000000000  0.83333333333333   6.941 > 1   36 Li  0.75000000000000  0.00000000000000  0.83333333333333   6.941 > 1   37 Li  0.00000000000000  0.25000000000000  0.83333333333333   6.941 > 1   38 Li  0.25000000000000  0.25000000000000  0.83333333333333   6.941 > 1   39 Li  0.50000000000000  0.25000000000000  0.83333333333333   6.941 > 1   40 Li  0.75000000000000  0.25000000000000  0.83333333333333   6.941 > 1   41 Li  0.00000000000000  0.50000000000000  0.83333333333333   6.941 > 1   42 Li  0.25000000000000  0.50000000000000  0.83333333333333   6.941 > 1   43 Li  0.50000000000000  0.50000000000000  0.83333333333333   6.941 > 1   44 Li  0.75000000000000  0.50000000000000  0.83333333333333   6.941 > 1   45 Li  0.00000000000000  0.75000000000000  0.83333333333333   6.941 > 1   46 Li  0.25000000000000  0.75000000000000  0.83333333333333   6.941 > 1   47 Li  0.50000000000000  0.75000000000000  0.83333333333333   6.941 > 1   48 Li  0.75000000000000  0.75000000000000  0.83333333333333   6.941 > 1  *49 As  0.08333333333333  0.16666666666667  0.16666666666667  74.922 > 2   50 As  0.33333333333333  0.16666666666667  0.16666666666667  74.922 > 2   51 As  0.58333333333333  0.16666666666667  0.16666666666667  74.922 > 2   52 As  0.83333333333333  0.16666666666667  0.16666666666667  74.922 > 2   53 As  0.08333333333333  0.41666666666667  0.16666666666667  74.922 > 2   54 As  0.33333333333333  0.41666666666667  0.16666666666667  74.922 > 2   55 As  0.58333333333333  0.41666666666667  0.16666666666667  74.922 > 2   56 As  0.83333333333333  0.41666666666667  0.16666666666667  74.922 > 2   57 As  0.08333333333333  0.66666666666667  0.16666666666667  74.922 > 2   58 As  0.33333333333333  0.66666666666667  0.16666666666667  74.922 > 2   59 As  0.58333333333333  0.66666666666667  0.16666666666667  74.922 > 2   60 As  0.83333333333333  0.66666666666667  0.16666666666667  74.922 > 2   61 As  0.08333333333333  0.91666666666667  0.16666666666667  74.922 > 2   62 As  0.33333333333333  0.91666666666667  0.16666666666667  74.922 > 2   63 As  0.58333333333333  0.91666666666667  0.16666666666667  74.922 > 2   64 As  0.83333333333333  0.91666666666667  0.16666666666667  74.922 > 2   65 As  0.08333333333333  0.16666666666667  0.50000000000000  74.922 > 2   66 As  0.33333333333333  0.16666666666667  0.50000000000000  74.922 > 2   67 As  0.58333333333333  0.16666666666667  0.50000000000000  74.922 > 2   68 As  0.83333333333333  0.16666666666667  0.50000000000000  74.922 > 2   69 As  0.08333333333333  0.41666666666667  0.50000000000000  74.922 > 2   70 As  0.33333333333333  0.41666666666667  0.50000000000000  74.922 > 2   71 As  0.58333333333333  0.41666666666667  0.50000000000000  74.922 > 2   72 As  0.83333333333333  0.41666666666667  0.50000000000000  74.922 > 2   73 As  0.08333333333333  0.66666666666667  0.50000000000000  74.922 > 2   74 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922 > 2   75 As  0.58333333333333  0.66666666666667  0.50000000000000  74.922 > 2   76 As  0.83333333333333  0.66666666666667  0.50000000000000  74.922 > 2   77 As  0.08333333333333  0.91666666666667  0.50000000000000  74.922 > 2   78 As  0.33333333333333  0.91666666666667  0.50000000000000  74.922 > 2   79 As  0.58333333333333  0.91666666666667  0.50000000000000  74.922 > 2   80 As  0.83333333333333  0.91666666666667  0.50000000000000  74.922 > 2   81 As  0.08333333333333  0.16666666666667  0.83333333333333  74.922 > 2   82 As  0.33333333333333  0.16666666666667  0.83333333333333  74.922 > 2   83 As  0.58333333333333  0.16666666666667  0.83333333333333  74.922 > 2   84 As  0.83333333333333  0.16666666666667  0.83333333333333  74.922 > 2   85 As  0.08333333333333  0.41666666666667  0.83333333333333  74.922 > 2   86 As  0.33333333333333  0.41666666666667  0.83333333333333  74.922 > 2   87 As  0.58333333333333  0.41666666666667  0.83333333333333  74.922 > 2   88 As  0.83333333333333  0.41666666666667  0.83333333333333  74.922 > 2   89 As  0.08333333333333  0.66666666666667  0.83333333333333  74.922 > 2   90 As  0.33333333333333  0.66666666666667  0.83333333333333  74.922 > 2   91 As  0.58333333333333  0.66666666666667  0.83333333333333  74.922 > 2   92 As  0.83333333333333  0.66666666666667  0.83333333333333  74.922 > 2   93 As  0.08333333333333  0.91666666666667  0.83333333333333  74.922 > 2   94 As  0.33333333333333  0.91666666666667  0.83333333333333  74.922 > 2   95 As  0.58333333333333  0.91666666666667  0.83333333333333  74.922 > 2   96 As  0.83333333333333  0.91666666666667  0.83333333333333  74.922 > 2  *97 Ba  0.16666666666667  0.08333333333333  0.00000000000000 137.327 > 3   98 Ba  0.41666666666667  0.08333333333333  0.00000000000000 137.327 > 3   99 Ba  0.66666666666667  0.08333333333333  0.00000000000000 137.327 > 3  100 Ba  0.91666666666667  0.08333333333333  0.00000000000000 137.327 > 3  101 Ba  0.16666666666667  0.33333333333333  0.00000000000000 137.327 > 3  102 Ba  0.41666666666667  0.33333333333333  0.00000000000000 137.327 > 3  103 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327 > 3  104 Ba  0.91666666666667  0.33333333333333  0.00000000000000 137.327 > 3  105 Ba  0.16666666666667  0.58333333333333  0.00000000000000 137.327 > 3  106 Ba  0.41666666666667  0.58333333333333  0.00000000000000 137.327 > 3  107 Ba  0.66666666666667  0.58333333333333  0.00000000000000 137.327 > 3  108 Ba  0.91666666666667  0.58333333333333  0.00000000000000 137.327 > 3  109 Ba  0.16666666666667  0.83333333333333  0.00000000000000 137.327 > 3  110 Ba  0.41666666666667  0.83333333333333  0.00000000000000 137.327 > 3  111 Ba  0.66666666666667  0.83333333333333  0.00000000000000 137.327 > 3  112 Ba  0.91666666666667  0.83333333333333  0.00000000000000 137.327 > 3  113 Ba  0.16666666666667  0.08333333333333  0.33333333333333 137.327 > 3  114 Ba  0.41666666666667  0.08333333333333  0.33333333333333 137.327 > 3  115 Ba  0.66666666666667  0.08333333333333  0.33333333333333 137.327 > 3  116 Ba  0.91666666666667  0.08333333333333  0.33333333333333 137.327 > 3  117 Ba  0.16666666666667  0.33333333333333  0.33333333333333 137.327 > 3  118 Ba  0.41666666666667  0.33333333333333  0.33333333333333 137.327 > 3  119 Ba  0.66666666666667  0.33333333333333  0.33333333333333 137.327 > 3  120 Ba  0.91666666666667  0.33333333333333  0.33333333333333 137.327 > 3  121 Ba  0.16666666666667  0.58333333333333  0.33333333333333 137.327 > 3  122 Ba  0.41666666666667  0.58333333333333  0.33333333333333 137.327 > 3  123 Ba  0.66666666666667  0.58333333333333  0.33333333333333 137.327 > 3  124 Ba  0.91666666666667  0.58333333333333  0.33333333333333 137.327 > 3  125 Ba  0.16666666666667  0.83333333333333  0.33333333333333 137.327 > 3  126 Ba  0.41666666666667  0.83333333333333  0.33333333333333 137.327 > 3  127 Ba  0.66666666666667  0.83333333333333  0.33333333333333 137.327 > 3  128 Ba  0.91666666666667  0.83333333333333  0.33333333333333 137.327 > 3  129 Ba  0.16666666666667  0.08333333333333  0.66666666666667 137.327 > 3  130 Ba  0.41666666666667  0.08333333333333  0.66666666666667 137.327 > 3  131 Ba  0.66666666666667  0.08333333333333  0.66666666666667 137.327 > 3  132 Ba  0.91666666666667  0.08333333333333  0.66666666666667 137.327 > 3  133 Ba  0.16666666666667  0.33333333333333  0.66666666666667 137.327 > 3  134 Ba  0.41666666666667  0.33333333333333  0.66666666666667 137.327 > 3  135 Ba  0.66666666666667  0.33333333333333  0.66666666666667 137.327 > 3  136 Ba  0.91666666666667  0.33333333333333  0.66666666666667 137.327 > 3  137 Ba  0.16666666666667  0.58333333333333  0.66666666666667 137.327 > 3  138 Ba  0.41666666666667  0.58333333333333  0.66666666666667 137.327 > 3  139 Ba  0.66666666666667  0.58333333333333  0.66666666666667 137.327 > 3  140 Ba  0.91666666666667  0.58333333333333  0.66666666666667 137.327 > 3  141 Ba  0.16666666666667  0.83333333333333  0.66666666666667 137.327 > 3  142 Ba  0.41666666666667  0.83333333333333  0.66666666666667 137.327 > 3  143 Ba  0.66666666666667  0.83333333333333  0.66666666666667 137.327 > 3  144 Ba  0.91666666666667  0.83333333333333  0.66666666666667 137.327 > 3----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            9.6531051   -0.0000000    0.0000000            0.0000000    9.6531051    0.0000000            0.0000000    0.0000000   45.3917432-------------------------- Born effective charges --------------------------    1 Li    0.9957698    0.0000000    0.0000000            0.0000000    0.9957698    0.0000000            0.0000000    0.0000000   -0.0040015    2 As   -3.1893956   -0.0000000    0.0000000           -0.0000000   -3.1893956    0.0000000            0.0000000    0.0000000   -5.8823552    3 Ba    2.1936258    0.0000000    0.0000000            0.0000000    2.1936258    0.0000000            0.0000000    0.0000000    5.8863567----------------------------------------------------------------------------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: 2574/2574Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 213Number of blocks in projector: 147Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 114Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 54Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 45Use standard eigh solver.Tree of FC basis block matrices:- (213, 207), data: False|-- (45, 45), data: True|-- (54, 51), data: True|-- (114, 111), data: True-----Solver_atoms: 1 -- 144 / 144Time (Solver_compr_matrix_reshape): 0.009Solver_block: 80 / 80 - Time: 0.192Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.223--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 432/432Permutation basis: 2574/2574Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 213Number of blocks in projector: 147Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 114Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 54Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 45Use standard eigh solver.Tree of FC basis block matrices:- (213, 207), data: False|-- (45, 45), data: True|-- (54, 51), data: True|-- (114, 111), data: TrueMax drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-07 23:44:39]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-07 23:44:39]-------------------------Compiled with OpenMP support (max 128 threads).Running in phono3py.load mode.Python version 3.14.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: force constantsHDF5 data compression filter: gzipCrystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".Supercell (dim): [4 4 3]Primitive matrix:  [1. 0. 0.]  [0. 1. 0.]  [0. 0. 1.]Spacegroup: P-6m2 (187)------------------------------ primitive cell ------------------------------Lattice vectors:  a    4.488603559793336    0.000000000000000    0.000000000000000  b   -2.244301779896668    3.887244710298293    0.000000000000000  c    0.000000000000000    0.000000000000000    4.539274940000000Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941    2 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922    3 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327-------------------------------- supercell ---------------------------------Lattice vectors:  a   17.954414239173342    0.000000000000000    0.000000000000000  b   -8.977207119586671   15.548978841193172    0.000000000000000  c    0.000000000000000    0.000000000000000   13.617824820000001Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.16666666666667   6.941 > 1    2 Li  0.25000000000000  0.00000000000000  0.16666666666667   6.941 > 1    3 Li  0.50000000000000  0.00000000000000  0.16666666666667   6.941 > 1    4 Li  0.75000000000000  0.00000000000000  0.16666666666667   6.941 > 1    5 Li  0.00000000000000  0.25000000000000  0.16666666666667   6.941 > 1    6 Li  0.25000000000000  0.25000000000000  0.16666666666667   6.941 > 1    7 Li  0.50000000000000  0.25000000000000  0.16666666666667   6.941 > 1    8 Li  0.75000000000000  0.25000000000000  0.16666666666667   6.941 > 1    9 Li  0.00000000000000  0.50000000000000  0.16666666666667   6.941 > 1   10 Li  0.25000000000000  0.50000000000000  0.16666666666667   6.941 > 1   11 Li  0.50000000000000  0.50000000000000  0.16666666666667   6.941 > 1   12 Li  0.75000000000000  0.50000000000000  0.16666666666667   6.941 > 1   13 Li  0.00000000000000  0.75000000000000  0.16666666666667   6.941 > 1   14 Li  0.25000000000000  0.75000000000000  0.16666666666667   6.941 > 1   15 Li  0.50000000000000  0.75000000000000  0.16666666666667   6.941 > 1   16 Li  0.75000000000000  0.75000000000000  0.16666666666667   6.941 > 1   17 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1   18 Li  0.25000000000000  0.00000000000000  0.50000000000000   6.941 > 1   19 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1   20 Li  0.75000000000000  0.00000000000000  0.50000000000000   6.941 > 1   21 Li  0.00000000000000  0.25000000000000  0.50000000000000   6.941 > 1   22 Li  0.25000000000000  0.25000000000000  0.50000000000000   6.941 > 1   23 Li  0.50000000000000  0.25000000000000  0.50000000000000   6.941 > 1   24 Li  0.75000000000000  0.25000000000000  0.50000000000000   6.941 > 1   25 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1   26 Li  0.25000000000000  0.50000000000000  0.50000000000000   6.941 > 1   27 Li  0.50000000000000  0.50000000000000  0.50000000000000   6.941 > 1   28 Li  0.75000000000000  0.50000000000000  0.50000000000000   6.941 > 1   29 Li  0.00000000000000  0.75000000000000  0.50000000000000   6.941 > 1   30 Li  0.25000000000000  0.75000000000000  0.50000000000000   6.941 > 1   31 Li  0.50000000000000  0.75000000000000  0.50000000000000   6.941 > 1   32 Li  0.75000000000000  0.75000000000000  0.50000000000000   6.941 > 1   33 Li  0.00000000000000  0.00000000000000  0.83333333333333   6.941 > 1   34 Li  0.25000000000000  0.00000000000000  0.83333333333333   6.941 > 1   35 Li  0.50000000000000  0.00000000000000  0.83333333333333   6.941 > 1   36 Li  0.75000000000000  0.00000000000000  0.83333333333333   6.941 > 1   37 Li  0.00000000000000  0.25000000000000  0.83333333333333   6.941 > 1   38 Li  0.25000000000000  0.25000000000000  0.83333333333333   6.941 > 1   39 Li  0.50000000000000  0.25000000000000  0.83333333333333   6.941 > 1   40 Li  0.75000000000000  0.25000000000000  0.83333333333333   6.941 > 1   41 Li  0.00000000000000  0.50000000000000  0.83333333333333   6.941 > 1   42 Li  0.25000000000000  0.50000000000000  0.83333333333333   6.941 > 1   43 Li  0.50000000000000  0.50000000000000  0.83333333333333   6.941 > 1   44 Li  0.75000000000000  0.50000000000000  0.83333333333333   6.941 > 1   45 Li  0.00000000000000  0.75000000000000  0.83333333333333   6.941 > 1   46 Li  0.25000000000000  0.75000000000000  0.83333333333333   6.941 > 1   47 Li  0.50000000000000  0.75000000000000  0.83333333333333   6.941 > 1   48 Li  0.75000000000000  0.75000000000000  0.83333333333333   6.941 > 1   49 As  0.08333333333333  0.16666666666667  0.16666666666667  74.922 > 49   50 As  0.33333333333333  0.16666666666667  0.16666666666667  74.922 > 49   51 As  0.58333333333333  0.16666666666667  0.16666666666667  74.922 > 49   52 As  0.83333333333333  0.16666666666667  0.16666666666667  74.922 > 49   53 As  0.08333333333333  0.41666666666667  0.16666666666667  74.922 > 49   54 As  0.33333333333333  0.41666666666667  0.16666666666667  74.922 > 49   55 As  0.58333333333333  0.41666666666667  0.16666666666667  74.922 > 49   56 As  0.83333333333333  0.41666666666667  0.16666666666667  74.922 > 49   57 As  0.08333333333333  0.66666666666667  0.16666666666667  74.922 > 49   58 As  0.33333333333333  0.66666666666667  0.16666666666667  74.922 > 49   59 As  0.58333333333333  0.66666666666667  0.16666666666667  74.922 > 49   60 As  0.83333333333333  0.66666666666667  0.16666666666667  74.922 > 49   61 As  0.08333333333333  0.91666666666667  0.16666666666667  74.922 > 49   62 As  0.33333333333333  0.91666666666667  0.16666666666667  74.922 > 49   63 As  0.58333333333333  0.91666666666667  0.16666666666667  74.922 > 49   64 As  0.83333333333333  0.91666666666667  0.16666666666667  74.922 > 49   65 As  0.08333333333333  0.16666666666667  0.50000000000000  74.922 > 49   66 As  0.33333333333333  0.16666666666667  0.50000000000000  74.922 > 49   67 As  0.58333333333333  0.16666666666667  0.50000000000000  74.922 > 49   68 As  0.83333333333333  0.16666666666667  0.50000000000000  74.922 > 49   69 As  0.08333333333333  0.41666666666667  0.50000000000000  74.922 > 49   70 As  0.33333333333333  0.41666666666667  0.50000000000000  74.922 > 49   71 As  0.58333333333333  0.41666666666667  0.50000000000000  74.922 > 49   72 As  0.83333333333333  0.41666666666667  0.50000000000000  74.922 > 49   73 As  0.08333333333333  0.66666666666667  0.50000000000000  74.922 > 49   74 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922 > 49   75 As  0.58333333333333  0.66666666666667  0.50000000000000  74.922 > 49   76 As  0.83333333333333  0.66666666666667  0.50000000000000  74.922 > 49   77 As  0.08333333333333  0.91666666666667  0.50000000000000  74.922 > 49   78 As  0.33333333333333  0.91666666666667  0.50000000000000  74.922 > 49   79 As  0.58333333333333  0.91666666666667  0.50000000000000  74.922 > 49   80 As  0.83333333333333  0.91666666666667  0.50000000000000  74.922 > 49   81 As  0.08333333333333  0.16666666666667  0.83333333333333  74.922 > 49   82 As  0.33333333333333  0.16666666666667  0.83333333333333  74.922 > 49   83 As  0.58333333333333  0.16666666666667  0.83333333333333  74.922 > 49   84 As  0.83333333333333  0.16666666666667  0.83333333333333  74.922 > 49   85 As  0.08333333333333  0.41666666666667  0.83333333333333  74.922 > 49   86 As  0.33333333333333  0.41666666666667  0.83333333333333  74.922 > 49   87 As  0.58333333333333  0.41666666666667  0.83333333333333  74.922 > 49   88 As  0.83333333333333  0.41666666666667  0.83333333333333  74.922 > 49   89 As  0.08333333333333  0.66666666666667  0.83333333333333  74.922 > 49   90 As  0.33333333333333  0.66666666666667  0.83333333333333  74.922 > 49   91 As  0.58333333333333  0.66666666666667  0.83333333333333  74.922 > 49   92 As  0.83333333333333  0.66666666666667  0.83333333333333  74.922 > 49   93 As  0.08333333333333  0.91666666666667  0.83333333333333  74.922 > 49   94 As  0.33333333333333  0.91666666666667  0.83333333333333  74.922 > 49   95 As  0.58333333333333  0.91666666666667  0.83333333333333  74.922 > 49   96 As  0.83333333333333  0.91666666666667  0.83333333333333  74.922 > 49   97 Ba  0.16666666666667  0.08333333333333  0.00000000000000 137.327 > 97   98 Ba  0.41666666666667  0.08333333333333  0.00000000000000 137.327 > 97   99 Ba  0.66666666666667  0.08333333333333  0.00000000000000 137.327 > 97  100 Ba  0.91666666666667  0.08333333333333  0.00000000000000 137.327 > 97  101 Ba  0.16666666666667  0.33333333333333  0.00000000000000 137.327 > 97  102 Ba  0.41666666666667  0.33333333333333  0.00000000000000 137.327 > 97  103 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327 > 97  104 Ba  0.91666666666667  0.33333333333333  0.00000000000000 137.327 > 97  105 Ba  0.16666666666667  0.58333333333333  0.00000000000000 137.327 > 97  106 Ba  0.41666666666667  0.58333333333333  0.00000000000000 137.327 > 97  107 Ba  0.66666666666667  0.58333333333333  0.00000000000000 137.327 > 97  108 Ba  0.91666666666667  0.58333333333333  0.00000000000000 137.327 > 97  109 Ba  0.16666666666667  0.83333333333333  0.00000000000000 137.327 > 97  110 Ba  0.41666666666667  0.83333333333333  0.00000000000000 137.327 > 97  111 Ba  0.66666666666667  0.83333333333333  0.00000000000000 137.327 > 97  112 Ba  0.91666666666667  0.83333333333333  0.00000000000000 137.327 > 97  113 Ba  0.16666666666667  0.08333333333333  0.33333333333333 137.327 > 97  114 Ba  0.41666666666667  0.08333333333333  0.33333333333333 137.327 > 97  115 Ba  0.66666666666667  0.08333333333333  0.33333333333333 137.327 > 97  116 Ba  0.91666666666667  0.08333333333333  0.33333333333333 137.327 > 97  117 Ba  0.16666666666667  0.33333333333333  0.33333333333333 137.327 > 97  118 Ba  0.41666666666667  0.33333333333333  0.33333333333333 137.327 > 97  119 Ba  0.66666666666667  0.33333333333333  0.33333333333333 137.327 > 97  120 Ba  0.91666666666667  0.33333333333333  0.33333333333333 137.327 > 97  121 Ba  0.16666666666667  0.58333333333333  0.33333333333333 137.327 > 97  122 Ba  0.41666666666667  0.58333333333333  0.33333333333333 137.327 > 97  123 Ba  0.66666666666667  0.58333333333333  0.33333333333333 137.327 > 97  124 Ba  0.91666666666667  0.58333333333333  0.33333333333333 137.327 > 97  125 Ba  0.16666666666667  0.83333333333333  0.33333333333333 137.327 > 97  126 Ba  0.41666666666667  0.83333333333333  0.33333333333333 137.327 > 97  127 Ba  0.66666666666667  0.83333333333333  0.33333333333333 137.327 > 97  128 Ba  0.91666666666667  0.83333333333333  0.33333333333333 137.327 > 97  129 Ba  0.16666666666667  0.08333333333333  0.66666666666667 137.327 > 97  130 Ba  0.41666666666667  0.08333333333333  0.66666666666667 137.327 > 97  131 Ba  0.66666666666667  0.08333333333333  0.66666666666667 137.327 > 97  132 Ba  0.91666666666667  0.08333333333333  0.66666666666667 137.327 > 97  133 Ba  0.16666666666667  0.33333333333333  0.66666666666667 137.327 > 97  134 Ba  0.41666666666667  0.33333333333333  0.66666666666667 137.327 > 97  135 Ba  0.66666666666667  0.33333333333333  0.66666666666667 137.327 > 97  136 Ba  0.91666666666667  0.33333333333333  0.66666666666667 137.327 > 97  137 Ba  0.16666666666667  0.58333333333333  0.66666666666667 137.327 > 97  138 Ba  0.41666666666667  0.58333333333333  0.66666666666667 137.327 > 97  139 Ba  0.66666666666667  0.58333333333333  0.66666666666667 137.327 > 97  140 Ba  0.91666666666667  0.58333333333333  0.66666666666667 137.327 > 97  141 Ba  0.16666666666667  0.83333333333333  0.66666666666667 137.327 > 97  142 Ba  0.41666666666667  0.83333333333333  0.66666666666667 137.327 > 97  143 Ba  0.66666666666667  0.83333333333333  0.66666666666667 137.327 > 97  144 Ba  0.91666666666667  0.83333333333333  0.66666666666667 137.327 > 97----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            9.6531051   -0.0000000    0.0000000            0.0000000    9.6531051    0.0000000            0.0000000    0.0000000   45.3917432-------------------------- Born effective charges --------------------------    1 Li    0.9957698    0.0000000    0.0000000            0.0000000    0.9957698    0.0000000            0.0000000    0.0000000   -0.0040015    2 As   -3.1893956   -0.0000000    0.0000000           -0.0000000   -3.1893956    0.0000000            0.0000000    0.0000000   -5.8823552    3 Ba    2.1936258    0.0000000    0.0000000            0.0000000    2.1936258    0.0000000            0.0000000    0.0000000    5.8863567----------------------------------------------------------------------------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[ 49, 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[ 97, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Expanding fc3.Symmetrizing fc3 by traditional approach (N=3).Symmetrizing fc2 by traditional approach (N=3).Max drift of fc3: -0.00000012 (xxy) -0.00000012 (xxy) -0.00000012 (xyx)fc3 was written into "fc3.hdf5".Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330----------- None of ph-ph interaction calculation was performed. -----------Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-07 23:44:45]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-07 23:44:46]-------------------------Compiled with OpenMP support (max 128 threads).Running in phono3py.load mode.Python version 3.14.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: conductivity-RTAHDF5 data compression filter: gzipCrystal structure was read from "phono3py.yaml".Supercell (dim): [4 4 3]Primitive matrix:  [1. 0. 0.]  [0. 1. 0.]  [0. 0. 1.]Spacegroup: P-6m2 (187)------------------------------ primitive cell ------------------------------Lattice vectors:  a    4.488603559793336    0.000000000000000    0.000000000000000  b   -2.244301779896668    3.887244710298293    0.000000000000000  c    0.000000000000000    0.000000000000000    4.539274940000000Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941    2 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922    3 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327-------------------------------- supercell ---------------------------------Lattice vectors:  a   17.954414239173342    0.000000000000000    0.000000000000000  b   -8.977207119586671   15.548978841193172    0.000000000000000  c    0.000000000000000    0.000000000000000   13.617824820000001Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.16666666666667   6.941 > 1    2 Li  0.25000000000000  0.00000000000000  0.16666666666667   6.941 > 1    3 Li  0.50000000000000  0.00000000000000  0.16666666666667   6.941 > 1    4 Li  0.75000000000000  0.00000000000000  0.16666666666667   6.941 > 1    5 Li  0.00000000000000  0.25000000000000  0.16666666666667   6.941 > 1    6 Li  0.25000000000000  0.25000000000000  0.16666666666667   6.941 > 1    7 Li  0.50000000000000  0.25000000000000  0.16666666666667   6.941 > 1    8 Li  0.75000000000000  0.25000000000000  0.16666666666667   6.941 > 1    9 Li  0.00000000000000  0.50000000000000  0.16666666666667   6.941 > 1   10 Li  0.25000000000000  0.50000000000000  0.16666666666667   6.941 > 1   11 Li  0.50000000000000  0.50000000000000  0.16666666666667   6.941 > 1   12 Li  0.75000000000000  0.50000000000000  0.16666666666667   6.941 > 1   13 Li  0.00000000000000  0.75000000000000  0.16666666666667   6.941 > 1   14 Li  0.25000000000000  0.75000000000000  0.16666666666667   6.941 > 1   15 Li  0.50000000000000  0.75000000000000  0.16666666666667   6.941 > 1   16 Li  0.75000000000000  0.75000000000000  0.16666666666667   6.941 > 1   17 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1   18 Li  0.25000000000000  0.00000000000000  0.50000000000000   6.941 > 1   19 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1   20 Li  0.75000000000000  0.00000000000000  0.50000000000000   6.941 > 1   21 Li  0.00000000000000  0.25000000000000  0.50000000000000   6.941 > 1   22 Li  0.25000000000000  0.25000000000000  0.50000000000000   6.941 > 1   23 Li  0.50000000000000  0.25000000000000  0.50000000000000   6.941 > 1   24 Li  0.75000000000000  0.25000000000000  0.50000000000000   6.941 > 1   25 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1   26 Li  0.25000000000000  0.50000000000000  0.50000000000000   6.941 > 1   27 Li  0.50000000000000  0.50000000000000  0.50000000000000   6.941 > 1   28 Li  0.75000000000000  0.50000000000000  0.50000000000000   6.941 > 1   29 Li  0.00000000000000  0.75000000000000  0.50000000000000   6.941 > 1   30 Li  0.25000000000000  0.75000000000000  0.50000000000000   6.941 > 1   31 Li  0.50000000000000  0.75000000000000  0.50000000000000   6.941 > 1   32 Li  0.75000000000000  0.75000000000000  0.50000000000000   6.941 > 1   33 Li  0.00000000000000  0.00000000000000  0.83333333333333   6.941 > 1   34 Li  0.25000000000000  0.00000000000000  0.83333333333333   6.941 > 1   35 Li  0.50000000000000  0.00000000000000  0.83333333333333   6.941 > 1   36 Li  0.75000000000000  0.00000000000000  0.83333333333333   6.941 > 1   37 Li  0.00000000000000  0.25000000000000  0.83333333333333   6.941 > 1   38 Li  0.25000000000000  0.25000000000000  0.83333333333333   6.941 > 1   39 Li  0.50000000000000  0.25000000000000  0.83333333333333   6.941 > 1   40 Li  0.75000000000000  0.25000000000000  0.83333333333333   6.941 > 1   41 Li  0.00000000000000  0.50000000000000  0.83333333333333   6.941 > 1   42 Li  0.25000000000000  0.50000000000000  0.83333333333333   6.941 > 1   43 Li  0.50000000000000  0.50000000000000  0.83333333333333   6.941 > 1   44 Li  0.75000000000000  0.50000000000000  0.83333333333333   6.941 > 1   45 Li  0.00000000000000  0.75000000000000  0.83333333333333   6.941 > 1   46 Li  0.25000000000000  0.75000000000000  0.83333333333333   6.941 > 1   47 Li  0.50000000000000  0.75000000000000  0.83333333333333   6.941 > 1   48 Li  0.75000000000000  0.75000000000000  0.83333333333333   6.941 > 1   49 As  0.08333333333333  0.16666666666667  0.16666666666667  74.922 > 49   50 As  0.33333333333333  0.16666666666667  0.16666666666667  74.922 > 49   51 As  0.58333333333333  0.16666666666667  0.16666666666667  74.922 > 49   52 As  0.83333333333333  0.16666666666667  0.16666666666667  74.922 > 49   53 As  0.08333333333333  0.41666666666667  0.16666666666667  74.922 > 49   54 As  0.33333333333333  0.41666666666667  0.16666666666667  74.922 > 49   55 As  0.58333333333333  0.41666666666667  0.16666666666667  74.922 > 49   56 As  0.83333333333333  0.41666666666667  0.16666666666667  74.922 > 49   57 As  0.08333333333333  0.66666666666667  0.16666666666667  74.922 > 49   58 As  0.33333333333333  0.66666666666667  0.16666666666667  74.922 > 49   59 As  0.58333333333333  0.66666666666667  0.16666666666667  74.922 > 49   60 As  0.83333333333333  0.66666666666667  0.16666666666667  74.922 > 49   61 As  0.08333333333333  0.91666666666667  0.16666666666667  74.922 > 49   62 As  0.33333333333333  0.91666666666667  0.16666666666667  74.922 > 49   63 As  0.58333333333333  0.91666666666667  0.16666666666667  74.922 > 49   64 As  0.83333333333333  0.91666666666667  0.16666666666667  74.922 > 49   65 As  0.08333333333333  0.16666666666667  0.50000000000000  74.922 > 49   66 As  0.33333333333333  0.16666666666667  0.50000000000000  74.922 > 49   67 As  0.58333333333333  0.16666666666667  0.50000000000000  74.922 > 49   68 As  0.83333333333333  0.16666666666667  0.50000000000000  74.922 > 49   69 As  0.08333333333333  0.41666666666667  0.50000000000000  74.922 > 49   70 As  0.33333333333333  0.41666666666667  0.50000000000000  74.922 > 49   71 As  0.58333333333333  0.41666666666667  0.50000000000000  74.922 > 49   72 As  0.83333333333333  0.41666666666667  0.50000000000000  74.922 > 49   73 As  0.08333333333333  0.66666666666667  0.50000000000000  74.922 > 49   74 As  0.33333333333333  0.66666666666667  0.50000000000000  74.922 > 49   75 As  0.58333333333333  0.66666666666667  0.50000000000000  74.922 > 49   76 As  0.83333333333333  0.66666666666667  0.50000000000000  74.922 > 49   77 As  0.08333333333333  0.91666666666667  0.50000000000000  74.922 > 49   78 As  0.33333333333333  0.91666666666667  0.50000000000000  74.922 > 49   79 As  0.58333333333333  0.91666666666667  0.50000000000000  74.922 > 49   80 As  0.83333333333333  0.91666666666667  0.50000000000000  74.922 > 49   81 As  0.08333333333333  0.16666666666667  0.83333333333333  74.922 > 49   82 As  0.33333333333333  0.16666666666667  0.83333333333333  74.922 > 49   83 As  0.58333333333333  0.16666666666667  0.83333333333333  74.922 > 49   84 As  0.83333333333333  0.16666666666667  0.83333333333333  74.922 > 49   85 As  0.08333333333333  0.41666666666667  0.83333333333333  74.922 > 49   86 As  0.33333333333333  0.41666666666667  0.83333333333333  74.922 > 49   87 As  0.58333333333333  0.41666666666667  0.83333333333333  74.922 > 49   88 As  0.83333333333333  0.41666666666667  0.83333333333333  74.922 > 49   89 As  0.08333333333333  0.66666666666667  0.83333333333333  74.922 > 49   90 As  0.33333333333333  0.66666666666667  0.83333333333333  74.922 > 49   91 As  0.58333333333333  0.66666666666667  0.83333333333333  74.922 > 49   92 As  0.83333333333333  0.66666666666667  0.83333333333333  74.922 > 49   93 As  0.08333333333333  0.91666666666667  0.83333333333333  74.922 > 49   94 As  0.33333333333333  0.91666666666667  0.83333333333333  74.922 > 49   95 As  0.58333333333333  0.91666666666667  0.83333333333333  74.922 > 49   96 As  0.83333333333333  0.91666666666667  0.83333333333333  74.922 > 49   97 Ba  0.16666666666667  0.08333333333333  0.00000000000000 137.327 > 97   98 Ba  0.41666666666667  0.08333333333333  0.00000000000000 137.327 > 97   99 Ba  0.66666666666667  0.08333333333333  0.00000000000000 137.327 > 97  100 Ba  0.91666666666667  0.08333333333333  0.00000000000000 137.327 > 97  101 Ba  0.16666666666667  0.33333333333333  0.00000000000000 137.327 > 97  102 Ba  0.41666666666667  0.33333333333333  0.00000000000000 137.327 > 97  103 Ba  0.66666666666667  0.33333333333333  0.00000000000000 137.327 > 97  104 Ba  0.91666666666667  0.33333333333333  0.00000000000000 137.327 > 97  105 Ba  0.16666666666667  0.58333333333333  0.00000000000000 137.327 > 97  106 Ba  0.41666666666667  0.58333333333333  0.00000000000000 137.327 > 97  107 Ba  0.66666666666667  0.58333333333333  0.00000000000000 137.327 > 97  108 Ba  0.91666666666667  0.58333333333333  0.00000000000000 137.327 > 97  109 Ba  0.16666666666667  0.83333333333333  0.00000000000000 137.327 > 97  110 Ba  0.41666666666667  0.83333333333333  0.00000000000000 137.327 > 97  111 Ba  0.66666666666667  0.83333333333333  0.00000000000000 137.327 > 97  112 Ba  0.91666666666667  0.83333333333333  0.00000000000000 137.327 > 97  113 Ba  0.16666666666667  0.08333333333333  0.33333333333333 137.327 > 97  114 Ba  0.41666666666667  0.08333333333333  0.33333333333333 137.327 > 97  115 Ba  0.66666666666667  0.08333333333333  0.33333333333333 137.327 > 97  116 Ba  0.91666666666667  0.08333333333333  0.33333333333333 137.327 > 97  117 Ba  0.16666666666667  0.33333333333333  0.33333333333333 137.327 > 97  118 Ba  0.41666666666667  0.33333333333333  0.33333333333333 137.327 > 97  119 Ba  0.66666666666667  0.33333333333333  0.33333333333333 137.327 > 97  120 Ba  0.91666666666667  0.33333333333333  0.33333333333333 137.327 > 97  121 Ba  0.16666666666667  0.58333333333333  0.33333333333333 137.327 > 97  122 Ba  0.41666666666667  0.58333333333333  0.33333333333333 137.327 > 97  123 Ba  0.66666666666667  0.58333333333333  0.33333333333333 137.327 > 97  124 Ba  0.91666666666667  0.58333333333333  0.33333333333333 137.327 > 97  125 Ba  0.16666666666667  0.83333333333333  0.33333333333333 137.327 > 97  126 Ba  0.41666666666667  0.83333333333333  0.33333333333333 137.327 > 97  127 Ba  0.66666666666667  0.83333333333333  0.33333333333333 137.327 > 97  128 Ba  0.91666666666667  0.83333333333333  0.33333333333333 137.327 > 97  129 Ba  0.16666666666667  0.08333333333333  0.66666666666667 137.327 > 97  130 Ba  0.41666666666667  0.08333333333333  0.66666666666667 137.327 > 97  131 Ba  0.66666666666667  0.08333333333333  0.66666666666667 137.327 > 97  132 Ba  0.91666666666667  0.08333333333333  0.66666666666667 137.327 > 97  133 Ba  0.16666666666667  0.33333333333333  0.66666666666667 137.327 > 97  134 Ba  0.41666666666667  0.33333333333333  0.66666666666667 137.327 > 97  135 Ba  0.66666666666667  0.33333333333333  0.66666666666667 137.327 > 97  136 Ba  0.91666666666667  0.33333333333333  0.66666666666667 137.327 > 97  137 Ba  0.16666666666667  0.58333333333333  0.66666666666667 137.327 > 97  138 Ba  0.41666666666667  0.58333333333333  0.66666666666667 137.327 > 97  139 Ba  0.66666666666667  0.58333333333333  0.66666666666667 137.327 > 97  140 Ba  0.91666666666667  0.58333333333333  0.66666666666667 137.327 > 97  141 Ba  0.16666666666667  0.83333333333333  0.66666666666667 137.327 > 97  142 Ba  0.41666666666667  0.83333333333333  0.66666666666667 137.327 > 97  143 Ba  0.66666666666667  0.83333333333333  0.66666666666667 137.327 > 97  144 Ba  0.91666666666667  0.83333333333333  0.66666666666667 137.327 > 97----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            9.6531051   -0.0000000    0.0000000            0.0000000    9.6531051    0.0000000            0.0000000    0.0000000   45.3917432-------------------------- Born effective charges --------------------------    1 Li    0.9957698    0.0000000    0.0000000            0.0000000    0.9957698    0.0000000            0.0000000    0.0000000   -0.0040015    2 As   -3.1893956   -0.0000000    0.0000000           -0.0000000   -3.1893956    0.0000000            0.0000000    0.0000000   -5.8823552    3 Ba    2.1936258    0.0000000    0.0000000            0.0000000    2.1936258    0.0000000            0.0000000    0.0000000    5.8863567----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: -0.00000012 (xxy) -0.00000012 (xxy) -0.00000012 (xyx)Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330Length for sampling mesh generation: 50.00Generating grid system ... [ 13 13 11 ]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.97, Number of G-points: 315, Lambda: 0.47Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/126) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126Calculating 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.927   (   0.000    0.000    0.000)    0.000   3.032   (   0.000    0.000    0.000)    0.000   3.032   (   0.000    0.000    0.000)    0.000   9.134   (   0.000    0.000    0.000)    0.000  10.763   (   0.000   -0.000    0.000)    0.000  10.763   (   0.000    0.000    0.000)    0.000======================= Grid point 1 (2/126) =======================q-point: ( 0.08  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.391   (  16.779    9.688    0.000)   19.375   0.485   (  20.587   11.886    0.000)   23.771   0.812   (  34.575   19.962    0.000)   39.924   1.951   (   1.974    1.140    0.000)    2.280   3.035   (   0.289    0.167    0.000)    0.334   3.947   (   2.713    1.566    0.000)    3.133   9.134   (   0.050    0.029    0.000)    0.057  10.759   (  -0.321   -0.185    0.000)    0.371  11.306   (  -0.151   -0.087    0.000)    0.175======================= Grid point 2 (3/126) =======================q-point: ( 0.15  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.760   (  15.403    8.893    0.000)   17.786   0.928   (  17.984   10.383    0.000)   20.766   1.558   (  30.283   17.484    0.000)   34.968   2.007   (   2.661    1.536    0.000)    3.073   3.047   (   0.735    0.425    0.000)    0.849   4.032   (   4.661    2.691    0.000)    5.382   9.137   (   0.230    0.133    0.000)    0.266  10.748   (  -0.574   -0.332    0.000)    0.663  11.296   (  -0.863   -0.498    0.000)    0.996======================= Grid point 3 (4/126) =======================q-point: ( 0.23  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.089   (  13.243    7.646    0.000)   15.292   1.303   (  14.819    8.556    0.000)   17.111   2.061   (   1.833    1.058    0.000)    2.117   2.177   (  23.536   13.589    0.000)   27.178   3.063   (   0.459    0.265    0.000)    0.530   4.155   (   5.960    3.441    0.000)    6.882   9.146   (   0.521    0.301    0.000)    0.602  10.735   (  -0.554   -0.320    0.000)    0.640  11.262   (  -2.147   -1.239    0.000)    2.479======================= Grid point 4 (5/126) =======================q-point: ( 0.31  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.360   (  10.396    6.002    0.000)   12.004   1.609   (  12.102    6.987    0.000)   13.974   2.085   (   0.291    0.168    0.000)    0.336   2.625   (  15.647    9.034    0.000)   18.068   3.057   (  -1.256   -0.725    0.000)    1.450   4.300   (   6.583    3.801    0.000)    7.602   9.160   (   0.736    0.425    0.000)    0.850  10.726   (  -0.220   -0.127    0.000)    0.254  11.199   (  -3.267   -1.886    0.000)    3.772======================= Grid point 5 (6/126) =======================q-point: ( 0.38  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.559   (   6.897    3.982    0.000)    7.964   1.857   (   9.401    5.428    0.000)   10.855   2.077   (  -0.827   -0.477    0.000)    0.954   2.898   (   8.505    4.910    0.000)    9.820   3.004   (  -3.186   -1.839    0.000)    3.679   4.442   (   5.432    3.136    0.000)    6.272   9.177   (   0.667    0.385    0.000)    0.770  10.725   (   0.108    0.063    0.000)    0.125  11.122   (  -3.141   -1.814    0.000)    3.627======================= Grid point 6 (7/126) =======================q-point: ( 0.46  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.669   (   2.535    1.464    0.000)    2.928   2.019   (   4.112    2.374    0.000)    4.748   2.058   (  -0.614   -0.355    0.000)    0.709   2.934   (  -2.129   -1.229    0.000)    2.458   3.024   (   2.668    1.540    0.000)    3.081   4.531   (   2.114    1.220    0.000)    2.440   9.188   (   0.271    0.156    0.000)    0.312  10.728   (   0.105    0.061    0.000)    0.122  11.068   (  -1.318   -0.761    0.000)    1.522======================= Grid point 14 (8/126) =======================q-point: ( 0.08  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.665   (   9.174   15.889    0.000)   18.347   0.841   (  11.841   20.509    0.000)   23.682   1.355   (  17.819   30.864    0.000)   35.638   1.991   (   1.536    2.660    0.000)    3.072   3.051   (   0.824    1.427    0.000)    1.648   3.999   (   2.066    3.578    0.000)    4.132   9.136   (   0.104    0.180    0.000)    0.208  10.757   (  -0.046   -0.079    0.000)    0.092  11.295   (  -0.619   -1.072    0.000)    1.237======================= Grid point 15 (9/126) =======================q-point: ( 0.15  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.979   (  10.653   12.469    0.000)   16.400   1.232   (  10.632   18.354    0.000)   21.211   1.951   (  20.566   20.806    0.000)   29.255   2.045   (   1.725    2.117    0.000)    2.731   3.088   (   0.293    3.719    0.000)    3.730   4.085   (   4.481    2.558    0.000)    5.160   9.142   (   0.323    0.451    0.000)    0.555  10.753   (  -0.770    0.689    0.000)    1.033  11.262   (  -1.000   -2.959    0.000)    3.124======================= Grid point 16 (10/126) =======================q-point: ( 0.23  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.264   (   9.560    9.654    0.000)   13.587   1.561   (   7.862   16.126    0.000)   17.941   2.083   (   0.708    1.146    0.000)    1.347   2.442   (  16.983   12.539    0.000)   21.111   3.122   (  -2.171    5.772    0.000)    6.167   4.203   (   6.911    1.366    0.000)    7.045   9.156   (   0.529    0.803    0.000)    0.962  10.744   (  -1.217    1.026    0.000)    1.592  11.202   (  -1.438   -4.679    0.000)    4.895======================= Grid point 17 (11/126) =======================q-point: ( 0.31  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.491   (   7.149    6.928    0.000)    9.956   1.823   (   5.704   13.421    0.000)   14.583   2.091   (  -0.635    0.479    0.000)    0.796   2.775   (  11.383    5.777    0.000)   12.765   3.117   (  -5.854    6.975    0.000)    9.106   4.342   (   7.889    0.518    0.000)    7.906   9.175   (   0.502    1.065    0.000)    1.177  10.734   (  -0.888    0.671    0.000)    1.113  11.125   (  -1.537   -5.306    0.000)    5.524======================= Grid point 18 (12/126) =======================q-point: ( 0.38  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.640   (   3.788    4.014    0.000)    5.519   2.016   (   2.999    9.633    0.000)   10.089   2.078   (  -1.317    0.608    0.000)    1.451   2.955   (   6.228    0.736    0.000)    6.272   3.069   (  -8.082    8.538    0.000)   11.756   4.464   (   5.887   -0.827    0.000)    5.945   9.192   (   0.127    1.077    0.000)    1.085  10.728   (  -0.223    0.068    0.000)    0.233  11.056   (  -0.474   -4.658    0.000)    4.682======================= Grid point 19 (13/126) =======================q-point: ( 0.46  0.08  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.694   (  -0.592    1.025    0.000)    1.184   2.068   (  -0.751    1.300    0.000)    1.502   2.094   (  -2.838    4.914    0.000)    5.674   3.010   (   1.135   -1.967    0.000)    2.271   3.038   (  -5.925   10.260    0.000)   11.848   4.512   (   1.787   -3.101    0.000)    3.579   9.198   (  -0.473    0.820    0.000)    0.947  10.727   (   0.109   -0.189    0.000)    0.218  11.027   (   1.865   -3.233    0.000)    3.732======================= Grid point 29 (14/126) =======================q-point: ( 0.15  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.229   (   7.126   12.342    0.000)   14.252   1.590   (   9.758   16.902    0.000)   19.517   2.083   (   0.901    1.560    0.000)    1.802   2.346   (  10.611   18.379    0.000)   21.222   3.183   (   3.149    5.455    0.000)    6.299   4.138   (   1.757    3.044    0.000)    3.515   9.156   (   0.527    0.912    0.000)    1.053  10.766   (   0.206    0.356    0.000)    0.411  11.184   (  -2.666   -4.617    0.000)    5.332======================= Grid point 30 (15/126) =======================q-point: ( 0.23  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.454   (   6.389    9.179    0.000)   11.184   1.884   (   5.733   15.448    0.000)   16.478   2.104   (   0.086    0.874    0.000)    0.878   2.656   (   9.563    8.768    0.000)   12.974   3.281   (  -1.185    9.667    0.000)    9.739   4.216   (   5.536    0.339    0.000)    5.547   9.178   (   0.645    1.338    0.000)    1.486  10.762   (  -1.142    0.412    0.000)    1.214  11.089   (  -1.845   -6.312    0.000)    6.576======================= Grid point 31 (16/126) =======================q-point: ( 0.31  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.621   (   4.162    5.892    0.000)    7.214   2.093   (   2.150   13.000    0.000)   13.177   2.106   (  -0.785    1.031    0.000)    1.296   2.844   (   6.903    1.106    0.000)    6.991   3.323   (  -6.907   13.280    0.000)   14.969   4.326   (   7.612   -1.979    0.000)    7.865   9.201   (   0.359    1.553    0.000)    1.594  10.740   (  -1.214   -0.277    0.000)    1.246  11.004   (  -0.492   -6.403    0.000)    6.422======================= Grid point 32 (17/126) =======================q-point: ( 0.38  0.15  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.709   (   1.010    2.735    0.000)    2.915   2.102   (  -1.172    1.729    0.000)    2.088   2.209   (  -1.928    9.554    0.000)    9.747   2.925   (   4.286   -3.687    0.000)    5.654   3.322   (  -9.564   15.928    0.000)   18.578   4.412   (   5.382   -4.400    0.000)    6.951   9.217   (  -0.311    1.412    0.000)    1.446  10.722   (  -0.174   -0.696    0.000)    0.717  10.953   (   1.418   -5.300    0.000)    5.487======================= Grid point 43 (18/126) =======================q-point: ( 0.23  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.615   (   3.897    6.751    0.000)    7.795   2.119   (   0.392    0.678    0.000)    0.783   2.150   (   6.477   11.218    0.000)   12.953   2.781   (   2.249    3.896    0.000)    4.498   3.473   (   5.134    8.892    0.000)   10.268   4.230   (   0.777    1.345    0.000)    1.554   9.207   (   0.859    1.488    0.000)    1.718  10.753   (  -0.738   -1.278    0.000)    1.476  10.971   (  -3.014   -5.220    0.000)    6.028======================= Grid point 44 (19/126) =======================q-point: ( 0.31  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.715   (   1.835    3.467    0.000)    3.923   2.131   (  -0.178    1.398    0.000)    1.409   2.324   (   1.268   10.346    0.000)   10.424   2.821   (   2.503   -3.220    0.000)    4.078   3.605   (  -3.257   14.103    0.000)   14.475   4.269   (   4.656   -3.771    0.000)    5.992   9.232   (   0.323    1.422    0.000)    1.458  10.721   (  -0.923   -1.498    0.000)    1.760  10.891   (  -0.521   -4.629    0.000)    4.658======================= Grid point 45 (20/126) =======================q-point: ( 0.38  0.23  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.747   (  -0.692    1.199    0.000)    1.384   2.137   (  -0.896    1.551    0.000)    1.791   2.388   (  -4.828    8.359    0.000)    9.653   2.823   (   3.505   -6.071    0.000)    7.010   3.645   (  -8.819   15.271    0.000)   17.635   4.294   (   4.146   -7.186    0.000)    8.296   9.242   (  -0.579    1.004    0.000)    1.159  10.704   (   0.542   -0.941    0.000)    1.086  10.863   (   1.928   -3.343    0.000)    3.859======================= Grid point 58 (21/126) =======================q-point: ( 0.31  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.759   (   0.603    1.044    0.000)    1.206   2.154   (   0.404    0.700    0.000)    0.808   2.495   (   3.589    6.216    0.000)    7.178   2.737   (  -2.405   -4.165    0.000)    4.809   3.836   (   4.447    7.703    0.000)    8.895   4.187   (  -2.087   -3.615    0.000)    4.174   9.253   (   0.334    0.578    0.000)    0.668  10.693   (  -0.537   -0.931    0.000)    1.075  10.830   (  -0.905   -1.568    0.000)    1.810======================= Grid point 181 (22/126) =======================q-point: ( 0.00  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.360   (   0.000    0.000   18.732)   18.732   0.360   (   0.000    0.000   18.732)   18.732   0.792   (   0.000    0.000   39.166)   39.166   3.002   (   0.000    0.000   -2.938)    2.938   3.002   (   0.000    0.000   -2.938)    2.938   3.650   (   0.000    0.000   10.673)   10.673   9.133   (   0.000    0.000   -0.084)    0.084  10.760   (   0.000    0.000   -0.201)    0.201  10.760   (   0.000    0.000   -0.201)    0.201======================= Grid point 182 (23/126) =======================q-point: ( 0.08  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.605   (  16.544    9.552   11.174)   22.131   0.680   (  15.390    8.885   15.527)   23.598   1.039   (  21.967   12.683   24.882)   35.532   2.844   ( -10.859   -6.269   18.103)   22.021   3.006   (   0.387    0.224   -2.810)    2.845   3.724   (   5.171    2.985   -1.758)    6.224   9.135   (   0.200    0.115    0.089)    0.247  10.756   (  -0.357   -0.206   -0.243)    0.479  10.864   (   7.357    4.247   -8.010)   11.676======================= Grid point 183 (24/126) =======================q-point: ( 0.15  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.983   (  12.335    7.121   17.113)   22.265   0.996   (  16.777    9.686    6.800)   20.531   1.631   (  26.440   15.265    9.121)   31.863   2.613   (  -8.343   -4.817   30.579)   32.061   3.021   (   0.845    0.488   -2.541)    2.721   3.856   (   6.270    3.620   -9.921)   12.282   9.143   (   0.497    0.287    0.560)    0.802  10.745   (  -0.621   -0.358   -0.338)    0.793  11.021   (   5.454    3.149  -12.743)   14.215======================= Grid point 184 (25/126) =======================q-point: ( 0.23  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.252   (  11.099    6.408   13.578)   18.671   1.352   (  14.233    8.217    4.961)   17.167   2.169   (  19.211   11.092    4.141)   22.567   2.477   (  -2.996   -1.730   25.899)   26.129   3.039   (   0.563    0.325   -2.270)    2.362   4.009   (   7.138    4.121  -11.148)   13.864   9.158   (   0.853    0.492    1.230)    1.576  10.730   (  -0.597   -0.344   -0.435)    0.815  11.101   (   1.652    0.954  -10.504)   10.676======================= Grid point 185 (26/126) =======================q-point: ( 0.31  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.482   (   8.893    5.134   10.336)   14.570   1.647   (  11.589    6.691    3.750)   13.897   2.335   (  -1.591   -0.919   19.529)   19.615   2.623   (  12.400    7.159    0.438)   14.325   3.036   (  -1.056   -0.610   -1.901)    2.258   4.178   (   7.406    4.276  -10.197)   13.308   9.181   (   1.060    0.612    1.966)    2.316  10.720   (  -0.262   -0.151   -0.526)    0.607  11.105   (  -1.035   -0.597   -7.084)    7.184======================= Grid point 186 (27/126) =======================q-point: ( 0.38  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.653   (   5.902    3.407    8.058)   10.553   1.882   (   8.816    5.090    2.488)   10.479   2.290   (  -1.883   -1.087   17.765)   17.897   2.851   (   7.270    4.197   -4.293)    9.428   2.989   (  -2.826   -1.631   -1.220)    3.484   4.333   (   5.814    3.357   -9.296)   11.467   9.204   (   0.895    0.517    2.584)    2.783  10.718   (   0.070    0.041   -0.616)    0.621  11.068   (  -1.881   -1.086   -4.381)    4.890======================= Grid point 187 (28/126) =======================q-point: ( 0.46  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.746   (   2.127    1.228    6.852)    7.279   2.032   (   3.765    2.174    1.328)    4.545   2.258   (  -0.782   -0.452   16.642)   16.667   2.928   (  -1.870   -1.080   -0.404)    2.197   2.959   (   2.288    1.321   -5.923)    6.485   4.428   (   2.202    1.271   -8.971)    9.324   9.218   (   0.349    0.201    2.922)    2.949  10.720   (   0.089    0.051   -0.679)    0.687  11.033   (  -0.892   -0.515   -2.885)    3.064======================= Grid point 195 (29/126) =======================q-point: ( 0.08  0.08  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.904   (   7.263   12.579   17.421)   22.682   0.914   (  10.831   18.760    7.271)   22.850   1.457   (  15.149   26.238   12.123)   32.633   2.663   (  -5.675   -9.829   29.016)   31.156   3.026   (   1.005    1.741   -2.297)    3.053   3.813   (   3.151    5.457   -8.581)   10.646   9.140   (   0.237    0.411    0.404)    0.623  10.753   (  -0.077   -0.134   -0.321)    0.357  10.979   (   3.431    5.943  -12.456)   14.222======================= Grid point 196 (30/126) =======================q-point: ( 0.15  0.08  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.159   (   8.892   10.194   14.563)   19.876   1.281   (  10.336   17.270    5.013)   20.742   1.978   (  17.940   18.366    5.151)   26.185   2.501   (  -3.821   -5.529   28.137)   28.929   3.069   (   0.232    4.266   -1.661)    4.584   3.932   (   5.377    3.779  -11.020)   12.831   9.152   (   0.560    0.681    0.941)    1.289  10.748   (  -0.808    0.636   -0.430)    1.114  11.067   (   2.491    1.066  -11.610)   11.922======================= Grid point 197 (31/126) =======================q-point: ( 0.23  0.08  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.398   (   8.225    7.967   11.173)   15.999   1.597   (   7.749   15.217    3.704)   17.474   2.329   (   3.393    2.279   14.638)   15.198   2.489   (   9.211    5.559    8.279)   13.575   3.106   (  -2.385    6.319   -1.369)    6.892   4.074   (   7.676    2.290  -10.539)   13.238   9.172   (   0.834    0.969    1.587)    2.038  10.738   (  -1.238    0.980   -0.523)    1.663  11.087   (   0.841   -2.426   -8.331)    8.717======================= Grid point 198 (32/126) =======================q-point: ( 0.31  0.08  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.590   (   6.277    5.529    8.429)   11.875   1.848   (   5.531   12.560    2.569)   13.962   2.309   (  -1.603   -1.349   18.022)   18.143   2.743   (   9.828    4.634   -2.871)   11.239   3.103   (  -5.878    7.377   -1.144)    9.501   4.232   (   8.289    1.201   -9.392)   12.584   9.198   (   0.807    1.114    2.202)    2.597  10.728   (  -0.904    0.662   -0.577)    1.260  11.058   (  -0.155   -4.006   -5.310)    6.653======================= Grid point 199 (33/126) =======================q-point: ( 0.38  0.08  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.717   (   3.454    2.882    6.671)    8.046   2.029   (   2.775    8.889    1.386)    9.414   2.275   (  -1.339   -0.184   16.487)   16.543   2.897   (   5.452    0.348   -5.325)    7.629   3.061   (  -7.830    8.777   -0.602)   11.777   4.362   (   5.893   -0.362   -8.774)   10.576   9.219   (   0.356    0.991    2.629)    2.832  10.722   (  -0.260    0.120   -0.599)    0.664  11.015   (   0.196   -4.003   -3.289)    5.185======================= Grid point 200 (34/126) =======================q-point: ( 0.46  0.08  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.761   (  -0.192    0.333    6.058)    6.070   2.101   (  -2.599    4.501    0.773)    5.255   2.262   (  -0.500    0.865   16.105)   16.136   2.944   (   1.413   -2.449   -6.018)    6.649   3.033   (  -6.230   10.788   -0.258)   12.460   4.413   (   1.588   -2.756   -8.641)    9.208   9.227   (  -0.364    0.631    2.775)    2.869  10.720   (   0.060   -0.105   -0.606)    0.618  10.996   (   1.786   -3.096   -2.569)    4.401======================= Grid point 210 (35/126) =======================q-point: ( 0.15  0.15  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.366   (   5.957   10.317   11.406)   16.493   1.621   (   9.320   16.144    3.168)   18.908   2.295   (   5.607    9.712    8.482)   14.061   2.441   (   1.643    2.846   15.318)   15.667   3.173   (   3.310    5.733   -0.725)    6.659   4.008   (   2.381    4.124  -10.640)   11.657   9.170   (   0.666    1.153    1.398)    1.930  10.760   (   0.194    0.336   -0.489)    0.624  11.060   (  -0.937   -1.623   -8.821)    9.018======================= Grid point 211 (36/126) =======================q-point: ( 0.23  0.15  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.554   (   5.412    7.478    8.414)   12.490   1.902   (   5.580   14.547    1.854)   15.691   2.318   (  -0.979   -1.223   17.655)   17.725   2.637   (   8.138    7.023   -1.736)   10.888   3.272   (  -1.368    9.821   -0.587)    9.933   4.106   (   6.101    1.210   -9.407)   11.278   9.196   (   0.833    1.422    1.827)    2.461  10.756   (  -1.178    0.481   -0.494)    1.365  11.013   (  -0.536   -4.589   -5.882)    7.480======================= Grid point 212 (37/126) =======================q-point: ( 0.31  0.15  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.692   (   3.614    4.519    6.166)    8.456   2.099   (   2.090   11.999    0.752)   12.203   2.295   (  -0.950    0.089   15.951)   15.980   2.796   (   6.091    0.545   -4.489)    7.586   3.313   (  -6.894   13.182   -0.762)   14.895   4.230   (   7.769   -1.268   -8.310)   11.447   9.224   (   0.565    1.452    2.144)    2.650  10.735   (  -1.243   -0.125   -0.434)    1.323  10.960   (   0.286   -5.470   -3.647)    6.581======================= Grid point 213 (38/126) =======================q-point: ( 0.38  0.15  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.764   (   1.047    1.853    5.043)    5.474   2.205   (  -1.742    8.359    0.297)    8.544   2.287   (  -0.937    1.554   14.945)   15.055   2.864   (   3.895   -3.585   -5.529)    7.655   3.312   (  -9.306   15.553   -0.735)   18.140   4.321   (   5.164   -3.791   -7.894)   10.166   9.241   (  -0.128    1.176    2.299)    2.585  10.718   (  -0.250   -0.497   -0.346)    0.655  10.923   (   1.635   -4.931   -2.468)    5.751======================= Grid point 224 (39/126) =======================q-point: ( 0.23  0.23  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.683   (   3.126    5.414    5.962)    8.639   2.150   (   5.962   10.326    0.292)   11.927   2.305   (  -0.010   -0.017   15.741)   15.741   2.737   (   1.776    3.075   -4.110)    5.431   3.463   (   4.977    8.621   -0.713)    9.980   4.135   (   1.191    2.063   -8.209)    8.548   9.226   (   0.863    1.495    1.906)    2.572  10.750   (  -0.655   -1.135   -0.239)    1.333  10.923   (  -2.368   -4.101   -3.901)    6.135======================= Grid point 225 (40/126) =======================q-point: ( 0.31  0.23  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.764   (   1.558    2.708    4.513)    5.488   2.285   (   0.365    3.614    7.548)    8.377   2.335   (   0.598    6.912    6.342)    9.400   2.767   (   2.234   -3.271   -4.936)    6.329   3.589   (  -3.344   13.663   -1.294)   14.126   4.185   (   4.804   -3.184   -7.280)    9.285   9.252   (   0.385    1.267    1.878)    2.298  10.720   (  -0.941   -1.241    0.028)    1.557  10.860   (  -0.155   -4.159   -2.617)    4.917======================= Grid point 226 (41/126) =======================q-point: ( 0.38  0.23  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.790   (  -0.504    0.872    4.110)    4.232   2.298   (  -0.829    1.436   11.890)   12.005   2.385   (  -4.340    7.516    1.192)    8.760   2.766   (   3.356   -5.814   -5.045)    8.397   3.626   (  -8.546   14.797   -1.579)   17.160   4.215   (   3.859   -6.689   -6.943)   10.385   9.262   (  -0.462    0.801    1.848)    2.067  10.705   (   0.408   -0.708    0.137)    0.829  10.838   (   1.878   -3.256   -2.167)    4.339======================= Grid point 239 (42/126) =======================q-point: ( 0.31  0.31  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.799   (   0.501    0.867    3.863)    3.991   2.312   (   0.309    0.536   12.558)   12.573   2.487   (   3.578    6.197   -0.203)    7.158   2.682   (  -2.545   -4.408   -4.906)    7.070   3.813   (   4.391    7.605   -1.769)    8.958   4.112   (  -2.024   -3.505   -6.770)    7.887   9.270   (   0.291    0.504    1.639)    1.739  10.696   (  -0.477   -0.827    0.429)    1.047  10.804   (  -0.816   -1.413   -2.169)    2.714======================= Grid point 362 (43/126) =======================q-point: ( 0.00  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.757   (   0.000    0.000   20.887)   20.887   0.757   (   0.000    0.000   20.887)   20.887   1.563   (   0.000   -0.000   37.709)   37.709   2.918   (   0.000    0.000   -5.435)    5.435   2.918   (   0.000   -0.000   -5.435)    5.435   3.908   (   0.000    0.000   13.412)   13.412   9.132   (   0.000    0.000   -0.010)    0.010  10.756   (   0.000    0.000   -0.158)    0.158  10.756   (   0.000    0.000   -0.158)    0.158======================= Grid point 363 (44/126) =======================q-point: ( 0.08  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.899   (  11.074    6.393   17.497)   21.672   0.972   (  14.867    8.584   15.824)   23.348   1.662   (   8.781    5.070   33.623)   35.119   2.926   (   0.723    0.417   -5.143)    5.210   2.952   (   2.534    1.463   -1.709)    3.388   3.850   (  -4.438   -2.562   10.312)   11.515   9.138   (   0.599    0.346    0.268)    0.742  10.751   (  -0.447   -0.258   -0.221)    0.561  10.793   (   2.958    1.708   -1.416)    3.697======================= Grid point 364 (45/126) =======================q-point: ( 0.15  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.196   (  13.842    7.992   12.881)   20.528   1.291   (  11.860    6.847   13.271)   19.069   1.954   (  15.748    9.092   22.032)   28.567   2.950   (   1.301    0.751   -4.436)    4.684   2.993   (  -0.053   -0.031    9.479)    9.479   3.759   (  -1.953   -1.128   -0.121)    2.258   9.159   (   1.196    0.691    1.032)    1.724  10.737   (  -0.737   -0.425   -0.363)    0.925  10.874   (   3.777    2.181   -3.633)    5.677======================= Grid point 365 (46/126) =======================q-point: ( 0.23  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.501   (  12.461    7.194    9.720)   17.364   1.519   (   8.529    4.924   11.515)   15.152   2.313   (  14.402    8.315   10.608)   19.726   2.921   (  -5.826   -3.363   17.466)   18.717   2.979   (   1.088    0.628   -3.530)    3.747   3.804   (   5.683    3.281   -8.626)   10.838   9.192   (   1.712    0.989    2.112)    2.893  10.720   (  -0.704   -0.406   -0.512)    0.960  10.949   (   2.559    1.477   -4.649)    5.508======================= Grid point 366 (47/126) =======================q-point: ( 0.31  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.691   (   6.651    3.840    9.062)   11.879   1.758   (  10.005    5.777    7.232)   13.630   2.579   (   8.193    4.730    2.710)    9.841   2.777   (  -5.321   -3.072   15.493)   16.667   2.991   (  -0.253   -0.146   -2.342)    2.360   3.975   (   8.243    4.759   -8.708)   12.900   9.235   (   1.886    1.089    3.273)    3.932  10.707   (  -0.367   -0.212   -0.653)    0.778  10.987   (   0.777    0.449   -4.059)    4.157======================= Grid point 367 (48/126) =======================q-point: ( 0.38  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.822   (   4.676    2.700    7.546)    9.279   1.955   (   7.137    4.121    4.777)    9.526   2.605   (  -3.527   -2.036    6.942)    8.049   2.786   (   4.088    2.360    2.831)    5.504   2.967   (  -1.659   -0.958   -0.694)    2.038   4.146   (   6.239    3.602   -7.804)   10.621   9.274   (   1.471    0.849    4.232)    4.560  10.703   (  -0.026   -0.015   -0.790)    0.791  10.991   (  -0.268   -0.155   -2.658)    2.676======================= Grid point 368 (49/126) =======================q-point: ( 0.46  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.898   (   1.783    1.029    7.141)    7.431   2.073   (   2.828    1.633    2.751)    4.270   2.542   (  -1.531   -0.884    6.262)    6.507   2.853   (   1.505    0.869   -0.850)    1.935   2.930   (  -1.125   -0.649    0.909)    1.585   4.245   (   2.247    1.297   -7.637)    8.066   9.298   (   0.548    0.317    4.760)    4.802  10.704   (   0.048    0.027   -0.887)    0.889  10.983   (  -0.246   -0.142   -1.551)    1.577======================= Grid point 376 (50/126) =======================q-point: ( 0.08  0.08  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.126   (   8.558   14.823   13.615)   21.871   1.208   (   7.434   12.876   13.544)   20.112   1.858   (   8.214   14.227   25.696)   30.499   2.930   (  -0.161   -0.280   -1.673)    1.704   3.004   (   1.889    3.271    2.908)    4.767   3.774   (  -2.265   -3.923    3.368)    5.645   9.152   (   0.591    1.023    0.781)    1.416  10.746   (  -0.156   -0.269   -0.338)    0.459  10.846   (   1.936    3.353   -3.077)    4.945======================= Grid point 377 (51/126) =======================q-point: ( 0.15  0.08  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.424   (   8.119   12.397   10.694)   18.275   1.446   (   8.532    9.912   11.396)   17.347   2.178   (  12.065   12.533   14.692)   22.771   2.895   (  -1.392   -5.063    8.270)    9.796   3.060   (  -0.683    5.406    3.888)    6.693   3.749   (   2.433    0.620   -6.752)    7.203   9.179   (   1.175    1.290    1.658)    2.407  10.738   (  -0.891    0.514   -0.494)    1.141  10.912   (   2.655    1.302   -4.439)    5.334======================= Grid point 378 (52/126) =======================q-point: ( 0.23  0.08  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.620   (   6.406    5.545    9.631)   12.827   1.710   (   7.120   12.559    7.435)   16.238   2.474   (  10.075    7.449    5.574)   13.714   2.789   (  -3.121   -7.011   13.843)   15.828   3.088   (  -3.740    7.905    0.562)    8.763   3.866   (   8.163    2.758   -9.170)   12.583   9.216   (   1.607    1.404    2.690)    3.434  10.726   (  -1.296    0.905   -0.611)    1.695  10.955   (   2.145   -1.031   -4.378)    4.983======================= Grid point 379 (53/126) =======================q-point: ( 0.31  0.08  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.764   (   5.174    3.512    7.691)    9.913   1.926   (   4.834   10.140    5.035)   12.310   2.619   (  -0.169   -1.574    3.733)    4.055   2.722   (   3.787   -0.319    8.167)    9.008   3.085   (  -6.081    8.615   -0.261)   10.548   4.044   (   8.812    2.157   -7.913)   12.038   9.258   (   1.574    1.237    3.645)    4.158  10.714   (  -0.980    0.697   -0.683)    1.383  10.966   (   1.347   -2.537   -3.211)    4.309======================= Grid point 380 (54/126) =======================q-point: ( 0.38  0.08  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.863   (   3.303    1.453    6.855)    7.747   2.072   (   2.143    6.907    2.841)    7.770   2.560   (  -1.971   -1.893    6.054)    6.642   2.806   (   3.858   -0.232    0.444)    3.891   3.057   (  -7.212    9.347    0.530)   11.818   4.185   (   5.736    0.341   -7.371)    9.347   9.290   (   0.936    0.766    4.298)    4.465  10.707   (  -0.392    0.289   -0.739)    0.885  10.958   (   1.166   -3.043   -1.863)    3.754======================= Grid point 381 (55/126) =======================q-point: ( 0.46  0.08  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.899   (   0.380   -0.658    6.733)    6.776   2.126   (  -2.049    3.548    1.803)    4.476   2.533   (   0.003   -0.005    5.383)    5.383   2.837   (   1.168   -2.024   -0.653)    2.427   3.039   (  -6.128   10.612    1.063)   12.300   4.237   (   1.235   -2.145   -7.343)    7.749   9.302   (  -0.083    0.143    4.524)    4.527  10.706   (  -0.076    0.132   -0.768)    0.783  10.952   (   1.655   -2.871   -1.273)    3.550======================= Grid point 391 (56/126) =======================q-point: ( 0.15  0.15  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.592   (   4.287    7.425    9.801)   13.022   1.720   (   7.988   13.836    6.703)   17.326   2.417   (   6.341   10.983    7.137)   14.553   2.778   (  -3.585   -6.209   13.102)   14.935   3.177   (   3.394    5.878    1.981)    7.071   3.796   (   2.458    4.258   -9.521)   10.715   9.209   (   1.026    1.777    2.406)    3.162  10.749   (   0.195    0.337   -0.550)    0.674  10.923   (  -0.064   -0.111   -4.469)    4.471======================= Grid point 392 (57/126) =======================q-point: ( 0.23  0.15  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.726   (   4.059    4.991    7.568)    9.933   1.961   (   4.989   11.945    4.053)   13.564   2.590   (   3.817    3.472   -0.814)    5.224   2.676   (  -0.909   -3.523   14.304)   14.760   3.271   (  -2.070   10.047    0.674)   10.280   3.915   (   7.142    2.386   -8.218)   11.146   9.247   (   1.302    1.631    3.076)    3.717  10.746   (  -1.275    0.762   -0.481)    1.562  10.913   (   0.816   -2.977   -3.545)    4.701======================= Grid point 393 (58/126) =======================q-point: ( 0.31  0.15  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.827   (   3.134    2.746    6.330)    7.578   2.125   (   1.808    9.458    1.864)    9.808   2.572   (  -1.752   -2.115    4.062)    4.903   2.729   (   4.337   -0.027    3.869)    5.812   3.305   (  -6.895   12.955    0.104)   14.675   4.063   (   7.991   -0.184   -6.878)   10.545   9.282   (   1.076    1.183    3.543)    3.887  10.726   (  -1.398    0.344   -0.395)    1.493  10.895   (   1.426   -4.283   -2.182)    5.014======================= Grid point 394 (59/126) =======================q-point: ( 0.38  0.15  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.883   (   1.321    0.709    6.032)    6.216   2.209   (  -1.425    6.732    0.390)    6.892   2.537   (  -0.654   -0.633    3.292)    3.415   2.779   (   2.809   -2.335    1.999)    4.164   3.304   (  -8.735   14.700    0.143)   17.101   4.161   (   4.714   -2.796   -6.575)    8.560   9.304   (   0.337    0.561    3.764)    3.821  10.711   (  -0.512    0.043   -0.351)    0.622  10.881   (   2.034   -4.402   -1.190)    4.993======================= Grid point 405 (60/126) =======================q-point: ( 0.23  0.23  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.813   (   2.176    3.769    6.087)    7.483   2.164   (   4.824    8.355    1.116)    9.712   2.577   (  -1.357   -2.351    1.273)    2.998   2.685   (   1.160    2.010    7.491)    7.842   3.457   (   4.604    7.974    0.131)    9.209   3.969   (   1.869    3.237   -6.868)    7.819   9.279   (   0.863    1.495    3.178)    3.616  10.748   (  -0.375   -0.649    0.081)    0.754  10.852   (  -1.592   -2.757   -2.616)    4.121======================= Grid point 406 (61/126) =======================q-point: ( 0.31  0.23  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.872   (   1.387    1.852    5.575)    6.036   2.289   (   0.909    7.009   -1.069)    7.148   2.540   (  -0.810   -1.354   -0.279)    1.602   2.714   (   1.847   -1.263    6.608)    6.976   3.568   (  -3.613   12.809   -0.724)   13.328   4.038   (   5.095   -2.347   -5.891)    8.134   9.303   (   0.535    0.856    3.099)    3.259  10.725   (  -1.209   -0.393    0.421)    1.339  10.812   (   0.665   -3.708   -1.605)    4.095======================= Grid point 407 (62/126) =======================q-point: ( 0.38  0.23  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.894   (  -0.231    0.399    5.546)    5.565   2.332   (  -3.152    5.457   -1.939)    6.593   2.525   (   0.438   -0.760   -1.408)    1.659   2.724   (   1.613   -2.795    7.095)    7.794   3.598   (  -8.029   13.902   -1.086)   16.090   4.074   (   3.427   -5.940   -5.601)    8.855   9.312   (  -0.159    0.275    3.038)    3.054  10.711   (  -0.064    0.111    0.436)    0.454  10.800   (   1.969   -3.414   -1.077)    4.086======================= Grid point 420 (63/126) =======================q-point: ( 0.31  0.31  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.898   (   0.399    0.691    5.410)    5.468   2.398   (   2.113    3.660   -2.774)    5.055   2.511   (  -0.861   -1.492   -5.639)    5.896   2.691   (  -0.438   -0.758   11.556)   11.589   3.783   (   4.651    8.055   -1.497)    9.421   3.971   (  -2.342   -4.057   -5.346)    7.108   9.315   (   0.177    0.307    2.701)    2.725  10.715   (  -0.271   -0.469    1.402)    1.503  10.761   (  -0.737   -1.276   -1.687)    2.240======================= Grid point 543 (64/126) =======================q-point: ( 0.00  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.189   (   0.000    0.000   21.914)   21.914   1.189   (   0.000    0.000   21.914)   21.914   2.298   (   0.000    0.000   35.560)   35.560   2.786   (   0.000    0.000   -7.769)    7.769   2.786   (   0.000    0.000   -7.769)    7.769   4.128   (   0.000    0.000    7.442)    7.442   9.133   (   0.000    0.000    0.207)    0.207  10.755   (   0.000    0.000    0.058)    0.058  10.755   (   0.000    0.000    0.058)    0.058======================= Grid point 544 (65/126) =======================q-point: ( 0.08  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.281   (   7.461    4.308   20.043)   21.817   1.327   (  10.218    5.900   19.203)   22.538   2.328   (   2.837    1.638   32.387)   32.552   2.802   (   1.377    0.795   -7.277)    7.449   2.867   (   6.827    3.942   -6.004)    9.909   4.031   (  -8.167   -4.715    6.230)   11.302   9.146   (   1.093    0.631    0.487)    1.353  10.749   (  -0.548   -0.316    0.001)    0.632  10.778   (   1.881    1.086   -0.287)    2.191======================= Grid point 545 (66/126) =======================q-point: ( 0.15  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.494   (  10.387    5.997   16.234)   20.184   1.561   (   8.929    5.155   14.458)   17.757   2.430   (   5.716    3.300   23.951)   24.844   2.846   (   2.385    1.377   -5.922)    6.531   3.070   (   9.790    5.652   -0.291)   11.308   3.794   ( -11.071   -6.392    1.913)   12.926   9.182   (   2.068    1.194    1.265)    2.703  10.732   (  -0.873   -0.504   -0.132)    1.017  10.833   (   2.758    1.592   -0.974)    3.330======================= Grid point 546 (67/126) =======================q-point: ( 0.23  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.723   (   5.545    3.201    9.867)   11.763   1.727   (   9.632    5.561   12.349)   16.619   2.556   (   4.505    2.601   12.059)   13.133   2.904   (   2.486    1.435   -3.925)    4.863   3.191   (  -2.628   -1.517   10.085)   10.532   3.662   (   3.177    1.834   -5.853)    6.907   9.238   (   2.767    1.598    2.360)    3.972  10.712   (  -0.833   -0.481   -0.280)    1.002  10.895   (   2.454    1.417   -1.297)    3.117======================= Grid point 547 (68/126) =======================q-point: ( 0.31  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.833   (   4.396    2.538    5.766)    7.682   1.923   (   7.402    4.274    8.753)   12.234   2.606   (  -0.623   -0.360   -0.103)    0.727   2.951   (   1.507    0.870   -1.479)    2.284   3.048   (  -6.024   -3.478   10.890)   12.922   3.854   (  10.049    5.802   -3.330)   12.072   9.304   (   2.880    1.663    3.515)    4.839  10.696   (  -0.495   -0.286   -0.423)    0.711  10.941   (   1.571    0.907   -0.852)    2.005======================= Grid point 548 (69/126) =======================q-point: ( 0.38  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.928   (   3.897    2.250    3.231)    5.539   2.063   (   4.757    2.746    5.569)    7.822   2.539   (  -4.218   -2.435   -7.254)    8.737   2.971   (   0.264    0.152    1.106)    1.148   2.972   (  -1.177   -0.679    9.511)    9.607   4.051   (   6.761    3.904   -1.701)    7.990   9.364   (   2.167    1.251    4.468)    5.121  10.689   (  -0.140   -0.081   -0.562)    0.584  10.967   (   0.770    0.444    0.119)    0.897======================= Grid point 549 (70/126) =======================q-point: ( 0.46  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.999   (   1.856    1.072    2.669)    3.423   2.137   (   1.689    0.975    3.439)    3.954   2.459   (  -2.130   -1.230   -9.486)    9.800   2.963   (  -0.038   -0.022    7.291)    7.291   2.970   (  -0.133   -0.077    3.000)    3.004   4.155   (   2.256    1.302   -1.442)    2.977   9.398   (   0.797    0.460    5.008)    5.092  10.688   (  -0.001   -0.001   -0.657)    0.657  10.978   (   0.224    0.129    0.908)    0.944======================= Grid point 557 (71/126) =======================q-point: ( 0.08  0.08  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.439   (   6.199   10.737   17.071)   21.098   1.498   (   5.560    9.631   15.843)   19.356   2.396   (   2.932    5.078   26.719)   27.354   2.836   (   1.442    2.497   -6.142)    6.785   3.008   (   5.589    9.681   -2.638)   11.485   3.862   (  -6.580  -11.397    3.557)   13.632   9.170   (   1.035    1.793    1.012)    2.305  10.741   (  -0.249   -0.432   -0.116)    0.512  10.812   (   1.290    2.235   -0.773)    2.694======================= Grid point 558 (72/126) =======================q-point: ( 0.15  0.08  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.659   (   5.632    7.597   12.315)   15.527   1.675   (   6.800    8.198   12.351)   16.310   2.510   (   4.151    4.166   16.656)   17.663   2.892   (   1.651    2.828   -3.601)    4.867   3.199   (   4.056    6.238    5.633)    9.332   3.651   (  -3.817   -6.156   -3.429)    8.015   9.215   (   1.935    2.069    1.917)    3.420  10.731   (  -0.969    0.362   -0.272)    1.069  10.860   (   2.281    0.962   -1.211)    2.756======================= Grid point 559 (73/126) =======================q-point: ( 0.23  0.08  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.782   (   4.267    3.031    7.309)    8.990   1.881   (   5.693    9.049    9.288)   14.162   2.589   (   1.759    0.132    4.656)    4.979   2.932   (   0.744    0.287    1.620)    1.805   3.176   (  -7.389    2.033    7.208)   10.520   3.725   (  10.091    4.110   -4.730)   11.878   9.275   (   2.527    1.960    2.961)    4.358  10.716   (  -1.362    0.846   -0.378)    1.647  10.904   (   2.503   -0.565   -1.072)    2.781======================= Grid point 560 (74/126) =======================q-point: ( 0.31  0.08  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.878   (   4.111    2.140    4.034)    6.145   2.040   (   3.561    6.815    6.000)    9.753   2.561   (  -2.212   -3.881   -4.947)    6.665   2.931   (   0.714   -1.384    7.612)    7.770   3.108   (  -6.383    7.978    2.808)   10.596   3.945   (   9.669    3.394   -1.879)   10.419   9.336   (   2.474    1.393    3.896)    4.821  10.703   (  -1.105    0.775   -0.448)    1.422  10.933   (   2.210   -1.689   -0.275)    2.795======================= Grid point 561 (75/126) =======================q-point: ( 0.38  0.08  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.961   (   3.348    1.065    2.759)    4.467   2.137   (   1.338    4.304    3.424)    5.660   2.477   (  -2.423   -3.311   -8.708)    9.626   2.939   (   1.083   -0.526    7.648)    7.742   3.088   (  -6.151    9.106    2.520)   11.274   4.099   (   5.546    1.006   -1.174)    5.757   9.381   (   1.638    0.484    4.536)    4.847  10.694   (  -0.574    0.483   -0.515)    0.910  10.948   (   1.856   -2.340    0.692)    3.066======================= Grid point 562 (76/126) =======================q-point: ( 0.46  0.08  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.998   (   0.528   -0.916    2.768)    2.963   2.170   (  -1.267    2.193    2.386)    3.479   2.436   (   0.544   -0.943   -9.876)    9.936   2.947   (   0.388   -0.674    7.124)    7.166   3.079   (  -5.571    9.647    2.861)   11.502   4.152   (   0.901   -1.566   -1.134)    2.133   9.398   (   0.269   -0.467    4.764)    4.794  10.692   (  -0.227    0.393   -0.548)    0.712  10.952   (   1.535   -2.663    1.126)    3.273======================= Grid point 572 (77/126) =======================q-point: ( 0.15  0.15  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.759   (   2.529    4.381    7.770)    9.271   1.882   (   5.906   10.230    9.056)   14.885   2.571   (   0.802    1.389    6.696)    6.886   2.928   (   0.112    0.194    3.128)    3.136   3.270   (   0.671    1.162    6.279)    6.421   3.639   (   3.274    5.670   -5.906)    8.817   9.262   (   1.475    2.555    2.703)    4.001  10.740   (   0.204    0.354   -0.379)    0.557  10.871   (   0.131    0.227   -1.090)    1.121======================= Grid point 573 (78/126) =======================q-point: ( 0.23  0.15  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.842   (   2.925    3.036    4.462)    6.139   2.058   (   3.732    8.205    5.296)   10.455   2.552   (  -1.485   -3.584   -2.831)    4.803   2.914   (  -0.137   -1.100    9.016)    9.084   3.294   (  -3.253    8.497    1.125)    9.168   3.811   (   8.841    4.470   -1.831)   10.074   9.313   (   1.841    1.879    3.367)    4.273  10.737   (  -1.352    1.074   -0.337)    1.759  10.875   (   1.546   -2.198   -0.456)    2.726======================= Grid point 574 (79/126) =======================q-point: ( 0.31  0.15  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.919   (   3.036    2.028    2.815)    4.610   2.170   (   1.218    5.904    2.343)    6.468   2.476   (  -2.001   -4.217   -7.629)    8.943   2.916   (   0.982   -0.351    9.069)    9.129   3.315   (  -6.307   12.112    0.708)   13.675   3.989   (   8.018    1.121   -0.366)    8.105   9.358   (   1.666    0.832    3.789)    4.222  10.720   (  -1.595    0.793   -0.253)    1.799  10.880   (   2.257   -3.443    0.512)    4.148======================= Grid point 575 (80/126) =======================q-point: ( 0.38  0.15  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.975   (   1.630    0.501    2.749)    3.235   2.220   (  -0.865    3.870    0.563)    4.005   2.419   (  -0.445   -2.428   -9.598)    9.910   2.930   (   0.890   -0.524    8.129)    8.194   3.319   (  -7.750   13.376    1.246)   15.509   4.092   (   4.115   -1.717   -0.203)    4.464   9.383   (   0.902   -0.213    3.979)    4.086  10.705   (  -0.792    0.530   -0.219)    0.978  10.883   (   2.333   -3.902    1.260)    4.718======================= Grid point 586 (81/126) =======================q-point: ( 0.23  0.23  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.903   (   1.781    3.086    2.793)    4.527   2.192   (   3.031    5.249    1.554)    6.258   2.470   (  -2.357   -4.082   -6.975)    8.418   2.904   (   0.062    0.108    9.664)    9.665   3.457   (   4.045    7.006   -0.479)    8.104   3.902   (   2.774    4.805    0.497)    5.571   9.347   (   0.836    1.449    3.436)    3.822  10.749   (  -0.115   -0.199   -0.118)    0.259  10.830   (  -1.066   -1.847    0.434)    2.177======================= Grid point 587 (82/126) =======================q-point: ( 0.31  0.23  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.958   (   1.599    1.855    2.592)    3.566   2.265   (   0.414    3.671   -1.367)    3.940   2.405   (  -1.070   -2.939   -8.748)    9.290   2.912   (   0.625   -0.005    9.033)    9.054   3.553   (  -3.170   11.119   -0.820)   11.591   3.992   (   4.787   -0.758    1.459)    5.062   9.369   (   0.698    0.318    3.304)    3.392  10.730   (  -1.550    0.276    0.051)    1.575  10.810   (   1.431   -3.234    1.327)    3.777======================= Grid point 588 (83/126) =======================q-point: ( 0.38  0.23  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.983   (  -0.182    0.315    2.837)    2.860   2.285   (  -1.544    2.672   -2.625)    4.052   2.380   (   0.916   -1.587   -9.202)    9.383   2.919   (   0.280   -0.486    8.767)    8.785   3.579   (  -6.846   11.853   -0.772)   13.710   4.032   (   2.417   -4.192    1.570)    5.087   9.376   (   0.220   -0.382    3.221)    3.251  10.717   (  -0.431    0.745    0.138)    0.871  10.807   (   1.958   -3.394    1.679)    4.263======================= Grid point 601 (84/126) =======================q-point: ( 0.31  0.31  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.987   (   0.473    0.820    2.843)    2.997   2.317   (   0.980    1.697   -4.787)    5.172   2.360   (  -0.903   -1.564   -8.070)    8.270   2.914   (   0.039    0.068    9.133)    9.134   3.724   (   3.053    5.288   -3.562)    7.069   3.966   (  -0.818   -1.417    4.242)    4.546   9.372   (   0.031    0.053    2.874)    2.874  10.731   (  -0.101   -0.175   -0.091)    0.221  10.765   (  -0.642   -1.112    2.219)    2.564======================= Grid point 724 (85/126) =======================q-point: ( 0.00  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.620   (   0.000    0.000   20.660)   20.660   1.620   (  -0.000    0.000   20.660)   20.660   2.604   (   0.000    0.000  -10.468)   10.468   2.604   (   0.000    0.000  -10.468)   10.468   2.977   (   0.000    0.000   31.745)   31.745   4.173   (   0.000    0.000   -3.470)    3.470   9.139   (  -0.000   -0.000    0.335)    0.335  10.758   (   0.000    0.000    0.188)    0.188  10.758   (   0.000    0.000    0.188)    0.188======================= Grid point 725 (86/126) =======================q-point: ( 0.08  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.680   (   4.927    2.845   19.277)   20.099   1.716   (   6.920    3.995   18.882)   20.503   2.632   (   2.460    1.420   -9.683)   10.091   2.712   (   9.541    5.509   -8.274)   13.778   2.951   (  -2.202   -1.271   28.273)   28.387   4.064   (  -9.128   -5.270   -3.322)   11.051   9.157   (   1.511    0.872    0.528)    1.823  10.751   (  -0.626   -0.361    0.152)    0.738  10.776   (   1.556    0.899    0.077)    1.799======================= Grid point 726 (87/126) =======================q-point: ( 0.15  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.820   (   6.793    3.922   15.673)   17.526   1.862   (   4.922    2.842   14.698)   15.759   2.713   (   4.430    2.558   -7.322)    8.932   2.874   (  -4.523   -2.611   19.153)   19.853   3.024   (  16.552    9.556   -3.070)   19.358   3.791   ( -13.528   -7.810   -2.597)   15.835   9.207   (   2.814    1.625    1.073)    3.422  10.731   (  -0.985   -0.569    0.063)    1.139  10.824   (   2.476    1.430   -0.110)    2.861======================= Grid point 727 (88/126) =======================q-point: ( 0.23  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.934   (   1.707    0.986   11.006)   11.181   1.969   (   5.966    3.444   11.108)   13.071   2.739   (  -7.296   -4.213    5.318)    9.963   2.825   (   5.087    2.937   -3.868)    7.033   3.344   (   3.030    1.749    5.692)    6.681   3.567   (   2.228    1.286   -3.640)    4.458   9.282   (   3.659    2.112    1.842)    4.609  10.709   (  -0.946   -0.546   -0.043)    1.093  10.884   (   2.621    1.513   -0.052)    3.027======================= Grid point 728 (89/126) =======================q-point: ( 0.31  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.957   (   0.653    0.377    7.106)    7.146   2.087   (   4.246    2.451    7.018)    8.561   2.552   (  -8.535   -4.928   -5.361)   11.219   2.932   (   4.097    2.365   -0.453)    4.752   3.217   (  -6.201   -3.580    6.146)    9.437   3.830   (  12.471    7.200    0.434)   14.407   9.368   (   3.711    2.142    2.646)    5.036  10.691   (  -0.607   -0.350   -0.149)    0.716  10.940   (   2.241    1.294    0.445)    2.625======================= Grid point 729 (90/126) =======================q-point: ( 0.38  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.979   (   1.476    0.852    2.698)    3.191   2.164   (   2.505    1.446    4.186)    5.088   2.375   (  -6.565   -3.791   -8.802)   11.617   3.006   (   2.303    1.330    2.115)    3.398   3.116   (  -2.902   -1.675    4.960)    5.986   4.058   (   7.306    4.218    1.646)    8.596   9.444   (   2.766    1.597    3.319)    4.607  10.681   (  -0.238   -0.137   -0.251)    0.372  10.984   (   1.542    0.890    1.194)    2.144======================= Grid point 730 (91/126) =======================q-point: ( 0.46  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.017   (   1.366    0.789   -0.126)    1.583   2.202   (   0.830    0.479    2.759)    2.920   2.264   (  -2.831   -1.634   -9.069)    9.640   3.039   (   0.685    0.395    3.501)    3.589   3.075   (  -0.841   -0.485    3.798)    3.920   4.167   (   2.291    1.323    1.861)    3.235   9.488   (   1.020    0.589    3.718)    3.900  10.678   (  -0.043   -0.025   -0.319)    0.323  11.008   (   0.561    0.324    1.745)    1.861======================= Grid point 738 (92/126) =======================q-point: ( 0.08  0.08  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.784   (   4.024    6.969   16.647)   18.490   1.826   (   3.315    5.742   15.954)   17.277   2.693   (   2.643    4.577   -7.979)    9.571   2.899   (  -1.524   -2.639   22.866)   23.068   2.921   (   7.930   13.735   -5.484)   16.782   3.873   (  -7.583  -13.134   -2.845)   15.431   9.190   (   1.419    2.459    0.897)    2.977  10.741   (  -0.335   -0.580    0.065)    0.672  10.805   (   1.117    1.934   -0.047)    2.234======================= Grid point 739 (93/126) =======================q-point: ( 0.15  0.08  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.907   (   2.389    2.212   12.203)   12.630   1.934   (   4.372    6.201   12.489)   14.613   2.774   (  -1.937   -2.753    3.791)    5.070   2.825   (   0.664    3.649    2.035)    4.230   3.259   (  12.033   13.092    1.169)   17.820   3.598   (  -8.177  -10.348   -2.165)   13.365   9.251   (   2.579    2.764    1.546)    4.084  10.727   (  -1.012    0.182   -0.060)    1.030  10.850   (   2.226    1.029   -0.069)    2.453======================= Grid point 740 (94/126) =======================q-point: ( 0.23  0.08  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.941   (   1.031   -0.154    8.766)    8.828   2.061   (   3.378    5.184    7.994)   10.109   2.618   (  -6.584   -7.574   -2.004)   10.234   2.936   (   2.235    7.537   -0.366)    7.870   3.302   (  -7.433   -3.390    4.795)    9.473   3.676   (  13.480    7.724   -0.405)   15.542   9.329   (   3.275    2.466    2.284)    4.693  10.711   (  -1.386    0.745   -0.144)    1.580  10.899   (   2.854   -0.096    0.242)    2.866======================= Grid point 741 (95/126) =======================q-point: ( 0.31  0.08  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.957   (   1.399   -0.177    4.544)    4.758   2.149   (   1.975    3.473    4.560)    6.063   2.434   (  -6.123   -6.500   -7.775)   11.840   3.027   (   0.562    5.420    3.052)    6.246   3.180   (  -5.238    1.313    3.400)    6.381   3.951   (  10.277    4.896    1.694)   11.509   9.406   (   3.206    1.542    2.928)    4.608  10.696   (  -1.188    0.813   -0.187)    1.451  10.943   (   2.870   -1.031    0.904)    3.181======================= Grid point 742 (96/126) =======================q-point: ( 0.38  0.08  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.990   (   2.180    0.338    0.917)    2.389   2.199   (   0.769    1.806    2.486)    3.168   2.295   (  -3.802   -3.957   -8.901)   10.457   3.060   (  -0.120    0.956    4.687)    4.785   3.142   (  -4.587    7.049    2.272)    8.712   4.117   (   5.294    1.765    2.132)    5.973   9.463   (   2.237    0.241    3.375)    4.056  10.687   (  -0.716    0.642   -0.216)    0.985  10.974   (   2.356   -1.820    1.593)    3.377======================= Grid point 743 (97/126) =======================q-point: ( 0.46  0.08  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.017   (   0.418   -0.725   -0.359)    0.910   2.214   (  -0.384    0.664    1.738)    1.899   2.238   (   0.493   -0.857   -8.905)    8.960   3.060   (   0.018   -0.032    4.126)    4.126   3.141   (  -4.845    8.389    2.905)   10.114   4.170   (   0.534   -0.931    2.167)    2.418   9.483   (   0.580   -1.005    3.540)    3.725  10.684   (  -0.356    0.616   -0.230)    0.747  10.985   (   1.444   -2.504    1.885)    3.451======================= Grid point 753 (98/126) =======================q-point: ( 0.15  0.15  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.930   (   0.246    0.426    9.285)    9.298   2.060   (   3.387    5.867    8.086)   10.549   2.636   (  -5.032   -8.716   -0.556)   10.079   2.962   (   4.105    7.110    0.856)    8.254   3.370   (  -3.251   -5.630    2.883)    7.112   3.574   (   7.060   12.229   -0.388)   14.126   9.312   (   1.868    3.235    2.129)    4.300  10.734   (   0.161    0.278   -0.219)    0.389  10.865   (   0.396    0.686    0.258)    0.833======================= Grid point 754 (99/126) =======================q-point: ( 0.23  0.15  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.936   (   0.655   -0.119    5.566)    5.606   2.156   (   1.968    4.155    4.152)    6.195   2.459   (  -5.294   -8.007   -6.607)   11.653   3.058   (   0.724    2.720    5.890)    6.528   3.300   (  -3.094    4.654   -0.610)    5.622   3.827   (   9.595    7.251    2.550)   12.294   9.375   (   2.262    2.088    2.593)    4.025  10.732   (  -1.271    1.124   -0.218)    1.711  10.882   (   1.995   -1.401    0.841)    2.579======================= Grid point 755 (100/126) =======================q-point: ( 0.31  0.15  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.959   (   1.781    0.558    1.910)    2.671   2.209   (   0.548    2.292    1.396)    2.739   2.313   (  -3.744   -5.329   -8.405)   10.633   3.064   (  -0.251   -0.397    5.968)    5.986   3.325   (  -4.395   11.160    0.153)   11.996   4.026   (   7.430    2.764    3.086)    8.507   9.426   (   2.130    0.512    2.856)    3.599  10.716   (  -1.607    1.037   -0.080)    1.915  10.904   (   2.756   -2.697    1.500)    4.137======================= Grid point 756 (101/126) =======================q-point: ( 0.38  0.15  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.997   (   1.488    0.384   -0.068)    1.539   2.226   (  -0.213    0.810   -0.459)    0.955   2.230   (  -1.268   -2.393   -8.189)    8.625   3.058   (   0.016   -0.288    4.697)    4.706   3.350   (  -6.294   12.125    1.651)   13.761   4.130   (   3.301   -0.381    3.084)    4.533   9.455   (   1.384   -0.901    2.966)    3.395  10.703   (  -0.957    0.885    0.020)    1.304  10.919   (   2.495   -3.479    1.964)    4.710======================= Grid point 767 (102/126) =======================q-point: ( 0.23  0.23  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.945   (   0.704    1.219    2.272)    2.673   2.217   (   1.104    1.913    0.748)    2.332   2.319   (  -3.389   -5.871   -8.094)   10.558   3.063   (  -0.364   -0.630    6.497)    6.537   3.437   (   4.078    7.064   -1.300)    8.260   3.961   (   3.571    6.185    4.145)    8.258   9.409   (   0.785    1.360    2.618)    3.053  10.745   (  -0.023   -0.040   -0.156)    0.162  10.855   (  -0.547   -0.947    1.602)    1.940======================= Grid point 768 (103/126) =======================q-point: ( 0.31  0.23  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.980   (   1.403    1.464    0.152)    2.033   2.230   (  -1.491   -2.494   -6.633)    7.241   2.236   (  -0.401   -0.121   -3.188)    3.215   3.055   (  -0.047   -0.232    5.361)    5.367   3.542   (  -1.679    9.673   -0.147)    9.819   4.064   (   3.847    1.258    4.596)    6.124   9.428   (   0.814   -0.178    2.483)    2.619  10.732   (  -1.323    0.460    0.159)    1.409  10.848   (   1.533   -2.522    2.059)    3.598======================= Grid point 769 (104/126) =======================q-point: ( 0.38  0.23  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.004   (  -0.135    0.234   -0.450)    0.525   2.200   (   0.320   -0.556   -7.938)    7.963   2.233   (   0.083   -0.145   -2.704)    2.709   3.054   (   0.034   -0.060    4.899)    4.899   3.574   (  -5.490    9.504    0.392)   10.983   4.105   (   1.080   -1.874    4.602)    5.084   9.434   (   0.554   -0.959    2.408)    2.650  10.721   (  -0.550    0.951    0.331)    1.147  10.850   (   1.777   -3.081    2.199)    4.181======================= Grid point 782 (105/126) =======================q-point: ( 0.31  0.31  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.006   (   0.471    0.816   -0.603)    1.119   2.198   (  -0.299   -0.517   -7.189)    7.214   2.229   (  -0.332   -0.576   -4.202)    4.255   3.054   (   0.017    0.029    5.074)    5.075   3.679   (   2.147    3.719   -0.830)    4.374   4.079   (   0.183    0.316    5.888)    5.900   9.424   (  -0.103   -0.178    2.155)    2.165  10.734   (  -0.066   -0.114    0.407)    0.427  10.815   (  -0.470   -0.814    2.340)    2.522======================= Grid point 905 (106/126) =======================q-point: ( 0.00  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.997   (   0.000    0.000   15.929)   15.929   1.997   (  -0.000    0.000   15.929)   15.929   2.370   (   0.000    0.000  -12.200)   12.200   2.370   (   0.000    0.000  -12.200)   12.200   3.546   (   0.000    0.000   23.970)   23.970   3.989   (   0.000    0.000  -14.145)   14.145   9.145   (   0.000    0.000    0.168)    0.168  10.761   (   0.000    0.000    0.099)    0.099  10.761   (   0.000    0.000    0.099)    0.099======================= Grid point 906 (107/126) =======================q-point: ( 0.08  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.026   (   2.245    1.296   13.794)   14.036   2.044   (   2.903    1.676   11.981)   12.441   2.422   (   4.555    2.630  -10.235)   11.507   2.529   (  13.857    8.000   -8.174)   17.967   3.457   (  -7.743   -4.471   20.404)   22.277   3.902   (  -7.064   -4.078  -11.590)   14.172   9.165   (   1.738    1.004    0.234)    2.021  10.753   (  -0.666   -0.385    0.087)    0.774  10.778   (   1.460    0.843    0.072)    1.688======================= Grid point 907 (108/126) =======================q-point: ( 0.15  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.084   (   2.586    1.493    8.944)    9.429   2.098   (   1.755    1.013    7.262)    7.540   2.569   (   7.919    4.572   -5.758)   10.806   2.929   (  16.707    9.646   -3.003)   19.524   3.213   ( -10.686   -6.170   10.711)   16.340   3.702   (  -9.646   -5.569   -4.822)   12.137   9.222   (   3.227    1.863    0.424)    3.750  10.733   (  -1.046   -0.604    0.057)    1.209  10.824   (   2.439    1.408    0.039)    2.817======================= Grid point 908 (109/126) =======================q-point: ( 0.23  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.123   (   0.305    0.176    6.150)    6.160   2.139   (   2.270    1.310    4.930)    5.584   2.760   (   8.262    4.770   -2.098)    9.769   2.773   ( -16.911   -9.763   -0.588)   19.536   3.438   (   9.025    5.211    3.839)   11.106   3.507   (   0.946    0.546   -2.707)    2.919   9.308   (   4.151    2.396    0.695)    4.843  10.709   (  -1.012   -0.584    0.019)    1.168  10.885   (   2.824    1.630    0.102)    3.262======================= Grid point 909 (110/126) =======================q-point: ( 0.31  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.109   (  -1.524   -0.880    6.797)    7.021   2.186   (   1.835    1.059    2.647)    3.391   2.417   ( -13.089   -7.557   -6.496)   16.451   2.929   (   6.265    3.617   -0.023)    7.235   3.298   (  -6.648   -3.838    1.997)    7.932   3.846   (  13.786    7.959    0.653)   15.932   9.405   (   4.169    2.407    0.977)    4.912  10.689   (  -0.672   -0.388   -0.021)    0.777  10.949   (   2.684    1.550    0.318)    3.115======================= Grid point 910 (111/126) =======================q-point: ( 0.38  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.060   (  -2.400   -1.386    5.557)    6.210   2.201   (  -5.984   -3.455   -8.216)   10.736   2.222   (   1.219    0.704    1.505)    2.061   3.041   (   3.547    2.048    1.059)    4.231   3.178   (  -3.828   -2.210    1.463)    4.656   4.090   (   7.611    4.394    1.046)    8.850   9.490   (   3.108    1.794    1.217)    3.789  10.678   (  -0.293   -0.169   -0.058)    0.344  11.004   (   1.994    1.151    0.609)    2.381======================= Grid point 911 (112/126) =======================q-point: ( 0.46  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.021   (  -0.816   -0.471    0.776)    1.220   2.119   (  -1.621   -0.936   -4.788)    5.141   2.241   (   0.429    0.248    1.040)    1.152   3.093   (   1.097    0.634    1.536)    1.991   3.121   (  -1.215   -0.701    1.061)    1.759   4.202   (   2.324    1.342    1.129)    2.911   9.540   (   1.152    0.665    1.363)    1.904  10.674   (  -0.066   -0.038   -0.083)    0.112  11.036   (   0.748    0.432    0.811)    1.185======================= Grid point 919 (113/126) =======================q-point: ( 0.08  0.08  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.068   (   1.395    2.416    9.822)   10.211   2.087   (   1.344    2.328    8.495)    8.910   2.531   (   4.669    8.086   -6.860)   11.587   2.810   (  10.709   18.549   -4.176)   21.822   3.285   (  -7.093  -12.286   13.666)   19.699   3.761   (  -5.438   -9.419   -6.727)   12.789   9.203   (   1.634    2.830    0.364)    3.288  10.743   (  -0.387   -0.671    0.055)    0.776  10.806   (   1.091    1.890    0.053)    2.183======================= Grid point 920 (114/126) =======================q-point: ( 0.15  0.08  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.103   (   1.023    0.288    5.900)    5.995   2.136   (   1.402    2.592    6.319)    6.972   2.713   (   5.158    8.265   -2.950)   10.179   2.916   ( -10.702  -11.129    1.562)   15.519   3.284   (  14.361   14.290    1.434)   20.310   3.561   (  -6.961   -8.195   -1.348)   10.836   9.273   (   2.935    3.167    0.595)    4.359  10.727   (  -1.027    0.045    0.005)    1.028  10.851   (   2.274    1.190    0.085)    2.568======================= Grid point 921 (115/126) =======================q-point: ( 0.23  0.08  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.107   (  -0.025   -1.440    6.049)    6.218   2.181   (   0.978    2.239    3.646)    4.388   2.536   ( -12.066  -12.797   -4.470)   18.147   2.928   (   3.379   11.150   -0.281)   11.655   3.365   (  -6.480   -5.776    1.542)    8.817   3.684   (  14.887   10.464    0.556)   18.205   9.361   (   3.673    2.768    0.856)    4.678  10.710   (  -1.376    0.646   -0.027)    1.521  10.905   (   3.107    0.255    0.236)    3.126======================= Grid point 922 (116/126) =======================q-point: ( 0.31  0.08  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.073   (  -1.454   -2.518    6.736)    7.337   2.212   (   0.573    0.861    1.465)    1.794   2.265   (  -7.086   -6.807   -8.129)   12.753   3.071   (   0.757    9.977    1.275)   10.087   3.225   (  -4.558   -3.101    1.030)    5.608   3.984   (  10.389    5.863    1.107)   11.981   9.447   (   3.597    1.639    1.082)    4.098  10.694   (  -1.211    0.817   -0.040)    1.461  10.959   (   3.264   -0.630    0.499)    3.361======================= Grid point 923 (117/126) =======================q-point: ( 0.38  0.08  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.025   (  -1.018   -1.866    2.819)    3.530   2.136   (  -2.459   -2.815   -6.308)    7.333   2.234   (   0.551    0.202    0.848)    1.031   3.132   (  -1.516    0.715    2.449)    2.967   3.166   (  -2.898    7.049    0.065)    7.622   4.156   (   5.087    2.298    1.241)    5.719   9.510   (   2.574    0.105    1.239)    2.858  10.685   (  -0.787    0.741   -0.045)    1.082  11.000   (   2.637   -1.535    0.754)    3.143======================= Grid point 924 (118/126) =======================q-point: ( 0.46  0.08  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.008   (   0.407   -0.706   -0.157)    0.829   2.102   (   0.400   -0.697   -3.969)    4.049   2.239   (   0.262   -0.455    0.648)    0.834   3.114   (  -0.045    0.077    1.281)    1.284   3.185   (  -4.610    7.982    1.259)    9.304   4.210   (   0.277   -0.484    1.255)    1.374   9.533   (   0.760   -1.316    1.298)    1.999  10.682   (  -0.438    0.759   -0.047)    0.877  11.015   (   1.400   -2.428    0.858)    2.931======================= Grid point 934 (119/126) =======================q-point: ( 0.15  0.15  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.098   (  -0.491   -0.851    5.731)    5.815   2.185   (   1.241    2.149    3.991)    4.700   2.575   (  -8.956  -15.512   -3.723)   18.294   2.970   (   6.303   10.916    0.073)   12.605   3.396   (  -4.504   -7.801    0.388)    9.016   3.594   (   8.942   15.489    1.101)   17.918   9.342   (   2.094    3.626    0.809)    4.265  10.731   (   0.111    0.192   -0.070)    0.232  10.872   (   0.613    1.062    0.246)    1.250======================= Grid point 935 (120/126) =======================q-point: ( 0.23  0.15  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.068   (  -1.094   -2.445    6.775)    7.285   2.215   (   0.270    0.847    1.737)    1.951   2.304   (  -6.790   -9.702   -7.805)   14.183   3.158   (   2.022    9.219    3.754)   10.157   3.271   (  -3.147   -1.521   -1.917)    3.987   3.875   (   9.454    8.721    1.527)   12.952   9.411   (   2.477    2.206    0.973)    3.457  10.729   (  -1.163    1.101   -0.070)    1.603  10.897   (   2.229   -0.866    0.482)    2.440======================= Grid point 936 (121/126) =======================q-point: ( 0.31  0.15  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.021   (  -1.072   -2.378    4.579)    5.270   2.147   (  -2.820   -4.068   -7.423)    8.922   2.229   (  -0.025   -0.282    0.285)    0.402   3.154   (  -1.920   -1.753    2.723)    3.765   3.320   (  -2.070   12.119   -0.465)   12.303   4.079   (   6.899    3.800    1.632)    8.043   9.466   (   2.371    0.323    1.057)    2.616  10.716   (  -1.554    1.187   -0.000)    1.956  10.928   (   3.003   -2.274    0.720)    3.835======================= Grid point 937 (122/126) =======================q-point: ( 0.38  0.15  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.996   (   0.003   -0.891    0.340)    0.953   2.094   (  -0.625   -1.423   -4.268)    4.542   2.224   (   0.446   -1.170   -0.152)    1.261   3.120   (  -0.757   -0.364    1.568)    1.779   3.377   (  -5.157   11.732    0.810)   12.841   4.183   (   2.750    0.503    1.619)    3.230   9.496   (   1.654   -1.300    1.091)    2.369  10.704   (  -1.044    1.128    0.053)    1.537  10.949   (   2.583   -3.278    0.878)    4.265======================= Grid point 948 (123/126) =======================q-point: ( 0.23  0.23  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.018   (  -1.301   -2.253    5.193)    5.808   2.151   (  -2.679   -4.641   -7.765)    9.435   2.228   (  -0.368   -0.638   -0.007)    0.737   3.162   (  -1.405   -2.434    2.953)    4.077   3.412   (   5.150    8.920   -0.984)   10.346   4.029   (   3.920    6.790    2.010)    8.094   9.446   (   0.742    1.285    0.976)    1.776  10.743   (   0.078    0.134   -0.029)    0.158  10.881   (  -0.285   -0.494    0.753)    0.945======================= Grid point 949 (124/126) =======================q-point: ( 0.31  0.23  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.990   (  -0.217   -0.670    1.163)    1.360   2.094   (  -0.918   -1.414   -4.870)    5.154   2.207   (  -0.394   -1.782   -0.834)    2.006   3.129   (  -0.729   -0.739    1.895)    2.161   3.544   (  -0.507    9.606    0.174)    9.621   4.138   (   3.342    2.249    2.144)    4.564   9.463   (   0.867   -0.475    0.918)    1.349  10.735   (  -1.165    0.691    0.118)    1.360  10.880   (   1.570   -2.239    0.896)    2.878======================= Grid point 950 (125/126) =======================q-point: ( 0.38  0.23  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.987   (   0.078   -0.136   -0.664)    0.683   2.081   (  -0.055    0.093   -3.337)    3.339   2.192   (   1.041   -1.806   -1.170)    2.390   3.119   (  -0.081    0.139    1.617)    1.625   3.587   (  -4.917    8.512    0.504)    9.844   4.179   (   0.408   -0.711    2.148)    2.299   9.468   (   0.745   -1.290    0.887)    1.734  10.727   (  -0.668    1.155    0.194)    1.348  10.883   (   1.726   -2.992    0.937)    3.579======================= Grid point 963 (126/126) =======================q-point: ( 0.31  0.31  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.986   (   0.020    0.035   -0.737)    0.738   2.086   (   0.186    0.322   -3.261)    3.282   2.169   (  -0.894   -1.549   -1.587)    2.391   3.121   (  -0.063   -0.109    1.681)    1.686   3.678   (   2.036    3.526    0.211)    4.077   4.169   (   0.460    0.796    2.500)    2.664   9.454   (  -0.181   -0.314    0.796)    0.874  10.742   (   0.032    0.055    0.256)    0.264  10.849   (  -0.450   -0.779    0.953)    1.310=================== 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/16731   10.0    711.432    711.432    465.950      0.000      0.000      0.000 3/16731   20.0     82.945     82.945     76.669      0.000      0.000     -0.000 3/16731   30.0     36.588     36.588     39.702      0.000      0.000     -0.000 3/16731   40.0     24.536     24.536     28.518     -0.000      0.000     -0.000 3/16731   50.0     19.100     19.100     22.764     -0.000      0.000     -0.000 3/16731   60.0     15.878     15.878     19.087     -0.000      0.000     -0.000 3/16731   70.0     13.670     13.670     16.475     -0.000      0.000     -0.000 3/16731   80.0     12.031     12.031     14.502     -0.000      0.000     -0.000 3/16731   90.0     10.752     10.752     12.954     -0.000      0.000     -0.000 3/16731  100.0      9.722      9.722     11.703     -0.000      0.000     -0.000 3/16731  110.0      8.873      8.873     10.672     -0.000      0.000     -0.000 3/16731  120.0      8.160      8.160      9.807     -0.000      0.000     -0.000 3/16731  130.0      7.553      7.553      9.071     -0.000      0.000     -0.000 3/16731  140.0      7.029      7.029      8.437     -0.000      0.000     -0.000 3/16731  150.0      6.573      6.573      7.886     -0.000      0.000     -0.000 3/16731  160.0      6.172      6.172      7.402     -0.000      0.000     -0.000 3/16731  170.0      5.818      5.818      6.974     -0.000      0.000     -0.000 3/16731  180.0      5.501      5.501      6.593     -0.000      0.000     -0.000 3/16731  190.0      5.217      5.217      6.251     -0.000      0.000     -0.000 3/16731  200.0      4.961      4.961      5.943     -0.000      0.000     -0.000 3/16731  210.0      4.729      4.729      5.663     -0.000      0.000     -0.000 3/16731  220.0      4.518      4.518      5.409     -0.000      0.000     -0.000 3/16731  230.0      4.324      4.324      5.177     -0.000      0.000     -0.000 3/16731  240.0      4.147      4.147      4.964     -0.000      0.000     -0.000 3/16731  250.0      3.983      3.983      4.768     -0.000      0.000     -0.000 3/16731  260.0      3.832      3.832      4.586     -0.000      0.000     -0.000 3/16731  270.0      3.692      3.692      4.418     -0.000      0.000     -0.000 3/16731  280.0      3.562      3.562      4.262     -0.000      0.000     -0.000 3/16731  290.0      3.440      3.440      4.117     -0.000      0.000     -0.000 3/16731  300.0      3.327      3.327      3.981     -0.000      0.000     -0.000 3/16731  310.0      3.221      3.221      3.854     -0.000      0.000     -0.000 3/16731  320.0      3.121      3.121      3.735     -0.000      0.000     -0.000 3/16731  330.0      3.028      3.028      3.622     -0.000      0.000     -0.000 3/16731  340.0      2.939      2.939      3.517     -0.000      0.000     -0.000 3/16731  350.0      2.856      2.856      3.417     -0.000      0.000     -0.000 3/16731  360.0      2.778      2.778      3.323     -0.000      0.000     -0.000 3/16731  370.0      2.703      2.703      3.234     -0.000      0.000     -0.000 3/16731  380.0      2.633      2.633      3.150     -0.000      0.000     -0.000 3/16731  390.0      2.566      2.566      3.070     -0.000      0.000     -0.000 3/16731  400.0      2.502      2.502      2.994     -0.000      0.000     -0.000 3/16731  410.0      2.442      2.442      2.921     -0.000      0.000     -0.000 3/16731  420.0      2.384      2.384      2.852     -0.000      0.000     -0.000 3/16731  430.0      2.329      2.329      2.786     -0.000      0.000     -0.000 3/16731  440.0      2.277      2.277      2.724     -0.000      0.000     -0.000 3/16731  450.0      2.226      2.226      2.664     -0.000      0.000     -0.000 3/16731  460.0      2.178      2.178      2.606     -0.000      0.000     -0.000 3/16731  470.0      2.132      2.132      2.551     -0.000      0.000     -0.000 3/16731  480.0      2.088      2.088      2.498     -0.000      0.000     -0.000 3/16731  490.0      2.046      2.046      2.448     -0.000      0.000     -0.000 3/16731  500.0      2.005      2.005      2.399     -0.000      0.000     -0.000 3/16731  510.0      1.966      1.966      2.352     -0.000      0.000     -0.000 3/16731  520.0      1.929      1.929      2.307     -0.000      0.000     -0.000 3/16731  530.0      1.893      1.893      2.264     -0.000      0.000     -0.000 3/16731  540.0      1.858      1.858      2.223     -0.000      0.000     -0.000 3/16731  550.0      1.824      1.824      2.182     -0.000      0.000     -0.000 3/16731  560.0      1.792      1.792      2.144     -0.000      0.000     -0.000 3/16731  570.0      1.761      1.761      2.106     -0.000      0.000     -0.000 3/16731  580.0      1.730      1.730      2.070     -0.000      0.000     -0.000 3/16731  590.0      1.701      1.701      2.035     -0.000      0.000     -0.000 3/16731  600.0      1.673      1.673      2.002     -0.000      0.000     -0.000 3/16731  610.0      1.646      1.646      1.969     -0.000      0.000     -0.000 3/16731  620.0      1.619      1.619      1.938     -0.000      0.000     -0.000 3/16731  630.0      1.594      1.594      1.907     -0.000      0.000     -0.000 3/16731  640.0      1.569      1.569      1.877     -0.000      0.000     -0.000 3/16731  650.0      1.545      1.545      1.849     -0.000      0.000     -0.000 3/16731  660.0      1.522      1.522      1.821     -0.000      0.000     -0.000 3/16731  670.0      1.499      1.499      1.794     -0.000      0.000     -0.000 3/16731  680.0      1.477      1.477      1.768     -0.000      0.000     -0.000 3/16731  690.0      1.456      1.456      1.742     -0.000      0.000     -0.000 3/16731  700.0      1.435      1.435      1.717     -0.000      0.000     -0.000 3/16731  710.0      1.415      1.415      1.693     -0.000      0.000     -0.000 3/16731  720.0      1.396      1.396      1.670     -0.000      0.000     -0.000 3/16731  730.0      1.377      1.377      1.647     -0.000      0.000     -0.000 3/16731  740.0      1.358      1.358      1.625     -0.000      0.000     -0.000 3/16731  750.0      1.340      1.340      1.604     -0.000      0.000     -0.000 3/16731  760.0      1.323      1.323      1.583     -0.000      0.000     -0.000 3/16731  770.0      1.306      1.306      1.562     -0.000      0.000     -0.000 3/16731  780.0      1.289      1.289      1.542     -0.000      0.000     -0.000 3/16731  790.0      1.273      1.273      1.523     -0.000      0.000     -0.000 3/16731  800.0      1.257      1.257      1.504     -0.000      0.000     -0.000 3/16731  810.0      1.241      1.241      1.485     -0.000      0.000     -0.000 3/16731  820.0      1.226      1.226      1.467     -0.000      0.000     -0.000 3/16731  830.0      1.212      1.212      1.450     -0.000      0.000     -0.000 3/16731  840.0      1.197      1.197      1.433     -0.000      0.000     -0.000 3/16731  850.0      1.183      1.183      1.416     -0.000      0.000     -0.000 3/16731  860.0      1.170      1.170      1.400     -0.000      0.000     -0.000 3/16731  870.0      1.156      1.156      1.384     -0.000      0.000     -0.000 3/16731  880.0      1.143      1.143      1.368     -0.000      0.000     -0.000 3/16731  890.0      1.130      1.130      1.353     -0.000      0.000     -0.000 3/16731  900.0      1.118      1.118      1.338     -0.000      0.000     -0.000 3/16731  910.0      1.106      1.106      1.323     -0.000      0.000     -0.000 3/16731  920.0      1.094      1.094      1.309     -0.000      0.000     -0.000 3/16731  930.0      1.082      1.082      1.295     -0.000      0.000     -0.000 3/16731  940.0      1.070      1.070      1.281     -0.000      0.000     -0.000 3/16731  950.0      1.059      1.059      1.268     -0.000      0.000     -0.000 3/16731  960.0      1.048      1.048      1.254     -0.000      0.000     -0.000 3/16731  970.0      1.037      1.037      1.242     -0.000      0.000     -0.000 3/16731  980.0      1.027      1.027      1.229     -0.000      0.000     -0.000 3/16731  990.0      1.017      1.017      1.217     -0.000      0.000     -0.000 3/16731 1000.0      1.006      1.006      1.205     -0.000      0.000     -0.000 3/16731Thermal conductivity related properties were written into "kappa-m131311.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-07 23:45:42]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|