# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/abd14b19-1c6d-4ef8-8d25-bcb20182354f/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 CsMgCl3 / P6_3/mmc (194) / materials id 23004](https://mdr.nims.go.jp/datasets/4e4b6c16-d121-43ae-97b9-624979b6c7ab)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-08 05:15:34]-------------------------Compiled with OpenMP support (max 128 threads).Running in phonopy.load mode.Python version 3.14.2Spglib version 2.6.1Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".Unit of length: angstromSettings:  Supercell: [2 2 3]  Primitive matrix:    [ 1. -0.  0.]    [0. 1. 0.]    [0. 0. 1.]Spacegroup: P6_3/mmc (194)Number of symmetry operations in supercell: 288------------------------------ primitive cell ------------------------------Lattice vectors:  a    7.199645617763331    0.000000000000000    0.000000000000000  b   -3.599822808881666    6.235076003228353    0.000000000000000  c    0.000000000000000    0.000000000000000    6.147624349999998Atomic positions (fractional):   *1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305    2 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305   *3 Cl  0.15739105246350  0.31478210492700  0.25000000000000  35.453    4 Cl  0.68521789507300  0.84260894753650  0.25000000000000  35.453    5 Cl  0.84260894753650  0.68521789507300  0.75000000000000  35.453    6 Cl  0.15739105246350  0.84260894753650  0.25000000000000  35.453    7 Cl  0.84260894753650  0.15739105246350  0.75000000000000  35.453    8 Cl  0.31478210492700  0.15739105246350  0.75000000000000  35.453   *9 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905   10 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905-------------------------------- unit cell ---------------------------------Lattice vectors:  a    7.199645617763331    0.000000000000000    0.000000000000000  b   -3.599822808881666    6.235076003228353    0.000000000000000  c    0.000000000000000    0.000000000000000    6.147624349999999Atomic positions (fractional):   *1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305 > 1    2 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305 > 2   *3 Cl  0.15739105246350  0.31478210492700  0.25000000000000  35.453 > 3    4 Cl  0.68521789507300  0.84260894753650  0.25000000000000  35.453 > 4    5 Cl  0.84260894753650  0.68521789507300  0.75000000000000  35.453 > 5    6 Cl  0.15739105246350  0.84260894753650  0.25000000000000  35.453 > 6    7 Cl  0.84260894753650  0.15739105246350  0.75000000000000  35.453 > 7    8 Cl  0.31478210492700  0.15739105246350  0.75000000000000  35.453 > 8   *9 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905 > 9   10 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905 > 10-------------------------------- super cell --------------------------------Lattice vectors:  a   14.399291235526661    0.000000000000000    0.000000000000000  b   -7.199645617763332   12.470152006456706    0.000000000000000  c    0.000000000000000    0.000000000000000   18.442873049999996Atomic positions (fractional):   *1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305 > 1    2 Mg  0.50000000000000  0.00000000000000  0.00000000000000  24.305 > 1    3 Mg  0.00000000000000  0.50000000000000  0.00000000000000  24.305 > 1    4 Mg  0.50000000000000  0.50000000000000  0.00000000000000  24.305 > 1    5 Mg  0.00000000000000  0.00000000000000  0.33333333333333  24.305 > 1    6 Mg  0.50000000000000  0.00000000000000  0.33333333333333  24.305 > 1    7 Mg  0.00000000000000  0.50000000000000  0.33333333333333  24.305 > 1    8 Mg  0.50000000000000  0.50000000000000  0.33333333333333  24.305 > 1    9 Mg  0.00000000000000  0.00000000000000  0.66666666666667  24.305 > 1   10 Mg  0.50000000000000  0.00000000000000  0.66666666666667  24.305 > 1   11 Mg  0.00000000000000  0.50000000000000  0.66666666666667  24.305 > 1   12 Mg  0.50000000000000  0.50000000000000  0.66666666666667  24.305 > 1   13 Mg  0.00000000000000  0.00000000000000  0.16666666666667  24.305 > 2   14 Mg  0.50000000000000  0.00000000000000  0.16666666666667  24.305 > 2   15 Mg  0.00000000000000  0.50000000000000  0.16666666666667  24.305 > 2   16 Mg  0.50000000000000  0.50000000000000  0.16666666666667  24.305 > 2   17 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305 > 2   18 Mg  0.50000000000000  0.00000000000000  0.50000000000000  24.305 > 2   19 Mg  0.00000000000000  0.50000000000000  0.50000000000000  24.305 > 2   20 Mg  0.50000000000000  0.50000000000000  0.50000000000000  24.305 > 2   21 Mg  0.00000000000000  0.00000000000000  0.83333333333333  24.305 > 2   22 Mg  0.50000000000000  0.00000000000000  0.83333333333333  24.305 > 2   23 Mg  0.00000000000000  0.50000000000000  0.83333333333333  24.305 > 2   24 Mg  0.50000000000000  0.50000000000000  0.83333333333333  24.305 > 2  *25 Cl  0.07869552623175  0.15739105246350  0.08333333333333  35.453 > 3   26 Cl  0.57869552623175  0.15739105246350  0.08333333333333  35.453 > 3   27 Cl  0.07869552623175  0.65739105246350  0.08333333333333  35.453 > 3   28 Cl  0.57869552623175  0.65739105246350  0.08333333333333  35.453 > 3   29 Cl  0.07869552623175  0.15739105246350  0.41666666666667  35.453 > 3   30 Cl  0.57869552623175  0.15739105246350  0.41666666666667  35.453 > 3   31 Cl  0.07869552623175  0.65739105246350  0.41666666666667  35.453 > 3   32 Cl  0.57869552623175  0.65739105246350  0.41666666666667  35.453 > 3   33 Cl  0.07869552623175  0.15739105246350  0.75000000000000  35.453 > 3   34 Cl  0.57869552623175  0.15739105246350  0.75000000000000  35.453 > 3   35 Cl  0.07869552623175  0.65739105246350  0.75000000000000  35.453 > 3   36 Cl  0.57869552623175  0.65739105246350  0.75000000000000  35.453 > 3   37 Cl  0.34260894753650  0.42130447376825  0.08333333333333  35.453 > 4   38 Cl  0.84260894753650  0.42130447376825  0.08333333333333  35.453 > 4   39 Cl  0.34260894753650  0.92130447376825  0.08333333333333  35.453 > 4   40 Cl  0.84260894753650  0.92130447376825  0.08333333333333  35.453 > 4   41 Cl  0.34260894753650  0.42130447376825  0.41666666666667  35.453 > 4   42 Cl  0.84260894753650  0.42130447376825  0.41666666666667  35.453 > 4   43 Cl  0.34260894753650  0.92130447376825  0.41666666666667  35.453 > 4   44 Cl  0.84260894753650  0.92130447376825  0.41666666666667  35.453 > 4   45 Cl  0.34260894753650  0.42130447376825  0.75000000000000  35.453 > 4   46 Cl  0.84260894753650  0.42130447376825  0.75000000000000  35.453 > 4   47 Cl  0.34260894753650  0.92130447376825  0.75000000000000  35.453 > 4   48 Cl  0.84260894753650  0.92130447376825  0.75000000000000  35.453 > 4   49 Cl  0.42130447376825  0.34260894753650  0.25000000000000  35.453 > 5   50 Cl  0.92130447376825  0.34260894753650  0.25000000000000  35.453 > 5   51 Cl  0.42130447376825  0.84260894753650  0.25000000000000  35.453 > 5   52 Cl  0.92130447376825  0.84260894753650  0.25000000000000  35.453 > 5   53 Cl  0.42130447376825  0.34260894753650  0.58333333333333  35.453 > 5   54 Cl  0.92130447376825  0.34260894753650  0.58333333333333  35.453 > 5   55 Cl  0.42130447376825  0.84260894753650  0.58333333333333  35.453 > 5   56 Cl  0.92130447376825  0.84260894753650  0.58333333333333  35.453 > 5   57 Cl  0.42130447376825  0.34260894753650  0.91666666666667  35.453 > 5   58 Cl  0.92130447376825  0.34260894753650  0.91666666666667  35.453 > 5   59 Cl  0.42130447376825  0.84260894753650  0.91666666666667  35.453 > 5   60 Cl  0.92130447376825  0.84260894753650  0.91666666666667  35.453 > 5   61 Cl  0.07869552623175  0.42130447376825  0.08333333333333  35.453 > 6   62 Cl  0.57869552623175  0.42130447376825  0.08333333333333  35.453 > 6   63 Cl  0.07869552623175  0.92130447376825  0.08333333333333  35.453 > 6   64 Cl  0.57869552623175  0.92130447376825  0.08333333333333  35.453 > 6   65 Cl  0.07869552623175  0.42130447376825  0.41666666666667  35.453 > 6   66 Cl  0.57869552623175  0.42130447376825  0.41666666666667  35.453 > 6   67 Cl  0.07869552623175  0.92130447376825  0.41666666666667  35.453 > 6   68 Cl  0.57869552623175  0.92130447376825  0.41666666666667  35.453 > 6   69 Cl  0.07869552623175  0.42130447376825  0.75000000000000  35.453 > 6   70 Cl  0.57869552623175  0.42130447376825  0.75000000000000  35.453 > 6   71 Cl  0.07869552623175  0.92130447376825  0.75000000000000  35.453 > 6   72 Cl  0.57869552623175  0.92130447376825  0.75000000000000  35.453 > 6   73 Cl  0.42130447376825  0.07869552623175  0.25000000000000  35.453 > 7   74 Cl  0.92130447376825  0.07869552623175  0.25000000000000  35.453 > 7   75 Cl  0.42130447376825  0.57869552623175  0.25000000000000  35.453 > 7   76 Cl  0.92130447376825  0.57869552623175  0.25000000000000  35.453 > 7   77 Cl  0.42130447376825  0.07869552623175  0.58333333333333  35.453 > 7   78 Cl  0.92130447376825  0.07869552623175  0.58333333333333  35.453 > 7   79 Cl  0.42130447376825  0.57869552623175  0.58333333333333  35.453 > 7   80 Cl  0.92130447376825  0.57869552623175  0.58333333333333  35.453 > 7   81 Cl  0.42130447376825  0.07869552623175  0.91666666666667  35.453 > 7   82 Cl  0.92130447376825  0.07869552623175  0.91666666666667  35.453 > 7   83 Cl  0.42130447376825  0.57869552623175  0.91666666666667  35.453 > 7   84 Cl  0.92130447376825  0.57869552623175  0.91666666666667  35.453 > 7   85 Cl  0.15739105246350  0.07869552623175  0.25000000000000  35.453 > 8   86 Cl  0.65739105246350  0.07869552623175  0.25000000000000  35.453 > 8   87 Cl  0.15739105246350  0.57869552623175  0.25000000000000  35.453 > 8   88 Cl  0.65739105246350  0.57869552623175  0.25000000000000  35.453 > 8   89 Cl  0.15739105246350  0.07869552623175  0.58333333333333  35.453 > 8   90 Cl  0.65739105246350  0.07869552623175  0.58333333333333  35.453 > 8   91 Cl  0.15739105246350  0.57869552623175  0.58333333333333  35.453 > 8   92 Cl  0.65739105246350  0.57869552623175  0.58333333333333  35.453 > 8   93 Cl  0.15739105246350  0.07869552623175  0.91666666666667  35.453 > 8   94 Cl  0.65739105246350  0.07869552623175  0.91666666666667  35.453 > 8   95 Cl  0.15739105246350  0.57869552623175  0.91666666666667  35.453 > 8   96 Cl  0.65739105246350  0.57869552623175  0.91666666666667  35.453 > 8  *97 Cs  0.33333333333333  0.16666666666667  0.08333333333333 132.905 > 9   98 Cs  0.83333333333333  0.16666666666667  0.08333333333333 132.905 > 9   99 Cs  0.33333333333333  0.66666666666667  0.08333333333333 132.905 > 9  100 Cs  0.83333333333333  0.66666666666667  0.08333333333333 132.905 > 9  101 Cs  0.33333333333333  0.16666666666667  0.41666666666667 132.905 > 9  102 Cs  0.83333333333333  0.16666666666667  0.41666666666667 132.905 > 9  103 Cs  0.33333333333333  0.66666666666667  0.41666666666667 132.905 > 9  104 Cs  0.83333333333333  0.66666666666667  0.41666666666667 132.905 > 9  105 Cs  0.33333333333333  0.16666666666667  0.75000000000000 132.905 > 9  106 Cs  0.83333333333333  0.16666666666667  0.75000000000000 132.905 > 9  107 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905 > 9  108 Cs  0.83333333333333  0.66666666666667  0.75000000000000 132.905 > 9  109 Cs  0.16666666666667  0.33333333333333  0.25000000000000 132.905 > 10  110 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905 > 10  111 Cs  0.16666666666667  0.83333333333333  0.25000000000000 132.905 > 10  112 Cs  0.66666666666667  0.83333333333333  0.25000000000000 132.905 > 10  113 Cs  0.16666666666667  0.33333333333333  0.58333333333333 132.905 > 10  114 Cs  0.66666666666667  0.33333333333333  0.58333333333333 132.905 > 10  115 Cs  0.16666666666667  0.83333333333333  0.58333333333333 132.905 > 10  116 Cs  0.66666666666667  0.83333333333333  0.58333333333333 132.905 > 10  117 Cs  0.16666666666667  0.33333333333333  0.91666666666667 132.905 > 10  118 Cs  0.66666666666667  0.33333333333333  0.91666666666667 132.905 > 10  119 Cs  0.16666666666667  0.83333333333333  0.91666666666667 132.905 > 10  120 Cs  0.66666666666667  0.83333333333333  0.91666666666667 132.905 > 10----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            2.8226675    0.0000000    0.0000000           -0.0000000    2.8226675    0.0000000            0.0000000    0.0000000    2.8614575-------------------------- Born effective charges --------------------------    1 Mg    1.7761930    0.0000000    0.0000000            0.0000000    1.7761930    0.0000000            0.0000000    0.0000000    2.3440364    2 Mg    1.7761930    0.0000000    0.0000000            0.0000000    1.7761930    0.0000000            0.0000000    0.0000000    2.3440364    3 Cl   -0.8982808    0.0000000    0.0000000           -0.0000000   -1.2463529    0.0000000            0.0000000    0.0000000   -1.2059511    4 Cl   -1.1593348   -0.1507196    0.0000000           -0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    5 Cl   -0.8982808    0.0000000    0.0000000           -0.0000000   -1.2463529    0.0000000            0.0000000    0.0000000   -1.2059511    6 Cl   -1.1593348    0.1507196    0.0000000            0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    7 Cl   -1.1593348    0.1507196    0.0000000            0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    8 Cl   -1.1593348   -0.1507196    0.0000000           -0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    9 Cs    1.4407575    0.0000000    0.0000000            0.0000000    1.4407575    0.0000000            0.0000000    0.0000000    1.2738168   10 Cs    1.4407575    0.0000000    0.0000000            0.0000000    1.4407575    0.0000000            0.0000000    0.0000000    1.2738168----------------------------------------------------------------------------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: 360/360Permutation basis: 5880/5880Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 300Number of blocks in projector: 202Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 157Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 76Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 67Use standard eigh solver.Tree of FC basis block matrices:- (300, 293), data: False|-- (67, 64), data: True|-- (76, 76), data: True|-- (157, 153), data: True-----Solver_atoms: 1 -- 120 / 120Time (Solver_compr_matrix_reshape): 0.011Solver_block: 80 / 80 - Time: 0.222Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.238--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 360/360Permutation basis: 5880/5880Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 300Number of blocks in projector: 202Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 157Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 76Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 67Use standard eigh solver.Tree of FC basis block matrices:- (300, 293), data: False|-- (67, 64), data: True|-- (76, 76), data: True|-- (157, 153), data: TrueMax drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-08 05:15:39]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 05:15: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): [2 2 3]Primitive matrix:  [ 1. -0.  0.]  [0. 1. 0.]  [0. 0. 1.]Spacegroup: P6_3/mmc (194)------------------------------ primitive cell ------------------------------Lattice vectors:  a    7.199645617763331    0.000000000000000    0.000000000000000  b   -3.599822808881666    6.235076003228353    0.000000000000000  c    0.000000000000000    0.000000000000000    6.147624349999998Atomic positions (fractional):    1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305    2 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305    3 Cl  0.15739105246350  0.31478210492700  0.25000000000000  35.453    4 Cl  0.68521789507300  0.84260894753650  0.25000000000000  35.453    5 Cl  0.84260894753650  0.68521789507300  0.75000000000000  35.453    6 Cl  0.15739105246350  0.84260894753650  0.25000000000000  35.453    7 Cl  0.84260894753650  0.15739105246350  0.75000000000000  35.453    8 Cl  0.31478210492700  0.15739105246350  0.75000000000000  35.453    9 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905   10 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905-------------------------------- supercell ---------------------------------Lattice vectors:  a   14.399291235526661    0.000000000000000    0.000000000000000  b   -7.199645617763332   12.470152006456706    0.000000000000000  c    0.000000000000000    0.000000000000000   18.442873049999996Atomic positions (fractional):    1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305 > 1    2 Mg  0.50000000000000  0.00000000000000  0.00000000000000  24.305 > 1    3 Mg  0.00000000000000  0.50000000000000  0.00000000000000  24.305 > 1    4 Mg  0.50000000000000  0.50000000000000  0.00000000000000  24.305 > 1    5 Mg  0.00000000000000  0.00000000000000  0.33333333333333  24.305 > 1    6 Mg  0.50000000000000  0.00000000000000  0.33333333333333  24.305 > 1    7 Mg  0.00000000000000  0.50000000000000  0.33333333333333  24.305 > 1    8 Mg  0.50000000000000  0.50000000000000  0.33333333333333  24.305 > 1    9 Mg  0.00000000000000  0.00000000000000  0.66666666666667  24.305 > 1   10 Mg  0.50000000000000  0.00000000000000  0.66666666666667  24.305 > 1   11 Mg  0.00000000000000  0.50000000000000  0.66666666666667  24.305 > 1   12 Mg  0.50000000000000  0.50000000000000  0.66666666666667  24.305 > 1   13 Mg  0.00000000000000  0.00000000000000  0.16666666666667  24.305 > 13   14 Mg  0.50000000000000  0.00000000000000  0.16666666666667  24.305 > 13   15 Mg  0.00000000000000  0.50000000000000  0.16666666666667  24.305 > 13   16 Mg  0.50000000000000  0.50000000000000  0.16666666666667  24.305 > 13   17 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305 > 13   18 Mg  0.50000000000000  0.00000000000000  0.50000000000000  24.305 > 13   19 Mg  0.00000000000000  0.50000000000000  0.50000000000000  24.305 > 13   20 Mg  0.50000000000000  0.50000000000000  0.50000000000000  24.305 > 13   21 Mg  0.00000000000000  0.00000000000000  0.83333333333333  24.305 > 13   22 Mg  0.50000000000000  0.00000000000000  0.83333333333333  24.305 > 13   23 Mg  0.00000000000000  0.50000000000000  0.83333333333333  24.305 > 13   24 Mg  0.50000000000000  0.50000000000000  0.83333333333333  24.305 > 13   25 Cl  0.07869552623175  0.15739105246350  0.08333333333333  35.453 > 25   26 Cl  0.57869552623175  0.15739105246350  0.08333333333333  35.453 > 25   27 Cl  0.07869552623175  0.65739105246350  0.08333333333333  35.453 > 25   28 Cl  0.57869552623175  0.65739105246350  0.08333333333333  35.453 > 25   29 Cl  0.07869552623175  0.15739105246350  0.41666666666667  35.453 > 25   30 Cl  0.57869552623175  0.15739105246350  0.41666666666667  35.453 > 25   31 Cl  0.07869552623175  0.65739105246350  0.41666666666667  35.453 > 25   32 Cl  0.57869552623175  0.65739105246350  0.41666666666667  35.453 > 25   33 Cl  0.07869552623175  0.15739105246350  0.75000000000000  35.453 > 25   34 Cl  0.57869552623175  0.15739105246350  0.75000000000000  35.453 > 25   35 Cl  0.07869552623175  0.65739105246350  0.75000000000000  35.453 > 25   36 Cl  0.57869552623175  0.65739105246350  0.75000000000000  35.453 > 25   37 Cl  0.34260894753650  0.42130447376825  0.08333333333333  35.453 > 37   38 Cl  0.84260894753650  0.42130447376825  0.08333333333333  35.453 > 37   39 Cl  0.34260894753650  0.92130447376825  0.08333333333333  35.453 > 37   40 Cl  0.84260894753650  0.92130447376825  0.08333333333333  35.453 > 37   41 Cl  0.34260894753650  0.42130447376825  0.41666666666667  35.453 > 37   42 Cl  0.84260894753650  0.42130447376825  0.41666666666667  35.453 > 37   43 Cl  0.34260894753650  0.92130447376825  0.41666666666667  35.453 > 37   44 Cl  0.84260894753650  0.92130447376825  0.41666666666667  35.453 > 37   45 Cl  0.34260894753650  0.42130447376825  0.75000000000000  35.453 > 37   46 Cl  0.84260894753650  0.42130447376825  0.75000000000000  35.453 > 37   47 Cl  0.34260894753650  0.92130447376825  0.75000000000000  35.453 > 37   48 Cl  0.84260894753650  0.92130447376825  0.75000000000000  35.453 > 37   49 Cl  0.42130447376825  0.34260894753650  0.25000000000000  35.453 > 49   50 Cl  0.92130447376825  0.34260894753650  0.25000000000000  35.453 > 49   51 Cl  0.42130447376825  0.84260894753650  0.25000000000000  35.453 > 49   52 Cl  0.92130447376825  0.84260894753650  0.25000000000000  35.453 > 49   53 Cl  0.42130447376825  0.34260894753650  0.58333333333333  35.453 > 49   54 Cl  0.92130447376825  0.34260894753650  0.58333333333333  35.453 > 49   55 Cl  0.42130447376825  0.84260894753650  0.58333333333333  35.453 > 49   56 Cl  0.92130447376825  0.84260894753650  0.58333333333333  35.453 > 49   57 Cl  0.42130447376825  0.34260894753650  0.91666666666667  35.453 > 49   58 Cl  0.92130447376825  0.34260894753650  0.91666666666667  35.453 > 49   59 Cl  0.42130447376825  0.84260894753650  0.91666666666667  35.453 > 49   60 Cl  0.92130447376825  0.84260894753650  0.91666666666667  35.453 > 49   61 Cl  0.07869552623175  0.42130447376825  0.08333333333333  35.453 > 61   62 Cl  0.57869552623175  0.42130447376825  0.08333333333333  35.453 > 61   63 Cl  0.07869552623175  0.92130447376825  0.08333333333333  35.453 > 61   64 Cl  0.57869552623175  0.92130447376825  0.08333333333333  35.453 > 61   65 Cl  0.07869552623175  0.42130447376825  0.41666666666667  35.453 > 61   66 Cl  0.57869552623175  0.42130447376825  0.41666666666667  35.453 > 61   67 Cl  0.07869552623175  0.92130447376825  0.41666666666667  35.453 > 61   68 Cl  0.57869552623175  0.92130447376825  0.41666666666667  35.453 > 61   69 Cl  0.07869552623175  0.42130447376825  0.75000000000000  35.453 > 61   70 Cl  0.57869552623175  0.42130447376825  0.75000000000000  35.453 > 61   71 Cl  0.07869552623175  0.92130447376825  0.75000000000000  35.453 > 61   72 Cl  0.57869552623175  0.92130447376825  0.75000000000000  35.453 > 61   73 Cl  0.42130447376825  0.07869552623175  0.25000000000000  35.453 > 73   74 Cl  0.92130447376825  0.07869552623175  0.25000000000000  35.453 > 73   75 Cl  0.42130447376825  0.57869552623175  0.25000000000000  35.453 > 73   76 Cl  0.92130447376825  0.57869552623175  0.25000000000000  35.453 > 73   77 Cl  0.42130447376825  0.07869552623175  0.58333333333333  35.453 > 73   78 Cl  0.92130447376825  0.07869552623175  0.58333333333333  35.453 > 73   79 Cl  0.42130447376825  0.57869552623175  0.58333333333333  35.453 > 73   80 Cl  0.92130447376825  0.57869552623175  0.58333333333333  35.453 > 73   81 Cl  0.42130447376825  0.07869552623175  0.91666666666667  35.453 > 73   82 Cl  0.92130447376825  0.07869552623175  0.91666666666667  35.453 > 73   83 Cl  0.42130447376825  0.57869552623175  0.91666666666667  35.453 > 73   84 Cl  0.92130447376825  0.57869552623175  0.91666666666667  35.453 > 73   85 Cl  0.15739105246350  0.07869552623175  0.25000000000000  35.453 > 85   86 Cl  0.65739105246350  0.07869552623175  0.25000000000000  35.453 > 85   87 Cl  0.15739105246350  0.57869552623175  0.25000000000000  35.453 > 85   88 Cl  0.65739105246350  0.57869552623175  0.25000000000000  35.453 > 85   89 Cl  0.15739105246350  0.07869552623175  0.58333333333333  35.453 > 85   90 Cl  0.65739105246350  0.07869552623175  0.58333333333333  35.453 > 85   91 Cl  0.15739105246350  0.57869552623175  0.58333333333333  35.453 > 85   92 Cl  0.65739105246350  0.57869552623175  0.58333333333333  35.453 > 85   93 Cl  0.15739105246350  0.07869552623175  0.91666666666667  35.453 > 85   94 Cl  0.65739105246350  0.07869552623175  0.91666666666667  35.453 > 85   95 Cl  0.15739105246350  0.57869552623175  0.91666666666667  35.453 > 85   96 Cl  0.65739105246350  0.57869552623175  0.91666666666667  35.453 > 85   97 Cs  0.33333333333333  0.16666666666667  0.08333333333333 132.905 > 97   98 Cs  0.83333333333333  0.16666666666667  0.08333333333333 132.905 > 97   99 Cs  0.33333333333333  0.66666666666667  0.08333333333333 132.905 > 97  100 Cs  0.83333333333333  0.66666666666667  0.08333333333333 132.905 > 97  101 Cs  0.33333333333333  0.16666666666667  0.41666666666667 132.905 > 97  102 Cs  0.83333333333333  0.16666666666667  0.41666666666667 132.905 > 97  103 Cs  0.33333333333333  0.66666666666667  0.41666666666667 132.905 > 97  104 Cs  0.83333333333333  0.66666666666667  0.41666666666667 132.905 > 97  105 Cs  0.33333333333333  0.16666666666667  0.75000000000000 132.905 > 97  106 Cs  0.83333333333333  0.16666666666667  0.75000000000000 132.905 > 97  107 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905 > 97  108 Cs  0.83333333333333  0.66666666666667  0.75000000000000 132.905 > 97  109 Cs  0.16666666666667  0.33333333333333  0.25000000000000 132.905 > 109  110 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905 > 109  111 Cs  0.16666666666667  0.83333333333333  0.25000000000000 132.905 > 109  112 Cs  0.66666666666667  0.83333333333333  0.25000000000000 132.905 > 109  113 Cs  0.16666666666667  0.33333333333333  0.58333333333333 132.905 > 109  114 Cs  0.66666666666667  0.33333333333333  0.58333333333333 132.905 > 109  115 Cs  0.16666666666667  0.83333333333333  0.58333333333333 132.905 > 109  116 Cs  0.66666666666667  0.83333333333333  0.58333333333333 132.905 > 109  117 Cs  0.16666666666667  0.33333333333333  0.91666666666667 132.905 > 109  118 Cs  0.66666666666667  0.33333333333333  0.91666666666667 132.905 > 109  119 Cs  0.16666666666667  0.83333333333333  0.91666666666667 132.905 > 109  120 Cs  0.66666666666667  0.83333333333333  0.91666666666667 132.905 > 109----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            2.8226675    0.0000000    0.0000000           -0.0000000    2.8226675    0.0000000            0.0000000    0.0000000    2.8614575-------------------------- Born effective charges --------------------------    1 Mg    1.7761930    0.0000000    0.0000000            0.0000000    1.7761930    0.0000000            0.0000000    0.0000000    2.3440364    2 Mg    1.7761930    0.0000000    0.0000000            0.0000000    1.7761930    0.0000000            0.0000000    0.0000000    2.3440364    3 Cl   -0.8982808    0.0000000    0.0000000           -0.0000000   -1.2463529    0.0000000            0.0000000    0.0000000   -1.2059511    4 Cl   -1.1593348   -0.1507196    0.0000000           -0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    5 Cl   -0.8982808    0.0000000    0.0000000           -0.0000000   -1.2463529    0.0000000            0.0000000    0.0000000   -1.2059511    6 Cl   -1.1593348    0.1507196    0.0000000            0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    7 Cl   -1.1593348    0.1507196    0.0000000            0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    8 Cl   -1.1593348   -0.1507196    0.0000000           -0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    9 Cs    1.4407575    0.0000000    0.0000000            0.0000000    1.4407575    0.0000000            0.0000000    0.0000000    1.2738168   10 Cs    1.4407575    0.0000000    0.0000000            0.0000000    1.4407575    0.0000000            0.0000000    0.0000000    1.2738168----------------------------------------------------------------------------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[ 25, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [-0.0050  0.0087  0.0000]    [ 0.0050 -0.0087  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.00000001 (yyy) 0.00000001 (yyy) 0.00000001 (yyy)fc3 was written into "fc3.hdf5".Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330----------- None of ph-ph interaction calculation was performed. -----------Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 05:15:46]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 05:15: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): [2 2 3]Primitive matrix:  [ 1. -0.  0.]  [0. 1. 0.]  [0. 0. 1.]Spacegroup: P6_3/mmc (194)------------------------------ primitive cell ------------------------------Lattice vectors:  a    7.199645617763331    0.000000000000000    0.000000000000000  b   -3.599822808881666    6.235076003228353    0.000000000000000  c    0.000000000000000    0.000000000000000    6.147624349999998Atomic positions (fractional):    1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305    2 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305    3 Cl  0.15739105246350  0.31478210492700  0.25000000000000  35.453    4 Cl  0.68521789507300  0.84260894753650  0.25000000000000  35.453    5 Cl  0.84260894753650  0.68521789507300  0.75000000000000  35.453    6 Cl  0.15739105246350  0.84260894753650  0.25000000000000  35.453    7 Cl  0.84260894753650  0.15739105246350  0.75000000000000  35.453    8 Cl  0.31478210492700  0.15739105246350  0.75000000000000  35.453    9 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905   10 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905-------------------------------- supercell ---------------------------------Lattice vectors:  a   14.399291235526661    0.000000000000000    0.000000000000000  b   -7.199645617763332   12.470152006456706    0.000000000000000  c    0.000000000000000    0.000000000000000   18.442873049999996Atomic positions (fractional):    1 Mg  0.00000000000000  0.00000000000000  0.00000000000000  24.305 > 1    2 Mg  0.50000000000000  0.00000000000000  0.00000000000000  24.305 > 1    3 Mg  0.00000000000000  0.50000000000000  0.00000000000000  24.305 > 1    4 Mg  0.50000000000000  0.50000000000000  0.00000000000000  24.305 > 1    5 Mg  0.00000000000000  0.00000000000000  0.33333333333333  24.305 > 1    6 Mg  0.50000000000000  0.00000000000000  0.33333333333333  24.305 > 1    7 Mg  0.00000000000000  0.50000000000000  0.33333333333333  24.305 > 1    8 Mg  0.50000000000000  0.50000000000000  0.33333333333333  24.305 > 1    9 Mg  0.00000000000000  0.00000000000000  0.66666666666667  24.305 > 1   10 Mg  0.50000000000000  0.00000000000000  0.66666666666667  24.305 > 1   11 Mg  0.00000000000000  0.50000000000000  0.66666666666667  24.305 > 1   12 Mg  0.50000000000000  0.50000000000000  0.66666666666667  24.305 > 1   13 Mg  0.00000000000000  0.00000000000000  0.16666666666667  24.305 > 13   14 Mg  0.50000000000000  0.00000000000000  0.16666666666667  24.305 > 13   15 Mg  0.00000000000000  0.50000000000000  0.16666666666667  24.305 > 13   16 Mg  0.50000000000000  0.50000000000000  0.16666666666667  24.305 > 13   17 Mg  0.00000000000000  0.00000000000000  0.50000000000000  24.305 > 13   18 Mg  0.50000000000000  0.00000000000000  0.50000000000000  24.305 > 13   19 Mg  0.00000000000000  0.50000000000000  0.50000000000000  24.305 > 13   20 Mg  0.50000000000000  0.50000000000000  0.50000000000000  24.305 > 13   21 Mg  0.00000000000000  0.00000000000000  0.83333333333333  24.305 > 13   22 Mg  0.50000000000000  0.00000000000000  0.83333333333333  24.305 > 13   23 Mg  0.00000000000000  0.50000000000000  0.83333333333333  24.305 > 13   24 Mg  0.50000000000000  0.50000000000000  0.83333333333333  24.305 > 13   25 Cl  0.07869552623175  0.15739105246350  0.08333333333333  35.453 > 25   26 Cl  0.57869552623175  0.15739105246350  0.08333333333333  35.453 > 25   27 Cl  0.07869552623175  0.65739105246350  0.08333333333333  35.453 > 25   28 Cl  0.57869552623175  0.65739105246350  0.08333333333333  35.453 > 25   29 Cl  0.07869552623175  0.15739105246350  0.41666666666667  35.453 > 25   30 Cl  0.57869552623175  0.15739105246350  0.41666666666667  35.453 > 25   31 Cl  0.07869552623175  0.65739105246350  0.41666666666667  35.453 > 25   32 Cl  0.57869552623175  0.65739105246350  0.41666666666667  35.453 > 25   33 Cl  0.07869552623175  0.15739105246350  0.75000000000000  35.453 > 25   34 Cl  0.57869552623175  0.15739105246350  0.75000000000000  35.453 > 25   35 Cl  0.07869552623175  0.65739105246350  0.75000000000000  35.453 > 25   36 Cl  0.57869552623175  0.65739105246350  0.75000000000000  35.453 > 25   37 Cl  0.34260894753650  0.42130447376825  0.08333333333333  35.453 > 37   38 Cl  0.84260894753650  0.42130447376825  0.08333333333333  35.453 > 37   39 Cl  0.34260894753650  0.92130447376825  0.08333333333333  35.453 > 37   40 Cl  0.84260894753650  0.92130447376825  0.08333333333333  35.453 > 37   41 Cl  0.34260894753650  0.42130447376825  0.41666666666667  35.453 > 37   42 Cl  0.84260894753650  0.42130447376825  0.41666666666667  35.453 > 37   43 Cl  0.34260894753650  0.92130447376825  0.41666666666667  35.453 > 37   44 Cl  0.84260894753650  0.92130447376825  0.41666666666667  35.453 > 37   45 Cl  0.34260894753650  0.42130447376825  0.75000000000000  35.453 > 37   46 Cl  0.84260894753650  0.42130447376825  0.75000000000000  35.453 > 37   47 Cl  0.34260894753650  0.92130447376825  0.75000000000000  35.453 > 37   48 Cl  0.84260894753650  0.92130447376825  0.75000000000000  35.453 > 37   49 Cl  0.42130447376825  0.34260894753650  0.25000000000000  35.453 > 49   50 Cl  0.92130447376825  0.34260894753650  0.25000000000000  35.453 > 49   51 Cl  0.42130447376825  0.84260894753650  0.25000000000000  35.453 > 49   52 Cl  0.92130447376825  0.84260894753650  0.25000000000000  35.453 > 49   53 Cl  0.42130447376825  0.34260894753650  0.58333333333333  35.453 > 49   54 Cl  0.92130447376825  0.34260894753650  0.58333333333333  35.453 > 49   55 Cl  0.42130447376825  0.84260894753650  0.58333333333333  35.453 > 49   56 Cl  0.92130447376825  0.84260894753650  0.58333333333333  35.453 > 49   57 Cl  0.42130447376825  0.34260894753650  0.91666666666667  35.453 > 49   58 Cl  0.92130447376825  0.34260894753650  0.91666666666667  35.453 > 49   59 Cl  0.42130447376825  0.84260894753650  0.91666666666667  35.453 > 49   60 Cl  0.92130447376825  0.84260894753650  0.91666666666667  35.453 > 49   61 Cl  0.07869552623175  0.42130447376825  0.08333333333333  35.453 > 61   62 Cl  0.57869552623175  0.42130447376825  0.08333333333333  35.453 > 61   63 Cl  0.07869552623175  0.92130447376825  0.08333333333333  35.453 > 61   64 Cl  0.57869552623175  0.92130447376825  0.08333333333333  35.453 > 61   65 Cl  0.07869552623175  0.42130447376825  0.41666666666667  35.453 > 61   66 Cl  0.57869552623175  0.42130447376825  0.41666666666667  35.453 > 61   67 Cl  0.07869552623175  0.92130447376825  0.41666666666667  35.453 > 61   68 Cl  0.57869552623175  0.92130447376825  0.41666666666667  35.453 > 61   69 Cl  0.07869552623175  0.42130447376825  0.75000000000000  35.453 > 61   70 Cl  0.57869552623175  0.42130447376825  0.75000000000000  35.453 > 61   71 Cl  0.07869552623175  0.92130447376825  0.75000000000000  35.453 > 61   72 Cl  0.57869552623175  0.92130447376825  0.75000000000000  35.453 > 61   73 Cl  0.42130447376825  0.07869552623175  0.25000000000000  35.453 > 73   74 Cl  0.92130447376825  0.07869552623175  0.25000000000000  35.453 > 73   75 Cl  0.42130447376825  0.57869552623175  0.25000000000000  35.453 > 73   76 Cl  0.92130447376825  0.57869552623175  0.25000000000000  35.453 > 73   77 Cl  0.42130447376825  0.07869552623175  0.58333333333333  35.453 > 73   78 Cl  0.92130447376825  0.07869552623175  0.58333333333333  35.453 > 73   79 Cl  0.42130447376825  0.57869552623175  0.58333333333333  35.453 > 73   80 Cl  0.92130447376825  0.57869552623175  0.58333333333333  35.453 > 73   81 Cl  0.42130447376825  0.07869552623175  0.91666666666667  35.453 > 73   82 Cl  0.92130447376825  0.07869552623175  0.91666666666667  35.453 > 73   83 Cl  0.42130447376825  0.57869552623175  0.91666666666667  35.453 > 73   84 Cl  0.92130447376825  0.57869552623175  0.91666666666667  35.453 > 73   85 Cl  0.15739105246350  0.07869552623175  0.25000000000000  35.453 > 85   86 Cl  0.65739105246350  0.07869552623175  0.25000000000000  35.453 > 85   87 Cl  0.15739105246350  0.57869552623175  0.25000000000000  35.453 > 85   88 Cl  0.65739105246350  0.57869552623175  0.25000000000000  35.453 > 85   89 Cl  0.15739105246350  0.07869552623175  0.58333333333333  35.453 > 85   90 Cl  0.65739105246350  0.07869552623175  0.58333333333333  35.453 > 85   91 Cl  0.15739105246350  0.57869552623175  0.58333333333333  35.453 > 85   92 Cl  0.65739105246350  0.57869552623175  0.58333333333333  35.453 > 85   93 Cl  0.15739105246350  0.07869552623175  0.91666666666667  35.453 > 85   94 Cl  0.65739105246350  0.07869552623175  0.91666666666667  35.453 > 85   95 Cl  0.15739105246350  0.57869552623175  0.91666666666667  35.453 > 85   96 Cl  0.65739105246350  0.57869552623175  0.91666666666667  35.453 > 85   97 Cs  0.33333333333333  0.16666666666667  0.08333333333333 132.905 > 97   98 Cs  0.83333333333333  0.16666666666667  0.08333333333333 132.905 > 97   99 Cs  0.33333333333333  0.66666666666667  0.08333333333333 132.905 > 97  100 Cs  0.83333333333333  0.66666666666667  0.08333333333333 132.905 > 97  101 Cs  0.33333333333333  0.16666666666667  0.41666666666667 132.905 > 97  102 Cs  0.83333333333333  0.16666666666667  0.41666666666667 132.905 > 97  103 Cs  0.33333333333333  0.66666666666667  0.41666666666667 132.905 > 97  104 Cs  0.83333333333333  0.66666666666667  0.41666666666667 132.905 > 97  105 Cs  0.33333333333333  0.16666666666667  0.75000000000000 132.905 > 97  106 Cs  0.83333333333333  0.16666666666667  0.75000000000000 132.905 > 97  107 Cs  0.33333333333333  0.66666666666667  0.75000000000000 132.905 > 97  108 Cs  0.83333333333333  0.66666666666667  0.75000000000000 132.905 > 97  109 Cs  0.16666666666667  0.33333333333333  0.25000000000000 132.905 > 109  110 Cs  0.66666666666667  0.33333333333333  0.25000000000000 132.905 > 109  111 Cs  0.16666666666667  0.83333333333333  0.25000000000000 132.905 > 109  112 Cs  0.66666666666667  0.83333333333333  0.25000000000000 132.905 > 109  113 Cs  0.16666666666667  0.33333333333333  0.58333333333333 132.905 > 109  114 Cs  0.66666666666667  0.33333333333333  0.58333333333333 132.905 > 109  115 Cs  0.16666666666667  0.83333333333333  0.58333333333333 132.905 > 109  116 Cs  0.66666666666667  0.83333333333333  0.58333333333333 132.905 > 109  117 Cs  0.16666666666667  0.33333333333333  0.91666666666667 132.905 > 109  118 Cs  0.66666666666667  0.33333333333333  0.91666666666667 132.905 > 109  119 Cs  0.16666666666667  0.83333333333333  0.91666666666667 132.905 > 109  120 Cs  0.66666666666667  0.83333333333333  0.91666666666667 132.905 > 109----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            2.8226675    0.0000000    0.0000000           -0.0000000    2.8226675    0.0000000            0.0000000    0.0000000    2.8614575-------------------------- Born effective charges --------------------------    1 Mg    1.7761930    0.0000000    0.0000000            0.0000000    1.7761930    0.0000000            0.0000000    0.0000000    2.3440364    2 Mg    1.7761930    0.0000000    0.0000000            0.0000000    1.7761930    0.0000000            0.0000000    0.0000000    2.3440364    3 Cl   -0.8982808    0.0000000    0.0000000           -0.0000000   -1.2463529    0.0000000            0.0000000    0.0000000   -1.2059511    4 Cl   -1.1593348   -0.1507196    0.0000000           -0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    5 Cl   -0.8982808    0.0000000    0.0000000           -0.0000000   -1.2463529    0.0000000            0.0000000    0.0000000   -1.2059511    6 Cl   -1.1593348    0.1507196    0.0000000            0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    7 Cl   -1.1593348    0.1507196    0.0000000            0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    8 Cl   -1.1593348   -0.1507196    0.0000000           -0.1507196   -0.9852988    0.0000000            0.0000000    0.0000000   -1.2059511    9 Cs    1.4407575    0.0000000    0.0000000            0.0000000    1.4407575    0.0000000            0.0000000    0.0000000    1.2738168   10 Cs    1.4407575    0.0000000    0.0000000            0.0000000    1.4407575    0.0000000            0.0000000    0.0000000    1.2738168----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: 0.00000001 (yyy) 0.00000001 (yyy) 0.00000001 (yyy)Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330Length for sampling mesh generation: 50.00Generating grid system ... [ 8 8 8 ]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.64, Number of G-points: 277, Lambda: 0.11Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/50) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 50Calculating 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   0.951   (   0.000    0.000    0.000)    0.000   1.019   (   0.000    0.000    0.000)    0.000   1.431   (   0.000    0.000    0.000)    0.000   1.431   (   0.000    0.000    0.000)    0.000   1.850   (   0.000    0.000    0.000)    0.000   1.961   (   0.000    0.000    0.000)    0.000   1.961   (   0.000    0.000    0.000)    0.000   3.005   (   0.000    0.000    0.000)    0.000   3.723   (   0.000    0.000    0.000)    0.000   3.723   (   0.000    0.000    0.000)    0.000   3.750   (   0.000    0.000    0.000)    0.000   3.750   (   0.000    0.000    0.000)    0.000   3.831   (   0.000    0.000    0.000)    0.000   3.831   (   0.000    0.000    0.000)    0.000   4.805   (   0.000    0.000    0.000)    0.000   4.805   (   0.000    0.000    0.000)    0.000   5.501   (   0.000    0.000    0.000)    0.000   5.501   (   0.000    0.000    0.000)    0.000   6.089   (   0.000    0.000    0.000)    0.000   6.472   (   0.000    0.000    0.000)    0.000   6.784   (   0.000    0.000    0.000)    0.000   6.784   (   0.000    0.000    0.000)    0.000   7.411   (   0.000    0.000    0.000)    0.000   7.461   (   0.000    0.000    0.000)    0.000   9.314   (   0.000    0.000    0.000)    0.000   9.314   (   0.000    0.000    0.000)    0.000  11.451   (   0.000    0.000    0.000)    0.000======================= Grid point 1 (2/50) =======================q-point: ( 0.12  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.165   (   7.042    4.066    0.000)    8.132   0.337   (  14.156    8.173    0.000)   16.346   0.621   (  26.009   15.016    0.000)   30.033   0.946   (  -0.452   -0.261    0.000)    0.522   1.000   (  -1.565   -0.904    0.000)    1.808   1.424   (  -0.849   -0.490    0.000)    0.980   1.439   (   0.577    0.333    0.000)    0.666   1.841   (   0.036    0.021    0.000)    0.042   2.017   (   3.836    2.215    0.000)    4.430   2.599   (  -2.752   -1.589    0.000)    3.178   3.033   (   2.336    1.349    0.000)    2.697   3.720   (  -0.248   -0.143    0.000)    0.286   3.727   (   0.308    0.178    0.000)    0.355   3.757   (   0.535    0.309    0.000)    0.618   3.792   (   3.401    1.964    0.000)    3.927   3.834   (   0.244    0.141    0.000)    0.282   3.838   (   0.530    0.306    0.000)    0.612   4.796   (  -0.800   -0.462    0.000)    0.924   4.903   (   0.363    0.210    0.000)    0.419   5.492   (  -0.813   -0.469    0.000)    0.939   5.555   (   4.460    2.575    0.000)    5.150   6.086   (  -0.238   -0.137    0.000)    0.275   6.463   (  -0.854   -0.493    0.000)    0.987   6.784   (   0.005    0.003    0.000)    0.006   6.786   (   0.155    0.089    0.000)    0.178   7.399   (  -0.999   -0.577    0.000)    1.153   7.460   (  -0.078   -0.045    0.000)    0.090   9.316   (   0.160    0.092    0.000)    0.185  10.501   (  -1.387   -0.801    0.000)    1.602  11.407   (  -3.559   -2.055    0.000)    4.109======================= Grid point 2 (3/50) =======================q-point: ( 0.25  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 105Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.323   (   6.477    3.740    0.000)    7.479   0.647   (  12.491    7.212    0.000)   14.423   0.924   (  -1.846   -1.066    0.000)    2.131   0.957   (  -1.602   -0.925    0.000)    1.849   1.185   (  22.092   12.755    0.000)   25.510   1.383   (  -2.760   -1.594    0.000)    3.187   1.457   (   1.133    0.654    0.000)    1.308   1.868   (   2.148    1.240    0.000)    2.480   2.106   (   3.447    1.990    0.000)    3.980   2.503   (  -5.553   -3.206    0.000)    6.412   3.105   (   3.571    2.062    0.000)    4.123   3.713   (  -0.325   -0.188    0.000)    0.376   3.736   (   0.462    0.267    0.000)    0.533   3.773   (   0.747    0.431    0.000)    0.862   3.841   (   0.321    0.185    0.000)    0.370   3.852   (   0.689    0.398    0.000)    0.796   3.892   (   4.791    2.766    0.000)    5.532   4.771   (  -1.191   -0.687    0.000)    1.375   4.913   (   0.455    0.263    0.000)    0.525   5.467   (  -1.185   -0.684    0.000)    1.368   5.693   (   6.969    4.023    0.000)    8.047   6.079   (  -0.320   -0.185    0.000)    0.369   6.433   (  -1.649   -0.952    0.000)    1.904   6.784   (   0.007    0.004    0.000)    0.008   6.790   (   0.217    0.125    0.000)    0.250   7.369   (  -1.449   -0.837    0.000)    1.674   7.457   (  -0.111   -0.064    0.000)    0.128   9.320   (   0.230    0.133    0.000)    0.266  10.456   (  -2.367   -1.366    0.000)    2.733  11.304   (  -4.876   -2.815    0.000)    5.630======================= Grid point 3 (4/50) =======================q-point: ( 0.38  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.458   (   4.861    2.806    0.000)    5.613   0.860   (  -3.044   -1.758    0.000)    3.515   0.908   (   9.872    5.700    0.000)   11.399   0.937   (  -0.427   -0.246    0.000)    0.493   1.295   (  -4.820   -2.783    0.000)    5.565   1.449   (   0.971    0.561    0.000)    1.122   1.670   (  16.705    9.645    0.000)   19.289   1.923   (   2.203    1.272    0.000)    2.543   2.166   (   1.684    0.972    0.000)    1.944   2.342   (  -8.316   -4.801    0.000)    9.602   3.182   (   2.732    1.577    0.000)    3.154   3.706   (  -0.217   -0.126    0.000)    0.251   3.746   (   0.341    0.197    0.000)    0.394   3.788   (   0.503    0.290    0.000)    0.581   3.847   (   0.214    0.124    0.000)    0.247   3.866   (   0.449    0.259    0.000)    0.519   3.992   (   3.441    1.987    0.000)    3.973   4.746   (  -0.886   -0.512    0.000)    1.023   4.921   (   0.220    0.127    0.000)    0.255   5.442   (  -0.867   -0.501    0.000)    1.001   5.849   (   5.805    3.352    0.000)    6.704   6.073   (  -0.214   -0.124    0.000)    0.248   6.390   (  -1.861   -1.075    0.000)    2.149   6.784   (   0.005    0.003    0.000)    0.006   6.795   (   0.152    0.088    0.000)    0.176   7.339   (  -1.047   -0.604    0.000)    1.209   7.455   (  -0.078   -0.045    0.000)    0.091   9.325   (   0.166    0.096    0.000)    0.191  10.402   (  -2.044   -1.180    0.000)    2.360  11.205   (  -3.277   -1.892    0.000)    3.784======================= Grid point 4 (5/50) =======================q-point: (-0.50  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 65Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.522   (   0.000    0.000    0.000)    0.000   0.816   (  -0.000   -0.000    0.000)    0.000   0.932   (   0.000    0.000    0.000)    0.000   1.064   (   0.000    0.000    0.000)    0.000   1.198   (  -0.000   -0.000    0.000)    0.000   1.457   (   0.000    0.000    0.000)    0.000   1.933   (   0.000    0.000    0.000)    0.000   1.951   (   0.000    0.000    0.000)    0.000   2.180   (   0.000    0.000    0.000)    0.000   2.185   (   0.000    0.000    0.000)    0.000   3.216   (   0.000    0.000    0.000)    0.000   3.704   (   0.000    0.000    0.000)    0.000   3.750   (   0.000    0.000    0.000)    0.000   3.794   (   0.000    0.000    0.000)    0.000   3.850   (   0.000    0.000    0.000)    0.000   3.872   (   0.000    0.000    0.000)    0.000   4.035   (   0.000    0.000    0.000)    0.000   4.735   (   0.000    0.000    0.000)    0.000   4.923   (   0.000    0.000    0.000)    0.000   5.431   (   0.000    0.000    0.000)    0.000   5.924   (   0.000    0.000    0.000)    0.000   6.070   (   0.000    0.000    0.000)    0.000   6.364   (   0.000    0.000    0.000)    0.000   6.784   (   0.000    0.000    0.000)    0.000   6.796   (   0.000    0.000    0.000)    0.000   7.326   (   0.000    0.000    0.000)    0.000   7.454   (   0.000    0.000    0.000)    0.000   9.327   (   0.000    0.000    0.000)    0.000  10.376   (   0.000    0.000    0.000)    0.000  11.165   (   0.000    0.000    0.000)    0.000======================= Grid point 10 (6/50) =======================q-point: ( 0.12  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.282   (   3.867    6.697    0.000)    7.733   0.584   (   8.153   14.121    0.000)   16.305   0.936   (  -0.427   -0.740    0.000)    0.854   0.966   (  -1.397   -2.420    0.000)    2.794   1.034   (  13.238   22.928    0.000)   26.475   1.420   (  -0.738   -1.278    0.000)    1.476   1.431   (   0.084    0.146    0.000)    0.169   1.860   (   1.204    2.085    0.000)    2.407   2.088   (   2.332    4.039    0.000)    4.663   2.529   (  -3.049   -5.281    0.000)    6.098   3.083   (   1.973    3.418    0.000)    3.946   3.716   (  -0.126   -0.217    0.000)    0.251   3.732   (   0.186    0.322    0.000)    0.372   3.772   (   0.630    1.091    0.000)    1.260   3.840   (   0.225    0.389    0.000)    0.450   3.848   (   0.371    0.643    0.000)    0.743   3.858   (   2.465    4.269    0.000)    4.929   4.782   (  -0.494   -0.855    0.000)    0.987   4.909   (   0.178    0.308    0.000)    0.356   5.514   (   0.359    0.621    0.000)    0.717   5.607   (   2.558    4.431    0.000)    5.116   6.081   (  -0.184   -0.319    0.000)    0.368   6.443   (  -0.853   -1.478    0.000)    1.706   6.784   (   0.019    0.032    0.000)    0.037   6.788   (   0.109    0.189    0.000)    0.219   7.378   (  -0.822   -1.423    0.000)    1.644   7.458   (  -0.063   -0.109    0.000)    0.126   9.336   (   1.105    1.913    0.000)    2.209  10.460   (  -1.866   -3.231    0.000)    3.731  11.332   (  -2.920   -5.057    0.000)    5.840======================= Grid point 11 (7/50) =======================q-point: ( 0.25  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 160Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.415   (   4.297    5.043    0.000)    6.626   0.857   (   7.488   12.337    0.000)   14.432   0.898   (  -2.970   -1.497    0.000)    3.326   0.928   (   0.140   -1.809    0.000)    1.815   1.356   (  -4.140   -0.203    0.000)    4.145   1.444   (   1.770    2.820    0.000)    3.329   1.488   (  14.783   10.933    0.000)   18.387   1.915   (   1.652    3.191    0.000)    3.594   2.158   (   1.541    2.747    0.000)    3.149   2.393   (  -5.149   -7.515    0.000)    9.109   3.158   (   2.505    2.948    0.000)    3.868   3.712   (  -0.293    0.077    0.000)    0.303   3.739   (   0.408    0.002    0.000)    0.408   3.796   (   0.292    1.886    0.000)    1.908   3.848   (   0.155    0.447    0.000)    0.473   3.861   (   0.478    0.382    0.000)    0.611   3.946   (   3.409    2.165    0.000)    4.038   4.765   (  -0.913    0.228    0.000)    0.941   4.912   (   0.284   -0.448    0.000)    0.531   5.502   (  -3.305    3.368    0.000)    4.719   5.734   (   7.844    0.533    0.000)    7.863   6.075   (  -0.205   -0.245    0.000)    0.320   6.405   (  -1.470   -1.814    0.000)    2.335   6.785   (  -0.011    0.092    0.000)    0.093   6.792   (   0.159    0.092    0.000)    0.184   7.347   (  -0.982   -1.227    0.000)    1.571   7.456   (  -0.074   -0.091    0.000)    0.117   9.385   (  -0.547    6.452    0.000)    6.475  10.380   (  -1.406   -6.311    0.000)    6.466  11.226   (  -3.056   -4.497    0.000)    5.438======================= Grid point 12 (8/50) =======================q-point: ( 0.38  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 160Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.516   (   2.103    2.690    0.000)    3.415   0.835   (  -2.253   -0.561    0.000)    2.322   0.919   (   0.587   -1.510    0.000)    1.620   1.081   (   5.053   10.288    0.000)   11.462   1.250   (  -6.921   -1.152    0.000)    7.016   1.467   (   0.105    1.002    0.000)    1.008   1.834   (  11.616    5.791    0.000)   12.980   1.967   (   0.761    2.974    0.000)    3.070   2.186   (  -1.721   -1.243    0.000)    2.123   2.214   (  -4.450   -5.175    0.000)    6.825   3.214   (   1.027    1.425    0.000)    1.756   3.708   (  -0.297    0.334    0.000)    0.447   3.744   (   0.376   -0.335    0.000)    0.504   3.814   (  -0.751    2.325    0.000)    2.443   3.853   (  -0.061    0.385    0.000)    0.390   3.870   (   0.247    0.051    0.000)    0.252   4.005   (   2.386   -0.906    0.000)    2.552   4.757   (  -1.154    1.636    0.000)    2.002   4.910   (   0.587   -1.273    0.000)    1.402   5.482   (  -3.074    4.168    0.000)    5.178   5.849   (   5.869   -3.129    0.000)    6.651   6.070   (  -0.074   -0.110    0.000)    0.133   6.365   (  -0.833   -1.262    0.000)    1.513   6.785   (  -0.058    0.132    0.000)    0.144   6.795   (   0.109   -0.038    0.000)    0.115   7.325   (  -0.336   -0.651    0.000)    0.732   7.454   (  -0.026   -0.046    0.000)    0.053   9.432   (  -4.302   10.411    0.000)   11.264  10.307   (   2.219   -8.416    0.000)    8.703  11.154   (  -0.650   -2.773    0.000)    2.848======================= Grid point 19 (9/50) =======================q-point: ( 0.25  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 105Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.510   (   2.311    4.003    0.000)    4.622   0.865   (  -1.027   -1.780    0.000)    2.055   0.894   (  -0.747   -1.293    0.000)    1.493   1.086   (   5.378    9.316    0.000)   10.757   1.349   (  -0.383   -0.663    0.000)    0.766   1.453   (   0.374    0.647    0.000)    0.748   1.741   (   6.524   11.299    0.000)   13.047   1.977   (   1.540    2.667    0.000)    3.079   2.203   (   0.874    1.514    0.000)    1.749   2.220   (  -5.293   -9.169    0.000)   10.587   3.212   (   1.238    2.144    0.000)    2.475   3.715   (   0.138    0.238    0.000)    0.275   3.737   (  -0.105   -0.182    0.000)    0.210   3.843   (   1.439    2.493    0.000)    2.878   3.857   (   0.224    0.388    0.000)    0.447   3.867   (   0.103    0.178    0.000)    0.206   3.968   (   0.007    0.012    0.000)    0.013   4.782   (   0.800    1.386    0.000)    1.601   4.893   (  -0.797   -1.380    0.000)    1.593   5.568   (   1.468    2.543    0.000)    2.936   5.739   (   0.237    0.411    0.000)    0.475   6.071   (  -0.093   -0.162    0.000)    0.187   6.366   (  -1.028   -1.781    0.000)    2.057   6.787   (   0.073    0.126    0.000)    0.146   6.793   (   0.003    0.004    0.000)    0.005   7.325   (  -0.516   -0.894    0.000)    1.032   7.454   (  -0.037   -0.065    0.000)    0.075   9.561   (   5.820   10.081    0.000)   11.641  10.214   (  -5.525   -9.570    0.000)   11.050  11.144   (  -1.861   -3.222    0.000)    3.721======================= Grid point 20 (10/50) =======================q-point: ( 0.38  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 100Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.555   (  -0.600    1.039    0.000)    1.200   0.842   (  -0.662    1.146    0.000)    1.323   0.881   (   1.138   -1.972    0.000)    2.277   1.189   (  -0.140    0.242    0.000)    0.279   1.335   (  -4.845    8.391    0.000)    9.689   1.466   (   1.111   -1.924    0.000)    2.222   1.887   (  -0.063    0.109    0.000)    0.125   2.021   (  -0.689    1.193    0.000)    1.377   2.076   (   1.821   -3.153    0.000)    3.641   2.218   (  -0.349    0.604    0.000)    0.697   3.235   (  -0.320    0.554    0.000)    0.639   3.718   (  -0.344    0.596    0.000)    0.688   3.734   (   0.350   -0.607    0.000)    0.701   3.861   (  -0.213    0.369    0.000)    0.426   3.868   (   0.106   -0.184    0.000)    0.212   3.875   (  -1.944    3.367    0.000)    3.888   3.961   (   1.745   -3.023    0.000)    3.491   4.804   (  -1.537    2.662    0.000)    3.074   4.872   (   1.336   -2.314    0.000)    2.672   5.592   (  -3.585    6.209    0.000)    7.169   5.747   (   3.561   -6.167    0.000)    7.121   6.069   (   0.024   -0.042    0.000)    0.048   6.346   (   0.314   -0.544    0.000)    0.629   6.789   (  -0.106    0.184    0.000)    0.212   6.793   (   0.083   -0.144    0.000)    0.166   7.315   (   0.161   -0.279    0.000)    0.322   7.453   (   0.011   -0.019    0.000)    0.022   9.702   (  -8.380   14.514    0.000)   16.759  10.078   (   7.562  -13.098    0.000)   15.124  11.110   (   0.741   -1.284    0.000)    1.483======================= Grid point 73 (11/50) =======================q-point: ( 0.00  0.00  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 40Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.187   (  -0.000    0.000   10.988)   10.988   0.187   (   0.000    0.000   10.988)   10.988   0.736   (   0.000    0.000   26.261)   26.261   0.962   (  -0.000   -0.000    1.113)    1.113   1.419   (   0.000   -0.000   -1.305)    1.305   1.419   (   0.000    0.000   -1.305)    1.305   1.848   (  -0.000   -0.000   23.371)   23.371   1.913   (   0.000   -0.000    5.996)    5.996   1.959   (   0.000    0.000   -0.020)    0.020   1.959   (   0.000    0.000   -0.020)    0.020   2.973   (  -0.000    0.000   -3.118)    3.118   3.342   (   0.000    0.000  -22.107)   22.107   3.342   (   0.000    0.000  -22.107)   22.107   3.770   (  -0.000   -0.000    4.610)    4.610   3.770   (   0.000    0.000    4.610)    4.610   4.204   (   0.000    0.000   18.548)   18.548   4.204   (   0.000   -0.000   18.548)   18.548   4.826   (   0.000   -0.000    1.930)    1.930   4.826   (   0.000   -0.000    1.930)    1.930   5.389   (   0.000    0.000  -10.300)   10.300   5.389   (   0.000   -0.000  -10.300)   10.300   5.860   (   0.000   -0.000  -18.015)   18.015   6.625   (   0.000   -0.000   10.486)   10.486   6.906   (   0.000    0.000   11.372)   11.372   6.906   (   0.000    0.000   11.372)   11.372   7.393   (   0.000   -0.000   -1.752)    1.752   9.223   (   0.000    0.000   -8.850)    8.850   9.223   (   0.000    0.000   -8.850)    8.850   9.778   (   0.000    0.000    5.489)    5.489  11.376   (   0.000   -0.000   -7.258)    7.258======================= Grid point 74 (12/50) =======================q-point: ( 0.12  0.00  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.382   (  12.171    7.027    5.000)   14.917   0.575   (  15.584    8.997    5.310)   18.762   0.770   (   8.082    4.666   18.741)   20.936   0.956   (  -0.560   -0.323    0.939)    1.140   1.413   (  -0.672   -0.388   -1.122)    1.364   1.424   (   0.066    0.038    1.091)    1.094   1.502   (  -8.812   -5.087   34.718)   36.178   1.879   (  -0.894   -0.516    3.160)    3.324   2.038   (   4.520    2.609    2.574)    5.819   2.296   (   7.806    4.507  -12.619)   15.508   3.001   (   2.297    1.326   -3.194)    4.152   3.349   (   0.575    0.332  -21.947)   21.957   3.370   (   2.247    1.297  -22.009)   22.162   3.768   (  -0.149   -0.086    4.663)    4.666   3.802   (   1.691    0.976    5.568)    5.901   4.206   (   0.178    0.103   18.438)   18.439   4.225   (   1.766    1.020   18.568)   18.680   4.817   (  -0.780   -0.450    1.946)    2.144   4.881   (   3.046    1.759    0.490)    3.551   5.381   (  -0.639   -0.369  -10.111)   10.138   5.428   (   3.419    1.974  -10.955)   11.645   5.857   (  -0.251   -0.145  -18.024)   18.026   6.611   (  -1.121   -0.647   10.237)   10.318   6.905   (  -0.105   -0.061   11.268)   11.269   6.917   (   0.929    0.537   12.048)   12.096   7.381   (  -1.020   -0.589   -1.773)    2.128   8.365   ( -31.030  -17.915   45.632)   58.018   9.224   (   0.105    0.061   -8.913)    8.913  10.158   (   8.236    4.755  -22.483)   24.412  11.331   (  -3.700   -2.136   -7.418)    8.561======================= Grid point 75 (13/50) =======================q-point: ( 0.25  0.00  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.665   (  11.786    6.804    2.205)   13.786   0.708   (   2.842    1.641   11.427)   11.889   0.937   (  -1.020   -0.589    0.641)    1.341   1.190   (  19.966   11.527    2.536)   23.194   1.366   (  -2.506   -1.447   37.920)   38.030   1.378   (  -2.434   -1.405   -0.523)    2.858   1.455   (   0.891    0.515    1.199)    1.580   1.886   (   1.316    0.760    1.561)    2.178   2.139   (   3.875    2.237    3.446)    5.648   2.347   (  -1.816   -1.049  -10.536)   10.743   3.072   (   3.553    2.051   -3.320)    5.278   3.366   (   0.798    0.461  -21.592)   21.612   3.434   (   2.975    1.718  -21.589)   21.860   3.764   (  -0.214   -0.124    4.815)    4.821   3.830   (   0.766    0.442    7.515)    7.567   4.211   (   0.243    0.140   18.152)   18.154   4.279   (   2.721    1.571   18.497)   18.762   4.794   (  -1.134   -0.655    2.027)    2.413   4.932   (   1.399    0.808    1.915)    2.505   5.362   (  -0.958   -0.553   -9.682)    9.745   5.539   (   5.715    3.299  -13.007)   14.585   5.849   (  -0.417   -0.241  -18.068)   18.075   6.578   (  -1.637   -0.945    9.977)   10.154   6.902   (  -0.145   -0.084   11.019)   11.020   6.945   (   1.307    0.754   13.546)   13.630   7.350   (  -1.525   -0.880   -1.821)    2.533   7.940   (  -9.691   -5.595   38.062)   39.673   9.227   (   0.153    0.088   -9.064)    9.066  10.227   (  -0.026   -0.015  -20.077)   20.077  11.222   (  -5.132   -2.963   -7.877)    9.857======================= Grid point 76 (14/50) =======================q-point: ( 0.38  0.00  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.777   (   2.961    1.709    5.638)    6.593   0.908   (  -1.524   -0.880    1.132)    2.093   0.913   (   9.334    5.389    0.764)   10.805   1.298   (  -4.388   -2.534    0.336)    5.078   1.310   (  -1.291   -0.745   43.630)   43.655   1.459   (   0.841    0.486    1.181)    1.529   1.639   (  15.495    8.946   -2.466)   18.061   1.924   (   1.596    0.921    0.133)    1.847   2.208   (   2.017    1.164    4.090)    4.707   2.249   (  -6.223   -3.593   -7.069)   10.079   3.150   (   2.781    1.605   -3.248)    4.567   3.382   (   0.508    0.293  -21.416)   21.424   3.494   (   1.978    1.142  -20.936)   21.060   3.759   (  -0.157   -0.090    5.025)    5.028   3.841   (   0.287    0.166    8.258)    8.265   4.216   (   0.156    0.090   17.820)   17.821   4.339   (   2.138    1.234   18.288)   18.454   4.770   (  -0.820   -0.473    2.184)    2.380   4.952   (   0.409    0.236    2.808)    2.848   5.341   (  -0.723   -0.418   -9.297)    9.334   5.666   (   4.632    2.674  -15.146)   16.063   5.838   (  -0.520   -0.300  -18.189)   18.199   6.542   (  -1.288   -0.744   10.240)   10.348   6.899   (  -0.101   -0.058   10.772)   10.772   6.972   (   0.945    0.545   14.972)   15.012   7.317   (  -1.153   -0.666   -1.879)    2.303   7.806   (  -3.023   -1.746   30.746)   30.944   9.230   (   0.111    0.064   -9.215)    9.216  10.206   (  -1.243   -0.717  -18.319)   18.376  11.117   (  -3.504   -2.023   -8.447)    9.366======================= Grid point 77 (15/50) =======================q-point: (-0.50  0.00  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 84Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.824   (   0.000    0.000    0.735)    0.735   0.878   (   0.000    0.000    3.671)    3.671   1.055   (   0.000    0.000   -0.520)    0.520   1.214   (   0.000    0.000    1.481)    1.481   1.295   (   0.000    0.000   46.003)   46.003   1.466   (   0.000    0.000    0.828)    0.828   1.875   (   0.000    0.000   -4.954)    4.954   1.945   (   0.000    0.000   -0.481)    0.481   2.124   (   0.000    0.000   -4.277)    4.277   2.231   (   0.000    0.000    4.384)    4.384   3.184   (   0.000    0.000   -3.069)    3.069   3.388   (   0.000    0.000  -21.495)   21.495   3.518   (   0.000    0.000  -20.602)   20.602   3.757   (   0.000    0.000    5.140)    5.140   3.844   (   0.000    0.000    8.399)    8.399   4.218   (   0.000    0.000   17.664)   17.664   4.366   (   0.000    0.000   18.171)   18.171   4.760   (   0.000    0.000    2.278)    2.278   4.956   (   0.000    0.000    2.976)    2.976   5.332   (   0.000    0.000   -9.152)    9.152   5.725   (   0.000    0.000  -16.108)   16.108   5.830   (   0.000    0.000  -18.333)   18.333   6.526   (   0.000    0.000   10.563)   10.563   6.897   (   0.000    0.000   10.670)   10.670   6.984   (   0.000    0.000   15.554)   15.554   7.303   (   0.000    0.000   -1.914)    1.914   7.773   (   0.000    0.000   28.600)   28.600   9.232   (   0.000    0.000   -9.278)    9.278  10.189   (   0.000    0.000  -17.740)   17.740  11.075   (   0.000    0.000   -8.720)    8.720======================= Grid point 83 (16/50) =======================q-point: ( 0.12  0.12  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.606   (   7.651   13.251    2.644)   15.528   0.690   (   1.742    3.017   12.470)   12.947   0.943   (  -0.500   -0.866    0.699)    1.220   1.054   (  11.668   20.209    3.471)   23.592   1.379   (  -1.851   -3.206   30.548)   30.771   1.419   (  -0.477   -0.826   -0.146)    0.965   1.430   (  -0.401   -0.694    7.654)    7.696   1.884   (   0.757    1.311    2.077)    2.570   2.116   (   2.569    4.450    3.179)    6.042   2.348   (  -0.305   -0.528  -11.417)   11.433   3.050   (   1.953    3.383   -3.305)    5.117   3.363   (   0.575    0.996  -21.687)   21.717   3.413   (   1.604    2.778  -21.745)   21.980   3.766   (  -0.059   -0.102    4.786)    4.787   3.823   (   0.477    0.827    7.076)    7.140   4.212   (   0.259    0.449   18.251)   18.258   4.260   (   1.325    2.294   18.519)   18.707   4.804   (  -0.460   -0.797    2.012)    2.212   4.920   (   0.943    1.634    1.471)    2.392   5.389   (   0.306    0.529  -10.798)   10.815   5.481   (   2.192    3.797  -11.288)   12.109   5.852   (  -0.227   -0.392  -18.051)   18.057   6.587   (  -0.926   -1.603    9.992)   10.162   6.905   (   0.026    0.046   11.186)   11.186   6.934   (   0.598    1.036   13.034)   13.089   7.359   (  -0.854   -1.479   -1.808)    2.487   8.007   (  -7.731  -13.391   40.747)   43.583   9.242   (   0.991    1.717   -9.133)    9.346  10.215   (   0.058    0.101  -20.652)   20.652  11.253   (  -3.020   -5.231   -7.679)    9.770======================= Grid point 84 (17/50) =======================q-point: ( 0.25  0.12  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 256Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.752   (   2.301    2.675    7.758)    8.523   0.865   (   7.027   11.715    1.126)   13.707   0.922   (  -0.896   -0.678    0.649)    1.297   1.320   (  -1.171   -1.517   41.616)   41.660   1.358   (  -3.879    0.294    0.131)    3.893   1.447   (   1.593    0.001    0.703)    1.741   1.472   (  13.795   11.864   -0.525)   18.203   1.923   (   1.009    2.661    0.733)    2.939   2.196   (   1.864    3.024    3.806)    5.206   2.283   (  -3.291   -5.258   -8.141)   10.235   3.125   (   2.539    2.985   -3.310)    5.130   3.385   (   0.406    1.386  -21.421)   21.470   3.469   (   2.000    1.524  -21.210)   21.358   3.764   (  -0.193    0.240    5.048)    5.058   3.833   (   0.382   -0.154    7.913)    7.924   4.221   (   0.015    0.889   17.961)   17.983   4.309   (   2.064    1.097   18.280)   18.429   4.789   (  -0.851    0.363    2.193)    2.381   4.941   (   0.569   -0.128    2.602)    2.667   5.387   (  -2.231    2.651  -10.339)   10.904   5.581   (   5.787    0.997  -13.112)   14.367   5.842   (  -0.368   -0.469  -18.139)   18.148   6.552   (  -1.164   -1.500   10.048)   10.225   6.905   (  -0.255    0.428   11.168)   11.180   6.956   (   1.050    0.279   14.226)   14.267   7.327   (  -1.066   -1.323   -1.866)    2.523   7.831   (  -3.386   -4.581   32.209)   32.709   9.286   (  -0.545    5.927   -9.581)   11.280  10.180   (  -0.556   -4.947  -18.458)   19.117  11.143   (  -3.296   -4.577   -8.119)    9.886======================= Grid point 85 (18/50) =======================q-point: ( 0.38  0.12  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 256Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.815   (   1.833    1.989    2.435)    3.639   0.889   (  -1.379   -0.855    1.669)    2.327   1.081   (   4.165   10.433    0.198)   11.236   1.257   (  -6.077   -1.326    1.287)    6.352   1.293   (  -0.427   -0.902   45.844)   45.855   1.469   (   0.547    0.361    0.149)    0.672   1.790   (  10.422    5.340   -3.719)   12.287   1.960   (   0.124    2.379   -0.474)    2.429   2.142   (  -4.777   -6.216   -4.735)    9.158   2.240   (   0.223    1.664    4.258)    4.578   3.182   (   1.056    1.480   -3.106)    3.598   3.400   (  -0.392    1.540  -21.417)   21.476   3.505   (   1.203   -0.180  -20.738)   20.773   3.764   (  -0.351    0.621    5.309)    5.357   3.835   (   0.456   -0.745    8.031)    8.078   4.228   (  -0.448    1.148   17.706)   17.748   4.344   (   1.616   -0.864   17.997)   18.090   4.783   (  -1.163    1.842    2.459)    3.286   4.942   (   0.620   -1.219    2.969)    3.269   5.375   (  -2.498    3.560   -9.771)   10.695   5.671   (   4.386   -2.072  -14.813)   15.588   5.831   (  -0.304   -0.360  -18.293)   18.299   6.524   (  -0.415   -0.908   10.466)   10.514   6.905   (  -0.406    0.680   11.151)   11.179   6.972   (   0.855   -0.598   15.018)   15.054   7.303   (  -0.384   -0.709   -1.925)    2.087   7.769   (  -0.587   -1.695   28.209)   28.266   9.329   (  -3.990    9.598   -9.979)   14.409  10.125   (   2.230   -7.502  -17.230)   18.925  11.067   (  -0.837   -2.700   -8.511)    8.968======================= Grid point 92 (19/50) =======================q-point: ( 0.25  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.809   (   1.599    2.769    3.434)    4.692   0.898   (  -0.723   -1.252    0.364)    1.491   1.096   (   5.649    9.785    1.004)   11.343   1.290   (  -0.701   -1.214   46.474)   46.495   1.351   (  -0.773   -1.339    0.476)    1.618   1.440   (   0.290    0.502   -1.082)    1.227   1.709   (   5.874   10.174   -2.654)   12.044   1.974   (   1.182    2.047   -0.173)    2.370   2.152   (  -4.140   -7.171   -5.254)    9.807   2.247   (   1.074    1.861    4.371)    4.871   3.180   (   1.279    2.215   -3.136)    4.046   3.417   (   0.972    1.683  -21.309)   21.398   3.488   (   0.159    0.275  -20.912)   20.914   3.773   (   0.357    0.619    5.551)    5.597   3.824   (  -0.384   -0.665    7.594)    7.633   4.244   (   0.700    1.212   17.758)   17.813   4.317   (  -0.117   -0.203   17.869)   17.871   4.810   (   0.921    1.595    2.578)    3.168   4.925   (  -0.747   -1.293    2.967)    3.322   5.444   (   1.398    2.421  -11.082)   11.429   5.591   (   0.177    0.307  -12.815)   12.820   5.832   (  -0.282   -0.489  -18.268)   18.277   6.523   (  -0.700   -1.212   10.370)   10.464   6.916   (   0.338    0.585   11.822)   11.841   6.956   (  -0.116   -0.201   14.283)   14.285   7.302   (  -0.568   -0.984   -1.934)    2.244   7.765   (  -1.186   -2.054   27.855)   27.956   9.447   (   5.332    9.236  -10.927)   15.269  10.039   (  -4.846   -8.393  -16.490)   19.127  11.059   (  -1.910   -3.308   -8.316)    9.151======================= Grid point 93 (20/50) =======================q-point: ( 0.38  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.848   (  -0.624    1.081    0.594)    1.382   0.875   (   0.289   -0.500    1.017)    1.169   1.216   (  -0.283    0.491    2.472)    2.536   1.277   (   0.161   -0.280   48.097)   48.098   1.318   (  -4.661    8.073   -0.804)    9.356   1.454   (   1.469   -2.545   -1.109)    3.141   1.834   (   0.205   -0.355   -4.714)    4.731   2.003   (  -0.720    1.246   -1.202)    1.875   2.045   (   1.346   -2.332   -2.215)    3.486   2.266   (  -0.446    0.773    4.667)    4.751   3.204   (  -0.332    0.574   -3.019)    3.091   3.439   (  -1.175    2.034  -21.232)   21.362   3.487   (   0.850   -1.472  -20.853)   20.922   3.783   (  -0.642    1.112    5.967)    6.103   3.814   (   0.693   -1.201    7.159)    7.292   4.260   (  -1.059    1.834   17.638)   17.764   4.311   (   1.157   -2.004   17.677)   17.828   4.834   (  -1.620    2.806    2.849)    4.315   4.904   (   1.344   -2.327    3.025)    4.046   5.469   (  -3.005    5.204  -11.044)   12.573   5.596   (   2.776   -4.808  -13.056)   14.188   5.826   (   0.068   -0.117  -18.341)   18.342   6.510   (   0.244   -0.423   10.609)   10.621   6.924   (  -0.619    1.071   12.317)   12.379   6.952   (   0.651   -1.127   14.036)   14.096   7.292   (   0.174   -0.302   -1.970)    2.001   7.745   (   0.363   -0.628   26.422)   26.432   9.577   (  -7.661   13.268  -12.044)   19.489   9.918   (   6.855  -11.873  -15.186)   20.459  11.024   (   0.718   -1.243   -8.421)    8.543======================= Grid point 146 (21/50) =======================q-point: ( 0.00  0.00  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 50Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.458   (   0.000    0.000   15.380)   15.380   0.458   (  -0.000    0.000   15.380)   15.380   0.995   (   0.000    0.000    2.028)    2.028   1.036   (   0.000   -0.000    6.430)    6.430   1.375   (   0.000   -0.000   -3.188)    3.188   1.375   (  -0.000    0.000   -3.188)    3.188   1.970   (   0.000    0.000    1.368)    1.368   1.970   (   0.000    0.000    1.368)    1.368   2.080   (  -0.000    0.000    9.956)    9.956   2.549   (   0.000   -0.000   41.410)   41.410   2.879   (  -0.000   -0.000   -6.111)    6.111   2.899   (   0.000    0.000  -21.139)   21.139   2.899   (  -0.000    0.000  -21.139)   21.139   3.910   (  -0.000    0.000    9.144)    9.144   3.910   (   0.000    0.000    9.144)    9.144   4.536   (   0.000    0.000   13.777)   13.777   4.536   (  -0.000    0.000   13.777)   13.777   4.873   (  -0.000    0.000    2.261)    2.261   4.873   (   0.000    0.000    2.261)    2.261   5.120   (   0.000    0.000  -15.410)   15.410   5.120   (   0.000    0.000  -15.410)   15.410   5.427   (   0.000   -0.000  -24.311)   24.311   6.829   (   0.000   -0.000    9.302)    9.302   7.215   (   0.000   -0.000   18.456)   18.456   7.215   (   0.000    0.000   18.456)   18.456   7.340   (   0.000   -0.000   -3.419)    3.419   8.963   (   0.000    0.000  -16.363)   16.363   8.963   (   0.000    0.000  -16.363)   16.363   9.942   (   0.000    0.000   10.527)   10.527  11.166   (   0.000   -0.000  -13.066)   13.066======================= Grid point 147 (22/50) =======================q-point: ( 0.12  0.00  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 148Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.557   (   7.569    4.370   12.018)   14.860   0.682   (  13.225    7.635    7.560)   17.040   0.982   (  -1.065   -0.615    1.532)    1.965   1.026   (   1.452    0.838    7.085)    7.281   1.374   (  -0.175   -0.101   -2.835)    2.842   1.396   (   1.712    0.988   -2.235)    2.983   1.940   (  -1.384   -0.799    2.825)    3.246   2.007   (   2.709    1.564    3.462)    4.666   2.148   (   4.502    2.599    8.231)    9.735   2.503   (  -3.373   -1.948   38.838)   39.033   2.880   (   0.527    0.304  -14.388)   14.401   2.929   (   2.465    1.423  -21.049)   21.241   2.935   (   2.500    1.443  -12.781)   13.103   3.909   (  -0.067   -0.039    9.171)    9.172   3.951   (   2.769    1.599    9.184)    9.724   4.535   (  -0.041   -0.024   13.625)   13.625   4.557   (   1.737    1.003   13.715)   13.861   4.863   (  -0.786   -0.454    2.139)    2.324   4.913   (   2.954    1.706    1.875)    3.893   5.117   (  -0.215   -0.124  -15.006)   15.008   5.150   (   2.552    1.473  -15.507)   15.784   5.424   (  -0.240   -0.138  -24.230)   24.231   6.813   (  -1.278   -0.738    9.321)    9.437   7.212   (  -0.257   -0.149   18.389)   18.392   7.238   (   1.814    1.047   18.852)   18.968   7.328   (  -1.040   -0.600   -3.413)    3.618   8.857   (  -7.984   -4.609    2.939)    9.676   8.963   (  -0.031   -0.018  -16.456)   16.456   9.875   (  -4.418   -2.551   -1.323)    5.270  11.115   (  -4.172   -2.409  -13.492)   14.326======================= Grid point 148 (23/50) =======================q-point: ( 0.25  0.00  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 150Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.760   (   9.210    5.317    7.457)   12.988   0.827   (   0.787    0.455    3.401)    3.520   0.952   (  -1.400   -0.808    0.690)    1.757   1.248   (  15.270    8.816    3.602)   17.997   1.358   (  -1.423   -0.822   -1.692)    2.359   1.441   (   1.932    1.115   -0.599)    2.310   1.925   (   0.088    0.051    2.474)    2.476   2.069   (   2.189    1.264    5.998)    6.509   2.253   (   4.107    2.371    7.576)    8.938   2.417   (  -3.811   -2.200   36.880)   37.141   2.906   (   1.564    0.903  -18.297)   18.386   3.002   (   3.478    2.008  -20.688)   21.075   3.002   (   3.133    1.809   -8.477)    9.216   3.907   (  -0.068   -0.039    9.246)    9.247   4.009   (   2.047    1.182   10.050)   10.324   4.533   (  -0.144   -0.083   13.166)   13.167   4.609   (   2.573    1.486   13.460)   13.784   4.841   (  -0.971   -0.561    2.148)    2.423   4.981   (   2.504    1.446    2.351)    3.726   5.109   (  -0.460   -0.266  -14.239)   14.249   5.229   (   4.013    2.317  -16.489)   17.128   5.416   (  -0.408   -0.235  -23.990)   23.995   6.777   (  -1.641   -0.947    9.347)    9.537   7.205   (  -0.364   -0.210   18.223)   18.228   7.285   (   0.871    0.503   11.877)   11.919   7.302   (   0.058    0.033    4.557)    4.557   8.662   (  -7.721   -4.458   26.845)   28.287   8.962   (  -0.040   -0.023  -16.682)   16.682   9.782   (  -3.323   -1.919  -18.784)   19.172  10.992   (  -5.894   -3.403  -14.559)   16.072======================= Grid point 149 (24/50) =======================q-point: ( 0.38  0.00  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 148Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.835   (   0.394    0.227    1.843)    1.898   0.921   (  -1.169   -0.675    0.329)    1.390   0.959   (   7.522    4.343    4.195)    9.646   1.304   (  -3.071   -1.773    0.087)    3.547   1.474   (   1.134    0.654    0.836)    1.553   1.584   (  12.021    6.941   -2.070)   14.035   1.936   (   0.622    0.359    1.406)    1.578   2.092   (  -0.446   -0.258    2.749)    2.797   2.327   (   2.180    1.259    7.162)    7.591   2.337   (  -2.694   -1.556   43.073)   43.185   2.943   (   1.411    0.815  -19.710)   19.777   3.067   (   2.212    1.277   -6.328)    6.824   3.074   (   2.453    1.416  -20.138)   20.336   3.906   (  -0.019   -0.011    9.338)    9.338   4.043   (   0.893    0.515   11.078)   11.126   4.529   (  -0.202   -0.117   12.524)   12.526   4.665   (   2.009    1.160   13.127)   13.331   4.823   (  -0.548   -0.317    2.607)    2.682   5.020   (   0.921    0.532    3.352)    3.517   5.098   (  -0.448   -0.259  -13.736)   13.746   5.319   (   3.376    1.949  -18.104)   18.519   5.405   (  -0.612   -0.353  -23.584)   23.595   6.744   (  -1.102   -0.636    9.389)    9.475   7.197   (  -0.257   -0.148   18.050)   18.053   7.260   (  -1.171   -0.676   -2.881)    3.182   7.345   (   1.744    1.007   20.326)   20.426   8.526   (  -3.914   -2.260   36.057)   36.339   8.961   (  -0.026   -0.015  -16.909)   16.909   9.722   (  -1.885   -1.088  -27.014)   27.102  10.869   (  -4.132   -2.385  -15.720)   16.428======================= Grid point 150 (25/50) =======================q-point: (-0.50  0.00  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 89Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.842   (   0.000    0.000    0.966)    0.966   0.905   (   0.000    0.000    0.434)    0.434   1.067   (   0.000    0.000    2.273)    2.273   1.249   (   0.000    0.000    1.634)    1.634   1.486   (   0.000    0.000    1.106)    1.106   1.746   (   0.000    0.000   -5.789)    5.789   1.945   (   0.000    0.000    0.891)    0.891   2.069   (   0.000    0.000    0.206)    0.206   2.306   (   0.000    0.000   48.900)   48.900   2.352   (   0.000    0.000    6.945)    6.945   2.961   (   0.000    0.000  -19.940)   19.940   3.095   (   0.000    0.000   -5.689)    5.689   3.104   (   0.000    0.000  -19.861)   19.861   3.906   (   0.000    0.000    9.382)    9.382   4.053   (   0.000    0.000   11.488)   11.488   4.526   (   0.000    0.000   12.188)   12.188   4.690   (   0.000    0.000   12.990)   12.990   4.816   (   0.000    0.000    2.946)    2.946   5.030   (   0.000    0.000    3.708)    3.708   5.092   (   0.000    0.000  -13.608)   13.608   5.369   (   0.000    0.000  -17.531)   17.531   5.389   (   0.000    0.000  -24.700)   24.700   6.731   (   0.000    0.000    9.421)    9.421   7.194   (   0.000    0.000   17.977)   17.977   7.246   (   0.000    0.000   -3.059)    3.059   7.366   (   0.000    0.000   20.923)   20.923   8.481   (   0.000    0.000   37.996)   37.996   8.961   (   0.000    0.000  -17.004)   17.004   9.699   (   0.000    0.000  -29.051)   29.051  10.819   (   0.000    0.000  -16.238)   16.238======================= Grid point 156 (26/50) =======================q-point: ( 0.12  0.12  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 148Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.715   (   5.852   10.135    8.335)   14.368   0.816   (   1.486    2.574    3.787)    4.814   0.961   (  -0.761   -1.317    0.923)    1.779   1.148   (   7.394   12.807    5.103)   15.644   1.384   (   0.420    0.728   -2.380)    2.524   1.412   (   0.348    0.603   -0.630)    0.940   1.933   (   0.185    0.321    2.789)    2.814   2.048   (   1.223    2.119    5.440)    5.965   2.228   (   2.693    4.664    7.708)    9.403   2.440   (  -2.250   -3.897   37.012)   37.284   2.897   (   0.841    1.457  -17.511)   17.592   2.978   (   1.860    3.222  -20.765)   21.095   2.984   (   1.826    3.162   -9.512)   10.189   3.910   (   0.052    0.090    9.255)    9.255   3.994   (   1.237    2.142    9.720)   10.030   4.536   (   0.062    0.108   13.322)   13.323   4.590   (   1.224    2.120   13.516)   13.736   4.852   (  -0.344   -0.596    2.300)    2.401   4.961   (   1.498    2.595    1.721)    3.456   5.121   (   0.231    0.399  -14.508)   14.516   5.196   (   1.758    3.044  -15.804)   16.191   5.418   (  -0.228   -0.395  -24.086)   24.090   6.787   (  -0.971   -1.682    9.343)    9.542   7.210   (  -0.026   -0.045   18.247)   18.247   7.272   (   1.213    2.102   19.548)   19.698   7.306   (  -0.865   -1.498   -3.380)    3.796   8.705   (  -5.156   -8.930   21.775)   24.093   8.975   (   0.709    1.229  -16.788)   16.848   9.800   (  -2.482   -4.298  -14.584)   15.405  11.027   (  -3.392   -5.874  -14.214)   15.750======================= Grid point 157 (27/50) =======================q-point: ( 0.25  0.12  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 256Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.832   (   0.249    0.340    2.387)    2.423   0.915   (   3.398    5.352    2.955)    6.994   0.942   (   1.364    3.724    2.377)    4.623   1.353   (  -2.624    1.346   -0.715)    3.034   1.438   (   3.721   -1.166   -0.112)    3.902   1.462   (   9.209   10.182    0.409)   13.734   1.948   (  -0.071    1.928    1.991)    2.772   2.079   (   0.929   -0.499    4.097)    4.231   2.314   (   2.149    3.409    7.544)    8.553   2.362   (  -2.436   -2.874   40.324)   40.499   2.935   (   1.267    1.997  -19.486)   19.629   3.044   (   2.201    2.002  -18.338)   18.578   3.050   (   2.201    2.527   -9.035)    9.636   3.913   (  -0.100    0.616    9.414)    9.435   4.026   (   1.049    0.223   10.633)   10.687   4.538   (  -0.297    0.572   12.776)   12.792   4.633   (   1.881    0.586   12.954)   13.103   4.843   (  -0.655    0.963    2.817)    3.048   5.002   (   1.251    0.432    2.418)    2.756   5.126   (  -1.023    1.883  -13.880)   14.044   5.268   (   3.602    1.342  -16.792)   17.227   5.408   (  -0.395   -0.579  -23.891)   23.902   6.753   (  -1.045   -1.395    9.382)    9.543   7.210   (  -0.434    0.751   18.108)   18.129   7.272   (  -1.207   -1.183   -2.750)    3.228   7.316   (   1.954    0.612   19.910)   20.015   8.547   (  -3.747   -6.000   33.311)   34.054   9.008   (  -0.459    4.641  -17.348)   17.964   9.715   (  -1.752   -4.851  -24.074)   24.620  10.903   (  -3.897   -5.025  -15.268)   16.539======================= Grid point 158 (28/50) =======================q-point: ( 0.38  0.12  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 256Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.844   (   0.266    0.740    1.121)    1.369   0.905   (  -0.508   -0.842    0.336)    1.039   1.102   (   2.367    9.270    2.342)    9.850   1.283   (  -4.510   -0.414    1.412)    4.744   1.471   (   1.621   -0.909    0.261)    1.877   1.693   (   6.638    3.203   -4.641)    8.710   1.961   (  -0.648    1.832    1.185)    2.276   2.064   (  -0.481   -2.154    0.690)    2.312   2.311   (  -0.984   -0.881   47.708)   47.726   2.363   (   0.207    2.147    7.395)    7.704   2.970   (   0.161    1.843  -20.113)   20.198   3.085   (   1.250    0.046  -16.415)   16.462   3.098   (   0.906    1.100   -9.390)    9.497   3.917   (  -0.466    1.104    9.550)    9.624   4.037   (   0.910   -1.148   11.220)   11.315   4.537   (  -0.622    0.995   12.252)   12.308   4.660   (   1.705   -1.575   12.271)   12.488   4.848   (  -1.165    2.616    3.617)    4.613   5.014   (   0.811   -1.080    3.144)    3.422   5.125   (  -1.711    2.777  -13.467)   13.857   5.333   (   3.069   -0.769  -17.410)   17.695   5.393   (  -0.548   -0.691  -24.183)   24.199   6.729   (  -0.293   -0.799    9.438)    9.476   7.210   (  -0.810    1.445   18.230)   18.305   7.246   (  -0.481   -0.649   -3.045)    3.150   7.345   (   1.520   -1.001   20.682)   20.762   8.464   (  -0.502   -3.268   36.797)   36.946   9.042   (  -3.196    7.725  -17.612)   19.496   9.653   (   1.258   -5.879  -27.641)   28.287  10.815   (  -1.215   -2.747  -16.074)   16.353======================= Grid point 165 (29/50) =======================q-point: ( 0.25  0.25  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 150Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.846   (   0.565    0.979    1.272)    1.702   0.905   (  -0.682   -1.182    0.237)    1.385   1.132   (   5.400    9.353    2.766)   11.148   1.376   (  -1.074   -1.860    1.981)    2.922   1.406   (   0.387    0.670   -1.877)    2.030   1.640   (   3.813    6.604   -3.463)    8.375   1.984   (   0.780    1.352    1.659)    2.278   2.051   (  -1.083   -1.876    0.908)    2.349   2.316   (  -0.861   -1.491   46.734)   46.766   2.375   (   1.343    2.326    7.868)    8.314   2.981   (   1.399    2.423  -20.216)   20.408   3.069   (   0.227    0.393  -20.005)   20.010   3.096   (   1.044    1.809   -6.145)    6.490   3.932   (   0.681    1.179    9.787)    9.882   4.017   (  -0.538   -0.932   10.887)   10.940   4.554   (   0.537    0.930   12.321)   12.368   4.629   (  -0.421   -0.730   12.083)   12.112   4.878   (   1.275    2.208    3.794)    4.571   4.996   (  -0.521   -0.902    3.024)    3.198   5.172   (   1.259    2.181  -14.148)   14.370   5.284   (   0.244    0.423  -16.319)   16.326   5.394   (  -0.436   -0.755  -24.307)   24.323   6.727   (  -0.588   -1.019    9.454)    9.527   7.231   (   0.638    1.105   15.997)   16.048   7.250   (  -0.461   -0.798   -0.674)    1.142   7.317   (  -0.221   -0.383   20.572)   20.577   8.450   (  -1.961   -3.397   36.139)   36.351   9.133   (   4.164    7.212  -19.371)   21.085   9.593   (  -3.974   -6.882  -25.621)   26.825  10.811   (  -2.102   -3.641  -15.916)   16.462======================= Grid point 166 (30/50) =======================q-point: ( 0.38  0.25  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 148Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.862   (  -0.519    0.899    0.678)    1.239   0.888   (   0.477   -0.826    0.393)    1.032   1.279   (  -0.753    1.305    3.402)    3.720   1.295   (  -4.226    7.319   -0.179)    8.454   1.430   (   2.069   -3.584   -0.871)    4.229   1.708   (   0.474   -0.820   -6.262)    6.333   1.997   (  -0.855    1.481    1.414)    2.219   2.028   (   0.726   -1.257    1.343)    1.977   2.302   (   0.085   -0.146   49.057)   49.057   2.399   (  -0.626    1.084    8.071)    8.168   3.013   (  -1.294    2.242  -20.378)   20.542   3.067   (   0.927   -1.605  -19.914)   20.000   3.115   (  -0.279    0.483   -5.856)    5.882   3.949   (  -1.082    1.873   10.062)   10.292   4.002   (   1.170   -2.026   10.823)   11.073   4.567   (  -1.011    1.751   12.044)   12.212   4.618   (   1.262   -2.185   11.850)   12.116   4.909   (  -1.743    3.018    4.029)    5.327   4.980   (   1.271   -2.201    3.575)    4.387   5.196   (  -2.231    3.864  -14.347)   15.025   5.288   (   1.871   -3.241  -15.771)   16.209   5.385   (   0.088   -0.152  -24.914)   24.914   6.716   (   0.214   -0.371    9.498)    9.508   7.236   (   0.149   -0.258   -3.331)    3.344   7.251   (  -1.314    2.277   19.181)   19.360   7.310   (   1.276   -2.211   20.302)   20.462   8.416   (   0.742   -1.285   36.780)   36.810   9.235   (  -5.859   10.148  -20.750)   23.830   9.495   (   5.178   -8.969  -25.343)   27.377  10.772   (   0.707   -1.225  -16.190)   16.252======================= Grid point 219 (31/50) =======================q-point: ( 0.00  0.00  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 40Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.789   (   0.000    0.000   16.355)   16.355   0.789   (   0.000   -0.000   16.355)   16.355   1.041   (   0.000    0.000    2.337)    2.337   1.096   (   0.000   -0.000    0.583)    0.583   1.278   (   0.000    0.000   -6.682)    6.682   1.278   (   0.000   -0.000   -6.682)    6.682   2.031   (   0.000    0.000    5.086)    5.086   2.031   (   0.000    0.000    5.086)    5.086   2.301   (   0.000    0.000   11.348)   11.348   2.498   (   0.000    0.000  -17.811)   17.811   2.498   (   0.000    0.000  -17.811)   17.811   2.727   (   0.000   -0.000   -8.766)    8.766   3.406   (   0.000    0.000   41.354)   41.354   4.138   (   0.000    0.000   13.162)   13.162   4.138   (   0.000    0.000   13.162)   13.162   4.734   (   0.000   -0.000   -1.527)    1.527   4.734   (   0.000   -0.000   -1.527)    1.527   4.810   (  -0.000   -0.000   -7.278)    7.278   4.810   (   0.000   -0.000   -7.278)    7.278   4.873   (   0.000   -0.000  -30.162)   30.162   4.904   (   0.000    0.000    0.510)    0.510   4.904   (   0.000    0.000    0.510)    0.510   7.001   (   0.000   -0.000    7.642)    7.642   7.255   (   0.000    0.000   -4.914)    4.914   7.636   (  -0.000   -0.000   22.482)   22.482   7.636   (   0.000   -0.000   22.482)   22.482   8.573   (   0.000    0.000  -21.484)   21.484   8.573   (   0.000    0.000  -21.484)   21.484  10.199   (   0.000   -0.000   14.532)   14.532  10.862   (   0.000   -0.000  -16.354)   16.354======================= Grid point 220 (32/50) =======================q-point: ( 0.12  0.00  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.839   (   4.064    2.346   14.685)   15.416   0.873   (   4.121    2.379    9.813)   10.906   1.011   (  -2.125   -1.227    1.114)    2.695   1.121   (   4.235    2.445    3.503)    6.016   1.287   (   0.644    0.372   -6.096)    6.141   1.330   (   3.843    2.219   -4.447)    6.282   2.018   (  -0.793   -0.458    5.594)    5.669   2.058   (   2.057    1.188    4.490)    5.080   2.346   (   3.228    1.864   10.313)   10.966   2.501   (   0.550    0.318  -17.443)   17.455   2.528   (   2.488    1.437  -17.943)   18.171   2.763   (   2.887    1.667   -8.706)    9.322   3.366   (  -2.854   -1.648   42.089)   42.218   4.137   (  -0.087   -0.050   13.106)   13.106   4.176   (   2.704    1.561   12.870)   13.244   4.724   (  -0.854   -0.493   -2.694)    2.869   4.742   (   0.734    0.424   -1.747)    1.941   4.804   (  -0.192   -0.111  -19.279)   19.281   4.810   (   0.024    0.014   -5.175)    5.175   4.905   (   0.014    0.008   -0.286)    0.286   4.911   (   2.783    1.607  -16.992)   17.293   4.950   (   3.423    1.976   -0.493)    3.983   6.986   (  -1.175   -0.678    7.734)    7.852   7.243   (  -1.058   -0.611   -4.914)    5.064   7.632   (  -0.320   -0.185   22.474)   22.477   7.667   (   2.502    1.444   22.936)   23.118   8.571   (  -0.182   -0.105  -21.566)   21.567   8.622   (   4.016    2.319  -20.486)   21.004  10.068   ( -10.503   -6.064   16.033)   20.103  10.797   (  -5.271   -3.043  -17.304)   18.343======================= Grid point 221 (33/50) =======================q-point: ( 0.25  0.00  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.892   (  -1.099   -0.635    3.064)    3.317   0.957   (   5.671    3.274   11.072)   12.864   0.959   (  -2.096   -1.210   -0.072)    2.421   1.299   (   0.228    0.132   -4.369)    4.377   1.325   (  11.040    6.374    3.730)   13.282   1.424   (   3.758    2.169   -1.189)    4.499   2.004   (  -0.359   -0.207    5.856)    5.871   2.109   (   2.038    1.177    4.393)    4.984   2.407   (   1.422    0.821    2.242)    2.779   2.546   (   3.433    1.982  -10.304)   11.040   2.604   (   3.734    2.156  -18.107)   18.613   2.845   (   3.836    2.215   -8.367)    9.467   3.303   (  -2.180   -1.259   43.376)   43.449   4.135   (  -0.098   -0.056   12.962)   12.962   4.240   (   2.529    1.460   12.638)   12.971   4.698   (  -1.307   -0.754   -3.919)    4.200   4.767   (   1.272    0.734   -3.803)    4.076   4.806   (   0.314    0.181  -26.904)   26.907   4.815   (   0.470    0.272   -2.632)    2.688   4.901   (  -0.412   -0.238   -1.256)    1.343   4.988   (   3.611    2.085   -7.450)    8.537   5.034   (   3.333    1.924   -0.367)    3.866   6.953   (  -1.584   -0.915    7.862)    8.072   7.211   (  -1.556   -0.898   -4.891)    5.211   7.623   (  -0.457   -0.264   22.453)   22.459   7.740   (   3.525    2.035   24.109)   24.450   8.566   (  -0.257   -0.148  -21.767)   21.769   8.746   (   6.450    3.724  -17.821)   19.315   9.767   ( -14.346   -8.283   17.509)   24.104  10.642   (  -7.478   -4.317  -19.570)   21.390======================= Grid point 222 (34/50) =======================q-point: ( 0.38  0.00  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.869   (  -0.586   -0.338    1.675)    1.807   0.920   (  -1.263   -0.729   -0.475)    1.534   1.082   (   4.736    2.734    7.320)    9.137   1.291   (  -1.008   -0.582   -1.783)    2.129   1.491   (   1.921    1.109    0.963)    2.418   1.549   (   7.275    4.200   -1.648)    8.560   2.001   (   0.020    0.012    5.538)    5.538   2.144   (   0.965    0.557    6.395)    6.491   2.421   (   0.201    0.116   -4.888)    4.893   2.629   (   3.072    1.774   -3.737)    5.153   2.684   (   2.778    1.604  -18.008)   18.291   2.924   (   2.647    1.528   -7.956)    8.523   3.269   (  -0.813   -0.469   44.196)   44.206   4.133   (  -0.043   -0.025   12.800)   12.800   4.286   (   1.373    0.793   12.822)   12.920   4.669   (  -0.979   -0.565   -3.976)    4.133   4.794   (   0.974    0.562  -15.078)   15.120   4.814   (   0.361    0.209  -19.353)   19.357   4.832   (   0.950    0.549   -2.243)    2.497   4.885   (  -0.932   -0.538   -1.226)    1.631   5.064   (   2.603    1.503   -4.310)    5.254   5.092   (   1.554    0.897    0.211)    1.807   6.920   (  -1.093   -0.631    7.891)    7.992   7.178   (  -1.158   -0.669   -4.832)    5.013   7.613   (  -0.326   -0.188   22.430)   22.433   7.814   (   2.521    1.455   25.390)   25.557   8.561   (  -0.182   -0.105  -21.972)   21.973   8.913   (   7.852    4.533  -16.455)   18.788   9.445   ( -12.816   -7.399   18.590)   23.761  10.486   (  -5.310   -3.066  -21.862)   22.706======================= Grid point 223 (35/50) =======================q-point: (-0.50  0.00  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 84Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.864   (   0.000    0.000    1.171)    1.171   0.904   (   0.000    0.000   -0.456)    0.456   1.148   (   0.000    0.000    5.130)    5.130   1.269   (   0.000    0.000    0.078)    0.078   1.513   (   0.000    0.000    1.659)    1.659   1.636   (   0.000    0.000   -4.670)    4.670   2.002   (   0.000    0.000    5.303)    5.303   2.155   (   0.000    0.000    8.146)    8.146   2.423   (   0.000    0.000   -6.393)    6.393   2.667   (   0.000    0.000   -2.403)    2.403   2.718   (   0.000    0.000  -17.887)   17.887   2.956   (   0.000    0.000   -7.784)    7.784   3.261   (   0.000    0.000   44.434)   44.434   4.133   (   0.000    0.000   12.725)   12.725   4.302   (   0.000    0.000   12.979)   12.979   4.657   (   0.000    0.000   -3.834)    3.834   4.806   (   0.000    0.000  -30.952)   30.952   4.819   (   0.000    0.000   -4.397)    4.397   4.852   (  -0.000   -0.000   -4.963)    4.963   4.865   (   0.000    0.000    1.508)    1.508   5.098   (  -0.000   -0.000   -4.836)    4.836   5.106   (   0.000    0.000    1.628)    1.628   6.907   (   0.000    0.000    7.881)    7.881   7.163   (   0.000    0.000   -4.796)    4.796   7.609   (   0.000    0.000   22.419)   22.419   7.845   (   0.000    0.000   25.961)   25.961   8.558   (   0.000    0.000  -22.058)   22.058   9.063   (   0.000    0.000  -32.381)   32.381   9.236   (  -0.000   -0.000   35.217)   35.217  10.421   (   0.000    0.000  -22.842)   22.842======================= Grid point 229 (36/50) =======================q-point: ( 0.12  0.12  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.896   (  -0.352   -0.609    3.942)    4.004   0.927   (   3.133    5.427   11.051)   12.704   0.977   (  -0.790   -1.368    1.065)    1.905   1.245   (   5.026    8.705    3.772)   10.737   1.315   (   1.620    2.807   -4.222)    5.323   1.389   (   1.828    3.166   -1.721)    4.040   2.016   (   0.140    0.243    5.958)    5.965   2.088   (   0.920    1.593    4.141)    4.531   2.407   (   2.088    3.617    8.401)    9.382   2.510   (   0.585    1.013  -17.999)   18.037   2.587   (   2.291    3.967  -16.329)   16.959   2.822   (   2.194    3.800   -8.487)    9.554   3.318   (  -1.462   -2.533   43.050)   43.149   4.138   (   0.068    0.117   13.037)   13.038   4.221   (   1.415    2.451   12.613)   12.927   4.712   (  -0.446   -0.773   -3.054)    3.182   4.752   (   0.386    0.669   -2.851)    2.954   4.798   (  -0.043   -0.075  -24.054)   24.054   4.821   (   0.462    0.801   -5.417)    5.495   4.903   (   0.113    0.196   -1.562)    1.578   4.969   (   1.875    3.247   -9.801)   10.494   5.008   (   1.973    3.417    0.319)    3.958   6.962   (  -0.921   -1.595    7.841)    8.054   7.220   (  -0.874   -1.514   -4.903)    5.205   7.630   (  -0.016   -0.027   22.499)   22.499   7.714   (   1.719    2.978   23.682)   23.930   8.577   (   0.399    0.691  -21.852)   21.867   8.699   (   2.958    5.124  -18.484)   19.408   9.855   (  -8.082  -13.999   17.325)   23.695  10.687   (  -4.260   -7.378  -18.893)   20.725======================= Grid point 230 (37/50) =======================q-point: ( 0.25  0.12  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 256Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.876   (  -0.632   -0.711    2.022)    2.235   0.930   (  -1.095   -1.426   -0.295)    1.822   1.068   (   3.474    7.030    8.120)   11.288   1.319   (  -0.782    2.021   -2.937)    3.650   1.447   (   3.104    0.096    0.776)    3.200   1.470   (   6.602    6.515    0.095)    9.276   2.022   (  -0.744    1.874    5.938)    6.271   2.119   (   1.735   -0.100    5.350)    5.625   2.456   (  -1.132    3.924   -0.121)    4.086   2.565   (   4.738    0.490  -11.648)   12.584   2.665   (   2.082    2.964  -14.742)   15.180   2.900   (   2.491    2.871   -8.076)    8.925   3.278   (  -0.931   -0.985   43.908)   43.929   4.141   (  -0.164    0.754   12.912)   12.935   4.262   (   1.507    0.512   12.629)   12.729   4.692   (  -1.372    0.153   -3.353)    3.626   4.770   (   1.441   -0.418   -3.305)    3.630   4.801   (   0.337   -0.012  -27.383)   27.385   4.840   (   0.338    1.689   -5.682)    5.938   4.905   (  -0.563    0.728   -1.724)    1.954   5.033   (   2.329    1.942   -5.294)    6.101   5.069   (   1.921    1.393    0.764)    2.493   6.929   (  -1.034   -1.337    7.914)    8.092   7.188   (  -1.074   -1.318   -4.862)    5.151   7.631   (  -0.571    1.065   22.576)   22.608   7.775   (   2.615    1.078   24.647)   24.809   8.595   (  -0.723    3.261  -22.251)   22.500   8.822   (   6.156    3.564  -15.961)   17.474   9.554   ( -10.431  -12.141   17.909)   24.019  10.531   (  -4.998   -6.193  -21.109)   22.559======================= Grid point 231 (38/50) =======================q-point: ( 0.38  0.12  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 256Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.867   (  -0.171    0.162    1.154)    1.178   0.904   (  -0.441   -0.737   -0.394)    0.945   1.180   (   0.882    6.568    4.848)    8.211   1.301   (  -2.755    1.455   -0.077)    3.117   1.486   (   2.196   -1.545    1.366)    3.013   1.601   (   3.611    0.689   -4.169)    5.558   2.026   (  -1.198    2.222    5.767)    6.295   2.141   (   1.137   -0.713    7.685)    7.801   2.464   (  -2.103    3.914   -5.091)    6.757   2.646   (   3.074    0.322   -5.144)    6.001   2.705   (   0.798   -0.097  -16.332)   16.351   2.953   (   1.009    1.128   -7.672)    7.820   3.264   (  -0.301   -0.049   44.156)   44.157   4.146   (  -0.588    1.306   12.774)   12.854   4.282   (   1.234   -1.204   12.827)   12.943   4.674   (  -1.162    0.990   -3.267)    3.606   4.787   (   1.537   -1.603   -2.675)    3.476   4.806   (   0.118    0.083  -30.315)   30.315   4.863   (  -0.054    1.946   -5.470)    5.806   4.897   (  -1.458    1.453   -1.740)    2.695   5.080   (   1.680    0.065   -3.434)    3.823   5.097   (   0.749   -0.356    1.244)    1.495   6.906   (  -0.319   -0.737    7.920)    7.961   7.163   (  -0.394   -0.687   -4.814)    4.879   7.631   (  -1.109    1.971   22.642)   22.755   7.816   (   2.080   -1.365   25.260)   25.382   8.611   (  -2.440    5.169  -22.253)   22.976   8.979   (   8.091    1.328  -21.402)   22.919   9.282   (  -7.810   -7.979   24.645)   27.056  10.420   (  -1.729   -3.188  -22.701)   22.989======================= Grid point 238 (39/50) =======================q-point: ( 0.25  0.25  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.870   (   0.076    0.131    1.138)    1.148   0.905   (  -0.614   -1.064   -0.338)    1.274   1.212   (   3.915    6.781    4.816)    9.193   1.366   (   0.927    1.606   -2.203)    2.880   1.427   (  -0.658   -1.140    2.923)    3.206   1.565   (   1.506    2.608   -3.786)    4.837   2.062   (   0.879    1.523    6.438)    6.674   2.115   (   0.029    0.051    7.032)    7.032   2.538   (   1.660    2.875    5.150)    6.128   2.585   (   1.426    2.470  -18.146)   18.369   2.698   (   0.061    0.106  -13.463)   13.464   2.949   (   1.032    1.788   -7.585)    7.861   3.267   (  -0.133   -0.230   43.899)   43.900   4.164   (   0.781    1.352   12.861)   12.955   4.258   (  -0.477   -0.826   12.671)   12.707   4.699   (   0.237    0.410   -2.436)    2.481   4.757   (  -0.385   -0.667   -2.258)    2.386   4.806   (   0.169    0.293  -30.405)   30.407   4.877   (   1.078    1.867   -6.531)    6.878   4.921   (   0.312    0.540   -3.165)    3.226   5.061   (   0.541    0.937   -2.535)    2.756   5.085   (   0.120    0.207    1.830)    1.846   6.904   (  -0.561   -0.972    7.959)    8.038   7.163   (  -0.560   -0.969   -4.832)    4.960   7.662   (   1.051    1.821   22.961)   23.057   7.779   (  -0.323   -0.559   24.569)   24.577   8.687   (   2.962    5.130  -23.731)   24.459   8.887   (   1.726    2.989  -16.948)   17.296   9.313   (  -6.382  -11.053   21.776)   25.241  10.418   (  -2.563   -4.440  -22.560)   23.135======================= Grid point 239 (40/50) =======================q-point: ( 0.38  0.25  0.38)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 144Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.875   (  -0.298    0.516    0.584)    0.834   0.891   (   0.342   -0.592   -0.076)    0.687   1.316   (  -3.067    5.312    1.852)    6.408   1.339   (  -1.563    2.708    2.037)    3.732   1.432   (   2.237   -3.874    1.364)    4.677   1.585   (   0.799   -1.384   -5.499)    5.727   2.075   (  -1.335    2.312    6.749)    7.258   2.120   (   0.802   -1.389    7.838)    8.000   2.557   (  -2.623    4.544   -0.070)    5.247   2.645   (   0.469   -0.812   -8.273)    8.326   2.678   (   1.068   -1.850  -17.967)   18.094   2.967   (  -0.187    0.324   -7.235)    7.245   3.265   (  -0.060    0.104   43.650)   43.650   4.183   (  -1.255    2.173   12.778)   13.022   4.243   (   1.323   -2.292   12.759)   13.030   4.704   (  -0.999    1.731   -2.173)    2.953   4.749   (   1.108   -1.920   -1.707)    2.798   4.809   (  -0.056    0.097  -31.349)   31.349   4.903   (  -0.967    1.675   -6.779)    7.050   4.923   (  -0.344    0.596   -4.785)    4.834   5.073   (   0.206   -0.357   -0.553)    0.690   5.083   (   0.450   -0.780    2.066)    2.254   6.894   (   0.192   -0.332    7.975)    7.984   7.153   (   0.167   -0.289   -4.818)    4.830   7.686   (  -1.802    3.121   23.262)   23.539   7.768   (   1.739   -3.012   24.323)   24.570   8.750   (  -4.486    7.770  -24.131)   25.745   8.933   (   3.357   -5.815  -28.897)   29.667   9.173   (   1.272   -2.203   34.630)   34.723  10.371   (   0.783   -1.357  -23.149)   23.202======================= Grid point 292 (41/50) =======================q-point: ( 0.00  0.00 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 30Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.083   (   0.000    0.000   -1.604)    1.604   1.083   (   0.000    0.000    1.604)    1.604   1.087   (   0.000    0.000  -12.291)   12.291   1.087   (   0.000    0.000  -12.291)   12.291   1.087   (   0.000   -0.000   12.291)   12.291   1.087   (   0.000   -0.000   12.291)   12.291   2.196   (  -0.000    0.000  -11.515)   11.515   2.196   (   0.000    0.000  -11.515)   11.515   2.196   (   0.000   -0.000   11.515)   11.515   2.196   (   0.000    0.000   11.515)   11.515   2.528   (   0.000   -0.000  -10.706)   10.706   2.528   (   0.000    0.000   10.706)   10.706   4.198   (   0.000    0.000  -36.237)   36.237   4.198   (   0.000    0.000   36.237)   36.237   4.437   (  -0.000   -0.000  -15.935)   15.935   4.437   (   0.000    0.000  -15.935)   15.935   4.437   (   0.000   -0.000   15.935)   15.935   4.437   (   0.000   -0.000   15.935)   15.935   4.881   (   0.000   -0.000   -3.036)    3.036   4.881   (   0.000    0.000   -3.036)    3.036   4.881   (  -0.000   -0.000    3.036)    3.036   4.881   (   0.000    0.000    3.036)    3.036   7.142   (   0.000    0.000   -6.252)    6.252   7.142   (   0.000   -0.000    6.252)    6.252   8.110   (   0.000   -0.000  -23.580)   23.580   8.110   (   0.000    0.000  -23.580)   23.580   8.110   (   0.000   -0.000   23.580)   23.580   8.110   (   0.000    0.000   23.580)   23.580  10.521   (   0.000   -0.000  -16.729)   16.729  10.521   (  -0.000    0.000   16.729)   16.729======================= Grid point 294 (42/50) =======================q-point: ( 0.12  0.00 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 92Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.007   (  -3.017   -1.742   -2.508)    4.292   1.007   (  -3.017   -1.742    2.508)    4.292   1.112   (   1.964    1.134  -11.314)   11.539   1.112   (   1.964    1.134   11.314)   11.539   1.217   (   7.143    4.124   -6.117)   10.268   1.217   (   7.143    4.124    6.117)   10.268   2.191   (  -0.335   -0.194  -11.777)   11.784   2.191   (  -0.335   -0.194   11.777)   11.784   2.221   (   2.079    1.201  -11.756)   11.999   2.221   (   2.079    1.201   11.756)   11.999   2.568   (   3.226    1.863  -10.379)   11.027   2.568   (   3.226    1.863   10.379)   11.027   4.169   (  -2.010   -1.160  -36.456)   36.529   4.169   (  -2.010   -1.160   36.456)   36.529   4.434   (  -0.280   -0.162  -15.743)   15.746   4.434   (  -0.280   -0.162   15.743)   15.746   4.467   (   2.175    1.256  -15.564)   15.765   4.467   (   2.175    1.256   15.564)   15.765   4.879   (  -0.137   -0.079   -2.963)    2.967   4.879   (  -0.137   -0.079    2.963)    2.967   4.923   (   3.302    1.906   -2.935)    4.811   4.923   (   3.302    1.906    2.935)    4.811   7.129   (  -1.077   -0.622   -6.291)    6.413   7.129   (  -1.077   -0.622    6.291)    6.413   8.106   (  -0.289   -0.167  -23.623)   23.626   8.106   (  -0.289   -0.167   23.623)   23.626   8.153   (   3.525    2.035  -24.305)   24.644   8.153   (   3.525    2.035   24.305)   24.644  10.429   (  -7.490   -4.324  -18.505)   20.426  10.429   (  -7.490   -4.324   18.505)   20.426======================= Grid point 296 (43/50) =======================q-point: ( 0.25  0.00 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 93Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.943   (  -2.385   -1.377   -1.698)    3.236   0.943   (  -2.385   -1.377    1.698)    3.236   1.167   (   2.580    1.489   -8.750)    9.243   1.167   (   2.580    1.489    8.750)    9.243   1.389   (   6.993    4.037   -2.392)    8.421   1.389   (   6.993    4.037    2.392)    8.421   2.184   (  -0.284   -0.164  -11.991)   11.996   2.184   (  -0.284   -0.164   11.991)   11.996   2.285   (   3.197    1.846  -12.682)   13.209   2.285   (   3.197    1.846   12.682)   13.209   2.660   (   4.355    2.514   -9.577)   10.817   2.660   (   4.355    2.514    9.577)   10.817   4.127   (  -1.375   -0.794  -37.309)   37.343   4.127   (  -1.375   -0.794   37.309)   37.343   4.425   (  -0.390   -0.225  -15.265)   15.272   4.425   (  -0.390   -0.225   15.265)   15.272   4.521   (   2.224    1.284  -14.688)   14.911   4.521   (   2.224    1.284   14.688)   14.911   4.875   (  -0.234   -0.135   -2.638)    2.652   4.875   (  -0.234   -0.135    2.638)    2.652   5.011   (   3.828    2.210   -2.583)    5.119   5.011   (   3.828    2.210    2.583)    5.119   7.097   (  -1.536   -0.887   -6.337)    6.580   7.097   (  -1.536   -0.887    6.337)    6.580   8.098   (  -0.413   -0.238  -23.730)   23.735   8.098   (  -0.413   -0.238   23.730)   23.735   8.258   (   5.153    2.975  -26.237)   26.903   8.258   (   5.153    2.975   26.237)   26.903  10.209   ( -10.561   -6.098  -22.623)   25.701  10.209   ( -10.561   -6.098   22.623)   25.701======================= Grid point 298 (44/50) =======================q-point: ( 0.38  0.00 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 92Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.901   (  -1.167   -0.674   -1.336)    1.897   0.901   (  -1.167   -0.674    1.336)    1.897   1.219   (   1.697    0.980   -5.387)    5.732   1.219   (   1.697    0.980    5.387)    5.732   1.516   (   3.730    2.154   -1.482)    4.555   1.516   (   3.730    2.154    1.482)    4.555   2.179   (  -0.118   -0.068  -11.979)   11.980   2.179   (  -0.118   -0.068   11.979)   11.980   2.355   (   2.546    1.470  -13.819)   14.128   2.355   (   2.546    1.470   13.819)   14.128   2.750   (   2.996    1.730   -8.863)    9.514   2.750   (   2.996    1.730    8.863)    9.514   4.106   (  -0.478   -0.276  -37.922)   37.926   4.106   (  -0.478   -0.276   37.922)   37.926   4.417   (  -0.274   -0.158  -14.749)   14.752   4.417   (  -0.274   -0.158   14.749)   14.752   4.564   (   1.386    0.800  -14.282)   14.372   4.564   (   1.386    0.800   14.282)   14.372   4.870   (  -0.195   -0.113   -1.942)    1.955   4.870   (  -0.195   -0.113    1.942)    1.955   5.084   (   2.267    1.309   -1.891)    3.230   5.084   (   2.267    1.309    1.891)    3.230   7.064   (  -1.111   -0.641   -6.319)    6.448   7.064   (  -1.111   -0.641    6.319)    6.448   8.089   (  -0.295   -0.170  -23.840)   23.842   8.089   (  -0.295   -0.170   23.840)   23.842   8.369   (   3.950    2.281  -28.822)   29.181   8.369   (   3.950    2.281   28.822)   29.181   9.987   (  -7.651   -4.417  -27.003)   28.412   9.987   (  -7.651   -4.417   27.003)   28.412======================= Grid point 300 (45/50) =======================q-point: (-0.50  0.00 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 55Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.888   (   0.000    0.000   -1.092)    1.092   0.888   (   0.000    0.000    1.092)    1.092   1.240   (   0.000    0.000   -3.161)    3.161   1.240   (   0.000    0.000    3.161)    3.161   1.559   (   0.000    0.000   -2.955)    2.955   1.559   (   0.000    0.000    2.955)    2.955   2.178   (   0.000    0.000  -11.945)   11.945   2.178   (   0.000    0.000   11.945)   11.945   2.387   (   0.000    0.000  -14.270)   14.270   2.387   (   0.000    0.000   14.270)   14.270   2.786   (   0.000    0.000   -8.598)    8.598   2.786   (   0.000    0.000    8.598)    8.598   4.101   (   0.000    0.000  -38.096)   38.096   4.101   (   0.000    0.000   38.096)   38.096   4.414   (   0.000    0.000  -14.517)   14.517   4.414   (   0.000    0.000   14.517)   14.517   4.581   (   0.000    0.000  -14.218)   14.218   4.581   (   0.000    0.000   14.218)   14.218   4.867   (   0.000    0.000   -1.397)    1.397   4.867   (   0.000    0.000    1.397)    1.397   5.111   (   0.000    0.000   -1.413)    1.413   5.111   (   0.000    0.000    1.413)    1.413   7.051   (   0.000    0.000   -6.294)    6.294   7.051   (   0.000    0.000    6.294)    6.294   8.086   (   0.000    0.000  -23.886)   23.886   8.086   (   0.000    0.000   23.886)   23.886   8.419   (   0.000    0.000  -30.309)   30.309   8.419   (   0.000    0.000   30.309)   30.309   9.892   (   0.000    0.000  -29.143)   29.143   9.892   (   0.000    0.000   29.143)   29.143======================= Grid point 312 (46/50) =======================q-point: ( 0.12  0.12 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 92Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.958   (  -1.534   -2.656   -1.787)    3.550   0.958   (  -1.534   -2.656    1.787)    3.550   1.158   (   1.844    3.194   -9.352)   10.053   1.158   (   1.844    3.194    9.352)   10.053   1.338   (   4.013    6.950   -3.284)    8.671   1.338   (   4.013    6.950    3.284)    8.671   2.197   (   0.385    0.667  -12.180)   12.204   2.197   (   0.385    0.667   12.180)   12.204   2.255   (   1.276    2.210  -12.006)   12.275   2.255   (   1.276    2.210   12.006)   12.275   2.633   (   2.481    4.297   -9.779)   10.966   2.633   (   2.481    4.297    9.779)   10.966   4.136   (  -0.948   -1.641  -37.066)   37.114   4.136   (  -0.948   -1.641   37.066)   37.114   4.430   (  -0.102   -0.177  -14.891)   14.893   4.430   (  -0.102   -0.177   14.891)   14.893   4.504   (   1.210    2.096  -14.165)   14.371   4.504   (   1.210    2.096   14.165)   14.371   4.878   (  -0.033   -0.057   -2.794)    2.795   4.878   (  -0.033   -0.057    2.794)    2.795   4.986   (   2.209    3.827   -2.699)    5.178   4.986   (   2.209    3.827    2.699)    5.178   7.106   (  -0.874   -1.514   -6.333)    6.570   7.106   (  -0.874   -1.514    6.333)    6.570   8.104   (  -0.034   -0.058  -22.678)   22.679   8.104   (  -0.034   -0.058   22.678)   22.679   8.223   (   2.692    4.662  -24.384)   24.971   8.223   (   2.692    4.662   24.384)   24.971  10.273   (  -5.984  -10.364  -21.399)   24.518  10.273   (  -5.984  -10.364   21.399)   24.518======================= Grid point 314 (47/50) =======================q-point: ( 0.25  0.12 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 160Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.912   (  -1.186   -1.410   -1.428)    2.331   0.912   (  -1.186   -1.410    1.428)    2.331   1.224   (   1.257    3.888   -6.484)    7.664   1.224   (   1.257    3.888    6.484)    7.664   1.463   (   4.142    2.774   -0.708)    5.035   1.463   (   4.142    2.774    0.708)    5.035   2.208   (  -1.054    2.412  -12.399)   12.675   2.208   (  -1.054    2.412   12.399)   12.675   2.313   (   3.263    1.042  -13.106)   13.546   2.313   (   3.263    1.042   13.106)   13.546   2.721   (   2.808    3.145   -9.120)   10.047   2.721   (   2.808    3.145    9.120)   10.047   4.111   (  -0.560   -0.617  -37.743)   37.753   4.111   (  -0.560   -0.617   37.743)   37.753   4.426   (  -0.394    0.327  -14.036)   14.045   4.426   (  -0.394    0.327   14.036)   14.045   4.544   (   1.481    0.865  -13.163)   13.274   4.544   (   1.481    0.865   13.163)   13.274   4.877   (  -0.256    0.321   -2.305)    2.341   4.877   (  -0.256    0.321    2.305)    2.341   5.058   (   2.323    2.196   -2.163)    3.860   5.058   (   2.323    2.196    2.163)    3.860   7.074   (  -1.039   -1.301   -6.345)    6.560   7.074   (  -1.039   -1.301    6.345)    6.560   8.104   (  -0.464    0.823  -21.898)   21.918   8.104   (  -0.464    0.823   21.898)   21.918   8.324   (   3.792    3.236  -25.213)   25.701   8.324   (   3.792    3.236   25.213)   25.701  10.054   (  -7.129   -8.571  -25.485)   27.816  10.054   (  -7.129   -8.571   25.485)   27.816======================= Grid point 316 (48/50) =======================q-point: ( 0.38  0.12 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 160Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.890   (  -0.387   -0.393   -1.009)    1.150   0.890   (  -0.387   -0.393    1.009)    1.150   1.269   (  -0.749    3.700   -3.204)    4.951   1.269   (  -0.749    3.700    3.204)    4.951   1.529   (   2.406   -0.919   -2.872)    3.858   1.529   (   2.406   -0.919    2.872)    3.858   2.210   (  -1.647    2.917  -12.240)   12.690   2.210   (  -1.647    2.917   12.240)   12.690   2.369   (   2.008   -0.073  -13.950)   14.094   2.369   (   2.008   -0.073   13.950)   14.094   2.779   (   1.213    1.055   -8.872)    9.016   2.779   (   1.213    1.055    8.872)    9.016   4.103   (  -0.158   -0.075  -38.054)   38.055   4.103   (  -0.158   -0.075   38.054)   38.055   4.424   (  -0.528    0.766  -13.400)   13.433   4.424   (  -0.528    0.766   13.400)   13.433   4.568   (   0.973   -0.364  -12.953)   12.995   4.568   (   0.973   -0.364   12.953)   12.995   4.876   (  -0.424    0.677   -1.506)    1.704   4.876   (  -0.424    0.677    1.506)    1.704   5.102   (   1.074    0.383   -1.433)    1.831   5.102   (   1.074    0.383    1.433)    1.831   7.050   (  -0.357   -0.694   -6.322)    6.370   7.050   (  -0.357   -0.694    6.322)    6.370   8.104   (  -0.876    1.547  -21.585)   21.658   8.104   (  -0.876    1.547   21.585)   21.658   8.400   (   2.186    0.601  -26.581)   26.678   8.400   (   2.186    0.601   26.581)   26.678   9.895   (  -2.717   -4.196  -28.754)   29.185   9.895   (  -2.717   -4.196   28.754)   29.185======================= Grid point 330 (49/50) =======================q-point: ( 0.25  0.25 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 93Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.891   (  -0.367   -0.635   -0.926)    1.181   0.891   (  -0.367   -0.635    0.926)    1.181   1.306   (   2.143    3.712   -3.894)    5.791   1.306   (   2.143    3.712    3.894)    5.791   1.492   (   0.134    0.231   -3.396)    3.407   1.492   (   0.134    0.231    3.396)    3.407   2.257   (   1.040    1.802  -11.761)   11.944   2.257   (   1.040    1.802   11.761)   11.944   2.336   (   0.834    1.445  -12.206)   12.319   2.336   (   0.834    1.445   12.206)   12.319   2.773   (   1.044    1.808   -9.139)    9.375   2.773   (   1.044    1.808    9.139)    9.375   4.104   (  -0.102   -0.177  -38.011)   38.012   4.104   (  -0.102   -0.177   38.011)   38.012   4.436   (   0.286    0.495  -10.796)   10.811   4.436   (   0.286    0.495   10.796)   10.811   4.554   (   0.122    0.211  -10.003)   10.006   4.554   (   0.122    0.211   10.003)   10.006   4.885   (   0.245    0.424   -1.586)    1.660   4.885   (   0.245    0.424    1.586)    1.660   5.092   (   0.632    1.095   -1.420)    1.901   5.092   (   0.632    1.095    1.420)    1.901   7.050   (  -0.551   -0.954   -6.350)    6.445   7.050   (  -0.551   -0.954    6.350)    6.445   8.125   (   0.620    1.075  -16.386)   16.433   8.125   (   0.620    1.075   16.386)   16.433   8.378   (   1.146    1.986  -19.407)   19.542   8.378   (   1.146    1.986   19.407)   19.542   9.899   (  -3.520   -6.097  -28.369)   29.229   9.899   (  -3.520   -6.097   28.369)   29.229======================= Grid point 332 (50/50) =======================q-point: ( 0.38  0.25 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 7.40e-04 7.40e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 0.00e+00 0.00e+00 Number of triplets: 92Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.885   (   0.044   -0.076   -0.436)    0.445   0.885   (   0.044   -0.076    0.436)    0.445   1.349   (  -2.010    3.482   -0.910)    4.122   1.349   (  -2.010    3.482    0.910)    4.122   1.489   (   1.401   -2.426   -3.962)    4.852   1.489   (   1.401   -2.426    3.962)    4.852   2.273   (  -1.474    2.552  -10.295)   10.709   2.273   (  -1.474    2.552   10.295)   10.709   2.356   (   0.547   -0.948  -11.011)   11.065   2.356   (   0.547   -0.948   11.011)   11.065   2.791   (  -0.104    0.180   -9.375)    9.378   2.791   (  -0.104    0.180    9.375)    9.378   4.102   (  -0.013    0.023  -38.081)   38.081   4.102   (  -0.013    0.023   38.081)   38.081   4.442   (  -0.419    0.725   -7.537)    7.583   4.442   (  -0.419    0.725    7.537)    7.583   4.555   (   0.358   -0.620   -7.268)    7.303   4.555   (   0.358   -0.620    7.268)    7.303   4.890   (  -0.299    0.518   -0.749)    0.959   4.890   (  -0.299    0.518    0.749)    0.959   5.103   (   0.069   -0.120   -0.688)    0.702   5.103   (   0.069   -0.120    0.688)    0.702   7.040   (   0.174   -0.302   -6.350)    6.360   7.040   (   0.174   -0.302    6.350)    6.360   8.137   (  -0.723    1.252  -11.385)   11.476   8.137   (  -0.723    1.252   11.385)   11.476   8.399   (   0.210   -0.364  -13.736)   13.743   8.399   (   0.210   -0.364   13.736)   13.743   9.834   (   0.961   -1.664  -29.625)   29.688   9.834   (   0.961   -1.664   29.625)   29.688=================== 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/15360   10.0      7.856      7.856      3.676     -0.000     -0.000      0.000 3/15360   20.0      3.499      3.499      4.313     -0.000     -0.000      0.000 3/15360   30.0      2.205      2.205      5.284     -0.000     -0.000      0.000 3/15360   40.0      1.613      1.613      6.229     -0.000     -0.000      0.000 3/15360   50.0      1.278      1.278      6.811     -0.000     -0.000      0.000 3/15360   60.0      1.064      1.064      7.036     -0.000     -0.000      0.000 3/15360   70.0      0.915      0.915      7.020     -0.000     -0.000      0.000 3/15360   80.0      0.805      0.805      6.867     -0.000     -0.000      0.000 3/15360   90.0      0.722      0.722      6.643     -0.000     -0.000      0.000 3/15360  100.0      0.655      0.655      6.388     -0.000     -0.000      0.000 3/15360  110.0      0.601      0.601      6.123     -0.000     -0.000      0.000 3/15360  120.0      0.556      0.556      5.861     -0.000     -0.000      0.000 3/15360  130.0      0.517      0.517      5.608     -0.000     -0.000      0.000 3/15360  140.0      0.484      0.484      5.367     -0.000     -0.000      0.000 3/15360  150.0      0.455      0.455      5.139     -0.000     -0.000      0.000 3/15360  160.0      0.429      0.429      4.925     -0.000     -0.000      0.000 3/15360  170.0      0.407      0.407      4.724     -0.000     -0.000      0.000 3/15360  180.0      0.386      0.386      4.535     -0.000     -0.000      0.000 3/15360  190.0      0.368      0.368      4.359     -0.000     -0.000      0.000 3/15360  200.0      0.351      0.351      4.194     -0.000     -0.000      0.000 3/15360  210.0      0.336      0.336      4.039     -0.000     -0.000      0.000 3/15360  220.0      0.322      0.322      3.895     -0.000     -0.000      0.000 3/15360  230.0      0.309      0.309      3.759     -0.000     -0.000      0.000 3/15360  240.0      0.297      0.297      3.631     -0.000     -0.000      0.000 3/15360  250.0      0.286      0.286      3.512     -0.000     -0.000      0.000 3/15360  260.0      0.275      0.275      3.399     -0.000     -0.000      0.000 3/15360  270.0      0.266      0.266      3.293     -0.000     -0.000      0.000 3/15360  280.0      0.257      0.257      3.192     -0.000     -0.000      0.000 3/15360  290.0      0.249      0.249      3.098     -0.000     -0.000      0.000 3/15360  300.0      0.241      0.241      3.008     -0.000     -0.000      0.000 3/15360  310.0      0.233      0.233      2.924     -0.000     -0.000      0.000 3/15360  320.0      0.227      0.227      2.843     -0.000     -0.000      0.000 3/15360  330.0      0.220      0.220      2.767     -0.000     -0.000      0.000 3/15360  340.0      0.214      0.214      2.695     -0.000     -0.000      0.000 3/15360  350.0      0.208      0.208      2.626     -0.000     -0.000      0.000 3/15360  360.0      0.202      0.202      2.560     -0.000     -0.000      0.000 3/15360  370.0      0.197      0.197      2.498     -0.000     -0.000      0.000 3/15360  380.0      0.192      0.192      2.438     -0.000     -0.000      0.000 3/15360  390.0      0.188      0.188      2.381     -0.000     -0.000      0.000 3/15360  400.0      0.183      0.183      2.327     -0.000     -0.000      0.000 3/15360  410.0      0.179      0.179      2.275     -0.000     -0.000      0.000 3/15360  420.0      0.175      0.175      2.225     -0.000     -0.000      0.000 3/15360  430.0      0.171      0.171      2.177     -0.000     -0.000      0.000 3/15360  440.0      0.167      0.167      2.131     -0.000     -0.000      0.000 3/15360  450.0      0.163      0.163      2.088     -0.000     -0.000      0.000 3/15360  460.0      0.160      0.160      2.045     -0.000     -0.000      0.000 3/15360  470.0      0.157      0.157      2.005     -0.000     -0.000      0.000 3/15360  480.0      0.153      0.153      1.966     -0.000     -0.000      0.000 3/15360  490.0      0.150      0.150      1.928     -0.000     -0.000      0.000 3/15360  500.0      0.147      0.147      1.892     -0.000     -0.000      0.000 3/15360  510.0      0.145      0.145      1.857     -0.000     -0.000      0.000 3/15360  520.0      0.142      0.142      1.823     -0.000     -0.000      0.000 3/15360  530.0      0.139      0.139      1.791     -0.000     -0.000      0.000 3/15360  540.0      0.137      0.137      1.760     -0.000     -0.000      0.000 3/15360  550.0      0.134      0.134      1.729     -0.000     -0.000      0.000 3/15360  560.0      0.132      0.132      1.700     -0.000     -0.000      0.000 3/15360  570.0      0.130      0.130      1.672     -0.000     -0.000      0.000 3/15360  580.0      0.128      0.128      1.644     -0.000     -0.000      0.000 3/15360  590.0      0.125      0.125      1.618     -0.000     -0.000      0.000 3/15360  600.0      0.123      0.123      1.592     -0.000     -0.000      0.000 3/15360  610.0      0.121      0.121      1.567     -0.000     -0.000      0.000 3/15360  620.0      0.119      0.119      1.543     -0.000     -0.000      0.000 3/15360  630.0      0.118      0.118      1.520     -0.000     -0.000      0.000 3/15360  640.0      0.116      0.116      1.497     -0.000     -0.000      0.000 3/15360  650.0      0.114      0.114      1.475     -0.000     -0.000      0.000 3/15360  660.0      0.112      0.112      1.454     -0.000     -0.000      0.000 3/15360  670.0      0.111      0.111      1.433     -0.000     -0.000      0.000 3/15360  680.0      0.109      0.109      1.412     -0.000     -0.000      0.000 3/15360  690.0      0.108      0.108      1.393     -0.000     -0.000      0.000 3/15360  700.0      0.106      0.106      1.374     -0.000     -0.000      0.000 3/15360  710.0      0.105      0.105      1.355     -0.000     -0.000      0.000 3/15360  720.0      0.103      0.103      1.337     -0.000     -0.000      0.000 3/15360  730.0      0.102      0.102      1.319     -0.000     -0.000      0.000 3/15360  740.0      0.100      0.100      1.302     -0.000     -0.000      0.000 3/15360  750.0      0.099      0.099      1.285     -0.000     -0.000      0.000 3/15360  760.0      0.098      0.098      1.269     -0.000     -0.000      0.000 3/15360  770.0      0.097      0.097      1.253     -0.000     -0.000      0.000 3/15360  780.0      0.095      0.095      1.237     -0.000     -0.000      0.000 3/15360  790.0      0.094      0.094      1.222     -0.000     -0.000      0.000 3/15360  800.0      0.093      0.093      1.207     -0.000     -0.000      0.000 3/15360  810.0      0.092      0.092      1.193     -0.000     -0.000      0.000 3/15360  820.0      0.091      0.091      1.179     -0.000     -0.000      0.000 3/15360  830.0      0.090      0.090      1.165     -0.000     -0.000      0.000 3/15360  840.0      0.089      0.089      1.152     -0.000     -0.000      0.000 3/15360  850.0      0.088      0.088      1.138     -0.000     -0.000      0.000 3/15360  860.0      0.087      0.087      1.126     -0.000     -0.000      0.000 3/15360  870.0      0.086      0.086      1.113     -0.000     -0.000      0.000 3/15360  880.0      0.085      0.085      1.101     -0.000     -0.000      0.000 3/15360  890.0      0.084      0.084      1.089     -0.000     -0.000      0.000 3/15360  900.0      0.083      0.083      1.077     -0.000     -0.000      0.000 3/15360  910.0      0.082      0.082      1.065     -0.000     -0.000      0.000 3/15360  920.0      0.081      0.081      1.054     -0.000     -0.000      0.000 3/15360  930.0      0.080      0.080      1.043     -0.000     -0.000      0.000 3/15360  940.0      0.079      0.079      1.032     -0.000     -0.000      0.000 3/15360  950.0      0.078      0.078      1.022     -0.000     -0.000      0.000 3/15360  960.0      0.078      0.078      1.011     -0.000     -0.000      0.000 3/15360  970.0      0.077      0.077      1.001     -0.000     -0.000      0.000 3/15360  980.0      0.076      0.076      0.991     -0.000     -0.000      0.000 3/15360  990.0      0.075      0.075      0.981     -0.000     -0.000      0.000 3/15360 1000.0      0.075      0.075      0.972     -0.000     -0.000      0.000 3/15360Thermal conductivity related properties were written into "kappa-m888.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 05:16:16]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|