# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/d7bbf03f-2767-4b4f-90de-ad601204c16a/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 ZnO / F-43m (216) / materials id 1986](https://mdr.nims.go.jp/datasets/44c57ae3-f305-45e9-8fc0-dae7c5e91708)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-07 19:01:49]-------------------------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: F-43m (216)Number of symmetry operations in supercell: 768------------------------------ primitive cell ------------------------------Lattice vectors:  a    0.000000000000000    2.279496510000000    2.279496510000000  b    2.279496510000000    0.000000000000000    2.279496510000000  c    2.279496510000000    2.279496510000000    0.000000000000000Atomic positions (fractional):   *1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380   *2 O   0.75000000000000  0.75000000000000  0.75000000000000  15.999-------------------------------- unit cell ---------------------------------Lattice vectors:  a    4.558993020000000    0.000000000000000    0.000000000000000  b    0.000000000000000    4.558993020000000    0.000000000000000  c    0.000000000000000    0.000000000000000    4.558993020000000Atomic positions (fractional):   *1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380 > 1    2 Zn  0.00000000000000  0.50000000000000  0.50000000000000  65.380 > 1    3 Zn  0.50000000000000  0.00000000000000  0.50000000000000  65.380 > 1    4 Zn  0.50000000000000  0.50000000000000  0.00000000000000  65.380 > 1   *5 O   0.25000000000000  0.25000000000000  0.75000000000000  15.999 > 2    6 O   0.25000000000000  0.75000000000000  0.25000000000000  15.999 > 2    7 O   0.75000000000000  0.25000000000000  0.25000000000000  15.999 > 2    8 O   0.75000000000000  0.75000000000000  0.75000000000000  15.999 > 2-------------------------------- super cell --------------------------------Lattice vectors:  a    9.117986040000000    0.000000000000000    0.000000000000000  b    0.000000000000000    9.117986040000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.117986040000000Atomic positions (fractional):   *1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380 > 1    2 Zn  0.50000000000000  0.00000000000000  0.00000000000000  65.380 > 1    3 Zn  0.00000000000000  0.50000000000000  0.00000000000000  65.380 > 1    4 Zn  0.50000000000000  0.50000000000000  0.00000000000000  65.380 > 1    5 Zn  0.00000000000000  0.00000000000000  0.50000000000000  65.380 > 1    6 Zn  0.50000000000000  0.00000000000000  0.50000000000000  65.380 > 1    7 Zn  0.00000000000000  0.50000000000000  0.50000000000000  65.380 > 1    8 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 1    9 Zn  0.00000000000000  0.25000000000000  0.25000000000000  65.380 > 1   10 Zn  0.50000000000000  0.25000000000000  0.25000000000000  65.380 > 1   11 Zn  0.00000000000000  0.75000000000000  0.25000000000000  65.380 > 1   12 Zn  0.50000000000000  0.75000000000000  0.25000000000000  65.380 > 1   13 Zn  0.00000000000000  0.25000000000000  0.75000000000000  65.380 > 1   14 Zn  0.50000000000000  0.25000000000000  0.75000000000000  65.380 > 1   15 Zn  0.00000000000000  0.75000000000000  0.75000000000000  65.380 > 1   16 Zn  0.50000000000000  0.75000000000000  0.75000000000000  65.380 > 1   17 Zn  0.25000000000000  0.00000000000000  0.25000000000000  65.380 > 1   18 Zn  0.75000000000000  0.00000000000000  0.25000000000000  65.380 > 1   19 Zn  0.25000000000000  0.50000000000000  0.25000000000000  65.380 > 1   20 Zn  0.75000000000000  0.50000000000000  0.25000000000000  65.380 > 1   21 Zn  0.25000000000000  0.00000000000000  0.75000000000000  65.380 > 1   22 Zn  0.75000000000000  0.00000000000000  0.75000000000000  65.380 > 1   23 Zn  0.25000000000000  0.50000000000000  0.75000000000000  65.380 > 1   24 Zn  0.75000000000000  0.50000000000000  0.75000000000000  65.380 > 1   25 Zn  0.25000000000000  0.25000000000000  0.00000000000000  65.380 > 1   26 Zn  0.75000000000000  0.25000000000000  0.00000000000000  65.380 > 1   27 Zn  0.25000000000000  0.75000000000000  0.00000000000000  65.380 > 1   28 Zn  0.75000000000000  0.75000000000000  0.00000000000000  65.380 > 1   29 Zn  0.25000000000000  0.25000000000000  0.50000000000000  65.380 > 1   30 Zn  0.75000000000000  0.25000000000000  0.50000000000000  65.380 > 1   31 Zn  0.25000000000000  0.75000000000000  0.50000000000000  65.380 > 1   32 Zn  0.75000000000000  0.75000000000000  0.50000000000000  65.380 > 1  *33 O   0.12500000000000  0.12500000000000  0.37500000000000  15.999 > 2   34 O   0.62500000000000  0.12500000000000  0.37500000000000  15.999 > 2   35 O   0.12500000000000  0.62500000000000  0.37500000000000  15.999 > 2   36 O   0.62500000000000  0.62500000000000  0.37500000000000  15.999 > 2   37 O   0.12500000000000  0.12500000000000  0.87500000000000  15.999 > 2   38 O   0.62500000000000  0.12500000000000  0.87500000000000  15.999 > 2   39 O   0.12500000000000  0.62500000000000  0.87500000000000  15.999 > 2   40 O   0.62500000000000  0.62500000000000  0.87500000000000  15.999 > 2   41 O   0.12500000000000  0.37500000000000  0.12500000000000  15.999 > 2   42 O   0.62500000000000  0.37500000000000  0.12500000000000  15.999 > 2   43 O   0.12500000000000  0.87500000000000  0.12500000000000  15.999 > 2   44 O   0.62500000000000  0.87500000000000  0.12500000000000  15.999 > 2   45 O   0.12500000000000  0.37500000000000  0.62500000000000  15.999 > 2   46 O   0.62500000000000  0.37500000000000  0.62500000000000  15.999 > 2   47 O   0.12500000000000  0.87500000000000  0.62500000000000  15.999 > 2   48 O   0.62500000000000  0.87500000000000  0.62500000000000  15.999 > 2   49 O   0.37500000000000  0.12500000000000  0.12500000000000  15.999 > 2   50 O   0.87500000000000  0.12500000000000  0.12500000000000  15.999 > 2   51 O   0.37500000000000  0.62500000000000  0.12500000000000  15.999 > 2   52 O   0.87500000000000  0.62500000000000  0.12500000000000  15.999 > 2   53 O   0.37500000000000  0.12500000000000  0.62500000000000  15.999 > 2   54 O   0.87500000000000  0.12500000000000  0.62500000000000  15.999 > 2   55 O   0.37500000000000  0.62500000000000  0.62500000000000  15.999 > 2   56 O   0.87500000000000  0.62500000000000  0.62500000000000  15.999 > 2   57 O   0.37500000000000  0.37500000000000  0.37500000000000  15.999 > 2   58 O   0.87500000000000  0.37500000000000  0.37500000000000  15.999 > 2   59 O   0.37500000000000  0.87500000000000  0.37500000000000  15.999 > 2   60 O   0.87500000000000  0.87500000000000  0.37500000000000  15.999 > 2   61 O   0.37500000000000  0.37500000000000  0.87500000000000  15.999 > 2   62 O   0.87500000000000  0.37500000000000  0.87500000000000  15.999 > 2   63 O   0.37500000000000  0.87500000000000  0.87500000000000  15.999 > 2   64 O   0.87500000000000  0.87500000000000  0.87500000000000  15.999 > 2----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            6.4261854    0.0000000    0.0000000            0.0000000    6.4261854    0.0000000            0.0000000    0.0000000    6.4261854-------------------------- Born effective charges --------------------------    1 Zn    2.1536386    0.0000000    0.0000000            0.0000000    2.1536386    0.0000000            0.0000000    0.0000000    2.1536386    2 O    -2.1536386    0.0000000    0.0000000            0.0000000   -2.1536386    0.0000000            0.0000000    0.0000000   -2.1536386----------------------------------------------------------------------------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: 42Number of blocks in projector: 42Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 2--- Eigsh_solver_block: 1 / 2 ---Block_size: 26Use standard eigh solver.--- Eigsh_solver_block: 2 / 2 ---Block_size: 16Use standard eigh solver.Tree of FC basis block matrices:- (42, 40), data: False|-- (16, 16), data: True|-- (26, 24), data: True-----Solver_atoms: 1 -- 64 / 64Time (Solver_compr_matrix_reshape): 0.002Solver_block: 80 / 80 - Time: 0.009Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.011--------------------------------- Symfc end --------------------------------Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 192/192Permutation basis: 852/852Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 42Number of blocks in projector: 42Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 2--- Eigsh_solver_block: 1 / 2 ---Block_size: 26Use standard eigh solver.--- Eigsh_solver_block: 2 / 2 ---Block_size: 16Use standard eigh solver.Tree of FC basis block matrices:- (42, 40), data: False|-- (16, 16), data: True|-- (26, 24), data: TrueMax drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-07 19:01:51]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-07 19:01:52]-------------------------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: F-43m (216)------------------------------ primitive cell ------------------------------Lattice vectors:  a    0.000000000000000    2.279496510000000    2.279496510000000  b    2.279496510000000    0.000000000000000    2.279496510000000  c    2.279496510000000    2.279496510000000    0.000000000000000Atomic positions (fractional):    1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380    2 O   0.75000000000000  0.75000000000000  0.75000000000000  15.999-------------------------------- supercell ---------------------------------Lattice vectors:  a    9.117986040000000    0.000000000000000    0.000000000000000  b    0.000000000000000    9.117986040000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.117986040000000Atomic positions (fractional):    1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380 > 1    2 Zn  0.50000000000000  0.00000000000000  0.00000000000000  65.380 > 1    3 Zn  0.00000000000000  0.50000000000000  0.00000000000000  65.380 > 1    4 Zn  0.50000000000000  0.50000000000000  0.00000000000000  65.380 > 1    5 Zn  0.00000000000000  0.00000000000000  0.50000000000000  65.380 > 1    6 Zn  0.50000000000000  0.00000000000000  0.50000000000000  65.380 > 1    7 Zn  0.00000000000000  0.50000000000000  0.50000000000000  65.380 > 1    8 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 1    9 Zn  0.00000000000000  0.25000000000000  0.25000000000000  65.380 > 1   10 Zn  0.50000000000000  0.25000000000000  0.25000000000000  65.380 > 1   11 Zn  0.00000000000000  0.75000000000000  0.25000000000000  65.380 > 1   12 Zn  0.50000000000000  0.75000000000000  0.25000000000000  65.380 > 1   13 Zn  0.00000000000000  0.25000000000000  0.75000000000000  65.380 > 1   14 Zn  0.50000000000000  0.25000000000000  0.75000000000000  65.380 > 1   15 Zn  0.00000000000000  0.75000000000000  0.75000000000000  65.380 > 1   16 Zn  0.50000000000000  0.75000000000000  0.75000000000000  65.380 > 1   17 Zn  0.25000000000000  0.00000000000000  0.25000000000000  65.380 > 1   18 Zn  0.75000000000000  0.00000000000000  0.25000000000000  65.380 > 1   19 Zn  0.25000000000000  0.50000000000000  0.25000000000000  65.380 > 1   20 Zn  0.75000000000000  0.50000000000000  0.25000000000000  65.380 > 1   21 Zn  0.25000000000000  0.00000000000000  0.75000000000000  65.380 > 1   22 Zn  0.75000000000000  0.00000000000000  0.75000000000000  65.380 > 1   23 Zn  0.25000000000000  0.50000000000000  0.75000000000000  65.380 > 1   24 Zn  0.75000000000000  0.50000000000000  0.75000000000000  65.380 > 1   25 Zn  0.25000000000000  0.25000000000000  0.00000000000000  65.380 > 1   26 Zn  0.75000000000000  0.25000000000000  0.00000000000000  65.380 > 1   27 Zn  0.25000000000000  0.75000000000000  0.00000000000000  65.380 > 1   28 Zn  0.75000000000000  0.75000000000000  0.00000000000000  65.380 > 1   29 Zn  0.25000000000000  0.25000000000000  0.50000000000000  65.380 > 1   30 Zn  0.75000000000000  0.25000000000000  0.50000000000000  65.380 > 1   31 Zn  0.25000000000000  0.75000000000000  0.50000000000000  65.380 > 1   32 Zn  0.75000000000000  0.75000000000000  0.50000000000000  65.380 > 1   33 O   0.12500000000000  0.12500000000000  0.37500000000000  15.999 > 33   34 O   0.62500000000000  0.12500000000000  0.37500000000000  15.999 > 33   35 O   0.12500000000000  0.62500000000000  0.37500000000000  15.999 > 33   36 O   0.62500000000000  0.62500000000000  0.37500000000000  15.999 > 33   37 O   0.12500000000000  0.12500000000000  0.87500000000000  15.999 > 33   38 O   0.62500000000000  0.12500000000000  0.87500000000000  15.999 > 33   39 O   0.12500000000000  0.62500000000000  0.87500000000000  15.999 > 33   40 O   0.62500000000000  0.62500000000000  0.87500000000000  15.999 > 33   41 O   0.12500000000000  0.37500000000000  0.12500000000000  15.999 > 33   42 O   0.62500000000000  0.37500000000000  0.12500000000000  15.999 > 33   43 O   0.12500000000000  0.87500000000000  0.12500000000000  15.999 > 33   44 O   0.62500000000000  0.87500000000000  0.12500000000000  15.999 > 33   45 O   0.12500000000000  0.37500000000000  0.62500000000000  15.999 > 33   46 O   0.62500000000000  0.37500000000000  0.62500000000000  15.999 > 33   47 O   0.12500000000000  0.87500000000000  0.62500000000000  15.999 > 33   48 O   0.62500000000000  0.87500000000000  0.62500000000000  15.999 > 33   49 O   0.37500000000000  0.12500000000000  0.12500000000000  15.999 > 33   50 O   0.87500000000000  0.12500000000000  0.12500000000000  15.999 > 33   51 O   0.37500000000000  0.62500000000000  0.12500000000000  15.999 > 33   52 O   0.87500000000000  0.62500000000000  0.12500000000000  15.999 > 33   53 O   0.37500000000000  0.12500000000000  0.62500000000000  15.999 > 33   54 O   0.87500000000000  0.12500000000000  0.62500000000000  15.999 > 33   55 O   0.37500000000000  0.62500000000000  0.62500000000000  15.999 > 33   56 O   0.87500000000000  0.62500000000000  0.62500000000000  15.999 > 33   57 O   0.37500000000000  0.37500000000000  0.37500000000000  15.999 > 33   58 O   0.87500000000000  0.37500000000000  0.37500000000000  15.999 > 33   59 O   0.37500000000000  0.87500000000000  0.37500000000000  15.999 > 33   60 O   0.87500000000000  0.87500000000000  0.37500000000000  15.999 > 33   61 O   0.37500000000000  0.37500000000000  0.87500000000000  15.999 > 33   62 O   0.87500000000000  0.37500000000000  0.87500000000000  15.999 > 33   63 O   0.37500000000000  0.87500000000000  0.87500000000000  15.999 > 33   64 O   0.87500000000000  0.87500000000000  0.87500000000000  15.999 > 33----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            6.4261854    0.0000000    0.0000000            0.0000000    6.4261854    0.0000000            0.0000000    0.0000000    6.4261854-------------------------- Born effective charges --------------------------    1 Zn    2.1536386    0.0000000    0.0000000            0.0000000    2.1536386    0.0000000            0.0000000    0.0000000    2.1536386    2 O    -2.1536386    0.0000000    0.0000000            0.0000000   -2.1536386    0.0000000            0.0000000    0.0000000   -2.1536386----------------------------------------------------------------------------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.00000025 (xyz) 0.00000025 (xyz) 0.00000025 (xzy)fc3 was written into "fc3.hdf5".Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330----------- None of ph-ph interaction calculation was performed. -----------Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-07 19:01:52]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-07 19:01:53]-------------------------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: F-43m (216)------------------------------ primitive cell ------------------------------Lattice vectors:  a    0.000000000000000    2.279496510000000    2.279496510000000  b    2.279496510000000    0.000000000000000    2.279496510000000  c    2.279496510000000    2.279496510000000    0.000000000000000Atomic positions (fractional):    1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380    2 O   0.75000000000000  0.75000000000000  0.75000000000000  15.999-------------------------------- supercell ---------------------------------Lattice vectors:  a    9.117986040000000    0.000000000000000    0.000000000000000  b    0.000000000000000    9.117986040000000    0.000000000000000  c    0.000000000000000    0.000000000000000    9.117986040000000Atomic positions (fractional):    1 Zn  0.00000000000000  0.00000000000000  0.00000000000000  65.380 > 1    2 Zn  0.50000000000000  0.00000000000000  0.00000000000000  65.380 > 1    3 Zn  0.00000000000000  0.50000000000000  0.00000000000000  65.380 > 1    4 Zn  0.50000000000000  0.50000000000000  0.00000000000000  65.380 > 1    5 Zn  0.00000000000000  0.00000000000000  0.50000000000000  65.380 > 1    6 Zn  0.50000000000000  0.00000000000000  0.50000000000000  65.380 > 1    7 Zn  0.00000000000000  0.50000000000000  0.50000000000000  65.380 > 1    8 Zn  0.50000000000000  0.50000000000000  0.50000000000000  65.380 > 1    9 Zn  0.00000000000000  0.25000000000000  0.25000000000000  65.380 > 1   10 Zn  0.50000000000000  0.25000000000000  0.25000000000000  65.380 > 1   11 Zn  0.00000000000000  0.75000000000000  0.25000000000000  65.380 > 1   12 Zn  0.50000000000000  0.75000000000000  0.25000000000000  65.380 > 1   13 Zn  0.00000000000000  0.25000000000000  0.75000000000000  65.380 > 1   14 Zn  0.50000000000000  0.25000000000000  0.75000000000000  65.380 > 1   15 Zn  0.00000000000000  0.75000000000000  0.75000000000000  65.380 > 1   16 Zn  0.50000000000000  0.75000000000000  0.75000000000000  65.380 > 1   17 Zn  0.25000000000000  0.00000000000000  0.25000000000000  65.380 > 1   18 Zn  0.75000000000000  0.00000000000000  0.25000000000000  65.380 > 1   19 Zn  0.25000000000000  0.50000000000000  0.25000000000000  65.380 > 1   20 Zn  0.75000000000000  0.50000000000000  0.25000000000000  65.380 > 1   21 Zn  0.25000000000000  0.00000000000000  0.75000000000000  65.380 > 1   22 Zn  0.75000000000000  0.00000000000000  0.75000000000000  65.380 > 1   23 Zn  0.25000000000000  0.50000000000000  0.75000000000000  65.380 > 1   24 Zn  0.75000000000000  0.50000000000000  0.75000000000000  65.380 > 1   25 Zn  0.25000000000000  0.25000000000000  0.00000000000000  65.380 > 1   26 Zn  0.75000000000000  0.25000000000000  0.00000000000000  65.380 > 1   27 Zn  0.25000000000000  0.75000000000000  0.00000000000000  65.380 > 1   28 Zn  0.75000000000000  0.75000000000000  0.00000000000000  65.380 > 1   29 Zn  0.25000000000000  0.25000000000000  0.50000000000000  65.380 > 1   30 Zn  0.75000000000000  0.25000000000000  0.50000000000000  65.380 > 1   31 Zn  0.25000000000000  0.75000000000000  0.50000000000000  65.380 > 1   32 Zn  0.75000000000000  0.75000000000000  0.50000000000000  65.380 > 1   33 O   0.12500000000000  0.12500000000000  0.37500000000000  15.999 > 33   34 O   0.62500000000000  0.12500000000000  0.37500000000000  15.999 > 33   35 O   0.12500000000000  0.62500000000000  0.37500000000000  15.999 > 33   36 O   0.62500000000000  0.62500000000000  0.37500000000000  15.999 > 33   37 O   0.12500000000000  0.12500000000000  0.87500000000000  15.999 > 33   38 O   0.62500000000000  0.12500000000000  0.87500000000000  15.999 > 33   39 O   0.12500000000000  0.62500000000000  0.87500000000000  15.999 > 33   40 O   0.62500000000000  0.62500000000000  0.87500000000000  15.999 > 33   41 O   0.12500000000000  0.37500000000000  0.12500000000000  15.999 > 33   42 O   0.62500000000000  0.37500000000000  0.12500000000000  15.999 > 33   43 O   0.12500000000000  0.87500000000000  0.12500000000000  15.999 > 33   44 O   0.62500000000000  0.87500000000000  0.12500000000000  15.999 > 33   45 O   0.12500000000000  0.37500000000000  0.62500000000000  15.999 > 33   46 O   0.62500000000000  0.37500000000000  0.62500000000000  15.999 > 33   47 O   0.12500000000000  0.87500000000000  0.62500000000000  15.999 > 33   48 O   0.62500000000000  0.87500000000000  0.62500000000000  15.999 > 33   49 O   0.37500000000000  0.12500000000000  0.12500000000000  15.999 > 33   50 O   0.87500000000000  0.12500000000000  0.12500000000000  15.999 > 33   51 O   0.37500000000000  0.62500000000000  0.12500000000000  15.999 > 33   52 O   0.87500000000000  0.62500000000000  0.12500000000000  15.999 > 33   53 O   0.37500000000000  0.12500000000000  0.62500000000000  15.999 > 33   54 O   0.87500000000000  0.12500000000000  0.62500000000000  15.999 > 33   55 O   0.37500000000000  0.62500000000000  0.62500000000000  15.999 > 33   56 O   0.87500000000000  0.62500000000000  0.62500000000000  15.999 > 33   57 O   0.37500000000000  0.37500000000000  0.37500000000000  15.999 > 33   58 O   0.87500000000000  0.37500000000000  0.37500000000000  15.999 > 33   59 O   0.37500000000000  0.87500000000000  0.37500000000000  15.999 > 33   60 O   0.87500000000000  0.87500000000000  0.37500000000000  15.999 > 33   61 O   0.37500000000000  0.37500000000000  0.87500000000000  15.999 > 33   62 O   0.87500000000000  0.37500000000000  0.87500000000000  15.999 > 33   63 O   0.37500000000000  0.87500000000000  0.87500000000000  15.999 > 33   64 O   0.87500000000000  0.87500000000000  0.87500000000000  15.999 > 33----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            6.4261854    0.0000000    0.0000000            0.0000000    6.4261854    0.0000000            0.0000000    0.0000000    6.4261854-------------------------- Born effective charges --------------------------    1 Zn    2.1536386    0.0000000    0.0000000            0.0000000    2.1536386    0.0000000            0.0000000    0.0000000    2.1536386    2 O    -2.1536386    0.0000000    0.0000000            0.0000000   -2.1536386    0.0000000            0.0000000    0.0000000   -2.1536386----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: 0.00000025 (xyz) 0.00000025 (xyz) 0.00000025 (xzy)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 ... [ 19 19 19 ]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.45, Number of G-points: 307, Lambda: 0.38Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/220) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 220Calculating 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  11.587   (   0.000    0.000    0.000)    0.000  11.587   (   0.000    0.000    0.000)    0.000  11.587   (  -0.000    0.000    0.000)    0.000======================= Grid point 1 (2/220) =======================q-point: ( 0.05  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.579   ( -16.253   16.253   16.253)   28.150   0.579   ( -16.253   16.253   16.253)   28.150   1.310   ( -37.301   37.301   37.301)   64.608  11.646   (  -3.353    3.353    3.353)    5.807  11.646   (  -3.353    3.353    3.353)    5.807  15.512   (  -2.806    2.806    2.806)    4.860======================= Grid point 2 (3/220) =======================q-point: ( 0.11  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.111   ( -14.314   14.314   14.314)   24.792   1.111   ( -14.314   14.314   14.314)   24.792   2.568   ( -35.144   35.144   35.144)   60.872  11.809   (  -5.820    5.820    5.820)   10.081  11.809   (  -5.820    5.820    5.820)   10.081  15.647   (  -4.798    4.798    4.798)    8.310======================= Grid point 3 (4/220) =======================q-point: ( 0.16  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.562   ( -11.599   11.599   11.599)   20.090   1.562   ( -11.599   11.599   11.599)   20.090   3.732   ( -31.976   31.976   31.976)   55.383  12.033   (  -6.881    6.881    6.881)   11.918  12.033   (  -6.881    6.881    6.881)   11.918  15.829   (  -5.463    5.463    5.463)    9.463======================= Grid point 4 (5/220) =======================q-point: ( 0.21  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.912   (  -8.654    8.654    8.654)   14.989   1.912   (  -8.654    8.654    8.654)   14.989   4.776   ( -28.203   28.203   28.203)   48.849  12.269   (  -6.538    6.538    6.538)   11.323  12.269   (  -6.538    6.538    6.538)   11.323  16.010   (  -4.765    4.765    4.765)    8.253======================= Grid point 5 (6/220) =======================q-point: ( 0.26  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.164   (  -5.944    5.944    5.944)   10.294   2.164   (  -5.944    5.944    5.944)   10.294   5.682   ( -24.086   24.086   24.086)   41.718  12.475   (  -5.217    5.217    5.217)    9.036  12.475   (  -5.217    5.217    5.217)    9.036  16.148   (  -3.119    3.119    3.119)    5.403======================= Grid point 6 (7/220) =======================q-point: ( 0.32  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.330   (  -3.752    3.752    3.752)    6.499   2.330   (  -3.752    3.752    3.752)    6.500   6.440   ( -19.658   19.658   19.658)   34.049  12.626   (  -3.524    3.524    3.524)    6.104  12.626   (  -3.524    3.524    3.524)    6.104  16.223   (  -1.220    1.220    1.220)    2.113======================= Grid point 7 (8/220) =======================q-point: ( 0.37  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.431   (  -2.161    2.161    2.161)    3.743   2.431   (  -2.161    2.161    2.161)    3.743   7.038   ( -14.773   14.773   14.773)   25.588  12.721   (  -1.991    1.991    1.991)    3.449  12.721   (  -1.991    1.991    1.991)    3.449  16.238   (   0.221   -0.221   -0.221)    0.383======================= Grid point 8 (9/220) =======================q-point: ( 0.42  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.486   (  -1.080    1.080    1.080)    1.871   2.486   (  -1.080    1.080    1.080)    1.871   7.457   (  -9.267    9.267    9.267)   16.051  12.769   (  -0.902    0.902    0.902)    1.563  12.769   (  -0.902    0.902    0.902)    1.563  16.218   (   0.748   -0.748   -0.748)    1.295======================= Grid point 9 (10/220) =======================q-point: ( 0.47  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 670Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.509   (  -0.324    0.324    0.324)    0.561   2.509   (  -0.324    0.324    0.324)    0.561   7.673   (  -3.172    3.172    3.172)    5.494  12.788   (  -0.243    0.243    0.243)    0.421  12.788   (  -0.243    0.243    0.243)    0.421  16.197   (   0.375   -0.375   -0.375)    0.649======================= Grid point 20 (11/220) =======================q-point: ( 0.05  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.758   (  -0.000    0.000   32.039)   32.039   0.758   (  -0.000    0.000   32.039)   32.039   1.310   (  -0.000    0.000   55.938)   55.938  11.672   (  -0.000    0.000    7.289)    7.289  11.672   (  -0.000    0.000    7.289)    7.289  15.531   (  -0.000    0.000    5.732)    5.732======================= Grid point 21 (12/220) =======================q-point: ( 0.11  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.187   (  -5.534    5.534   28.035)   29.107   1.211   (  -3.876    3.876   30.086)   30.581   2.258   ( -25.757   25.757   44.996)   57.892  11.802   (  -2.373    2.373   10.173)   10.712  11.809   (  -2.839    2.839   10.271)   11.028  15.635   (  -2.235    2.235    7.535)    8.171======================= Grid point 22 (13/220) =======================q-point: ( 0.16  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.611   (  -2.295    2.295   27.084)   27.277   1.620   (  -5.892    5.892   23.197)   24.648   3.385   ( -29.988   29.988   36.020)   55.642  12.003   (  -3.829    3.829   11.553)   12.760  12.026   (  -4.460    4.460   12.046)   13.598  15.795   (  -3.698    3.698    7.648)    9.265======================= Grid point 23 (14/220) =======================q-point: ( 0.21  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.939   (  -1.274    1.274   22.559)   22.631   1.983   (  -4.326    4.326   18.870)   19.837   4.445   ( -28.540   28.540   29.907)   50.234  12.230   (  -3.954    3.954   11.457)   12.749  12.276   (  -4.526    4.526   12.455)   14.003  15.968   (  -3.996    3.996    6.201)    8.390======================= Grid point 24 (15/220) =======================q-point: ( 0.26  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.190   (   0.054   -0.054   18.553)   18.553   2.254   (  -2.149    2.149   15.164)   15.466   5.385   ( -25.487   25.487   24.619)   43.649  12.439   (  -3.034    3.034   10.396)   11.247  12.508   (  -3.269    3.269   11.673)   12.555  16.111   (  -3.225    3.225    3.775)    5.921======================= Grid point 25 (16/220) =======================q-point: ( 0.32  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.370   (   1.446   -1.446   15.780)   15.913   2.435   (  -0.047    0.047   12.075)   12.075   6.186   ( -21.706   21.706   19.552)   36.394  12.605   (  -1.580    1.580    9.139)    9.408  12.688   (  -1.308    1.308   10.063)   10.232  16.198   (  -1.860    1.860    1.180)    2.883======================= Grid point 26 (17/220) =======================q-point: ( 0.37  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.494   (   2.710   -2.710   14.382)   14.884   2.540   (   1.615   -1.615    9.470)    9.742   6.833   ( -17.328   17.328   14.353)   28.400  12.724   (  -0.047    0.047    8.385)    8.386  12.799   (   0.635   -0.635    8.004)    8.054  16.227   (  -0.546    0.546   -0.821)    1.127======================= Grid point 27 (18/220) =======================q-point: ( 0.42  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.582   (   3.846   -3.846   14.158)   15.167   2.585   (   2.689   -2.689    7.133)    8.084   7.310   ( -12.263   12.263    8.748)   19.424  12.806   (   1.327   -1.327    8.500)    8.704  12.847   (   1.998   -1.998    5.783)    6.437  16.214   (   0.155   -0.155   -1.768)    1.782======================= Grid point 28 (19/220) =======================q-point: ( 0.47  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.585   (   3.179   -3.179    4.853)    6.615   2.647   (   4.973   -4.973   14.700)   16.295   7.594   (  -6.490    6.490    2.660)    9.556  12.847   (   2.550   -2.550    3.572)    5.076  12.863   (   2.546   -2.546    9.376)   10.044  16.189   (   0.017   -0.017   -1.688)    1.688======================= Grid point 29 (20/220) =======================q-point: (-0.47  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.553   (   3.202   -3.202    2.519)    5.183   2.692   (   6.225   -6.225   15.538)   17.859   7.667   (  -0.252    0.252   -3.666)    3.683  12.817   (   2.437   -2.437    1.451)    3.740  12.905   (   3.789   -3.789   10.484)   11.774  16.180   (  -0.694    0.694   -1.088)    1.465======================= Grid point 30 (21/220) =======================q-point: (-0.42  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.497   (   2.994   -2.994    0.141)    4.236   2.716   (   7.712   -7.712   16.279)   19.595   7.522   (   5.979   -5.979   -9.752)   12.908  12.770   (   2.098   -2.098   -0.543)    3.016  12.927   (   5.251   -5.251   11.090)   13.347  16.195   (  -1.329    1.329   -0.691)    2.002======================= Grid point 31 (22/220) =======================q-point: (-0.37  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.418   (   2.905   -2.905   -2.208)    4.664   2.708   (   9.533   -9.533   16.666)   21.436   7.171   (  11.781  -11.781  -15.191)   22.547  12.706   (   2.057   -2.057   -2.363)    3.749  12.912   (   7.005   -7.005   10.572)   14.489  16.217   (  -1.234    1.234   -0.992)    2.008======================= Grid point 32 (23/220) =======================q-point: (-0.32  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.308   (   3.347   -3.347   -4.456)    6.501   2.655   (  11.798  -11.798   16.615)   23.546   6.635   (  17.010  -17.010  -19.854)   31.190  12.616   (   2.686   -2.686   -3.909)    5.450  12.841   (   8.920   -8.920    8.708)   15.329  16.218   (  -0.128    0.128   -1.982)    1.991======================= Grid point 33 (24/220) =======================q-point: (-0.26  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.154   (   4.661   -4.661   -6.522)    9.273   2.541   (  14.609  -14.609   16.247)   26.283   5.934   (  21.787  -21.787  -23.807)   38.937  12.488   (   4.003   -4.003   -4.995)    7.550  12.699   (  10.658  -10.658    5.804)   16.151  16.170   (   1.779   -1.779   -3.203)    4.073======================= Grid point 34 (25/220) =======================q-point: (-0.21  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.936   (   6.976   -6.976   -8.204)   12.831   2.351   (  17.973  -17.973   15.933)   29.999   5.084   (  26.325  -26.325  -27.032)   46.008  12.316   (   5.629   -5.629   -5.359)    9.596  12.487   (  11.740  -11.740    2.628)   16.810  16.062   (   3.905   -3.905   -4.002)    6.820======================= Grid point 35 (26/220) =======================q-point: (-0.16  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.639   (  10.153  -10.153   -8.965)   16.928   2.079   (  21.680  -21.680   16.164)   34.660   4.099   (  30.796  -30.796  -28.947)   52.295  12.111   (   6.891   -6.891   -4.725)   10.831  12.229   (  11.682  -11.682    0.155)   16.522  15.905   (   5.513   -5.513   -3.816)    8.680======================= Grid point 36 (27/220) =======================q-point: (-0.11  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.265   (  13.811  -13.811   -7.471)   20.912   1.727   (  24.955  -24.955   16.111)   38.795   3.007   (  35.438  -35.438  -26.590)   56.734  11.902   (   7.090   -7.090   -2.915)   10.442  11.968   (  10.174  -10.174   -0.813)   14.411  15.733   (   5.979   -5.979   -2.402)    8.790======================= Grid point 37 (28/220) =======================q-point: (-0.05  0.05  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.854   (  17.022  -17.022   -0.000)   24.073   1.262   (  26.234  -26.234   -0.000)   37.100   1.929   (  40.213  -40.213   -0.000)   56.869  11.732   (   5.818   -5.818   -0.000)    8.228  11.758   (   7.287   -7.287   -0.000)   10.305  15.589   (   5.024   -5.024   -0.000)    7.105======================= Grid point 40 (29/220) =======================q-point: ( 0.11  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.464   (  -0.000    0.000   28.798)   28.798   1.464   (  -0.000    0.000   28.798)   28.798   2.567   (  -0.000    0.000   52.609)   52.609  11.914   (  -0.000    0.000   13.368)   13.368  11.914   (  -0.000    0.000   13.368)   13.368  15.711   (  -0.000    0.000    9.462)    9.462======================= Grid point 41 (30/220) =======================q-point: ( 0.16  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.812   (  -2.349    2.349   25.557)   25.772   1.859   (  -4.504    4.504   26.006)   26.774   3.324   ( -14.428   14.428   46.154)   50.463  12.101   (  -1.576    1.576   15.321)   15.482  12.108   (  -2.167    2.167   15.254)   15.559  15.840   (  -1.276    1.276    9.625)    9.793======================= Grid point 42 (31/220) =======================q-point: ( 0.21  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.141   (  -2.357    2.357   21.423)   21.681   2.220   (  -1.173    1.173   24.757)   24.812   4.229   ( -21.323   21.323   36.736)   47.527  12.327   (  -1.968    1.968   16.044)   16.284  12.351   (  -3.041    3.041   15.697)   16.276  15.980   (  -1.898    1.898    7.927)    8.369======================= Grid point 43 (32/220) =======================q-point: ( 0.26  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.407   (  -1.117    1.117   17.541)   17.612   2.488   (   1.999   -1.999   23.712)   23.880   5.120   ( -22.637   22.637   28.570)   42.908  12.549   (  -1.162    1.162   15.886)   15.971  12.594   (  -2.492    2.492   14.839)   15.251  16.100   (  -1.860    1.860    4.966)    5.620======================= Grid point 44 (33/220) =======================q-point: ( 0.32  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.596   (   0.609   -0.609   14.265)   14.291   2.683   (   4.260   -4.260   22.804)   23.587   5.918   ( -21.155   21.155   21.696)   36.956  12.741   (   0.343   -0.343   15.462)   15.470  12.796   (  -0.864    0.864   13.053)   13.110  16.176   (  -1.309    1.309    1.700)    2.514======================= Grid point 45 (34/220) =======================q-point: ( 0.37  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.711   (   2.348   -2.348   11.598)   12.064   2.826   (   5.847   -5.847   22.434)   23.909   6.588   ( -18.162   18.162   15.446)   29.972  12.891   (   1.982   -1.982   15.369)   15.622  12.930   (   1.227   -1.227   10.763)   10.901  16.201   (  -0.605    0.605   -0.983)    1.303======================= Grid point 46 (35/220) =======================q-point: ( 0.42  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.759   (   3.835   -3.835    9.356)   10.814   2.936   (   6.985   -6.985   22.626)   24.689   7.105   ( -14.066   14.066    9.286)   21.953  12.989   (   3.127   -3.127    8.306)    9.410  13.008   (   3.418   -3.418   15.939)   16.656  16.187   (  -0.172    0.172   -2.562)    2.573======================= Grid point 47 (36/220) =======================q-point: ( 0.47  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.753   (   4.931   -4.931    7.299)   10.095   3.028   (   7.903   -7.903   23.113)   25.674   7.444   (  -9.004    9.004    2.956)   13.072  12.983   (   4.361   -4.361    5.897)    8.533  13.105   (   4.637   -4.637   17.086)   18.302  16.159   (  -0.309    0.309   -3.032)    3.064======================= Grid point 48 (37/220) =======================q-point: ( 0.53  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.704   (   5.570   -5.570    5.244)    9.463   3.105   (   8.812   -8.812   23.574)   26.665   7.583   (  -3.206    3.206   -3.422)    5.680  12.931   (   4.763   -4.763    3.645)    7.659  13.189   (   5.875   -5.875   18.292)   20.091  16.138   (  -0.991    0.991   -2.857)    3.182======================= Grid point 49 (38/220) =======================q-point: (-0.42  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.620   (   5.766   -5.766    3.127)    8.733   3.163   (   9.879   -9.879   23.794)   27.593   7.512   (   2.904   -2.904   -9.419)   10.275  12.853   (   4.517   -4.517    1.609)    6.587  13.251   (   7.416   -7.416   18.769)   21.501  16.140   (  -1.818    1.818   -2.663)    3.702======================= Grid point 50 (39/220) =======================q-point: (-0.37  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.512   (   5.674   -5.674    1.002)    8.086   3.195   (  11.237  -11.237   23.675)   28.513   7.236   (   8.866   -8.866  -14.569)   19.222  12.762   (   4.067   -4.067   -0.155)    5.754  13.271   (   9.370   -9.370   17.788)   22.181  16.156   (  -2.200    2.200   -2.833)    4.208======================= Grid point 51 (40/220) =======================q-point: (-0.32  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.382   (   5.630   -5.630   -0.990)    8.023   3.187   (  12.979  -12.979   23.157)   29.549   6.774   (  14.431  -14.431  -18.588)   27.605  12.661   (   3.919   -3.919   -1.558)    5.758  13.220   (  11.551  -11.551   15.038)   22.203  16.168   (  -1.677    1.677   -3.280)    4.048======================= Grid point 52 (41/220) =======================q-point: (-0.26  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.226   (   6.080   -6.080   -2.631)    8.992   3.125   (  15.128  -15.128   22.106)   30.764   6.149   (  19.637  -19.637  -21.240)   34.962  12.542   (   4.394   -4.394   -2.452)    6.681  13.080   (  13.524  -13.524   10.839)   21.984  16.151   (  -0.162    0.162   -3.579)    3.587======================= Grid point 53 (42/220) =======================q-point: (-0.21  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.036   (   7.376   -7.376   -3.553)   11.020   2.993   (  17.544  -17.544   20.008)   31.873   5.386   (  24.702  -24.702  -21.760)   41.157  12.399   (   5.454   -5.454   -2.629)    8.149  12.850   (  14.743  -14.743    6.112)   21.727  16.090   (   2.013   -2.013   -3.259)    4.327======================= Grid point 54 (43/220) =======================q-point: (-0.16  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.802   (   9.555   -9.555   -3.081)   13.859   2.766   (  19.852  -19.852   14.753)   31.716   4.523   (  29.831  -29.831  -17.322)   45.605  12.232   (   6.679   -6.679   -1.863)    9.627  12.553   (  14.711  -14.711    2.099)   20.910  15.986   (   4.207   -4.207   -2.043)    6.291======================= Grid point 55 (44/220) =======================q-point: (-0.11  0.11  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.530   (  12.148  -12.148   -0.000)   17.180   2.382   (  21.964  -21.964   -0.000)   31.061   3.661   (  34.426  -34.426   -0.000)   48.685  12.056   (   7.421   -7.421   -0.000)   10.495  12.238   (  13.155  -13.155   -0.000)   18.604  15.857   (   5.706   -5.706   -0.000)    8.069======================= Grid point 60 (45/220) =======================q-point: ( 0.16  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.079   (  -0.000    0.000   24.345)   24.345   2.079   (  -0.000    0.000   24.345)   24.345   3.728   (  -0.000    0.000   47.825)   47.825  12.273   (  -0.000    0.000   17.309)   17.309  12.273   (  -0.000    0.000   17.309)   17.309  15.943   (  -0.000    0.000    9.986)    9.986======================= Grid point 61 (46/220) =======================q-point: ( 0.21  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.355   (  -0.785    0.785   21.378)   21.407   2.409   (  -4.463    4.463   21.585)   22.489   4.356   (  -8.176    8.176   42.825)   44.359  12.492   (  -1.003    1.003   18.122)   18.178  12.495   (  -1.349    1.349   17.904)   18.005  16.055   (  -0.227    0.227    8.499)    8.505======================= Grid point 62 (47/220) =======================q-point: ( 0.26  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.598   (  -0.357    0.357   18.005)   18.012   2.741   (  -2.859    2.859   20.333)   20.731   5.062   ( -13.425   13.425   35.044)   39.857  12.726   (  -0.929    0.929   18.215)   18.263  12.734   (  -1.488    1.488   17.116)   17.244  16.141   (  -0.066    0.066    5.578)    5.579======================= Grid point 63 (48/220) =======================q-point: ( 0.32  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.783   (   0.797   -0.797   14.915)   14.957   3.005   (   0.464   -0.464   20.732)   20.742   5.761   ( -15.700   15.700   26.816)   34.815  12.946   (  -0.440    0.440   15.297)   15.310  12.947   (   0.145   -0.145   18.199)   18.200  16.184   (   0.153   -0.153    2.063)    2.074======================= Grid point 64 (49/220) =======================q-point: ( 0.37  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.905   (   2.283   -2.283   12.333)   12.749   3.204   (   3.369   -3.369   21.867)   22.380   6.389   ( -15.461   15.461   19.080)   29.019  13.099   (   1.421   -1.421   12.892)   13.048  13.138   (   1.718   -1.718   18.600)   18.758  16.184   (   0.264   -0.264   -1.120)    1.181======================= Grid point 65 (50/220) =======================q-point: ( 0.42  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.964   (   3.834   -3.834   10.196)   11.548   3.361   (   5.527   -5.527   23.188)   24.470   6.902   ( -13.331   13.331   11.790)   22.235  13.175   (   3.530   -3.530   10.277)   11.426  13.300   (   3.305   -3.305   19.664)   20.212  16.151   (   0.136   -0.136   -3.350)    3.355======================= Grid point 66 (51/220) =======================q-point: ( 0.47  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.965   (   5.272   -5.272    8.313)   11.167   3.488   (   7.110   -7.110   24.285)   26.284   7.266   (  -9.723    9.723    4.750)   14.547  13.176   (   5.322   -5.322    7.704)   10.770  13.443   (   4.724   -4.724   21.238)   22.264  16.107   (  -0.336    0.336   -4.527)    4.551======================= Grid point 67 (52/220) =======================q-point: ( 0.53  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.914   (   6.454   -6.454    6.496)   11.203   3.594   (   8.340   -8.340   24.851)   27.509   7.453   (  -4.979    4.979   -1.953)    7.307  13.114   (   6.389   -6.389    5.312)   10.482  13.572   (   6.146   -6.146   22.728)   24.333  16.068   (  -1.172    1.172   -5.001)    5.269======================= Grid point 68 (53/220) =======================q-point: (-0.42  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.820   (   7.248   -7.248    4.649)   11.255   3.676   (   9.398   -9.398   24.735)   28.080   7.448   (   0.484   -0.484   -7.890)    7.920  13.011   (   6.608   -6.608    3.184)    9.873  13.677   (   7.906   -7.906   23.191)   25.745  16.048   (  -2.209    2.209   -5.207)    6.073======================= Grid point 69 (54/220) =======================q-point: (-0.37  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.691   (   7.588   -7.588    2.793)   11.088   3.728   (  10.436  -10.436   23.868)   28.062   7.252   (   6.214   -6.214  -12.512)   15.289  12.889   (   6.175   -6.175    1.393)    8.844  13.730   (  10.170  -10.170   21.644)   25.987  16.049   (  -3.059    3.059   -5.243)    6.797======================= Grid point 70 (55/220) =======================q-point: (-0.32  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.537   (   7.562   -7.562    1.060)   10.746   3.740   (  11.597  -11.597   22.060)   27.489   6.880   (  11.877  -11.877  -15.299)   22.720  12.761   (   5.530   -5.530    0.034)    7.821  13.695   (  12.694  -12.694   17.596)   25.138  16.064   (  -3.203    3.203   -4.825)    6.619======================= Grid point 71 (56/220) =======================q-point: (-0.26  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.368   (   7.460   -7.460   -0.296)   10.554   3.694   (  12.974  -12.974   18.688)   26.190   6.362   (  17.353  -17.353  -15.517)   29.035  12.634   (   5.152   -5.152   -0.764)    7.326  13.545   (  14.907  -14.907   11.576)   24.050  16.080   (  -2.276    2.276   -3.719)    4.918======================= Grid point 72 (57/220) =======================q-point: (-0.21  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.186   (   7.675   -7.675   -0.858)   10.888   3.559   (  14.602  -14.602   12.172)   23.971   5.739   (  22.603  -22.603  -11.443)   33.952  12.504   (   5.321   -5.321   -0.834)    7.571  13.279   (  16.224  -16.224    5.094)   23.503  16.078   (  -0.352    0.352   -2.021)    2.081======================= Grid point 73 (58/220) =======================q-point: (-0.16  0.16  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.995   (   8.445   -8.445   -0.000)   11.943   3.275   (  16.711  -16.711   -0.000)   23.633   5.087   (  27.182  -27.182   -0.000)   38.442  12.369   (   5.967   -5.967   -0.000)    8.439  12.929   (  16.232  -16.232   -0.000)   22.955  16.047   (   2.044   -2.044   -0.000)    2.891======================= Grid point 80 (59/220) =======================q-point: ( 0.21  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.585   (  -0.000    0.000   19.552)   19.552   2.585   (  -0.000    0.000   19.552)   19.552   4.770   (  -0.000    0.000   42.351)   42.351  12.694   (  -0.000    0.000   18.709)   18.709  12.694   (  -0.000    0.000   18.709)   18.709  16.148   (  -0.000    0.000    7.302)    7.302======================= Grid point 81 (60/220) =======================q-point: ( 0.26  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.797   (   0.026   -0.026   16.942)   16.943   2.856   (  -4.483    4.483   17.163)   18.297   5.292   (  -4.426    4.426   38.075)   38.586  12.916   (  -0.540    0.540   18.110)   18.127  12.919   (  -0.782    0.782   18.495)   18.528  16.211   (   0.642   -0.642    4.712)    4.799======================= Grid point 82 (61/220) =======================q-point: ( 0.32  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.970   (   0.691   -0.691   14.275)   14.309   3.161   (  -4.840    4.840   16.119)   17.512   5.835   (  -7.416    7.416   31.673)   33.364  13.125   (  -0.063    0.063   16.407)   16.408  13.151   (  -0.838    0.838   18.205)   18.243  16.224   (   1.303   -1.303    1.447)    2.343======================= Grid point 83 (62/220) =======================q-point: ( 0.37  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.093   (   1.788   -1.788   11.913)   12.179   3.440   (  -2.714    2.714   16.933)   17.363   6.351   (  -9.074    9.074   24.128)   27.329  13.288   (   1.336   -1.336   14.029)   14.156  13.374   (  -0.130    0.130   18.495)   18.496  16.189   (   1.509   -1.509   -1.725)    2.744======================= Grid point 84 (63/220) =======================q-point: ( 0.42  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.162   (   3.147   -3.147    9.917)   10.870   3.676   (   0.076   -0.076   18.745)   18.745   6.797   (  -9.138    9.138   16.178)   20.706  13.381   (   3.312   -3.312   11.374)   12.301  13.584   (   1.063   -1.063   19.685)   19.743  16.122   (   1.148   -1.148   -4.260)    4.559======================= Grid point 85 (64/220) =======================q-point: ( 0.47  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.177   (   4.632   -4.632    8.167)   10.470   3.870   (   2.724   -2.724   20.526)   20.885   7.134   (  -7.593    7.593    8.301)   13.572  13.397   (   5.376   -5.376    8.722)   11.570  13.781   (   2.490   -2.490   21.619)   21.904  16.044   (   0.316   -0.316   -6.046)    6.062======================= Grid point 86 (65/220) =======================q-point: ( 0.53  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.137   (   6.100   -6.100    6.507)   10.805   4.024   (   4.990   -4.990   21.535)   22.662   7.331   (  -4.638    4.638    1.000)    6.636  13.338   (   7.033   -7.033    6.239)   11.741  13.966   (   4.181   -4.181   23.602)   24.331  15.972   (  -0.859    0.859   -7.326)    7.426======================= Grid point 87 (66/220) =======================q-point: (-0.42  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.047   (   7.373   -7.373    4.841)   11.496   4.136   (   6.811   -6.811   21.348)   23.421   7.368   (  -0.600    0.600   -5.033)    5.104  13.223   (   7.918   -7.918    4.029)   11.901  14.124   (   6.430   -6.430   24.389)   26.029  15.918   (  -2.300    2.300   -8.220)    8.840======================= Grid point 88 (67/220) =======================q-point: (-0.37  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.912   (   8.263   -8.263    3.177)   12.110   4.200   (   8.211   -8.211   19.700)   22.868   7.244   (   4.153   -4.153   -9.019)   10.763  13.074   (   7.923   -7.923    2.183)   11.415  14.216   (   9.399   -9.399   22.323)   25.981  15.893   (  -3.810    3.810   -8.124)    9.749======================= Grid point 89 (68/220) =======================q-point: (-0.32  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.743   (   8.640   -8.640    1.630)   12.327   4.207   (   9.333   -9.333   16.180)   20.881   6.976   (   9.271   -9.271  -10.176)   16.597  12.915   (   7.244   -7.244    0.810)   10.276  14.188   (  12.565  -12.565   16.360)   24.154  15.909   (  -4.733    4.733   -6.161)    9.097======================= Grid point 90 (69/220) =======================q-point: (-0.26  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.555   (   8.518   -8.518    0.439)   12.054   4.133   (  10.464  -10.464    9.936)   17.825   6.596   (  14.418  -14.418   -7.525)   21.734  12.762   (   6.311   -6.311    0.036)    8.925  14.009   (  14.998  -14.998    7.911)   22.638  15.964   (  -4.331    4.331   -2.934)    6.792======================= Grid point 91 (70/220) =======================q-point: (-0.21  0.21  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.365   (   8.100   -8.100   -0.000)   11.455   3.935   (  12.057  -12.057   -0.000)   17.051   6.159   (  19.106  -19.106   -0.000)   27.020  12.625   (   5.583   -5.583   -0.000)    7.896  13.690   (  16.256  -16.256   -0.000)   22.989  16.030   (  -2.669    2.669    0.000)    3.774======================= Grid point 100 (71/220) =======================q-point: ( 0.26  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.983   (  -0.000    0.000   15.006)   15.006   2.983   (  -0.000    0.000   15.006)   15.006   5.683   (  -0.000    0.000   36.703)   36.703  13.119   (  -0.000    0.000   17.750)   17.750  13.119   (  -0.000    0.000   17.750)   17.750  16.263   (  -0.000    0.000    2.518)    2.518======================= Grid point 101 (72/220) =======================q-point: ( 0.32  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.140   (   0.405   -0.405   12.797)   12.810   3.204   (  -4.785    4.785   13.059)   14.708   6.111   (  -1.938    1.938   32.827)   32.941  13.316   (   0.141   -0.141   16.261)   16.262  13.330   (  -1.019    1.019   16.799)   16.860  16.263   (   1.180   -1.180   -0.285)    1.693======================= Grid point 102 (73/220) =======================q-point: ( 0.37  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.258   (   1.146   -1.146   10.708)   10.830   3.488   (  -6.764    6.764   12.280)   15.566   6.517   (  -2.954    2.954   27.226)   27.544  13.479   (   1.067   -1.067   14.030)   14.111  13.553   (  -1.725    1.725   16.394)   16.575  16.205   (   2.055   -2.055   -2.995)    4.173======================= Grid point 103 (74/220) =======================q-point: ( 0.42  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.333   (   2.190   -2.190    8.876)    9.401   3.785   (  -6.119    6.119   12.962)   15.585   6.866   (  -3.344    3.344   20.282)   20.826  13.584   (   2.678   -2.678   11.448)   12.058  13.788   (  -1.854    1.854   17.044)   17.244  16.106   (   2.215   -2.215   -5.229)    6.096======================= Grid point 104 (75/220) =======================q-point: ( 0.47  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.361   (   3.503   -3.503    7.246)    8.778   4.061   (  -3.884    3.884   14.380)   15.393   7.130   (  -2.971    2.971   12.561)   13.245  13.616   (   4.686   -4.686    8.812)   11.025  14.032   (  -1.332    1.332   18.806)   18.900  15.989   (   1.664   -1.664   -7.077)    7.458======================= Grid point 105 (76/220) =======================q-point: ( 0.53  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.337   (   5.000   -5.000    5.695)    9.079   4.290   (  -0.939    0.939   15.464)   15.521   7.287   (  -1.732    1.732    5.000)    5.568  13.572   (   6.670   -6.670    6.309)   11.348  14.280   (  -0.005    0.005   21.187)   21.187  15.870   (   0.546   -0.546   -8.990)    9.023======================= Grid point 106 (77/220) =======================q-point: ( 0.58  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.261   (   6.522   -6.522    4.135)   10.108   4.456   (   2.064   -2.064   15.246)   15.523   7.325   (   0.359   -0.359   -1.148)    1.255  13.458   (   8.177   -8.177    4.060)   12.256  14.509   (   2.611   -2.611   22.942)   23.237  15.761   (  -1.263    1.263  -11.107)   11.250======================= Grid point 107 (78/220) =======================q-point: (-0.37  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.133   (   7.847   -7.847    2.577)   11.392   4.542   (   4.669   -4.669   12.957)   14.542   7.248   (   3.250   -3.250   -4.600)    6.502  13.295   (   8.858   -8.858    2.175)   12.715  14.660   (   7.077   -7.077   20.987)   23.251  15.687   (  -4.117    4.117  -11.385)   12.787======================= Grid point 108 (79/220) =======================q-point: (-0.32  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.962   (   8.731   -8.731    1.125)   12.399   4.533   (   6.716   -6.716    7.967)   12.397   7.077   (   6.822   -6.822   -4.352)   10.584  13.109   (   8.604   -8.604    0.778)   12.193  14.634   (  11.727  -11.727   11.525)   20.197  15.707   (  -6.429    6.429   -6.157)   10.981======================= Grid point 109 (80/220) =======================q-point: (-0.26  0.26  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.762   (   8.980   -8.980   -0.000)   12.699   4.407   (   8.461   -8.461   -0.000)   11.966   6.846   (  10.743  -10.743   -0.000)   15.192  12.924   (   7.613   -7.613   -0.000)   10.766  14.396   (  14.083  -14.083   -0.000)   19.916  15.823   (  -6.051    6.051    0.000)    8.557======================= Grid point 120 (81/220) =======================q-point: ( 0.32  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.282   (  -0.000    0.000   10.952)   10.952   3.282   (  -0.000    0.000   10.952)   10.952   6.464   (  -0.000    0.000   30.888)   30.888  13.500   (  -0.000    0.000   15.013)   15.013  13.500   (  -0.000    0.000   15.013)   15.013  16.262   (   0.000   -0.000   -2.543)    2.543======================= Grid point 121 (82/220) =======================q-point: ( 0.37  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.392   (   0.531   -0.531    9.100)    9.131   3.462   (  -5.297    5.297    9.369)   11.995   6.803   (  -0.152    0.152   26.997)   26.998  13.656   (   0.640   -0.640   12.992)   13.024  13.683   (  -1.711    1.711   13.639)   13.852  16.203   (   1.404   -1.404   -4.589)    5.001======================= Grid point 122 (83/220) =======================q-point: ( 0.42  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.467   (   1.260   -1.260    7.414)    7.625   3.730   (  -8.490    8.490    8.684)   14.818   7.083   (   0.400   -0.400   21.573)   21.581  13.764   (   1.849   -1.849   10.560)   10.879  13.898   (  -3.451    3.451   13.233)   14.104  16.097   (   2.393   -2.393   -6.042)    6.925======================= Grid point 123 (84/220) =======================q-point: ( 0.47  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.503   (   2.274   -2.274    5.885)    6.707   4.037   (  -9.015    9.015    8.824)   15.505   7.276   (   1.229   -1.229   14.978)   15.079  13.810   (   3.562   -3.562    8.012)    9.463  14.149   (  -4.714    4.714   13.981)   15.489  15.960   (   2.804   -2.804   -7.088)    8.122======================= Grid point 124 (85/220) =======================q-point: ( 0.53  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.496   (   3.598   -3.598    4.416)    6.737   4.333   (  -7.304    7.304    9.079)   13.752   7.370   (   2.117   -2.117    8.048)    8.587  13.783   (   5.545   -5.545    5.547)    9.605  14.435   (  -5.117    5.117   15.640)   17.232  15.805   (   2.801   -2.801   -8.440)    9.323======================= Grid point 125 (86/220) =======================q-point: ( 0.58  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.437   (   5.139   -5.139    2.924)    7.834   4.569   (  -4.086    4.086    8.398)   10.193   7.369   (   2.923   -2.923    2.214)    4.689  13.683   (   7.418   -7.418    3.304)   10.999  14.737   (  -4.068    4.068   17.748)   18.657  15.634   (   2.299   -2.299  -11.052)   11.520======================= Grid point 126 (87/220) =======================q-point: (-0.37  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.326   (   6.694   -6.694    1.417)    9.572   4.702   (  -0.190    0.190    5.629)    5.636   7.297   (   3.736   -3.736   -0.849)    5.352  13.522   (   8.744   -8.744    1.409)   12.446  15.000   (   1.207   -1.207   18.587)   18.665  15.461   (  -1.304    1.304  -14.259)   14.377======================= Grid point 127 (88/220) =======================q-point: (-0.32  0.32  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.163   (   7.991   -7.991   -0.000)   11.301   4.698   (   3.569   -3.569   -0.000)    5.047   7.189   (   4.931   -4.931   -0.000)    6.974  13.321   (   9.174   -9.174   -0.000)   12.975  14.970   (  10.442  -10.442   -0.000)   14.768  15.491   (  -7.998    7.998    0.000)   11.311======================= Grid point 140 (89/220) =======================q-point: ( 0.37  0.37  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.493   (  -0.000    0.000    7.407)    7.407   3.493   (  -0.000    0.000    7.407)    7.407   7.103   (  -0.000    0.000   24.303)   24.303  13.804   (  -0.000    0.000   11.185)   11.185  13.804   (  -0.000    0.000   11.185)   11.185  16.159   (   0.000   -0.000   -5.863)    5.863======================= Grid point 141 (90/220) =======================q-point: ( 0.42  0.37  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.563   (   0.527   -0.527    5.816)    5.864   3.640   (  -5.851    5.851    6.038)   10.244   7.348   (   1.173   -1.173   19.876)   19.945  13.910   (   0.957   -0.957    8.929)    9.031  13.953   (  -2.722    2.722    9.585)   10.329  16.071   (   1.479   -1.479   -6.368)    6.702======================= Grid point 142 (91/220) =======================q-point: ( 0.47  0.37  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.602   (   1.223   -1.223    4.352)    4.683   3.890   (  -9.819    9.819    5.193)   14.825   7.500   (   2.810   -2.810   14.201)   14.747  13.962   (   2.313   -2.313    6.501)    7.278  14.156   (  -5.623    5.623    8.918)   11.948  15.949   (   2.725   -2.725   -6.161)    7.267======================= Grid point 143 (92/220) =======================q-point: ( 0.53  0.37  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.605   (   2.240   -2.240    2.934)    4.318   4.191   ( -10.944   10.944    4.481)   16.112   7.545   (   4.295   -4.295    8.062)   10.094  13.950   (   4.044   -4.044    4.105)    7.039  14.416   (  -8.036    8.036    8.655)   14.285  15.804   (   3.905   -3.905   -5.683)    7.925======================= Grid point 144 (93/220) =======================q-point: ( 0.58  0.37  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.564   (   3.602   -3.602    1.484)    5.306   4.475   (  -9.281    9.281    3.109)   13.488   7.496   (   5.057   -5.057    2.782)    7.673  13.869   (   5.948   -5.948    1.890)    8.622  14.716   (  -9.604    9.604    7.398)   15.466  15.635   (   5.430   -5.430   -4.988)    9.157======================= Grid point 145 (94/220) =======================q-point: (-0.37  0.37  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.471   (   5.172   -5.172   -0.000)    7.314   4.667   (  -5.355    5.355    0.000)    7.573   7.394   (   4.745   -4.745   -0.000)    6.711  13.722   (   7.661   -7.661   -0.000)   10.834  15.001   (  -9.818    9.818    0.000)   13.885  15.449   (   7.093   -7.093   -0.000)   10.031======================= Grid point 160 (95/220) =======================q-point: ( 0.42  0.42  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.627   (  -0.000    0.000    4.271)    4.271   3.627   (  -0.000    0.000    4.271)    4.271   7.573   (  -0.000    0.000   16.024)   16.024  14.012   (  -0.000    0.000    6.830)    6.830  14.012   (  -0.000    0.000    6.830)    6.830  16.017   (   0.000   -0.000   -5.892)    5.892======================= Grid point 161 (96/220) =======================q-point: ( 0.47  0.42  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.663   (   0.485   -0.485    2.831)    2.913   3.743   (  -6.271    6.271    2.956)    9.348   7.705   (   2.056   -2.056   10.768)   11.154  14.065   (   1.124   -1.124    4.518)    4.789  14.122   (  -3.723    3.723    4.976)    7.245  15.937   (   1.571   -1.571   -4.586)    5.096======================= Grid point 162 (97/220) =======================q-point: ( 0.53  0.42  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.669   (   1.169   -1.169    1.434)    2.188   3.970   ( -10.546   10.546    1.730)   15.014   7.724   (   4.132   -4.132    5.008)    7.696  14.062   (   2.523   -2.523    2.189)    4.185  14.299   (  -7.346    7.346    3.253)   10.886  15.840   (   3.158   -3.158   -2.680)    5.209======================= Grid point 163 (98/220) =======================q-point: (-0.42 -0.58  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.639   (   2.214   -2.214   -0.000)    3.131   4.243   ( -11.620   11.620    0.000)   16.433   7.639   (   5.394   -5.394   -0.000)    7.628  13.998   (   4.194   -4.194   -0.000)    5.932  14.522   (  -9.826    9.826    0.000)   13.897  15.732   (   4.779   -4.779   -0.000)    6.758======================= Grid point 181 (99/220) =======================q-point: ( 0.47  0.47  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 550Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.692   (  -0.000    0.000    1.394)    1.394   3.692   (  -0.000    0.000    1.394)    1.394   7.827   (  -0.000    0.000    5.672)    5.672  14.118   (  -0.000    0.000    2.289)    2.289  14.118   (  -0.000    0.000    2.289)    2.289  15.915   (   0.000   -0.000   -2.472)    2.472======================= Grid point 182 (100/220) =======================q-point: (-0.47 -0.53  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1000Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.695   (   0.464   -0.464   -0.000)    0.656   3.777   (  -6.427    6.427    0.000)    9.089   7.832   (   2.377   -2.377   -0.000)    3.361  14.117   (   1.175   -1.175   -0.000)    1.661  14.180   (  -4.170    4.170    0.000)    5.897  15.881   (   1.640   -1.640   -0.000)    2.319======================= Grid point 406 (101/220) =======================q-point: ( 0.16  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.403   (  -0.000   -1.064   26.666)   26.688   1.829   (  -0.000   24.546   23.350)   33.879   2.928   (  -0.000   27.096   44.624)   52.206  11.941   (  -0.000    2.708   11.918)   12.222  11.999   (  -0.000    7.241   13.326)   15.167  15.750   (  -0.000    3.294    8.472)    9.090======================= Grid point 407 (102/220) =======================q-point: ( 0.21  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.714   (  -5.925   -1.894   22.566)   23.408   2.273   (  13.176   21.961   21.953)   33.732   3.848   ( -18.861   27.782   36.932)   49.916  12.129   (  -2.595    2.447   12.298)   12.805  12.264   (   0.624    9.992   15.102)   18.119  15.895   (  -1.861    3.380    7.617)    8.538======================= Grid point 408 (103/220) =======================q-point: ( 0.26  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.997   (  -6.977   -3.876   18.964)   20.575   2.569   (  17.583   18.487   20.232)   32.562   4.796   ( -23.676   25.984   29.509)   45.896  12.328   (  -3.985    0.992   11.685)   12.386  12.545   (   3.137   11.286   15.416)   19.361  16.035   (  -2.944    2.573    5.350)    6.626======================= Grid point 409 (104/220) =======================q-point: ( 0.32  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.221   (  -6.055   -5.471   16.229)   18.165   2.770   (  19.230   15.195   17.893)   30.345   5.656   ( -23.222   22.820   23.267)   40.017  12.505   (  -4.235   -0.935   10.632)   11.483  12.793   (   6.782   11.095   14.514)   19.488  16.138   (  -3.107    1.206    2.471)    4.149======================= Grid point 410 (105/220) =======================q-point: ( 0.37  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.389   (  -4.289   -6.181   14.350)   16.202   2.894   (  20.077   12.120   15.844)   28.302   6.387   ( -20.701   18.886   17.548)   33.063  12.648   (  -3.609   -2.578    9.683)   10.650  12.972   (  10.551    9.758   12.907)   19.317  16.190   (  -2.595   -0.184   -0.147)    2.606======================= Grid point 411 (106/220) =======================q-point: ( 0.42  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.513   (  -2.084   -5.867   12.933)   14.354   2.957   (  20.400    9.190   14.592)   26.712   6.968   ( -16.973   14.280   11.846)   25.146  12.755   (  -2.376   -3.428    9.062)    9.976  13.073   (  13.481    7.748   11.228)   19.179  16.195   (  -1.892   -1.073   -1.856)    2.860======================= Grid point 412 (107/220) =======================q-point: ( 0.47  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.601   (   0.595   -4.455   11.167)   12.037   2.979   (  19.996    6.086   14.538)   25.460   7.375   ( -12.223    8.973    5.864)   16.258  12.834   (  -0.578   -3.207    8.341)    8.955  13.107   (  14.834    5.331   10.299)   18.829  16.174   (  -1.499   -1.153   -2.439)    3.086======================= Grid point 413 (108/220) =======================q-point: ( 0.53  0.11  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.651   (   4.095   -1.748    8.127)    9.266   2.980   (  18.443    2.234   16.054)   24.554   7.586   (  -6.583    3.071   -0.393)    7.275  12.883   (   2.299   -1.475    6.165)    6.743  13.105   (  13.945    2.089   11.315)   18.079  16.155   (  -1.660   -0.486   -2.169)    2.774======================= Grid point 426 (109/220) =======================q-point: ( 0.21  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.004   (  -0.000   -5.182   24.681)   25.219   2.339   (  -0.000   19.527   20.591)   28.378   3.946   (  -0.000   17.594   42.991)   46.452  12.267   (  -0.000   -0.249   15.988)   15.990  12.356   (  -0.000    7.083   17.196)   18.597  15.957   (  -0.000    1.224    8.935)    9.018======================= Grid point 427 (110/220) =======================q-point: ( 0.26  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.239   (  -4.590   -6.589   22.157)   23.567   2.734   (   8.314   18.585   18.004)   27.178   4.693   ( -11.638   19.590   35.841)   42.471  12.456   (  -2.514   -1.713   15.704)   15.996  12.647   (   2.002    9.636   17.698)   20.250  16.064   (  -0.967    0.535    6.600)    6.692======================= Grid point 428 (111/220) =======================q-point: ( 0.32  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.455   (  -5.139   -8.077   19.809)   22.001   3.001   (  13.870   15.144   17.133)   26.744   5.459   ( -16.952   19.711   27.832)   38.086  12.639   (  -3.760   -3.806   14.950)   15.878  12.924   (   5.689   10.683   17.062)   20.919  16.138   (  -1.624   -0.382    3.344)    3.737======================= Grid point 429 (112/220) =======================q-point: ( 0.37  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.626   (  -3.768   -8.735   17.825)   20.204   3.171   (  16.781   11.412   16.896)   26.407   6.160   ( -17.891   17.973   20.469)   32.590  12.796   (  -3.809   -5.769   14.139)   15.739  13.146   (  10.027   10.194   15.828)   21.330  16.169   (  -1.933   -1.192    0.129)    2.275======================= Grid point 430 (113/220) =======================q-point: ( 0.42  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.752   (  -1.253   -8.180   15.785)   17.823   3.276   (  17.900    7.562   17.440)   26.110   6.747   ( -16.302   14.848   13.650)   25.933  12.919   (  -2.788   -6.896   13.275)   15.217  13.292   (  13.764    8.400   14.776)   21.871  16.162   (  -1.992   -1.583   -2.266)    3.407======================= Grid point 431 (114/220) =======================q-point: ( 0.47  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.829   (   2.249   -6.369   13.101)   14.740   3.343   (  17.512    3.413   18.972)   26.043   7.185   ( -13.005   10.623    7.019)   18.200  13.005   (  -0.449   -6.604   11.629)   13.380  13.371   (  15.526    5.429   15.000)   22.261  16.133   (  -2.026   -1.354   -3.528)    4.287======================= Grid point 432 (115/220) =======================q-point: ( 0.53  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.853   (   6.396   -3.786    9.734)   12.247   3.395   (  15.991   -0.847   21.050)   26.448   7.446   (  -8.372    5.527    0.470)   10.043  13.038   (   3.641   -4.577    8.304)   10.157  13.420   (  14.369    1.250   17.252)   22.487  16.103   (  -2.265   -0.507   -3.868)    4.511======================= Grid point 433 (116/220) =======================q-point: (-0.42  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.820   (  10.162   -1.340    6.388)   12.078   3.441   (  14.337   -4.579   22.639)   27.185   7.512   (  -2.767   -0.096   -5.733)    6.366  13.007   (   7.693   -2.006    4.493)    9.132  13.471   (  12.219   -3.078   19.748)   23.426  16.088   (  -2.745    0.714   -3.800)    4.742======================= Grid point 434 (117/220) =======================q-point: (-0.37  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.739   (  12.747    0.512    3.528)   13.236   3.477   (  13.359   -7.561   23.177)   27.800   7.377   (   3.344   -5.787  -11.101)   12.958  12.928   (   9.617   -0.208    1.685)    9.765  13.511   (  11.701   -6.602   20.218)   24.275  16.095   (  -3.163    1.881   -3.751)    5.254======================= Grid point 435 (118/220) =======================q-point: (-0.32  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.623   (  14.216    1.831    1.142)   14.379   3.488   (  13.167  -10.145   22.714)   28.146   7.051   (   9.553  -11.136  -15.177)   21.110  12.825   (  10.030    0.808   -0.219)   10.064  13.501   (  12.748   -9.611   18.216)   24.222  16.114   (  -2.992    2.521   -3.767)    5.431======================= Grid point 436 (119/220) =======================q-point: (-0.26  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.481   (  15.101    2.617   -0.782)   15.346   3.457   (  13.597  -12.717   21.246)   28.249   6.557   (  15.720  -15.912  -17.564)   28.439  12.711   (  10.065    1.057   -1.434)   10.222  13.406   (  14.492  -12.134   14.071)   23.563  16.128   (  -1.793    2.313   -3.582)    4.626======================= Grid point 437 (120/220) =======================q-point: (-0.21  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.316   (  16.035    2.697   -2.041)   16.388   3.363   (  14.439  -15.481   18.257)   27.954   5.926   (  21.999  -20.038  -17.504)   34.524  12.583   (  10.370    0.494   -1.920)   10.558  13.210   (  16.153  -13.770    8.707)   22.943  16.120   (   0.474    1.234   -2.922)    3.207======================= Grid point 438 (121/220) =======================q-point: (-0.16  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.126   (  17.518    2.027   -2.115)   17.761   3.174   (  15.485  -18.442   12.111)   26.955   5.201   (  28.568  -23.422  -12.937)   39.142  12.438   (  11.070   -0.668   -1.501)   11.191  12.925   (  17.097  -14.071    3.523)   22.421  16.076   (   3.345   -0.410   -1.691)    3.771======================= Grid point 439 (122/220) =======================q-point: (-0.11  0.16  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.914   (  19.753    1.136   -0.000)   19.786   2.834   (  17.275  -21.249   -0.000)   27.385   4.465   (  34.587  -25.872   -0.000)   43.193  12.281   (  11.883   -1.876   -0.000)   12.030  12.589   (  16.841  -12.782   -0.000)   21.142  15.997   (   6.007   -2.118   -0.000)    6.369======================= Grid point 446 (123/220) =======================q-point: ( 0.26  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.529   (  -0.000   -4.664   20.673)   21.193   2.774   (  -0.000   14.881   16.979)   22.578   4.897   (  -0.000   10.721   39.150)   40.592  12.662   (  -0.000   -2.553   17.824)   18.006  12.773   (  -0.000    6.785   18.526)   19.729  16.139   (  -0.000   -0.695    6.441)    6.478======================= Grid point 447 (124/220) =======================q-point: ( 0.32  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.718   (  -2.544   -5.653   19.050)   20.033   3.108   (   3.835   15.031   14.419)   21.179   5.489   (  -6.658   12.287   32.897)   35.743  12.838   (  -2.410   -4.472   17.095)   17.834  13.063   (   2.874    9.110   17.920)   20.307  16.182   (  -0.134   -1.833    3.346)    3.817======================= Grid point 448 (125/220) =======================q-point: ( 0.37  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.885   (  -1.708   -6.008   17.091)   18.197   3.364   (   7.220   11.746   14.475)   19.991   6.075   ( -10.686   12.843   25.347)   30.358  13.002   (  -3.386   -6.394   16.199)   17.741  13.319   (   6.836    9.852   16.880)   20.706  16.179   (  -0.482   -2.605    0.021)    2.649======================= Grid point 449 (126/220) =======================q-point: ( 0.42  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.009   (   0.902   -5.433   14.805)   15.796   3.553   (   9.304    7.209   16.254)   20.068   6.600   ( -11.997   11.941   17.577)   24.402  13.138   (  -2.659   -7.599   15.030)   17.050  13.518   (  10.576    8.777   16.340)   21.352  16.139   (  -1.061   -2.733   -2.782)    4.042======================= Grid point 450 (127/220) =======================q-point: ( 0.47  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.079   (   4.401   -4.223   12.178)   13.620   3.698   (  10.168    2.537   18.732)   21.465   7.020   ( -11.025    9.588    9.986)   17.698  13.229   (   0.195   -7.394   12.869)   14.843  13.661   (  12.410    5.812   17.361)   22.118  16.083   (  -1.754   -2.157   -4.683)    5.446======================= Grid point 451 (128/220) =======================q-point: ( 0.53  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.090   (   8.062   -2.946    9.447)   12.764   3.813   (  10.328   -1.616   20.875)   23.346   7.296   (  -8.196    6.027    2.777)   10.545  13.253   (   4.851   -5.852    9.408)   12.095  13.778   (  11.742    1.478   20.041)   23.274  16.027   (  -2.472   -0.980   -5.773)    6.356======================= Grid point 452 (129/220) =======================q-point: ( 0.58  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.044   (  11.269   -1.953    6.879)   13.347   3.904   (  10.258   -4.967   21.938)   24.721   7.406   (  -3.902    1.611   -3.658)    5.586  13.201   (   8.983   -4.163    5.849)   11.500  13.887   (  10.618   -2.881   22.439)   24.992  15.987   (  -3.204    0.608   -6.375)    7.160======================= Grid point 453 (130/220) =======================q-point: (-0.37  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.947   (  13.682   -1.188    4.563)   14.472   3.966   (  10.151   -7.618   21.659)   25.104   7.339   (   1.418   -3.247   -8.726)    9.418  13.094   (  11.052   -2.846    3.095)   11.825  13.972   (  10.829   -6.798   22.576)   25.946  15.970   (  -3.881    2.343   -6.561)    7.975======================= Grid point 454 (131/220) =======================q-point: (-0.32  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.810   (  15.217   -0.454    2.517)   15.430   3.990   (  10.064   -9.863   19.954)   24.428   7.102   (   7.330   -8.156  -11.804)   16.112  12.961   (  11.441   -1.698    1.155)   11.624  13.985   (  12.334  -10.425   19.518)   25.333  15.980   (  -4.167    3.776   -5.951)    8.188======================= Grid point 455 (132/220) =======================q-point: (-0.26  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.647   (  16.002    0.324    0.847)   16.028   3.960   (  10.100  -11.999   16.461)   22.736   6.719   (  13.509  -12.813  -12.194)   22.256  12.823   (  11.048   -0.747   -0.050)   11.073  13.886   (  14.362  -13.409   13.523)   23.853  16.013   (  -3.546    4.297   -4.305)    7.041======================= Grid point 456 (133/220) =======================q-point: (-0.21  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.472   (  16.360    1.021   -0.165)   16.393   3.847   (  10.511  -14.219   10.205)   20.416   6.229   (  19.679  -17.007   -8.757)   27.444  12.689   (  10.628   -0.234   -0.476)   10.641  13.657   (  16.161  -15.174    6.278)   23.041  16.055   (  -1.799    3.573   -2.100)    4.518======================= Grid point 457 (134/220) =======================q-point: (-0.16  0.21  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.293   (  16.732    1.496   -0.000)   16.799   3.604   (  11.947  -16.522   -0.000)   20.389   5.695   (  25.172  -20.563   -0.000)   32.503  12.561   (  10.618   -0.274   -0.000)   10.622  13.320   (  17.214  -15.413   -0.000)   23.106  16.081   (   0.727    1.829    0.000)    1.968======================= Grid point 466 (135/220) =======================q-point: ( 0.32  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.956   (  -0.000   -2.486   16.384)   16.571   3.123   (  -0.000   11.330   13.293)   17.466   5.747   (  -0.000    5.670   34.441)   34.905  13.073   (  -0.000   -3.740   17.503)   17.898  13.194   (  -0.000    6.379   17.525)   18.650  16.239   (   0.000   -2.067    1.996)    2.873======================= Grid point 467 (136/220) =======================q-point: ( 0.37  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.112   (   0.265   -2.391   14.967)   15.159   3.406   (  -0.663   11.984   11.512)   16.631   6.204   (  -3.340    6.172   28.891)   29.731  13.230   (  -1.986   -5.372   16.518)   17.483  13.460   (   2.807    8.469   16.218)   18.510  16.210   (   0.431   -3.312   -0.985)    3.482======================= Grid point 468 (137/220) =======================q-point: ( 0.42  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.231   (   2.762   -1.754   12.780)   13.193   3.675   (  -0.456    9.305   12.466)   15.563   6.622   (  -5.524    6.169   21.926)   23.437  13.369   (  -1.866   -6.411   15.193)   16.596  13.695   (   5.485    8.791   15.538)   18.676  16.136   (   0.242   -3.821   -3.645)    5.286======================= Grid point 469 (138/220) =======================q-point: ( 0.47  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.300   (   5.728   -1.162   10.447)   11.971   3.915   (   0.830    5.320   14.743)   15.695   6.968   (  -6.275    5.361   14.196)   16.421  13.466   (   0.940   -6.343   12.896)   14.402  13.899   (   6.713    6.904   16.645)   19.230  16.040   (  -0.544   -3.430   -5.748)    6.715======================= Grid point 470 (139/220) =======================q-point: ( 0.53  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.313   (   8.643   -0.929    8.242)   11.979   4.116   (   2.639    1.231   16.844)   17.093   7.209   (  -5.453    3.617    6.469)    9.201  13.491   (   5.531   -5.510    9.598)   12.372  14.090   (   6.494    3.358   19.456)   20.784  15.943   (  -1.656   -2.219   -7.458)    7.955======================= Grid point 471 (140/220) =======================q-point: ( 0.58  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.269   (  11.333   -1.013    6.178)   12.947   4.268   (   4.423   -2.489   17.765)   18.476   7.322   (  -3.108    1.041   -0.331)    3.295  13.434   (   9.644   -4.825    6.283)   12.480  14.272   (   6.523   -0.792   22.174)   23.127  15.859   (  -2.892   -0.298   -9.003)    9.461======================= Grid point 472 (141/220) =======================q-point: (-0.37  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.172   (  13.627   -1.202    4.225)   14.318   4.363   (   5.800   -5.670   16.941)   18.783   7.297   (   0.565   -2.171   -5.197)    5.661  13.310   (  12.029   -4.334    3.593)   13.281  14.416   (   7.970   -5.412   22.466)   24.444  15.801   (  -4.246    2.295   -9.817)   10.939======================= Grid point 473 (142/220) =======================q-point: (-0.32  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.028   (  15.339   -1.237    2.422)   15.578   4.391   (   6.634   -8.294   13.986)   17.561   7.144   (   5.266   -5.823   -7.217)   10.664  13.149   (  12.731   -3.611    1.634)   13.334  14.457   (  10.614  -10.241   17.944)   23.227  15.793   (  -5.385    4.977   -8.106)   10.930======================= Grid point 474 (143/220) =======================q-point: (-0.26  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.851   (  16.329   -0.912    0.915)   16.380   4.334   (   7.176  -10.482    8.445)   15.254   6.885   (  10.542   -9.698   -5.630)   15.391  12.979   (  12.295   -2.521    0.417)   12.558  14.337   (  13.195  -13.751    9.047)   21.096  15.853   (  -5.335    6.167   -3.839)    9.013======================= Grid point 475 (144/220) =======================q-point: (-0.21  0.26  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.658   (  16.600   -0.150   -0.000)   16.601   4.166   (   8.144  -12.429   -0.000)   14.860   6.563   (  15.675  -13.523   -0.000)   20.702  12.822   (  11.392   -1.289   -0.000)   11.465  14.055   (  15.060  -15.248   -0.000)   21.431  15.955   (  -3.921    5.422    0.000)    6.691======================= Grid point 486 (145/220) =======================q-point: ( 0.37  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.287   (  -0.000    0.114   12.249)   12.249   3.390   (  -0.000    8.947    9.942)   13.375   6.483   (  -0.000    1.841   29.198)   29.256  13.456   (  -0.000   -3.625   15.413)   15.833  13.569   (  -0.000    5.929   14.782)   15.927  16.229   (   0.000   -2.765   -2.743)    3.895======================= Grid point 487 (146/220) =======================q-point: ( 0.42  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.406   (   3.502    1.111   10.575)   11.195   3.643   (  -5.193   10.113    9.099)   14.562   6.816   (  -1.047    1.139   23.943)   23.993  13.587   (  -0.794   -4.431   14.173)   14.870  13.802   (   1.317    7.799   13.302)   15.476  16.141   (   0.741   -3.949   -4.819)    6.275======================= Grid point 488 (147/220) =======================q-point: ( 0.47  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.476   (   6.131    1.525    8.498)   10.589   3.936   (  -6.524    9.034    9.975)   14.955   7.078   (  -1.492    0.213   17.352)   17.418  13.687   (   1.475   -4.381   11.985)   12.846  14.033   (   1.024    7.837   13.618)   15.745  16.017   (   0.766   -4.371   -6.417)    7.803======================= Grid point 489 (148/220) =======================q-point: ( 0.53  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.496   (   8.386    1.270    6.625)   10.762   4.218   (  -5.214    6.401   11.179)   13.897   7.250   (  -1.323   -0.851   10.096)   10.218  13.719   (   5.623   -4.029    8.892)   11.266  14.278   (  -0.019    6.257   15.875)   17.063  15.879   (   0.211   -4.002   -7.975)    8.926======================= Grid point 490 (149/220) =======================q-point: ( 0.58  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.463   (  10.607    0.604    4.851)   11.679   4.450   (  -2.590    2.742   11.664)   12.259   7.321   (  -0.468   -1.978    3.372)    3.937  13.669   (   9.516   -4.145    5.785)   11.883  14.538   (   0.020    3.556   18.855)   19.187  15.738   (  -0.775   -2.758  -10.209)   10.603======================= Grid point 491 (150/220) =======================q-point: (-0.37  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.377   (  12.738   -0.226    3.106)   13.113   4.601   (   0.352   -1.319   10.411)   10.500   7.298   (   1.181   -3.197   -1.301)    3.648  13.543   (  12.155   -4.529    3.179)   13.355  14.774   (   2.592   -1.290   20.675)   20.877  15.606   (  -2.822    0.333  -12.875)   13.184======================= Grid point 492 (151/220) =======================q-point: (-0.32  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.238   (  14.563   -0.945    1.437)   14.664   4.645   (   2.929   -5.138    6.600)    8.862   7.199   (   3.731   -4.745   -2.595)    6.570  13.366   (  13.344   -4.576    1.225)   14.159  14.868   (   8.236   -9.256   14.242)   18.877  15.557   (  -6.540    6.139   -9.480)   13.051======================= Grid point 493 (152/220) =======================q-point: (-0.26  0.32  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.056   (  15.830   -1.264   -0.000)   15.880   4.561   (   5.017   -8.234   -0.000)    9.642   7.052   (   7.045   -6.906   -0.000)    9.866  13.168   (  13.252   -3.901   -0.000)   13.814  14.695   (  11.570  -13.170   -0.000)   17.530  15.673   (  -6.980    7.709    0.000)   10.399======================= Grid point 506 (153/220) =======================q-point: ( 0.42  0.37  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.524   (  -0.000    2.399    8.383)    8.720   3.585   (  -0.000    7.648    6.958)   10.339   7.089   (  -0.000   -1.125   23.017)   23.045  13.775   (  -0.000   -2.381   12.054)   12.287  13.868   (  -0.000    5.505   10.985)   12.287  16.125   (   0.000   -2.977   -5.905)    6.613======================= Grid point 507 (154/220) =======================q-point: ( 0.47  0.37  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.603   (   5.861    3.469    6.567)    9.462   3.824   (  -8.538    9.766    6.485)   14.503   7.299   (   0.575   -2.845   17.618)   17.856  13.872   (   1.609   -2.220   10.282)   10.641  14.070   (  -1.939    7.287    9.797)   12.363  16.004   (   0.993   -4.207   -6.525)    7.827======================= Grid point 508 (155/220) =======================q-point: ( 0.53  0.37  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.629   (   7.711    3.191    4.869)    9.662   4.127   ( -10.048   10.197    6.430)   15.694   7.413   (   1.408   -4.445   11.343)   12.263  13.912   (   5.250   -2.005    7.395)    9.288  14.315   (  -4.907    8.060   10.518)   14.131  15.858   (   1.629   -5.040   -6.846)    8.656======================= Grid point 509 (156/220) =======================q-point: ( 0.58  0.37  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.609   (   9.471    2.346    3.252)   10.285   4.418   (  -8.676    8.273    5.992)   13.402   7.429   (   2.157   -5.441    5.276)    7.880  13.873   (   8.693   -2.571    4.455)   10.101  14.607   (  -6.749    8.245   11.759)   15.868  15.691   (   2.286   -5.688   -7.596)    9.761======================= Grid point 510 (157/220) =======================q-point: ( 0.63  0.37  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.538   (  11.339    1.193    1.615)   11.515   4.636   (  -5.422    4.393    4.183)    8.136   7.371   (   2.620   -5.523    0.993)    6.193  13.758   (  11.324   -3.626    1.934)   12.047  14.921   (  -6.962    7.484   12.393)   16.064  15.494   (   3.113   -6.151   -9.255)   11.541======================= Grid point 511 (158/220) =======================q-point: (-0.32  0.37  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.413   (  13.130   -0.010   -0.000)   13.130   4.725   (  -1.238   -0.535    0.000)    1.349   7.280   (   2.963   -4.956   -0.000)    5.775  13.580   (  12.909   -4.490   -0.000)   13.668  15.142   (   1.815   -2.855   -0.000)    3.383  15.331   (  -3.413    2.150    0.000)    4.034======================= Grid point 526 (159/220) =======================q-point: ( 0.47  0.42  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.676   (  -0.000    3.956    4.824)    6.239   3.713   (  -0.000    7.179    4.175)    8.305   7.533   (  -0.000   -3.315   15.152)   15.510  14.006   (  -0.000   -0.545    7.784)    7.803  14.073   (  -0.000    5.173    6.692)    8.458  15.981   (   0.000   -3.077   -5.954)    6.702======================= Grid point 527 (160/220) =======================q-point: ( 0.53  0.42  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.714   (   6.900    4.591    3.149)    8.866   3.939   ( -10.211   10.220    3.366)   14.834   7.615   (   1.612   -5.542    9.504)   11.120  14.054   (   4.397   -0.049    5.299)    6.886  14.249   (  -5.808    7.316    5.566)   10.874  15.866   (   1.364   -4.579   -4.843)    6.803======================= Grid point 528 (161/220) =======================q-point: ( 0.58  0.42  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.704   (   8.269    3.864    1.597)    9.266   4.227   ( -11.625   11.202    2.211)   16.295   7.591   (   2.976   -7.067    3.974)    8.637  14.026   (   7.374   -0.715    2.465)    7.809  14.494   (  -9.298    9.444    4.350)   13.948  15.733   (   2.893   -6.127   -3.251)    7.515======================= Grid point 529 (162/220) =======================q-point: (-0.37 -0.58  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.647   (   9.760    2.656   -0.000)   10.115   4.488   (  -9.823    9.104    0.000)   13.393   7.490   (   3.429   -7.153   -0.000)    7.933  13.925   (   9.666   -2.133   -0.000)    9.898  14.760   ( -10.736   10.661    0.000)   15.130  15.587   (   4.611   -7.733   -0.000)    9.004======================= Grid point 547 (163/220) =======================q-point: ( 0.53  0.47  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.749   (   0.000    4.693    1.565)    4.947   3.777   (  -0.000    7.126    1.401)    7.262   7.773   (  -0.000   -4.541    5.352)    7.018  14.129   (   0.000    0.913    2.730)    2.879  14.176   (  -0.000    4.986    2.238)    5.466  15.877   (   0.000   -3.240   -2.530)    4.111======================= Grid point 548 (164/220) =======================q-point: (-0.42 -0.53  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.750   (   7.167    4.916   -0.000)    8.691   3.978   ( -10.698   10.486    0.000)   14.980   7.727   (   1.980   -6.519   -0.000)    6.813  14.115   (   5.559    0.783   -0.000)    5.614  14.316   (  -7.524    7.580    0.000)   10.680  15.806   (   1.609   -4.862   -0.000)    5.121======================= Grid point 813 (165/220) =======================q-point: ( 0.32  0.21  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.394   (  -0.000   -6.152   21.092)   21.971   3.186   (  -0.000   19.126   13.349)   23.324   5.239   (  -0.000   18.150   32.190)   36.954  12.595   (  -0.000   -2.726   15.050)   15.295  12.996   (  -0.000   12.199   18.004)   21.748  16.117   (  -0.000   -1.254    4.047)    4.237======================= Grid point 814 (166/220) =======================q-point: ( 0.37  0.21  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.569   (  -1.207   -6.370   20.085)   21.105   3.485   (   8.341   16.701   11.639)   21.999   5.832   (  -8.581   16.752   25.927)   32.038  12.736   (  -1.321   -3.809   14.487)   15.038  13.325   (   4.469   13.241   17.167)   22.136  16.128   (  -1.172   -2.802    0.973)    3.189======================= Grid point 815 (167/220) =======================q-point: ( 0.42  0.21  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.733   (  -0.392   -6.177   19.136)   20.112   3.662   (  13.835   13.408   11.664)   22.522   6.398   ( -12.543   14.525   18.821)   26.880  12.870   (  -1.438   -4.440   14.220)   14.966  13.582   (  10.126   12.591   15.785)   22.588  16.106   (  -2.456   -3.656   -1.810)    4.761======================= Grid point 816 (168/220) =======================q-point: ( 0.47  0.21  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.874   (   1.645   -5.246   17.797)   18.627   3.757   (  16.854    9.914   12.594)   23.258   6.875   ( -12.958   11.280   11.758)   20.818  12.994   (  -0.411   -4.283   13.993)   14.639  13.742   (  15.297   10.383   14.548)   23.526  16.066   (  -3.598   -3.554   -3.713)    6.274======================= Grid point 817 (169/220) =======================q-point: ( 0.53  0.21  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.981   (   4.594   -3.575   15.547)   16.601   3.801   (  18.007    6.197   14.152)   23.726   7.217   ( -10.850    7.164    4.950)   13.912  13.099   (   1.873   -3.224   12.970)   13.496  13.812   (  18.296    6.921   14.323)   24.244  16.027   (  -4.353   -2.533   -4.618)    6.833======================= Grid point 818 (170/220) =======================q-point: ( 0.58  0.21  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.040   (   8.280   -1.282   12.209)   14.808   3.817   (  17.589    2.131   16.039)   23.898   7.397   (  -6.826    2.458   -1.387)    7.387  13.166   (   5.652   -1.243   10.069)   11.614  13.826   (  18.127    2.469   15.853)   24.207  16.004   (  -4.717   -0.910   -4.831)    6.813======================= Grid point 833 (171/220) =======================q-point: ( 0.37  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.864   (  -0.000   -5.406   19.224)   19.970   3.456   (  -0.000   16.426   10.247)   19.360   5.943   (  -0.000   11.072   28.651)   30.716  12.963   (  -0.000   -5.331   16.547)   17.384  13.403   (  -0.000   11.471   16.909)   20.433  16.171   (   0.000   -3.697    0.526)    3.735======================= Grid point 834 (172/220) =======================q-point: ( 0.42  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.023   (   0.351   -5.098   18.703)   19.389   3.722   (   4.235   14.665    9.065)   17.753   6.392   (  -5.082    9.891   22.465)   25.067  13.094   (  -1.003   -5.825   16.253)   17.294  13.704   (   4.557   12.316   15.405)   20.242  16.115   (  -0.755   -4.865   -2.236)    5.407======================= Grid point 835 (173/220) =======================q-point: ( 0.47  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.162   (   2.627   -4.184   17.265)   17.958   3.913   (   7.471   11.057   10.274)   16.841   6.792   (  -7.788    8.213   15.250)   18.991  13.223   (  -0.385   -5.690   15.903)   16.895  13.931   (   9.514   11.287   14.325)   20.570  16.034   (  -2.136   -4.972   -4.470)    7.019======================= Grid point 836 (174/220) =======================q-point: ( 0.53  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.258   (   6.055   -2.602   14.753)   16.158   4.048   (   9.362    6.458   12.688)   17.039   7.101   (  -8.049    5.792    7.884)   12.668  13.331   (   2.142   -4.746   14.561)   15.464  14.079   (  13.020    8.243   14.667)   21.275  15.954   (  -3.711   -3.909   -6.038)    8.094======================= Grid point 837 (175/220) =======================q-point: ( 0.58  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.298   (  10.010   -0.706   11.478)   15.246   4.142   (  10.045    1.527   14.991)   18.110   7.285   (  -6.085    2.722    1.115)    6.758  13.388   (   6.614   -3.017   11.273)   13.414  14.171   (  13.632    3.423   16.911)   21.989  15.894   (  -4.922   -1.872   -7.077)    8.821======================= Grid point 838 (176/220) =======================q-point: (-0.37  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.275   (  13.817    1.059    8.057)   16.030   4.197   (   9.848   -3.145   16.061)   19.101   7.326   (  -2.222   -0.757   -4.310)    4.908  13.364   (  11.431   -1.052    6.831)   13.358  14.238   (  12.425   -2.225   19.279)   23.043  15.863   (  -5.571    0.745   -7.576)    9.434======================= Grid point 839 (177/220) =======================q-point: (-0.32  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.193   (  16.897    2.488    4.967)   17.787   4.210   (   9.106   -7.226   15.295)   19.212   7.221   (   3.102   -4.382   -7.690)    9.378  13.266   (  14.360    0.567    3.158)   14.714  14.264   (  11.926   -7.556   18.587)   23.341  15.867   (  -5.718    3.423   -6.987)    9.656======================= Grid point 840 (178/220) =======================q-point: (-0.26  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.062   (  18.947    3.634    2.435)   19.445   4.168   (   8.181  -10.690   12.472)   18.351   6.983   (   9.348   -7.906   -8.424)   14.862  13.131   (  15.195    1.873    0.883)   15.335  14.200   (  12.791  -11.723   13.795)   22.167  15.913   (  -5.204    5.328   -4.892)    8.911======================= Grid point 841 (179/220) =======================q-point: (-0.21  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.902   (  19.971    4.595    0.651)   20.503   4.054   (   7.642  -13.606    7.387)   17.265   6.641   (  15.872  -11.127   -5.968)   20.282  12.987   (  14.937    2.913   -0.140)   15.219  14.012   (  14.378  -14.063    6.707)   21.201  15.991   (  -3.791    5.752   -2.148)    7.216======================= Grid point 842 (180/220) =======================q-point: (-0.16  0.26  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.729   (  20.223    5.431   -0.000)   20.940   3.837   (   8.433  -15.823   -0.000)   17.930   6.239   (  21.782  -13.998   -0.000)   25.892  12.853   (  14.521    3.568   -0.000)   14.953  13.707   (  16.039  -14.484   -0.000)   21.611  16.073   (  -1.590    4.755    0.000)    5.014======================= Grid point 853 (181/220) =======================q-point: ( 0.42  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.269   (  -0.000   -1.949   15.689)   15.810   3.665   (  -0.000   14.064    7.907)   16.134   6.557   (  -0.000    4.644   24.373)   24.812  13.345   (  -0.000   -5.622   16.262)   17.206  13.764   (  -0.000   10.699   14.164)   17.751  16.138   (   0.000   -4.990   -3.338)    6.004======================= Grid point 854 (182/220) =======================q-point: ( 0.47  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.412   (   3.131   -1.029   14.696)   15.061   3.913   (  -0.531   12.778    7.597)   14.875   6.864   (  -2.420    3.156   18.281)   18.708  13.470   (  -0.096   -5.317   15.952)   16.815  14.027   (   3.334   11.404   12.538)   17.273  16.026   (  -0.451   -5.828   -5.287)    7.882======================= Grid point 855 (183/220) =======================q-point: ( 0.53  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.506   (   7.178    0.400   12.232)   14.188   4.143   (  -0.105    8.915    9.580)   13.086   7.099   (  -3.645    1.653   11.266)   11.956  13.581   (   2.352   -4.414   14.566)   15.401  14.242   (   5.803    9.980   12.708)   17.168  15.901   (  -1.741   -5.559   -6.930)    9.053======================= Grid point 856 (184/220) =======================q-point: ( 0.58  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.537   (  11.010    1.649    9.312)   14.515   4.334   (   1.292    3.921   11.627)   12.338   7.241   (  -3.424    0.115    4.391)    5.569  13.635   (   7.165   -3.032   11.193)   13.631  14.425   (   6.493    6.124   15.392)   17.793  15.784   (  -3.410   -3.978   -8.665)   10.126======================= Grid point 857 (185/220) =======================q-point: ( 0.63  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.505   (  14.369    2.422    6.477)   15.946   4.464   (   2.914   -1.215   12.168)   12.570   7.282   (  -1.605   -1.464   -1.066)    2.420  13.599   (  12.188   -1.790    6.853)   14.097  14.585   (   6.726    0.231   18.511)   19.697  15.693   (  -5.054   -0.793  -10.396)   11.587======================= Grid point 858 (186/220) =======================q-point: (-0.32  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.413   (  17.110    2.841    3.887)   17.775   4.515   (   4.108   -5.893   10.456)   12.685   7.223   (   1.794   -3.218   -3.957)    5.407  13.480   (  15.289   -0.912    3.288)   15.665  14.675   (   8.270   -6.928   17.346)   20.428  15.657   (  -6.484    3.862   -9.819)   12.384======================= Grid point 859 (187/220) =======================q-point: (-0.26  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.269   (  19.042    3.172    1.655)   19.375   4.471   (   4.778   -9.684    6.299)   12.501   7.081   (   6.424   -5.360   -3.602)    9.109  13.314   (  16.259    0.049    1.027)   16.291  14.611   (  10.441  -12.423    9.354)   18.731  15.714   (  -6.744    7.158   -4.843)   10.962======================= Grid point 860 (188/220) =======================q-point: (-0.21  0.32  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.091   (  20.003    3.682   -0.000)   20.339   4.322   (   5.580  -12.311   -0.000)   13.517   6.877   (  11.481   -8.012   -0.000)   14.000  13.141   (  15.852    1.364   -0.000)   15.911  14.374   (  12.439  -14.304   -0.000)   18.956  15.846   (  -5.570    7.237    0.000)    9.132======================= Grid point 873 (189/220) =======================q-point: ( 0.47  0.37  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.582   (  -0.000    2.048   11.293)   11.478   3.825   (  -0.000   12.545    5.964)   13.891   7.061   (  -0.000   -0.970   19.139)   19.163  13.700   (  -0.000   -3.903   14.228)   14.753  14.050   (  -0.000    9.990   10.455)   14.460  16.026   (   0.000   -5.467   -6.086)    8.181======================= Grid point 874 (190/220) =======================q-point: ( 0.53  0.37  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.691   (   7.025    3.431    9.394)   12.222   4.074   (  -5.847   11.526    6.268)   14.364   7.229   (  -0.457   -2.844   13.173)   13.484  13.812   (   2.134   -2.875   13.183)   13.661  14.280   (   0.092   10.638    9.302)   14.132  15.881   (  -0.007   -6.338   -6.942)    9.400======================= Grid point 875 (191/220) =======================q-point: ( 0.58  0.37  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.724   (  11.088    4.299    6.781)   13.690   4.342   (  -6.292    8.201    7.305)   12.657   7.310   (  -0.435   -4.087    6.918)    8.047  13.869   (   7.119   -1.687    9.898)   12.308  14.519   (  -1.397    9.293   11.101)   14.544  15.723   (  -0.497   -6.431   -8.078)   10.337======================= Grid point 876 (192/220) =======================q-point: ( 0.63  0.37  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.693   (  13.965    4.330    4.362)   15.258   4.557   (  -4.061    3.522    7.090)    8.898   7.310   (   0.097   -4.450    1.761)    4.787  13.831   (  12.116   -1.102    5.680)   13.427  14.779   (  -1.961    5.982   14.591)   15.892  15.558   (  -1.336   -5.173  -10.509)   11.789======================= Grid point 877 (193/220) =======================q-point: (-0.32  0.37  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.603   (  16.327    3.911    2.082)   16.917   4.667   (  -1.008   -1.814    4.747)    5.181   7.255   (   1.348   -4.183   -0.809)    4.469  13.703   (  15.275   -1.109    2.232)   15.477  14.996   (   2.299   -2.782   15.900)   16.304  15.418   (  -4.864    1.420  -12.720)   13.692======================= Grid point 878 (194/220) =======================q-point: (-0.26  0.37  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.457   (  18.100    3.403   -0.000)   18.417   4.641   (   1.942   -6.600   -0.000)    6.880   7.173   (   3.531   -4.175   -0.000)    5.468  13.518   (  16.497   -0.918   -0.000)   16.523  14.925   (   8.269  -12.231   -0.000)   14.764  15.498   (  -8.236    8.742    0.000)   12.011======================= Grid point 893 (195/220) =======================q-point: ( 0.53  0.42  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.789   (   0.000    5.242    6.636)    8.457   3.939   (  -0.000   11.865    3.911)   12.493   7.430   (  -0.000   -5.345   12.492)   13.588  13.987   (   0.000   -1.035   10.288)   10.340  14.244   (  -0.000    9.445    6.330)   11.370  15.877   (   0.000   -5.808   -6.238)    8.523======================= Grid point 894 (196/220) =======================q-point: ( 0.58  0.42  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.849   (  10.036    6.371    4.376)   12.667   4.192   (  -9.721   11.253    3.724)   15.330   7.464   (   0.805   -7.141    6.981)   10.018  14.056   (   5.917    0.183    7.528)    9.576  14.455   (  -5.226   10.296    5.693)   12.873  15.729   (   0.875   -7.176   -5.560)    9.120======================= Grid point 895 (197/220) =======================q-point: ( 0.63  0.42  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.826   (  12.874    6.112    2.150)   14.413   4.462   (  -9.318    8.551    2.732)   12.939   7.417   (   1.376   -7.533    2.333)    8.005  14.024   (  11.013    0.384    3.373)   11.525  14.727   (  -8.935   10.498    5.727)   14.928  15.562   (   2.245   -8.520   -4.682)    9.978======================= Grid point 896 (198/220) =======================q-point: (-0.32  0.42  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.743   (  14.835    5.139   -0.000)   15.700   4.642   (  -5.863    3.755    0.000)    6.962   7.330   (   1.457   -6.252   -0.000)    6.419  13.897   (  13.962   -0.396   -0.000)   13.968  15.005   (  -9.684   10.135    0.000)   14.018  15.381   (   3.928   -9.663   -0.000)   10.430======================= Grid point 913 (199/220) =======================q-point: ( 0.58  0.47  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.890   (  -0.000    6.947    2.156)    7.274   4.001   (   0.000   11.694    1.396)   11.777   7.627   (  -0.000   -7.816    4.381)    8.960  14.156   (   0.000    1.403    3.934)    4.177  14.341   (  -0.000    9.147    2.108)    9.387  15.766   (   0.000   -6.313   -2.776)    6.897======================= Grid point 914 (200/220) =======================q-point: (-0.37 -0.53  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.899   (  10.996    7.293   -0.000)   13.194   4.237   ( -10.971   11.355    0.000)   15.790   7.545   (   1.250   -8.722   -0.000)    8.811  14.145   (   8.135    1.627   -0.000)    8.296  14.525   (  -8.407   10.413    0.000)   13.383  15.658   (   1.557   -7.869   -0.000)    8.022======================= Grid point 1219 (201/220) =======================q-point: ( 0.47  0.32  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.193   (  -0.000   -4.414   19.315)   19.813   4.005   (  -0.000   14.750    6.327)   16.049   6.696   (  -0.000    7.255   17.748)   19.174  13.218   (  -0.000   -4.895   16.707)   17.410  14.048   (  -0.000   13.559   13.328)   19.012  16.005   (   0.000   -6.440   -4.302)    7.745======================= Grid point 1220 (202/220) =======================q-point: ( 0.53  0.32  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.355   (   2.178   -3.632   18.624)   19.099   4.210   (   3.970   12.525    6.088)   14.481   6.959   (  -3.496    4.918   11.706)   13.170  13.361   (   0.599   -3.730   16.999)   17.414  14.313   (   5.948   12.997   11.599)   18.408  15.878   (  -2.171   -6.897   -5.987)    9.388======================= Grid point 1221 (203/220) =======================q-point: ( 0.58  0.32  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.481   (   5.455   -2.449   16.488)   17.539   4.345   (   6.902    8.402    7.211)   13.047   7.152   (  -4.763    2.817    5.406)    7.737  13.499   (   3.080   -2.413   16.038)   16.508  14.476   (  11.387    9.840   11.173)   18.743  15.770   (  -5.268   -5.574   -7.285)   10.579======================= Grid point 1222 (204/220) =======================q-point: ( 0.63  0.32  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.548   (   9.301   -0.870   13.187)   16.161   4.412   (   8.375    2.962    8.323)   12.174   7.254   (  -3.558    0.911    0.007)    3.673  13.588   (   7.555   -0.865   12.702)   14.804  14.544   (  13.615    3.747   12.458)   18.831  15.710   (  -7.593   -2.176   -8.059)   11.284======================= Grid point 1239 (205/220) =======================q-point: ( 0.53  0.37  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.598   (   0.000   -0.862   15.417)   15.441   4.135   (  -0.000   13.676    4.936)   14.539   7.055   (  -0.000    0.612   13.248)   13.262  13.607   (   0.000   -3.900   16.563)   17.016  14.316   (  -0.000   12.749    9.760)   16.056  15.878   (  -0.000   -7.239   -6.531)    9.750======================= Grid point 1240 (206/220) =======================q-point: ( 0.58  0.37  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.735   (   5.062    0.008   13.770)   14.671   4.338   (  -0.407   10.935    4.973)   12.020   7.181   (  -1.478   -1.081    7.531)    7.750  13.748   (   2.463   -2.514   15.850)   16.236  14.546   (   3.195   12.081    8.607)   15.174  15.718   (  -1.816   -7.625   -7.682)   10.975======================= Grid point 1241 (207/220) =======================q-point: ( 0.63  0.37  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.796   (   9.970    1.451   10.479)   14.536   4.499   (   0.256    5.273    5.895)    7.913   7.239   (  -1.677   -1.732    2.341)    3.360  13.836   (   7.644   -1.105   12.383)   14.594  14.724   (   5.044    7.891   10.552)   14.109  15.576   (  -4.600   -5.749   -9.342)   11.895======================= Grid point 1242 (208/220) =======================q-point: (-0.32  0.37  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.781   (  13.993    2.659    6.740)   15.757   4.579   (   1.496   -1.473    5.644)    6.022   7.238   (  -0.232   -1.630   -1.014)    1.934  13.821   (  13.381    0.162    6.996)   15.101  14.842   (   6.255   -1.326   13.112)   14.588  15.494   (  -7.403    0.479  -10.406)   12.779======================= Grid point 1243 (209/220) =======================q-point: (-0.26  0.37  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.693   (  16.937    3.329    3.135)   17.544   4.557   (   2.610   -7.502    3.522)    8.689   7.187   (   2.941   -1.631   -1.756)    3.794  13.704   (  16.823    1.093    2.479)   17.039  14.818   (   7.489  -10.838    8.163)   15.497  15.540   (  -8.310    7.417   -5.634)   12.483======================= Grid point 1244 (210/220) =======================q-point: (-0.21  0.37  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.544   (  18.548    3.540   -0.000)   18.883   4.430   (   3.812  -11.395   -0.000)   12.016   7.093   (   7.279   -2.624   -0.000)    7.738  13.534   (  17.662    2.159   -0.000)   17.793  14.623   (   9.054  -13.515   -0.000)   16.267  15.697   (  -7.319    8.340    0.000)   11.096======================= Grid point 1259 (211/220) =======================q-point: ( 0.58  0.42  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.892   (  -0.000    3.180    9.807)   10.310   4.231   (  -0.000   12.736    3.347)   13.169   7.304   (  -0.000   -5.040    8.257)    9.673  13.959   (   0.000   -1.355   13.358)   13.427  14.497   (  -0.000   12.150    5.871)   13.494  15.717   (  -0.000   -7.964   -6.933)   10.559======================= Grid point 1260 (212/220) =======================q-point: ( 0.63  0.42  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.977   (   8.896    4.111    7.046)   12.070   4.434   (  -5.220    9.172    3.131)   11.008   7.309   (  -0.086   -5.682    3.621)    6.738  14.058   (   6.195   -0.064   10.210)   11.942  14.712   (  -2.628   11.624    5.678)   13.201  15.541   (  -0.320   -8.934   -6.986)   11.345======================= Grid point 1261 (213/220) =======================q-point: ( 0.68  0.42  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.957   (  13.474    4.689    3.447)   14.677   4.600   (  -4.523    3.031    2.447)    5.969   7.271   (   0.188   -4.551    0.559)    4.589  14.038   (  12.661    0.708    4.577)   13.481  14.958   (  -4.815    8.306    8.522)   12.838  15.358   (  -0.580   -8.022   -8.277)   11.541======================= Grid point 1262 (214/220) =======================q-point: (-0.26  0.42  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.858   (  15.930    4.295   -0.000)   16.499   4.648   (  -1.352   -3.274   -0.000)    3.542   7.223   (   1.109   -2.581   -0.000)    2.809  13.902   (  16.064    0.774   -0.000)   16.082  15.069   (   3.860   -9.670   -0.000)   10.412  15.300   (  -8.260    7.813    0.000)   11.370======================= Grid point 1279 (215/220) =======================q-point: ( 0.63  0.47  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   4.045   (  -0.000    5.976    3.366)    6.858   4.284   (   0.000   12.012    1.192)   12.071   7.433   (  -0.000   -8.381    2.816)    8.841  14.188   (   0.000    1.232    5.583)    5.718  14.587   (  -0.000   11.836    1.948)   11.996  15.588   (  -0.000   -9.001   -3.382)    9.615======================= Grid point 1280 (216/220) =======================q-point: (-0.32 -0.53  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   4.058   (  10.817    5.928   -0.000)   12.335   4.472   (  -7.488    8.366    0.000)   11.228   7.351   (   0.424   -7.444   -0.000)    7.456  14.182   (   8.959    1.458   -0.000)    9.077  14.786   (  -7.436   11.929    0.000)   14.057  15.446   (   1.471  -10.346   -0.000)   10.450======================= Grid point 1626 (217/220) =======================q-point: ( 0.63  0.42  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   3.921   (  -0.000   -0.817   13.566)   13.591   4.496   (  -0.000    9.568    1.860)    9.747   7.218   (  -0.000   -2.031    3.710)    4.230  13.928   (   0.000   -1.326   15.677)   15.733  14.790   (  -0.000   12.932    5.414)   14.020  15.515   (  -0.000   -9.385   -8.075)   12.381======================= Grid point 1627 (218/220) =======================q-point: ( 0.68  0.42  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 3430Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   4.020   (   5.867   -0.501   10.745)   12.253   4.590   (  -0.152    3.858    0.389)    3.880   7.227   (  -0.387   -0.987    0.830)    1.347  14.052   (   5.645   -0.481   12.226)   13.475  14.963   (   3.683    8.577    5.275)   10.722  15.339   (  -4.888   -7.145   -8.606)   12.207======================= Grid point 1646 (219/220) =======================q-point: ( 0.68  0.47  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   4.146   (  -0.000    2.517    5.403)    5.961   4.517   (   0.000    7.242    0.134)    7.244   7.272   (  -0.000   -4.855    1.080)    4.973  14.209   (   0.000    0.452    7.256)    7.270  14.873   (  -0.000   12.568    1.783)   12.694  15.357   (  -0.000  -10.816   -4.459)   11.699======================= Grid point 1647 (220/220) =======================q-point: (-0.26  0.47  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 5.98e-04 3.36e-05 Number of triplets: 1810Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   4.153   (   8.212    2.022   -0.000)    8.457   4.590   (  -2.854    1.176    0.000)    3.087   7.235   (   0.136   -1.944   -0.000)    1.948  14.207   (   8.518    0.607   -0.000)    8.540  15.049   (  -4.268    7.907    0.000)    8.985  15.202   (  -0.047   -7.949   -0.000)    7.949=================== 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/41154   10.0   1788.838   1788.838   1788.838      0.000     -0.000     -0.000 3/41154   20.0    972.959    972.959    972.959      0.000     -0.000     -0.000 3/41154   30.0    553.451    553.451    553.451      0.000     -0.000     -0.000 3/41154   40.0    409.226    409.226    409.226      0.000     -0.000     -0.000 3/41154   50.0    338.282    338.282    338.282      0.000     -0.000     -0.000 3/41154   60.0    294.611    294.611    294.611      0.000     -0.000     -0.000 3/41154   70.0    263.110    263.110    263.110      0.000     -0.000     -0.000 3/41154   80.0    237.827    237.827    237.827      0.000     -0.000     -0.000 3/41154   90.0    216.176    216.176    216.176      0.000     -0.000     -0.000 3/41154  100.0    197.036    197.036    197.036      0.000     -0.000     -0.000 3/41154  110.0    179.957    179.957    179.957      0.000     -0.000     -0.000 3/41154  120.0    164.746    164.746    164.746      0.000     -0.000     -0.000 3/41154  130.0    151.269    151.269    151.269      0.000     -0.000     -0.000 3/41154  140.0    139.380    139.380    139.380      0.000     -0.000     -0.000 3/41154  150.0    128.917    128.917    128.917      0.000     -0.000     -0.000 3/41154  160.0    119.711    119.711    119.711      0.000     -0.000     -0.000 3/41154  170.0    111.598    111.598    111.598      0.000     -0.000     -0.000 3/41154  180.0    104.430    104.430    104.430      0.000     -0.000     -0.000 3/41154  190.0     98.075     98.075     98.075      0.000     -0.000     -0.000 3/41154  200.0     92.418     92.418     92.418      0.000     -0.000     -0.000 3/41154  210.0     87.362     87.362     87.362      0.000     -0.000     -0.000 3/41154  220.0     82.823     82.823     82.823      0.000     -0.000     -0.000 3/41154  230.0     78.733     78.733     78.733      0.000     -0.000     -0.000 3/41154  240.0     75.030     75.030     75.030      0.000     -0.000     -0.000 3/41154  250.0     71.665     71.665     71.665      0.000     -0.000     -0.000 3/41154  260.0     68.595     68.595     68.595      0.000     -0.000     -0.000 3/41154  270.0     65.784     65.784     65.784      0.000     -0.000     -0.000 3/41154  280.0     63.202     63.202     63.202      0.000     -0.000     -0.000 3/41154  290.0     60.822     60.822     60.822      0.000     -0.000     -0.000 3/41154  300.0     58.621     58.621     58.621      0.000     -0.000     -0.000 3/41154  310.0     56.580     56.580     56.580      0.000     -0.000     -0.000 3/41154  320.0     54.682     54.682     54.682      0.000     -0.000     -0.000 3/41154  330.0     52.913     52.913     52.913      0.000     -0.000     -0.000 3/41154  340.0     51.259     51.259     51.259      0.000     -0.000     -0.000 3/41154  350.0     49.710     49.710     49.710      0.000     -0.000     -0.000 3/41154  360.0     48.257     48.257     48.257      0.000     -0.000     -0.000 3/41154  370.0     46.889     46.889     46.889      0.000     -0.000     -0.000 3/41154  380.0     45.601     45.601     45.601      0.000     -0.000     -0.000 3/41154  390.0     44.384     44.384     44.384      0.000     -0.000     -0.000 3/41154  400.0     43.233     43.233     43.233      0.000     -0.000     -0.000 3/41154  410.0     42.144     42.144     42.144      0.000     -0.000     -0.000 3/41154  420.0     41.110     41.110     41.110      0.000     -0.000     -0.000 3/41154  430.0     40.127     40.127     40.127      0.000     -0.000     -0.000 3/41154  440.0     39.193     39.193     39.193      0.000     -0.000     -0.000 3/41154  450.0     38.303     38.303     38.303      0.000     -0.000     -0.000 3/41154  460.0     37.454     37.454     37.454      0.000     -0.000     -0.000 3/41154  470.0     36.643     36.643     36.643      0.000     -0.000     -0.000 3/41154  480.0     35.868     35.868     35.868      0.000     -0.000     -0.000 3/41154  490.0     35.126     35.126     35.126      0.000     -0.000     -0.000 3/41154  500.0     34.415     34.415     34.415      0.000     -0.000     -0.000 3/41154  510.0     33.734     33.734     33.734      0.000     -0.000     -0.000 3/41154  520.0     33.080     33.080     33.080      0.000     -0.000     -0.000 3/41154  530.0     32.452     32.452     32.452      0.000     -0.000     -0.000 3/41154  540.0     31.848     31.848     31.848      0.000     -0.000     -0.000 3/41154  550.0     31.267     31.267     31.267      0.000     -0.000     -0.000 3/41154  560.0     30.707     30.707     30.707      0.000     -0.000     -0.000 3/41154  570.0     30.168     30.168     30.168      0.000     -0.000     -0.000 3/41154  580.0     29.648     29.648     29.648      0.000     -0.000     -0.000 3/41154  590.0     29.146     29.146     29.146      0.000     -0.000     -0.000 3/41154  600.0     28.661     28.661     28.661      0.000     -0.000     -0.000 3/41154  610.0     28.193     28.193     28.193      0.000     -0.000     -0.000 3/41154  620.0     27.740     27.740     27.740      0.000     -0.000     -0.000 3/41154  630.0     27.302     27.302     27.302      0.000     -0.000     -0.000 3/41154  640.0     26.878     26.878     26.878      0.000     -0.000     -0.000 3/41154  650.0     26.467     26.467     26.467      0.000     -0.000     -0.000 3/41154  660.0     26.069     26.069     26.069      0.000     -0.000     -0.000 3/41154  670.0     25.683     25.683     25.683      0.000     -0.000     -0.000 3/41154  680.0     25.308     25.308     25.308      0.000     -0.000     -0.000 3/41154  690.0     24.945     24.945     24.945      0.000     -0.000     -0.000 3/41154  700.0     24.592     24.592     24.592      0.000     -0.000     -0.000 3/41154  710.0     24.249     24.249     24.249      0.000     -0.000     -0.000 3/41154  720.0     23.916     23.916     23.916      0.000     -0.000     -0.000 3/41154  730.0     23.592     23.592     23.592      0.000     -0.000     -0.000 3/41154  740.0     23.276     23.276     23.276      0.000     -0.000     -0.000 3/41154  750.0     22.970     22.970     22.970      0.000     -0.000     -0.000 3/41154  760.0     22.671     22.671     22.671      0.000     -0.000     -0.000 3/41154  770.0     22.381     22.381     22.381      0.000     -0.000     -0.000 3/41154  780.0     22.097     22.097     22.097      0.000     -0.000     -0.000 3/41154  790.0     21.821     21.821     21.821      0.000     -0.000     -0.000 3/41154  800.0     21.552     21.552     21.552      0.000     -0.000     -0.000 3/41154  810.0     21.290     21.290     21.290      0.000     -0.000     -0.000 3/41154  820.0     21.034     21.034     21.034      0.000     -0.000     -0.000 3/41154  830.0     20.784     20.784     20.784      0.000     -0.000     -0.000 3/41154  840.0     20.541     20.541     20.541      0.000     -0.000     -0.000 3/41154  850.0     20.303     20.303     20.303      0.000     -0.000     -0.000 3/41154  860.0     20.070     20.070     20.070      0.000     -0.000     -0.000 3/41154  870.0     19.843     19.843     19.843      0.000     -0.000     -0.000 3/41154  880.0     19.621     19.621     19.621      0.000     -0.000     -0.000 3/41154  890.0     19.404     19.404     19.404      0.000     -0.000     -0.000 3/41154  900.0     19.192     19.192     19.192      0.000     -0.000     -0.000 3/41154  910.0     18.985     18.985     18.985      0.000     -0.000     -0.000 3/41154  920.0     18.782     18.782     18.782      0.000     -0.000     -0.000 3/41154  930.0     18.583     18.583     18.583      0.000     -0.000     -0.000 3/41154  940.0     18.389     18.389     18.389      0.000     -0.000     -0.000 3/41154  950.0     18.198     18.198     18.198      0.000     -0.000     -0.000 3/41154  960.0     18.012     18.012     18.012      0.000     -0.000     -0.000 3/41154  970.0     17.829     17.829     17.829      0.000     -0.000     -0.000 3/41154  980.0     17.651     17.651     17.651      0.000     -0.000     -0.000 3/41154  990.0     17.475     17.475     17.475      0.000     -0.000     -0.000 3/41154 1000.0     17.304     17.304     17.304      0.000     -0.000     -0.000 3/41154Thermal conductivity related properties were written into "kappa-m191919.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-07 19:03:21]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|