# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/d52c6d5d-a642-44ab-b962-f685561c72f1/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 KNa2CuO2 / I4mm (107) / materials id 545359](https://mdr.nims.go.jp/datasets/44455db8-a1f4-412b-a0d2-42b3868cb41d)

## Fulltext

------------------------------------ calculate fc2 ------------------------------------        _  _ __ | |__   ___  _ __   ___   _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_|                            |_|    |___/                                      2.47.1-------------------------[time 2026-01-08 10:41: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: [3 3 1]  Primitive matrix:    [-0.5  0.5  0.5]    [ 0.5 -0.5  0.5]    [ 0.5  0.5 -0.5]Spacegroup: I4mm (107)Number of symmetry operations in supercell: 144------------------------------ primitive cell ------------------------------Lattice vectors:  a   -2.188854930000001    2.188854930000001    5.256536895000000  b    2.188854930000001   -2.188854930000001    5.256536895000000  c    2.188854930000001    2.188854930000001   -5.256536895000000Atomic positions (fractional):   *1 K   0.56334326478607  0.56334326478607  0.00000000000000  39.098   *2 O   0.16918583876024  0.16918583876024  0.00000000000000  15.999   *3 O   0.82931053362736  0.82931053362736  0.00000000000000  15.999   *4 Na  0.77203019469644  0.27203019469644  0.50000000000000  22.990    5 Na  0.27203019469644  0.77203019469644  0.50000000000000  22.990   *6 Cu  0.99903880343342  0.99903880343342  0.00000000000000  63.546-------------------------------- unit cell ---------------------------------Lattice vectors:  a    4.377709860000001    0.000000000000000    0.000000000000000  b    0.000000000000000    4.377709860000001    0.000000000000000  c    0.000000000000000    0.000000000000000   10.513073790000000Atomic positions (fractional):   *1 K   0.50000000000000  0.50000000000000  0.06334326478607  39.098 > 1    2 K   0.00000000000000  0.00000000000000  0.56334326478607  39.098 > 1   *3 O   0.00000000000000  0.00000000000000  0.16918583876024  15.999 > 2   *4 O   0.50000000000000  0.50000000000000  0.32931053362736  15.999 > 3    5 O   0.50000000000000  0.50000000000000  0.66918583876024  15.999 > 2    6 O   0.00000000000000  0.00000000000000  0.82931053362736  15.999 > 3   *7 Na  0.00000000000000  0.50000000000000  0.27203019469644  22.990 > 4    8 Na  0.50000000000000  0.00000000000000  0.27203019469644  22.990 > 5    9 Na  0.50000000000000  0.00000000000000  0.77203019469644  22.990 > 4   10 Na  0.00000000000000  0.50000000000000  0.77203019469644  22.990 > 5  *11 Cu  0.50000000000000  0.50000000000000  0.49903880343342  63.546 > 6   12 Cu  0.00000000000000  0.00000000000000  0.99903880343342  63.546 > 6-------------------------------- super cell --------------------------------Lattice vectors:  a   13.133129580000004    0.000000000000000    0.000000000000000  b    0.000000000000000   13.133129580000004    0.000000000000000  c    0.000000000000000    0.000000000000000   10.513073790000000Atomic positions (fractional):   *1 K   0.16666666666667  0.16666666666667  0.06334326478607  39.098 > 1    2 K   0.50000000000000  0.16666666666667  0.06334326478607  39.098 > 1    3 K   0.83333333333333  0.16666666666667  0.06334326478607  39.098 > 1    4 K   0.16666666666667  0.50000000000000  0.06334326478607  39.098 > 1    5 K   0.50000000000000  0.50000000000000  0.06334326478607  39.098 > 1    6 K   0.83333333333333  0.50000000000000  0.06334326478607  39.098 > 1    7 K   0.16666666666667  0.83333333333333  0.06334326478607  39.098 > 1    8 K   0.50000000000000  0.83333333333333  0.06334326478607  39.098 > 1    9 K   0.83333333333333  0.83333333333333  0.06334326478607  39.098 > 1   10 K   0.00000000000000  0.00000000000000  0.56334326478607  39.098 > 1   11 K   0.33333333333333  0.00000000000000  0.56334326478607  39.098 > 1   12 K   0.66666666666667  0.00000000000000  0.56334326478607  39.098 > 1   13 K   0.00000000000000  0.33333333333333  0.56334326478607  39.098 > 1   14 K   0.33333333333333  0.33333333333333  0.56334326478607  39.098 > 1   15 K   0.66666666666667  0.33333333333333  0.56334326478607  39.098 > 1   16 K   0.00000000000000  0.66666666666667  0.56334326478607  39.098 > 1   17 K   0.33333333333333  0.66666666666667  0.56334326478607  39.098 > 1   18 K   0.66666666666667  0.66666666666667  0.56334326478607  39.098 > 1  *19 O   0.00000000000000  0.00000000000000  0.16918583876024  15.999 > 2   20 O   0.33333333333333  0.00000000000000  0.16918583876024  15.999 > 2   21 O   0.66666666666667  0.00000000000000  0.16918583876024  15.999 > 2   22 O   0.00000000000000  0.33333333333333  0.16918583876024  15.999 > 2   23 O   0.33333333333333  0.33333333333333  0.16918583876024  15.999 > 2   24 O   0.66666666666667  0.33333333333333  0.16918583876024  15.999 > 2   25 O   0.00000000000000  0.66666666666667  0.16918583876024  15.999 > 2   26 O   0.33333333333333  0.66666666666667  0.16918583876024  15.999 > 2   27 O   0.66666666666667  0.66666666666667  0.16918583876024  15.999 > 2  *28 O   0.16666666666667  0.16666666666667  0.32931053362736  15.999 > 3   29 O   0.50000000000000  0.16666666666667  0.32931053362736  15.999 > 3   30 O   0.83333333333333  0.16666666666667  0.32931053362736  15.999 > 3   31 O   0.16666666666667  0.50000000000000  0.32931053362736  15.999 > 3   32 O   0.50000000000000  0.50000000000000  0.32931053362736  15.999 > 3   33 O   0.83333333333333  0.50000000000000  0.32931053362736  15.999 > 3   34 O   0.16666666666667  0.83333333333333  0.32931053362736  15.999 > 3   35 O   0.50000000000000  0.83333333333333  0.32931053362736  15.999 > 3   36 O   0.83333333333333  0.83333333333333  0.32931053362736  15.999 > 3   37 O   0.16666666666667  0.16666666666667  0.66918583876024  15.999 > 2   38 O   0.50000000000000  0.16666666666667  0.66918583876024  15.999 > 2   39 O   0.83333333333333  0.16666666666667  0.66918583876024  15.999 > 2   40 O   0.16666666666667  0.50000000000000  0.66918583876024  15.999 > 2   41 O   0.50000000000000  0.50000000000000  0.66918583876024  15.999 > 2   42 O   0.83333333333333  0.50000000000000  0.66918583876024  15.999 > 2   43 O   0.16666666666667  0.83333333333333  0.66918583876024  15.999 > 2   44 O   0.50000000000000  0.83333333333333  0.66918583876024  15.999 > 2   45 O   0.83333333333333  0.83333333333333  0.66918583876024  15.999 > 2   46 O   0.00000000000000  0.00000000000000  0.82931053362736  15.999 > 3   47 O   0.33333333333333  0.00000000000000  0.82931053362736  15.999 > 3   48 O   0.66666666666667  0.00000000000000  0.82931053362736  15.999 > 3   49 O   0.00000000000000  0.33333333333333  0.82931053362736  15.999 > 3   50 O   0.33333333333333  0.33333333333333  0.82931053362736  15.999 > 3   51 O   0.66666666666667  0.33333333333333  0.82931053362736  15.999 > 3   52 O   0.00000000000000  0.66666666666667  0.82931053362736  15.999 > 3   53 O   0.33333333333333  0.66666666666667  0.82931053362736  15.999 > 3   54 O   0.66666666666667  0.66666666666667  0.82931053362736  15.999 > 3  *55 Na  0.00000000000000  0.16666666666667  0.27203019469644  22.990 > 4   56 Na  0.33333333333333  0.16666666666667  0.27203019469644  22.990 > 4   57 Na  0.66666666666667  0.16666666666667  0.27203019469644  22.990 > 4   58 Na  0.00000000000000  0.50000000000000  0.27203019469644  22.990 > 4   59 Na  0.33333333333333  0.50000000000000  0.27203019469644  22.990 > 4   60 Na  0.66666666666667  0.50000000000000  0.27203019469644  22.990 > 4   61 Na  0.00000000000000  0.83333333333333  0.27203019469644  22.990 > 4   62 Na  0.33333333333333  0.83333333333333  0.27203019469644  22.990 > 4   63 Na  0.66666666666667  0.83333333333333  0.27203019469644  22.990 > 4   64 Na  0.16666666666667  0.00000000000000  0.27203019469644  22.990 > 5   65 Na  0.50000000000000  0.00000000000000  0.27203019469644  22.990 > 5   66 Na  0.83333333333333  0.00000000000000  0.27203019469644  22.990 > 5   67 Na  0.16666666666667  0.33333333333333  0.27203019469644  22.990 > 5   68 Na  0.50000000000000  0.33333333333333  0.27203019469644  22.990 > 5   69 Na  0.83333333333333  0.33333333333333  0.27203019469644  22.990 > 5   70 Na  0.16666666666667  0.66666666666667  0.27203019469644  22.990 > 5   71 Na  0.50000000000000  0.66666666666667  0.27203019469644  22.990 > 5   72 Na  0.83333333333333  0.66666666666667  0.27203019469644  22.990 > 5   73 Na  0.16666666666667  0.00000000000000  0.77203019469644  22.990 > 4   74 Na  0.50000000000000  0.00000000000000  0.77203019469644  22.990 > 4   75 Na  0.83333333333333  0.00000000000000  0.77203019469644  22.990 > 4   76 Na  0.16666666666667  0.33333333333333  0.77203019469644  22.990 > 4   77 Na  0.50000000000000  0.33333333333333  0.77203019469644  22.990 > 4   78 Na  0.83333333333333  0.33333333333333  0.77203019469644  22.990 > 4   79 Na  0.16666666666667  0.66666666666667  0.77203019469644  22.990 > 4   80 Na  0.50000000000000  0.66666666666667  0.77203019469644  22.990 > 4   81 Na  0.83333333333333  0.66666666666667  0.77203019469644  22.990 > 4   82 Na  0.00000000000000  0.16666666666667  0.77203019469644  22.990 > 5   83 Na  0.33333333333333  0.16666666666667  0.77203019469644  22.990 > 5   84 Na  0.66666666666667  0.16666666666667  0.77203019469644  22.990 > 5   85 Na  0.00000000000000  0.50000000000000  0.77203019469644  22.990 > 5   86 Na  0.33333333333333  0.50000000000000  0.77203019469644  22.990 > 5   87 Na  0.66666666666667  0.50000000000000  0.77203019469644  22.990 > 5   88 Na  0.00000000000000  0.83333333333333  0.77203019469644  22.990 > 5   89 Na  0.33333333333333  0.83333333333333  0.77203019469644  22.990 > 5   90 Na  0.66666666666667  0.83333333333333  0.77203019469644  22.990 > 5  *91 Cu  0.16666666666667  0.16666666666667  0.49903880343342  63.546 > 6   92 Cu  0.50000000000000  0.16666666666667  0.49903880343342  63.546 > 6   93 Cu  0.83333333333333  0.16666666666667  0.49903880343342  63.546 > 6   94 Cu  0.16666666666667  0.50000000000000  0.49903880343342  63.546 > 6   95 Cu  0.50000000000000  0.50000000000000  0.49903880343342  63.546 > 6   96 Cu  0.83333333333333  0.50000000000000  0.49903880343342  63.546 > 6   97 Cu  0.16666666666667  0.83333333333333  0.49903880343342  63.546 > 6   98 Cu  0.50000000000000  0.83333333333333  0.49903880343342  63.546 > 6   99 Cu  0.83333333333333  0.83333333333333  0.49903880343342  63.546 > 6  100 Cu  0.00000000000000  0.00000000000000  0.99903880343342  63.546 > 6  101 Cu  0.33333333333333  0.00000000000000  0.99903880343342  63.546 > 6  102 Cu  0.66666666666667  0.00000000000000  0.99903880343342  63.546 > 6  103 Cu  0.00000000000000  0.33333333333333  0.99903880343342  63.546 > 6  104 Cu  0.33333333333333  0.33333333333333  0.99903880343342  63.546 > 6  105 Cu  0.66666666666667  0.33333333333333  0.99903880343342  63.546 > 6  106 Cu  0.00000000000000  0.66666666666667  0.99903880343342  63.546 > 6  107 Cu  0.33333333333333  0.66666666666667  0.99903880343342  63.546 > 6  108 Cu  0.66666666666667  0.66666666666667  0.99903880343342  63.546 > 6----------------------------------------------------------------------------NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".--------------------------- Dielectric constant ----------------------------            3.3400940    0.0000000    0.0000000            0.0000000    3.3400940    0.0000000            0.0000000    0.0000000    3.4803409-------------------------- Born effective charges --------------------------    1 K     1.1098283    0.0000000    0.0000000            0.0000000    1.1098283    0.0000000            0.0000000    0.0000000    1.1499791    2 O    -1.2017887    0.0000000    0.0000000            0.0000000   -1.2017887    0.0000000            0.0000000    0.0000000   -2.3213243    3 O    -0.9895916    0.0000000    0.0000000            0.0000000   -0.9895916    0.0000000            0.0000000    0.0000000   -2.5544378    4 Na    1.0151837    0.0000000    0.0000000            0.0000000    0.9405108    0.0000000            0.0000000    0.0000000    1.0156978    5 Na    0.9405108    0.0000000    0.0000000            0.0000000    1.0151837    0.0000000            0.0000000    0.0000000    1.0156978    6 Cu   -0.8741425    0.0000000    0.0000000            0.0000000   -0.8741425    0.0000000            0.0000000    0.0000000    1.6943874----------------------------------------------------------------------------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: 324/324Permutation basis: 3528/3528Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 460Number of blocks in projector: 460Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 12--- Eigsh_solver_block: 1 / 12 ---Block_size: 168Use standard eigh solver.--- Eigsh_solver_block: 2 / 12 ---Block_size: 144Use standard eigh solver.--- Eigsh_solver_block: 3 / 12 ---Block_size: 106Use standard eigh solver.--- Eigsh_solver_block: 4 / 12 ---Block_size: 34Use standard eigh solver.--- Eigsh_solver_block: 5 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 6 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 7 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 8 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 9 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 10 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 11 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 12 / 12 ---Block_size: 1Use standard eigh solver.Tree of FC basis block matrices:- (460, 449), data: False|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (34, 34), data: True|-- (106, 101), data: True|-- (144, 144), data: True|-- (168, 162), data: True-----Solver_atoms: 1 -- 108 / 108Time (Solver_compr_matrix_reshape): 0.003Solver_block: 100 / 160 - Time: 0.334Solver_block: 160 / 160 - Time: 0.246Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.588--------------------------------- Symfc end --------------------------------Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 324/324Permutation basis: 3528/3528Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 460Number of blocks in projector: 460Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 12--- Eigsh_solver_block: 1 / 12 ---Block_size: 168Use standard eigh solver.--- Eigsh_solver_block: 2 / 12 ---Block_size: 144Use standard eigh solver.--- Eigsh_solver_block: 3 / 12 ---Block_size: 106Use standard eigh solver.--- Eigsh_solver_block: 4 / 12 ---Block_size: 34Use standard eigh solver.--- Eigsh_solver_block: 5 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 6 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 7 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 8 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 9 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 10 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 11 / 12 ---Block_size: 1Use standard eigh solver.--- Eigsh_solver_block: 12 / 12 ---Block_size: 1Use standard eigh solver.Tree of FC basis block matrices:- (460, 449), data: False|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (1, 1), data: True|-- (34, 34), data: True|-- (106, 101), data: True|-- (144, 144), data: True|-- (168, 162), data: TrueMax drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) Force constants are written into "force_constants.hdf5".---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d)----------------------------------------------------------------------------Summary of calculation was written in "phonopy.yaml".-------------------------[time 2026-01-08 10:41:56]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate fc3 -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 10:41:56]-------------------------Compiled with OpenMP support (max 128 threads).Running in phono3py.load mode.Python version 3.14.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: force constantsHDF5 data compression filter: gzipCrystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".Supercell (dim): [3 3 1]Primitive matrix:  [-0.5  0.5  0.5]  [ 0.5 -0.5  0.5]  [ 0.5  0.5 -0.5]Spacegroup: I4mm (107)------------------------------ primitive cell ------------------------------Lattice vectors:  a   -2.188854930000001    2.188854930000001    5.256536895000000  b    2.188854930000001   -2.188854930000001    5.256536895000000  c    2.188854930000001    2.188854930000001   -5.256536895000000Atomic positions (fractional):    1 K   0.56334326478607  0.56334326478607  0.00000000000000  39.098    2 O   0.16918583876024  0.16918583876024  0.00000000000000  15.999    3 O   0.82931053362736  0.82931053362736  0.00000000000000  15.999    4 Na  0.77203019469644  0.27203019469644  0.50000000000000  22.990    5 Na  0.27203019469644  0.77203019469644  0.50000000000000  22.990    6 Cu  0.99903880343342  0.99903880343342  0.00000000000000  63.546-------------------------------- supercell ---------------------------------Lattice vectors:  a   13.133129580000004    0.000000000000000    0.000000000000000  b    0.000000000000000   13.133129580000004    0.000000000000000  c    0.000000000000000    0.000000000000000   10.513073790000000Atomic positions (fractional):    1 K   0.16666666666667  0.16666666666667  0.06334326478607  39.098 > 1    2 K   0.50000000000000  0.16666666666667  0.06334326478607  39.098 > 1    3 K   0.83333333333333  0.16666666666667  0.06334326478607  39.098 > 1    4 K   0.16666666666667  0.50000000000000  0.06334326478607  39.098 > 1    5 K   0.50000000000000  0.50000000000000  0.06334326478607  39.098 > 1    6 K   0.83333333333333  0.50000000000000  0.06334326478607  39.098 > 1    7 K   0.16666666666667  0.83333333333333  0.06334326478607  39.098 > 1    8 K   0.50000000000000  0.83333333333333  0.06334326478607  39.098 > 1    9 K   0.83333333333333  0.83333333333333  0.06334326478607  39.098 > 1   10 K   0.00000000000000  0.00000000000000  0.56334326478607  39.098 > 1   11 K   0.33333333333333  0.00000000000000  0.56334326478607  39.098 > 1   12 K   0.66666666666667  0.00000000000000  0.56334326478607  39.098 > 1   13 K   0.00000000000000  0.33333333333333  0.56334326478607  39.098 > 1   14 K   0.33333333333333  0.33333333333333  0.56334326478607  39.098 > 1   15 K   0.66666666666667  0.33333333333333  0.56334326478607  39.098 > 1   16 K   0.00000000000000  0.66666666666667  0.56334326478607  39.098 > 1   17 K   0.33333333333333  0.66666666666667  0.56334326478607  39.098 > 1   18 K   0.66666666666667  0.66666666666667  0.56334326478607  39.098 > 1   19 O   0.00000000000000  0.00000000000000  0.16918583876024  15.999 > 19   20 O   0.33333333333333  0.00000000000000  0.16918583876024  15.999 > 19   21 O   0.66666666666667  0.00000000000000  0.16918583876024  15.999 > 19   22 O   0.00000000000000  0.33333333333333  0.16918583876024  15.999 > 19   23 O   0.33333333333333  0.33333333333333  0.16918583876024  15.999 > 19   24 O   0.66666666666667  0.33333333333333  0.16918583876024  15.999 > 19   25 O   0.00000000000000  0.66666666666667  0.16918583876024  15.999 > 19   26 O   0.33333333333333  0.66666666666667  0.16918583876024  15.999 > 19   27 O   0.66666666666667  0.66666666666667  0.16918583876024  15.999 > 19   28 O   0.16666666666667  0.16666666666667  0.32931053362736  15.999 > 28   29 O   0.50000000000000  0.16666666666667  0.32931053362736  15.999 > 28   30 O   0.83333333333333  0.16666666666667  0.32931053362736  15.999 > 28   31 O   0.16666666666667  0.50000000000000  0.32931053362736  15.999 > 28   32 O   0.50000000000000  0.50000000000000  0.32931053362736  15.999 > 28   33 O   0.83333333333333  0.50000000000000  0.32931053362736  15.999 > 28   34 O   0.16666666666667  0.83333333333333  0.32931053362736  15.999 > 28   35 O   0.50000000000000  0.83333333333333  0.32931053362736  15.999 > 28   36 O   0.83333333333333  0.83333333333333  0.32931053362736  15.999 > 28   37 O   0.16666666666667  0.16666666666667  0.66918583876024  15.999 > 19   38 O   0.50000000000000  0.16666666666667  0.66918583876024  15.999 > 19   39 O   0.83333333333333  0.16666666666667  0.66918583876024  15.999 > 19   40 O   0.16666666666667  0.50000000000000  0.66918583876024  15.999 > 19   41 O   0.50000000000000  0.50000000000000  0.66918583876024  15.999 > 19   42 O   0.83333333333333  0.50000000000000  0.66918583876024  15.999 > 19   43 O   0.16666666666667  0.83333333333333  0.66918583876024  15.999 > 19   44 O   0.50000000000000  0.83333333333333  0.66918583876024  15.999 > 19   45 O   0.83333333333333  0.83333333333333  0.66918583876024  15.999 > 19   46 O   0.00000000000000  0.00000000000000  0.82931053362736  15.999 > 28   47 O   0.33333333333333  0.00000000000000  0.82931053362736  15.999 > 28   48 O   0.66666666666667  0.00000000000000  0.82931053362736  15.999 > 28   49 O   0.00000000000000  0.33333333333333  0.82931053362736  15.999 > 28   50 O   0.33333333333333  0.33333333333333  0.82931053362736  15.999 > 28   51 O   0.66666666666667  0.33333333333333  0.82931053362736  15.999 > 28   52 O   0.00000000000000  0.66666666666667  0.82931053362736  15.999 > 28   53 O   0.33333333333333  0.66666666666667  0.82931053362736  15.999 > 28   54 O   0.66666666666667  0.66666666666667  0.82931053362736  15.999 > 28   55 Na  0.00000000000000  0.16666666666667  0.27203019469644  22.990 > 55   56 Na  0.33333333333333  0.16666666666667  0.27203019469644  22.990 > 55   57 Na  0.66666666666667  0.16666666666667  0.27203019469644  22.990 > 55   58 Na  0.00000000000000  0.50000000000000  0.27203019469644  22.990 > 55   59 Na  0.33333333333333  0.50000000000000  0.27203019469644  22.990 > 55   60 Na  0.66666666666667  0.50000000000000  0.27203019469644  22.990 > 55   61 Na  0.00000000000000  0.83333333333333  0.27203019469644  22.990 > 55   62 Na  0.33333333333333  0.83333333333333  0.27203019469644  22.990 > 55   63 Na  0.66666666666667  0.83333333333333  0.27203019469644  22.990 > 55   64 Na  0.16666666666667  0.00000000000000  0.27203019469644  22.990 > 64   65 Na  0.50000000000000  0.00000000000000  0.27203019469644  22.990 > 64   66 Na  0.83333333333333  0.00000000000000  0.27203019469644  22.990 > 64   67 Na  0.16666666666667  0.33333333333333  0.27203019469644  22.990 > 64   68 Na  0.50000000000000  0.33333333333333  0.27203019469644  22.990 > 64   69 Na  0.83333333333333  0.33333333333333  0.27203019469644  22.990 > 64   70 Na  0.16666666666667  0.66666666666667  0.27203019469644  22.990 > 64   71 Na  0.50000000000000  0.66666666666667  0.27203019469644  22.990 > 64   72 Na  0.83333333333333  0.66666666666667  0.27203019469644  22.990 > 64   73 Na  0.16666666666667  0.00000000000000  0.77203019469644  22.990 > 55   74 Na  0.50000000000000  0.00000000000000  0.77203019469644  22.990 > 55   75 Na  0.83333333333333  0.00000000000000  0.77203019469644  22.990 > 55   76 Na  0.16666666666667  0.33333333333333  0.77203019469644  22.990 > 55   77 Na  0.50000000000000  0.33333333333333  0.77203019469644  22.990 > 55   78 Na  0.83333333333333  0.33333333333333  0.77203019469644  22.990 > 55   79 Na  0.16666666666667  0.66666666666667  0.77203019469644  22.990 > 55   80 Na  0.50000000000000  0.66666666666667  0.77203019469644  22.990 > 55   81 Na  0.83333333333333  0.66666666666667  0.77203019469644  22.990 > 55   82 Na  0.00000000000000  0.16666666666667  0.77203019469644  22.990 > 64   83 Na  0.33333333333333  0.16666666666667  0.77203019469644  22.990 > 64   84 Na  0.66666666666667  0.16666666666667  0.77203019469644  22.990 > 64   85 Na  0.00000000000000  0.50000000000000  0.77203019469644  22.990 > 64   86 Na  0.33333333333333  0.50000000000000  0.77203019469644  22.990 > 64   87 Na  0.66666666666667  0.50000000000000  0.77203019469644  22.990 > 64   88 Na  0.00000000000000  0.83333333333333  0.77203019469644  22.990 > 64   89 Na  0.33333333333333  0.83333333333333  0.77203019469644  22.990 > 64   90 Na  0.66666666666667  0.83333333333333  0.77203019469644  22.990 > 64   91 Cu  0.16666666666667  0.16666666666667  0.49903880343342  63.546 > 91   92 Cu  0.50000000000000  0.16666666666667  0.49903880343342  63.546 > 91   93 Cu  0.83333333333333  0.16666666666667  0.49903880343342  63.546 > 91   94 Cu  0.16666666666667  0.50000000000000  0.49903880343342  63.546 > 91   95 Cu  0.50000000000000  0.50000000000000  0.49903880343342  63.546 > 91   96 Cu  0.83333333333333  0.50000000000000  0.49903880343342  63.546 > 91   97 Cu  0.16666666666667  0.83333333333333  0.49903880343342  63.546 > 91   98 Cu  0.50000000000000  0.83333333333333  0.49903880343342  63.546 > 91   99 Cu  0.83333333333333  0.83333333333333  0.49903880343342  63.546 > 91  100 Cu  0.00000000000000  0.00000000000000  0.99903880343342  63.546 > 91  101 Cu  0.33333333333333  0.00000000000000  0.99903880343342  63.546 > 91  102 Cu  0.66666666666667  0.00000000000000  0.99903880343342  63.546 > 91  103 Cu  0.00000000000000  0.33333333333333  0.99903880343342  63.546 > 91  104 Cu  0.33333333333333  0.33333333333333  0.99903880343342  63.546 > 91  105 Cu  0.66666666666667  0.33333333333333  0.99903880343342  63.546 > 91  106 Cu  0.00000000000000  0.66666666666667  0.99903880343342  63.546 > 91  107 Cu  0.33333333333333  0.66666666666667  0.99903880343342  63.546 > 91  108 Cu  0.66666666666667  0.66666666666667  0.99903880343342  63.546 > 91----------------------------------------------------------------------------NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz".--------------------------- Dielectric constant ----------------------------            3.3400940    0.0000000    0.0000000            0.0000000    3.3400940    0.0000000            0.0000000    0.0000000    3.4803409-------------------------- Born effective charges --------------------------    1 K     1.1098283    0.0000000    0.0000000            0.0000000    1.1098283    0.0000000            0.0000000    0.0000000    1.1499791    2 O    -1.2017887    0.0000000    0.0000000            0.0000000   -1.2017887    0.0000000            0.0000000    0.0000000   -2.3213243    3 O    -0.9895916    0.0000000    0.0000000            0.0000000   -0.9895916    0.0000000            0.0000000    0.0000000   -2.5544378    4 Na    1.0151837    0.0000000    0.0000000            0.0000000    0.9405108    0.0000000            0.0000000    0.0000000    1.0156978    5 Na    0.9405108    0.0000000    0.0000000            0.0000000    1.0151837    0.0000000            0.0000000    0.0000000    1.0156978    6 Cu   -0.8741425    0.0000000    0.0000000            0.0000000   -0.8741425    0.0000000            0.0000000    0.0000000    1.6943874----------------------------------------------------------------------------Sets of supercell forces were read from "FORCES_FC3.xz".Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz".----------------------------- Force constants ------------------------------Computing fc3[ 1, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 19, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 28, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 55, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0100  0.0000]    [ 0.0000 -0.0100  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Computing fc3[ 91, x, x ] using numpy.linalg.pinv.Displacements (in Angstrom):    [ 0.0100  0.0000  0.0000]    [-0.0100  0.0000  0.0000]    [ 0.0000  0.0000  0.0100]    [ 0.0000  0.0000 -0.0100]Expanding fc3.Symmetrizing fc3 by traditional approach (N=3).Symmetrizing fc2 by traditional approach (N=3).Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz)fc3 was written into "fc3.hdf5".Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) fc2 was written into "fc2.hdf5".--------------------------- Calculation settings ---------------------------Non-analytical term correction (NAC): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330----------- None of ph-ph interaction calculation was performed. -----------Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 10:42:03]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|---------------------------------- calculate LTC -----------------------------------        _                      _____  _ __ | |__   ___  _ __   ___|___ / _ __  _   _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_|                                |_|    |___/                                       3.23.0-------------------------[time 2026-01-08 10:42:04]-------------------------Compiled with OpenMP support (max 128 threads).Running in phono3py.load mode.Python version 3.14.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: conductivity-RTAHDF5 data compression filter: gzipCrystal structure was read from "phono3py.yaml".Supercell (dim): [3 3 1]Primitive matrix:  [-0.5  0.5  0.5]  [ 0.5 -0.5  0.5]  [ 0.5  0.5 -0.5]Spacegroup: I4mm (107)------------------------------ primitive cell ------------------------------Lattice vectors:  a   -2.188854930000001    2.188854930000001    5.256536895000000  b    2.188854930000001   -2.188854930000001    5.256536895000000  c    2.188854930000001    2.188854930000001   -5.256536895000000Atomic positions (fractional):    1 K   0.56334326478607  0.56334326478607  0.00000000000000  39.098    2 O   0.16918583876024  0.16918583876024  0.00000000000000  15.999    3 O   0.82931053362736  0.82931053362736  0.00000000000000  15.999    4 Na  0.77203019469644  0.27203019469644  0.50000000000000  22.990    5 Na  0.27203019469644  0.77203019469644  0.50000000000000  22.990    6 Cu  0.99903880343342  0.99903880343342  0.00000000000000  63.546-------------------------------- supercell ---------------------------------Lattice vectors:  a   13.133129580000004    0.000000000000000    0.000000000000000  b    0.000000000000000   13.133129580000004    0.000000000000000  c    0.000000000000000    0.000000000000000   10.513073790000000Atomic positions (fractional):    1 K   0.16666666666667  0.16666666666667  0.06334326478607  39.098 > 1    2 K   0.50000000000000  0.16666666666667  0.06334326478607  39.098 > 1    3 K   0.83333333333333  0.16666666666667  0.06334326478607  39.098 > 1    4 K   0.16666666666667  0.50000000000000  0.06334326478607  39.098 > 1    5 K   0.50000000000000  0.50000000000000  0.06334326478607  39.098 > 1    6 K   0.83333333333333  0.50000000000000  0.06334326478607  39.098 > 1    7 K   0.16666666666667  0.83333333333333  0.06334326478607  39.098 > 1    8 K   0.50000000000000  0.83333333333333  0.06334326478607  39.098 > 1    9 K   0.83333333333333  0.83333333333333  0.06334326478607  39.098 > 1   10 K   0.00000000000000  0.00000000000000  0.56334326478607  39.098 > 1   11 K   0.33333333333333  0.00000000000000  0.56334326478607  39.098 > 1   12 K   0.66666666666667  0.00000000000000  0.56334326478607  39.098 > 1   13 K   0.00000000000000  0.33333333333333  0.56334326478607  39.098 > 1   14 K   0.33333333333333  0.33333333333333  0.56334326478607  39.098 > 1   15 K   0.66666666666667  0.33333333333333  0.56334326478607  39.098 > 1   16 K   0.00000000000000  0.66666666666667  0.56334326478607  39.098 > 1   17 K   0.33333333333333  0.66666666666667  0.56334326478607  39.098 > 1   18 K   0.66666666666667  0.66666666666667  0.56334326478607  39.098 > 1   19 O   0.00000000000000  0.00000000000000  0.16918583876024  15.999 > 19   20 O   0.33333333333333  0.00000000000000  0.16918583876024  15.999 > 19   21 O   0.66666666666667  0.00000000000000  0.16918583876024  15.999 > 19   22 O   0.00000000000000  0.33333333333333  0.16918583876024  15.999 > 19   23 O   0.33333333333333  0.33333333333333  0.16918583876024  15.999 > 19   24 O   0.66666666666667  0.33333333333333  0.16918583876024  15.999 > 19   25 O   0.00000000000000  0.66666666666667  0.16918583876024  15.999 > 19   26 O   0.33333333333333  0.66666666666667  0.16918583876024  15.999 > 19   27 O   0.66666666666667  0.66666666666667  0.16918583876024  15.999 > 19   28 O   0.16666666666667  0.16666666666667  0.32931053362736  15.999 > 28   29 O   0.50000000000000  0.16666666666667  0.32931053362736  15.999 > 28   30 O   0.83333333333333  0.16666666666667  0.32931053362736  15.999 > 28   31 O   0.16666666666667  0.50000000000000  0.32931053362736  15.999 > 28   32 O   0.50000000000000  0.50000000000000  0.32931053362736  15.999 > 28   33 O   0.83333333333333  0.50000000000000  0.32931053362736  15.999 > 28   34 O   0.16666666666667  0.83333333333333  0.32931053362736  15.999 > 28   35 O   0.50000000000000  0.83333333333333  0.32931053362736  15.999 > 28   36 O   0.83333333333333  0.83333333333333  0.32931053362736  15.999 > 28   37 O   0.16666666666667  0.16666666666667  0.66918583876024  15.999 > 19   38 O   0.50000000000000  0.16666666666667  0.66918583876024  15.999 > 19   39 O   0.83333333333333  0.16666666666667  0.66918583876024  15.999 > 19   40 O   0.16666666666667  0.50000000000000  0.66918583876024  15.999 > 19   41 O   0.50000000000000  0.50000000000000  0.66918583876024  15.999 > 19   42 O   0.83333333333333  0.50000000000000  0.66918583876024  15.999 > 19   43 O   0.16666666666667  0.83333333333333  0.66918583876024  15.999 > 19   44 O   0.50000000000000  0.83333333333333  0.66918583876024  15.999 > 19   45 O   0.83333333333333  0.83333333333333  0.66918583876024  15.999 > 19   46 O   0.00000000000000  0.00000000000000  0.82931053362736  15.999 > 28   47 O   0.33333333333333  0.00000000000000  0.82931053362736  15.999 > 28   48 O   0.66666666666667  0.00000000000000  0.82931053362736  15.999 > 28   49 O   0.00000000000000  0.33333333333333  0.82931053362736  15.999 > 28   50 O   0.33333333333333  0.33333333333333  0.82931053362736  15.999 > 28   51 O   0.66666666666667  0.33333333333333  0.82931053362736  15.999 > 28   52 O   0.00000000000000  0.66666666666667  0.82931053362736  15.999 > 28   53 O   0.33333333333333  0.66666666666667  0.82931053362736  15.999 > 28   54 O   0.66666666666667  0.66666666666667  0.82931053362736  15.999 > 28   55 Na  0.00000000000000  0.16666666666667  0.27203019469644  22.990 > 55   56 Na  0.33333333333333  0.16666666666667  0.27203019469644  22.990 > 55   57 Na  0.66666666666667  0.16666666666667  0.27203019469644  22.990 > 55   58 Na  0.00000000000000  0.50000000000000  0.27203019469644  22.990 > 55   59 Na  0.33333333333333  0.50000000000000  0.27203019469644  22.990 > 55   60 Na  0.66666666666667  0.50000000000000  0.27203019469644  22.990 > 55   61 Na  0.00000000000000  0.83333333333333  0.27203019469644  22.990 > 55   62 Na  0.33333333333333  0.83333333333333  0.27203019469644  22.990 > 55   63 Na  0.66666666666667  0.83333333333333  0.27203019469644  22.990 > 55   64 Na  0.16666666666667  0.00000000000000  0.27203019469644  22.990 > 64   65 Na  0.50000000000000  0.00000000000000  0.27203019469644  22.990 > 64   66 Na  0.83333333333333  0.00000000000000  0.27203019469644  22.990 > 64   67 Na  0.16666666666667  0.33333333333333  0.27203019469644  22.990 > 64   68 Na  0.50000000000000  0.33333333333333  0.27203019469644  22.990 > 64   69 Na  0.83333333333333  0.33333333333333  0.27203019469644  22.990 > 64   70 Na  0.16666666666667  0.66666666666667  0.27203019469644  22.990 > 64   71 Na  0.50000000000000  0.66666666666667  0.27203019469644  22.990 > 64   72 Na  0.83333333333333  0.66666666666667  0.27203019469644  22.990 > 64   73 Na  0.16666666666667  0.00000000000000  0.77203019469644  22.990 > 55   74 Na  0.50000000000000  0.00000000000000  0.77203019469644  22.990 > 55   75 Na  0.83333333333333  0.00000000000000  0.77203019469644  22.990 > 55   76 Na  0.16666666666667  0.33333333333333  0.77203019469644  22.990 > 55   77 Na  0.50000000000000  0.33333333333333  0.77203019469644  22.990 > 55   78 Na  0.83333333333333  0.33333333333333  0.77203019469644  22.990 > 55   79 Na  0.16666666666667  0.66666666666667  0.77203019469644  22.990 > 55   80 Na  0.50000000000000  0.66666666666667  0.77203019469644  22.990 > 55   81 Na  0.83333333333333  0.66666666666667  0.77203019469644  22.990 > 55   82 Na  0.00000000000000  0.16666666666667  0.77203019469644  22.990 > 64   83 Na  0.33333333333333  0.16666666666667  0.77203019469644  22.990 > 64   84 Na  0.66666666666667  0.16666666666667  0.77203019469644  22.990 > 64   85 Na  0.00000000000000  0.50000000000000  0.77203019469644  22.990 > 64   86 Na  0.33333333333333  0.50000000000000  0.77203019469644  22.990 > 64   87 Na  0.66666666666667  0.50000000000000  0.77203019469644  22.990 > 64   88 Na  0.00000000000000  0.83333333333333  0.77203019469644  22.990 > 64   89 Na  0.33333333333333  0.83333333333333  0.77203019469644  22.990 > 64   90 Na  0.66666666666667  0.83333333333333  0.77203019469644  22.990 > 64   91 Cu  0.16666666666667  0.16666666666667  0.49903880343342  63.546 > 91   92 Cu  0.50000000000000  0.16666666666667  0.49903880343342  63.546 > 91   93 Cu  0.83333333333333  0.16666666666667  0.49903880343342  63.546 > 91   94 Cu  0.16666666666667  0.50000000000000  0.49903880343342  63.546 > 91   95 Cu  0.50000000000000  0.50000000000000  0.49903880343342  63.546 > 91   96 Cu  0.83333333333333  0.50000000000000  0.49903880343342  63.546 > 91   97 Cu  0.16666666666667  0.83333333333333  0.49903880343342  63.546 > 91   98 Cu  0.50000000000000  0.83333333333333  0.49903880343342  63.546 > 91   99 Cu  0.83333333333333  0.83333333333333  0.49903880343342  63.546 > 91  100 Cu  0.00000000000000  0.00000000000000  0.99903880343342  63.546 > 91  101 Cu  0.33333333333333  0.00000000000000  0.99903880343342  63.546 > 91  102 Cu  0.66666666666667  0.00000000000000  0.99903880343342  63.546 > 91  103 Cu  0.00000000000000  0.33333333333333  0.99903880343342  63.546 > 91  104 Cu  0.33333333333333  0.33333333333333  0.99903880343342  63.546 > 91  105 Cu  0.66666666666667  0.33333333333333  0.99903880343342  63.546 > 91  106 Cu  0.00000000000000  0.66666666666667  0.99903880343342  63.546 > 91  107 Cu  0.33333333333333  0.66666666666667  0.99903880343342  63.546 > 91  108 Cu  0.66666666666667  0.66666666666667  0.99903880343342  63.546 > 91----------------------------------------------------------------------------NAC parameters were read from "phono3py.yaml".--------------------------- Dielectric constant ----------------------------            3.3400940    0.0000000    0.0000000            0.0000000    3.3400940    0.0000000            0.0000000    0.0000000    3.4803409-------------------------- Born effective charges --------------------------    1 K     1.1098283    0.0000000    0.0000000            0.0000000    1.1098283    0.0000000            0.0000000    0.0000000    1.1499791    2 O    -1.2017887    0.0000000    0.0000000            0.0000000   -1.2017887    0.0000000            0.0000000    0.0000000   -2.3213243    3 O    -0.9895916    0.0000000    0.0000000            0.0000000   -0.9895916    0.0000000            0.0000000    0.0000000   -2.5544378    4 Na    1.0151837    0.0000000    0.0000000            0.0000000    0.9405108    0.0000000            0.0000000    0.0000000    1.0156978    5 Na    0.9405108    0.0000000    0.0000000            0.0000000    1.0151837    0.0000000            0.0000000    0.0000000    1.0156978    6 Cu   -0.8741425    0.0000000    0.0000000            0.0000000   -0.8741425    0.0000000            0.0000000    0.0000000    1.6943874----------------------------------------------------------------------------fc3 was read from "fc3.hdf5".fc2 was read from "fc2.hdf5".----------------------------- Force constants ------------------------------Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz)Max drift of fc2: 0.00000000 (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 ... Generalized regular grid: [ 1 11 110 ]Grid generation matrix:  [ 0 11 11 ]  [ 11 0 11 ]  [ 5 5 0 ]fc3-r2q-transformation over three atoms: True--------------------------- Phonon calculations ----------------------------Use NAC by Gonze et al. (no real space sum in current implementation)  PRB 50, 13035(R) (1994), PRB 55, 10355 (1997)  G-cutoff distance: 0.89, Number of G-points: 309, Lambda: 0.17Running harmonic phonon calculations...-------------------- Lattice thermal conductivity (RTA) --------------------======================= Grid point 0 (1/126) =======================q-point: ( 0.00  0.00  0.00)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|  -0.000   (   0.000    0.000    0.000)    0.000  -0.000   (   0.000    0.000    0.000)    0.000  -0.000   (   0.000    0.000    0.000)    0.000   2.888   (   0.000   -0.000    0.000)    0.000   2.888   (   0.000   -0.000    0.000)    0.000   2.915   (   0.000    0.000   -0.000)    0.000   4.561   (  -0.000   -0.000   -0.000)    0.000   4.561   (   0.000    0.000    0.000)    0.000   4.620   (   0.000   -0.000    0.000)    0.000   6.076   (   0.000   -0.000   -0.000)    0.000   6.076   (  -0.000    0.000   -0.000)    0.000   6.529   (   0.000    0.000    0.000)    0.000   6.529   (   0.000    0.000   -0.000)    0.000   6.553   (   0.000   -0.000    0.000)    0.000  11.193   (   0.000    0.000    0.000)    0.000  11.193   (   0.000    0.000   -0.000)    0.000  18.389   (   0.000   -0.000   -0.000)    0.000  18.862   (   0.000   -0.000   -0.000)    0.000======================= Grid point 1 (2/126) =======================q-point: ( 0.00  0.00  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.275   (   8.439    8.439    0.000)   11.934   0.920   (  20.993   20.993    0.000)   29.689   1.398   (  31.861   31.861    0.000)   45.058   2.755   (  -5.580   -5.580    0.000)    7.891   2.824   (  -2.550   -2.550    0.000)    3.606   3.414   (  12.834   12.834    0.000)   18.150   4.466   (  -2.471   -2.471    0.000)    3.494   4.790   (   5.345    5.345    0.000)    7.558   5.294   (  -4.156   -4.156    0.000)    5.877   6.174   (   3.847    3.847    0.000)    5.441   6.204   (   5.719    5.719    0.000)    8.088   6.485   (  -2.884   -2.884    0.000)    4.079   6.495   (  -1.228   -1.228    0.000)    1.737   7.262   (   0.213    0.213    0.000)    0.301  11.097   (  -4.547   -4.547    0.000)    6.431  11.872   (  -2.678   -2.678    0.000)    3.787  18.353   (  -1.643   -1.643    0.000)    2.324  18.890   (   1.356    1.356    0.000)    1.917======================= Grid point 2 (3/126) =======================q-point: ( 0.00  0.00  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.731   (  13.568   13.568    0.000)   19.187   1.668   (  13.266   13.266    0.000)   18.760   2.532   (  -3.192   -3.192    0.000)    4.514   2.566   (  21.151   21.151    0.000)   29.912   2.753   (   1.949    1.949    0.000)    2.756   3.896   (   5.338    5.338    0.000)    7.549   4.334   (  -4.400   -4.400    0.000)    6.222   5.009   (   4.953    4.953    0.000)    7.005   5.183   (   2.591    2.591    0.000)    3.664   6.173   (  -5.899   -5.899    0.000)    8.343   6.373   (   0.254    0.254    0.000)    0.360   6.489   (   3.978    3.978    0.000)    5.625   6.604   (   8.610    8.610    0.000)   12.176   7.335   (   4.522    4.522    0.000)    6.396  10.831   (  -8.005   -8.005    0.000)   11.320  11.703   (  -5.516   -5.516    0.000)    7.800  18.266   (  -2.380   -2.380    0.000)    3.366  18.970   (   2.406    2.406    0.000)    3.403======================= Grid point 3 (4/126) =======================q-point: ( 0.00  0.00  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.386   (  17.530   17.530    0.000)   24.791   1.854   (  -2.644   -2.644    0.000)    3.740   2.671   (   8.365    8.365    0.000)   11.829   2.674   (  -0.475   -0.475    0.000)    0.672   3.363   (   5.308    5.308    0.000)    7.507   3.749   (  -2.780   -2.780    0.000)    3.931   4.074   (  -8.111   -8.111    0.000)   11.471   5.193   (   3.861    3.861    0.000)    5.460   5.565   (   9.401    9.401    0.000)   13.295   5.854   (  -6.800   -6.800    0.000)    9.616   6.425   (   3.243    3.243    0.000)    4.586   6.615   (   2.164    2.164    0.000)    3.060   7.003   (   9.285    9.285    0.000)   13.131   7.730   (  14.618   14.618    0.000)   20.672  10.459   (  -9.495   -9.495    0.000)   13.429  11.406   (  -8.949   -8.949    0.000)   12.656  18.169   (  -2.191   -2.191    0.000)    3.099  19.082   (   2.833    2.833    0.000)    4.007======================= Grid point 4 (5/126) =======================q-point: ( 0.00  0.00  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.702   (  -3.304   -3.304    0.000)    4.672   2.107   (  14.559   14.559    0.000)   20.589   2.667   (   0.131    0.131    0.000)    0.186   2.926   ( -12.387  -12.387    0.000)   17.517   2.982   (   5.239    5.239    0.000)    7.409   3.706   (  -8.439   -8.439    0.000)   11.935   4.045   (   6.374    6.374    0.000)    9.014   5.328   (   2.623    2.623    0.000)    3.710   5.345   ( -10.313  -10.313    0.000)   14.585   5.984   (   3.662    3.662    0.000)    5.179   6.650   (   7.125    7.125    0.000)   10.076   6.755   (   4.674    4.674    0.000)    6.611   7.310   (   4.970    4.970    0.000)    7.029   8.489   (  20.842   20.842    0.000)   29.475  10.082   (  -8.134   -8.134    0.000)   11.503  10.950   ( -13.076  -13.076    0.000)   18.492  18.091   (  -1.511   -1.511    0.000)    2.137  19.193   (   2.341    2.341    0.000)    3.311======================= Grid point 5 (6/126) =======================q-point: ( 0.00  0.00  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.615   (  -0.924   -0.924    0.000)    1.307   2.138   (  -6.256   -6.256    0.000)    8.847   2.673   (  -0.061   -0.061    0.000)    0.087   3.029   (  10.568   10.568    0.000)   14.945   3.066   (  -0.487   -0.487    0.000)    0.688   3.482   (  -1.953   -1.953    0.000)    2.761   4.075   (  -3.837   -3.837    0.000)    5.427   5.028   (  -3.693   -3.693    0.000)    5.222   5.405   (   0.982    0.982    0.000)    1.388   6.047   (   0.347    0.347    0.000)    0.491   6.920   (   2.348    2.348    0.000)    3.321   6.936   (   5.094    5.094    0.000)    7.203   7.402   (  -0.006   -0.006    0.000)    0.009   9.349   (  18.494   18.494    0.000)   26.154   9.833   (  -3.295   -3.295    0.000)    4.659  10.339   ( -15.064  -15.064    0.000)   21.304  18.047   (  -0.548   -0.548    0.000)    0.775  19.263   (   0.917    0.917    0.000)    1.297======================= Grid point 11 (7/126) =======================q-point: (-0.19 -0.01  0.01)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.653   (  -0.000   -5.319  -22.501)   23.121   1.197   (  -0.000  -26.986   -2.590)   27.110   2.023   (  -0.000  -38.716   -6.692)   39.291   2.838   (  -0.000    2.293   -4.387)    4.950   3.075   (   0.000   -9.701    4.141)   10.548   3.342   (  -0.000    0.702   -9.619)    9.644   4.467   (   0.000    0.960    7.770)    7.829   4.722   (   0.000   -5.931    2.327)    6.372   5.194   (   0.000   -5.331   16.374)   17.220   5.821   (   0.000    9.152    1.609)    9.292   6.301   (   0.000    1.769    6.948)    7.169   6.570   (  -0.000    2.427   -6.358)    6.805   6.895   (  -0.000  -18.820  -12.940)   22.839   7.144   (   0.000    3.259    2.231)    3.950  11.189   (  -0.000    0.201   -0.046)    0.206  11.454   (   0.000    7.490   15.114)   16.868  18.264   (   0.000    3.995    0.322)    4.008  19.418   (  -0.000   14.139  -37.115)   39.717======================= Grid point 12 (8/126) =======================q-point: (-0.19 -0.01  0.10)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.601   (   1.203   -1.203  -27.167)   27.221   0.934   (  21.064  -21.064   -1.427)   29.823   1.486   (  30.448  -30.448  -10.019)   44.210   2.721   (  -4.982    4.982    3.596)    7.911   3.010   (   0.746   -0.746   -7.564)    7.637   3.496   (   1.493   -1.493   -9.835)   10.059   4.425   (  -2.030    2.030    4.260)    5.137   4.772   (   5.190   -5.190    1.779)    7.552   5.000   (   1.289   -1.289   22.716)   22.789   6.072   (   2.566   -2.566    9.444)   10.118   6.156   (   2.787   -2.787    3.834)    5.499   6.635   (  -0.621    0.621  -12.822)   12.852   6.647   (   1.229   -1.229  -11.588)   11.718   7.225   (  -0.455    0.455    5.267)    5.306  11.098   (  -4.515    4.515   -0.121)    6.386  11.584   (   5.479   -5.479   19.054)   20.569  18.321   (  -3.889    3.889   -2.576)    6.074  19.677   ( -21.844   21.844  -50.422)   59.132======================= Grid point 14 (9/126) =======================q-point: (-0.19 -0.01  0.28)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.893   (  12.985    5.384   -8.180)   16.264   1.820   (  20.054    6.748   -2.819)   21.346   2.597   (   5.391   -1.269   -5.083)    7.517   2.860   (   8.897   -0.563   -2.783)    9.339   3.087   (  16.994   -4.908   -7.104)   19.062   3.673   (   0.088   12.014    1.464)   12.103   4.409   (   0.793   -5.279    4.832)    7.200   4.943   (   8.358    6.216    5.548)   11.802   5.186   (   3.484   -5.158    9.364)   11.245   5.783   (  -9.111   13.906    1.202)   16.668   6.145   (  -8.359   -2.229    3.560)    9.355   6.492   (  -2.453   -0.984   -4.601)    5.306   7.044   (  -4.374    3.956   -0.022)    5.898   7.471   (  28.850   -0.386   -7.163)   29.728  10.914   (  -7.973  -12.866    2.008)   15.269  11.397   (  -9.611    6.905    6.707)   13.603  18.194   (  -1.763   -1.005    3.012)    3.632  19.242   (  -2.185   -1.362  -19.753)   19.920======================= Grid point 15 (10/126) =======================q-point: (-0.19 -0.01  0.37)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.383   (  15.126   12.005    2.192)   19.435   2.087   (   8.880   -9.045   -3.100)   13.048   2.703   (   4.102    3.171   -1.038)    5.287   2.891   (   4.461   -1.369  -10.269)   11.280   3.400   (  -3.274   -4.820    3.176)    6.636   3.754   (   9.151    5.089   -3.431)   11.019   4.208   (  -2.705  -12.840   -0.120)   13.122   5.109   (  -1.425    3.883    4.646)    6.220   5.419   (   6.023    9.178    8.452)   13.854   5.783   (  -7.970    0.359   -2.529)    8.370   6.109   (  -1.244   13.104    2.892)   13.476   6.512   (   0.549    6.087    1.322)    6.253   7.074   (  -3.355    8.026   -3.728)    9.465   8.081   (  29.457   -0.437   -3.734)   29.696  10.448   ( -13.199   -9.373    1.033)   16.221  11.263   (  -8.361   -2.044    3.617)    9.336  18.140   (  -0.951   -1.538    2.203)    2.851  19.232   (   2.534   -0.708   -8.835)    9.218======================= Grid point 16 (11/126) =======================q-point: (-0.19 -0.01  0.46)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.825   (   5.099    1.644    4.798)    7.192   2.140   (  13.573   -1.567   -0.567)   13.675   2.650   ( -11.070   -3.646    2.412)   11.902   2.894   (   0.474   -1.219   -8.970)    9.065   3.233   (   6.741   -4.570   -1.024)    8.208   3.835   (  -3.846  -13.130   -2.144)   13.848   4.075   (   6.106    2.478    0.432)    6.604   5.114   (  -8.864    3.944    1.760)    9.860   5.398   (  -1.085   -3.582    2.718)    4.625   5.853   (  -3.684   10.976   -4.626)   12.467   6.366   (  -2.356   12.533    1.890)   12.892   6.750   (   5.892    8.139    6.662)   12.056   7.268   (   1.360   11.177   -4.448)   12.106   8.626   (  16.280    5.273   -0.815)   17.133  10.055   (  -7.957   -4.411   -0.518)    9.113  10.966   (  -6.741  -11.282    1.736)   13.257  18.092   (  -0.236   -1.828    1.253)    2.229  19.291   (   4.330   -1.147   -2.702)    5.231======================= Grid point 17 (12/126) =======================q-point: (-0.19 -0.01  0.55)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.841   (   5.963   -5.937   -1.260)    8.508   2.053   (  -7.642   -5.619    2.213)    9.740   2.616   (   2.806    2.409    5.761)    6.846   3.055   (   5.773    6.467   -6.564)   10.874   3.214   (   2.564   -6.598    0.331)    7.086   3.649   (   4.872   -5.224   -3.942)    8.159   4.027   (  -3.271   -6.695    6.846)   10.119   4.990   (  -5.821    2.690   -6.195)    8.916   5.405   (   1.836   -0.034    2.957)    3.481   5.939   (  -6.766    8.321   -7.870)   13.303   6.512   ( -10.343   13.949    2.695)   17.573   7.015   (   5.592    3.530    5.389)    8.530   7.463   (   3.345    0.932   -0.504)    3.509   8.886   ( -13.742   17.559   -0.356)   22.300  10.062   (  11.753    0.820   -0.056)   11.781  10.541   (  -1.742  -19.332   -0.624)   19.420  18.061   (   0.720   -1.541    1.003)    1.975  19.343   (   4.049   -2.586   -0.055)    4.805======================= Grid point 22 (13/126) =======================q-point: (-0.38 -0.02  0.02)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.263   (  -0.000   -6.485  -12.260)   13.869   2.287   (  -0.000  -20.172   -6.207)   21.105   2.751   (  -0.000    1.125   -0.631)    1.290   3.072   (  -0.000   -3.233   -4.192)    5.294   3.328   (   0.000   -2.051    3.397)    3.968   3.830   (  -0.000  -11.365   -9.133)   14.580   4.258   (   0.000   -3.939   14.720)   15.238   5.004   (   0.000  -12.276    7.512)   14.392   5.269   (   0.000   -9.157    4.290)   10.112   5.410   (  -0.000   11.143   -5.581)   12.463   5.872   (   0.000   10.474    6.251)   12.198   6.542   (  -0.000    2.743   -3.081)    4.125   6.755   (   0.000    8.047    2.804)    8.521   8.302   (  -0.000  -30.191   -8.666)   31.410  10.612   (   0.000   20.648    8.565)   22.354  11.181   (  -0.000    0.172   -0.033)    0.176  18.156   (  -0.000    1.342   -0.744)    1.534  19.490   (  -0.000    5.033  -11.703)   12.740======================= Grid point 23 (14/126) =======================q-point: (-0.38 -0.02  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.172   (   1.497   -4.284  -19.669)   20.186   1.898   (   5.284  -20.712   -5.241)   22.008   2.624   (  -4.672   -5.828    1.171)    7.561   2.991   (   3.097   -9.799   -8.543)   13.364   3.205   (   1.501  -10.028   -3.689)   10.789   3.664   (   6.099   -4.052   -1.530)    7.481   4.268   (   2.217   -0.242    9.666)    9.920   4.777   (   2.124  -10.515   11.479)   15.711   5.025   (  -3.548   -5.775    6.388)    9.313   5.748   (   8.349    5.499    2.337)   10.267   6.030   (  -3.040    6.633    7.494)   10.459   6.615   (   1.953    6.738   -7.112)    9.990   7.001   (   7.648    3.418    4.548)    9.532   7.648   (  -3.239  -27.720  -10.363)   29.770  10.855   (  -8.963   10.791    4.058)   14.603  11.270   (   4.330    4.352    5.608)    8.315  18.192   (  -1.641    3.600   -2.880)    4.894  19.618   (  -3.694   10.078  -17.556)   20.577======================= Grid point 24 (15/126) =======================q-point: (-0.38 -0.02  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.142   (   3.062   -3.062  -21.004)   21.446   1.713   (  12.220  -12.220   -1.947)   17.390   2.523   (   0.834   -0.834    0.684)    1.363   2.807   (  19.887  -19.887  -10.848)   30.144   3.168   (  -0.554    0.554  -10.628)   10.657   3.737   (   3.924   -3.924    4.998)    7.468   4.271   (  -1.887    1.887    3.024)    4.033   4.706   (   7.785   -7.785   17.617)   20.774   4.907   (   3.895   -3.895    5.071)    7.488   5.909   (   0.370   -0.370   11.998)   12.010   6.043   (  -3.823    3.823    4.901)    7.298   6.889   (  -0.802    0.802   -8.567)    8.642   7.030   (   5.436   -5.436  -10.900)   13.338   7.237   (   8.721   -8.721    2.708)   12.627  10.844   (  -7.689    7.689   -0.575)   10.889  11.387   (  -1.461    1.461   10.302)   10.507  18.217   (  -3.879    3.879   -4.034)    6.810  19.697   (  -9.660    9.660  -19.675)   23.953======================= Grid point 27 (16/126) =======================q-point: (-0.38 -0.02  0.47)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.447   (   6.572    4.088   -0.599)    7.763   2.566   (   9.604   -6.423   -4.663)   12.459   2.710   (  -0.487   -1.018    1.094)    1.572   3.116   (   1.018   -0.430   -3.794)    3.951   3.293   (  -3.597    0.312    2.613)    4.456   4.002   (   6.171    1.676   -4.763)    7.973   4.310   (  -1.846   -2.864    3.651)    4.995   5.005   (  -1.318   -9.827    0.955)    9.961   5.298   (   6.964    4.389    2.905)    8.730   5.604   (  -7.356    3.667   -1.027)    8.283   5.777   (   1.191   15.140    3.346)   15.551   6.422   (  -1.666   -2.562    0.682)    3.132   6.790   (  -3.461   12.415   -1.043)   12.930   8.808   (  19.776   -5.089   -4.694)   20.953  10.165   ( -15.564   -1.907    4.332)   16.268  11.169   (  -0.800   -1.310    0.480)    1.608  18.133   (   0.013   -0.967    0.562)    1.119  19.429   (   0.545   -1.309   -4.375)    4.599======================= Grid point 28 (17/126) =======================q-point: ( 0.62 -0.02 -0.44)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.649   (   1.131    6.205    8.499)   10.584   2.443   ( -11.119  -10.646    1.909)   15.512   2.649   (   9.382   -2.904   -2.380)   10.106   3.048   (   2.200   -3.735   -1.922)    4.741   3.336   (  -1.656    2.267    0.638)    2.879   4.034   (  -6.296   -5.151   -3.148)    8.722   4.180   (   3.298   -2.467   -0.018)    4.118   4.894   (  -1.178    1.008   -0.861)    1.773   5.338   (  -5.621    1.993    0.118)    5.965   5.662   (   1.165    7.747   -3.862)    8.734   6.029   (  -9.122   10.961   -4.087)   14.834   6.529   (   7.001    8.685    7.092)   13.219   7.072   (   2.015   14.735   -1.759)   14.976   8.737   ( -13.230   -6.878    2.223)   15.076  10.084   (   9.655   -0.815   -1.372)    9.786  11.083   (   0.181   -6.152   -0.757)    6.201  18.113   (   0.835   -1.751    1.278)    2.323  19.446   (   4.403   -2.550    1.695)    5.363======================= Grid point 29 (18/126) =======================q-point: ( 0.62 -0.02 -0.35)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.748   (  -3.230    6.576    9.613)   12.086   1.990   (  -8.477   -7.608    2.219)   11.605   2.857   (  12.899   -3.149   -3.038)   13.621   3.068   (  -0.314    2.042    1.247)    2.413   3.224   (  -5.114   -2.657    3.390)    6.686   3.903   (   0.486   -0.441   -3.394)    3.457   4.073   (   1.127   -9.787    4.889)   10.998   4.956   (  -2.477    5.634   -9.090)   10.977   5.317   (  -4.100    5.237    0.865)    6.707   5.814   (  -3.093    8.076   -9.270)   12.678   5.970   ( -12.908    8.599   -4.374)   16.115   6.888   (   6.197    7.028    7.109)   11.761   7.396   (   1.925   10.339    1.728)   10.657   8.238   ( -24.809    4.003    2.190)   25.225  10.312   (  12.117   -3.580    0.956)   12.670  10.951   (   4.850   -9.997   -4.677)   12.056  18.099   (   1.718   -2.100    1.795)    3.253  19.482   (   5.512   -4.356    4.899)    8.565======================= Grid point 33 (19/126) =======================q-point: ( 0.43 -0.03  0.03)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.406   (   0.000    6.601    1.321)    6.732   2.637   (   0.000   12.343    3.808)   12.917   2.718   (   0.000   -1.375    0.445)    1.445   3.114   (   0.000    1.060    3.342)    3.506   3.287   (  -0.000   -4.315   -3.104)    5.316   3.993   (   0.000    4.666    4.805)    6.698   4.323   (  -0.000   -0.603   -5.657)    5.689   5.234   (  -0.000    8.499   -2.103)    8.755   5.235   (   0.000   -3.427    2.771)    4.407   5.425   (   0.000    3.472   -0.305)    3.486   5.659   (  -0.000   -8.628   -5.157)   10.052   6.508   (   0.000   -0.454    1.058)    1.151   6.602   (  -0.000   -5.268   -0.514)    5.293   8.864   (   0.000   20.791    5.355)   21.469  10.187   (  -0.000  -17.333   -5.488)   18.181  11.179   (   0.000   -0.058    0.012)    0.059  18.144   (  -0.000    0.037   -0.547)    0.548  19.443   (   0.000    0.478    4.597)    4.622======================= Grid point 34 (20/126) =======================q-point: (-0.57 -0.03  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.549   (   3.235    3.006   -9.582)   10.550   2.635   (  -6.925   -6.038   -1.707)    9.345   2.691   (  -1.757    1.022    0.313)    2.057   3.122   (  -1.510    0.540    2.848)    3.268   3.272   (   2.483    1.066   -0.192)    2.709   4.049   (   1.428   -2.636    0.325)    3.015   4.200   (   0.577   -6.447    7.318)    9.770   5.020   ( -11.832    1.273   -2.231)   12.108   5.270   (   3.321  -10.297   -0.256)   10.823   5.566   (   2.794    5.029    4.834)    7.514   5.748   (  15.054   -4.154   -0.437)   15.623   6.447   (  -0.989    4.083   -3.077)    5.207   6.783   (  12.698    2.514    1.790)   13.067   8.883   (  -5.926  -14.141   -2.869)   15.599  10.099   (  -0.507   10.826    2.380)   11.096  11.162   (  -1.915    0.084    0.234)    1.931  18.144   (  -1.088    1.088   -1.567)    2.196  19.490   (  -1.552    5.214   -1.860)    5.749======================= Grid point 35 (21/126) =======================q-point: (-0.57 -0.03  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.634   (   4.113   -2.022  -15.132)   15.810   2.193   ( -11.259  -11.803   -1.718)   16.402   2.785   (   6.372    1.293   -2.015)    6.807   3.108   (  -2.213    1.796    3.330)    4.383   3.282   (  -2.635   -0.641   -0.121)    2.715   3.898   (   0.683   -6.343    0.452)    6.395   4.214   (  -1.572    0.277   -0.206)    1.609   4.871   (   0.545   -1.943    5.120)    5.503   5.141   (   9.480  -13.407    2.983)   16.689   5.715   (   2.015    2.795    7.651)    8.391   5.832   (   5.557   -6.965    5.048)   10.241   6.699   (   9.519    8.822   -7.274)   14.878   7.125   (  13.837    2.331    1.647)   14.128   8.283   (  -5.762  -26.724   -5.234)   27.835  10.382   (  -3.179   15.627    1.954)   16.066  11.109   (  -5.099    2.262    3.021)    6.344  18.145   (  -2.241    2.399   -2.682)    4.239  19.572   (  -4.049    7.591   -6.078)   10.534======================= Grid point 36 (22/126) =======================q-point: (-0.57 -0.03  0.30)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.658   (   4.641   -4.641  -15.574)   16.901   1.905   (  -0.121    0.121   -0.623)    0.646   2.896   (  10.138  -10.138   -5.741)   15.444   3.086   (  -4.506    4.506    5.156)    8.197   3.239   (  -4.835    4.835   -0.829)    6.888   3.865   (   3.827   -3.827   -1.732)    5.682   4.141   (  -4.775    4.775   -3.170)    7.459   4.898   (   3.085   -3.085    9.803)   10.730   5.106   (  16.251  -16.251    3.092)   23.189   5.743   (   3.408   -3.408   11.394)   12.372   5.815   (   1.115   -1.115    4.737)    4.992   6.959   (   0.308   -0.308   -4.276)    4.298   7.404   (   4.012   -4.012   -5.546)    7.934   7.775   (  15.033  -15.033   -2.036)   21.357  10.522   (  -8.138    8.138   -1.472)   11.602  11.072   (  -6.017    6.017    6.512)   10.716  18.148   (  -2.957    2.957   -3.205)    5.270  19.612   (  -6.811    6.811   -7.466)   12.187======================= Grid point 40 (23/126) =======================q-point: ( 0.43 -0.03 -0.34)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.297   (  -7.279    3.494   10.734)   13.432   2.287   ( -16.520   -2.029    5.492)   17.526   2.703   (  -1.701   -2.294   -0.111)    2.858   3.067   (  -2.667    0.179    5.949)    6.522   3.366   (   0.149    1.542   -3.451)    3.783   3.811   ( -10.593   -0.600    7.116)   12.775   4.296   (  -1.867    1.098   -9.169)    9.422   4.996   (  -6.521   -2.830   -5.093)    8.744   5.130   (  -6.222   -0.046   -5.655)    8.408   5.669   (   1.210   13.990    0.710)   14.060   5.862   (   8.603   -0.199   -6.730)   10.924   6.494   (   4.823   -1.123    3.308)    5.955   6.899   (   5.951   10.491   -1.128)   12.114   8.259   ( -29.821   -3.770    7.995)   31.103  10.548   (  18.175   -5.406   -5.968)   19.879  11.199   (   2.299    0.903   -1.933)    3.136  18.144   (   1.204   -1.032    0.547)    1.678  19.469   (   4.469   -1.856   11.160)   12.164======================= Grid point 41 (24/126) =======================q-point: ( 0.43 -0.03 -0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.249   (  -7.084    6.200   14.631)   17.398   1.939   (  -9.908   -1.858    1.647)   10.214   2.618   (  -7.323    5.072    2.734)    9.318   3.046   (  -1.317    0.807    8.090)    8.236   3.282   ( -10.832    5.222    5.149)   13.081   3.735   (   1.294    0.890   -0.906)    1.813   4.236   (   1.260   -5.479   -1.396)    5.793   4.931   (  -9.362   11.360  -11.961)   18.967   5.018   (  -8.612    4.035   -8.073)   12.474   5.892   (   1.913    4.225   -7.480)    8.801   5.955   (   3.508   -3.107   -8.466)    9.676   6.739   (   9.028    8.712    9.245)   15.584   7.149   (  -1.641    6.057    2.001)    6.586   7.598   ( -25.185    0.212    5.778)   25.841  10.670   (   9.694   -8.604   -0.114)   12.962  11.304   (   5.837   -1.450   -8.113)   10.099  18.150   (   2.420   -2.142    1.625)    3.617  19.524   (   6.812   -4.967   15.200)   17.381======================= Grid point 44 (25/126) =======================q-point: ( 0.24 -0.04  0.04)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.840   (   0.000   12.230   10.261)   15.964   1.730   (   0.000   24.272    3.082)   24.467   2.595   (   0.000   11.278    9.036)   14.452   2.877   (   0.000    9.565    0.457)    9.576   3.273   (  -0.000    8.375   -2.944)    8.877   3.418   (   0.000    7.955    6.989)   10.589   4.462   (  -0.000    0.851   -9.482)    9.520   4.881   (  -0.000    9.166   -2.411)    9.478   5.289   (  -0.000    4.551  -11.263)   12.147   5.586   (   0.000  -13.131    1.988)   13.281   6.172   (  -0.000   -9.162   -3.157)    9.690   6.522   (   0.000   -2.072    4.026)    4.528   6.976   (  -0.000   -7.301   -2.213)    7.629   7.486   (   0.000   30.466    8.843)   31.724  11.185   (   0.000   -0.215    0.035)    0.218  11.188   (  -0.000  -17.275   -9.628)   19.776  18.205   (  -0.000   -1.927   -2.891)    3.475  19.258   (   0.000   -3.081   21.467)   21.687======================= Grid point 45 (26/126) =======================q-point: ( 0.24 -0.04  0.13)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.194   (   6.051   15.338    0.618)   16.500   2.199   (  -1.434   16.105    3.571)   16.559   2.668   (   0.446    3.167    3.890)    5.036   2.958   (  -5.164    2.992    3.793)    7.072   3.412   (   1.836    5.081   -1.066)    5.506   3.693   (   2.130    6.939    4.336)    8.456   4.426   (  -7.157   -0.106   -4.233)    8.316   5.079   (  -0.599    1.552   -4.176)    4.496   5.249   (   3.801    3.742   -4.906)    7.247   5.652   (  13.735   -2.024    0.319)   13.887   5.947   (   0.773   -9.777   -1.600)    9.937   6.461   (  -0.758   -0.472   -0.049)    0.895   6.900   (   8.224   -8.187    0.672)   11.624   8.123   (  -2.511   31.893    5.624)   32.482  10.647   (  -9.149  -17.649   -3.966)   20.272  11.241   (   3.521   -4.883   -2.088)    6.372  18.167   (  -0.932   -0.868   -2.461)    2.771  19.248   (  -0.640    2.476   10.739)   11.040======================= Grid point 46 (27/126) =======================q-point: ( 0.24 -0.04  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.670   (   9.518   11.875   -7.480)   16.958   2.292   (  -9.778   10.804    5.068)   15.428   2.710   (  -2.245   -3.977   -2.200)    5.069   2.954   (  -3.115    0.338    8.515)    9.073   3.313   (  -3.158   -0.198   -1.572)    3.533   3.983   (   4.565    8.928    3.404)   10.590   4.132   ( -13.000   -4.426    1.103)   13.777   5.021   (   3.119   -6.183   -3.790)    7.894   5.455   (   0.657   -1.695   -4.535)    4.886   5.641   (   6.918   -5.324    4.017)    9.609   6.099   (  13.606   -1.774    1.196)   13.773   6.569   (   7.888    4.469   -8.115)   12.167   7.000   (  13.250   -2.347    4.910)   14.324   8.631   (  -3.885   20.432    2.293)   20.924  10.151   (  -4.777  -13.065   -1.388)   13.980  11.131   (  -4.734   -4.340   -1.231)    6.539  18.130   (  -1.716    0.001   -1.523)    2.294  19.317   (  -1.280    5.565    3.607)    6.754======================= Grid point 47 (28/126) =======================q-point: ( 0.24 -0.04  0.31)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.915   (  -0.476    3.386   -7.494)    8.237   2.226   ( -10.142   -8.305    0.926)   13.141   2.581   (   0.277    8.473   -3.619)    9.218   2.954   (   2.223    4.616    7.650)    9.206   3.328   (  -3.404    0.934   -0.012)    3.529   3.790   (  -8.452   -0.034    3.602)    9.188   4.155   (  -4.292    1.841   -3.784)    6.011   4.919   (   3.748   -8.565    3.986)   10.163   5.395   (   1.028   -0.041   -2.304)    2.523   5.757   (   8.862   -4.902    5.430)   11.491   6.216   (  13.629  -12.153    0.122)   18.261   6.887   (   8.513    6.739   -7.076)   12.960   7.319   (  12.765    3.559    2.757)   13.536   8.685   (   1.381  -11.467   -0.617)   11.567  10.059   (  -0.255    9.171    0.449)    9.185  10.876   ( -12.801   -1.923    0.389)   12.951  18.100   (  -2.137    1.063   -1.298)    2.717  19.402   (  -2.948    6.276   -0.107)    6.934======================= Grid point 48 (29/126) =======================q-point: (-0.76 -0.04  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.945   (   1.251   -1.251   -1.996)    2.667   1.953   (  -3.426    3.426   -3.629)    6.054   2.692   (  -4.538    4.538   -2.193)    6.782   3.130   (  -0.617    0.617    5.867)    5.932   3.234   (   0.177   -0.177   -2.472)    2.484   3.765   (  -4.610    4.610    3.996)    7.647   4.022   (  -4.523    4.523   -7.311)    9.714   4.880   (   5.160   -5.160    8.618)   11.292   5.430   (  -0.129    0.129   -3.335)    3.340   5.788   (   7.268   -7.268    7.682)   12.831   6.226   (  16.209  -16.209   -0.462)   22.928   7.075   (   0.267   -0.267   -4.680)    4.695   7.525   (  -1.151    1.151   -1.115)    1.974   8.501   (  18.826  -18.826    0.187)   26.625  10.268   (  -5.327    5.327   -1.812)    7.748  10.635   ( -10.675   10.675    3.219)   15.436  18.088   (  -1.920    1.920   -1.335)    3.026  19.438   (  -5.043    5.043   -1.230)    7.238======================= Grid point 53 (30/126) =======================q-point: ( 0.24 -0.04 -0.15)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.680   (  -7.058    3.475   18.709)   20.296   1.382   ( -21.748   14.457    2.126)   26.201   2.195   ( -32.942   16.166    6.806)   37.321   2.686   (  -1.857  -10.497   -1.340)   10.744   2.974   (   1.323   -4.506   12.747)   13.584   3.605   (  -6.131   16.016   -1.765)   17.240   4.405   (   0.310   -4.443   -4.181)    6.109   4.823   (  -2.735    7.903   -5.188)    9.841   5.092   (  -5.401   -7.105  -14.340)   16.891   5.968   (   7.913   11.431   -5.224)   14.852   6.252   (   0.454   -4.267   -7.497)    8.638   6.561   (   4.006    0.868    9.978)   10.787   6.861   ( -18.241   -3.979   13.687)   23.150   7.201   (   1.117    4.842   -3.729)    6.213  11.022   (   3.576  -10.257   -0.382)   10.869  11.574   (   5.777    5.375  -12.864)   15.092  18.244   (   3.143   -1.731   -1.298)    3.815  19.355   (   9.433   -4.937   32.705)   34.394======================= Grid point 55 (31/126) =======================q-point: ( 0.05 -0.05  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.181   (   0.000   10.183    0.000)   10.183   0.598   (   0.000   28.406    0.000)   28.406   1.037   (   0.000   48.159    0.000)   48.159   2.842   (  -0.000   -5.061   -0.000)    5.061   2.951   (   0.000    5.853    0.000)    5.853   3.153   (   0.000    4.628    0.000)    4.628   4.554   (  -0.000   -0.542   -0.000)    0.542   4.649   (   0.000    2.832    0.000)    2.832   5.387   (   0.000    0.057    0.000)    0.057   6.017   (  -0.000   -5.725   -0.000)    5.725   6.247   (   0.000   14.752    0.000)   14.752   6.520   (  -0.000   -2.593   -0.000)    2.593   6.522   (  -0.000   -0.674   -0.000)    0.674   7.236   (  -0.000   -2.520   -0.000)    2.520  11.192   (  -0.000   -0.118   -0.000)    0.118  11.855   (  -0.000   -6.984   -0.000)    6.984  18.371   (  -0.000   -1.720   -0.000)    1.720  18.876   (   0.000    1.436    0.000)    1.436======================= Grid point 56 (32/126) =======================q-point: ( 0.05 -0.05  0.14)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.510   (   7.328   13.870    0.000)   15.687   1.362   (  14.856   21.392    0.000)   26.044   2.124   (  16.163   32.954    0.000)   36.704   2.696   (  -7.796    0.842    0.000)    7.841   2.818   (  -9.108    2.076    0.000)    9.341   3.627   (  18.755    7.210   -0.000)   20.093   4.445   (  -6.167   -0.115    0.000)    6.168   4.867   (   8.331    3.850   -0.000)    9.177   5.268   (  -8.469    0.802    0.000)    8.507   6.022   (  15.073  -10.052   -0.000)   18.117   6.335   (  -2.830   -5.875   -0.000)    6.521   6.465   (  -2.233   -0.588    0.000)    2.309   6.686   (  -5.179   24.449    0.000)   24.991   7.233   (   5.424   -2.284   -0.000)    5.885  11.025   ( -10.709   -3.151    0.000)   11.163  11.726   (   3.452  -10.569   -0.000)   11.118  18.307   (  -1.452   -2.668   -0.000)    3.038  18.934   (   1.118    2.877    0.000)    3.086======================= Grid point 57 (33/126) =======================q-point: ( 0.05 -0.05  0.23)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.056   (  13.705   17.215    0.000)   22.005   1.892   (   0.904    8.526    0.000)    8.574   2.552   (   2.485    5.606    0.000)    6.132   2.714   (  -4.063    1.609    0.000)    4.370   3.126   (  12.242   22.802    0.000)   25.880   3.838   (  -0.433   -8.832   -0.000)    8.842   4.263   (  -9.672   -2.639    0.000)   10.026   5.114   (   3.480    4.558    0.000)    5.735   5.332   (   9.627   10.054    0.000)   13.919   5.984   (  -2.901  -11.177   -0.000)   11.547   6.252   (  13.273   -8.109   -0.000)   15.554   6.485   (   5.113    1.324   -0.000)    5.281   6.976   (  -4.606   12.337    0.000)   13.169   7.510   (   7.832   17.058    0.000)   18.770  10.670   ( -11.196   -8.240    0.000)   13.902  11.525   (  -1.916  -10.757   -0.000)   10.927  18.214   (  -2.056   -2.479   -0.000)    3.220  19.036   (   1.632    3.975    0.000)    4.297======================= Grid point 58 (34/126) =======================q-point: ( 0.05 -0.05  0.32)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.739   (  13.719   13.003   -0.000)   18.902   1.862   (  -5.831    6.825    0.000)    8.977   2.683   (  -1.691    1.197    0.000)    2.072   2.835   (   6.157    4.545   -0.000)    7.653   3.243   ( -13.588   -9.494    0.000)   16.576   3.878   (   8.907   10.709    0.000)   13.928   3.927   ( -12.507   -6.198    0.000)   13.958   5.245   (   4.647    0.185   -0.000)    4.651   5.550   (  -7.971   -8.472   -0.000)   11.632   5.830   (  10.065    1.825   -0.000)   10.229   6.421   (  13.134   -1.787   -0.000)   13.255   6.687   (   4.812    6.312    0.000)    7.937   7.166   (   7.761    4.092   -0.000)    8.774   8.150   (  10.503   24.387    0.000)   26.553  10.257   (  -8.892  -10.296   -0.000)   13.605  11.187   (  -8.964  -11.118   -0.000)   14.282  18.128   (  -1.993   -1.608    0.000)    2.561  19.156   (   1.277    4.322    0.000)    4.507======================= Grid point 59 (35/126) =======================q-point: ( 0.05 -0.05  0.41)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.706   (  -7.445    3.528    0.000)    8.239   2.245   (  -1.451   -3.906   -0.000)    4.167   2.683   (  -0.913    0.603    0.000)    1.094   2.789   (   9.110    4.193   -0.000)   10.028   3.121   (  -4.472    6.191    0.000)    7.637   3.587   (  -8.053   -2.605    0.000)    8.464   4.124   (  -0.179    1.026    0.000)    1.042   5.114   (  -3.776  -10.521   -0.000)   11.178   5.389   (   0.216    2.437    0.000)    2.446   5.977   (   6.263   -3.719   -0.000)    7.284   6.679   (  13.639   -2.687   -0.000)   13.901   6.923   (   2.370    8.732    0.000)    9.048   7.381   (   3.545    2.853   -0.000)    4.551   8.893   (  20.101   13.548   -0.000)   24.241   9.943   (  -5.220   -3.115    0.000)    6.079  10.690   ( -16.688  -10.810    0.000)   19.883  18.069   (  -1.497   -0.513    0.000)    1.583  19.254   (  -0.034    3.595    0.000)    3.595======================= Grid point 60 (36/126) =======================q-point: ( 0.05 -0.05 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.667   (  -5.682    5.682    0.000)    8.036   2.048   (   1.730   -1.730   -0.000)    2.447   2.682   (  -0.879    0.879    0.000)    1.244   3.102   (  -4.009    4.009    0.000)    5.670   3.115   (   3.238   -3.238   -0.000)    4.579   3.518   (  -4.465    4.465    0.000)    6.314   4.022   (  -0.357    0.357    0.000)    0.505   4.968   (   2.272   -2.272   -0.000)    3.213   5.424   (  -0.809    0.809    0.000)    1.144   5.996   (   4.927   -4.927   -0.000)    6.968   6.798   (   9.773   -9.773   -0.000)   13.821   7.072   (  -3.372    3.372    0.000)    4.768   7.427   (  -2.507    2.507    0.000)    3.545   9.295   (  18.710  -18.710   -0.000)   26.460  10.062   ( -15.159   15.159    0.000)   21.438  10.146   (  -1.639    1.639    0.000)    2.319  18.048   (  -0.605    0.605    0.000)    0.855  19.293   (  -1.929    1.929    0.000)    2.729======================= Grid point 67 (37/126) =======================q-point: (-0.15 -0.05  0.05)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.602   (  -0.000   -2.860  -27.061)   27.212   0.617   (  -0.000  -28.586   -1.956)   28.653   1.143   (  -0.000  -41.591  -14.096)   43.915   2.864   (   0.000   -0.747    2.940)    3.033   2.905   (   0.000   -6.053    4.830)    7.744   3.532   (  -0.000   21.404  -23.706)   31.939   4.489   (   0.000    0.900    6.653)    6.713   4.636   (   0.000   -2.528    1.712)    3.053   5.019   (   0.000  -13.910   23.685)   27.467   5.962   (   0.000    4.449    4.389)    6.249   6.184   (   0.000  -11.359    7.665)   13.703   6.615   (  -0.000    1.683   -8.343)    8.511   6.669   (  -0.000   -3.654  -12.650)   13.167   7.225   (   0.000    3.795    4.809)    6.126  11.193   (  -0.000    0.120   -0.053)    0.131  11.443   (   0.000  -10.052   19.774)   22.182  18.370   (  -0.000    5.573   -4.228)    6.995  20.029   (  -0.000   52.120  -56.587)   76.932======================= Grid point 70 (38/126) =======================q-point: (-0.15 -0.05  0.33)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.071   (  14.263   11.640   -3.730)   18.784   1.908   (   9.094    0.165   -1.418)    9.205   2.575   (   6.070    2.065   -1.816)    6.664   2.843   (   0.422   -1.610  -11.508)   11.628   3.171   (  21.043   10.723   -4.652)   24.071   3.796   (  -7.387   -0.547    3.843)    8.344   4.256   (  -2.264   -8.773    0.734)    9.090   5.076   (   4.150    4.218    4.033)    7.161   5.218   (  11.326   10.061   10.624)   18.503   5.964   (  -8.938   -1.101    2.040)    9.234   6.161   (  -2.484    6.195    9.839)   11.889   6.557   (  -4.775    8.024   -8.212)   12.435   7.060   (   7.651   -1.657   -6.554)   10.210   7.526   (  20.018    6.105   -1.888)   21.013  10.670   (  -8.607  -10.521   -0.023)   13.593  11.454   (  -9.038   -1.456    6.777)   11.390  18.167   (  -1.539   -1.519    3.021)    3.715  19.213   (  -0.657   -1.234  -15.652)   15.714======================= Grid point 71 (39/126) =======================q-point: (-0.15 -0.05  0.42)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.670   (  12.008   14.356    4.749)   19.309   1.889   (   6.702   -7.944   -2.355)   10.657   2.778   (  -1.255    2.936   -6.430)    7.180   2.896   (   5.000    2.359   -7.403)    9.240   3.200   (  -8.031  -13.293    4.073)   16.056   3.900   (   8.613    6.598   -1.912)   11.017   3.941   (  -4.604  -10.915   -1.652)   11.961   5.209   (   0.054    4.959    3.692)    6.182   5.476   (  -5.211   -4.534    6.337)    9.374   5.875   (   1.764    9.033   -3.294)    9.775   6.350   (   2.789    9.842    8.344)   13.201   6.673   (   0.804    7.109   -1.011)    7.225   7.235   (   3.104    7.423   -5.705)    9.863   8.163   (  24.550    9.654   -1.287)   26.411  10.267   ( -10.459   -8.139   -0.955)   13.287  11.148   ( -10.137   -8.888    3.929)   14.043  18.106   (  -1.094   -1.649    1.770)    2.655  19.223   (   2.158   -0.053   -6.400)    6.754======================= Grid point 72 (40/126) =======================q-point: (-0.15 -0.05  0.51)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.738   (   3.177   -5.945   -3.197)    7.460   2.175   (  -0.132   -0.442    6.044)    6.062   2.606   (  -4.445    1.019    4.098)    6.131   2.931   (   4.872    4.734   -9.376)   11.578   3.132   (   5.380   -4.588   -1.241)    7.178   3.614   (  -1.937   -7.718   -2.912)    8.473   4.081   (  -0.825   -2.026    4.440)    4.950   5.135   (  -8.169   -1.779   -2.289)    8.668   5.360   (   2.484    0.077    2.907)    3.825   6.049   (  -3.181    7.201   -7.025)   10.551   6.626   (  -0.253   12.229    6.173)   13.701   6.877   (   6.574    3.913    3.195)    8.290   7.408   (   2.037    2.193   -2.219)    3.726   8.895   (  13.524   19.945   -0.122)   24.098   9.965   (  -3.550   -4.457   -2.146)    6.089  10.668   (  -9.858  -16.886    2.145)   19.670  18.058   (  -0.408   -1.375    1.002)    1.749  19.273   (   2.635   -0.464   -1.800)    3.224======================= Grid point 73 (41/126) =======================q-point: (-0.15 -0.05 -0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.713   (   5.187   -5.187   -4.590)    8.652   2.012   (  -1.439    1.439    3.715)    4.236   2.656   (  -1.063    1.063    2.147)    2.621   3.099   (   3.734   -3.734    0.304)    5.289   3.155   (  -1.258    1.258   -3.086)    3.562   3.551   (   4.248   -4.248   -3.725)    7.069   3.941   (   0.262   -0.262    7.579)    7.588   5.035   (  -1.941    1.941   -6.241)    6.818   5.396   (   0.771   -0.771    2.974)    3.167   6.083   (  -5.111    5.111   -8.883)   11.452   6.761   (  -9.255    9.255    4.517)   13.846   7.020   (   3.245   -3.245    4.508)    6.433   7.426   (   2.420   -2.420    0.085)    3.423   9.305   ( -18.688   18.688   -0.943)   26.445  10.090   (  15.101  -15.101   -2.867)   21.548  10.120   (   1.743   -1.743    2.593)    3.578  18.038   (   0.566   -0.566    0.844)    1.163  19.298   (   1.701   -1.701   -0.433)    2.445======================= Grid point 79 (42/126) =======================q-point: (-0.34 -0.06  0.06)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.163   (  -0.000   -2.765  -21.682)   21.857   1.820   (  -0.000  -24.642   -6.378)   25.454   2.694   (  -0.000  -10.585   -0.602)   10.602   2.970   (  -0.000   -6.023   -8.664)   10.552   3.189   (   0.000  -10.304    5.874)   11.861   3.571   (  -0.000   -9.930   -8.562)   13.112   4.228   (   0.000    1.040   14.184)   14.222   4.771   (   0.000   -9.400    9.932)   13.675   5.089   (   0.000   -8.174    7.718)   11.242   5.634   (  -0.000    9.662   -2.761)   10.049   6.064   (   0.000    7.391    7.588)   10.592   6.617   (  -0.000    4.339   -5.532)    7.031   6.899   (   0.000    5.364    6.336)    8.302   7.687   (  -0.000  -27.763  -11.173)   29.927  10.986   (   0.000   14.646   10.261)   17.883  11.186   (  -0.000    0.230   -0.053)    0.236  18.210   (  -0.000    4.214   -3.414)    5.424  19.660   (  -0.000   11.724  -18.364)   21.788======================= Grid point 80 (43/126) =======================q-point: (-0.34 -0.06  0.15)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.117   (  -0.336   -1.804  -23.719)   23.790   1.442   (  12.841  -22.754   -4.108)   26.448   2.373   (  14.533  -24.212  -11.204)   30.380   2.654   (  -8.457   -9.269    0.972)   12.585   3.138   (   5.733   -2.637   -1.758)    6.550   3.601   (   6.197   -0.321   -2.351)    6.635   4.288   (   0.535    1.541    7.974)    8.139   4.627   (   1.196   -3.697   14.588)   15.097   4.880   (  -1.591   -7.958    6.836)   10.611   5.858   (   4.173    4.519    5.623)    8.334   6.111   (  -2.845    0.637    6.895)    7.486   6.759   (   5.268    6.619   -9.528)   12.742   7.026   (   3.059   -1.944    5.722)    6.774   7.165   (  -1.954  -16.686   -8.706)   18.922  11.009   (  -8.529    5.058    0.933)    9.960  11.349   (   5.106    2.098    9.396)   10.898  18.298   (  -3.371    6.421   -6.198)    9.540  19.901   (  -8.837   17.446  -21.867)   29.336======================= Grid point 84 (44/126) =======================q-point: (-0.34 -0.06  0.52)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.575   (   5.923    8.017    1.070)   10.025   2.383   (  11.683  -10.406   -4.178)   16.194   2.683   (  -7.056   -1.965    1.187)    7.420   3.057   (  -2.909   -4.375   -1.618)    5.497   3.309   (   3.138    0.658   -1.500)    3.540   4.047   (   7.588    0.847   -2.072)    7.911   4.159   (  -3.450  -10.596   -1.604)   11.259   4.941   (  -3.175    3.506    3.901)    6.132   5.373   (   2.190    1.147    3.618)    4.382   5.694   (  -4.145    7.768   -0.836)    8.844   6.101   (   3.505   12.895    2.117)   13.530   6.460   (  -1.664    6.975    1.281)    7.284   7.070   (  -1.925   14.082   -2.045)   14.359   8.685   (  17.433   -5.574   -3.090)   18.561  10.117   ( -11.862   -2.530    1.938)   12.283  11.104   (  -2.291   -5.472    1.391)    6.093  18.106   (  -0.066   -1.579    0.583)    1.685  19.393   (   0.658   -2.080   -3.762)    4.349======================= Grid point 85 (45/126) =======================q-point: ( 0.66 -0.06 -0.39)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.784   (  -0.158    5.904    4.984)    7.728   2.220   ( -11.664  -10.474    1.628)   15.761   2.563   (  13.152   -4.154   -1.339)   13.857   3.032   (   2.072    2.349   -0.852)    3.246   3.327   (  -2.201   -3.437    0.760)    4.151   3.879   (   0.782   -7.591   -4.519)    8.869   4.070   (  -0.039   -7.823    4.146)    8.854   5.008   (  -2.824    7.182   -4.640)    9.005   5.346   (  -0.070   -0.619    1.962)    2.059   5.881   (  -3.091   12.450   -6.166)   14.233   6.228   (  -9.680    9.208   -1.634)   13.460   6.745   (   6.558   10.221    6.740)   13.889   7.356   (   1.650   10.892   -0.684)   11.037   8.668   ( -13.013    2.385    1.079)   13.274  10.067   (   8.794   -0.855   -0.219)    8.838  10.892   (   0.666  -12.443   -1.317)   12.530  18.075   (   0.614   -1.803    0.942)    2.125  19.389   (   3.306   -2.735    1.347)    4.498======================= Grid point 86 (46/126) =======================q-point: ( 0.66 -0.06 -0.30)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.889   (  -0.015    0.015    3.760)    3.761   1.890   (  -4.355    4.355    2.032)    6.485   2.694   (   9.245   -9.245   -2.377)   13.289   3.168   (  -1.901    1.901    4.324)    5.091   3.171   (  -0.179    0.179    2.187)    2.201   3.880   (   3.467   -3.467   -7.267)    8.766   3.880   (   5.610   -5.610    6.794)   10.445   5.056   (  -3.819    3.819   -8.964)   10.465   5.355   (   0.445   -0.445    3.896)    3.947   5.995   (  -8.056    8.056  -13.490)   17.657   6.204   ( -12.802   12.802    2.113)   18.228   6.961   (   0.033   -0.033    6.711)    6.711   7.492   (   0.074   -0.074    2.392)    2.394   8.503   ( -19.541   19.541   -0.016)   27.636  10.220   (   5.975   -5.975    3.082)    8.995  10.719   (  11.072  -11.072   -5.333)   16.542  18.062   (   1.398   -1.398    1.080)    2.253  19.398   (   3.604   -3.604    3.104)    5.968======================= Grid point 90 (47/126) =======================q-point: (-0.53 -0.07  0.07)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.516   (  -0.000    3.777   -9.675)   10.386   2.709   (   0.000    0.549    0.974)    1.118   2.713   (  -0.000   -6.409   -3.891)    7.498   3.122   (   0.000   -0.137    3.078)    3.081   3.250   (   0.000    1.974    0.391)    2.013   4.042   (  -0.000    0.199   -0.112)    0.228   4.176   (   0.000  -11.858    9.317)   15.081   5.201   (  -0.000  -12.718    0.107)   12.718   5.304   (  -0.000    8.565   -4.106)    9.498   5.413   (   0.000   -5.196    3.078)    6.039   5.592   (   0.000    3.876    0.285)    3.887   6.525   (  -0.000    2.168   -0.711)    2.281   6.584   (   0.000    3.367    1.318)    3.615   8.949   (  -0.000  -14.706   -3.113)   15.032  10.104   (   0.000   11.434    2.840)   11.782  11.179   (  -0.000    0.078   -0.010)    0.079  18.156   (  -0.000    1.203   -1.682)    2.068  19.506   (  -0.000    5.518   -1.903)    5.837======================= Grid point 91 (48/126) =======================q-point: (-0.53 -0.07  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.580   (   1.334    0.363  -17.826)   17.879   2.397   (  -6.009  -15.453   -5.546)   17.484   2.691   (   0.714   -2.206    0.373)    2.348   3.121   (   0.282   -0.289    1.342)    1.402   3.299   (   1.494   -0.508    3.034)    3.420   3.896   (  -1.254   -8.408    0.990)    8.558   4.126   (   7.816   -2.497    6.072)   10.208   4.890   (   1.853  -15.610    4.139)   16.256   5.075   (  -6.397   -6.114    1.082)    8.915   5.644   (   7.296    1.351    3.865)    8.366   5.748   (   1.066    5.947    4.007)    7.250   6.568   (   1.893    7.046   -4.099)    8.369   6.855   (  11.773    3.448    3.249)   12.691   8.404   (  -4.830  -27.657   -6.696)   28.864  10.441   (  -2.295   18.153    4.851)   18.930  11.171   (  -1.334    0.815    1.026)    1.870  18.186   (  -1.491    3.189   -3.435)    4.919  19.642   (  -2.495    9.326   -6.675)   11.736======================= Grid point 92 (49/126) =======================q-point: (-0.53 -0.07  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.587   (   1.928   -2.140  -19.828)   20.036   1.964   (  -4.428  -10.413   -0.894)   11.351   2.676   (   9.249   -9.465   -3.280)   13.634   3.128   (  -0.214    0.363    0.273)    0.502   3.299   (   1.083   -0.987    2.358)    2.776   3.785   (   3.039   -4.174   -4.371)    6.765   4.209   (   0.404   -0.292    1.159)    1.262   4.719   (  14.623  -15.907    8.019)   23.047   4.881   (   0.281   -6.404    6.651)    9.237   5.708   (   0.534   -0.915    8.972)    9.034   5.840   (  -2.330    4.568    6.075)    7.950   6.866   (   7.742    5.983   -3.405)   10.360   7.160   (  14.357    1.618   -0.458)   14.455   7.723   (  -4.790  -25.716   -6.591)   26.975  10.666   (  -6.007   11.255    0.311)   12.762  11.170   (  -2.999    3.348    5.557)    7.147  18.218   (  -3.570    4.588   -4.849)    7.570  19.755   (  -6.661    9.916   -8.867)   14.877======================= Grid point 97 (50/126) =======================q-point: ( 0.47 -0.07 -0.29)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.418   (  -8.578    8.323    6.759)   13.731   2.148   ( -10.635  -10.160    2.996)   15.010   2.739   (  -2.194    5.180    2.431)    6.128   3.074   (  -2.278    0.005    7.092)    7.449   3.322   (   4.664   -6.302   -4.235)    8.911   3.845   ( -10.476    3.852    5.125)   12.282   4.210   (   1.528   -9.131    0.094)    9.258   4.992   (   0.818    4.038   -7.024)    8.143   5.253   (  -9.556    8.308   -8.122)   15.043   5.802   (   5.600    1.303   -0.992)    5.834   5.977   (  -5.267    9.114   -9.566)   14.224   6.565   (   7.276    7.972    6.116)   12.406   7.128   (   2.778   11.084    1.441)   11.517   8.175   ( -28.402   -3.220    5.734)   29.154  10.420   (  14.433   -6.421   -1.872)   15.907  11.181   (   5.166   -3.270   -4.429)    7.550  18.116   (   0.970   -1.576    0.107)    1.854  19.419   (   3.204   -2.767    9.636)   10.525======================= Grid point 98 (51/126) =======================q-point: ( 0.47 -0.07 -0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.425   ( -10.290   10.290    8.154)   16.681   1.887   (   1.277   -1.277    1.192)    2.164   2.786   (  -9.236    9.236    5.408)   14.137   3.039   (   0.565   -0.565    9.316)    9.350   3.285   (   6.801   -6.801   -4.755)   10.729   3.796   (  -6.368    6.368    4.343)    9.999   4.097   (   6.694   -6.694    1.512)    9.587   5.060   (  -3.653    3.653   -6.988)    8.690   5.248   ( -11.463   11.463  -10.707)   19.428   5.874   (   0.935   -0.935   -1.447)    1.960   5.999   (  -4.431    4.431  -14.837)   16.106   6.841   (   1.034   -1.034    7.835)    7.971   7.273   (  -5.556    5.556    8.111)   11.293   7.745   ( -15.531   15.531    1.007)   21.987  10.493   (   8.760   -8.760    1.508)   12.480  11.213   (   6.808   -6.808   -7.844)   12.418  18.110   (   1.664   -1.664    0.408)    2.389  19.427   (   4.078   -4.078   11.647)   12.997======================= Grid point 101 (52/126) =======================q-point: ( 0.28 -0.08  0.08)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.132   (   0.000   15.187    1.835)   15.297   2.199   (   0.000   20.730    3.081)   20.958   2.665   (   0.000    0.954    7.267)    7.329   3.039   (  -0.000    2.725   -1.050)    2.920   3.383   (  -0.000    1.670   -2.190)    2.754   3.661   (   0.000   13.175    7.820)   15.321   4.496   (  -0.000    1.061   -9.638)    9.696   5.079   (  -0.000    8.990   -0.994)    9.045   5.309   (   0.000  -11.884    4.536)   12.720   5.404   (  -0.000    6.605   -8.431)   10.710   5.952   (  -0.000  -11.548   -1.788)   11.686   6.490   (   0.000   -0.917    2.164)    2.350   6.793   (  -0.000   -9.750   -1.327)    9.840   8.154   (   0.000   32.387    6.195)   32.974  10.765   (  -0.000  -22.765   -6.645)   23.714  11.181   (   0.000   -0.156    0.022)    0.158  18.177   (  -0.000   -0.846   -2.537)    2.675  19.255   (   0.000    2.378   11.478)   11.722======================= Grid point 102 (53/126) =======================q-point: ( 0.28 -0.08  0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.506   (   5.638   14.202   -5.255)   16.158   2.475   (  -7.419   10.272    4.428)   13.422   2.734   (   0.902    2.360   -1.684)    3.036   2.999   (  -1.971    1.027    7.256)    7.589   3.356   (  -0.539   -4.850   -3.197)    5.834   3.903   (   2.134    9.018    4.641)   10.365   4.337   (  -5.689  -10.034    0.241)   11.536   5.053   (  -6.626   -0.294   -3.396)    7.452   5.332   (   7.359    3.004    0.260)    7.952   5.601   (   1.607   -6.884   -1.837)    7.304   5.815   (  11.612   -1.822    0.011)   11.754   6.472   (   1.134    1.454   -3.998)    4.403   6.748   (  10.533   -5.049    2.980)   12.055   8.723   (  -3.985   23.393    3.617)   24.004  10.244   (  -3.824  -18.366   -3.423)   19.070  11.179   (  -0.501   -1.619   -0.471)    1.759  18.159   (  -1.039    0.217   -1.729)    2.029  19.340   (  -0.833    6.295    4.309)    7.674======================= Grid point 103 (54/126) =======================q-point: ( 0.28 -0.08  0.26)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.862   (   4.719    6.501  -12.165)   14.578   2.445   ( -10.325   -4.690    1.876)   11.494   2.615   (  -3.345    4.410    0.522)    5.560   2.971   (  -3.040    1.964    6.158)    7.143   3.346   (   1.321    0.617   -0.914)    1.721   3.961   (  -6.465   -7.399    3.484)   10.425   4.158   (   4.203    4.200    1.723)    6.187   4.898   (  -3.650   -5.051   -0.859)    6.291   5.359   (   2.843   -7.329   -1.187)    7.950   5.585   (   6.826   -0.334    3.009)    7.467   5.968   (  10.492  -10.764    2.336)   15.212   6.684   (  10.026    6.501   -8.327)   14.564   7.008   (  15.305    2.658    4.740)   16.242   8.767   (  -7.171   -9.427   -0.900)   11.879  10.067   (  -0.442    7.380    0.494)    7.410  11.074   (  -6.475   -1.322    0.192)    6.612  18.145   (  -2.096    1.590   -1.654)    3.108  19.460   (  -2.597    8.154   -0.092)    8.558======================= Grid point 104 (55/126) =======================q-point: (-0.72 -0.08  0.35)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.954   (  -1.100   -0.108   -9.126)    9.193   2.032   (  -5.314   -7.098   -3.604)    9.571   2.789   (  -3.750    9.128    3.369)   10.428   3.062   (   5.222    4.698    2.916)    7.605   3.258   (  -2.146   -6.418   -1.172)    6.868   3.828   (  -2.044    2.776    3.356)    4.811   4.178   (  -7.479    1.179   -5.669)    9.458   4.778   (   3.760   -5.090    8.770)   10.815   5.370   (   7.305   -3.633   -3.514)    8.883   5.661   (   4.456   -4.411    5.545)    8.370   5.901   (  13.009  -15.969    3.313)   20.862   7.003   (   6.407    3.590   -4.600)    8.666   7.410   (  13.231    5.062    0.040)   14.166   8.274   (   1.314  -24.431   -1.428)   24.508  10.325   (  -1.507   13.325   -0.345)   13.415  10.883   ( -11.210    2.532    2.397)   11.739  18.138   (  -2.844    2.709   -1.924)    4.374  19.547   (  -5.276    7.614   -2.039)    9.485======================= Grid point 110 (56/126) =======================q-point: ( 0.28 -0.08 -0.10)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.817   (  -9.713    9.713   11.027)   17.615   1.681   ( -13.041   13.041    1.299)   18.488   2.530   (   2.180   -2.180   -0.188)    3.088   2.642   ( -23.113   23.113    6.565)   33.339   2.903   (   1.647   -1.647   14.109)   14.300   3.844   (  -4.909    4.909   -4.983)    8.546   4.317   (   3.746   -3.746   -1.714)    5.568   4.983   (  -4.708    4.708   -2.773)    7.213   5.026   (  -4.852    4.852  -14.521)   16.061   6.126   (   4.641   -4.641   -3.841)    7.604   6.165   (  -1.808    1.808  -14.354)   14.580   6.651   (  -1.485    1.485   13.551)   13.713   6.798   (  -6.234    6.234   13.111)   15.799   7.307   (  -5.375    5.375   -2.991)    8.169  10.835   (   7.917   -7.917    0.360)   11.202  11.594   (   3.482   -3.482  -10.073)   11.213  18.203   (   1.957   -1.957   -2.809)    3.943  19.259   (   3.703   -3.703   23.946)   24.512======================= Grid point 112 (57/126) =======================q-point: ( 0.09 -0.09  0.09)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.435   (   0.000   14.233    0.000)   14.233   1.172   (   0.000   26.688    0.000)   26.688   1.953   (   0.000   38.408    0.000)   38.408   2.760   (  -0.000   -1.590   -0.000)    1.590   3.114   (   0.000    9.246    0.000)    9.246   3.247   (   0.000    4.228    0.000)    4.228   4.544   (  -0.000   -0.215   -0.000)    0.215   4.741   (   0.000    6.120    0.000)    6.120   5.390   (   0.000    0.456    0.000)    0.456   5.841   (  -0.000  -11.064   -0.000)   11.064   6.369   (  -0.000   -6.238   -0.000)    6.238   6.503   (  -0.000   -1.066   -0.000)    1.066   6.740   (   0.000   25.876    0.000)   25.876   7.158   (  -0.000   -4.788   -0.000)    4.788  11.189   (  -0.000   -0.196   -0.000)    0.196  11.640   (  -0.000  -13.691   -0.000)   13.691  18.323   (  -0.000   -2.784   -0.000)    2.784  18.922   (   0.000    3.035    0.000)    3.035======================= Grid point 113 (58/126) =======================q-point: ( 0.09 -0.09  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.836   (   7.105   17.149    0.000)   18.562   1.793   (   7.009   19.706    0.000)   20.915   2.515   (   2.346    5.564    0.000)    6.038   2.862   (  -5.938    5.407    0.000)    8.031   2.999   (  -5.070   22.784    0.000)   23.341   3.691   (  14.062   -0.988   -0.000)   14.096   4.454   (  -8.231    0.907    0.000)    8.280   4.995   (   7.432    8.210    0.000)   11.074   5.295   (  -5.425    1.141    0.000)    5.544   5.792   (  15.709  -10.939   -0.000)   19.143   6.165   (  -2.106   -9.650   -0.000)    9.878   6.457   (  -1.063   -0.059    0.000)    1.064   7.033   (   0.783   -1.814   -0.000)    1.976   7.401   (   2.107   27.747    0.000)   27.827  10.932   ( -13.920   -6.890    0.000)   15.532  11.471   (   7.051  -12.707   -0.000)   14.532  18.249   (  -1.239   -2.728   -0.000)    2.996  19.012   (   0.725    4.676    0.000)    4.732======================= Grid point 114 (59/126) =======================q-point: ( 0.09 -0.09  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.425   (  13.576   17.807    0.000)   22.392   2.051   (  -8.651    7.543    0.000)   11.477   2.702   (   3.905    6.701    0.000)    7.756   2.768   (  -4.147    3.568    0.000)    5.470   3.430   (  -3.137   -0.451    0.000)    3.169   3.723   (   4.568    6.206    0.000)    7.706   4.207   ( -13.649   -3.143    0.000)   14.006   5.155   (   3.186   -1.150   -0.000)    3.387   5.520   (  11.648    5.733   -0.000)   12.982   5.741   (  -2.386  -11.163   -0.000)   11.415   6.121   (  14.777   -3.533   -0.000)   15.193   6.550   (   6.437    4.793   -0.000)    8.025   7.027   (   6.158   -3.676   -0.000)    7.171   8.044   (   0.806   29.727    0.000)   29.738  10.458   ( -10.516  -12.685   -0.000)   16.477  11.300   (  -1.717   -9.815   -0.000)    9.964  18.169   (  -1.844   -1.727    0.000)    2.527  19.137   (   0.554    5.822    0.000)    5.848======================= Grid point 115 (60/126) =======================q-point: ( 0.09 -0.09  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.878   (  -5.927    5.579    0.000)    8.140   2.146   (   6.673   13.152    0.000)   14.748   2.710   (  -1.814   -2.961   -0.000)    3.473   2.746   (  -3.338   -4.589   -0.000)    5.674   3.229   (  -5.373    5.884    0.000)    7.968   3.814   ( -13.348   -4.480    0.000)   14.080   4.077   (   4.296    7.755    0.000)    8.865   5.133   (   2.467  -11.067   -0.000)   11.338   5.428   (  -4.262   -1.900    0.000)    4.666   5.801   (  10.308   -4.030   -0.000)   11.067   6.382   (  14.526   -3.434   -0.000)   14.926   6.835   (   6.125    7.024    0.000)    9.320   7.211   (  12.534    1.910   -0.000)   12.679   8.618   (   5.846   16.828    0.000)   17.815  10.049   (  -5.259   -7.844   -0.000)    9.444  10.984   ( -11.099   -7.801    0.000)   13.566  18.106   (  -1.991   -0.420    0.000)    2.035  19.262   (  -0.666    5.942    0.000)    5.979======================= Grid point 116 (61/126) =======================q-point: ( 0.09 -0.09  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.828   (  -9.113    7.627    0.000)   11.884   2.076   (  -4.264   -8.474   -0.000)    9.486   2.705   (  -1.059    1.526    0.000)    1.857   2.951   (  10.439    7.237   -0.000)   12.702   3.220   (  -6.503    3.263    0.000)    7.276   3.614   (  -5.066    4.861    0.000)    7.021   4.100   (  -5.767   -2.285    0.000)    6.204   4.924   (   2.207   -7.693   -0.000)    8.003   5.434   (   0.047    1.946    0.000)    1.946   5.861   (   7.531   -6.813   -0.000)   10.155   6.536   (  14.829  -11.218   -0.000)   18.594   7.073   (   2.916    5.215    0.000)    5.975   7.456   (   1.869    4.389    0.000)    4.770   8.883   (  17.083  -13.843   -0.000)   21.987  10.060   (   1.037   13.076    0.000)   13.118  10.536   ( -19.403   -3.343    0.000)   19.689  18.071   (  -1.662    0.798    0.000)    1.844  19.341   (  -2.701    4.769    0.000)    5.480======================= Grid point 123 (62/126) =======================q-point: (-0.10 -0.10  0.10)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 105Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.094   (  -0.000    0.000   -5.089)    5.089   0.094   (  -0.000    0.000   -5.089)    5.089   0.806   (  -0.000    0.000  -42.231)   42.231   2.840   (   0.000   -0.000    5.030)    5.030   2.840   (   0.000   -0.000    5.030)    5.030   4.229   (   0.000   -0.000    1.669)    1.669   4.499   (   0.000   -0.000    6.294)    6.294   4.499   (   0.000   -0.000    6.294)    6.294   4.618   (   0.000   -0.000    0.244)    0.244   6.008   (   0.000   -0.000    5.407)    5.407   6.008   (   0.000   -0.000    5.407)    5.407   6.633   (  -0.000    0.000   -9.097)    9.097   6.633   (  -0.000    0.000   -9.097)    9.097   7.280   (   0.000   -0.000    8.655)    8.655  11.194   (  -0.000    0.000   -0.056)    0.056  11.194   (  -0.000    0.000   -0.056)    0.056  18.440   (  -0.000    0.000   -5.074)    5.074  21.168   (   0.000   -0.000    4.919)    4.919======================= Grid point 126 (63/126) =======================q-point: (-0.10 -0.10  0.37)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.360   (  15.563   15.563    0.703)   22.021   1.864   (  -2.255   -2.255   -1.021)    3.349   2.702   (   8.590    8.590   -3.161)   12.553   2.814   (  -0.643   -0.643  -12.275)   12.309   3.351   (   0.889    0.889    1.752)    2.157   3.746   (   1.238    1.238   -0.528)    1.829   4.078   (  -7.769   -7.769   -0.550)   11.000   5.160   (   3.818    3.818    3.502)    6.435   5.457   (   8.977    8.977    9.825)   16.053   5.875   (  -3.959   -3.959   -1.153)    5.717   6.276   (   4.266    4.266   13.144)   14.462   6.689   (   0.613    0.613   -7.039)    7.092   7.107   (   7.598    7.598   -8.696)   13.823   7.733   (  14.963   14.963   -0.364)   21.165  10.469   (  -9.288   -9.288   -0.951)   13.169  11.348   (  -8.188   -8.188    5.675)   12.896  18.134   (  -1.521   -1.521    2.531)    3.322  19.200   (   0.150    0.150  -10.862)   10.864======================= Grid point 127 (64/126) =======================q-point: (-0.10 -0.10  0.46)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.733   (  -2.809   -2.809   -3.111)    5.046   2.009   (  13.324   13.324    7.412)   20.248   2.781   (  -5.314   -5.314   -6.482)    9.924   2.915   (  -6.510   -6.510   -1.468)    9.322   3.009   (   4.875    4.875   -2.766)    7.429   3.723   (  -8.124   -8.124   -1.747)   11.621   4.041   (   4.667    4.667    0.749)    6.642   5.297   (   2.710    2.710    3.295)    5.054   5.321   (  -8.204   -8.204    1.823)   11.745   6.039   (   4.712    4.712   -4.992)    8.326   6.534   (   7.747    7.747   11.097)   15.595   6.774   (   3.717    3.717   -2.164)    5.684   7.356   (   3.769    3.769   -4.008)    6.670   8.493   (  20.674   20.674   -0.360)   29.239  10.103   (  -7.767   -7.767   -2.162)   11.195  10.915   ( -12.691  -12.691    3.530)   18.292  18.075   (  -1.246   -1.246    1.303)    2.192  19.234   (   1.156    1.156   -3.957)    4.282======================= Grid point 128 (65/126) =======================q-point: (-0.10 -0.10  0.55)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.662   (  -0.697   -0.697   -4.728)    4.830   2.092   (  -5.476   -5.476    4.636)    9.025   2.652   (   1.596    1.596    1.802)    2.888   3.070   (  -0.924   -0.924   -0.443)    1.380   3.074   (   7.644    7.644   -4.053)   11.545   3.516   (  -1.609   -1.609   -3.614)    4.270   4.003   (  -4.714   -4.714    6.932)    9.618   5.081   (  -2.611   -2.611   -5.087)    6.285   5.376   (   1.018    1.018    2.981)    3.311   6.133   (   0.707    0.707   -8.633)    8.691   6.840   (   5.390    5.390    9.577)   12.240   6.920   (   2.225    2.225   -0.071)    3.147   7.408   (  -0.599   -0.599   -0.529)    0.999   9.348   (  18.483   18.483    0.205)   26.139   9.869   (  -3.055   -3.055   -3.560)    5.598  10.315   ( -15.005  -15.005    2.522)   21.370  18.038   (  -0.499   -0.499    0.848)    1.103  19.274   (   0.580    0.580   -1.024)    1.312======================= Grid point 134 (66/126) =======================q-point: (-0.29 -0.11  0.11)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.136   (  -0.000   -1.372  -24.450)   24.488   1.272   (  -0.000  -27.906   -5.375)   28.419   2.210   (  -0.000  -32.328  -13.223)   34.928   2.862   (   0.000   -4.970    1.134)    5.098   2.956   (   0.000  -10.954    8.286)   13.735   3.543   (  -0.000    7.154   -9.553)   11.935   4.273   (   0.000    2.665   11.715)   12.015   4.619   (   0.000   -5.462    9.654)   11.092   4.915   (   0.000   -8.990   10.345)   13.705   5.798   (   0.000    6.089    0.652)    6.124   6.146   (   0.000   -0.359    8.378)    8.386   6.712   (  -0.000    4.441   -7.782)    8.960   6.971   (   0.000    1.343    9.265)    9.362   7.202   (  -0.000  -17.548  -11.369)   20.909  11.190   (  -0.000    0.212   -0.071)    0.224  11.200   (   0.000    5.690   10.342)   11.804  18.336   (  -0.000    7.719   -6.988)   10.412  20.005   (  -0.000   22.474  -21.998)   31.448======================= Grid point 135 (67/126) =======================q-point: (-0.29 -0.11  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.974   (  21.185  -21.185   -2.710)   30.082   1.032   (   8.483   -8.483  -15.844)   19.873   1.841   (  19.288  -19.288  -27.371)   38.642   2.619   (  -3.330    3.330    7.139)    8.553   2.989   (  10.639  -10.639    7.059)   16.619   3.702   (  -8.458    8.458   -8.796)   14.848   4.312   (  -0.737    0.737    7.282)    7.356   4.589   (   2.397   -2.397   19.505)   19.797   4.728   (   4.401   -4.401    2.631)    6.757   5.915   (  -0.216    0.216    6.578)    6.585   6.066   (   3.140   -3.140    4.579)    6.378   6.864   (   2.397   -2.397  -10.605)   11.134   6.869   (   1.783   -1.783  -10.306)   10.610   7.043   (   0.904   -0.904   12.751)   12.815  11.101   (  -4.430    4.430   -0.192)    6.268  11.329   (   3.933   -3.933    8.258)    9.956  18.434   (  -5.783    5.783   -8.423)   11.740  20.327   ( -21.424   21.424  -18.563)   35.533======================= Grid point 139 (68/126) =======================q-point: (-0.29 -0.11  0.56)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.751   (   3.369    7.481    1.632)    8.366   2.175   (  13.173   -8.498   -2.622)   15.894   2.605   ( -14.242   -5.081    2.479)   15.323   3.011   (  -0.295    0.757   -3.085)    3.190   3.283   (   6.101   -3.611   -3.534)    7.922   3.890   (  -1.943  -11.216   -3.000)   11.771   4.052   (   2.234   -2.922    1.736)    4.067   5.075   (  -3.212    7.305    2.434)    8.343   5.346   (  -0.277   -2.734    1.993)    3.394   5.928   (  -1.035   13.601   -1.777)   13.755   6.297   (   3.956    7.466    4.996)    9.816   6.672   (  -0.808   10.910    0.065)   10.940   7.339   (   0.005   10.308   -2.488)   10.604   8.646   (  14.841    3.881   -1.202)   15.387  10.065   (  -8.305   -2.492   -0.390)    8.680  10.926   (  -3.992  -11.814    2.161)   12.656  18.071   (  -0.199   -1.659    0.574)    1.767  19.348   (   0.640   -2.077   -2.890)    3.616======================= Grid point 140 (69/126) =======================q-point: ( 0.71 -0.11 -0.35)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.861   (   1.643    1.360    0.176)    2.140   2.015   (  -5.323   -8.967    0.579)   10.444   2.543   (   4.493    2.592    1.694)    5.457   3.112   (   2.782    4.060    0.439)    4.941   3.213   (   0.309   -6.037    0.423)    6.059   3.757   (   4.532   -4.322   -5.354)    8.239   3.886   (  -2.576   -8.874    5.472)   10.739   5.123   (  -2.992    3.594   -6.854)    8.298   5.330   (   1.667   -0.793    4.521)    4.883   6.136   (  -5.404   11.128  -11.643)   16.988   6.428   (  -7.216    9.842    5.701)   13.470   6.908   (   3.989    5.007    4.874)    8.046   7.466   (   1.473   -0.077    0.357)    1.517   8.893   ( -13.570   18.557   -0.295)   22.991  10.053   (   8.389   -0.353    1.283)    8.494  10.566   (   2.479  -18.709   -2.063)   18.985  18.042   (   0.511   -1.314    0.675)    1.563  19.337   (   2.143   -2.164    0.783)    3.144======================= Grid point 145 (70/126) =======================q-point: (-0.48 -0.12  0.12)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.568   (  -0.000    1.618  -18.493)   18.563   2.433   (  -0.000  -18.700   -8.904)   20.712   2.724   (   0.000   -0.008    2.407)    2.407   3.119   (   0.000    0.163    0.905)    0.920   3.261   (   0.000   -4.477    3.520)    5.695   3.960   (   0.000   -5.774   15.862)   16.880   3.982   (  -0.000   -7.038   -6.099)    9.314   4.857   (   0.000  -17.930    5.562)   18.773   5.214   (   0.000  -11.255    2.698)   11.574   5.509   (  -0.000    9.826   -4.642)   10.867   5.745   (   0.000    8.516    6.303)   10.595   6.598   (  -0.000    4.805   -2.676)    5.499   6.683   (   0.000    4.642    4.576)    6.518   8.457   (  -0.000  -28.076   -7.113)   28.963  10.467   (   0.000   19.247    6.248)   20.236  11.182   (  -0.000    0.192   -0.031)    0.195  18.202   (  -0.000    3.551   -3.732)    5.151  19.668   (  -0.000   10.025   -6.840)   12.136======================= Grid point 146 (71/126) =======================q-point: (-0.48 -0.12  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.577   (  -0.673   -0.642  -21.712)   21.732   2.013   (   1.186  -21.297   -6.589)   22.325   2.589   (  -2.199   -6.779    2.305)    7.490   3.113   (   1.413   -4.080   -3.452)    5.528   3.189   (   8.149   -8.033    5.258)   12.593   3.745   (   0.916   -7.381   -5.702)    9.372   4.079   (   9.358   -0.519    8.572)   12.701   4.560   (   1.760  -13.146    9.697)   16.430   4.934   (  -4.125   -7.425    3.720)    9.273   5.698   (   0.747    3.494    2.676)    4.464   5.887   (  -1.920    6.408    7.310)    9.909   6.719   (   4.279    6.713   -2.397)    8.314   6.901   (  10.636    1.047    4.151)   11.465   7.832   (  -4.393  -26.135   -8.364)   27.790  10.775   (  -4.329   13.347    3.832)   14.545  11.196   (   0.229    1.467    2.415)    2.835  18.285   (  -2.572    6.432   -6.151)    9.264  19.874   (  -4.357   12.965   -9.220)   16.495======================= Grid point 147 (72/126) =======================q-point: (-0.48 -0.12  0.30)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.552   (   1.375   -1.375  -21.258)   21.347   1.749   (  10.593  -10.593   -1.684)   15.075   2.503   (   5.494   -5.494    1.296)    7.877   3.034   (  10.894  -10.894  -11.508)   19.229   3.223   (   5.225   -5.225    5.193)    9.031   3.711   (   3.370   -3.370   -2.836)    5.545   4.209   (   0.677   -0.677    3.242)    3.380   4.408   (   9.766   -9.766   12.413)   18.569   4.801   (   2.347   -2.347    5.664)    6.564   5.722   (  -1.856    1.856    7.559)    8.002   5.934   (  -3.909    3.909    6.675)    8.667   6.891   (   3.600   -3.600    6.356)    8.144   7.206   (   4.909   -4.909   -7.500)   10.220   7.285   (   8.134   -8.134   -6.438)   13.183  10.855   (  -7.411    7.411   -0.566)   10.495  11.227   (  -1.595    1.595    6.377)    6.764  18.331   (  -6.026    6.026   -7.139)   11.118  19.973   ( -10.714   10.714   -9.487)   17.877======================= Grid point 152 (73/126) =======================q-point: ( 0.52 -0.12 -0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.634   (  -8.377   11.932    1.916)   14.705   1.941   (  -7.387   -8.506    2.850)   11.621   2.862   (   6.133    5.014    2.641)    8.351   3.056   (   2.420   -2.313    3.938)    5.168   3.136   (  -2.714   -8.935    1.328)    9.432   3.928   (  -2.539    0.983    0.578)    2.784   3.994   (   0.784   -8.045    3.613)    8.853   5.108   (  -0.542    5.770   -6.802)    8.936   5.359   (  -0.158    0.373   -4.980)    4.996   5.911   (  -2.640   10.127   -0.346)   10.471   6.133   (  -8.883    6.450  -12.376)   16.543   6.752   (   5.507    7.720    6.481)   11.486   7.337   (  -0.858    8.071    4.470)    9.266   8.196   ( -24.837    7.249    2.092)   25.957  10.290   (  11.037   -6.133    1.322)   12.695  11.052   (   7.711   -9.070   -5.610)   13.160  18.083   (   0.915   -1.528   -0.283)    1.803  19.363   (   1.996   -2.500    7.385)    8.048======================= Grid point 156 (74/126) =======================q-point: ( 0.33 -0.13  0.13)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.447   (  -0.000   14.648   -4.212)   15.242   2.578   (   0.000   15.187    2.156)   15.339   2.688   (   0.000    1.229    2.407)    2.703   3.034   (   0.000   -1.875    3.937)    4.361   3.349   (  -0.000   -4.061   -2.883)    4.980   3.889   (   0.000    8.348    5.099)    9.782   4.387   (  -0.000  -13.528   -1.400)   13.601   5.205   (   0.000    3.378    0.522)    3.418   5.226   (   0.000    4.776    2.009)    5.182   5.507   (  -0.000    0.091   -7.221)    7.221   5.720   (  -0.000   -9.345   -0.416)    9.355   6.489   (   0.000    0.950    0.850)    1.274   6.609   (  -0.000   -6.621   -0.296)    6.628   8.768   (   0.000   24.457    4.134)   24.804  10.288   (  -0.000  -21.011   -4.441)   21.475  11.179   (   0.000   -0.041    0.008)    0.042  18.170   (  -0.000    0.305   -1.805)    1.831  19.349   (   0.000    6.532    4.578)    7.976======================= Grid point 157 (75/126) =======================q-point: ( 0.33 -0.13  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.764   (   3.836   10.250  -11.129)   15.608   2.622   (  -6.357    2.486    1.505)    6.990   2.716   (  -5.572   -3.182   -1.976)    6.714   3.016   (   0.005    0.397    7.697)    7.708   3.301   (   2.144   -0.301   -2.636)    3.411   4.003   (   1.958   -7.811    3.594)    8.819   4.065   (   2.681   -4.941    5.749)    8.040   5.076   ( -11.838   -0.177   -3.049)   12.226   5.278   (   3.455  -10.580    0.362)   11.136   5.498   (   2.519    1.016    1.022)    2.902   5.750   (  10.813   -5.086    0.679)   11.968   6.528   (   3.760    4.266   -5.096)    7.636   6.725   (  11.491    2.160    4.089)   12.387   8.915   (  -5.963   -8.553   -0.710)   10.451  10.075   (  -0.301    5.948    0.334)    5.965  11.160   (  -2.147   -0.389    0.048)    2.182  18.181   (  -1.311    2.108   -1.952)    3.157  19.504   (  -1.530    9.433    0.005)    9.556======================= Grid point 158 (76/126) =======================q-point: (-0.67 -0.13  0.31)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.949   (   1.185    2.221  -16.157)   16.352   2.216   ( -11.711  -12.455   -0.899)   17.120   2.791   (   3.603    6.667    2.206)    7.892   3.003   (   0.082   -0.295    6.188)    6.196   3.293   (  -0.642   -4.053   -0.869)    4.195   3.847   (   0.712   -2.484    3.423)    4.289   4.208   (   5.224    0.667    0.884)    5.340   4.805   (  -6.902   -4.359    1.965)    8.396   5.138   (  12.280  -15.176   -1.784)   19.603   5.596   (   2.157    0.467    3.512)    4.147   5.738   (   3.560   -7.970    5.184)   10.152   6.825   (   9.993    6.179   -5.652)   13.038   7.076   (  16.169    3.681    3.108)   16.871   8.363   (  -7.157  -25.244   -2.947)   26.404  10.351   (  -1.155   16.266    1.160)   16.348  11.068   (  -6.493    0.578    1.328)    6.653  18.201   (  -2.989    3.901   -2.607)    5.563  19.650   (  -4.484    9.989   -2.406)   11.211======================= Grid point 159 (77/126) =======================q-point: (-0.67 -0.13  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.911   (   0.805   -0.805   -0.088)    1.142   1.979   (   0.627   -0.627  -16.674)   16.698   2.933   (  -4.626    4.626    9.665)   11.671   2.989   (  10.919  -10.919   -3.691)   15.877   3.236   (  -3.774    3.774    1.411)    5.521   3.862   (  -0.001    0.001    1.038)    1.038   4.218   (  -1.962    1.962   -4.702)    5.460   4.706   (   1.456   -1.456    9.454)    9.676   5.133   (  19.773  -19.773   -4.689)   28.354   5.578   (   2.225   -2.225    5.837)    6.631   5.703   (   2.401   -2.401    6.380)    7.228   7.008   (   3.273   -3.273   -1.215)    4.785   7.483   (   3.024   -3.024   -2.751)    5.085   7.815   (  13.593  -13.593   -1.838)   19.311  10.546   (  -7.657    7.657   -0.892)   10.866  10.976   (  -5.893    5.893    3.487)    9.034  18.212   (  -4.269    4.269   -2.938)    6.714  19.711   (  -7.948    7.948   -3.105)   11.661======================= Grid point 167 (78/126) =======================q-point: ( 0.14 -0.14  0.14)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.764   (   0.000   17.101    0.000)   17.101   1.701   (   0.000   24.140    0.000)   24.140   2.482   (   0.000    9.738    0.000)    9.738   2.885   (   0.000   14.397    0.000)   14.397   3.301   (   0.000    7.826    0.000)    7.826   3.342   (   0.000    6.347    0.000)    6.347   4.557   (   0.000    1.714    0.000)    1.714   4.900   (   0.000    8.871    0.000)    8.871   5.417   (   0.000    2.397    0.000)    2.397   5.566   (  -0.000  -15.026   -0.000)   15.026   6.192   (  -0.000  -10.110   -0.000)   10.110   6.481   (  -0.000   -0.957   -0.000)    0.957   6.998   (  -0.000   -7.450   -0.000)    7.450   7.393   (   0.000   32.361    0.000)   32.361  11.185   (  -0.000   -0.209   -0.000)    0.209  11.291   (  -0.000  -19.693   -0.000)   19.693  18.262   (  -0.000   -2.827   -0.000)    2.827  19.004   (   0.000    4.913    0.000)    4.913======================= Grid point 168 (79/126) =======================q-point: ( 0.14 -0.14  0.23)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.206   (   7.181   18.092    0.000)   19.465   2.162   (  -1.911   15.423    0.000)   15.541   2.612   (   3.446    5.742    0.000)    6.696   2.929   (  -9.638    1.455    0.000)    9.747   3.416   (   1.739   10.069    0.000)   10.218   3.655   (   4.245    2.747   -0.000)    5.057   4.464   ( -10.522   -1.105    0.000)   10.580   5.130   (   0.615    1.609    0.000)    1.723   5.291   (   7.430   -0.576   -0.000)    7.453   5.656   (  10.895   -0.420   -0.000)   10.903   5.950   (   0.549  -10.534   -0.000)   10.549   6.470   (   1.262    1.525    0.000)    1.980   6.890   (   6.569   -9.099   -0.000)   11.223   8.067   (  -1.907   32.811    0.000)   32.867  10.686   ( -10.837  -17.177   -0.000)   20.310  11.263   (   4.663   -6.460   -0.000)    7.967  18.200   (  -1.104   -1.819   -0.000)    2.128  19.130   (   0.172    6.761    0.000)    6.763======================= Grid point 169 (80/126) =======================q-point: ( 0.14 -0.14  0.32)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.768   (  11.090   14.538    0.000)   18.285   2.228   ( -10.032    9.717    0.000)   13.967   2.776   (  -4.713   -3.769    0.000)    6.035   2.811   (  -1.288    2.839    0.000)    3.118   3.334   (  -4.463   -1.686    0.000)    4.771   3.948   (   5.400    9.336    0.000)   10.785   4.121   ( -12.770   -4.971    0.000)   13.703   5.061   (   2.492   -7.185   -0.000)    7.605   5.519   (   0.588   -6.186   -0.000)    6.214   5.579   (   6.238   -2.649   -0.000)    6.777   6.085   (  14.036   -2.190   -0.000)   14.205   6.675   (   8.618    6.570   -0.000)   10.837   6.930   (  12.844   -3.713   -0.000)   13.370   8.607   (  -3.116   21.495    0.000)   21.720  10.165   (  -5.797  -13.379   -0.000)   14.581  11.143   (  -4.423   -5.231   -0.000)    6.850  18.147   (  -1.870   -0.255    0.000)    1.887  19.278   (  -0.835    7.792    0.000)    7.836======================= Grid point 170 (81/126) =======================q-point: ( 0.14 -0.14  0.41)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.003   (  -5.998    5.979    0.000)    8.469   2.208   (  -8.133   -7.462    0.000)   11.038   2.664   (   5.866    4.549   -0.000)    7.423   2.821   (  -0.791    7.002    0.000)    7.046   3.331   (  -3.430    1.965    0.000)    3.953   3.756   (  -8.586   -0.560    0.000)    8.604   4.193   (  -1.823    3.612    0.000)    4.046   4.878   (   1.637  -11.809   -0.000)   11.922   5.417   (   1.527   -0.031   -0.000)    1.527   5.703   (   7.436   -4.743   -0.000)    8.820   6.215   (  15.071  -13.100   -0.000)   19.969   6.966   (   7.058    5.580   -0.000)    8.998   7.283   (  14.546    4.925   -0.000)   15.358   8.691   (   1.053  -10.856   -0.000)   10.907  10.054   (  -0.215    9.299    0.000)    9.301  10.873   ( -12.887   -2.879    0.000)   13.205  18.114   (  -2.313    1.224    0.000)    2.617  19.401   (  -2.958    7.398    0.000)    7.967======================= Grid point 171 (82/126) =======================q-point: ( 0.14  0.86 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.964   (   1.719   -1.719   -0.000)    2.432   1.993   (  -7.589    7.589    0.000)   10.733   2.737   (  -1.305    1.305    0.000)    1.845   3.039   (  -1.381    1.381    0.000)    1.953   3.258   (  -0.498    0.498    0.000)    0.705   3.730   (  -4.618    4.618    0.000)    6.531   4.106   (  -4.554    4.554    0.000)    6.440   4.783   (   6.221   -6.221   -0.000)    8.798   5.464   (  -0.531    0.531    0.000)    0.751   5.714   (   6.845   -6.845   -0.000)    9.681   6.229   (  17.542  -17.542   -0.000)   24.808   7.121   (   0.711   -0.711   -0.000)    1.005   7.536   (  -1.510    1.510    0.000)    2.135   8.499   (  18.491  -18.491   -0.000)   26.150  10.285   (  -5.098    5.098    0.000)    7.210  10.603   ( -10.574   10.574    0.000)   14.954  18.102   (  -2.157    2.157    0.000)    3.050  19.449   (  -5.522    5.522    0.000)    7.809======================= Grid point 191 (83/126) =======================q-point: (-0.25 -0.15  0.15)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.674   (  -0.000  -29.176   -4.033)   29.454   0.892   (  -0.000  -29.500   -7.016)   30.323   1.664   (  -0.000  -11.734  -35.748)   37.624   2.756   (   0.000   -4.710    7.790)    9.103   2.769   (   0.000   -6.528    9.304)   11.366   3.823   (  -0.000   18.635   -8.191)   20.355   4.321   (   0.000    1.718   10.214)   10.358   4.524   (   0.000   -5.036   12.065)   13.074   4.712   (   0.000   -9.679    6.718)   11.782   5.889   (   0.000    2.816    2.886)    4.032   6.041   (   0.000   -8.460    6.395)   10.605   6.789   (  -0.000    2.761   -9.151)    9.558   6.895   (  -0.000   -6.776   -9.136)   11.375   7.040   (   0.000   -0.207   12.494)   12.496  11.194   (  -0.000    0.127   -0.082)    0.151  11.231   (   0.000   -1.768    5.464)    5.743  18.500   (  -0.000    6.886   -8.912)   11.262  20.600   (  -0.000   32.373  -11.059)   34.210======================= Grid point 196 (84/126) =======================q-point: (-0.25 -0.15  0.61)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.813   (   2.113   -1.510   -4.104)    4.857   2.062   (   2.897   -3.490    5.007)    6.756   2.548   (  -6.130    1.228    1.957)    6.551   3.067   (   1.225    3.771   -4.399)    5.922   3.175   (   3.335   -5.631   -3.083)    7.234   3.701   (   0.054   -6.754   -6.053)    9.069   3.961   (  -3.469   -5.565    7.726)   10.134   5.190   (  -3.380    3.104   -1.979)    4.997   5.292   (   1.798   -1.744    2.881)    3.817   6.205   (  -0.753   11.533   -7.337)   13.689   6.448   (   5.485    7.448   10.454)   13.958   6.852   (   0.713    6.023   -0.833)    6.122   7.447   (   0.370    0.381   -1.609)    1.694   8.896   (  13.499   19.435    0.006)   23.663  10.014   (  -5.187   -2.517   -2.630)    6.337  10.616   (  -6.982  -17.607    2.996)   19.176  18.041   (  -0.345   -1.231    0.532)    1.384  19.311   (   0.362   -1.373   -1.897)    2.370======================= Grid point 197 (85/126) =======================q-point: ( 0.75 -0.15 -0.30)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.830   (   4.408   -4.408   -7.360)    9.645   1.912   (  -0.723    0.723    6.466)    6.546   2.622   (  -3.234    3.234    1.003)    4.682   3.091   (   2.983   -2.983    0.512)    4.249   3.189   (  -0.451    0.451    0.027)    0.639   3.665   (   4.761   -4.761   -8.128)   10.555   3.774   (   0.367   -0.367    9.081)    9.095   5.170   (  -1.448    1.448   -7.218)    7.503   5.314   (   0.770   -0.770    5.421)    5.530   6.313   (  -5.457    5.457  -14.460)   16.391   6.603   (  -5.527    5.527   11.849)   14.195   6.950   (   0.533   -0.533    2.335)    2.453   7.423   (   2.122   -2.122    0.210)    3.009   9.325   ( -18.594   18.594   -1.003)   26.315  10.051   (   2.053   -2.053    4.562)    5.407  10.163   (  14.973  -14.973   -4.464)   21.641  18.023   (   0.484   -0.484    0.529)    0.865  19.303   (   1.075   -1.075    0.018)    1.520======================= Grid point 202 (86/126) =======================q-point: (-0.44 -0.16  0.16)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.591   (  -0.000    0.355  -22.090)   22.093   1.973   (  -0.000  -25.186   -9.714)   26.994   2.670   (   0.000   -5.971    2.595)    6.510   3.046   (   0.000  -13.911    9.295)   16.731   3.107   (  -0.000   -3.838   -4.094)    5.612   3.741   (  -0.000  -12.829   -9.385)   15.895   3.961   (   0.000    3.872   12.954)   13.520   4.553   (   0.000  -10.390   10.678)   14.898   4.990   (   0.000  -10.110    3.870)   10.826   5.685   (  -0.000    6.884   -2.463)    7.311   5.909   (   0.000    6.758    8.178)   10.609   6.717   (  -0.000    6.360   -4.697)    7.906   6.742   (   0.000    0.736    9.424)    9.452   7.880   (  -0.000  -26.243   -8.643)   27.629  10.826   (   0.000   14.608    6.381)   15.941  11.187   (  -0.000    0.247   -0.050)    0.253  18.313   (  -0.000    7.271   -6.639)    9.846  19.921   (  -0.000   14.340   -9.174)   17.023======================= Grid point 203 (87/126) =======================q-point: (-0.44 -0.16  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.505   (   6.073  -14.727  -12.891)   20.492   1.548   (   4.927  -15.966  -10.996)   20.003   2.464   (  -0.491   -4.126    3.434)    5.390   2.789   (   8.406  -21.240  -14.718)   27.174   3.053   (  10.834   -9.587    8.018)   16.541   3.670   (   0.276    2.590   -3.518)    4.377   4.123   (   5.555    3.523    8.512)   10.757   4.366   (  -2.156   -6.479   11.319)   13.219   4.784   (  -0.523   -6.620    3.991)    7.747   5.769   (  -2.241    3.051    3.640)    5.252   5.989   (  -1.327    2.849    6.218)    6.967   6.780   (   5.181   -1.627   11.272)   12.512   6.953   (   9.649    4.820   -5.384)   12.056   7.357   (  -4.202  -18.974   -8.802)   21.334  10.988   (  -5.513    7.434    1.181)    9.330  11.226   (   1.615    1.140    3.927)    4.396  18.447   (  -4.545    8.762   -8.709)   13.164  20.175   (  -7.883   15.753   -7.930)   19.318======================= Grid point 209 (88/126) =======================q-point: ( 0.56 -0.16 -0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.809   (   1.534   -1.534    4.446)    4.947   1.900   ( -10.297   10.297   -3.275)   14.926   2.764   (  12.308  -12.308   -4.245)   17.916   3.058   (  -1.071    1.071    8.890)    9.018   3.081   (  -3.728    3.728    4.502)    6.933   3.777   (   7.132   -7.132    4.061)   10.873   3.997   (  -0.313    0.313   -4.301)    4.323   5.204   (  -3.042    3.042   -6.427)    7.734   5.306   (   3.647   -3.647    0.824)    5.222   6.138   (  -8.295    8.295    4.725)   12.647   6.275   (  -8.324    8.324  -15.120)   19.162   6.843   (  -1.423    1.423    5.125)    5.506   7.434   (  -1.672    1.672    3.750)    4.434   8.500   ( -20.215   20.215    0.317)   28.589  10.157   (   6.898   -6.898    3.296)   10.297  10.825   (  11.837  -11.837   -5.525)   17.628  18.056   (   1.060   -1.060   -0.470)    1.570  19.322   (   1.324   -1.324    4.708)    5.067======================= Grid point 213 (89/126) =======================q-point: ( 0.37 -0.17  0.17)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.722   (  -0.000   11.466  -10.127)   15.298   2.707   (  -0.000    0.278   -1.653)    1.676   2.771   (  -0.000    1.173   -1.982)    2.303   3.000   (   0.000   -0.654    9.017)    9.040   3.277   (  -0.000   -1.305   -2.740)    3.035   3.998   (   0.000  -21.098    8.427)   22.719   4.013   (   0.000    3.870    2.236)    4.469   5.202   (   0.000  -10.820    1.227)   10.890   5.374   (  -0.000   10.448   -2.874)   10.836   5.376   (  -0.000   -8.455   -1.280)    8.552   5.604   (  -0.000   -2.221   -0.170)    2.227   6.535   (   0.000    3.580   -0.276)    3.591   6.558   (   0.000    1.258    1.257)    1.778   8.981   (  -0.000   -8.341   -0.571)    8.361  10.078   (   0.000    5.312    0.220)    5.317  11.179   (   0.000    0.095   -0.005)    0.095  18.196   (  -0.000    2.328   -2.068)    3.114  19.521   (   0.000    9.889    0.058)    9.889======================= Grid point 214 (90/126) =======================q-point: (-0.63 -0.17  0.26)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.922   (   1.085    4.938  -16.468)   17.227   2.478   ( -12.019  -12.601   -2.754)   17.630   2.705   (   3.409   -2.918   -1.594)    4.762   3.013   (  -0.601    0.477    8.791)    8.824   3.265   (   3.517   -5.638    0.504)    6.664   3.772   (   6.155   -4.569    8.258)   11.268   4.082   (   4.143   -0.376    0.967)    4.271   4.880   (   4.090  -19.402   -2.512)   19.987   5.054   (  -7.156   -8.738    1.412)   11.383   5.572   (   0.133    4.681    2.326)    5.229   5.686   (   2.391    0.306    3.154)    3.969   6.641   (   6.189    5.584   -3.446)    9.020   6.787   (  11.306    3.297    3.582)   12.310   8.502   (  -5.343  -25.997   -3.498)   26.770  10.372   (  -0.783   17.866    2.342)   18.036  11.158   (  -2.495    0.223    0.393)    2.535  18.254   (  -1.964    5.081   -3.159)    6.297  19.726   (  -2.612   11.738   -2.465)   12.275======================= Grid point 215 (91/126) =======================q-point: (-0.63 -0.17  0.35)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.963   (  -1.955   -2.580  -13.914)   14.286   1.976   (  -4.039   -8.487   -4.998)   10.646   2.732   (  12.149  -11.732   -2.177)   17.029   3.007   (  -2.193    1.966   10.944)   11.334   3.266   (   2.283    0.070    0.886)    2.450   3.836   (   2.029    0.213   -0.982)    2.264   4.194   (   5.562   -1.619    0.593)    5.823   4.663   (  13.472  -16.627   -1.338)   21.442   4.776   (  -0.018  -10.131    3.977)   10.884   5.581   (  -2.180   -0.179    4.146)    4.687   5.712   (  -1.955    3.120    6.712)    7.655   6.879   (   8.612   -1.900    0.901)    8.865   7.186   (  17.537    7.164   -1.704)   19.020   7.826   (  -6.703  -24.968   -3.839)   26.135  10.660   (  -3.880   12.613    0.285)   13.199  11.095   (  -4.888    1.916    2.489)    5.810  18.309   (  -4.792    6.416   -4.021)    8.961  19.867   (  -6.866   10.681   -3.333)   13.128======================= Grid point 224 (92/126) =======================q-point: ( 0.18 -0.18  0.18)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.132   (   0.000   18.035    0.000)   18.035   2.170   (   0.000   20.874    0.000)   20.874   2.576   (   0.000    3.172    0.000)    3.172   3.061   (   0.000    0.639    0.000)    0.639   3.405   (   0.000    1.720    0.000)    1.720   3.576   (   0.000   14.338    0.000)   14.338   4.594   (  -0.000   -0.374   -0.000)    0.374   5.083   (   0.000    7.912    0.000)    7.912   5.260   (  -0.000  -11.999   -0.000)   11.999   5.502   (   0.000    5.794    0.000)    5.794   5.959   (  -0.000  -12.205   -0.000)   12.205   6.468   (  -0.000   -0.168   -0.000)    0.168   6.806   (  -0.000  -10.171   -0.000)   10.171   8.091   (   0.000   33.448    0.000)   33.448  10.833   (  -0.000  -23.961   -0.000)   23.961  11.181   (  -0.000   -0.150   -0.000)    0.150  18.212   (  -0.000   -1.853   -0.000)    1.853  19.129   (   0.000    7.082    0.000)    7.082======================= Grid point 225 (93/126) =======================q-point: ( 0.18 -0.18  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.572   (   6.849   16.778    0.000)   18.122   2.426   (  -8.517   10.243    0.000)   13.321   2.764   (   5.142    7.860    0.000)    9.393   2.906   (  -6.672   -3.748    0.000)    7.653   3.390   (  -0.469   -4.727   -0.000)    4.750   3.855   (   2.561   10.055    0.000)   10.376   4.329   (  -6.507  -13.189   -0.000)   14.707   5.095   (  -5.834   -0.470    0.000)    5.853   5.317   (   9.647    1.339   -0.000)    9.740   5.636   (  -2.978   -6.946   -0.000)    7.557   5.802   (  12.422   -1.539   -0.000)   12.517   6.526   (   4.793    3.898   -0.000)    6.177   6.707   (   8.310   -7.217   -0.000)   11.006   8.687   (  -3.501   24.622    0.000)   24.870  10.279   (  -4.721  -19.232   -0.000)   19.803  11.184   (  -0.101   -2.029   -0.000)    2.031  18.179   (  -1.133   -0.084    0.000)    1.136  19.294   (  -0.555    8.985    0.000)    9.002======================= Grid point 226 (94/126) =======================q-point: ( 0.18 -0.18  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.009   (   4.858    8.648    0.000)    9.919   2.407   ( -10.418    2.345    0.000)   10.679   2.621   (  -2.999   -3.699   -0.000)    4.762   2.885   (  -3.566    3.580    0.000)    5.053   3.360   (   0.283    1.393    0.000)    1.421   3.925   (  -6.201   -8.355   -0.000)   10.405   4.137   (   6.876    4.464   -0.000)    8.198   4.909   (  -5.844   -6.603   -0.000)    8.817   5.368   (   3.801   -7.716   -0.000)    8.601   5.561   (   5.465    0.123   -0.000)    5.466   5.944   (  10.866  -11.787   -0.000)   16.031   6.794   (   9.028    5.163   -0.000)   10.400   6.932   (  16.948    3.574   -0.000)   17.321   8.774   (  -7.087   -7.892   -0.000)   10.607  10.063   (  -0.554    6.425    0.000)    6.449  11.073   (  -6.544   -1.873    0.000)    6.807  18.163   (  -2.275    1.852    0.000)    2.933  19.459   (  -2.539    9.573    0.000)    9.904======================= Grid point 227 (95/126) =======================q-point: (-0.82 -0.18  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.010   (  -5.864   -7.579   -0.000)    9.583   2.123   (  -4.595    3.777    0.000)    5.948   2.749   (   0.276    3.027    0.000)    3.040   3.012   (   3.433    9.511    0.000)   10.112   3.273   (  -1.594   -7.002   -0.000)    7.181   3.796   (  -2.192    3.457    0.000)    4.093   4.247   (  -5.223    2.155    0.000)    5.650   4.678   (   2.057   -7.211   -0.000)    7.499   5.408   (   7.547   -1.639   -0.000)    7.723   5.613   (   2.192   -4.137   -0.000)    4.682   5.862   (  15.666  -19.103   -0.000)   24.705   7.049   (   5.959    1.646   -0.000)    6.182   7.407   (  14.714    6.873   -0.000)   16.240   8.288   (   0.308  -24.206   -0.000)   24.208  10.327   (  -0.742   13.862    0.000)   13.882  10.860   ( -11.744    1.638    0.000)   11.858  18.159   (  -3.173    3.127    0.000)    4.455  19.565   (  -5.522    8.295    0.000)    9.964======================= Grid point 248 (96/126) =======================q-point: (-0.20 -0.20  0.20)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.197   (  -0.000    0.000   -5.803)    5.803   0.197   (  -0.000    0.000   -5.803)    5.803   1.602   (  -0.000    0.000  -41.211)   41.211   2.700   (   0.000   -0.000    9.546)    9.546   2.700   (   0.000   -0.000    9.546)    9.546   4.181   (   0.000   -0.000    3.456)    3.456   4.340   (   0.000   -0.000    9.741)    9.741   4.340   (   0.000   -0.000    9.741)    9.741   4.611   (   0.000   -0.000    0.396)    0.396   5.918   (   0.000   -0.000    3.620)    3.620   5.918   (   0.000   -0.000    3.620)    3.620   6.819   (  -0.000    0.000   -9.577)    9.577   6.819   (  -0.000    0.000   -9.577)    9.577   7.045   (   0.000   -0.000   15.647)   15.647  11.195   (  -0.000    0.000   -0.086)    0.086  11.195   (  -0.000    0.000   -0.086)    0.086  18.578   (  -0.000    0.000   -9.144)    9.144  21.034   (   0.000   -0.000    8.872)    8.872======================= Grid point 253 (97/126) =======================q-point: (-0.20 -0.20  0.65)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.779   (  -0.086   -0.086   -7.083)    7.084   1.975   (  -4.436   -4.436    7.327)    9.646   2.622   (   3.365    3.365    0.950)    4.853   3.080   (  -1.935   -1.935   -0.592)    2.800   3.153   (   3.164    3.164   -3.511)    5.687   3.618   (  -0.897   -0.897   -6.993)    7.107   3.840   (  -5.459   -5.459    9.535)   12.269   5.198   (  -1.230   -1.230   -6.553)    6.780   5.294   (   1.151    1.151    5.582)    5.815   6.351   (   2.196    2.196  -13.408)   13.763   6.600   (   5.624    5.624   14.770)   16.775   6.925   (   1.581    1.581   -0.575)    2.309   7.418   (  -1.520   -1.520   -0.405)    2.188   9.341   (  18.499   18.499    0.574)   26.168   9.956   (  -2.539   -2.539   -5.140)    6.269  10.248   ( -14.927  -14.927    4.244)   21.532  18.023   (  -0.445   -0.445    0.497)    0.802  19.294   (  -0.217   -0.217   -0.898)    0.949======================= Grid point 259 (98/126) =======================q-point: (-0.39 -0.21  0.21)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.401   (  -0.000  -29.577   -8.033)   30.648   1.519   (  -0.000  -10.958  -15.974)   19.371   2.462   (  -0.000  -12.639   -7.905)   14.907   2.763   (   0.000  -12.274   11.754)   16.995   2.872   (  -0.000  -14.906   -6.498)   16.261   3.675   (  -0.000    6.484   -4.232)    7.742   4.056   (   0.000    4.474   10.368)   11.292   4.395   (   0.000   -5.862   11.363)   12.786   4.796   (   0.000   -8.449    3.847)    9.283   5.796   (   0.000    3.929   -0.312)    3.942   5.999   (   0.000    1.139    7.037)    7.129   6.723   (   0.000   -1.786   14.713)   14.821   6.847   (  -0.000    5.785   -6.169)    8.458   7.403   (  -0.000  -19.158   -8.844)   21.101  11.061   (   0.000    8.076    4.716)    9.352  11.192   (  -0.000    0.226   -0.067)    0.236  18.498   (  -0.000    9.980   -9.226)   13.591  20.264   (  -0.000   18.414   -6.624)   19.569======================= Grid point 260 (99/126) =======================q-point: (-0.39 -0.21  0.30)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.035   (  21.146  -21.146   -3.614)   30.122   1.259   (  18.578  -18.578  -10.082)   28.141   2.416   (   9.175   -9.175  -27.903)   30.772   2.457   (  -1.090    1.090    9.428)    9.553   2.819   (  11.721  -11.721    7.521)   18.203   3.796   (  -7.276    7.276   -1.073)   10.346   4.168   (   0.730   -0.730    7.188)    7.261   4.276   (   1.388   -1.388   13.084)   13.230   4.680   (   3.060   -3.060    2.343)    4.921   5.823   (  -2.213    2.213    3.261)    4.520   5.992   (   2.291   -2.291    3.428)    4.717   6.724   (   2.311   -2.311   16.881)   17.194   7.039   (   2.823   -2.823   -7.587)    8.574   7.067   (   3.484   -3.484   -7.655)    9.104  11.105   (  -4.326    4.326   -0.186)    6.121  11.229   (   1.006   -1.006    3.065)    3.379  18.617   (  -6.640    6.640  -10.133)   13.815  20.493   ( -13.348   13.348   -1.671)   18.950======================= Grid point 270 (100/126) =======================q-point: (-0.58 -0.22  0.22)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.910   (  -0.000    6.428  -15.770)   17.030   2.611   (  -0.000  -15.042   -8.999)   17.528   2.685   (   0.000   -2.655    1.159)    2.898   3.027   (   0.000    3.392    7.682)    8.397   3.196   (   0.000  -12.414    3.219)   12.824   3.687   (   0.000   -3.537   11.778)   12.298   4.045   (  -0.000   -1.430   -0.933)    1.708   4.828   (  -0.000  -21.928   -1.655)   21.990   5.159   (   0.000  -11.897    2.625)   12.183   5.579   (  -0.000    8.617   -2.603)    9.001   5.646   (   0.000    4.936    3.963)    6.330   6.599   (   0.000    1.262    3.705)    3.914   6.639   (  -0.000    6.365   -1.490)    6.537   8.561   (  -0.000  -26.462   -3.636)   26.710  10.381   (   0.000   18.453    2.834)   18.670  11.183   (   0.000    0.208   -0.018)    0.208  18.276   (  -0.000    5.586   -3.371)    6.524  19.754   (  -0.000   12.393   -2.443)   12.632======================= Grid point 271 (101/126) =======================q-point: (-0.58 -0.22  0.31)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 606Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.966   (  -0.350   -0.991  -18.264)   18.294   2.133   (  -6.773  -21.025   -5.506)   22.765   2.540   (   3.848  -10.396    2.645)   11.396   3.002   (   5.366   -6.235   12.052)   14.592   3.152   (   5.256   -0.375   -0.241)    5.274   3.791   (   2.096    2.223    1.453)    3.383   4.008   (   9.386   -6.349   -0.617)   11.348   4.458   (   4.071  -17.254    0.647)   17.740   4.871   (  -4.650   -8.466    2.925)   10.092   5.657   (  -4.555    3.270    1.365)    5.771   5.754   (  -1.581    5.336    6.330)    8.429   6.684   (   8.855   -2.319    3.771)    9.900   6.893   (   9.898    7.191   -1.680)   12.349   7.954   (  -4.761  -25.392   -4.365)   26.201  10.720   (  -1.985   14.739    1.937)   14.997  11.167   (  -2.100    0.568    0.804)    2.319  18.396   (  -3.294    8.561   -4.836)   10.369  19.984   (  -4.092   12.913   -3.005)   13.875======================= Grid point 272 (102/126) =======================q-point: (-0.58 -0.22  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.773   (   8.577   -8.577   -0.845)   12.159   1.934   (   2.776   -2.776  -18.283)   18.700   2.478   (  10.290  -10.290    1.182)   14.600   2.999   (   4.833   -4.833   13.813)   15.411   3.221   (   4.310   -4.310   -5.237)    8.036   3.805   (   3.965   -3.965   -5.444)    7.815   4.157   (   2.919   -2.919    2.012)    4.593   4.263   (  13.424  -13.424    2.980)   19.217   4.708   (   0.514   -0.514    3.701)    3.772   5.619   (  -3.575    3.575    3.417)    6.103   5.796   (  -4.291    4.291    7.044)    9.297   6.764   (   8.314   -8.314    5.359)   12.921   7.313   (   4.553   -4.553   -3.791)    7.472   7.396   (   6.882   -6.882   -4.332)   10.653  10.864   (  -7.187    7.187   -0.345)   10.170  11.141   (  -2.328    2.328    2.808)    4.328  18.458   (  -7.539    7.539   -5.397)   11.949  20.082   (  -9.628    9.628   -2.755)   13.893======================= Grid point 281 (103/126) =======================q-point: ( 0.23 -0.23  0.23)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.500   (   0.000   17.113    0.000)   17.113   2.558   (   0.000   15.989    0.000)   15.989   2.657   (   0.000    4.781    0.000)    4.781   2.993   (  -0.000   -5.799   -0.000)    5.799   3.380   (  -0.000   -3.425   -0.000)    3.425   3.835   (   0.000    9.860    0.000)    9.860   4.397   (  -0.000  -19.083   -0.000)   19.083   5.203   (   0.000    3.287    0.000)    3.287   5.203   (   0.000    6.332    0.000)    6.332   5.614   (   0.000    2.915    0.000)    2.915   5.693   (  -0.000  -13.192   -0.000)   13.192   6.480   (   0.000    1.499    0.000)    1.499   6.612   (  -0.000   -7.227   -0.000)    7.227   8.726   (   0.000   25.677    0.000)   25.677  10.334   (  -0.000  -22.283   -0.000)   22.283  11.179   (  -0.000   -0.035   -0.000)    0.035  18.192   (  -0.000   -0.010   -0.000)    0.010  19.300   (   0.000    9.399    0.000)    9.399======================= Grid point 282 (104/126) =======================q-point: ( 0.23 -0.23  0.32)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.888   (   5.187   13.235    0.000)   14.215   2.608   (  -8.760    7.358    0.000)   11.440   2.746   (  -7.220   -8.868   -0.000)   11.435   2.910   (   2.269    2.524    0.000)    3.394   3.336   (   1.564   -0.278   -0.000)    1.589   3.955   (   4.144  -19.179   -0.000)   19.622   4.011   (   1.976    4.137    0.000)    4.584   5.110   ( -11.453   -0.943    0.000)   11.491   5.267   (   3.588   -8.592   -0.000)    9.312   5.502   (   1.213   -1.578   -0.000)    1.990   5.737   (   9.800   -5.779   -0.000)   11.377   6.618   (   7.117    4.448   -0.000)    8.392   6.645   (   9.714    1.541   -0.000)    9.835   8.920   (  -5.805   -6.243   -0.000)    8.525  10.074   (  -0.453    3.912    0.000)    3.939  11.159   (  -2.182   -0.562    0.000)    2.253  18.203   (  -1.426    2.453    0.000)    2.837  19.502   (  -1.471   11.029    0.000)   11.127======================= Grid point 283 (105/126) =======================q-point: (-0.77 -0.23  0.41)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.130   (  -1.403    2.011    0.000)    2.452   2.239   ( -10.521  -10.450    0.000)   14.829   2.730   (   1.064    6.911    0.000)    6.992   2.939   (   2.985    0.033   -0.000)    2.985   3.306   (  -0.618   -5.878   -0.000)    5.910   3.809   (   1.897   -1.089   -0.000)    2.188   4.195   (   8.013    1.084   -0.000)    8.086   4.788   ( -10.858   -5.294    0.000)   12.080   5.163   (  13.575  -14.245   -0.000)   19.677   5.579   (   2.240    1.577   -0.000)    2.739   5.668   (   2.871  -11.397   -0.000)   11.753   6.891   (   8.836    2.666   -0.000)    9.230   7.034   (  18.295    6.492   -0.000)   19.413   8.391   (  -7.546  -24.663   -0.000)   25.791  10.340   (  -0.553   16.429    0.000)   16.439  11.056   (  -6.965    0.057    0.000)    6.965  18.228   (  -3.309    4.504    0.000)    5.589  19.672   (  -4.541   10.722    0.000)   11.644======================= Grid point 285 (106/126) =======================q-point: ( 0.23  0.77 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.911   (   1.164   -1.164   -0.000)    1.646   2.181   (  -1.582    1.582    0.000)    2.237   2.802   (  -2.227    2.227    0.000)    3.150   3.025   (  11.217  -11.217   -0.000)   15.864   3.216   (  -4.729    4.729    0.000)    6.688   3.851   (  -1.176    1.176    0.000)    1.662   4.280   (   0.408   -0.408   -0.000)    0.577   4.596   (  -0.462    0.462    0.000)    0.653   5.200   (  21.263  -21.263   -0.000)   30.070   5.523   (   1.745   -1.745   -0.000)    2.468   5.625   (   2.239   -2.239   -0.000)    3.166   7.017   (   4.645   -4.645   -0.000)    6.569   7.509   (   2.682   -2.682   -0.000)    3.793   7.834   (  12.897  -12.897   -0.000)   18.239  10.554   (  -7.479    7.479    0.000)   10.577  10.942   (  -5.935    5.935    0.000)    8.393  18.242   (  -4.801    4.801    0.000)    6.790  19.739   (  -8.202    8.202    0.000)   11.600======================= Grid point 317 (107/126) =======================q-point: (-0.35 -0.25  0.25)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.769   (  -0.000  -30.296   -5.809)   30.848   0.997   (  -0.000  -37.571   -5.133)   37.920   2.316   (  -0.000    0.013  -24.918)   24.918   2.561   (   0.000   -6.803   12.016)   13.808   2.620   (  -0.000   -9.550   -0.313)    9.555   3.891   (   0.000   12.166    0.355)   12.171   4.130   (   0.000    2.516    9.177)    9.515   4.279   (   0.000   -5.851   11.503)   12.906   4.653   (   0.000   -4.879    1.702)    5.167   5.853   (   0.000    1.751    1.006)    2.020   5.947   (   0.000   -5.190    3.636)    6.337   6.702   (   0.000    0.095   17.847)   17.847   6.945   (  -0.000    3.354   -6.927)    7.696   7.094   (  -0.000  -10.480   -7.620)   12.957  11.170   (   0.000    2.855    1.671)    3.308  11.196   (  -0.000    0.134   -0.078)    0.155  18.690   (  -0.000    7.504  -10.526)   12.927  20.648   (   0.000   16.668    3.079)   16.950======================= Grid point 328 (108/126) =======================q-point: (-0.54 -0.26  0.26)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.974   (  -0.000   -1.117  -17.571)   17.606   2.181   (  -0.000  -26.011  -11.542)   28.457   2.587   (   0.000   -6.963    6.351)    9.425   2.840   (   0.000  -17.416   11.803)   21.038   3.121   (   0.000    4.473    1.447)    4.701   3.765   (   0.000    6.827    7.074)    9.831   3.903   (  -0.000  -12.564   -6.297)   14.054   4.427   (   0.000  -13.633    2.032)   13.784   4.923   (   0.000  -10.374    3.104)   10.829   5.723   (  -0.000    5.244   -1.427)    5.434   5.766   (   0.000    6.075    6.356)    8.793   6.578   (   0.000   -3.386    6.759)    7.560   6.788   (  -0.000    7.626   -2.597)    8.056   8.005   (  -0.000  -25.637   -4.437)   26.018  10.741   (   0.000   15.345    2.723)   15.585  11.188   (  -0.000    0.261   -0.030)    0.263  18.432   (  -0.000    9.486   -5.137)   10.788  20.028   (  -0.000   13.810   -2.762)   14.084======================= Grid point 329 (109/126) =======================q-point: (-0.54 -0.26  0.35)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 605Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.623   (   5.486  -26.627   -4.474)   27.552   1.847   (   5.255  -12.049  -15.010)   19.952   2.360   (   0.487   -6.218    5.708)    8.455   2.765   (  15.101  -14.701   13.154)   24.843   3.154   (   1.736   -2.372  -12.760)   13.094   3.724   (   2.171   -3.586   -2.779)    5.030   3.991   (   7.712    4.030    5.149)   10.111   4.218   (  -1.702   -7.812    4.318)    9.087   4.721   (  -1.253   -5.634    2.535)    6.303   5.719   (  -5.632    2.876    1.652)    6.536   5.868   (  -2.109    4.907    6.081)    8.093   6.575   (   8.532   -7.106    8.156)   13.777   7.047   (  10.186    6.875   -3.598)   12.805   7.484   (  -4.324  -19.344   -4.468)   20.319  10.969   (  -3.237    9.267    0.703)    9.841  11.180   (  -1.260    0.672    1.241)    1.892  18.598   (  -5.339   10.308   -6.347)   13.230  20.249   (  -6.027   12.231   -1.120)   13.681======================= Grid point 339 (110/126) =======================q-point: ( 0.27 -0.27  0.27)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.830   (   0.000   14.405    0.000)   14.405   2.773   (   0.000    6.069    0.000)    6.069   2.789   (   0.000    4.503    0.000)    4.503   2.851   (  -0.000   -7.323   -0.000)    7.323   3.314   (  -0.000   -1.492   -0.000)    1.492   3.905   (  -0.000  -25.766   -0.000)   25.766   3.988   (   0.000    5.124    0.000)    5.124   5.180   (  -0.000   -7.635   -0.000)    7.635   5.403   (   0.000   10.790    0.000)   10.790   5.411   (  -0.000  -13.528   -0.000)   13.528   5.598   (  -0.000   -2.576   -0.000)    2.576   6.537   (   0.000    4.134    0.000)    4.134   6.545   (   0.000    0.241    0.000)    0.241   8.984   (  -0.000   -5.685   -0.000)    5.685  10.079   (   0.000    2.758    0.000)    2.758  11.179   (   0.000    0.101    0.000)    0.101  18.218   (   0.000    2.706    0.000)    2.706  19.518   (   0.000   11.545    0.000)   11.545======================= Grid point 340 (111/126) =======================q-point: (-0.73 -0.27  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.107   (   1.976    7.528    0.000)    7.783   2.509   ( -13.260  -12.436    0.000)   18.179   2.723   (  -0.254   -0.863   -0.000)    0.900   2.894   (   0.919   -1.720   -0.000)    1.950   3.266   (   4.925   -9.759   -0.000)   10.931   3.680   (   9.151   -1.907   -0.000)    9.348   4.069   (   2.777    0.766   -0.000)    2.881   4.921   (   5.698  -19.121   -0.000)   19.952   5.033   (  -7.944   -9.972   -0.000)   12.749   5.557   (  -0.840    5.034    0.000)    5.104   5.646   (   0.954   -1.631   -0.000)    1.890   6.687   (  10.714    0.572   -0.000)   10.730   6.738   (   8.281    7.536   -0.000)   11.197   8.535   (  -5.461  -25.340   -0.000)   25.921  10.350   (  -0.395   17.687    0.000)   17.691  11.155   (  -2.827    0.067    0.000)    2.828  18.287   (  -2.163    5.817    0.000)    6.206  19.747   (  -2.584   12.424    0.000)   12.689======================= Grid point 342 (112/126) =======================q-point: (-0.73 -0.27  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.972   (  -7.020  -10.555   -0.000)   12.676   2.193   (   0.176   -0.136   -0.000)    0.222   2.749   (  12.512  -11.283   -0.000)   16.848   2.855   (  -1.076   -0.682    0.000)    1.274   3.258   (   1.906    0.999   -0.000)    2.152   3.843   (   1.909    2.673    0.000)    3.284   4.187   (   7.861   -1.753   -0.000)    8.055   4.691   (  -3.119   -5.818   -0.000)    6.601   4.740   (  15.168  -22.258   -0.000)   26.935   5.547   (  -3.363   -0.864    0.000)    3.472   5.633   (  -1.875    3.400    0.000)    3.883   6.866   (   9.592   -4.739   -0.000)   10.699   7.204   (  17.809    9.054   -0.000)   19.978   7.863   (  -7.148  -24.675   -0.000)   25.690  10.657   (  -3.106   13.129    0.000)   13.492  11.072   (  -5.662    1.454    0.000)    5.846  18.350   (  -5.270    7.118    0.000)    8.857  19.897   (  -6.757   10.702    0.000)   12.656======================= Grid point 376 (113/126) =======================q-point: (-0.30 -0.30  0.30)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 105Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.314   (  -0.000    0.000   -6.376)    6.376   0.314   (  -0.000    0.000   -6.376)    6.376   2.363   (  -0.000    0.000  -38.341)   38.341   2.490   (   0.000   -0.000   11.947)   11.947   2.490   (   0.000   -0.000   11.947)   11.947   4.092   (   0.000   -0.000    6.137)    6.137   4.156   (   0.000   -0.000    8.823)    8.823   4.156   (   0.000   -0.000    8.823)    8.823   4.604   (   0.000   -0.000    0.398)    0.398   5.871   (   0.000   -0.000    1.435)    1.435   5.871   (   0.000   -0.000    1.435)    1.435   6.711   (   0.000   -0.000   18.568)   18.568   6.980   (  -0.000    0.000   -7.147)    7.147   6.980   (  -0.000    0.000   -7.147)    7.147  11.197   (  -0.000    0.000   -0.082)    0.082  11.197   (  -0.000    0.000   -0.082)    0.082  18.773   (  -0.000    0.000  -10.873)   10.873  20.845   (   0.000   -0.000   10.576)   10.576======================= Grid point 387 (114/126) =======================q-point: (-0.49 -0.31  0.31)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.564   (  -0.000  -32.385   -8.300)   33.432   1.775   (  -0.000  -21.189  -11.460)   24.090   2.418   (   0.000   -7.640    8.429)   11.376   2.534   (   0.000  -11.939   11.152)   16.337   3.137   (  -0.000   -3.666  -15.654)   16.078   3.712   (  -0.000   -0.863   -0.112)    0.871   3.900   (   0.000    5.675    5.625)    7.990   4.240   (   0.000   -6.594    4.650)    8.069   4.738   (   0.000   -7.261    2.404)    7.649   5.804   (  -0.000    2.767   -0.422)    2.799   5.876   (   0.000    3.784    5.671)    6.818   6.473   (   0.000   -5.754   10.051)   11.581   6.939   (  -0.000    6.504   -3.303)    7.294   7.530   (  -0.000  -19.625   -4.467)   20.127  11.003   (   0.000    9.839    1.726)    9.989  11.193   (  -0.000    0.236   -0.039)    0.240  18.657   (  -0.000   11.503   -6.687)   13.306  20.315   (   0.000   13.388   -0.154)   13.389======================= Grid point 388 (115/126) =======================q-point: (-0.49 -0.31  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 336Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.103   (  20.888  -20.888   -3.159)   29.708   1.444   (  24.454  -24.454   -8.920)   35.715   2.286   (   0.898   -0.898    7.567)    7.673   2.484   (  10.584  -10.584   10.625)   18.357   3.086   (   2.750   -2.750  -25.212)   25.510   3.752   (  -4.199    4.199    5.287)    7.951   4.057   (   1.664   -1.664    4.172)    4.790   4.092   (   1.214   -1.214    6.551)    6.772   4.644   (   1.907   -1.907    1.337)    3.011   5.782   (  -3.357    3.357    1.205)    4.898   5.927   (  -0.118    0.118    3.403)    3.407   6.438   (   4.894   -4.894   11.468)   13.395   7.149   (   2.986   -2.986   -3.906)    5.752   7.179   (   3.999   -3.999   -4.000)    6.927  11.108   (  -4.243    4.243   -0.112)    6.002  11.193   (  -0.434    0.434    1.050)    1.216  18.791   (  -7.388    7.388   -7.367)   12.784  20.468   (  -8.124    8.124    2.721)   11.807======================= Grid point 398 (116/126) =======================q-point: (-0.68 -0.32  0.32)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.082   (   0.000    9.358    0.000)    9.358   2.682   (  -0.000   -8.483   -0.000)    8.483   2.716   (  -0.000  -10.472   -0.000)   10.472   2.921   (   0.000    7.743    0.000)    7.743   3.156   (  -0.000  -22.214   -0.000)   22.214   3.564   (   0.000    0.944    0.000)    0.944   4.048   (   0.000    0.354    0.000)    0.354   4.856   (  -0.000  -21.382   -0.000)   21.382   5.132   (  -0.000  -12.942   -0.000)   12.942   5.604   (   0.000    8.135    0.000)    8.135   5.605   (   0.000    3.094    0.000)    3.094   6.561   (  -0.000   -0.087   -0.000)    0.087   6.653   (   0.000    6.942    0.000)    6.942   8.595   (  -0.000  -25.801   -0.000)   25.801  10.355   (   0.000   18.130    0.000)   18.130  11.183   (   0.000    0.213    0.000)    0.213  18.311   (   0.000    6.370    0.000)    6.370  19.775   (   0.000   13.042    0.000)   13.042======================= Grid point 400 (117/126) =======================q-point: (-0.68 -0.32  0.41)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 363Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.176   (  -5.975   -8.162   -0.000)   10.115   2.197   (  -6.555  -12.572   -0.000)   14.179   2.508   (   9.117  -14.189   -0.000)   16.865   2.840   (   4.939   -7.800   -0.000)    9.232   3.153   (   5.769    2.904   -0.000)    6.459   3.766   (   5.075    6.292    0.000)    8.083   4.020   (   6.130   -6.424   -0.000)    8.880   4.475   (   5.673  -20.020   -0.000)   20.808   4.837   (  -4.835   -8.552   -0.000)    9.824   5.651   (  -5.915    2.150    0.000)    6.294   5.678   (  -2.140    5.629    0.000)    6.022   6.643   (  10.634   -4.872   -0.000)   11.696   6.913   (   9.560    8.722   -0.000)   12.941   7.996   (  -4.837  -25.175   -0.000)   25.636  10.701   (  -1.380   15.166    0.000)   15.228  11.160   (  -2.750    0.380    0.000)    2.776  18.446   (  -3.576    9.361    0.000)   10.021  20.009   (  -3.871   12.553    0.000)   13.136======================= Grid point 402 (118/126) =======================q-point: ( 0.32  0.68 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.781   (   7.589   -7.589   -0.000)   10.732   2.154   (   3.792   -3.792   -0.000)    5.363   2.465   (  12.429  -12.429   -0.000)   17.577   2.814   (   6.039   -6.039   -0.000)    8.540   3.267   (   0.209   -0.209   -0.000)    0.295   3.864   (   1.871   -1.871   -0.000)    2.646   4.138   (   3.810   -3.810   -0.000)    5.387   4.249   (  17.103  -17.103   -0.000)   24.187   4.671   (  -0.337    0.337    0.000)    0.476   5.587   (  -4.191    4.191    0.000)    5.927   5.717   (  -4.387    4.387    0.000)    6.204   6.711   (   9.661   -9.661   -0.000)   13.662   7.350   (   4.425   -4.425   -0.000)    6.258   7.438   (   6.720   -6.720   -0.000)    9.503  10.867   (  -7.103    7.103    0.000)   10.045  11.115   (  -2.649    2.649    0.000)    3.746  18.513   (  -8.117    8.117    0.000)   11.479  20.104   (  -8.959    8.959    0.000)   12.670======================= Grid point 447 (119/126) =======================q-point: (-0.45 -0.35  0.35)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 330Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.880   (  -0.000  -31.980   -5.283)   32.413   1.093   (  -0.000  -40.780   -4.578)   41.036   2.310   (   0.000   -3.087    8.353)    8.905   2.344   (   0.000   -6.262    9.647)   11.501   3.058   (  -0.000   -2.440  -28.612)   28.716   3.816   (   0.000    7.600    7.682)   10.806   3.992   (   0.000    2.998    5.072)    5.892   4.113   (   0.000   -5.771    5.881)    8.239   4.628   (   0.000   -3.241    0.916)    3.368   5.844   (   0.000    1.187    0.173)    1.199   5.894   (   0.000   -1.992    2.138)    2.922   6.391   (   0.000   -1.440   13.071)   13.150   7.046   (  -0.000    3.648   -3.641)    5.153   7.206   (  -0.000  -11.214   -3.993)   11.904  11.152   (   0.000    4.585    0.513)    4.613  11.197   (  -0.000    0.140   -0.046)    0.147  18.873   (  -0.000    8.255   -7.764)   11.332  20.558   (   0.000    9.151    4.978)   10.418======================= Grid point 458 (120/126) =======================q-point: (-0.64 -0.36  0.36)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   2.190   (  -0.000   -0.323   -0.000)    0.323   2.335   (  -0.000  -27.694   -0.000)   27.694   2.479   (  -0.000  -11.380   -0.000)   11.380   2.681   (  -0.000  -17.419   -0.000)   17.419   3.101   (   0.000    9.539    0.000)    9.539   3.698   (   0.000    7.942    0.000)    7.942   3.959   (  -0.000  -10.200   -0.000)   10.200   4.428   (  -0.000  -16.271   -0.000)   16.271   4.888   (  -0.000  -10.317   -0.000)   10.317   5.698   (   0.000    5.446    0.000)    5.446   5.736   (   0.000    4.707    0.000)    4.707   6.511   (  -0.000   -4.830   -0.000)    4.830   6.813   (   0.000    8.051    0.000)    8.051   8.048   (  -0.000  -25.458   -0.000)   25.458  10.716   (   0.000   15.658    0.000)   15.658  11.188   (   0.000    0.266    0.000)    0.266  18.485   (   0.000   10.318    0.000)   10.318  20.050   (   0.000   13.247    0.000)   13.247======================= Grid point 460 (121/126) =======================q-point: (-0.64 -0.36  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 366Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.671   (   2.640  -27.738   -0.000)   27.864   2.031   (   7.797  -15.447   -0.000)   17.303   2.295   (   3.324   -6.603   -0.000)    7.392   2.592   (  13.060  -13.932   -0.000)   19.097   3.281   (  -1.703    7.235    0.000)    7.433   3.767   (   5.458   -9.193   -0.000)   10.691   3.934   (   8.140    4.124   -0.000)    9.125   4.182   (  -1.478   -9.394   -0.000)    9.510   4.695   (  -1.412   -4.954   -0.000)    5.151   5.703   (  -6.442    2.901    0.000)    7.065   5.799   (  -2.999    5.571    0.000)    6.327   6.497   (   9.824   -8.575   -0.000)   13.040   7.082   (  10.416    7.162   -0.000)   12.641   7.526   (  -4.338  -19.472   -0.000)   19.949  10.963   (  -2.549    9.864    0.000)   10.188  11.170   (  -2.162    0.536    0.000)    2.228  18.663   (  -5.668   10.939    0.000)   12.320  20.255   (  -5.213   10.687    0.000)   11.891======================= Grid point 508 (122/126) =======================q-point: (-0.40 -0.40  0.40)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 126Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.427   (  -0.000    0.000   -4.963)    4.963   0.427   (  -0.000    0.000   -4.963)    4.963   2.278   (   0.000   -0.000    9.145)    9.145   2.278   (   0.000   -0.000    9.145)    9.145   3.036   (  -0.000    0.000  -31.468)   31.468   3.932   (   0.000   -0.000   11.228)   11.228   4.023   (   0.000   -0.000    4.913)    4.913   4.023   (   0.000   -0.000    4.913)    4.913   4.597   (   0.000   -0.000    0.247)    0.247   5.856   (   0.000   -0.000    0.366)    0.366   5.856   (   0.000   -0.000    0.366)    0.366   6.387   (   0.000   -0.000   13.800)   13.800   7.085   (  -0.000    0.000   -3.736)    3.736   7.085   (  -0.000    0.000   -3.736)    3.736  11.198   (  -0.000    0.000   -0.048)    0.048  11.198   (  -0.000    0.000   -0.048)    0.048  18.964   (  -0.000    0.000   -8.218)    8.218  20.659   (   0.000   -0.000    8.024)    8.024======================= Grid point 521 (123/126) =======================q-point: (-0.59 -0.41  0.41)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.659   (  -0.000  -34.991   -0.000)   34.991   1.914   (  -0.000  -28.180   -0.000)   28.180   2.297   (  -0.000   -4.846   -0.000)    4.846   2.409   (  -0.000  -10.246   -0.000)   10.246   3.307   (   0.000    9.658    0.000)    9.658   3.705   (  -0.000   -9.856   -0.000)    9.856   3.847   (   0.000    6.045    0.000)    6.045   4.202   (  -0.000   -7.776   -0.000)    7.776   4.712   (  -0.000   -6.495   -0.000)    6.495   5.808   (   0.000    2.405    0.000)    2.405   5.813   (   0.000    4.987    0.000)    4.987   6.375   (  -0.000   -7.354   -0.000)    7.354   6.970   (   0.000    6.728    0.000)    6.728   7.573   (  -0.000  -19.765   -0.000)   19.765  10.988   (   0.000   10.391    0.000)   10.391  11.193   (   0.000    0.240    0.000)    0.240  18.725   (   0.000   12.161    0.000)   12.161  20.310   (   0.000   11.360    0.000)   11.360======================= Grid point 523 (124/126) =======================q-point: ( 0.41  0.59 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 198Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   1.137   (  20.724  -20.724   -0.000)   29.308   1.546   (  27.995  -27.995   -0.000)   39.591   2.207   (   1.694   -1.694   -0.000)    2.396   2.356   (   8.327   -8.327   -0.000)   11.777   3.405   (  -3.860    3.860    0.000)    5.458   3.644   (   1.654   -1.654   -0.000)    2.338   4.017   (   1.965   -1.965   -0.000)    2.779   4.026   (   1.214   -1.214   -0.000)    1.717   4.631   (   1.469   -1.469   -0.000)    2.078   5.772   (  -3.721    3.721    0.000)    5.263   5.885   (  -1.869    1.869    0.000)    2.643   6.327   (   6.556   -6.556   -0.000)    9.271   7.186   (   3.028   -3.028   -0.000)    4.283   7.217   (   4.148   -4.148   -0.000)    5.866  11.109   (  -4.211    4.211    0.000)    5.956  11.183   (  -0.842    0.842    0.000)    1.190  18.867   (  -7.786    7.786    0.000)   11.011  20.436   (  -6.320    6.320    0.000)    8.938======================= Grid point 583 (125/126) =======================q-point: (-0.55 -0.45  0.45)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 216Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.937   (  -0.000  -33.091   -0.000)   33.091   1.144   (  -0.000  -42.398   -0.000)   42.398   2.218   (  -0.000   -2.995   -0.000)    2.995   2.243   (  -0.000   -5.585   -0.000)    5.585   3.439   (   0.000    1.927    0.000)    1.927   3.648   (   0.000    3.916    0.000)    3.916   3.943   (   0.000    3.142    0.000)    3.142   4.056   (  -0.000   -6.264   -0.000)    6.264   4.619   (  -0.000   -2.587   -0.000)    2.587   5.843   (   0.000    1.012    0.000)    1.012   5.870   (  -0.000   -0.041   -0.000)    0.041   6.260   (  -0.000   -2.653   -0.000)    2.653   7.081   (   0.000    3.737    0.000)    3.737   7.244   (  -0.000  -11.429   -0.000)   11.429  11.147   (   0.000    5.027    0.000)    5.027  11.197   (   0.000    0.142    0.000)    0.142  18.953   (   0.000    8.715    0.000)    8.715  20.503   (   0.000    6.765    0.000)    6.765======================= Grid point 646 (126/126) =======================q-point: ( 0.50  0.50 -0.50)Boundary mean free path (millimeter): 1000.000Mass variance parameters: 1.64e-04 3.36e-05 3.36e-05 0.00e+00 0.00e+00 2.11e-04 Number of triplets: 63Calculating Gamma of ph-isotope with tetrahedron methodFrequency     group velocity (x, y, z)     |gv|   0.479   (   0.000   -0.000    0.000)    0.000   0.479   (   0.000   -0.000    0.000)    0.000   2.184   (  -0.000    0.000   -0.000)    0.000   2.184   (  -0.000    0.000   -0.000)    0.000   3.448   (  -0.000    0.000   -0.000)    0.000   3.722   (  -0.000    0.000   -0.000)    0.000   3.976   (  -0.000    0.000   -0.000)    0.000   3.976   (  -0.000    0.000   -0.000)    0.000   4.595   (  -0.000    0.000   -0.000)    0.000   5.853   (  -0.000    0.000   -0.000)    0.000   5.853   (  -0.000    0.000   -0.000)    0.000   6.246   (  -0.000    0.000   -0.000)    0.000   7.120   (  -0.000    0.000   -0.000)    0.000   7.120   (  -0.000    0.000   -0.000)    0.000  11.199   (  -0.000    0.000   -0.000)    0.000  11.199   (  -0.000    0.000   -0.000)    0.000  19.049   (  -0.000    0.000   -0.000)    0.000  20.576   (  -0.000    0.000   -0.000)    0.000=================== End of collection of collisions ===================----------- Thermal conductivity (W/m-k) with tetrahedron method -----------#  T(K)        xx         yy         zz         yz         xz         xy        #ipm    0.0      0.000      0.000      0.000      0.000      0.000      0.000 3/21780   10.0     22.935     22.935     13.253     -0.000      0.000      0.000 3/21780   20.0     24.355     24.355     11.787     -0.000      0.000      0.000 3/21780   30.0     16.990     16.990      8.396     -0.000      0.000      0.000 3/21780   40.0     12.149     12.149      6.224     -0.000      0.000      0.000 3/21780   50.0      9.226      9.226      4.860     -0.000      0.000      0.000 3/21780   60.0      7.387      7.387      3.963     -0.000      0.000      0.000 3/21780   70.0      6.166      6.166      3.343     -0.000      0.000      0.000 3/21780   80.0      5.311      5.311      2.893     -0.000      0.000      0.000 3/21780   90.0      4.683      4.683      2.554     -0.000      0.000      0.000 3/21780  100.0      4.203      4.203      2.290     -0.000      0.000      0.000 3/21780  110.0      3.822      3.822      2.078     -0.000      0.000      0.000 3/21780  120.0      3.512      3.512      1.906     -0.000      0.000      0.000 3/21780  130.0      3.253      3.253      1.762     -0.000      0.000      0.000 3/21780  140.0      3.033      3.033      1.642     -0.000      0.000      0.000 3/21780  150.0      2.844      2.844      1.539     -0.000      0.000      0.000 3/21780  160.0      2.679      2.679      1.451     -0.000      0.000      0.000 3/21780  170.0      2.533      2.533      1.374     -0.000      0.000      0.000 3/21780  180.0      2.403      2.403      1.306     -0.000      0.000      0.000 3/21780  190.0      2.287      2.287      1.246     -0.000      0.000      0.000 3/21780  200.0      2.182      2.182      1.193     -0.000      0.000      0.000 3/21780  210.0      2.087      2.087      1.145     -0.000      0.000      0.000 3/21780  220.0      2.000      2.000      1.102     -0.000      0.000      0.000 3/21780  230.0      1.921      1.921      1.062     -0.000      0.000      0.000 3/21780  240.0      1.847      1.847      1.026     -0.000      0.000      0.000 3/21780  250.0      1.780      1.780      0.993     -0.000      0.000      0.000 3/21780  260.0      1.717      1.717      0.962     -0.000      0.000      0.000 3/21780  270.0      1.658      1.658      0.933     -0.000      0.000      0.000 3/21780  280.0      1.604      1.604      0.906     -0.000      0.000      0.000 3/21780  290.0      1.553      1.553      0.881     -0.000      0.000      0.000 3/21780  300.0      1.505      1.505      0.858     -0.000      0.000      0.000 3/21780  310.0      1.460      1.460      0.835     -0.000      0.000      0.000 3/21780  320.0      1.418      1.418      0.814     -0.000      0.000      0.000 3/21780  330.0      1.378      1.378      0.794     -0.000      0.000      0.000 3/21780  340.0      1.340      1.340      0.775     -0.000      0.000      0.000 3/21780  350.0      1.305      1.305      0.757     -0.000      0.000      0.000 3/21780  360.0      1.271      1.271      0.740     -0.000      0.000      0.000 3/21780  370.0      1.239      1.239      0.724     -0.000      0.000      0.000 3/21780  380.0      1.208      1.208      0.708     -0.000      0.000      0.000 3/21780  390.0      1.179      1.179      0.693     -0.000      0.000      0.000 3/21780  400.0      1.151      1.151      0.678     -0.000      0.000      0.000 3/21780  410.0      1.125      1.125      0.665     -0.000      0.000      0.000 3/21780  420.0      1.100      1.100      0.651     -0.000      0.000      0.000 3/21780  430.0      1.076      1.076      0.639     -0.000      0.000      0.000 3/21780  440.0      1.052      1.052      0.626     -0.000      0.000      0.000 3/21780  450.0      1.030      1.030      0.615     -0.000      0.000      0.000 3/21780  460.0      1.009      1.009      0.603     -0.000      0.000      0.000 3/21780  470.0      0.989      0.989      0.592     -0.000      0.000      0.000 3/21780  480.0      0.969      0.969      0.582     -0.000      0.000      0.000 3/21780  490.0      0.950      0.950      0.571     -0.000      0.000      0.000 3/21780  500.0      0.932      0.932      0.562     -0.000      0.000      0.000 3/21780  510.0      0.915      0.915      0.552     -0.000      0.000      0.000 3/21780  520.0      0.898      0.898      0.543     -0.000      0.000      0.000 3/21780  530.0      0.882      0.882      0.534     -0.000      0.000      0.000 3/21780  540.0      0.866      0.866      0.525     -0.000      0.000      0.000 3/21780  550.0      0.851      0.851      0.517     -0.000      0.000      0.000 3/21780  560.0      0.836      0.836      0.509     -0.000      0.000      0.000 3/21780  570.0      0.822      0.822      0.501     -0.000      0.000      0.000 3/21780  580.0      0.808      0.808      0.493     -0.000      0.000      0.000 3/21780  590.0      0.795      0.795      0.486     -0.000      0.000      0.000 3/21780  600.0      0.782      0.782      0.478     -0.000      0.000      0.000 3/21780  610.0      0.770      0.770      0.471     -0.000      0.000      0.000 3/21780  620.0      0.758      0.758      0.465     -0.000      0.000      0.000 3/21780  630.0      0.746      0.746      0.458     -0.000      0.000      0.000 3/21780  640.0      0.735      0.735      0.452     -0.000      0.000      0.000 3/21780  650.0      0.724      0.724      0.445     -0.000      0.000      0.000 3/21780  660.0      0.714      0.714      0.439     -0.000      0.000      0.000 3/21780  670.0      0.703      0.703      0.433     -0.000      0.000      0.000 3/21780  680.0      0.693      0.693      0.427     -0.000      0.000      0.000 3/21780  690.0      0.683      0.683      0.422     -0.000      0.000      0.000 3/21780  700.0      0.674      0.674      0.416     -0.000      0.000      0.000 3/21780  710.0      0.665      0.665      0.411     -0.000      0.000      0.000 3/21780  720.0      0.656      0.656      0.406     -0.000      0.000      0.000 3/21780  730.0      0.647      0.647      0.401     -0.000      0.000      0.000 3/21780  740.0      0.638      0.638      0.396     -0.000      0.000      0.000 3/21780  750.0      0.630      0.630      0.391     -0.000      0.000      0.000 3/21780  760.0      0.622      0.622      0.386     -0.000      0.000      0.000 3/21780  770.0      0.614      0.614      0.381     -0.000      0.000      0.000 3/21780  780.0      0.607      0.607      0.377     -0.000      0.000      0.000 3/21780  790.0      0.599      0.599      0.372     -0.000      0.000      0.000 3/21780  800.0      0.592      0.592      0.368     -0.000      0.000      0.000 3/21780  810.0      0.585      0.585      0.364     -0.000      0.000      0.000 3/21780  820.0      0.578      0.578      0.360     -0.000      0.000      0.000 3/21780  830.0      0.571      0.571      0.356     -0.000      0.000      0.000 3/21780  840.0      0.564      0.564      0.352     -0.000      0.000      0.000 3/21780  850.0      0.558      0.558      0.348     -0.000      0.000      0.000 3/21780  860.0      0.551      0.551      0.344     -0.000      0.000      0.000 3/21780  870.0      0.545      0.545      0.340     -0.000      0.000      0.000 3/21780  880.0      0.539      0.539      0.337     -0.000      0.000      0.000 3/21780  890.0      0.533      0.533      0.333     -0.000      0.000      0.000 3/21780  900.0      0.527      0.527      0.330     -0.000      0.000      0.000 3/21780  910.0      0.522      0.522      0.326     -0.000      0.000      0.000 3/21780  920.0      0.516      0.516      0.323     -0.000      0.000      0.000 3/21780  930.0      0.511      0.511      0.320     -0.000      0.000      0.000 3/21780  940.0      0.505      0.505      0.316     -0.000      0.000      0.000 3/21780  950.0      0.500      0.500      0.313     -0.000      0.000      0.000 3/21780  960.0      0.495      0.495      0.310     -0.000      0.000      0.000 3/21780  970.0      0.490      0.490      0.307     -0.000      0.000      0.000 3/21780  980.0      0.485      0.485      0.304     -0.000      0.000      0.000 3/21780  990.0      0.480      0.480      0.301     -0.000      0.000      0.000 3/21780 1000.0      0.475      0.475      0.298     -0.000      0.000      0.000 3/21780Thermal conductivity related properties were written into "kappa-m111110.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 10:42:26]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|