# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/7f92c858-33bb-4e73-97f0-573b0e85d485/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 LiBr / Fm-3m (225) / materials id 23259](https://mdr.nims.go.jp/datasets/d05e2793-a256-4fc3-b6ec-1647ea8e41b1)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-08 05:22:33]-------------------------Compiled with OpenMP support (max 128 threads).Running in phonopy.load mode.Python version 3.14.2Spglib version 2.6.1Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".Unit of length: angstromSettings:  Supercell: [2 2 2]  Primitive matrix:    [0.  0.5 0.5]    [0.5 0.  0.5]    [0.5 0.5 0. ]Spacegroup: Fm-3m (225)Number of symmetry operations in supercell: 1536------------------------------ primitive cell ------------------------------Lattice vectors:  a    0.000000000000000    2.703621060000000    2.703621060000000  b    2.703621060000000    0.000000000000000    2.703621060000000  c    2.703621060000000    2.703621060000000    0.000000000000000Atomic positions (fractional):   *1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941   *2 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904-------------------------------- unit cell ---------------------------------Lattice vectors:  a    5.407242120000000    0.000000000000000    0.000000000000000  b    0.000000000000000    5.407242120000000    0.000000000000000  c    0.000000000000000    0.000000000000000    5.407242120000000Atomic positions (fractional):   *1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941 > 1    2 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1    3 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1    4 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 1   *5 Br  0.50000000000000  0.00000000000000  0.00000000000000  79.904 > 2    6 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904 > 2    7 Br  0.00000000000000  0.00000000000000  0.50000000000000  79.904 > 2    8 Br  0.00000000000000  0.50000000000000  0.00000000000000  79.904 > 2-------------------------------- super cell --------------------------------Lattice vectors:  a   10.814484240000001    0.000000000000000    0.000000000000000  b    0.000000000000000   10.814484240000001    0.000000000000000  c    0.000000000000000    0.000000000000000   10.814484240000001Atomic positions (fractional):   *1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941 > 1    2 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 1    4 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 1    5 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1    6 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1    7 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1    8 Li  0.50000000000000  0.50000000000000  0.50000000000000   6.941 > 1    9 Li  0.00000000000000  0.25000000000000  0.25000000000000   6.941 > 1   10 Li  0.50000000000000  0.25000000000000  0.25000000000000   6.941 > 1   11 Li  0.00000000000000  0.75000000000000  0.25000000000000   6.941 > 1   12 Li  0.50000000000000  0.75000000000000  0.25000000000000   6.941 > 1   13 Li  0.00000000000000  0.25000000000000  0.75000000000000   6.941 > 1   14 Li  0.50000000000000  0.25000000000000  0.75000000000000   6.941 > 1   15 Li  0.00000000000000  0.75000000000000  0.75000000000000   6.941 > 1   16 Li  0.50000000000000  0.75000000000000  0.75000000000000   6.941 > 1   17 Li  0.25000000000000  0.00000000000000  0.25000000000000   6.941 > 1   18 Li  0.75000000000000  0.00000000000000  0.25000000000000   6.941 > 1   19 Li  0.25000000000000  0.50000000000000  0.25000000000000   6.941 > 1   20 Li  0.75000000000000  0.50000000000000  0.25000000000000   6.941 > 1   21 Li  0.25000000000000  0.00000000000000  0.75000000000000   6.941 > 1   22 Li  0.75000000000000  0.00000000000000  0.75000000000000   6.941 > 1   23 Li  0.25000000000000  0.50000000000000  0.75000000000000   6.941 > 1   24 Li  0.75000000000000  0.50000000000000  0.75000000000000   6.941 > 1   25 Li  0.25000000000000  0.25000000000000  0.00000000000000   6.941 > 1   26 Li  0.75000000000000  0.25000000000000  0.00000000000000   6.941 > 1   27 Li  0.25000000000000  0.75000000000000  0.00000000000000   6.941 > 1   28 Li  0.75000000000000  0.75000000000000  0.00000000000000   6.941 > 1   29 Li  0.25000000000000  0.25000000000000  0.50000000000000   6.941 > 1   30 Li  0.75000000000000  0.25000000000000  0.50000000000000   6.941 > 1   31 Li  0.25000000000000  0.75000000000000  0.50000000000000   6.941 > 1   32 Li  0.75000000000000  0.75000000000000  0.50000000000000   6.941 > 1  *33 Br  0.25000000000000  0.00000000000000  0.00000000000000  79.904 > 2   34 Br  0.75000000000000  0.00000000000000  0.00000000000000  79.904 > 2   35 Br  0.25000000000000  0.50000000000000  0.00000000000000  79.904 > 2   36 Br  0.75000000000000  0.50000000000000  0.00000000000000  79.904 > 2   37 Br  0.25000000000000  0.00000000000000  0.50000000000000  79.904 > 2   38 Br  0.75000000000000  0.00000000000000  0.50000000000000  79.904 > 2   39 Br  0.25000000000000  0.50000000000000  0.50000000000000  79.904 > 2   40 Br  0.75000000000000  0.50000000000000  0.50000000000000  79.904 > 2   41 Br  0.25000000000000  0.25000000000000  0.25000000000000  79.904 > 2   42 Br  0.75000000000000  0.25000000000000  0.25000000000000  79.904 > 2   43 Br  0.25000000000000  0.75000000000000  0.25000000000000  79.904 > 2   44 Br  0.75000000000000  0.75000000000000  0.25000000000000  79.904 > 2   45 Br  0.25000000000000  0.25000000000000  0.75000000000000  79.904 > 2   46 Br  0.75000000000000  0.25000000000000  0.75000000000000  79.904 > 2   47 Br  0.25000000000000  0.75000000000000  0.75000000000000  79.904 > 2   48 Br  0.75000000000000  0.75000000000000  0.75000000000000  79.904 > 2   49 Br  0.00000000000000  0.00000000000000  0.25000000000000  79.904 > 2   50 Br  0.50000000000000  0.00000000000000  0.25000000000000  79.904 > 2   51 Br  0.00000000000000  0.50000000000000  0.25000000000000  79.904 > 2   52 Br  0.50000000000000  0.50000000000000  0.25000000000000  79.904 > 2   53 Br  0.00000000000000  0.00000000000000  0.75000000000000  79.904 > 2   54 Br  0.50000000000000  0.00000000000000  0.75000000000000  79.904 > 2   55 Br  0.00000000000000  0.50000000000000  0.75000000000000  79.904 > 2   56 Br  0.50000000000000  0.50000000000000  0.75000000000000  79.904 > 2   57 Br  0.00000000000000  0.25000000000000  0.00000000000000  79.904 > 2   58 Br  0.50000000000000  0.25000000000000  0.00000000000000  79.904 > 2   59 Br  0.00000000000000  0.75000000000000  0.00000000000000  79.904 > 2   60 Br  0.50000000000000  0.75000000000000  0.00000000000000  79.904 > 2   61 Br  0.00000000000000  0.25000000000000  0.50000000000000  79.904 > 2   62 Br  0.50000000000000  0.25000000000000  0.50000000000000  79.904 > 2   63 Br  0.00000000000000  0.75000000000000  0.50000000000000  79.904 > 2   64 Br  0.50000000000000  0.75000000000000  0.50000000000000  79.904 > 2----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            3.5739202    0.0000000    0.0000000            0.0000000    3.5739202    0.0000000            0.0000000    0.0000000    3.5739202-------------------------- Born effective charges --------------------------    1 Li    1.2211001    0.0000000    0.0000000            0.0000000    1.2211001    0.0000000            0.0000000    0.0000000    1.2211001    2 Br   -1.2211001    0.0000000    0.0000000            0.0000000   -1.2211001    0.0000000            0.0000000    0.0000000   -1.2211001----------------------------------------------------------------------------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: 192/192Permutation basis: 852/852Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 33Number of blocks in projector: 33Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 28Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 3Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 2Use standard eigh solver.Tree of FC basis block matrices:- (33, 31), data: False|-- (2, 2), data: True|-- (3, 3), data: True|-- (28, 26), data: True-----Solver_atoms: 1 -- 64 / 64Time (Solver_compr_matrix_reshape): 0.001Solver_block: 80 / 80 - Time: 0.006Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.007--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 192/192Permutation basis: 852/852Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 33Number of blocks in projector: 33Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 28Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 3Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 2Use standard eigh solver.Tree of FC basis block matrices:- (33, 31), data: False|-- (2, 2), data: True|-- (3, 3), data: True|-- (28, 26), data: TrueMax drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-08 05:22:36]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 05:22:36]-------------------------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 2]Primitive matrix:  [0.  0.5 0.5]  [0.5 0.  0.5]  [0.5 0.5 0. ]Spacegroup: Fm-3m (225)------------------------------ primitive cell ------------------------------Lattice vectors:  a    0.000000000000000    2.703621060000000    2.703621060000000  b    2.703621060000000    0.000000000000000    2.703621060000000  c    2.703621060000000    2.703621060000000    0.000000000000000Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941    2 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904-------------------------------- supercell ---------------------------------Lattice vectors:  a   10.814484240000001    0.000000000000000    0.000000000000000  b    0.000000000000000   10.814484240000001    0.000000000000000  c    0.000000000000000    0.000000000000000   10.814484240000001Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941 > 1    2 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 1    4 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 1    5 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1    6 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1    7 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1    8 Li  0.50000000000000  0.50000000000000  0.50000000000000   6.941 > 1    9 Li  0.00000000000000  0.25000000000000  0.25000000000000   6.941 > 1   10 Li  0.50000000000000  0.25000000000000  0.25000000000000   6.941 > 1   11 Li  0.00000000000000  0.75000000000000  0.25000000000000   6.941 > 1   12 Li  0.50000000000000  0.75000000000000  0.25000000000000   6.941 > 1   13 Li  0.00000000000000  0.25000000000000  0.75000000000000   6.941 > 1   14 Li  0.50000000000000  0.25000000000000  0.75000000000000   6.941 > 1   15 Li  0.00000000000000  0.75000000000000  0.75000000000000   6.941 > 1   16 Li  0.50000000000000  0.75000000000000  0.75000000000000   6.941 > 1   17 Li  0.25000000000000  0.00000000000000  0.25000000000000   6.941 > 1   18 Li  0.75000000000000  0.00000000000000  0.25000000000000   6.941 > 1   19 Li  0.25000000000000  0.50000000000000  0.25000000000000   6.941 > 1   20 Li  0.75000000000000  0.50000000000000  0.25000000000000   6.941 > 1   21 Li  0.25000000000000  0.00000000000000  0.75000000000000   6.941 > 1   22 Li  0.75000000000000  0.00000000000000  0.75000000000000   6.941 > 1   23 Li  0.25000000000000  0.50000000000000  0.75000000000000   6.941 > 1   24 Li  0.75000000000000  0.50000000000000  0.75000000000000   6.941 > 1   25 Li  0.25000000000000  0.25000000000000  0.00000000000000   6.941 > 1   26 Li  0.75000000000000  0.25000000000000  0.00000000000000   6.941 > 1   27 Li  0.25000000000000  0.75000000000000  0.00000000000000   6.941 > 1   28 Li  0.75000000000000  0.75000000000000  0.00000000000000   6.941 > 1   29 Li  0.25000000000000  0.25000000000000  0.50000000000000   6.941 > 1   30 Li  0.75000000000000  0.25000000000000  0.50000000000000   6.941 > 1   31 Li  0.25000000000000  0.75000000000000  0.50000000000000   6.941 > 1   32 Li  0.75000000000000  0.75000000000000  0.50000000000000   6.941 > 1   33 Br  0.25000000000000  0.00000000000000  0.00000000000000  79.904 > 33   34 Br  0.75000000000000  0.00000000000000  0.00000000000000  79.904 > 33   35 Br  0.25000000000000  0.50000000000000  0.00000000000000  79.904 > 33   36 Br  0.75000000000000  0.50000000000000  0.00000000000000  79.904 > 33   37 Br  0.25000000000000  0.00000000000000  0.50000000000000  79.904 > 33   38 Br  0.75000000000000  0.00000000000000  0.50000000000000  79.904 > 33   39 Br  0.25000000000000  0.50000000000000  0.50000000000000  79.904 > 33   40 Br  0.75000000000000  0.50000000000000  0.50000000000000  79.904 > 33   41 Br  0.25000000000000  0.25000000000000  0.25000000000000  79.904 > 33   42 Br  0.75000000000000  0.25000000000000  0.25000000000000  79.904 > 33   43 Br  0.25000000000000  0.75000000000000  0.25000000000000  79.904 > 33   44 Br  0.75000000000000  0.75000000000000  0.25000000000000  79.904 > 33   45 Br  0.25000000000000  0.25000000000000  0.75000000000000  79.904 > 33   46 Br  0.75000000000000  0.25000000000000  0.75000000000000  79.904 > 33   47 Br  0.25000000000000  0.75000000000000  0.75000000000000  79.904 > 33   48 Br  0.75000000000000  0.75000000000000  0.75000000000000  79.904 > 33   49 Br  0.00000000000000  0.00000000000000  0.25000000000000  79.904 > 33   50 Br  0.50000000000000  0.00000000000000  0.25000000000000  79.904 > 33   51 Br  0.00000000000000  0.50000000000000  0.25000000000000  79.904 > 33   52 Br  0.50000000000000  0.50000000000000  0.25000000000000  79.904 > 33   53 Br  0.00000000000000  0.00000000000000  0.75000000000000  79.904 > 33   54 Br  0.50000000000000  0.00000000000000  0.75000000000000  79.904 > 33   55 Br  0.00000000000000  0.50000000000000  0.75000000000000  79.904 > 33   56 Br  0.50000000000000  0.50000000000000  0.75000000000000  79.904 > 33   57 Br  0.00000000000000  0.25000000000000  0.00000000000000  79.904 > 33   58 Br  0.50000000000000  0.25000000000000  0.00000000000000  79.904 > 33   59 Br  0.00000000000000  0.75000000000000  0.00000000000000  79.904 > 33   60 Br  0.50000000000000  0.75000000000000  0.00000000000000  79.904 > 33   61 Br  0.00000000000000  0.25000000000000  0.50000000000000  79.904 > 33   62 Br  0.50000000000000  0.25000000000000  0.50000000000000  79.904 > 33   63 Br  0.00000000000000  0.75000000000000  0.50000000000000  79.904 > 33   64 Br  0.50000000000000  0.75000000000000  0.50000000000000  79.904 > 33----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            3.5739202    0.0000000    0.0000000            0.0000000    3.5739202    0.0000000            0.0000000    0.0000000    3.5739202-------------------------- Born effective charges --------------------------    1 Li    1.2211001    0.0000000    0.0000000            0.0000000    1.2211001    0.0000000            0.0000000    0.0000000    1.2211001    2 Br   -1.2211001    0.0000000    0.0000000            0.0000000   -1.2211001    0.0000000            0.0000000    0.0000000   -1.2211001----------------------------------------------------------------------------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]Computing fc3[ 33, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]Expanding fc3.Symmetrizing fc3 by traditional approach (N=3).Symmetrizing fc2 by traditional approach (N=3).Max drift of fc3: 0.00000003 (yyy) 0.00000003 (yyy) 0.00000003 (yyy)fc3 was written into "fc3.hdf5".Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330----------- None of ph-ph interaction calculation was performed. -----------Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 05:22:37]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 05:22:38]-------------------------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 2]Primitive matrix:  [0.  0.5 0.5]  [0.5 0.  0.5]  [0.5 0.5 0. ]Spacegroup: Fm-3m (225)------------------------------ primitive cell ------------------------------Lattice vectors:  a    0.000000000000000    2.703621060000000    2.703621060000000  b    2.703621060000000    0.000000000000000    2.703621060000000  c    2.703621060000000    2.703621060000000    0.000000000000000Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941    2 Br  0.50000000000000  0.50000000000000  0.50000000000000  79.904-------------------------------- supercell ---------------------------------Lattice vectors:  a   10.814484240000001    0.000000000000000    0.000000000000000  b    0.000000000000000   10.814484240000001    0.000000000000000  c    0.000000000000000    0.000000000000000   10.814484240000001Atomic positions (fractional):    1 Li  0.00000000000000  0.00000000000000  0.00000000000000   6.941 > 1    2 Li  0.50000000000000  0.00000000000000  0.00000000000000   6.941 > 1    3 Li  0.00000000000000  0.50000000000000  0.00000000000000   6.941 > 1    4 Li  0.50000000000000  0.50000000000000  0.00000000000000   6.941 > 1    5 Li  0.00000000000000  0.00000000000000  0.50000000000000   6.941 > 1    6 Li  0.50000000000000  0.00000000000000  0.50000000000000   6.941 > 1    7 Li  0.00000000000000  0.50000000000000  0.50000000000000   6.941 > 1    8 Li  0.50000000000000  0.50000000000000  0.50000000000000   6.941 > 1    9 Li  0.00000000000000  0.25000000000000  0.25000000000000   6.941 > 1   10 Li  0.50000000000000  0.25000000000000  0.25000000000000   6.941 > 1   11 Li  0.00000000000000  0.75000000000000  0.25000000000000   6.941 > 1   12 Li  0.50000000000000  0.75000000000000  0.25000000000000   6.941 > 1   13 Li  0.00000000000000  0.25000000000000  0.75000000000000   6.941 > 1   14 Li  0.50000000000000  0.25000000000000  0.75000000000000   6.941 > 1   15 Li  0.00000000000000  0.75000000000000  0.75000000000000   6.941 > 1   16 Li  0.50000000000000  0.75000000000000  0.75000000000000   6.941 > 1   17 Li  0.25000000000000  0.00000000000000  0.25000000000000   6.941 > 1   18 Li  0.75000000000000  0.00000000000000  0.25000000000000   6.941 > 1   19 Li  0.25000000000000  0.50000000000000  0.25000000000000   6.941 > 1   20 Li  0.75000000000000  0.50000000000000  0.25000000000000   6.941 > 1   21 Li  0.25000000000000  0.00000000000000  0.75000000000000   6.941 > 1   22 Li  0.75000000000000  0.00000000000000  0.75000000000000   6.941 > 1   23 Li  0.25000000000000  0.50000000000000  0.75000000000000   6.941 > 1   24 Li  0.75000000000000  0.50000000000000  0.75000000000000   6.941 > 1   25 Li  0.25000000000000  0.25000000000000  0.00000000000000   6.941 > 1   26 Li  0.75000000000000  0.25000000000000  0.00000000000000   6.941 > 1   27 Li  0.25000000000000  0.75000000000000  0.00000000000000   6.941 > 1   28 Li  0.75000000000000  0.75000000000000  0.00000000000000   6.941 > 1   29 Li  0.25000000000000  0.25000000000000  0.50000000000000   6.941 > 1   30 Li  0.75000000000000  0.25000000000000  0.50000000000000   6.941 > 1   31 Li  0.25000000000000  0.75000000000000  0.50000000000000   6.941 > 1   32 Li  0.75000000000000  0.75000000000000  0.50000000000000   6.941 > 1   33 Br  0.25000000000000  0.00000000000000  0.00000000000000  79.904 > 33   34 Br  0.75000000000000  0.00000000000000  0.00000000000000  79.904 > 33   35 Br  0.25000000000000  0.50000000000000  0.00000000000000  79.904 > 33   36 Br  0.75000000000000  0.50000000000000  0.00000000000000  79.904 > 33   37 Br  0.25000000000000  0.00000000000000  0.50000000000000  79.904 > 33   38 Br  0.75000000000000  0.00000000000000  0.50000000000000  79.904 > 33   39 Br  0.25000000000000  0.50000000000000  0.50000000000000  79.904 > 33   40 Br  0.75000000000000  0.50000000000000  0.50000000000000  79.904 > 33   41 Br  0.25000000000000  0.25000000000000  0.25000000000000  79.904 > 33   42 Br  0.75000000000000  0.25000000000000  0.25000000000000  79.904 > 33   43 Br  0.25000000000000  0.75000000000000  0.25000000000000  79.904 > 33   44 Br  0.75000000000000  0.75000000000000  0.25000000000000  79.904 > 33   45 Br  0.25000000000000  0.25000000000000  0.75000000000000  79.904 > 33   46 Br  0.75000000000000  0.25000000000000  0.75000000000000  79.904 > 33   47 Br  0.25000000000000  0.75000000000000  0.75000000000000  79.904 > 33   48 Br  0.75000000000000  0.75000000000000  0.75000000000000  79.904 > 33   49 Br  0.00000000000000  0.00000000000000  0.25000000000000  79.904 > 33   50 Br  0.50000000000000  0.00000000000000  0.25000000000000  79.904 > 33   51 Br  0.00000000000000  0.50000000000000  0.25000000000000  79.904 > 33   52 Br  0.50000000000000  0.50000000000000  0.25000000000000  79.904 > 33   53 Br  0.00000000000000  0.00000000000000  0.75000000000000  79.904 > 33   54 Br  0.50000000000000  0.00000000000000  0.75000000000000  79.904 > 33   55 Br  0.00000000000000  0.50000000000000  0.75000000000000  79.904 > 33   56 Br  0.50000000000000  0.50000000000000  0.75000000000000  79.904 > 33   57 Br  0.00000000000000  0.25000000000000  0.00000000000000  79.904 > 33   58 Br  0.50000000000000  0.25000000000000  0.00000000000000  79.904 > 33   59 Br  0.00000000000000  0.75000000000000  0.00000000000000  79.904 > 33   60 Br  0.50000000000000  0.75000000000000  0.00000000000000  79.904 > 33   61 Br  0.00000000000000  0.25000000000000  0.50000000000000  79.904 > 33   62 Br  0.50000000000000  0.25000000000000  0.50000000000000  79.904 > 33   63 Br  0.00000000000000  0.75000000000000  0.50000000000000  79.904 > 33   64 Br  0.50000000000000  0.75000000000000  0.50000000000000  79.904 > 33----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            3.5739202    0.0000000    0.0000000            0.0000000    3.5739202    0.0000000            0.0000000    0.0000000    3.5739202-------------------------- Born effective charges --------------------------    1 Li    1.2211001    0.0000000    0.0000000            0.0000000    1.2211001    0.0000000            0.0000000    0.0000000    1.2211001    2 Br   -1.2211001    0.0000000    0.0000000            0.0000000   -1.2211001    0.0000000            0.0000000    0.0000000   -1.2211001----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: 0.00000003 (yyy) 0.00000003 (yyy) 0.00000003 (yyy)Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330Length for sampling mesh generation: 50.00Generating grid system ... [ 16 16 16 ]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: 1.22, Number of G-points: 307, Lambda: 0.24Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/145) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 145Calculating 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   5.472   (   0.000    0.000    0.000)    0.000   5.472   (   0.000    0.000    0.000)    0.000   5.472   (   0.000    0.000    0.000)    0.000======================= Grid point 1 (2/145) =======================q-point: ( 0.06  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 408Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.444   ( -12.600   12.600   12.600)   21.823   0.444   ( -12.600   12.600   12.600)   21.823   0.780   ( -22.196   22.196   22.196)   38.445   5.441   (   1.623   -1.623   -1.623)    2.811   5.441   (   1.623   -1.623   -1.623)    2.811  10.116   (   1.852   -1.852   -1.852)    3.209======================= Grid point 2 (3/145) =======================q-point: ( 0.12  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 417Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.867   ( -11.748   11.748   11.748)   20.348   0.867   ( -11.748   11.748   11.748)   20.348   1.528   ( -20.816   20.816   20.816)   36.054   5.378   (   1.620   -1.620   -1.620)    2.806   5.378   (   1.620   -1.620   -1.620)    2.806  10.044   (   1.906   -1.906   -1.906)    3.302======================= Grid point 3 (4/145) =======================q-point: ( 0.19  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 408Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.252   ( -10.334   10.334   10.334)   17.900   1.252   ( -10.334   10.334   10.334)   17.900   2.211   ( -18.343   18.343   18.343)   31.772   5.354   (  -0.541    0.541    0.541)    0.938   5.354   (  -0.541    0.541    0.541)    0.938  10.013   (  -0.526    0.526    0.526)    0.911======================= Grid point 4 (5/145) =======================q-point: ( 0.25  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 417Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.578   (  -8.414    8.414    8.414)   14.573   1.578   (  -8.414    8.414    8.414)   14.573   2.788   ( -14.821   14.821   14.821)   25.671   5.429   (  -3.771    3.771    3.771)    6.532   5.429   (  -3.771    3.771    3.771)    6.532  10.097   (  -4.406    4.406    4.406)    7.631======================= Grid point 5 (6/145) =======================q-point: ( 0.31  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 408Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.832   (  -6.222    6.222    6.222)   10.777   1.832   (  -6.222    6.222    6.222)   10.777   3.232   ( -10.738   10.738   10.738)   18.598   5.605   (  -6.063    6.063    6.063)   10.502   5.605   (  -6.063    6.063    6.063)   10.502  10.309   (  -7.459    7.459    7.459)   12.920======================= Grid point 6 (7/145) =======================q-point: ( 0.38  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 417Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.010   (  -4.046    4.046    4.046)    7.008   2.010   (  -4.046    4.046    4.046)    7.008   3.534   (  -6.761    6.761    6.761)   11.711   5.822   (  -6.035    6.035    6.035)   10.452   5.822   (  -6.035    6.035    6.035)   10.452  10.582   (  -7.731    7.731    7.731)   13.390======================= Grid point 7 (8/145) =======================q-point: ( 0.44  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 408Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.114   (  -1.987    1.987    1.987)    3.442   2.114   (  -1.987    1.987    1.987)    3.442   3.706   (  -3.219    3.219    3.219)    5.576   5.996   (  -3.687    3.687    3.687)    6.386   5.996   (  -3.687    3.687    3.687)    6.386  10.808   (  -4.850    4.850    4.850)    8.401======================= Grid point 8 (9/145) =======================q-point: (-0.50  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 213Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.148   (  -0.000    0.000    0.000)    0.000   2.148   (  -0.000    0.000    0.000)    0.000   3.761   (  -0.000    0.000    0.000)    0.000   6.062   (  -0.000    0.000    0.000)    0.000   6.062   (  -0.000    0.000    0.000)    0.000  10.895   (  -0.000    0.000    0.000)    0.000======================= Grid point 18 (10/145) =======================q-point: ( 0.06  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 344Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.523   (   0.000   -0.000   22.052)   22.052   0.523   (   0.000   -0.000   22.052)   22.052   0.881   (   0.000   -0.000   37.257)   37.257   5.491   (   0.000   -0.000    1.639)    1.639   5.491   (   0.000   -0.000    1.639)    1.639  10.030   (  -0.000    0.000  -10.368)   10.368======================= Grid point 19 (11/145) =======================q-point: ( 0.12  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.844   (  -5.876    5.876   18.845)   20.596   0.847   (  -5.970    5.970   19.520)   21.267   1.444   ( -11.537   11.537   31.340)   35.333   5.429   (   5.492   -5.492    2.269)    8.091   5.485   (   2.399   -2.399    2.648)    4.304   9.930   (  -3.193    3.193  -13.739)   14.462======================= Grid point 20 (12/145) =======================q-point: ( 0.19  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.197   (  -7.255    7.255   15.267)   18.395   1.213   (  -7.031    7.031   16.930)   19.634   2.073   ( -14.181   14.181   24.636)   31.767   5.351   (   3.032   -3.032    2.936)    5.197   5.462   (   2.570   -2.570    4.328)    5.652   9.871   (  -5.565    5.565  -12.977)   15.177======================= Grid point 21 (13/145) =======================q-point: ( 0.25  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.521   (  -6.906    6.906   12.192)   15.622   1.548   (  -5.961    5.961   14.543)   16.809   2.646   ( -13.492   13.492   18.401)   26.508   5.379   (  -1.596    1.596    5.066)    5.546   5.485   (   0.329   -0.329    6.132)    6.149   9.908   (  -9.333    9.333   -9.934)   16.519======================= Grid point 22 (14/145) =======================q-point: ( 0.31  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.792   (  -5.533    5.533    9.526)   12.327   1.819   (  -4.086    4.086   12.166)   13.468   3.110   ( -11.055   11.055   12.536)   20.039   5.533   (  -5.240    5.240    6.983)   10.182   5.592   (  -2.858    2.858    7.043)    8.121  10.073   ( -13.145   13.145   -7.095)   19.897======================= Grid point 23 (15/145) =======================q-point: ( 0.38  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.994   (  -3.552    3.552    7.154)    8.742   2.018   (  -2.068    2.068   10.075)   10.491   3.444   (  -7.910    7.910    7.390)   13.407   5.753   (  -6.258    6.258    6.600)   11.040   5.765   (  -4.898    4.898    6.421)    9.445  10.326   ( -14.857   14.857   -6.593)   22.021======================= Grid point 24 (16/145) =======================q-point: ( 0.44  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.120   (  -1.397    1.397    4.890)    5.274   2.149   (  -0.053    0.053    8.356)    8.356   3.650   (  -4.729    4.729    3.173)    7.402   5.943   (  -4.477    4.477    4.406)    7.714   5.945   (  -4.725    4.725    3.595)    7.588  10.571   ( -13.395   13.395   -9.041)   20.991======================= Grid point 25 (17/145) =======================q-point: ( 0.50  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.172   (   0.639   -0.639    2.600)    2.753   2.214   (   2.028   -2.028    6.759)    7.342   3.740   (  -1.676    1.676   -0.290)    2.388   6.038   (  -1.679    1.679   -0.676)    2.470   6.052   (  -1.645    1.645    1.584)    2.815  10.708   (  -9.294    9.294  -13.385)   18.759======================= Grid point 26 (18/145) =======================q-point: (-0.44  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.152   (   2.478   -2.478    0.327)    3.520   2.210   (   4.199   -4.199    4.983)    7.752   3.722   (   1.427   -1.427   -3.355)    3.915   6.007   (   1.708   -1.708   -4.167)    4.817   6.044   (   2.478   -2.478   -1.437)    3.788  10.683   (  -4.178    4.178  -17.677)   18.639======================= Grid point 27 (19/145) =======================q-point: (-0.38  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.066   (   4.242   -4.242   -1.761)    6.252   2.133   (   6.412   -6.412    2.893)    9.518   3.594   (   4.883   -4.883   -6.315)    9.358   5.878   (   4.451   -4.451   -5.032)    8.059   5.909   (   6.197   -6.197   -4.083)    9.668  10.508   (  -0.024    0.024  -19.942)   19.942======================= Grid point 28 (20/145) =======================q-point: (-0.31  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.916   (   6.191   -6.191   -3.506)    9.431   1.980   (   8.622   -8.622    0.662)   12.211   3.346   (   8.928   -8.928   -9.193)   15.618   5.684   (   7.741   -7.741   -5.714)   12.349   5.709   (   5.838   -5.838   -2.461)    8.615  10.256   (   1.588   -1.588  -18.916)   19.049======================= Grid point 29 (21/145) =======================q-point: (-0.25  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.698   (   8.526   -8.526   -4.797)   12.977   1.751   (  10.798  -10.798   -1.227)   15.319   2.967   (  13.492  -13.492  -11.491)   22.274   5.451   (   6.048   -6.048   -5.754)   10.309   5.573   (   5.638   -5.638    2.494)    8.354  10.033   (   0.295   -0.295  -14.735)   14.740======================= Grid point 30 (22/145) =======================q-point: (-0.19  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.411   (  11.073  -11.073   -5.377)   16.557   1.456   (  12.855  -12.855   -2.273)   18.321   2.463   (  18.124  -18.124  -12.107)   28.347   5.298   (   1.810   -1.810   -4.180)    4.902   5.512   (   4.324   -4.324    6.840)    9.175   9.920   (  -2.587    2.587   -9.052)    9.764======================= Grid point 31 (23/145) =======================q-point: (-0.12  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.069   (  13.349  -13.349   -4.476)   19.402   1.111   (  14.630  -14.630   -2.190)   20.805   1.865   (  22.335  -22.335   -9.388)   32.952   5.273   (  -2.447    2.447   -1.875)    3.936   5.516   (   2.800   -2.800    6.887)    7.944   9.935   (  -4.786    4.786   -3.893)    7.808======================= Grid point 32 (24/145) =======================q-point: (-0.06  0.06  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 620Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.709   (  14.874  -14.874    0.000)   21.034   0.742   (  15.745  -15.745    0.000)   22.267   1.245   (  25.574  -25.574    0.000)   36.166   5.347   (  -4.221    4.221    0.000)    5.970   5.512   (   1.694   -1.694    0.000)    2.395  10.029   (  -4.704    4.704    0.000)    6.653======================= Grid point 35 (25/145) =======================q-point: ( 0.12  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 349Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.008   (   0.000   -0.000   19.700)   19.700   1.008   (   0.000   -0.000   19.700)   19.700   1.704   (   0.000   -0.000   33.595)   33.595   5.543   (   0.000   -0.000    2.658)    2.658   5.543   (   0.000   -0.000    2.658)    2.658   9.675   (  -0.000    0.000  -20.225)   20.225======================= Grid point 36 (26/145) =======================q-point: ( 0.19  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.266   (  -3.302    3.302   17.226)   17.848   1.284   (  -4.534    4.534   17.845)   18.962   2.149   (  -6.151    6.151   28.855)   30.138   5.506   (   4.990   -4.990    3.952)    8.088   5.556   (   1.516   -1.516    3.317)    3.949   9.500   (  -6.738    6.738  -23.226)   25.104======================= Grid point 37 (27/145) =======================q-point: ( 0.25  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.541   (  -4.597    4.597   14.138)   15.561   1.602   (  -5.454    5.454   16.173)   17.918   2.619   (  -8.555    8.555   21.898)   25.017   5.464   (   2.740   -2.740    5.861)    7.026   5.567   (   1.066   -1.066    4.515)    4.760   9.455   ( -12.500   12.500  -22.733)   28.798======================= Grid point 38 (28/145) =======================q-point: ( 0.31  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.796   (  -4.553    4.553   11.295)   13.001   1.892   (  -3.931    3.931   14.833)   15.841   3.033   (  -8.711    8.711   14.610)   19.111   5.531   (  -1.979    1.979    6.973)    7.514   5.623   (  -1.039    1.039    5.501)    5.693   9.551   ( -17.368   17.368  -20.821)   32.200======================= Grid point 39 (29/145) =======================q-point: ( 0.38  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.006   (  -3.450    3.450    8.839)   10.096   2.122   (  -1.623    1.623   13.671)   13.862   3.352   (  -7.480    7.480    8.112)   13.330   5.694   (  -4.956    4.956    6.099)    9.291   5.740   (  -3.377    3.377    5.515)    7.295   9.758   ( -20.368   20.368  -20.102)   35.125======================= Grid point 40 (30/145) =======================q-point: ( 0.44  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.155   (  -1.577    1.577    6.635)    6.999   2.284   (   0.671   -0.671   12.529)   12.565   3.566   (  -5.511    5.511    2.964)    8.339   5.869   (  -4.946    4.946    3.152)    7.672   5.891   (  -4.347    4.347    4.360)    7.537   9.996   ( -20.417   20.417  -21.669)   36.101======================= Grid point 41 (31/145) =======================q-point: ( 0.50  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.230   (   0.624   -0.624    4.496)    4.581   2.380   (   2.890   -2.890   11.104)   11.832   3.677   (  -3.164    3.164   -0.898)    4.563   5.976   (  -2.783    2.783   -0.506)    3.969   6.021   (  -3.061    3.061    2.378)    4.939  10.172   ( -17.520   17.520  -24.740)   35.014======================= Grid point 42 (32/145) =======================q-point: (-0.44  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.231   (   2.732   -2.732    2.385)    4.540   2.405   (   5.116   -5.116    9.164)   11.676   3.692   (  -0.490    0.490   -3.826)    3.888   5.985   (   0.156   -0.156   -2.976)    2.984   6.070   (   0.305   -0.305    0.077)    0.438  10.221   ( -12.775   12.775  -27.437)   32.851======================= Grid point 43 (33/145) =======================q-point: (-0.38  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.162   (   4.553   -4.553    0.455)    6.455   2.354   (   7.327   -7.327    6.814)   12.402   3.610   (   2.642   -2.642   -6.113)    7.165   5.913   (   2.735   -2.735   -2.884)    4.824   6.002   (   4.614   -4.614   -2.045)    6.838  10.136   (  -7.845    7.845  -27.771)   29.905======================= Grid point 44 (34/145) =======================q-point: (-0.31  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.034   (   6.203   -6.203   -1.059)    8.837   2.224   (   9.453   -9.453    4.478)   14.099   3.427   (   6.394   -6.394   -7.718)   11.888   5.810   (   4.238   -4.238   -0.264)    5.999   5.819   (   8.099   -8.099   -3.445)   11.960   9.968   (  -4.352    4.352  -24.419)   25.182======================= Grid point 45 (35/145) =======================q-point: (-0.25  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.852   (   8.012   -8.012   -1.913)   11.491   2.023   (  11.429  -11.429    2.524)   16.359   3.136   (  10.772  -10.772   -8.017)   17.215   5.573   (   8.954   -8.954   -3.568)   13.156   5.724   (   4.453   -4.453    3.026)    6.987   9.805   (  -3.170    3.170  -17.385)   17.954======================= Grid point 46 (36/145) =======================q-point: (-0.19  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.620   (  10.182  -10.182   -1.775)   14.508   1.758   (  13.153  -13.153    0.938)   18.624   2.750   (  15.447  -15.447   -5.818)   22.607   5.355   (   6.372   -6.372   -2.215)    9.279   5.672   (   3.757   -3.757    3.914)    6.599   9.735   (  -3.852    3.852   -8.310)    9.936======================= Grid point 47 (37/145) =======================q-point: (-0.12  0.12  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 621Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.347   (  12.316  -12.316    0.000)   17.418   1.442   (  14.402  -14.402    0.000)   20.368   2.305   (  19.712  -19.712    0.000)   27.877   5.248   (   1.625   -1.625    0.000)    2.298   5.622   (   2.999   -2.999    0.000)    4.242   9.791   (  -4.886    4.886    0.000)    6.910======================= Grid point 52 (38/145) =======================q-point: ( 0.19  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 344Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.426   (   0.000   -0.000   16.319)   16.319   1.426   (   0.000   -0.000   16.319)   16.319   2.418   (   0.000   -0.000   27.872)   27.872   5.607   (   0.000   -0.000    2.747)    2.747   5.607   (   0.000   -0.000    2.747)    2.747   9.105   (  -0.000    0.000  -28.752)   28.752======================= Grid point 53 (39/145) =======================q-point: ( 0.25  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.628   (  -2.126    2.126   14.020)   14.339   1.660   (  -4.492    4.492   14.552)   15.878   2.753   (  -2.895    2.895   22.988)   23.349   5.596   (   3.064   -3.064    3.533)    5.591   5.630   (   0.575   -0.575    2.929)    3.040   8.870   (  -9.672    9.672  -30.829)   33.727======================= Grid point 54 (40/145) =======================q-point: ( 0.31  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.839   (  -3.184    3.184   11.485)   12.336   1.947   (  -5.778    5.778   13.490)   15.772   3.062   (  -4.053    4.053   16.085)   17.076   5.591   (   0.908   -0.908    4.779)    4.949   5.661   (  -0.314    0.314    3.479)    3.507   8.835   ( -17.618   17.618  -30.490)   39.375======================= Grid point 55 (41/145) =======================q-point: ( 0.38  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.034   (  -3.212    3.212    9.206)   10.266   2.219   (  -4.140    4.140   13.226)   14.464   3.306   (  -4.549    4.549    8.740)   10.852   5.670   (  -2.619    2.619    4.815)    6.074   5.731   (  -2.146    2.146    3.720)    4.801   8.958   ( -22.688   22.688  -29.992)   43.920======================= Grid point 56 (42/145) =======================q-point: ( 0.44  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.193   (  -2.158    2.158    7.206)    7.825   2.435   (  -1.242    1.242   13.141)   13.258   3.476   (  -4.607    4.607    2.413)    6.948   5.800   (  -4.046    4.046    3.236)    6.573   5.842   (  -3.751    3.751    3.256)    6.224   9.161   ( -24.620   24.620  -30.769)   46.465======================= Grid point 57 (43/145) =======================q-point: ( 0.50  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.294   (  -0.216    0.216    5.308)    5.316   2.577   (   1.797   -1.797   12.386)   12.644   3.576   (  -3.953    3.953   -2.131)    5.983   5.910   (  -3.001    3.001    1.065)    4.375   5.966   (  -3.927    3.927    2.159)    5.959   9.358   ( -23.490   23.490  -32.279)   46.320======================= Grid point 58 (44/145) =======================q-point: (-0.44  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.323   (   2.093   -2.093    3.424)    4.525   2.636   (   4.581   -4.581   10.535)   12.367   3.609   (  -2.387    2.387   -4.871)    5.927   5.957   (  -0.713    0.713   -0.184)    1.026   6.056   (  -2.016    2.016    0.754)    2.949   9.487   ( -19.964   19.964  -32.655)   43.168======================= Grid point 60 (45/145) =======================q-point: (-0.38  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.278   (   4.221   -4.221    1.674)    6.199   2.609   (   6.991   -6.991    7.931)   12.675   3.572   (   0.111   -0.111   -6.136)    6.138   5.945   (   1.544   -1.544    0.201)    2.192   6.062   (   1.726   -1.726   -0.592)    2.511   9.529   ( -15.195   15.195  -29.862)   36.790======================= Grid point 61 (46/145) =======================q-point: (-0.31  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.171   (   5.915   -5.915    0.317)    8.372   2.499   (   9.039   -9.039    5.213)   13.805   3.461   (   3.427   -3.427   -6.041)    7.745   5.900   (   3.007   -3.007    1.473)    4.500   5.959   (   6.075   -6.075   -1.455)    8.714   9.517   ( -10.504   10.504  -22.715)   27.142======================= Grid point 62 (47/145) =======================q-point: (-0.25  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.017   (   7.300   -7.300   -0.350)   10.329   2.314   (  10.794  -10.794    2.534)   15.474   3.276   (   7.360   -7.360   -4.139)   11.201   5.761   (   9.195   -9.195   -1.351)   13.073   5.844   (   3.513   -3.513    1.880)    5.312   9.514   (  -6.958    6.958  -11.804)   15.367======================= Grid point 63 (48/145) =======================q-point: (-0.19  0.19  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 620Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.830   (   8.699   -8.699    0.000)   12.302   2.061   (  12.244  -12.244    0.000)   17.316   3.031   (  11.569  -11.569    0.000)   16.362   5.531   (   9.409   -9.409    0.000)   13.307   5.776   (   3.502   -3.502    0.000)    4.953   9.581   (  -4.912    4.912    0.000)    6.946======================= Grid point 70 (49/145) =======================q-point: ( 0.25  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 349Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.759   (   0.000   -0.000   12.464)   12.464   1.759   (   0.000   -0.000   12.464)   12.464   2.979   (   0.000   -0.000   20.391)   20.391   5.664   (   0.000   -0.000    2.052)    2.052   5.664   (   0.000   -0.000    2.052)    2.052   8.367   (  -0.000    0.000  -34.481)   34.481======================= Grid point 71 (50/145) =======================q-point: ( 0.31  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.911   (  -1.617    1.617   10.376)   10.626   1.954   (  -5.008    5.008   10.856)   12.962   3.197   (  -0.649    0.649   15.074)   15.102   5.664   (   1.437   -1.437    2.318)    3.083   5.686   (  -0.188    0.188    1.884)    1.903   8.101   ( -11.637   11.637  -34.969)   38.648======================= Grid point 72 (51/145) =======================q-point: ( 0.38  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.069   (  -2.592    2.592    8.377)    9.144   2.222   (  -7.096    7.096   10.243)   14.340   3.348   (  -0.697    0.697    8.331)    8.389   5.681   (  -0.560    0.560    3.239)    3.335   5.724   (  -1.304    1.304    1.954)    2.687   8.075   ( -20.098   20.098  -34.458)   44.668======================= Grid point 73 (52/145) =======================q-point: ( 0.44  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.218   (  -2.629    2.629    6.643)    7.612   2.496   (  -5.797    5.797   10.573)   13.379   3.427   (  -1.378    1.378    1.655)    2.556   5.763   (  -2.644    2.644    3.819)    5.345   5.794   (  -2.862    2.862    1.811)    4.434   8.200   ( -24.587   24.587  -34.625)   49.071======================= Grid point 74 (53/145) =======================q-point: ( 0.50  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.335   (  -1.490    1.490    5.035)    5.459   2.717   (  -2.184    2.184   11.020)   11.445   3.460   (  -2.674    2.674   -3.680)    5.276   5.872   (  -2.572    2.572    3.841)    5.291   5.893   (  -3.970    3.970    1.339)    5.772   8.387   ( -25.915   25.915  -34.781)   50.526======================= Grid point 75 (54/145) =======================q-point: ( 0.56  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.396   (   0.596   -0.596    3.411)    3.514   2.840   (   2.253   -2.253   10.121)   10.610   3.472   (  -3.549    3.549   -6.465)    8.184   5.957   (  -0.882    0.882    3.671)    3.877   5.997   (  -3.660    3.660    0.676)    5.220   8.583   ( -24.674   24.674  -32.566)   47.729======================= Grid point 76 (55/145) =======================q-point: (-0.38  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.385   (   2.968   -2.968    1.845)    4.584   2.846   (   5.761   -5.761    7.248)   10.905   3.472   (  -2.619    2.619   -6.134)    7.165   5.995   (   1.076   -1.076    3.201)    3.544   6.063   (  -1.389    1.389    0.027)    1.964   8.776   ( -21.354   21.354  -25.740)   39.680======================= Grid point 77 (56/145) =======================q-point: (-0.31  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.305   (   4.986   -4.986    0.601)    7.077   2.749   (   7.959   -7.959    3.573)   11.809   3.445   (   0.124   -0.124   -3.684)    3.688   5.984   (   2.413   -2.413    1.939)    3.925   6.050   (   2.474   -2.474   -0.330)    3.514   8.984   ( -16.657   16.657  -14.012)   27.409======================= Grid point 78 (57/145) =======================q-point: (-0.25  0.25  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 621Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.175   (   6.388   -6.388    0.000)    9.034   2.566   (   9.541   -9.541    0.000)   13.493   3.380   (   3.675   -3.675    0.000)    5.197   5.931   (   3.041   -3.041    0.000)    4.301   5.941   (   6.597   -6.597    0.000)    9.329   9.229   ( -11.558   11.558    0.000)   16.345======================= Grid point 87 (58/145) =======================q-point: ( 0.31  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 344Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.002   (   0.000   -0.000    8.631)    8.631   2.002   (   0.000   -0.000    8.631)    8.631   3.348   (   0.000   -0.000   11.227)   11.227   5.700   (   0.000   -0.000    1.064)    1.064   5.700   (   0.000   -0.000    1.064)    1.064   7.551   (  -0.000    0.000  -34.971)   34.971======================= Grid point 88 (59/145) =======================q-point: ( 0.38  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.109   (  -1.493    1.493    6.854)    7.171   2.163   (  -5.795    5.795    7.260)   10.949   3.438   (   0.824   -0.824    5.612)    5.731   5.707   (   0.036   -0.036    1.593)    1.594   5.717   (  -0.667    0.667    0.802)    1.238   7.300   ( -11.737   11.737  -33.055)   36.989======================= Grid point 89 (60/145) =======================q-point: ( 0.44  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.227   (  -2.498    2.498    5.324)    6.389   2.419   (  -8.807    8.807    6.836)   14.207   3.443   (   1.333   -1.333    0.066)    1.887   5.754   (  -1.853    1.853    0.699)    2.712   5.759   (  -1.921    1.921    4.055)    4.881   7.290   ( -19.174   19.174  -32.345)   42.208======================= Grid point 90 (61/145) =======================q-point: ( 0.50  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.341   (  -2.495    2.495    3.963)    5.306   2.703   (  -8.186    8.186    7.207)   13.637   3.391   (   0.792   -0.792   -4.101)    4.251   5.820   (  -3.221    3.221    0.537)    4.587   5.878   (  -2.433    2.433    6.768)    7.592   7.414   ( -23.530   23.530  -31.978)   46.151======================= Grid point 91 (62/145) =======================q-point: ( 0.56  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.424   (  -1.215    1.215    2.579)    3.099   2.940   (  -3.976    3.976    8.362)   10.076   3.325   (  -1.120    1.120   -7.256)    7.427   5.911   (  -4.070    4.070    0.314)    5.764   6.002   (  -0.791    0.791    7.170)    7.257   7.640   ( -25.906   25.906  -27.395)   45.746======================= Grid point 92 (63/145) =======================q-point: (-0.38  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.449   (   0.968   -0.968    1.175)    1.804   3.037   (   3.806   -3.806    6.382)    8.349   3.311   (  -5.072    5.072   -6.158)    9.453   6.004   (  -3.535    3.535    0.093)    5.000   6.064   (   1.094   -1.094    4.031)    4.318   7.985   ( -25.513   25.513  -15.652)   39.330======================= Grid point 93 (64/145) =======================q-point: (-0.31  0.31  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 620Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.406   (   3.257   -3.257    0.000)    4.606   2.941   (   6.689   -6.689    0.000)    9.459   3.385   (  -3.362    3.362    0.000)    4.754   6.049   (   2.025   -2.025    0.000)    2.863   6.063   (  -1.182    1.182    0.000)    1.671   8.449   ( -22.066   22.066    0.000)   31.206======================= Grid point 104 (65/145) =======================q-point: ( 0.38  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 349Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.161   (   0.000   -0.000    5.205)    5.205   2.161   (   0.000   -0.000    5.205)    5.205   3.490   (   0.000   -0.000    1.155)    1.155   5.714   (   0.000   -0.000    0.293)    0.293   5.714   (   0.000   -0.000    0.293)    0.293   6.812   (  -0.000    0.000  -27.575)   27.575======================= Grid point 105 (66/145) =======================q-point: ( 0.44  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.231   (  -1.561    1.561    3.773)    4.371   2.293   (  -6.576    6.576    4.049)   10.143   3.463   (   1.318   -1.318   -2.731)    3.307   5.726   (  -0.887    0.887    0.141)    1.262   5.749   (  -2.501    2.501    2.299)    4.219   6.630   (  -8.748    8.748  -23.934)   26.942======================= Grid point 106 (67/145) =======================q-point: ( 0.50  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.317   (  -2.607    2.607    2.544)    4.479   2.538   ( -10.310   10.310    3.424)   14.978   3.383   (   1.930   -1.930   -4.003)    4.844   5.761   (  -2.071    2.071    0.096)    2.930   5.895   (  -4.394    4.394    8.267)   10.341   6.624   ( -14.957   14.957  -24.552)   32.407======================= Grid point 107 (68/145) =======================q-point: ( 0.56  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.401   (  -2.520    2.520    1.311)    3.798   2.823   ( -10.714   10.714    2.899)   15.427   3.289   (   2.519   -2.519   -3.311)    4.863   5.825   (  -3.348    3.348    0.062)    4.736   6.078   (  -0.932    0.932    9.051)    9.147   6.797   ( -23.124   23.124  -18.668)   37.656======================= Grid point 108 (69/145) =======================q-point: (-0.38  0.38  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 621Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.454   (  -1.147    1.147    0.000)    1.622   3.082   (  -9.450    9.450    0.000)   13.364   3.188   (   4.071   -4.071    0.000)    5.757   5.914   (  -4.096    4.096    0.000)    5.793   6.120   (   1.115   -1.115    0.000)    1.577   7.278   ( -26.843   26.843    0.000)   37.961======================= Grid point 123 (70/145) =======================q-point: ( 0.44  0.44  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 344Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.247   (   0.000   -0.000    2.376)    2.376   2.247   (   0.000   -0.000    2.376)    2.376   3.434   (  -0.000    0.000   -4.471)    4.471   5.717   (  -0.000    0.000   -0.022)    0.022   5.717   (  -0.000    0.000   -0.022)    0.022   6.327   (  -0.000    0.000  -13.896)   13.896======================= Grid point 124 (71/145) =======================q-point: ( 0.50  0.44  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.287   (  -1.652    1.652    1.185)    2.620   2.353   (  -7.061    7.061    1.283)   10.068   3.375   (   1.010   -1.010   -3.081)    3.396   5.727   (  -0.952    0.952   -0.023)    1.346   5.815   (  -7.965    7.965    2.857)   11.621   6.227   (  -2.093    2.093  -10.118)   10.542======================= Grid point 125 (72/145) =======================q-point: (-0.44 -0.56  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 620Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.346   (  -2.676    2.676    0.000)    3.784   2.577   ( -10.950   10.950    0.000)   15.486   3.324   (   1.884   -1.884    0.000)    2.664   5.762   (  -2.119    2.119    0.000)    2.997   6.157   (   0.513   -0.513    0.000)    0.726   6.171   ( -18.706   18.706    0.000)   26.454======================= Grid point 143 (73/145) =======================q-point: (-0.50 -0.50  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 205Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.274   (   0.000   -0.000    0.000)    0.000   2.274   (   0.000   -0.000    0.000)    0.000   3.367   (   0.000   -0.000    0.000)    0.000   5.716   (   0.000   -0.000    0.000)    0.000   5.716   (   0.000   -0.000    0.000)    0.000   6.169   (   0.000   -0.000    0.000)    0.000======================= Grid point 304 (74/145) =======================q-point: ( 0.19  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.103   (   0.000    7.969   17.733)   19.441   1.141   (   0.000   10.586   17.620)   20.555   1.890   (   0.000   14.699   28.470)   32.040   5.358   (  -0.000  -10.424    3.836)   11.107   5.565   (   0.000    1.840    2.751)    3.310   9.776   (  -0.000    5.208  -16.570)   17.369======================= Grid point 305 (75/145) =======================q-point: ( 0.25  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.406   (  -2.628    8.118   14.920)   17.187   1.470   (  -1.222   10.995   15.710)   19.215   2.398   (  -6.495   14.140   22.652)   27.481   5.329   (  -1.562   -4.880    5.960)    7.860   5.592   (   5.129    2.930    3.597)    6.916   9.692   (  -6.721    6.858  -16.067)   18.718======================= Grid point 306 (76/145) =======================q-point: ( 0.31  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.683   (  -3.221    6.817   11.991)   14.165   1.773   (  -0.528    9.530   14.038)   16.975   2.867   (  -8.977   11.966   16.219)   22.065   5.427   (  -3.570    2.354    7.997)    9.069   5.622   (   4.373    4.075    4.775)    7.650   9.729   ( -14.160    8.975  -14.024)   21.857======================= Grid point 307 (77/145) =======================q-point: ( 0.38  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.910   (  -2.810    4.804    9.557)   11.059   2.021   (   1.257    7.618   12.284)   14.509   3.245   (  -8.764    8.997   10.178)   16.166   5.629   (  -4.469    7.030    8.122)   11.635   5.704   (   0.374    4.691    5.354)    7.128   9.892   ( -20.563   10.235  -12.755)   26.273======================= Grid point 308 (78/145) =======================q-point: ( 0.44  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.079   (  -1.549    2.843    7.509)    8.178   2.202   (   3.510    5.479   10.572)   12.414   3.510   (  -7.098    5.847    5.098)   10.515   5.827   (  -3.457    4.793    4.212)    7.257   5.856   (  -3.672    6.968    6.147)    9.991  10.117   ( -24.118    9.136  -13.831)   29.265======================= Grid point 309 (79/145) =======================q-point: ( 0.50  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.184   (   0.445    1.272    5.376)    5.542   2.310   (   5.716    2.919    9.025)   11.074   3.664   (  -4.797    2.846    1.081)    5.681   5.952   (  -4.377    2.514    1.729)    5.335   6.013   (  -1.967    4.644    2.717)    5.729  10.308   ( -24.104    5.394  -17.145)   30.067======================= Grid point 310 (80/145) =======================q-point: (-0.44  0.12  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1032Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.220   (   2.914    0.000    2.914)    4.120   2.346   (   7.521   -0.000    7.521)   10.636   3.715   (  -2.116    0.000   -2.116)    2.993   6.006   (  -1.976    0.000   -1.976)    2.795   6.066   (   0.049    0.000    0.049)    0.069  10.382   ( -21.224    0.000  -21.224)   30.015======================= Grid point 322 (81/145) =======================q-point: ( 0.25  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.495   (   0.000    5.444   15.839)   16.749   1.527   (   0.000    8.260   15.426)   17.499   2.513   (   0.000    7.848   24.786)   25.999   5.474   (  -0.000   -8.095    5.480)    9.776   5.631   (   0.000    2.037    2.852)    3.505   9.283   (  -0.000   12.514  -25.758)   28.637======================= Grid point 323 (82/145) =======================q-point: ( 0.31  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.735   (  -0.630    5.956   13.283)   14.571   1.816   (  -2.508    8.758   13.908)   16.626   2.879   (  -3.807    7.613   18.455)   20.324   5.491   (  -1.212   -2.185    7.066)    7.494   5.673   (   3.093    3.082    3.193)    5.410   9.206   (  -9.537   16.127  -25.620)   31.740======================= Grid point 324 (83/145) =======================q-point: ( 0.38  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.949   (  -1.039    5.200   10.755)   11.991   2.090   (  -1.593    7.007   13.147)   14.982   3.191   (  -5.696    6.659   11.483)   14.444   5.611   (  -2.567    4.251    7.083)    8.650   5.722   (   1.806    3.879    3.731)    5.677   9.276   ( -18.492   17.823  -24.901)   35.773======================= Grid point 325 (84/145) =======================q-point: ( 0.44  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.124   (  -0.692    3.663    8.600)    9.373   2.311   (   0.870    4.584   12.580)   13.418   3.425   (  -6.095    5.294    5.213)    9.610   5.778   (  -2.317    6.089    4.246)    7.776   5.816   (  -1.918    5.011    4.392)    6.934   9.446   ( -24.810   17.244  -25.454)   39.507======================= Grid point 326 (85/145) =======================q-point: ( 0.50  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.244   (   0.699    1.865    6.547)    6.843   2.463   (   3.677    1.915   11.759)   12.469   3.575   (  -5.406    3.558    0.346)    6.481   5.892   (  -3.179    3.315    1.569)    4.854   5.958   (  -2.408    5.720    2.814)    6.815   9.630   ( -27.476   14.161  -27.528)   41.392======================= Grid point 327 (86/145) =======================q-point: ( 0.56  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.298   (   2.903    0.085    4.390)    5.264   2.540   (   6.186   -1.000   10.289)   12.048   3.642   (  -3.813    1.509   -3.057)    5.115   5.957   (  -2.275    1.029   -0.785)    2.617   6.050   (  -1.442    2.360    0.802)    2.880   9.750   ( -26.702    9.097  -29.822)   41.050======================= Grid point 328 (87/145) =======================q-point: (-0.38  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.280   (   5.324   -1.555    2.277)    5.996   2.534   (   8.151   -3.993    8.066)   12.143   3.627   (  -1.266   -0.815   -5.235)    5.447   5.948   (   0.276   -1.457   -1.887)    2.400   6.058   (   0.573   -2.226   -0.682)    2.398   9.763   ( -23.680    3.191  -30.380)   38.651======================= Grid point 329 (88/145) =======================q-point: (-0.31  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.195   (   7.474   -3.007    0.541)    8.075   2.446   (   9.644   -6.814    5.478)   13.017   3.530   (   2.353   -3.405   -6.388)    7.612   5.887   (   2.298   -3.477   -0.654)    4.219   5.965   (   4.153   -6.381   -1.944)    7.858   9.682   ( -20.085   -2.203  -27.508)   34.131======================= Grid point 330 (89/145) =======================q-point: (-0.25  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.055   (   9.334   -4.342   -0.548)   10.309   2.283   (  10.888   -9.312    3.072)   14.652   3.346   (   7.055   -6.225   -6.326)   11.338   5.775   (   7.251   -8.663   -2.200)   11.510   5.818   (   3.702   -4.203    1.369)    5.765   9.567   ( -17.457   -5.965  -20.522)   27.595======================= Grid point 331 (90/145) =======================q-point: (-0.19  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.869   (  11.217   -5.752   -0.822)   12.632   2.054   (  12.144  -11.305    1.106)   16.629   3.080   (  12.579   -9.165   -4.381)   16.169   5.554   (   9.338   -7.183   -1.747)   11.910   5.756   (   3.357   -3.929    2.505)    5.743   9.500   ( -16.484   -7.659  -10.410)   20.946======================= Grid point 332 (91/145) =======================q-point: (-0.12  0.19  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.647   (  13.253   -7.254   -0.000)   15.108   1.773   (  13.710  -12.201   -0.000)   18.353   2.752   (  18.178  -11.988   -0.000)   21.776   5.384   (   8.495   -1.706   -0.000)    8.664   5.695   (   3.134   -3.329    0.000)    4.572   9.542   ( -16.546   -7.705    0.000)   18.252======================= Grid point 339 (92/145) =======================q-point: ( 0.31  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.827   (   0.000    5.056   12.770)   13.735   1.845   (   0.000    7.136   12.031)   13.988   3.015   (   0.000    3.093   18.249)   18.509   5.590   (   0.000   -4.002    4.315)    5.885   5.690   (   0.000    2.236    2.148)    3.100   8.603   (  -0.000   17.740  -32.485)   37.014======================= Grid point 340 (93/145) =======================q-point: ( 0.38  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.010   (   1.403    5.907   10.383)   12.028   2.107   (  -4.976    7.674   11.256)   14.503   3.233   (  -1.618    2.619   11.826)   12.220   5.631   (  -0.980    1.034    4.927)    5.129   5.734   (   1.457    3.231    2.096)    4.117   8.527   ( -11.468   21.915  -32.490)   40.834======================= Grid point 341 (94/145) =======================q-point: ( 0.44  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.171   (   0.786    5.229    8.324)    9.861   2.374   (  -4.260    5.996   11.291)   13.475   3.385   (  -2.828    2.538    5.050)    6.320   5.739   (  -1.500    4.969    4.155)    6.649   5.792   (  -0.219    3.985    2.306)    4.609   8.602   ( -20.884   22.633  -32.702)   44.920======================= Grid point 342 (95/145) =======================q-point: ( 0.50  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.299   (   1.052    3.647    6.509)    7.535   2.592   (  -1.196    3.011   11.485)   11.933   3.478   (  -3.983    2.650   -0.796)    4.849   5.845   (  -1.592    3.701    2.299)    4.639   5.894   (  -2.329    5.613    2.314)    6.503   8.750   ( -26.678   20.719  -33.565)   47.620======================= Grid point 343 (96/145) =======================q-point: ( 0.56  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.375   (   2.499    1.561    4.728)    5.571   2.727   (   2.710   -0.544   10.691)   11.043   3.524   (  -4.600    2.351   -4.664)    6.960   5.921   (  -1.946    1.578    1.622)    2.984   5.999   (  -2.247    4.169    0.995)    4.839   8.902   ( -28.764   16.766  -33.769)   47.422======================= Grid point 344 (97/145) =======================q-point: (-0.38  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.384   (   4.723   -0.556    2.964)    5.604   2.760   (   5.876   -4.134    8.338)   11.007   3.531   (  -3.651    1.269   -6.064)    7.191   5.957   (  -0.424   -0.833    1.533)    1.795   6.057   (  -1.495    0.736   -0.050)    1.667   9.016   ( -27.868   11.396  -31.216)   43.370======================= Grid point 345 (98/145) =======================q-point: (-0.31  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.324   (   7.014   -2.307    1.411)    7.517   2.693   (   7.675   -7.290    5.255)   11.818   3.494   (  -0.737   -0.643   -5.490)    5.576   5.950   (   1.267   -2.779    1.876)    3.584   6.040   (   1.074   -3.599   -0.751)    3.830   9.095   ( -25.138    5.402  -24.280)   35.364======================= Grid point 346 (99/145) =======================q-point: (-0.25  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.207   (   8.869   -3.580    0.325)    9.569   2.542   (   8.825   -9.754    2.382)   13.368   3.405   (   3.472   -3.155   -3.444)    5.819   5.911   (   2.392   -4.135    1.618)    5.044   5.931   (   4.760   -7.102   -0.865)    8.593   9.176   ( -21.819   -0.288  -13.022)   25.411======================= Grid point 347 (100/145) =======================q-point: (-0.19  0.25  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.049   (  10.201   -4.612   -0.000)   11.195   2.323   (  10.224  -11.242    0.000)   15.196   3.263   (   8.132   -6.066   -0.000)   10.145   5.747   (   8.170   -9.186    0.000)   12.294   5.850   (   2.726   -3.862    0.000)    4.727   9.305   ( -18.854   -4.804    0.000)   19.456======================= Grid point 356 (101/145) =======================q-point: ( 0.38  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.082   (   0.000    6.570    8.419)   10.679   2.084   (   0.000    5.950    9.402)   11.126   3.342   (   0.000   -0.235    9.812)    9.815   5.671   (   0.000   -0.842    2.774)    2.899   5.728   (   0.000    2.399    1.136)    2.654   7.815   (  -0.000   20.128  -34.820)   40.219======================= Grid point 357 (102/145) =======================q-point: ( 0.44  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.212   (   2.597    6.508    7.053)    9.942   2.334   (  -7.226    7.958    8.374)   13.626   3.414   (  -0.033   -0.757    3.758)    3.834   5.726   (  -0.907    2.918    3.677)    4.781   5.769   (   0.382    3.365    0.977)    3.525   7.743   ( -11.940   23.473  -34.342)   43.277======================= Grid point 358 (103/145) =======================q-point: ( 0.50  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.331   (   1.733    5.648    5.504)    8.074   2.606   (  -7.169    6.665    8.638)   13.055   3.421   (  -0.650   -0.290   -1.809)    1.944   5.823   (  -0.786    3.364    2.243)    4.120   5.835   (  -1.572    4.997    3.180)    6.128   7.813   ( -20.925   23.496  -34.234)   46.496======================= Grid point 359 (104/145) =======================q-point: ( 0.56  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.422   (   2.032    3.868    4.052)    5.959   2.832   (  -3.585    2.941    9.217)   10.317   3.393   (  -2.573    0.899   -6.121)    6.701   5.901   (  -2.584    2.882    1.812)    4.274   5.947   (  -1.064    3.728    3.829)    5.449   7.964   ( -26.484   21.704  -32.567)   47.255======================= Grid point 360 (105/145) =======================q-point: (-0.38  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.460   (   3.629    1.597    2.566)    4.723   2.945   (   2.248   -2.751    8.016)    8.768   3.377   (  -5.133    2.679   -7.256)    9.283   5.982   (  -2.681    0.849    2.462)    3.738   6.024   (  -0.567    1.551    2.270)    2.807   8.173   ( -28.883   18.227  -26.400)   43.167======================= Grid point 362 (106/145) =======================q-point: (-0.31  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.432   (   5.848   -0.606    1.116)    5.984   2.906   (   5.563   -7.095    3.710)    9.750   3.401   (  -4.105    2.433   -3.917)    6.173   6.020   (   0.461   -2.353    2.787)    3.676   6.052   (  -1.381   -0.308   -0.342)    1.456   8.446   ( -28.420   12.911  -14.682)   34.496======================= Grid point 363 (107/145) =======================q-point: (-0.25  0.31  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.341   (   7.756   -2.368   -0.000)    8.110   2.758   (   7.010   -9.070    0.000)   11.464   3.420   (  -0.415   -0.069    0.000)    0.420   5.987   (   1.808   -3.303    0.000)    3.765   6.029   (   1.050   -4.200    0.000)    4.329   8.787   ( -25.735    6.343    0.000)   26.506======================= Grid point 374 (108/145) =======================q-point: ( 0.44  0.38  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.237   (   0.000    6.297    5.110)    8.110   2.262   (   0.000    7.463    6.075)    9.623   3.462   (  -0.000   -2.057    0.696)    2.172   5.731   (   0.000    1.977    2.834)    3.456   5.744   (   0.000    2.512    0.338)    2.534   7.046   (  -0.000   18.384  -30.484)   35.599======================= Grid point 375 (109/145) =======================q-point: ( 0.50  0.38  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.338   (   2.748    6.937    3.945)    8.440   2.492   (  -8.589    9.294    5.244)   13.698   3.415   (   0.793   -2.414   -3.017)    3.944   5.781   (  -0.145    3.581    0.256)    3.593   5.830   (  -1.552    4.217    5.935)    7.444   6.992   ( -10.621   20.769  -29.558)   37.654======================= Grid point 376 (110/145) =======================q-point: ( 0.56  0.38  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.426   (   2.000    6.035    2.723)    6.917   2.768   (  -9.810    8.572    5.220)   14.035   3.329   (   0.901   -2.440   -4.989)    5.626   5.839   (  -1.630    4.167    0.193)    4.479   5.978   (  -1.254    2.602    8.385)    8.868   7.085   ( -19.530   22.555  -27.084)   40.295======================= Grid point 377 (111/145) =======================q-point: ( 0.62  0.38  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.485   (   2.483    4.164    1.404)    5.047   3.018   (  -6.867    4.944    6.761)   10.831   3.229   (  -0.376   -1.733   -6.976)    7.198   5.917   (  -3.154    4.031    0.150)    5.120   6.079   (   0.227   -0.727    5.411)    5.464   7.359   ( -26.298   22.891  -16.399)   38.530======================= Grid point 378 (112/145) =======================q-point: (-0.31  0.38  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.491   (   4.085    1.753   -0.000)    4.446   3.066   (   4.203   -6.473    0.000)    7.718   3.258   (  -7.728    5.730    0.000)    9.620   5.999   (  -3.526    2.627    0.000)    4.397   6.082   (   0.943   -2.216    0.000)    2.408   7.834   ( -29.130   19.280    0.000)   34.933======================= Grid point 391 (113/145) =======================q-point: ( 0.50  0.44  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.322   (   0.000    6.184    2.343)    6.613   2.366   (   0.000    8.851    2.949)    9.329   3.407   (  -0.000   -2.209   -4.111)    4.667   5.747   (   0.000    2.574   -0.002)    2.574   5.820   (   0.000    7.069    5.168)    8.757   6.451   (  -0.000   11.393  -20.349)   23.322======================= Grid point 392 (114/145) =======================q-point: ( 0.56  0.44  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.398   (   2.555    7.137    1.253)    7.684   2.575   (  -9.251   10.606    1.822)   14.191   3.330   (   0.935   -2.856   -2.840)    4.134   5.783   (  -0.285    3.609    0.001)    3.620   6.029   (  -2.852    5.913   11.590)   13.320   6.411   (  -8.433   16.541  -20.226)   27.456======================= Grid point 393 (115/145) =======================q-point: (-0.38 -0.56  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.458   (   2.029    6.217   -0.000)    6.540   2.833   ( -11.435   10.280    0.000)   15.376   3.257   (   1.914   -3.999    0.000)    4.433   5.840   (  -1.729    4.226   -0.000)    4.566   6.137   (   0.419   -1.257    0.000)    1.325   6.702   ( -20.469   24.869   -0.000)   32.210======================= Grid point 410 (116/145) =======================q-point: (-0.44 -0.50  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 580Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.349   (   0.000    6.155    0.000)    6.155   2.400   (   0.000    9.419    0.000)    9.419   3.346   (  -0.000   -1.876    0.000)    1.876   5.747   (   0.000    2.594    0.000)    2.594   5.932   (   0.000   17.635    0.000)   17.635   6.163   (  -0.000   -0.551    0.000)    0.551======================= Grid point 608 (117/145) =======================q-point: ( 0.38  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.947   (   0.000    4.723   12.493)   13.356   2.062   (   0.000   10.910   11.113)   15.573   3.120   (   0.000    5.810   13.312)   14.525   5.580   (   0.000    3.933    7.602)    8.559   5.760   (   0.000    3.658    2.400)    4.375   9.038   (  -0.000   17.093  -26.700)   31.702======================= Grid point 609 (118/145) =======================q-point: ( 0.44  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.127   (   0.630    3.817   10.806)   11.478   2.310   (  -0.575    9.853    9.978)   14.035   3.315   (  -2.788    4.172    7.449)    8.982   5.739   (  -1.057    7.578    6.115)    9.795   5.822   (   1.343    4.119    2.280)    4.896   8.984   ( -12.028   15.499  -27.648)   33.901======================= Grid point 610 (119/145) =======================q-point: ( 0.50  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2052Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.268   (   1.376    2.926    8.775)    9.352   2.522   (   1.112    7.131    9.753)   12.133   3.451   (  -4.573    2.759    1.858)    5.655   5.867   (   0.822    4.544    1.882)    4.987   5.911   (  -1.491    6.578    3.273)    7.497   9.017   ( -22.489   11.707  -29.254)   38.712======================= Grid point 611 (120/145) =======================q-point: ( 0.56  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.362   (   2.607    1.628    6.613)    7.292   2.666   (   4.163    3.736    9.193)   10.762   3.529   (  -5.428    1.382   -2.541)    6.151   5.918   (  -0.227    2.508    0.717)    2.618   6.017   (  -1.337    4.083    0.772)    4.365   9.076   ( -29.071    6.287  -31.036)   42.987======================= Grid point 612 (121/145) =======================q-point: (-0.38  0.25  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1034Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.395   (   4.446    0.000    4.446)    6.288   2.718   (   7.296   -0.000    7.296)   10.318   3.553   (  -4.982    0.000   -4.982)    7.046   5.940   (  -0.322    0.000   -0.322)    0.455   6.052   (  -0.692    0.000   -0.692)    0.978   9.103   ( -31.539    0.000  -31.539)   44.603======================= Grid point 626 (122/145) =======================q-point: ( 0.44  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.216   (   0.000    4.602   10.639)   11.591   2.282   (   0.000   10.188    7.942)   12.918   3.352   (   0.000    1.287    6.524)    6.650   5.720   (   0.000    5.269    4.642)    7.022   5.804   (   0.000    3.925    1.325)    4.143   8.334   (  -0.000   22.023  -33.304)   39.928======================= Grid point 627 (123/145) =======================q-point: ( 0.50  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.354   (   3.336    4.562    8.479)   10.190   2.517   (  -4.350    8.294    7.972)   12.299   3.414   (  -1.362    0.616    0.969)    1.782   5.839   (   0.079    5.905    3.040)    6.642   5.863   (  -0.064    4.801    1.409)    5.003   8.260   ( -12.663   19.043  -33.853)   40.853======================= Grid point 628 (124/145) =======================q-point: ( 0.56  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.446   (   3.771    3.538    6.417)    8.241   2.733   (  -2.479    4.933    8.424)   10.072   3.426   (  -3.174    0.404   -4.026)    5.142   5.902   (   0.621    3.217    1.528)    3.615   5.959   (  -1.310    4.923    1.387)    5.280   8.271   ( -22.862   14.310  -33.701)   43.165======================= Grid point 629 (125/145) =======================q-point: ( 0.62  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.493   (   4.746    1.812    4.657)    6.892   2.861   (   2.195    0.192    7.420)    7.740   3.417   (  -5.513    0.745   -6.713)    8.718   5.943   (  -0.091    0.886    1.921)    2.118   6.023   (  -1.226    1.990    0.131)    2.341   8.330   ( -29.072    8.591  -31.259)   43.545======================= Grid point 630 (126/145) =======================q-point: (-0.31  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.483   (   6.591    0.069    3.094)    7.282   2.849   (   5.716   -5.001    3.740)    8.466   3.422   (  -5.419    0.991   -5.374)    7.696   5.963   (   0.163   -1.933    2.379)    3.069   6.035   (  -0.798   -1.344   -0.702)    1.714   8.418   ( -31.587    2.193  -24.744)   40.184======================= Grid point 631 (127/145) =======================q-point: (-0.25  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.414   (   8.788   -1.224    1.554)    9.008   2.719   (   6.525   -8.781    0.703)   10.962   3.426   (  -1.778   -0.043   -2.645)    3.187   5.948   (   0.979   -3.812    1.928)    4.383   5.990   (   0.811   -4.874   -0.935)    5.028   8.545   ( -31.471   -4.521  -13.641)   34.596======================= Grid point 632 (128/145) =======================q-point: (-0.19  0.31  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.291   (  10.407   -2.283   -0.000)   10.655   2.528   (   7.437  -10.617    0.000)   12.963   3.394   (   2.905   -1.891   -0.000)    3.467   5.883   (   3.624   -7.716    0.000)    8.525   5.900   (   1.567   -4.154    0.000)    4.439   8.732   ( -29.877  -10.774    0.000)   31.761======================= Grid point 643 (129/145) =======================q-point: ( 0.50  0.38  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.430   (   0.000    9.813    4.888)   10.963   2.432   (   0.000    6.007    7.932)    9.950   3.417   (  -0.000   -1.544   -0.826)    1.751   5.816   (   0.000    5.078    4.399)    6.718   5.823   (   0.000    4.110    0.456)    4.135   7.547   (  -0.000   22.469  -33.534)   40.365======================= Grid point 644 (130/145) =======================q-point: ( 0.56  0.38  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.509   (   5.789    6.721    4.885)   10.126   2.682   (  -8.253    7.123    6.276)   12.580   3.365   (  -0.263   -1.883   -4.727)    5.095   5.876   (   0.054    4.295    0.434)    4.317   5.929   (  -0.464    3.997    4.883)    6.327   7.484   ( -12.230   19.490  -31.622)   39.108======================= Grid point 645 (131/145) =======================q-point: ( 0.62  0.38  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2052Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.561   (   5.435    4.945    3.499)    8.138   2.910   (  -5.795    3.633    6.957)    9.756   3.282   (  -2.198   -1.460   -7.736)    8.174   5.933   (  -0.837    3.519    0.663)    3.677   6.017   (  -0.231    0.822    4.323)    4.407   7.552   ( -22.160   15.668  -26.060)   37.626======================= Grid point 646 (132/145) =======================q-point: (-0.31  0.38  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.572   (   6.320    2.933    1.954)    7.236   3.000   (   2.205   -5.206    4.235)    7.064   3.259   (  -7.669    3.137   -5.504)    9.947   5.987   (  -1.683    1.538    0.798)    2.415   6.047   (   0.240   -1.997    1.988)    2.828   7.760   ( -28.628   10.340  -14.802)   33.847======================= Grid point 647 (133/145) =======================q-point: (-0.25  0.38  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1098Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.524   (   7.879    0.885   -0.000)    7.928   2.889   (   4.395   -8.777    0.000)    9.816   3.352   (  -5.605    2.618    0.000)    6.186   6.015   (  -1.379   -1.386    0.000)    1.955   6.016   (   0.774   -3.500    0.000)    3.585   8.104   ( -31.531    3.260    0.000)   31.700======================= Grid point 660 (134/145) =======================q-point: ( 0.56  0.44  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.512   (   0.000    9.647    2.261)    9.909   2.577   (  -0.000    8.057    4.494)    9.226   3.343   (  -0.000   -3.177   -4.414)    5.439   5.828   (   0.000    4.214    0.053)    4.214   5.957   (  -0.000    4.111    8.296)    9.259   6.842   (  -0.000   20.443  -26.292)   33.304======================= Grid point 661 (135/145) =======================q-point: ( 0.62  0.44  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.583   (   6.169    8.094    1.557)   10.295   2.796   ( -10.156    8.105    3.088)   13.356   3.245   (   0.684   -4.440   -3.988)    6.008   5.880   (  -0.244    4.482    0.035)    4.488   6.078   (  -0.192   -0.113    6.570)    6.574   6.889   ( -11.858   20.708  -16.992)   29.294======================= Grid point 662 (136/145) =======================q-point: (-0.31  0.44  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.603   (   5.990    5.818   -0.000)    8.351   3.035   ( -11.019    7.162    0.000)   13.142   3.137   (   1.973   -6.356    0.000)    6.655   5.938   (  -1.209    3.835   -0.000)    4.021   6.096   (   0.310   -2.359    0.000)    2.379   7.212   ( -22.071   17.225    0.000)   27.997======================= Grid point 679 (137/145) =======================q-point: (-0.38 -0.50  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 589Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.537   (   0.000    9.602    0.000)    9.602   2.631   (   0.000    9.315    0.000)    9.315   3.278   (  -0.000   -3.978    0.000)    3.978   5.829   (   0.000    4.247    0.000)    4.247   6.142   (  -0.000   -1.314    0.000)    1.314   6.443   (   0.000   24.096   -0.000)   24.096======================= Grid point 913 (138/145) =======================q-point: ( 0.56  0.38  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.519   (   0.000    1.652    9.113)    9.261   2.667   (  -0.000   10.261    4.632)   11.258   3.390   (  -0.000   -1.029   -3.212)    3.373   5.924   (  -0.000    4.358    0.610)    4.400   5.939   (   0.000    4.858    2.062)    5.278   7.981   (  -0.000   13.258  -33.449)   35.981======================= Grid point 914 (139/145) =======================q-point: ( 0.62  0.38  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.601   (   3.989    0.955    7.189)    8.277   2.847   (  -2.561    7.307    4.615)    9.014   3.321   (  -2.101   -2.127   -7.044)    7.652   5.967   (   1.586    3.161    0.881)    3.645   6.001   (  -0.741    2.076    0.779)    2.337   7.802   ( -13.501    6.651  -30.646)   34.142======================= Grid point 915 (140/145) =======================q-point: (-0.31  0.38  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1032Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.624   (   5.486    0.000    5.486)    7.758   2.947   (   2.445    0.000    2.445)    3.458   3.265   (  -7.073   -0.000   -7.073)   10.003   5.981   (   1.608    0.000    1.608)    2.274   6.021   (  -0.271   -0.000   -0.271)    0.384   7.742   ( -24.206    0.000  -24.206)   34.232======================= Grid point 931 (141/145) =======================q-point: ( 0.62  0.44  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1096Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.705   (   0.000    2.742    6.786)    7.319   2.745   (  -0.000   10.088    2.160)   10.317   3.264   (  -0.000   -3.633   -6.744)    7.660   5.931   (   0.000    4.478    0.127)    4.480   6.014   (  -0.000    0.969    4.690)    4.789   7.272   (  -0.000   14.402  -25.669)   29.433======================= Grid point 932 (142/145) =======================q-point: ( 0.69  0.44  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 2048Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.733   (   7.613    3.939    3.556)    9.280   2.947   (  -7.354    4.877    4.472)    9.893   3.134   (  -1.442   -4.532   -7.836)    9.166   5.979   (   0.567    3.792    0.124)    3.836   6.050   (   0.037   -2.202    2.620)    3.423   7.245   ( -13.135    8.491  -14.600)   21.396======================= Grid point 933 (143/145) =======================q-point: (-0.25  0.44  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 1088Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.705   (   7.891    2.538   -0.000)    8.289   2.963   (   2.059   -8.623    0.000)    8.865   3.170   (  -9.577    4.718    0.000)   10.676   6.004   (   0.868    1.741   -0.000)    1.946   6.027   (   0.212   -3.629    0.000)    3.636   7.434   ( -23.843    1.251    0.000)   23.875======================= Grid point 948 (144/145) =======================q-point: (-0.31 -0.50  0.19)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 580Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.769   (   0.000   10.039    0.000)   10.039   2.800   (   0.000    5.015    0.000)    5.015   3.160   (  -0.000   -6.320    0.000)    6.320   5.933   (   0.000    4.518    0.000)    4.518   6.099   (  -0.000   -2.414    0.000)    2.414   6.934   (   0.000   16.292    0.000)   16.292======================= Grid point 1216 (145/145) =======================q-point: (-0.25 -0.50  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.46e-03 1.56e-04 Number of triplets: 299Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.858   (   0.000   -0.000    0.000)    0.000   2.987   (  -0.000   -8.576    0.000)    8.576   2.987   (   0.000    8.576    0.000)    8.576   6.029   (  -0.000   -3.675    0.000)    3.675   6.029   (   0.000    3.675    0.000)    3.675   7.130   (   0.000   -0.000   -0.000)    0.000=================== End of collection of collisions ===================----------- Thermal conductivity (W/m-k) with tetrahedron method -----------#  T(K)        xx         yy         zz         yz         xz         xy        #ipm    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/24576   10.0    694.825    694.825    694.825     -0.000     -0.000      0.000 3/24576   20.0    160.618    160.618    160.618      0.000      0.000      0.000 3/24576   30.0     67.538     67.538     67.538      0.000      0.000      0.000 3/24576   40.0     43.292     43.292     43.292      0.000      0.000      0.000 3/24576   50.0     30.705     30.705     30.705      0.000      0.000      0.000 3/24576   60.0     22.862     22.862     22.862      0.000      0.000      0.000 3/24576   70.0     17.779     17.779     17.779      0.000      0.000      0.000 3/24576   80.0     14.370     14.370     14.370      0.000      0.000      0.000 3/24576   90.0     11.995     11.995     11.995      0.000      0.000      0.000 3/24576  100.0     10.274     10.274     10.274      0.000      0.000      0.000 3/24576  110.0      8.983      8.983      8.983      0.000      0.000      0.000 3/24576  120.0      7.985      7.985      7.985      0.000      0.000      0.000 3/24576  130.0      7.192      7.192      7.192      0.000      0.000      0.000 3/24576  140.0      6.548      6.548      6.548      0.000      0.000      0.000 3/24576  150.0      6.014      6.014      6.014      0.000      0.000      0.000 3/24576  160.0      5.565      5.565      5.565      0.000      0.000      0.000 3/24576  170.0      5.181      5.181      5.181      0.000      0.000      0.000 3/24576  180.0      4.850      4.850      4.850      0.000      0.000      0.000 3/24576  190.0      4.560      4.560      4.560      0.000      0.000      0.000 3/24576  200.0      4.305      4.305      4.305      0.000      0.000      0.000 3/24576  210.0      4.078      4.078      4.078      0.000      0.000      0.000 3/24576  220.0      3.874      3.874      3.874      0.000      0.000      0.000 3/24576  230.0      3.691      3.691      3.691      0.000      0.000      0.000 3/24576  240.0      3.525      3.525      3.525      0.000      0.000      0.000 3/24576  250.0      3.374      3.374      3.374      0.000      0.000      0.000 3/24576  260.0      3.236      3.236      3.236      0.000      0.000      0.000 3/24576  270.0      3.109      3.109      3.109      0.000      0.000      0.000 3/24576  280.0      2.992      2.992      2.992      0.000      0.000      0.000 3/24576  290.0      2.883      2.883      2.883      0.000      0.000      0.000 3/24576  300.0      2.783      2.783      2.783      0.000      0.000      0.000 3/24576  310.0      2.689      2.689      2.689      0.000      0.000      0.000 3/24576  320.0      2.602      2.602      2.602      0.000      0.000      0.000 3/24576  330.0      2.520      2.520      2.520      0.000      0.000      0.000 3/24576  340.0      2.443      2.443      2.443      0.000      0.000      0.000 3/24576  350.0      2.371      2.371      2.371      0.000      0.000      0.000 3/24576  360.0      2.304      2.304      2.304      0.000      0.000      0.000 3/24576  370.0      2.240      2.240      2.240      0.000      0.000      0.000 3/24576  380.0      2.179      2.179      2.179      0.000      0.000      0.000 3/24576  390.0      2.122      2.122      2.122      0.000      0.000      0.000 3/24576  400.0      2.068      2.068      2.068      0.000      0.000      0.000 3/24576  410.0      2.016      2.016      2.016      0.000      0.000      0.000 3/24576  420.0      1.967      1.967      1.967      0.000      0.000      0.000 3/24576  430.0      1.921      1.921      1.921      0.000      0.000      0.000 3/24576  440.0      1.876      1.876      1.876      0.000      0.000      0.000 3/24576  450.0      1.834      1.834      1.834      0.000      0.000      0.000 3/24576  460.0      1.794      1.794      1.794      0.000      0.000      0.000 3/24576  470.0      1.755      1.755      1.755      0.000      0.000      0.000 3/24576  480.0      1.718      1.718      1.718      0.000      0.000      0.000 3/24576  490.0      1.682      1.682      1.682      0.000      0.000      0.000 3/24576  500.0      1.648      1.648      1.648      0.000      0.000      0.000 3/24576  510.0      1.616      1.616      1.616      0.000      0.000      0.000 3/24576  520.0      1.584      1.584      1.584      0.000      0.000      0.000 3/24576  530.0      1.554      1.554      1.554      0.000      0.000      0.000 3/24576  540.0      1.525      1.525      1.525      0.000      0.000      0.000 3/24576  550.0      1.497      1.497      1.497      0.000      0.000      0.000 3/24576  560.0      1.470      1.470      1.470      0.000      0.000      0.000 3/24576  570.0      1.444      1.444      1.444      0.000      0.000      0.000 3/24576  580.0      1.419      1.419      1.419      0.000      0.000      0.000 3/24576  590.0      1.395      1.395      1.395      0.000      0.000      0.000 3/24576  600.0      1.371      1.371      1.371      0.000      0.000      0.000 3/24576  610.0      1.349      1.349      1.349      0.000      0.000      0.000 3/24576  620.0      1.327      1.327      1.327      0.000      0.000      0.000 3/24576  630.0      1.306      1.306      1.306      0.000      0.000      0.000 3/24576  640.0      1.285      1.285      1.285      0.000      0.000      0.000 3/24576  650.0      1.265      1.265      1.265      0.000      0.000      0.000 3/24576  660.0      1.246      1.246      1.246      0.000      0.000      0.000 3/24576  670.0      1.227      1.227      1.227      0.000      0.000      0.000 3/24576  680.0      1.209      1.209      1.209      0.000      0.000      0.000 3/24576  690.0      1.192      1.192      1.192      0.000      0.000      0.000 3/24576  700.0      1.175      1.175      1.175      0.000      0.000      0.000 3/24576  710.0      1.158      1.158      1.158      0.000      0.000      0.000 3/24576  720.0      1.142      1.142      1.142      0.000      0.000      0.000 3/24576  730.0      1.126      1.126      1.126      0.000      0.000      0.000 3/24576  740.0      1.111      1.111      1.111      0.000      0.000      0.000 3/24576  750.0      1.096      1.096      1.096      0.000      0.000      0.000 3/24576  760.0      1.082      1.082      1.082      0.000      0.000      0.000 3/24576  770.0      1.067      1.067      1.067      0.000      0.000      0.000 3/24576  780.0      1.054      1.054      1.054      0.000      0.000      0.000 3/24576  790.0      1.040      1.040      1.040      0.000      0.000      0.000 3/24576  800.0      1.027      1.027      1.027      0.000      0.000      0.000 3/24576  810.0      1.015      1.015      1.015      0.000      0.000      0.000 3/24576  820.0      1.002      1.002      1.002      0.000      0.000      0.000 3/24576  830.0      0.990      0.990      0.990      0.000      0.000      0.000 3/24576  840.0      0.978      0.978      0.978      0.000      0.000      0.000 3/24576  850.0      0.967      0.967      0.967      0.000      0.000      0.000 3/24576  860.0      0.956      0.956      0.956      0.000      0.000      0.000 3/24576  870.0      0.945      0.945      0.945      0.000      0.000      0.000 3/24576  880.0      0.934      0.934      0.934      0.000      0.000      0.000 3/24576  890.0      0.923      0.923      0.923      0.000      0.000      0.000 3/24576  900.0      0.913      0.913      0.913      0.000      0.000      0.000 3/24576  910.0      0.903      0.903      0.903      0.000      0.000      0.000 3/24576  920.0      0.893      0.893      0.893      0.000      0.000      0.000 3/24576  930.0      0.884      0.884      0.884      0.000      0.000      0.000 3/24576  940.0      0.874      0.874      0.874      0.000      0.000      0.000 3/24576  950.0      0.865      0.865      0.865      0.000      0.000      0.000 3/24576  960.0      0.856      0.856      0.856      0.000      0.000      0.000 3/24576  970.0      0.847      0.847      0.847      0.000      0.000      0.000 3/24576  980.0      0.838      0.838      0.838      0.000      0.000      0.000 3/24576  990.0      0.830      0.830      0.830      0.000      0.000      0.000 3/24576 1000.0      0.822      0.822      0.822      0.000      0.000      0.000 3/24576Thermal conductivity related properties were written into "kappa-m161616.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 05:22:43]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|